CAM evolution is associated with gene family expansion in an explosive bromeliad radiation.

The subgenus Tillandsia (Bromeliaceae) belongs to one of the fastest radiating clades in the plant kingdom and is characterised by the repeated evolution of Crassulacean acid metabolism (CAM). Despite its complex genetic basis, this water-conserving trait has evolved independently across many plant families and is regarded as a key innovation trait and driver of ecological diversification in Bromeliaceae. By producing high-quality genome assemblies of a Tillandsia species pair displaying divergent photosynthetic phenotypes, and combining genome-wide investigations of synteny, transposable element (TE) dynamics, sequence evolution, gene family evolution and temporal differential expression, we were able to pinpoint the genomic drivers of CAM evolution in Tillandsia. Several large-scale rearrangements associated with karyotype changes between the two genomes and a highly dynamic TE landscape shaped the genomes of Tillandsia. However, our analyses show that rewiring of photosynthetic metabolism is mainly obtained through regulatory evolution rather than coding sequence evolution, as CAM-related genes are differentially expressed across a 24-hour cycle between the two species but are not candidates of positive selection. Gene orthology analyses reveal that CAM-related gene families manifesting differential expression underwent accelerated gene family expansion in the constitutive CAM species, further supporting the view of gene family evolution as a driver of CAM evolution.

Low-frequency somatic mutations are heritable in tropical trees Dicorynia guianensis and Sextonia rubra.

Somatic mutations potentially play a role in plant evolution, but common expectations pertaining to plant somatic mutations remain insufficiently tested. Unlike in most animals, the plant germline is assumed to be set aside late in development, leading to the expectation that plants accumulate somatic mutations along growth. Therefore, several predictions were made on the fate of somatic mutations: mutations have generally low frequency in plant tissues; mutations at high frequency have a higher chance of intergenerational transmission; branching topology of the tree dictates mutation distribution; and exposure to UV (ultraviolet) radiation increases mutagenesis. To provide insights into mutation accumulation and transmission in plants, we produced two high-quality reference genomes and a unique dataset of 60 high-coverage whole-genome sequences of two tropical tree species, Dicorynia guianensis (Fabaceae) and Sextonia rubra (Lauraceae). We identified 15,066 de novo somatic mutations in D. guianensis and 3,208 in S. rubra, surprisingly almost all found at low frequency. We demonstrate that 1) low-frequency mutations can be transmitted to the next generation; 2) mutation phylogenies deviate from the branching topology of the tree; and 3) mutation rates and mutation spectra are not demonstrably affected by differences in UV exposure. Altogether, our results suggest far more complex links between plant growth, aging, UV exposure, and mutation rates than commonly thought.

Allele surfing causes maladaptation in a Pacific salmon of conservation concern.

How various factors, including demography, recombination or genome duplication, may impact the efficacy of natural selection and the burden of deleterious mutations, is a central question in evolutionary biology and genetics. In this study, we show that key evolutionary processes, including variations in i) effective population size (Ne) ii) recombination rates and iii) chromosome inheritance, have influenced the genetic load and efficacy of selection in Coho salmon (Oncorhynchus kisutch), a widely distributed salmonid species on the west coast of North America. Using whole genome resequencing data from 14 populations at different migratory distances from their southern glacial refugium, we found evidence supporting gene surfing, wherein reduced Ne at the postglacial recolonization front, leads to a decrease in the efficacy of selection and a surf of deleterious alleles in the northernmost populations. Furthermore, our results indicate that recombination rates play a prime role in shaping the load along the genome. Additionally, we identified variation in polyploidy as a contributing factor to within-genome variation of the load. Overall, our results align remarkably well with expectations under the nearly neutral theory of molecular evolution. We discuss the fundamental and applied implications of these findings for evolutionary and conservation genomics.

Demographic responses of oceanic island birds to local and regional ecological disruptions revealed by whole-genome sequencing.

