Gene Expression Changes and Community Turnover Differentially Shape the Global Ocean Metatranscriptome.

Ocean microbial communities strongly influence the biogeochemistry, food webs, and climate of our planet. Despite recent advances in understanding their taxonomic and genomic compositions, little is known about how their transcriptomes vary globally. Here, we present a dataset of 187 metatranscriptomes and 370 metagenomes from 126 globally distributed sampling stations and establish a resource of 47 million genes to study community-level transcriptomes across depth layers from pole-to-pole. We examine gene expression changes and community turnover as the underlying mechanisms shaping community transcriptomes along these axes of environmental variation and show how their individual contributions differ for multiple biogeochemically relevant processes. Furthermore, we find the relative contribution of gene expression changes to be significantly lower in polar than in non-polar waters and hypothesize that in polar regions, alterations in community activity in response to ocean warming will be driven more strongly by changes in organismal composition than by gene regulatory mechanisms. VIDEO ABSTRACT.

Global Trends in Marine Plankton Diversity across Kingdoms of Life.

The ocean is home to myriad small planktonic organisms that underpin the functioning of marine ecosystems. However, their spatial patterns of diversity and the underlying drivers remain poorly known, precluding projections of their responses to global changes. Here we investigate the latitudinal gradients and global predictors of plankton diversity across archaea, bacteria, eukaryotes, and major virus clades using both molecular and imaging data from Tara Oceans. We show a decline of diversity for most planktonic groups toward the poles, mainly driven by decreasing ocean temperatures. Projections into the future suggest that severe warming of the surface ocean by the end of the 21st century could lead to tropicalization of the diversity of most planktonic groups in temperate and polar regions. These changes may have multiple consequences for marine ecosystem functioning and services and are expected to be particularly significant in key areas for carbon sequestration, fisheries, and marine conservation. VIDEO ABSTRACT.

Interspecific introgression patterns reveal the origins of worldwide cultivated bananas in New Guinea.

Hybridizations between Musa species and subspecies, enabled by their transport via human migration, were proposed to have played an important role in banana domestication. We exploited sequencing data of 226 Musaceae accessions, including wild and cultivated accessions, to characterize the inter(sub)specific hybridization pattern that gave rise to cultivated bananas. We identified 11 genetic pools that contributed to cultivars, including two contributors of unknown origin. Informative alleles for each of these genetic pools were pinpointed and used to obtain genome ancestry mosaics of accessions. Diploid and triploid cultivars had genome mosaics involving three up to possibly seven contributors. The simplest mosaics were found for some diploid cultivars from New Guinea, combining three contributors, i.e., banksii and zebrina representing Musa acuminata subspecies and, more unexpectedly, the New Guinean species Musa schizocarpa. Breakpoints of M. schizocarpa introgressions were found to be conserved between New Guinea cultivars and the other analyzed diploid and triploid cultivars. This suggests that plants bearing these M. schizocarpa introgressions were transported from New Guinea and gave rise to currently cultivated bananas. Many cultivars showed contrasted mosaics with predominant ancestry from their geographical origin across Southeast Asia to New Guinea. This revealed that further diversification occurred in different Southeast Asian regions through hybridization with other Musa (sub)species, including two unknown ancestors that we propose to be M. acuminata ssp. halabanensis and a yet to be characterized M. acuminata subspecies. These results highlighted a dynamic crop formation process that was initiated in New Guinea, with subsequent diversification throughout Southeast Asia.

Dynamics of Gene Loss following Ancient Whole-Genome Duplication in the Cryptic Paramecium Complex.

Whole-genome duplications (WGDs) have shaped the gene repertoire of many eukaryotic lineages. The redundancy created by WGDs typically results in a phase of massive gene loss. However, some WGD-derived paralogs are maintained over long evolutionary periods, and the relative contributions of different selective pressures to their maintenance are still debated. Previous studies have revealed a history of three successive WGDs in the lineage of the ciliate Paramecium tetraurelia and two of its sister species from the Paramecium aurelia complex. Here, we report the genome sequence and analysis of 10 additional P. aurelia species and 1 additional out group, revealing aspects of post-WGD evolution in 13 species sharing a common ancestral WGD. Contrary to the morphological radiation of vertebrates that putatively followed two WGD events, members of the cryptic P. aurelia complex have remained morphologically indistinguishable after hundreds of millions of years. Biases in gene retention compatible with dosage constraints appear to play a major role opposing post-WGD gene loss across all 13 species. In addition, post-WGD gene loss has been slower in Paramecium than in other species having experienced genome duplication, suggesting that the selective pressures against post-WGD gene loss are especially strong in Paramecium. A near complete lack of recent single-gene duplications in Paramecium provides additional evidence for strong selective pressures against gene dosage changes. This exceptional data set of 13 species sharing an ancestral WGD and 2 closely related out group species will be a useful resource for future studies on Paramecium as a major model organism in the evolutionary cell biology.

