The Chalcidoidea bush of life: evolutionary history of a massive radiation of minute wasps.

Chalcidoidea are mostly parasitoid wasps that include as many as 500 000 estimated species. Capturing phylogenetic signal from such a massive radiation can be daunting. Chalcidoidea is an excellent example of a hyperdiverse group that has remained recalcitrant to phylogenetic resolution. We combined 1007 exons obtained with Anchored Hybrid Enrichment with 1048 ultra-conserved elements (UCEs) for 433 taxa including all extant families, >95% of all subfamilies, and 356 genera chosen to represent the vast diversity of the superfamily. Going back and forth between the molecular results and our collective knowledge of morphology and biology, we detected bias in the analyses that was driven by the saturation of nucleotide data. Our final results are based on a concatenated analysis of the least saturated exons and UCE datasets (2054 loci, 284 106 sites). Our analyses support an expected sister relationship with Mymarommatoidea. Seven previously recognized families were not monophyletic, so support for a new classification is discussed. Natural history in some cases would appear to be more informative than morphology, as illustrated by the elucidation of a clade of plant gall associates and a clade of taxa with planidial first-instar larvae. The phylogeny suggests a transition from smaller soft-bodied wasps to larger and more heavily sclerotized wasps, with egg parasitism as potentially ancestral for the entire superfamily. Deep divergences in Chalcidoidea coincide with an increase in insect families in the fossil record, and an early shift to phytophagy corresponds with the beginning of the "Angiosperm Terrestrial Revolution". Our dating analyses suggest a middle Jurassic origin of 174 Ma (167.3-180.5 Ma) and a crown age of 162.2 Ma (153.9-169.8 Ma) for Chalcidoidea. During the Cretaceous, Chalcidoidea may have undergone a rapid radiation in southern Gondwana with subsequent dispersals to the Northern Hemisphere. This scenario is discussed with regard to knowledge about the host taxa of chalcid wasps, their fossil record and Earth's palaeogeographic history.

Genetic independence between traits separated by metamorphosis is widespread but varies with biological function.

Why is metamorphosis so pervasive? Does it facilitate the independent (micro)evolution of quantitative traits in distinct life stages, similarly to how it enables some limbs and organs to develop at specific life stages? We tested this hypothesis by measuring the expression of 6400 genes in 41 Drosophila melanogaster inbred lines at larval and adult stages. Only 30% of the genes showed significant genetic correlations between larval and adult expression. By contrast, 46% of the traits showed some level of genetic independence between stages. Gene ontology terms enrichment revealed that across stages correlated traits were often involved in proteins synthesis, insecticide resistance and innate immunity, while a vast number of genes expression traits associated with energy metabolism were independent between life stages. We compared our results to a similar case: genetic constraints between males and females in gonochoric species (i.e. sexual antagonism). We expected selection for the separation between males and females to be higher than between juvenile and adult functions, as gonochorism is a more common strategy in the animal kingdom than metamorphosis. Surprisingly, we found that inter-stage constraints were lower than inter-sexual genetic constraints. Overall, our results show that metamorphosis enables a large part of the transcriptome to evolve independently at different life stages.

Molecular complexity and gene expression controlling cell turnover during a digestive cycle of carnivorous sponge Lycopodina hypogea.

Lycopodina hypogea is a carnivorous sponge that tolerates laboratory husbandry very well. During a digestion cycle, performed without any digestive cavity, this species undergoes spectacular morphological changes leading to a total regression of long filaments that ensure the capture of prey and their reformation at the end of the cycle. This phenomenon is a unique opportunity to analyze the molecular and cellular determinants that ensure digestion in the sister group of all other metazoans. Using differential transcriptomic analysis coupled with cell biology studies of proliferation, differentiation, and programmed cell deaths (i.e., autophagy and the destructive/constructive function of apoptosis), we demonstrate that the molecular and cellular actors that ensure digestive homeostasis in a sister group of all remaining animals are similar in variety and complexity to those controlling tissue homeostasis in higher vertebrates. During a digestion cycle, most of these actors are finely tuned in a coordinated manner. Our data benefits from complementary approaches coupling in silico and cell biology studies and demonstrate that the nutritive function is provided by the coordination of molecular network that impacts the cells turnover in the entire organism.

Ultra-Conserved Elements and morphology reciprocally illuminate conflicting phylogenetic hypotheses in Chalcididae (Hymenoptera, Chalcidoidea).

