Seedling root traits impact plant establishment under challenging environments. Pearl millet is one of the most heat and drought tolerant cereal crops that provides a vital food source across the sub-Saharan Sahel region. Pearl millet’s early root system features a single fast-growing primary root which we hypothesize is an adaptation to the Sahelian climate. Using crop modeling, we demonstrate that early drought stress is an important constraint in agrosystems in the Sahel where pearl millet was domesticated. Furthermore, we show that increased pearl millet primary root growth is correlated with increased early water stress tolerance in field conditions. Genetics including genome-wide association study and quantitative trait loci (QTL) approaches identify genomic regions controlling this key root trait. Combining gene expression data, re-sequencing and re-annotation of one of these genomic regions identified a glutaredoxin-encoding gene PgGRXC9 as the candidate stress resilience root growth regulator. Functional characterization of its closest Arabidopsis homolog AtROXY19 revealed a novel role for this glutaredoxin (GRX) gene clade in regulating cell elongation. In summary, our study suggests a conserved function for GRX genes in conferring root cell elongation and enhancing resilience of pearl millet to its Sahelian environment.

The rhizosheath, the layer of soil that adheres strongly to roots, influences water and nutrients acquisition. Pearl millet is a cereal crop that plays a major role for food security in arid regions of sub-Saharan Africa and India. We previously showed that root-adhering soil mass is a heritable trait in pearl millet and that it correlates with changes in rhizosphere microbiota structure and functions. Here, we studied the correlation between root-adhering soil mass and root hair development, root architecture, and symbiosis with arbuscular mycorrhizal fungi and we analysed the genetic control of this trait using genome wide association (GWAS) combined with bulk segregant analysis and gene expression studies. Root-adhering soil mass was weakly correlated only to root hairs traits in pearl millet. Twelve QTLs for rhizosheath formation were identified by GWAS. Bulk segregant analysis on a biparental population validated five of these QTLs. Combining genetics with a comparison of global gene expression in the root tip of contrasted inbred lines revealed candidate genes that might control rhizosheath formation in pearl millet. Our study indicates that rhizosheath formation is under complex genetic control in pearl millet and suggests that it is mainly regulated by root exudation.

Ecosystems worldwide are suffering from climate change. Coral reef ecosystems are globally threatened by increasing sea surface temperatures. However, gene expression plasticity provides the potential for organisms to respond rapidly and effectively to environmental changes, and would be favored in variable environments. In this study, we investigated the thermal stress response in Pocillopora coral colonies from two contrasting environments by exposing them to heat stress. We compared the physiological state, bacterial and Symbiodiniaceae communities (using 16S and ITS2 metabarcoding), and gene expression levels (using RNA-Seq) between control conditions and heat stress (the temperature just below the first signs of compromised health). Colonies from both thermal regimes remained apparently normal and presented open and colored polyps during heat stress, with no change in bacterial and Symbiodiniaceae community composition. In contrast, they differed in their transcriptomic responses. The colonies from Oman displayed a more plastic transcriptome, but some genes had a higher basal expression level (frontloading) compared to the less thermotolerant colonies from New Caledonia. In terms of biological functions, we observed an increase in the expression of stress response genes (including induction of tumor necrosis factor receptors, heat shock proteins, and detoxification of reactive oxygen species), together with a decrease in the expression of genes involved in morpho-anatomical functions. Gene regulation (transcription factors, mobile elements, histone modifications and DNA methylation) appeared to be overrepresented in the Oman colonies, indicating possible epigenetic regulation. These results show that transcriptomic plasticity and frontloading can be co-occurring processes in corals confronted with highly variable thermal regimes.

Correctly delimiting species and populations is a prerequisite for studies of connectivity, adaptation and conservation. Genomic data are particularly useful to test species differentiation for organisms with few informative morphological characters or low discrimination of cytoplasmic markers, as in Scleractinians. Here we applied Restriction site Associated DNA sequencing (RAD-sequencing) to the study of species differentiation and genetic structure in populations of Pocillopora spp. from Oman and French Polynesia, with the objectives to test species hypotheses, and to study the genetic structure among sampling sites within species. We focused here on coral colonies morphologically similar to P. acuta (damicornis type β). We tested the impact of different filtering strategies on the stability of the results. The main genetic differentiation was observed between samples from Oman and French Polynesia. These samples corresponded to different previously defined primary species hypotheses (PSH), i.e., PSHs 12 and 13 in Oman, and PSH 5 in French Polynesia. In Oman, we did not observe any clear differentiation between the two putative species PSH 12 and 13, nor between sampling sites. In French Polynesia, where a single species hypothesis was studied, there was no differentiation between sites. Our analyses allowed the identification of clonal lineages in Oman and French Polynesia. The impact of clonality on genetic diversity is discussed in light of individual-based simulations.