Tracking population structure and phenology through time using ancient genomes from waterlogged white oak wood.

Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices.

Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes.

Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression.

Genome evolution across 1,011 Saccharomyces cerevisiae isolates.

Large-scale population genomic surveys are essential to explore the phenotypic diversity of natural populations. Here we report the whole-genome sequencing and phenotyping of 1,011 Saccharomyces cerevisiae isolates, which together provide an accurate evolutionary picture of the genomic variants that shape the species-wide phenotypic landscape of this yeast. Genomic analyses support a single 'out-of-China' origin for this species, followed by several independent domestication events. Although domesticated isolates exhibit high variation in ploidy, aneuploidy and genome content, genome evolution in wild isolates is mainly driven by the accumulation of single nucleotide polymorphisms. A common feature is the extensive loss of heterozygosity, which represents an essential source of inter-individual variation in this mainly asexual species. Most of the single nucleotide polymorphisms, including experimentally identified functional polymorphisms, are present at very low frequencies. The largest numbers of variants identified by genome-wide association are copy-number changes, which have a greater phenotypic effect than do single nucleotide polymorphisms. This resource will guide future population genomics and genotype-phenotype studies in this classic model system.

Oak stands along an elevation gradient have different molecular strategies for regulating bud phenology.

BACKGROUND: Global warming raises serious concerns about the persistence of species and populations locally adapted to their environment, simply because of the shift it produces in their adaptive landscape. For instance, the phenological cycle of tree species may be strongly affected by higher winter temperatures and late frost in spring. Given the variety of ecosystem services they provide, the question of forest tree adaptation has received increasing attention in the scientific community and catalyzed research efforts in ecology, evolutionary biology and functional genomics to study their adaptive capacity to respond to such perturbations. RESULTS: In the present study, we used an elevation gradient in the Pyrenees Mountains to explore the gene expression network underlying dormancy regulation in natural populations of sessile oak stands sampled along an elevation cline and potentially adapted to different climatic conditions mainly driven by temperature. By performing analyses of gene expression in terminal buds we identified genes displaying significant dormancy, elevation or dormancy-by-elevation interaction effects. Our Results highlighted that low- and high-altitude populations have evolved different molecular strategies for minimizing late frost damage and maximizing the growth period, thereby increasing potentially their respective fitness in these contrasting environmental conditions. More particularly, population from high elevation overexpressed genes involved in the inhibition of cell elongation and delaying flowering time while genes involved in cell division and flowering, enabling buds to flush earlier were identified in population from low elevation. CONCLUSION: Our study made it possible to identify key dormancy-by-elevation responsive genes revealing that the stands analyzed in this study have evolved distinct molecular strategies to adapt their bud phenology in response to temperature.

Gene expression and genetic divergence in oak species highlight adaptive genes to soil water constraints.

Drought and waterlogging impede tree growth and may even lead to tree death. Oaks, an emblematic group of tree species, have evolved a range of adaptations to cope with these constraints. The two most widely distributed European species, pedunculate (PO; Quercus robur L.) and sessile oak (SO; Quercus petraea Matt. Lieb), have overlapping ranges, but their respective distribution are highly constrained by local soil conditions. These contrasting ecological preferences between two closely related and frequently hybridizing species constitute a powerful model to explore the functional bases of the adaptive responses in oak. We exposed oak seedlings to waterlogging and drought, conditions typically encountered by the two species in their respective habitats, and studied changes in gene expression in roots using RNA-seq. We identified genes that change in expression between treatments differentially depending on species. These "species × environment"-responsive genes revealed adaptive molecular strategies involving adventitious and lateral root formation, aerenchyma formation in PO, and osmoregulation and ABA regulation in SO. With this experimental design, we also identified genes with different expression between species independently of water conditions imposed. Surprisingly, this category included genes with functions consistent with a role in intrinsic reproductive barriers. Finally, we compared our findings with those for a genome scan of species divergence and found that the expressional candidate genes included numerous highly differentiated genetic markers between the two species. By combining transcriptomic analysis, gene annotation, pathway analyses, as well as genome scan for genetic differentiation among species, we were able to highlight loci likely involved in adaptation of the two species to their respective ecological niches.