Disentangling the effects of ecological disruptions operating at different spatial and temporal scales in shaping past species' demography is particularly important in the current context of rapid environmental changes driven by both local and regional factors. We argue that volcanic oceanic islands provide useful settings to study the influence of past ecological disruptions operating at local and regional scales on population demographic histories. We investigate potential drivers of past population dynamics for three closely related species of passerine birds from two volcanic oceanic islands, Reunion and Mauritius (Mascarene archipelago), with distinct volcanic history. Using ABC and PSMC inferences from complete genomes, we reconstructed the demographic history of the Reunion Grey White-eye (Zosterops borbonicus (Pennant, 1781)), the Reunion Olive White-eye (Z. olivaceus (Linnaeus, 1766)) and the Mauritius Grey White-eye (Z. mauritianus (Gmelin, 1789)) and searched for possible causes underlying similarities or differences between species living on the same or different islands. Both demographic inferences strongly support ancient and long-term expansions in all species. They also reveal different trajectories between species inhabiting different islands, but consistent demographic trajectories in species or populations from the same island. Species from Reunion appear to have experienced synchronous reductions in population size during the Last Glacial Maximum, a trend not seen in the Mauritian species. Overall, this study suggests that local events may have played a role in shaping population trajectories of these island species. It also highlights the potential of our conceptual framework to disentangle the effects of local and regional drivers on past species' demography and long-term population processes.

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.

Correction: Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific Salmon.

[This corrects the article DOI: 10.1371/journal.pgen.1008348.].

Drivers of genomic landscapes of differentiation across a Populus divergence gradient.

Speciation, the continuous process by which new species form, is often investigated by looking at the variation of nucleotide diversity and differentiation across the genome (hereafter genomic landscapes). A key challenge lies in how to determine the main evolutionary forces at play shaping these patterns. One promising strategy, albeit little used to date, is to comparatively investigate these genomic landscapes as progression through time by using a series of species pairs along a divergence gradient. Here, we resequenced 201 whole-genomes from eight closely related Populus species, with pairs of species at different stages along the divergence gradient to learn more about speciation processes. Using population structure and ancestry analyses, we document extensive introgression between some species pairs, especially those with parapatric distributions. We further investigate genomic landscapes, focusing on within-species (i.e. nucleotide diversity and recombination rate) and among-species (i.e. relative and absolute divergence) summary statistics of diversity and divergence. We observe relatively conserved patterns of genomic divergence across species pairs. Independent of the stage across the divergence gradient, we find support for signatures of linked selection (i.e. the interaction between natural selection and genetic linkage) in shaping these genomic landscapes, along with gene flow and standing genetic variation. We highlight the importance of investigating genomic patterns on multiple species across a divergence gradient and discuss prospects to better understand the evolutionary forces shaping the genomic landscapes of diversity and differentiation.

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.

Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific Salmon.

A thorough reconstruction of historical processes is essential for a comprehensive understanding of the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded recently, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored.

Within-island diversification in a passerine bird.

The presence of congeneric taxa on the same island suggests the possibility of in situ divergence, but can also result from multiple colonizations of previously diverged lineages. Here, using genome-wide data from a large population sample, we test the hypothesis that intra-island divergence explains the occurrence of four geographical forms meeting at hybrid zones in the Reunion grey white-eye (Zosterops borbonicus), a species complex endemic to the small volcanic island of Reunion. Using population genomic and phylogenetic analyses, we reconstructed the population history of the different forms. We confirmed the monophyly of the complex and found that one of the lowland forms is paraphyletic and basal relative to others, a pattern highly consistent with in situ divergence. Our results suggest initial colonization of the island through the lowlands, followed by expansion into the highlands, which led to the evolution of a distinct geographical form, genetically and ecologically different from the lowland ones. Lowland forms seem to have experienced periods of geographical isolation, but they diverged from one another by sexual selection rather than niche change. Overall, low dispersal capabilities in this island bird combined with both geographical and ecological opportunities seem to explain how divergence occurred at such a small spatial scale.

Oak symbolism in the light of genomics.

Throughout the Northern Hemisphere, human societies, political systems, and religions have appropriated oaks in symbolic representations. In this review, we explore the possible associations between recent genetic and genomic findings and the symbolic representations of oaks. We first consider the ways in which evolutionary history during the Holocene has tightened links between humans and oaks in Europe, and how this may have led to symbolic representations. We then show how recent findings concerning the structure and evolution of the oak genome have provided additional knowledge about symbolic representations, such as longevity, cohesiveness, and robustness.

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.

Population structure of Venturia inaequalis, a causal agent of apple scab, in response to heterogeneous apple tree cultivation.