The freshwater snail Biomphalaria glabrata is an intermediate host of Schistosoma mansoni , the agent of human intestinal schistosomiasis. However, much is to be discovered about its innate immune system that appears as a complex black box, in which the immune cells (called hemocytes) play a major role in both cellular and humoral response towards pathogens. Until now, hemocyte classification has been based exclusively on cell morphology and ultrastructural description and depending on the authors considered from 2 to 5 hemocyte populations have been described. In this study, we proposed to evaluate the hemocyte heterogeneity at the transcriptomic level. To accomplish this objective, we used single cell RNA sequencing (scRNAseq) technology coupled to a droplet-based system to separate hemocytes and analyze their transcriptome at a unique cell level in naive Biomphalaria glabrata snails. We were able to demonstrate the presence of 7 hemocyte transcriptomic populations defined by the expression of specific marker genes. As a result, scRNAseq approach showed a high heterogeneity within hemocytes, but provides a detailed description of the different hemocyte transcriptomic populations in B. glabrata supported by distinct cellular functions and lineage trajectory. As a main result, scRNAseq revealed the 3 main population as a super-group of hemocyte diversity but, on the contrary, a great hemocytes plasticity with a probable capacity of hemocytes to engage to different activation pathways. This work opens a new field of research to understand the role of hemocytes particularly in response to pathogens, and towards S. mansoni parasites.

Rationale: Patients with colorectal cancer die mainly due to liver metastases (CRC-LM). Although the tumor microenvironment (TME) plays an important role in tumor development and therapeutic response, our understanding of the individual TME components, especially cancer-associated fibroblasts (CAFs), remains limited. Methods: We analyzed CRC-LM CAFs and cancer cells by single-cell transcriptomics and used bioinformatics for data analysis and integration with related available single-cell and bulk transcriptomic datasets. We validated key findings by RT-qPCR, western blotting, and immunofluorescence. Results: By single-cell transcriptomic analysis of 4,397 CAFs from six CRC-LM samples, we identified two main CAF populations, contractile CAFs and extracellular matrix (ECM)-remodeling/pro-angiogenic CAFs, and four subpopulations with distinct phenotypes. We found that ECM-remodeling/pro-angiogenic CAFs derive from portal resident fibroblasts. They associate with areas of strong desmoplastic reaction and Wnt signaling in low-proliferating tumor cells engulfed in a stiff extracellular matrix. By integrating public single-cell primary liver tumor data, we propose a model to explain how different liver malignancies recruit CAFs of different origins to this organ. Lastly, we found that LTBP2 plays an important role in modulating collagen biosynthesis, ECM organization, and adhesion pathways. We developed fully human antibodies against LTBP2 that depleted LTBP2+ CAFs in vitro. Conclusion: This study complements recent reports on CRC-LM CAF heterogeneity at the single-cell resolution. The number of sequenced CAFs was more than one order of magnitude larger compared to existing data. LTBP2 targeting by antibodies might create opportunities to deplete ECM-remodeling CAFs in CRC-LMs. This might be combined with other therapies, e.g., anti-angiogenic compounds as already done in CRC. Moreover, we showed that in intrahepatic cholangiocarcinoma, in which ECM-remodeling CAF proportion is similar to that of CRC-LM, several genes expressed by ECM-remodeling CAFs, such as LTBP2, were associated with survival.

Despite the growing interest in the rabbit model for developmental and stem cell biology, the characterization of embryos at the molecular level is still poorly documented. We conducted a transcriptome analysis of rabbit preimplantation embryos from E2.7 (morula stage) to E6.6 (early primitive streak stage) using bulk and single-cell RNA-sequencing. In parallel, we studied oxidative phosphorylation and glycolysis, and analysed active and repressive epigenetic modifications during blastocyst formation and expansion. We generated a transcriptomic, epigenetic and metabolic map of the pluripotency continuum in rabbit preimplantation embryos, and identified novel markers of naive pluripotency that might be instrumental for deriving naive pluripotent stem cell lines. Although the rabbit is evolutionarily closer to mice than to primates, we found that the transcriptome of rabbit epiblast cells shares common features with those of humans and non-human primates.