Recent technical advances combined with novel computational approaches have promised the acceleration of our understanding of the tree of life. However, when it comes to hyperdiverse and poorly known groups of invertebrates, studies are still scarce. As published phylogenies will be rarely challenged by future taxonomists, careful attention must be paid to potential analytical bias. We present the first molecular phylogenetic hypothesis for the family Chalcididae, a group of parasitoid wasps, with a representative sampling (144 ingroups and seven outgroups) that covers all described subfamilies and tribes, and 82% of the known genera. Analyses of 538 Ultra-Conserved Elements (UCEs) with supermatrix (RAxML and IQTREE) and gene tree reconciliation approaches (ASTRAL, ASTRID) resulted in highly supported topologies in overall agreement with morphology but reveal conflicting topologies for some of the deepest nodes. To resolve these conflicts, we explored the phylogenetic tree space with clustering and gene genealogy interrogation methods, analyzed marker and taxon properties that could bias inferences and performed a thorough morphological analysis (130 characters encoded for 40 taxa representative of the diversity). This joint analysis reveals that UCEs enable attainment of resolution between ancestry and convergent/divergent evolution when morphology is not informative enough, but also shows that a systematic exploration of bias with different analytical methods and a careful analysis of morphological features is required to prevent publication of artifactual results. We highlight a GC content bias for maximum-likelihood approaches, an artifactual mid-point rooting of the ASTRAL tree and a deleterious effect of high percentage of missing data (>85% missing UCEs) on gene tree reconciliation methods. Based on the results we propose a new classification of the family into eight subfamilies and ten tribes that lay the foundation for future studies on the evolutionary history of Chalcididae.

Argonaute proteins regulate HIV-1 multiply spliced RNA and viral production in a Dicer independent manner.

Argonaute (Ago) proteins associate with microRNAs (miRNAs) to form the core of the RNA-induced silencing complex (RISC) that mediates post-transcriptional gene silencing of target mRNAs. As key players in anti-viral defense, Ago proteins are thought to have the ability to interact with human immunodeficiency virus type 1 (HIV-1) RNA. However, the role of this interaction in regulating HIV-1 replication has been debated. Here, we used high throughput sequencing of RNA isolated by cross-linking immunoprecipitation (HITS-CLIP) to explore the interaction between Ago2 and HIV-1 RNA in infected cells. By only considering reads of 50 nucleotides length in our analysis, we identified more than 30 distinct binding sites for Ago2 along the viral RNA genome. Using reporter assays, we found four binding sites, located near splice donor sites, capable of repressing Luciferase gene expression in an Ago-dependent manner. Furthermore, inhibition of Ago1 and Ago2 levels in cells expressing HIV-1 led to an increase of viral multiply spliced transcripts and to a strong reduction in the extracellular CAp24 level. Depletion of Dicer did not affect these activities. Our results highlight a new role of Ago proteins in the control of multiply spliced HIV-1 transcript levels and viral production, independently of the miRNA pathway.

In Vitro Corticogenesis from Embryonic Stem Cells Recapitulates the In Vivo Epigenetic Control of Imprinted Gene Expression.

In vitro corticogenesis from embryonic stem cells (ESCs) is an attractive model of cortical development and a promising tool for cortical therapy. It is unknown to which extent epigenetic mechanisms crucial for cortex development and function, such as parental genomic imprinting, are recapitulated by in vitro corticogenesis. Here, using genome-wide transcriptomic and methylation analyses on hybrid mouse tissues and cells, we find a high concordance of imprinting status between in vivo and ESC-derived cortices. Notably, in vitro corticogenesis strictly reproduced the in vivo parent-of-origin-dependent expression of 41 imprinted genes (IGs), including Mest and Cdkn1c known to control corticogenesis. Parent-of-origin-dependent DNA methylation was also conserved at 14 of 18 imprinted differentially methylated regions. The least concordant imprinted locus was Gpr1-Zdbf2, where the aberrant bi-allelic expression of Zdbf2 and Adam23 was concomitant with a gain of methylation on the maternal allele in vitro. Combined, our data argue for a broad conservation of the epigenetic mechanisms at imprinted loci in cortical cells derived from ESCs. We propose that in vitro corticogenesis helps to define the still poorly understood mechanisms that regulate imprinting in the brain and the roles of IGs in cortical development.

Prostaglandin D2 acts through the Dp2 receptor to influence male germ cell differentiation in the foetal mouse testis.