BACKGROUND: Tracking newly emergent virulent populations in agroecosystems provides an opportunity to increase our understanding of the co-evolution dynamics of pathogens and their hosts. On the one hand host plants exert selective pressure on pathogen populations, thus dividing them into subpopulations of different virulence, while on the other hand they create an opportunity for secondary contact between the two divergent populations on one tree. The main objectives of the study were to explore whether the previously reported structure between two Venturia inaequalis population types, virulent or avirulent towards Malus x domestica cultivars carrying Rvi6 gene, is maintained or broken several years after the first emergence of new virulent strains in Poland, and to investigate the relationship between 'new' and 'native' populations derived from the same commercial orchards. For this purpose, we investigated the genetic structure of populations of the apple scab fungus, occurring on apple tree cultivars containing Rvi6, Rvi1 or Rvi17 resistance gene or no resistance at all, based on microsatellite data obtained from 606 strains sampled in 10 orchards composed of various host cultivars. RESULTS: Application of genetic distance inferring and clustering methods allowed us to observe clear genetic distinctness of the populations virulent towards cultivars carrying Rvi6 gene from the Rvi6-avirulent populations and substructures within the Rvi6-group as a consequence of independent immigration events followed by rare, long-distance dispersals. We did not observe such a structuring effect of other genes determining apple scab resistance on any other populations, which in turn were genetically homogenous. However, in two orchards the co-occurrence of strains of different virulence pattern on the same trees was detected, blurring the genetic boundaries between populations. CONCLUSIONS: Among several resistance genes studied, only Rvi6 exerted selective pressure on pathogens populations: those virulent toward Rvi6 hosts show unique and clear genetic and virulence pattern. For the first time in commercial Malus x domestica orchards, we reported secondary contacts between populations virulent and avirulent toward Rvi6 hosts. These two populations, first diverged in allopatry, second came into contact and subsequently began interbreeding, in such way that they show unambiguous footprints of gene flow today.

High-Throughput DNA sequencing of ancient wood.

Reconstructing the colonization and demographic dynamics that gave rise to extant forests is essential to forecasts of forest responses to environmental changes. Classical approaches to map how population of trees changed through space and time largely rely on pollen distribution patterns, with only a limited number of studies exploiting DNA molecules preserved in wooden tree archaeological and subfossil remains. Here, we advance such analyses by applying high-throughput (HTS) DNA sequencing to wood archaeological and subfossil material for the first time, using a comprehensive sample of 167 European white oak waterlogged remains spanning a large temporal (from 550 to 9,800 years) and geographical range across Europe. The successful characterization of the endogenous DNA and exogenous microbial DNA of 140 (~83%) samples helped the identification of environmental conditions favouring long-term DNA preservation in wood remains, and started to unveil the first trends in the DNA decay process in wood material. Additionally, the maternally inherited chloroplast haplotypes of 21 samples from three periods of forest human-induced use (Neolithic, Bronze Age and Middle Ages) were found to be consistent with those of modern populations growing in the same geographic areas. Our work paves the way for further studies aiming at using ancient DNA preserved in wood to reconstruct the micro-evolutionary response of trees to climate change and human forest management.

Extensive recent secondary contacts between four European white oak species.

Historical trajectories of tree species during the late Quaternary have been well reconstructed through genetic and palaeobotanical studies. However, many congeneric tree species are interfertile, and the timing and contribution of introgression to species divergence during their evolutionary history remains largely unknown. We quantified past and current gene flow events between four morphologically divergent oak species (Quercus petraea, Q. robur, Q. pyrenaica, Q. pubescens), by two independent inference methods: diffusion approximation to the joint frequency spectrum (∂a∂i) and approximate Bayesian computation (ABC). For each pair of species, alternative scenarios of speciation allowing gene flow over different timescales were evaluated. Analyses of 3524 single nucleotide polymorphisms (SNPs) randomly distributed in the genome, showed that these species evolved in complete isolation for most of their history, but recently came into secondary contact, probably facilitated by the most recent period of postglacial warming. We demonstrated that: there was sufficient genetic differentiation before secondary contact for the accumulation of barriers to gene flow; and current European white oak genomes are a mosaic of genes that have crossed species boundaries and genes impermeable to gene flow.

Back to America: tracking the origin of European introduced populations of Quercus rubra L.

Quercus rubra has been introduced in Europe since the end of the 17th century. It is widely distributed today across this continent and considered invasive in some countries. Here, we investigated the distribution of genetic diversity of both native and introduced populations with the aim of tracing the origin of introduced populations. A large sampling of 883 individuals from 73 native and 38 European locations were genotyped at 69 SNPs. In the natural range, we found a continuous geographic gradient of variation with a predominant latitudinal component. We explored the existence of ancestral populations by performing Bayesian clustering analysis and found support for two or three ancestral genetic clusters. Approximate Bayesian Computations analyses based on these two or three clusters support recent extensive secondary contacts between them, suggesting that present-day continuous genetic variation resulted from recent admixture. In the introduced range, one main genetic cluster was not recovered in Europe, suggesting that source populations were preferentially located in the northern part of the natural distribution. However, our results cannot refute the introduction of populations from the southern states that did not survive in Europe.