Bolder individuals have greater access to food sources and reproductive partners but are also at increased risk of predation. Boldness is believed to be consistent across time and contexts, but few studies have investigated the stability of this trait across variable environments, such as varying stress loads or long periods of time. Moreover, the underlying molecular components of boldness are poorly studied. Here, we report that boldness of 1154 European sea bass, evaluated using group risk-taking tests, is consistent over seven months and for individuals subjected to multiple environments, including a chronically stressful environment. Differences in risk-taking behaviour were further supported by differences observed in the responses to a novel environment test: shy individuals displayed more group dispersion, more thigmotaxic behaviour and lower activity levels. Transcriptomic analyses performed on extreme phenotypes revealed that bold individuals display greater expression for genes involved in social and exploration behaviours, and memory in the pituitary, and genes involved in immunity and responses to stimuli in the head kidney. This study demonstrates that personality traits come with an underpinning molecular signature, especially in organs involved in the endocrine and immune systems. As such, our results help to depict state–behaviour feedback mechanisms, previously proposed as key in shaping animal personality.

Fish species, such as zebrafish (Danio rerio), can regenerate their appendages after amputation through the formation of a heterogeneous cellular structure named blastema. Here, by combining live imaging of triple transgenic zebrafish embryos and single-cell RNA sequencing we established a detailed cell atlas of the regenerating caudal fin in zebrafish larvae. We confirmed the presence of macrophage subsets that govern zebrafish fin regeneration, and identified a foxd3-positive cell population within the regenerating fin. Genetic depletion of these foxd3-positive neural crest-derived cells (NCdC) showed that they are involved in blastema formation and caudal fin regeneration. Finally, chemical inhibition and transcriptomic analysis demonstrated that these foxd3-positive cells regulate macrophage recruitment and polarization through the NRG1/ErbB pathway. Here, we show the diversity of the cells required for blastema formation, identify a discrete foxd3-positive NCdC population, and reveal the critical function of the NRG1/ErbB pathway in controlling the dialogue between macrophages and NCdC.

Scleractinian corals are of great ecological interest as ecosystem engineer species. Accordingly, there is a wealth of studies on their adaptive abilities facing climate change. Such studies should rely on precise species and population delimitation. Nevertheless species delimitation in corals can be hindered by the lack of adequate genetic markers, by hybridization, and by morphological plasticity. Here we applied RAD sequencing to the study of species delimitation and genetic structure in populations of Pocillopora spp. from Oman and French Polynesia with the objectives to test primary species hypotheses based on mitochondrial DNA sequencing, and to study the genetic structure among sampling sites inside species. Regarding the varying levels of missing data observed among samples we tested different filtering strategy. The main genetic differentiation was observed between samples from Oman and French Polynesia, which also corresponded to different mitochondrial lineages and species hypotheses. In Oman, we did not observe any clear differentiation according to the main mitochondrial lineages considered here, nor between sampling sites. In French Polynesia where a single mitochondrial lineage was studied, we did not evidence any differentiation according to sampling sites. These results provide an additional example of the importance of using independent nuclear markers for the study of species delimitation. Our analyses also allowed the identification of clonal lineages in our samples, and to take them into account in our interpretations. We used simulations to study the impact of clonal reproduction on the distribution of statistics of genetic diversity and genetic structure among loci.

Despite the growing interest in the rabbit model for developmental and stem cell biology, the characterization of embryos at the molecular level is still poorly documented. We conducted a transcriptome analysis of rabbit pre-implantation embryos from E2.7 (morula stage) to E6.6 (early primitive streak stage) using bulk and single-cell RNA-sequencing, and single-cell Biomark qPCR. In parallel, we studied oxidative phosphorylation and glycolysis and analyzed active and repressive epigenetic modifications during blastocyst formation and expansion. We generated a transcriptomic, epigenetic, and metabolic map of the pluripotency continuum in rabbit preimplantation embryos and identified novel markers of naïve pluripotency that might be instrumental for deriving naïve pluripotent stem cell lines. Although the rabbit is evolutionarily closer to mice than to primates, we found that the transcriptome of rabbit epiblast cells shares common features with that of humans and non-human primates.