Shared pedigree relationships and transmission of unreduced gametes in cultivated banana.

BACKGROUND AND AIMS: Cultivated bananas resulted from inter(sub)specific hybridizations involving Musa species and subspecies (M. acuminata subspecies, M. schizocarpa, M. balbisiana) and the subsequent selection, centuries ago, of hybrids with parthenocarpic, seedless fruits. Cultivars have low fertility and are vegetatively propagated, forming groups of somaclones. Relatively few of them, mainly triploids, are grown on a large scale and characterization of their parental relationships may be useful for breeding strategies. Here we investigate parental relationships and gamete-type contributions among diploid and polyploid banana cultivars. METHODS: We used SNP genotyping data from whole-genome sequencing of 178 banana individuals, including 111 cultivars, 55 wild bananas and 12 synthetic F1 hybrids. We analysed the proportion of SNP sites in accordance with direct parentage with a global statistic and along chromosomes for selected individuals. KEY RESULTS: We characterized parentage relationships for 7 diploid cultivars, 11 triploid cultivars and 1 tetraploid cultivar. Results showed that both diploid and triploid cultivars could have contributed gametes to other banana cultivars. Diploids may have contributed 1x or 2x gametes and triploids 1x to 3x gametes. The Mchare diploid cultivar group, nowadays only found in East Africa, was found as parent of two diploid and eight triploid cultivars. In five of its identified triploid offspring, corresponding to main export or locally popular dessert bananas, Mchare contributed a 2x gamete with full genome restitution without recombination. Analyses of remaining haplotypes in these Mchare offspring suggested ancestral pedigree relationships between different interspecific banana cultivars. CONCLUSIONS: The current cultivated banana resulted from different pathways of formation, with implication of recombined or un-recombined unreduced gametes produced by diploid or triploid cultivars. Identification of dessert banana's parents and the types of gametes they contributed should support the design of breeding strategies.

Cyanorak v2.1: a scalable information system dedicated to the visualization and expert curation of marine and brackish picocyanobacteria genomes.