Through intercellular signalling, the somatic compartment of the foetal testis is able to program primordial germ cells to undergo spermatogenesis. Fibroblast growth factor 9 and several members of the transforming growth factor ß superfamily are involved in this process in the foetal testis, counteracting the induction of meiosis by retinoic acid and activating germinal mitotic arrest. Here, using in vitro and in vivo approaches, we show that prostaglandin D2 (PGD2), which is produced through both L-Pgds and H-Pgds enzymatic activities in the somatic and germ cell compartments of the foetal testis, plays a role in mitotic arrest in male germ cells by activating the expression and nuclear localization of the CDK inhibitor p21(Cip1) and by repressing pluripotency markers. We show that PGD2 acts through its Dp2 receptor, at least in part through direct effects in germ cells, and contributes to the proper differentiation of male germ cells through the upregulation of the master gene Nanos2. Our data identify PGD2 signalling as an early pathway that acts in both paracrine and autocrine manners, and contributes to the differentiation of germ cells in the foetal testis.

A duplicated NUCLEOLIN gene with antagonistic activity is required for chromatin organization of silent 45S rDNA in Arabidopsis.

In plants as well as in animals, hundreds to thousands of 45S rRNA gene copies localize in Nucleolus Organizer Regions (NORs), and the activation or repression of specific sets of rDNA depends on epigenetic mechanisms. Previously, we reported that the Arabidopsis thaliana nucleolin protein NUC1, an abundant and evolutionarily conserved nucleolar protein in eukaryotic organisms, is required for maintaining DNA methylation levels and for controlling the expression of specific rDNA variants in Arabidopsis. Interestingly, in contrast with animal or yeast cells, plants contain a second nucleolin gene. Here, we report that Arabidopsis NUC1 and NUC2 nucleolin genes are both required for plant growth and survival and that NUC2 disruption represses flowering. However, these genes seem to be functionally antagonistic. In contrast with NUC1, disruption of NUC2 induces CG hypermethylation of rDNA and NOR association with the nucleolus. Moreover, NUC2 loss of function triggers major changes in rDNA spatial organization, expression, and transgenerational stability. Our analyses indicate that silencing of specific rRNA genes is mostly determined by the active or repressed state of the NORs and that nucleolin proteins play a key role in the developmental control of this process.

Empirical assessment of RAD sequencing for interspecific phylogeny.

Next-generation sequencing opened up new possibilities in phylogenetics; however, choosing an appropriate method of sample preparation remains challenging. Here, we demonstrate that restriction-site-associated DNA sequencing (RAD-seq) generates useful data for phylogenomics. Analysis of our RAD library using current bioinformatic and phylogenetic tools produced 400× more sites than our Sanger approach (2,262,825 nt/species), fully resolving relationships between 18 species of ground beetles (divergences up to 17 My). This suggests that RAD-seq is promising to infer phylogeny of eukaryotic species, though potential biases need to be evaluated and new methodologies developed to take full advantage of such data.

Who is the puppet master? Replication of a parasitic wasp-associated virus correlates with host behaviour manipulation.

Many parasites modify their host behaviour to improve their own transmission and survival, but the proximate mechanisms remain poorly understood. An original model consists of the parasitoid Dinocampus coccinellae and its coccinellid host, Coleomegilla maculata; during the behaviour manipulation, the parasitoid is not in contact with its host anymore. We report herein the discovery and characterization of a new RNA virus of the parasitoid (D. coccinellae paralysis virus, DcPV). Using a combination of RT-qPCR and transmission electron microscopy, we demonstrate that DcPV is stored in the oviduct of parasitoid females, replicates in parasitoid larvae and is transmitted to the host during larval development. Next, DcPV replication in the host's nervous tissue induces a severe neuropathy and antiviral immune response that correlate with the paralytic symptoms characterizing the behaviour manipulation. Remarkably, virus clearance correlates with recovery of normal coccinellid behaviour. These results provide evidence that changes in ladybeetle behaviour most likely result from DcPV replication in the cerebral ganglia rather than by manipulation by the parasitoid. This offers stimulating prospects for research on parasitic manipulation by suggesting for the first time that behaviour manipulation could be symbiont-mediated.

De novo assembly, characterization and functional annotation of Senegalese sole (Solea senegalensis) and common sole (Solea solea) transcriptomes: integration in a database and design of a microarray.