Genomic data allow an in-depth and renewed study of local adaptation. The red coral (Corallium rubrum, Cnidaria) is a highly genetically structured species and a promising model for the study of adaptive processes along an environmental gradient. Here, we used RAD-Sequencing in order to study the vertical genetic structure of this species and to search for signals of local adaptation to depth and thermal regime in the red coral. Previous studies have shown different thermotolerance levels according to depth in this species which could correspond to genetic or environmental differences. We designed a sampling scheme with six pairs of shallow vs deep populations distributed in three geographical regions as replicates. Our results showed significant differentiation among locations and among sites separated by around 20 m depth. The tests of association between genetics and environment allowed the identification of candidate loci under selection but with a potentially high rate of false positive. We discuss the methodological obstacles and biases encountered for the detection of selected loci in such a strongly genetically structured species. On this basis, we also discuss the significance of the candidate loci for local adaptation detected in each geographical region and the evolution of red coral populations along environmental gradients.

In most animals, sex determination occurs at conception, when sex chromosomes are segregated following Mendelian laws. However, in multiple reptiles and fishes, this genetic sex can be overridden by external factors after fertilization or birth. In some species, the genetic sex may also be governed by multiple genes, further limiting our understanding of sex determination in such species. We used the European sea bass (Dicentrarchus labrax) as a model and combined genomic (using a single nucleotide polymorphism chip) and transcriptomic (RNA-Sequencing) approaches to thoroughly depict this polygenic sex determination system and its interaction with temperature. We estimated genetic sex tendency (eGST), defined as the estimated genetic liability to become a given sex under a liability threshold model for sex determination, which accurately predicts the future phenotypic sex. We found evidence that energetic pathways, concerning the regulation of lipids and glucose, are involved in sex determination and could explain why females tend to exhibit higher energy levels and improved growth compared to males. Besides, early exposure to hightemperature up-regulated sox3, followed by sox9a in individuals with intermediate eGST, but not in individuals showing highly female-biased eGST, providing the most parsimonious explanation for temperature-induced masculinization. This gonadal state was maintained likely by DNA methylation and the up-regulation of several genes involved in histone modifications, including jmjd1c. Overall, we describe a sex determination system resulting from continuous genetic and environmental influences in an animal. Our results provide significant progress in our understanding of the mechanisms underlying temperature-induced masculinization in fish.

Chromatin organization is essential for appropriate interpretation of the genetic information. Here, we demonstrated that the chromatin-associated proteins HP1 are dispensable for hepatocytes survival but are essential within hepatocytes to prevent liver tumor development in mice with HP1β being pivotal in these functions. Yet, we found that the loss of HP1 per se is not sufficient to induce cell transformation but renders cells more resistant to specific stress such as the expression of oncogenes and thus in fine, more prone to cell transformation. Molecular characterization of HP1-Triple KO premalignant livers and BMEL cells revealed that HP1 are essential for the maintenance of heterochromatin organization and for the regulation of specific genes with most of them having well characterized functions in liver functions and homeostasis. We further showed that some specific retrotransposons get reactivated upon loss of HP1, correlating with overexpression of genes in their neighborhood. Interestingly, we found that, although HP1-dependent genes are characterized by enrichment H3K9me3, this mark does not require HP1 for its maintenance and is not sufficient to maintain gene repression in absence of HP1. Finally, we demonstrated that the loss of TRIM28 association with HP1 recapitulated several phenotypes induced by the loss of HP1 including the reactivation of some retrotransposons and the increased incidence of liver cancer development. Altogether, our findings indicate that HP1 proteins act as guardians of liver homeostasis to prevent tumor development by modulating multiple chromatin-associated events within both the heterochromatic and euchromatic compartments, partly through regulation of the corepressor TRIM28 activity.

Gastrointestinal stromal tumours (GISTs), the most common mesenchymal neoplasm of the gastrointestinal tract, result from deregulated proliferation of transformed KIT-positive interstitial cells of Cajal that share mesenchymal progenitors with smooth muscle cells. Despite the identification of selective KIT inhibitors, primary resistance and relapse remain a major concern. Moreover, most patients develop resistance partly through reactivation of KIT and its downstream signalling pathways. We previously identified the Limb Expression 1 (LIX1) gene as a unique marker of digestive mesenchyme immaturity. We also demonstrated that LIX1 regulates mesenchymal progenitor proliferation and differentiation by controlling the Hippo effector YAP1, which is constitutively activated in many sarcomas. Therefore, we wanted to determine LIX1 role in GIST development. We found that LIX1 is strongly up-regulated in GIST samples and this is associated with unfavourable prognosis. Moreover, LIX1 controls GIST cell proliferation in vitro and in vivo. Upon LIX1 inactivation in GIST cells, YAP1/TAZ activity is reduced, KIT (the GIST signature) is down-regulated, and cells acquire smooth muscle lineage features. Our data highlight LIX1 role in digestive mesenchyme-derived cell-fate decisions and identify this novel regulator as a target for drug design for GIST treatment by influencing its differentiation status.