Emergence of new virulent populations of apple scab from nonagricultural disease reservoirs.

Plant pathogens adapt readily to new crop varieties in agrosystems, and it is crucial to understand the factors underlying the epidemic spread of new virulent strains if we are to develop more efficient strategies to control them. In this study we used multilocus microsatellite typing, molecular epidemiology tools and a large collection of isolates from cultivated, wild and ornamental apples to investigate the origin of new virulent populations of Venturia inaequalis, an ascomycete fungus causing apple scab on varieties carrying the Rvi6 resistance gene. We demonstrated a common origin at the European scale of populations infecting apples (Malus × domestica) carrying the Rvi6 resistance and Malus floribunda, the progenitor of the Rvi6 resistance. Demographic modeling indicated that the Rvi6-virulent lineage separated several thousands of years ago from populations infecting non-Rvi6 hosts, without detectable gene flow between the two lineages. These findings show that 'breakdowns' of plant resistance genes can be caused by the selection and migration of virulent genotypes from standing genetic variation maintained in environmental disease reservoirs, here ornamental crabapples. This work stresses the need to take better account of pathogen diversity in resistance screenings of breeding lines and in resistance deployment strategies, in order to enhance sustainable disease management.

When virulence originates from nonagricultural hosts: evolutionary and epidemiological consequences of introgressions following secondary contacts in Venturia inaequalis.

In pathogens, introgressions through secondary contacts between divergent populations from agricultural and nonagricultural disease reservoirs are expected to have crucial evolutionary and epidemiological implications. Despite the importance of this question for disease management, experimental demonstrations of these implications remain scarce. Recently, we identified a virulent population of the apple scab pathogen Venturia inaequalis that migrated from nonagricultural hosts to European domestic apple orchards. Here, we investigated the occurrence of gene flow between agricultural and nonagricultural populations sampled in two orchards, and thereafter its consequences on the pathogenicity of hybrids. Population genetic structure and demographic inferences based on the genotypes of 104 strains revealed a high amount of gene flow between the two populations in one orchard. In this site, mating between populations was made possible by the presence of a common host. Our results revealed an invasion of the virulent trait in the agricultural population; a main direction of introgression in hybrids from the agricultural to nonagricultural genetic backgrounds; and a population of hybrids with transgressive traits. We demonstrate a secondary contact with gene flow between divergent populations of pathogens. Our findings highlight evolutionary and epidemiological changes in pathogens and have concrete implications for sustainable disease management.

The genetic structure of a Venturia inaequalis population in a heterogeneous host population composed of different Malus species.

BACKGROUND: Adaptation, which induces differentiation between populations in relation to environmental conditions, can initiate divergence. The balance between gene flow and selection determines the maintenance of such a structure in sympatry. Studying these two antagonistic forces in plant pathogens is made possible because of the high ability of pathogens to disperse and of the strong selective pressures exerted by their hosts. In this article, we analysed the genetic structure of the population of the apple scab fungus, Venturia inaequalis, in a heterogeneous environment composed of various Malus species. Inferences were drawn from microsatellite and AFLP data obtained from 114 strains sampled in a single orchard on nine different Malus species to determine the forces that shape the genetic structure of the pathogen. RESULTS: Using clustering methods, we first identified two specialist subpopulations: (i) a virulent subpopulation sampled on Malus trees carrying the Rvi6 resistance gene; and (ii) a subpopulation infecting only Malus trees that did not carry this resistance gene. A genome scan of loci on these two subpopulations did not detect any locus under selection. Additionally, we did not detect any other particular substructure linked to different hosts. However, an isolation-by-distance (IBD) pattern at the orchard scale revealed free gene flow within each subpopulation. CONCLUSIONS: Our work shows a rare example of a very strong effect of a resistance gene on pathogen populations. Despite the high diversity of Malus hosts, the presence of Rvi6 seems sufficient to explain the observed genetic structure. Moreover, detection of an IBD pattern at the orchard scale revealed a very low average dispersal distance that is particularly significant for epidemiologists and landscape managers for the design of scab control strategies.