BACKGROUND: Senegalese sole (Solea senegalensis) and common sole (S. solea) are two economically and evolutionary important flatfish species both in fisheries and aquaculture. Although some genomic resources and tools were recently described in these species, further sequencing efforts are required to establish a complete transcriptome, and to identify new molecular markers. Moreover, the comparative analysis of transcriptomes will be useful to understand flatfish evolution. RESULTS: A comprehensive characterization of the transcriptome for each species was carried out using a large set of Illumina data (more than 1,800 millions reads for each sole species) and 454 reads (more than 5 millions reads only in S. senegalensis), providing coverages ranging from 1,384x to 2,543x. After a de novo assembly, 45,063 and 38,402 different transcripts were obtained, comprising 18,738 and 22,683 full-length cDNAs in S. senegalensis and S. solea, respectively. A reference transcriptome with the longest unique transcripts and putative non-redundant new transcripts was established for each species. A subset of 11,953 reference transcripts was qualified as highly reliable orthologs (>97% identity) between both species. A small subset of putative species-specific, lineage-specific and flatfish-specific transcripts were also identified. Furthermore, transcriptome data permitted the identification of single nucleotide polymorphisms and simple-sequence repeats confirmed by FISH to be used in further genetic and expression studies. Moreover, evidences on the retention of crystallins crybb1, crybb1-like and crybb3 in the two species of soles are also presented. Transcriptome information was applied to the design of a microarray tool in S. senegalensis that was successfully tested and validated by qPCR. Finally, transcriptomic data were hosted and structured at SoleaDB. CONCLUSIONS: Transcriptomes and molecular markers identified in this study represent a valuable source for future genomic studies in these economically important species. Orthology analysis provided new clues regarding sole genome evolution indicating a divergent evolution of crystallins in flatfish. The design of a microarray and establishment of a reference transcriptome will be useful for large-scale gene expression studies. Moreover, the integration of transcriptomic data in the SoleaDB will facilitate the management of genomic information in these important species.

Transcriptome population genomics reveals severe bottleneck and domestication cost in the African rice (Oryza glaberrima).

The African cultivated rice (Oryza glaberrima) was domesticated in West Africa 3000 years ago. Although less cultivated than the Asian rice (O. sativa), O. glaberrima landraces often display interesting adaptation to rustic environment (e.g. drought). Here, using RNA-seq technology, we were able to compare more than 12,000 transcripts between 9 O. glaberrima, 10 wild O. barthii and one O. meridionalis individuals. With a synonymous nucleotide diversity πs = 0.0006 per site, O. glaberrima appears as the least genetically diverse crop grass ever documented. Using approximate Bayesian computation, we estimated that O. glaberrima experienced a severe bottleneck during domestication. This demographic scenario almost fully accounts for the pattern of genetic diversity across O. glaberrima genome as we detected very few outliers regions where positive selection may have further impacted genetic diversity. Moreover, the large excess of derived nonsynonymous substitution that we detected suggests that the O. glaberrima population suffered from the 'cost of domestication'. In addition, we used this genome-scale data set to demonstrate that (i) O. barthii genetic diversity is positively correlated with recombination rate and negatively with gene density, (ii) expression level is negatively correlated with evolutionary constraint, and (iii) one region on chromosome 5 (position 4-6 Mb) exhibits a clear signature of introgression with a yet unidentified Oryza species. This work represents the first genome-wide survey of the African rice genetic diversity and paves the way for further comparison between the African and the Asian rice, notably regarding the genetics underlying domestication traits.

Complete mitogenome data from a European specimen of Ostrinia scapulalis (Walker, 1859) (Lepidoptera, Pyraloidea, Crambidae, Pyraustinae).

We present an assembly and annotation of the mitogenome of a European specimen of the Adzuki bean borer, Ostrinia scapulalis (Walker, 1859). The data were obtained by combining WGS data issue of a de novo and a previously published sequence library (Gschloessl et al., 2018). We also provide the phylogenetic positioning of the mitogenome within the Ostrinia genus, the Crambidae family and with more distant Lepidoptera species.

Optimized DNA extraction and library preparation for minute arthropods: Application to target enrichment in chalcid wasps used for biocontrol.