Aging-dependent changes in tissue function are associated with the development of metabolic diseases. However, the molecular connections linking aging, obesity, and diabetes remain unclear. Lamin A, lamin C, and progerin, products of the Lmna gene, have antagonistic functions on energy metabolism and life span. Lamin C, albeit promoting obesity, increases life span, suggesting that this isoform is crucial for maintaining healthy conditions under metabolic stresses. Because β-cell loss during obesity or aging leads to diabetes, we investigated the contribution of lamin C to β-cell function in physiopathological conditions. We demonstrate that aged lamin C only-expressing mice (Lmna LCS/LCS ) become obese but remain glucose tolerant due to adaptive mechanisms including increased β-cell mass and insulin secretion. Triggering diabetes in young mice revealed that Lmna LCS/LCS animals normalize their fasting glycemia by both increasing insulin secretion and regenerating β-cells. Genome-wide analyses combined to functional analyses revealed an increase of mitochondrial biogenesis and global translational rate in Lmna LCS/LCS islets, two major processes involved in insulin secretion. Altogether, our results demonstrate for the first time that the sole expression of lamin C protects from glucose intolerance through a β-cell-adaptive transcriptional program during metabolic stresses, highlighting Lmna gene processing as a new therapeutic target for diabetes treatment.

Changes in ocean water temperature associated with global climate change are bound to enormously affect fish populations, with potential major economic consequences in the aquaculture and fisheries industries. A link between temperature fluctuations and changes in fish stress response is well established. In this study, we aimed to assess the effects of a short- (4 days) or a long-term (4 months) exposure to warm temperature in the stress physiology of European sea bass (Dicentrarchus labrax) larvae and juveniles. First, cortisol (i.e. the main stress hormone in fishes) analysis was used to confirm that a steady and short-term elevation of temperature acts as a physiological stressful event in these fish, and cortisol release is indeed above a metabolic increase linked to temperature. Moreover, our results verified that measurement of cortisol released into the water can be reliably employed as a non-invasive indicator of acute thermal stress in experimental conditions. Secondly, the different effects on the genetic cascade underlying the stress response between long-term low or high thermal treatments were evaluated at two larval development stages via candidate-gene and whole-transcriptome approaches. Interestingly, opposite expression for some key stress genes (nr3c1, nr3c2 and hsd11b2) were observed between developmental stages, highlighting the distinct adaptive mechanisms controlling the primary and secondary responses to a stressor. Surprising expression patterns for some understudied genes involved in the stress axis were also revealed, including crhr1, mc2r, mc5r, trh or trhr, which should be further explored. Finally, evaluation of cortisol content in scales was successfully used as a biomarker of chronic thermal stress, with 10x more cortisol in fish kept at 21 °C vs 16 °C after 4 months, supporting the gene expression results observed. The use of such a method as a proxy of long-term stress, unprecedented in the literature, holds a vast array of applications in further research, in particular, in the investigation of the impact of global warming on wild fish populations.

Additional file 1: Table S0. REPLICATES. Table S1. Genes differentially expressed between Full moon and new Moon at ZT0. Table S2. Genes differentially expressed between Full moon and new Moon at ZT3. Table S3. Genes differentially expressed between Full moon and new Moon at ZT6. Table S4. Genes differentially expressed between Full moon and new Moon at ZT9. Table S5. Genes differentially expressed between Full moon and new Moon at ZT12. Table S6. Genes differentially expressed between Full moon and new Moon at ZT15. Table S7. Genes differentially expressed between Full moon and new Moon at ZT18. Table S8. Genes differentially expressed between Full moon and new Moon at ZT21. Table S9. Genes differentially expressed between Full moon and new Moon at ZT24. Figure S10. Rhythmic genes in common between Coffee and Arabidopsis thaliana. Figure S11. Pairwise phase plots of the core clock genes. Figure S12. Number of genes per phase under new moon and the full moon. Table S13. Correlation matrix between Redox homeostasis genes and LHY. Table S14. Correlation matrix between genes involved in lipids pathways and LHY. Table S14. Primers used for genes studied by qRT-PCR. Figure S16. RT-QPCR of LHY gene during new moon. Figure S17. RT-QPCR time-course of core clock genes, POR1A and GWD1.