Cyanorak v2.1 (http://www.sb-roscoff.fr/cyanorak) is an information system dedicated to visualizing, comparing and curating the genomes of Prochlorococcus, Synechococcus and Cyanobium, the most abundant photosynthetic microorganisms on Earth. The database encompasses sequences from 97 genomes, covering most of the wide genetic diversity known so far within these groups, and which were split into 25,834 clusters of likely orthologous groups (CLOGs). The user interface gives access to genomic characteristics, accession numbers as well as an interactive map showing strain isolation sites. The main entry to the database is through search for a term (gene name, product, etc.), resulting in a list of CLOGs and individual genes. Each CLOG benefits from a rich functional annotation including EggNOG, EC/K numbers, GO terms, TIGR Roles, custom-designed Cyanorak Roles as well as several protein motif predictions. Cyanorak also displays a phyletic profile, indicating the genotype and pigment type for each CLOG, and a genome viewer (Jbrowse) to visualize additional data on each genome such as predicted operons, genomic islands or transcriptomic data, when available. This information system also includes a BLAST search tool, comparative genomic context as well as various data export options. Altogether, Cyanorak v2.1 constitutes an invaluable, scalable tool for comparative genomics of ecologically relevant marine microorganisms.

Transcriptome of the synganglion in the tick Ixodes ricinus and evolution of the cys-loop ligand-gated ion channel family in ticks.

BACKGROUND: Ticks represent a major health issue for humans and domesticated animals. Exploring the expression landscape of the tick's central nervous system (CNS), known as the synganglion, would be an important step in understanding tick physiology and in managing tick-borne diseases, but studies on that topic are still relatively scarce. Neuron-specific genes like the cys-loop ligand-gated ion channels (cys-loop LGICs, or cysLGICs) are important pharmacological targets of acaricides. To date their sequence have not been well catalogued for ticks, and their phylogeny has not been fully studied. RESULTS: We carried out the sequencing of transcriptomes of the I. ricinus synganglion, for adult ticks in different conditions (unfed males, unfed females, and partially-fed females). The de novo assembly of these transcriptomes allowed us to obtain a large collection of cys-loop LGICs sequences. A reference meta-transcriptome based on synganglion and whole body transcriptomes was then produced, showing high completeness and allowing differential expression analyses between synganglion and whole body. Many of the genes upregulated in the synganglion were associated with neurotransmission and/or localized in neurons or the synaptic membrane. As the first step of a functional study of cysLGICs, we cloned the predicted sequence of the resistance to dieldrin (RDL) subunit homolog, and functionally reconstituted the first GABA-gated receptor of Ixodes ricinus. A phylogenetic study was performed for the nicotinic acetylcholine receptors (nAChRs) and other cys-loop LGICs respectively, revealing tick-specific expansions of some types of receptors (especially for Histamine-like subunits and GluCls). CONCLUSIONS: We established a large catalogue of genes preferentially expressed in the tick CNS, including the cysLGICs. We discovered tick-specific gene family expansion of some types of cysLGIC receptors, and a case of intragenic duplication, suggesting a complex pattern of gene expression among different copies or different alternative transcripts of tick neuro-receptors.

Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp.

BACKGROUND: Dinoflagellates are aquatic protists particularly widespread in the oceans worldwide. Some are responsible for toxic blooms while others live in symbiotic relationships, either as mutualistic symbionts in corals or as parasites infecting other protists and animals. Dinoflagellates harbor atypically large genomes (~ 3 to 250 Gb), with gene organization and gene expression patterns very different from closely related apicomplexan parasites. Here we sequenced and analyzed the genomes of two early-diverging and co-occurring parasitic dinoflagellate Amoebophrya strains, to shed light on the emergence of such atypical genomic features, dinoflagellate evolution, and host specialization. RESULTS: We sequenced, assembled, and annotated high-quality genomes for two Amoebophrya strains (A25 and A120), using a combination of Illumina paired-end short-read and Oxford Nanopore Technology (ONT) MinION long-read sequencing approaches. We found a small number of transposable elements, along with short introns and intergenic regions, and a limited number of gene families, together contribute to the compactness of the Amoebophrya genomes, a feature potentially linked with parasitism. While the majority of Amoebophrya proteins (63.7% of A25 and 59.3% of A120) had no functional assignment, we found many orthologs shared with Dinophyceae. Our analyses revealed a strong tendency for genes encoded by unidirectional clusters and high levels of synteny conservation between the two genomes despite low interspecific protein sequence similarity, suggesting rapid protein evolution. Most strikingly, we identified a large portion of non-canonical introns, including repeated introns, displaying a broad variability of associated splicing motifs never observed among eukaryotes. Those introner elements appear to have the capacity to spread over their respective genomes in a manner similar to transposable elements. Finally, we confirmed the reduction of organelles observed in Amoebophrya spp., i.e., loss of the plastid, potential loss of a mitochondrial genome and functions. CONCLUSION: These results expand the range of atypical genome features found in basal dinoflagellates and raise questions regarding speciation and the evolutionary mechanisms at play while parastitism was selected for in this particular unicellular lineage.

Chromosome reciprocal translocations have accompanied subspecies evolution in bananas.

Chromosome rearrangements and the way that they impact genetic differentiation and speciation have long raised questions from evolutionary biologists. They are also a major concern for breeders because of their bearing on chromosome recombination. Banana is a major crop that derives from inter(sub)specific hybridizations between various once geographically isolated Musa species and subspecies. We sequenced 155 accessions, including banana cultivars and representatives of Musa diversity, and genotyped-by-sequencing 1059 individuals from 11 progenies. We precisely characterized six large reciprocal translocations and showed that they emerged in different (sub)species of Musa acuminata, the main contributor to currently cultivated bananas. Most diploid and triploid cultivars analyzed were structurally heterozygous for 1 to 4 M. acuminata translocations, highlighting their complex origin. We showed that all translocations induced a recombination reduction of variable intensity and extent depending on the translocations, involving only the breakpoint regions, a chromosome arm, or an entire chromosome. The translocated chromosomes were found preferentially transmitted in many cases. We explore and discuss the possible mechanisms involved in this preferential transmission and its impact on translocation colonization.

Ion channel profiling of the Lymnaea stagnalis ganglia via transcriptome analysis.

BACKGROUND: The pond snail Lymnaea stagnalis (L. stagnalis) has been widely used as a model organism in neurobiology, ecotoxicology, and parasitology due to the relative simplicity of its central nervous system (CNS). However, its usefulness is restricted by a limited availability of transcriptome data. While sequence information for the L. stagnalis CNS transcripts has been obtained from EST libraries and a de novo RNA-seq assembly, the quality of these assemblies is limited by a combination of low coverage of EST libraries, the fragmented nature of de novo assemblies, and lack of reference genome. RESULTS: In this study, taking advantage of the recent availability of a preliminary L. stagnalis genome, we generated an RNA-seq library from the adult L. stagnalis CNS, using a combination of genome-guided and de novo assembly programs to identify 17,832 protein-coding L. stagnalis transcripts. We combined our library with existing resources to produce a transcript set with greater sequence length, completeness, and diversity than previously available ones. Using our assembly and functional domain analysis, we profiled L. stagnalis CNS transcripts encoding ion channels and ionotropic receptors, which are key proteins for CNS function, and compared their sequences to other vertebrate and invertebrate model organisms. Interestingly, L. stagnalis transcripts encoding numerous putative Ca2+ channels showed the most sequence similarity to those of Mus musculus, Danio rerio, Xenopus tropicalis, Drosophila melanogaster, and Caenorhabditis elegans, suggesting that many calcium channel-related signaling pathways may be evolutionarily conserved. CONCLUSIONS: Our study provides the most thorough characterization to date of the L. stagnalis transcriptome and provides insights into differences between vertebrates and invertebrates in CNS transcript diversity, according to function and protein class. Furthermore, this study provides a complete characterization of the ion channels of Lymnaea stagnalis, opening new avenues for future research on fundamental neurobiological processes in this model system.

Characterization of Mollivirus kamchatka, the First Modern Representative of the Proposed Molliviridae Family of Giant Viruses.

Microbes trapped in permanently frozen paleosoils (permafrost) are the focus of increasing research in the context of global warming. Our previous investigations led to the discovery and reactivation of two Acanthamoeba-infecting giant viruses, Mollivirus sibericum and Pithovirus sibericum, from a 30,000-year old permafrost layer. While several modern pithovirus strains have since been isolated, no contemporary mollivirus relative was found. We now describe Mollivirus kamchatka, a close relative to M. sibericum, isolated from surface soil sampled on the bank of the Kronotsky River in Kamchatka, Russian Federation. This discovery confirms that molliviruses have not gone extinct and are at least present in a distant subarctic continental location. This modern isolate exhibits a nucleocytoplasmic replication cycle identical to that of M. sibericum Its spherical particle (0.6 µm in diameter) encloses a 648-kb GC-rich double-stranded DNA genome coding for 480 proteins, of which 61% are unique to these two molliviruses. The 461 homologous proteins are highly conserved (92% identical residues, on average), despite the presumed stasis of M. sibericum for the last 30,000 years. Selection pressure analyses show that most of these proteins contribute to virus fitness. The comparison of these first two molliviruses clarify their evolutionary relationship with the pandoraviruses, supporting their provisional classification in a distinct family, the Molliviridae, pending the eventual discovery of intermediary missing links better demonstrating their common ancestry.IMPORTANCE Virology has long been viewed through the prism of human, cattle, or plant diseases, leading to a largely incomplete picture of the viral world. The serendipitous discovery of the first giant virus visible under a light microscope (i.e., >0.3 µm in diameter), mimivirus, opened a new era of environmental virology, now incorporating protozoan-infecting viruses. Planet-wide isolation studies and metagenome analyses have shown the presence of giant viruses in most terrestrial and aquatic environments, including upper Pleistocene frozen soils. Those systematic surveys have led authors to propose several new distinct families, including the Mimiviridae, Marseilleviridae, Faustoviridae, Pandoraviridae, and Pithoviridae We now propose to introduce one additional family, the Molliviridae, following the description of M. kamchatka, the first modern relative of M. sibericum, previously isolated from 30,000-year-old arctic permafrost.

Sequence diversity and evolution of a group of iflaviruses associated with ticks.

We studied a group of tick-associated viruses with characteristics of members of the family Iflaviridae, a family of viruses frequently found in arthropods. Our aim was to gain insight into the evolutionary dynamics of this group of viruses, which may be linked to the biology of ticks. We explored assembled RNA-Seq data sets for different species of ticks. We identified members of five different iflavirus species, four of them novel, and discovered nine new genome sequences, including variants. Five variants represented a virus species associated with Ixodes ricinus. Unexpectedly, a sequence found in the Ixodes scapularis cell line ISE6 was nearly identical to the sequences of I. ricinus variants, suggesting a contamination of this cell line by I. ricinus material. Analysing patterns of substitutions between these variants, we detected a strong excess of synonymous mutations, suggesting evolution under strong positive selection. The phylogenies of the viruses and of their tick hosts were not congruent, suggesting recurrent host changes across tick genera during their evolution. Overall, our work constitutes a step in the understanding of the interactions between this family of viruses and ticks.

Adaptive introgression as a driver of local adaptation to climate in European white oaks.

Latitudinal and elevational gradients provide valuable experimental settings for studies of the potential impact of global warming on forest tree species. The availability of long-term phenological surveys in common garden experiments for traits associated with climate, such as bud flushing for sessile oaks (Quercus petraea), provide an ideal opportunity to investigate this impact. We sequenced 18 sessile oak populations and used available sequencing data for three other closely related European white oak species (Quercus pyrenaica, Quercus pubescens, and Quercus robur) to explore the evolutionary processes responsible for shaping the genetic variation across latitudinal and elevational gradients in extant sessile oaks. We used phenotypic surveys in common garden experiments and climatic data for the population of origin to perform genome-wide scans for population differentiation and genotype-environment and genotype-phenotype associations. The inferred historical relationships between Q. petraea populations suggest that interspecific gene flow occurred between Q. robur and Q. petraea populations from cooler or wetter areas. A genome-wide scan of differentiation between Q. petraea populations identified single nucleotide polymorphisms (SNPs) displaying strong interspecific relative divergence between these two species. These SNPs followed genetic clines along climatic or phenotypic gradients, providing further support for the likely contribution of introgression to the adaptive divergence of Q. petraea populations. Overall, the results indicate that outliers and associated SNPs are Q. robur ancestry-informative. We discuss the results of this study in the framework of the postglacial colonization scenario, in which introgression and diversifying selection have been proposed as essential drivers of Q. petraea microevolution.

Massive postglacial gene flow between European white oaks uncovered genes underlying species barriers.

Oaks are dominant forest tree species widely distributed across the Northern Hemisphere, where they constitute natural resources of economic, ecological, social and historical value. Hybridisation and adaptive introgression have long been thought to be major drivers of their ecological success. Therefore, the maintenance of species barriers remains a key question, given the extent of interspecific gene flow. In this study, we made use of the tremendous genetic variation among four European white oak species (31 million single nucleotide polymorphisms (SNPs)) to infer the evolutionary history of these species, study patterns of genetic differentiation and identify reproductive barriers. We first analysed the ecological and historical relationships among these species and inferred a long-term strict isolation followed by a recent and extensive postglacial contact using approximate Bayesian computation. Assuming this demographic scenario, we then performed backward simulations to generate the expected distributions of differentiation under neutrality to scan their genomes for reproductive barriers. We finally identified important intrinsic and ecological functions driving the reproductive isolation. We discussed the importance of identifying the genetic basis for the ecological preferences between these oak species and its implications for the renewal of European forests under global warming.

Genome-defence small RNAs exapted for epigenetic mating-type inheritance.

In the ciliate Paramecium, transposable elements and their single-copy remnants are deleted during the development of somatic macronuclei from germline micronuclei, at each sexual generation. Deletions are targeted by scnRNAs, small RNAs produced from the germ line during meiosis that first scan the maternal macronuclear genome to identify missing sequences, and then allow the zygotic macronucleus to reproduce the same deletions. Here we show that this process accounts for the maternal inheritance of mating types in Paramecium tetraurelia, a long-standing problem in epigenetics. Mating type E depends on expression of the transmembrane protein mtA, and the default type O is determined during development by scnRNA-dependent excision of the mtA promoter. In the sibling species Paramecium septaurelia, mating type O is determined by coding-sequence deletions in a different gene, mtB, which is specifically required for mtA expression. These independently evolved mechanisms suggest frequent exaptation of the scnRNA pathway to regulate cellular genes and mediate transgenerational epigenetic inheritance of essential phenotypic polymorphisms.