Target enrichment is increasingly used for genotyping of plant and animal species or to better understand the evolutionary history of important lineages through the inference of statistically robust phylogenies. Limitations to routine target enrichment are both the complexity of current protocols and low input DNA quantity. Thus, working with tiny organisms such as microarthropods can be challenging. Here, we propose easy to set up optimizations for DNA extraction and library preparation prior to target enrichment. Prepared libraries were used to capture 1,432 ultraconserved elements (UCEs) from microhymenoptera (Chalcidoidea), which are among the tiniest insects on Earth and the most commercialized worldwide for biological control purposes. Results show no correlation between input DNA quantities (1.8-250 ng, 0.4 ng with an extra whole genome amplification step) and the number of sequenced UCEs on an Illumina MiSeq. Phylogenetic inferences highlight the potential of UCEs to solve relationships within the families of chalcid wasps, which has not been achieved so far. The protocol (library preparation + target enrichment) allows processing 96 specimens in five working days, by a single person, without requiring the use of expensive robotic molecular biology platforms, which could help to generalize the use of target enrichment for minute specimens.

New WGS data and annotation of the heterosomal vs. autosomal localization of Ostrinia scapulalis (Lepidoptera, Crambidae) nuclear genomic scaffolds.

Here, we introduce new whole-genome shotgun sequencing and annotation data describing the autosomal vs. Z-heterosomal localization of nuclear genomic scaffolds of the moth species Ostrinia scapulalis. Four WGS libraries (corresponding to 2 males and 2 females) were sequenced with an Illumina HiSeq2500 sequencing technology, and the so-called 'AD-ratio' method was applied to distinguish between autosomal and Z-heterosomal scaffolds based on sequencing depth comparisons between homogametic (male) and heterogametic (female) libraries. A total of 25,760 scaffolds (corresponding to 341.69 Mb) were labelled as autosomal and 1273 scaffolds (15.29 Mb) were labelled as Z-heterosomal, totaling about 357 Mb. Besides, 4874 scaffolds (29.07 Mb) remain ambiguous because of a lack of AD-ratio reproducibility between the two replicates. The annotation method was evaluated a posteriori, by comparing depth-based annotation with the exact localization of known genes. Raw genomic data have been deposited and made accessible via the EMBL ENA BioProject id PRJEB26557. Comprehensive annotation is made accessible via the LepidoDB database (http://bipaa.genouest.org/sp/ostrinia_scapulalis/download/genome/v1.2/).

Transcriptome Profiles of Nod Factor-independent Symbiosis in the Tropical Legume Aeschynomene evenia.

Nod factors (NF) were assumed to be indispensable for the establishment of a rhizobium-legume symbiosis until the discovery that certain Bradyrhizobium strains interacting with certain Aeschynomene species lack the canonical nodABC genes required for their synthesis. So far, the molecular dialogue between Aeschynomene and its symbionts remains an open question. Here we report a time course transcriptional analysis of Aeschynomene evenia in response to inoculation with Bradyrhizobium ORS278. The NF-independent symbiotic process was monitored at five time points between bacterial infection and nodule maturity. The five time points correspond to three specific events, root infection by crack entry, nodule organogenesis, and the establishment of the nitrogen fixing process. During the third stage, about 80 NCR-like genes and eight symbiotic genes known to be involved in signaling, bacterial infection or nodulation regulation were highly expressed. Comparative gene expression analyses at the five time points also enabled the selection of genes with an expression profile that makes them promising markers to monitor early plant responses to bacteria. Such markers could be used in bioassays to identify the nature of the bacterial signal(s). Our data represent valuable resources for investigation of this Nod factor-independent symbiosis.

Using insects to detect, monitor and predict the distribution of Xylella fastidiosa: a case study in Corsica.

We sampled ca 2500 specimens of Philaenus spumarius (Hemiptera: Aphrophoridae) throughout Corsica without a priori knowledge on the presence of symptoms on plants. We screened 448 specimens for the presence of Xylella fastidiosa (Xf) using qPCR and a custom nested PCR. qPCR appeared versatile and under-estimated the prevalence of Xf. Nested PCR showed that Xf was present in all populations. Molecular results were validated by prediction on the distribution of Xf made from tests conducted on plants, which shows the pertinence of using vectors in risk assessment studies. Xf was detected in tenerals and adults. Thus, P. spumarius could acquire Xf from its host plant, mostly Cistus monspeliensis in Corsica, which may act as reservoir for the next season. This contrasts with other observations and suggests that management strategies may have to be adapted on a case-by-case basis. At least two genetic entities and several variants of Xf not yet identified on plants were present in the insects, which suggests ancient introductions of Xf and a probable underestimation of the current diversity of the strains present in Corsica. Interestingly 6% of the specimens carried two subspecies of Xf. Studies are required to better characterize the strains present in Corsica and to determine how the disease was introduced, spread and why no sign of a potential epidemic was detected earlier. This study shows that, when sensitive enough methods are implemented, spittlebugs (and more specifically P. spumarius for which species distribution modelling shows it could be a good sentinel for Europe) can be used to predict and better assess the exact distribution of Xf. Furthermore, Xf multiply only in their foregut and does not become circulative, which facilitates its detection.

New Insights into the Anti-pathogenic Potential of Lactococcus garvieae against Staphylococcus aureus Based on RNA Sequencing Profiling.

The bio-preservation potential of Lactococcus garvieae lies in its capacity to inhibit the growth of staphylococci, especially Staphylococcus aureus, in dairy products and in vitro. In vitro, inhibition is modulated by the level of aeration, owing to hydrogen peroxide (H2O2) production by L. garvieae under aeration. The S. aureus response to this inhibition has already been studied. However, the molecular mechanisms of L. garvieae underlying the antagonism against S. aureus have never been explored. This study provides evidence of the presence of another extracellular inhibition effector in vitro. This effector was neither a protein, nor a lipid, nor a polysaccharide, nor related to an L-threonine deficiency. To better understand the H2O2-related inhibition mechanism at the transcriptome level and to identify other mechanisms potentially involved, we used RNA sequencing to determine the transcriptome response of L. garvieae to different aeration levels and to the presence or absence of S. aureus. The L. garvieae transcriptome differed radically between different aeration levels mainly in biological processes related to fundamental functions and nutritional adaptation. The transcriptomic response of L. garvieae to aeration level differed according to the presence or absence of S. aureus. The higher concentration of H2O2 with high aeration was not associated with a higher expression of L. garvieae H2O2-synthesis genes (pox, sodA, and spxA1) but rather with a repression of L. garvieae H2O2-degradation genes (trxB1, ahpC, ahpF, and gpx). We showed that L. garvieae displayed an original, previously undiscovered, H2O2 production regulation mechanism among bacteria. In addition to the key factor H2O2, the involvement of another extracellular effector in the antagonism against S. aureus was shown. Future studies should explore the relation between H2O2-metabolism, H2O2-producing LAB and the pathogen they inhibit. The nature of the other extracellular effector should also be determined.

Two genomes of highly polyphagous lepidopteran pests (Spodoptera frugiperda, Noctuidae) with different host-plant ranges.

Emergence of polyphagous herbivorous insects entails significant adaptation to recognize, detoxify and digest a variety of host-plants. Despite of its biological and practical importance - since insects eat 20% of crops - no exhaustive analysis of gene repertoires required for adaptations in generalist insect herbivores has previously been performed. The noctuid moth Spodoptera frugiperda ranks as one of the world's worst agricultural pests. This insect is polyphagous while the majority of other lepidopteran herbivores are specialist. It consists of two morphologically indistinguishable strains ("C" and "R") that have different host plant ranges. To describe the evolutionary mechanisms that both enable the emergence of polyphagous herbivory and lead to the shift in the host preference, we analyzed whole genome sequences from laboratory and natural populations of both strains. We observed huge expansions of genes associated with chemosensation and detoxification compared with specialist Lepidoptera. These expansions are largely due to tandem duplication, a possible adaptation mechanism enabling polyphagy. Individuals from natural C and R populations show significant genomic differentiation. We found signatures of positive selection in genes involved in chemoreception, detoxification and digestion, and copy number variation in the two latter gene families, suggesting an adaptive role for structural variation.

Differential root transcriptomics in a polyploid non-model crop: the importance of respiration during osmotic stress.

To explore the transcriptomic global response to osmotic stress in roots, 18 mRNA-seq libraries were generated from three triploid banana genotypes grown under mild osmotic stress (5% PEG) and control conditions. Illumina sequencing produced 568 million high quality reads, of which 70-84% were mapped to the banana diploid reference genome. Using different uni- and multivariate statistics, 92 genes were commonly identified as differentially expressed in the three genotypes. Using our in house workflow to analyze GO enriched and underlying biochemical pathways, we present the general processes affected by mild osmotic stress in the root and focus subsequently on the most significantly overrepresented classes associated with: respiration, glycolysis and fermentation. We hypothesize that in fast growing and oxygen demanding tissues, mild osmotic stress leads to a lower energy level, which induces a metabolic shift towards (i) a higher oxidative respiration, (ii) alternative respiration and (iii) fermentation. To confirm the mRNA-seq results, a subset of twenty up-regulated transcripts were further analysed by RT-qPCR in an independent experiment at three different time points. The identification and annotation of this set of genes provides a valuable resource to understand the importance of energy sensing during mild osmotic stress.