Adaptive introgression facilitate adaptation to high latitudes in European aspen (Populus tremula L.)

Ploidy as a leaky reproductive barrier: mechanisms, rates and evolutionary significance of interploidy gene flow

BACKGROUND

Whole-genome duplication (polyploidization) is a dominant force in sympatric speciation, particularly in plants. Genome doubling instantly poses a barrier to gene flow owing to the strong crossing incompatibilities between individuals differing in ploidy. The strength of the barrier, however, varies from species to species and recent genetic investigations revealed cases of rampant interploidy introgression in multiple ploidy-variable species.

SCOPE

Here, we review novel insights into the frequency of interploidy gene flow in natural systems and summarize the underlying mechanisms promoting interploidy gene flow. Field surveys, occasionally complemented by crossing experiments, suggest frequent opportunities for interploidy gene flow, particularly in the direction from diploid to tetraploid, and between (higher) polyploids. However, a scarcity of accompanying population genetic evidence and a virtual lack of integration of these approaches leave the underlying mechanisms and levels of realized interploidy gene flow in nature largely unknown. Finally, we discuss potential consequences of interploidy genome permeability on polyploid speciation and adaptation and highlight novel avenues that have just recently been opened by the very first genomic studies of ploidy-variable species. Standing in stark contrast with rapidly accumulating evidence for evolutionary importance of homoploid introgression, similar cases in ploidy-variable systems are yet to be documented.

CONCLUSIONS

The genomics era provides novel opportunity to re-evaluate the role of interploidy introgression in speciation and adaptation. To achieve this goal, interdisciplinary studies bordering ecology and population genetics and genomics are needed.

Patterns of genetic variation and local adaptation of a native herbivore to a lethal invasive plant

Understanding the evolutionary processes that influence fitness is critical to predicting species' responses to selection. Interactions among evolutionary processes including gene flow, drift and the strength of selection can lead to either local adaptation or maladaptation, especially in heterogenous landscapes. Populations experiencing novel environments or resources are ideal for understanding the mechanisms underlying adaptation or maladaptation, specifically in locally co-evolved interactions. We used the interaction between a native herbivore that oviposits on a patchily distributed introduced plant that in turn causes significant mortality to the larvae to test for signatures of local adaptation in areas where the two co-occurred. We used whole-genome sequencing to explore population structure, patterns of gene flow and signatures of local adaptation. We found signatures of local adaptation in response to the introduced plant in the absence of strong population structure with no genetic differentiation and low genetic variation. Additionally, we found localized allele frequency differences within a single population between habitats with and without the lethal plant, highlighting the effects of strong selection. Finally, we identified that selection was acting on larval ability to feed on the plant rather than on females' ability to avoid oviposition, thus uncovering the specific ontogenetic target of selection. Our work highlights the potential for adaptation to occur in a fine-grained landscape in the presence of gene flow and low genetic variation.

Evolution of the Correlated Genomic Variation Landscape Across a Divergence Continuum in the Genus Castanopsis

The heterogeneous landscape of genomic variation has been well documented in population genomic studies. However, disentangling the intricate interplay of evolutionary forces influencing the genetic variation landscape over time remains challenging. In this study, we assembled a chromosome-level genome for Castanopsis eyrei and sequenced the whole genomes of 276 individuals from 12 Castanopsis species, spanning a broad divergence continuum. We found highly correlated genomic variation landscapes across these species. Furthermore, variations in genetic diversity and differentiation along the genome were strongly associated with recombination rates and gene density. These results suggest that long-term linked selection and conserved genomic features have contributed to the formation of a common genomic variation landscape. By examining how correlations between population summary statistics change throughout the species divergence continuum, we determined that background selection alone does not fully explain the observed patterns of genomic variation; the effects of recurrent selective sweeps must be considered. We further revealed that extensive gene flow has significantly influenced patterns of genomic variation in Castanopsis species. The estimated admixture proportion correlated positively with recombination rate and negatively with gene density, supporting a scenario of selection against gene flow. Additionally, putative introgression regions exhibited strong signals of positive selection, an enrichment of functional genes, and reduced genetic burdens, indicating that adaptive introgression has played a role in shaping the genomes of hybridizing species. This study provides insights into how different evolutionary forces have interacted in driving the evolution of the genomic variation landscape.

An Improved Chromosome-scale Genome Assembly and Population Genetics resource for Populus tremula

Aspen (Populus tremula L.) is a keystone species and a model system for forest tree genomics. We present an updated resource comprising a chromosome-scale assembly, population genetics and genomics data. Using the resource, we explore the genetic basis of natural variation in leaf size and shape, traits with complex genetic architecture. We generated the genome assembly using long-read sequencing, optical and high-density genetic maps. We conducted whole-genome resequencing of the Umeå Aspen (UmAsp) collection. Using the assembly and re-sequencing data from the UmAsp, Swedish Aspen (SwAsp) and Scottish Aspen (ScotAsp) collections we performed genome-wide association analyses (GWAS) using Single Nucleotide Polymorphisms (SNPs) for 26 leaf physiognomy phenotypes. We conducted Assay of Transposase Accessible Chromatin sequencing (ATAC-Seq), identified genomic regions of accessible chromatin, and subset SNPs to these regions, improving the GWAS detection rate. We identified candidate long non-coding RNAs in leaf samples, quantified their expression in an updated co-expression network, and used this to explore the functions of candidate genes identified from the GWAS. A GWAS found SNP associations for seven traits. The associated SNPs were in or near genes annotated with developmental functions, which represent candidates for further study. Of particular interest was a ~177-kbp region harbouring associations with several leaf phenotypes in ScotAsp. We have incorporated the assembly, population genetics, genomics, and GWAS data into the PlantGenIE.org web resource, including updating existing genomics data to the new genome version, to enable easy exploration and visualisation. We provide all raw and processed data to facilitate reuse in future studies.

Why we thrive beneath a northern sky - genomic signals of selection in apple for adaptation to northern Sweden

Good understanding of the genomic regions underlying adaptation of apple to boreal climates is needed to facilitate efficient breeding of locally adapted apple cultivars. Proper infrastructure for phenotyping and evaluation is essential for identification of traits responsible for adaptation, and dissection of their genetic composition. However, such infrastructure is costly and currently not available for the boreal zone of northern Sweden. Therefore, we used historical pomological data on climate adaptation of 59 apple cultivars and whole genome sequencing to identify genomic regions that have undergone historical selection among apple cultivars recommended for cultivation in northern Sweden. We found the apple collection to be composed of two ancestral groups that are largely concordant with the grouping into 'hardy' and 'not hardy' cultivars based on the pomological literature. Using a number of genome-wide scans for signals of selection, we obtained strong evidence of positive selection at a genomic region around 29 MbHFTH1 of chromosome 1 among apple cultivars in the 'hardy' group. Using phased genotypic data from the 20 K apple Infinium® SNP array, we identified haplotypes associated with the two cultivar groups and traced transmission of these haplotypes through the pedigrees of some apple cultivars. This demonstrates that historical data from pomological literature can be analyzed by population genomic approaches as a step towards revealing the genomic control of a key property for a horticultural niche market. Such knowledge is needed to facilitate efficient breeding strategies for development of locally adapted apple cultivars in the future. The current study illustrates the response to a very strong selective pressure imposed on tree crops by climatic factors, and the importance of genetic research on this topic and feasibility of breeding efforts in the light of the ongoing climate change.

Genomic evidence for evolutionary history and local adaptation of two endemic apricots: Prunus hongpingensis and P. zhengheensis

Apricot, belonging to the Armeniaca section of Rosaceae, is one of the economically important crop fruits that has been extensively cultivated. The natural wild apricots offer valuable genetic resources for crop improvement. However, some of them are endemic, with small populations, and are even at risk of extinction. In this study we unveil chromosome-level genome assemblies for two southern China endemic apricots, Prunus hongpingensis (PHP) and P. zhengheensis (PZH). We also characterize their evolutionary history and the genomic basis of their local adaptation using whole-genome resequencing data. Our findings reveal that PHP and PZH are closely related to Prunus armeniaca and form a distinct lineage. Both species experienced a decline in effective population size following the Last Glacial Maximum (LGM), which likely contributed to their current small population sizes. Despite the observed decrease in genetic diversity and heterozygosity, we do not observe an increased accumulation of deleterious mutations in these two endemic apricots. This is likely due to the combined effects of a low inbreeding coefficient and strong purifying selection. Furthermore, we identify a set of genes that have undergone positive selection and are associated with local environmental adaptation in PHP and PZH, respectively. These candidate genes can serve as valuable genetic resources for targeted breeding and improvement of cultivated apricots. Overall, our study not only enriches our comprehension of the evolutionary history of apricot species but also offers crucial insights for the conservation and future breeding of other endemic species amidst rapid climate changes.

The role of introgressive hybridization in shaping the geographically isolated gene pools of wax palm populations (genus Ceroxylon)

The speciation continuum is the process by which genetic groups diverge until they reach reproductive isolation. It has become common in the literature to show that this process is gradual and flickering, with possibly many instances of secondary contact and introgression after divergence has started. The level of divergence might vary among genomic regions due to, among others, the different forces and roles of selection played by the shared regions. Through hybrid capture, we sequenced ca. 4,000 nuclear regions in populations of six species of wax palms, five of which form a monophyletic group (genus Ceroxylon, Arecaceae: Ceroxyloideae). We show that in this group, the different populations show varying degrees of introgressive hybridization, and two of them are backcrosses of the other three 'pure' species. This is particularly interesting because these three species are dioecious, have a shared main pollinator, and have slightly overlapping reproductive seasons but highly divergent morphologies. Our work supports shows wax palms diverge under positive and background selection in allopatry, and hybridize due to secondary contact and inefficient reproductive barriers, which sustain genetic diversity. Introgressed regions are generally not under positive selection. Peripheral populations are backcrosses of other species; thus, introgressive hybridization is likely modulated by demographic effects rather than selective pressures. In general, these species might function as an 'evolutionary syngameon' where expanding, peripheral, small, and isolated populations maintain diversity by crossing with available individuals of other wax palms. In the Andean context, species can benefit from gained variation from a second taxon or the enhancement of population sizes by recreating a common genetic pool.

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss)

BACKGROUND AND AIMS

New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae.

METHODS

We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally.

KEY RESULTS

Using sliding-window population genetic analyses we find that differentiation is concentrated within 'islands' of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap.

CONCLUSIONS

We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth's climate rapidly changes.

Introgression and persistence of cultivar alleles in wild carrot (Daucus carota) populations in the United States

PREMISE

Cultivated species and their wild relatives often hybridize in the wild, and the hybrids can survive and reproduce in some environments. However, it is unclear whether cultivar alleles are permanently incorporated into the wild genomes or whether they are purged by natural selection. This question is key to accurately assessing the risk of escape and spread of cultivar genes into wild populations.

METHODS

We used genomic data and population genomic methods to study hybridization and introgression between cultivated and wild carrot (Daucus carota) in the United States. We used single nucleotide polymorphisms (SNPs) obtained via genotyping by sequencing for 450 wild individuals from 29 wild georeferenced populations in seven states and 144 cultivars from the United States, Europe, and Asia.

RESULTS

Cultivated and wild carrot formed two genetically differentiated groups, and evidence of crop-wild admixture was detected in several but not all wild carrot populations in the United States. Two regions were identified where cultivar alleles were present in wild carrots: California and Nantucket Island (Massachusetts). Surprisingly, there was no evidence of introgression in some populations with a long-known history of sympatry with the crop, suggesting that post-hybridization barriers might prevent introgression in some areas.

CONCLUSIONS

Our results provide support for the introgression and long-term persistence of cultivar alleles in wild carrots populations. We thus anticipate that the release of genetically engineered (GE) cultivars would lead to the introduction and spread of GE alleles in wild carrot populations.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

Genetic markers and tree properties predicting wood biorefining potential in aspen (Populus tremula) bioenergy feedstock

BACKGROUND

Wood represents the majority of the biomass on land and constitutes a renewable source of biofuels and other bioproducts. However, wood is recalcitrant to bioconversion, raising a need for feedstock improvement in production of, for instance, biofuels. We investigated the properties of wood that affect bioconversion, as well as the underlying genetics, to help identify superior tree feedstocks for biorefining.

RESULTS

We recorded 65 wood-related and growth traits in a population of 113 natural aspen genotypes from Sweden ( https://doi.org/10.5061/dryad.gtht76hrd ). These traits included three growth and field performance traits, 20 traits for wood chemical composition, 17 traits for wood anatomy and structure, and 25 wood saccharification traits as indicators of bioconversion potential. Glucose release after saccharification with acidic pretreatment correlated positively with tree stem height and diameter and the carbohydrate content of the wood, and negatively with the content of lignin and the hemicellulose sugar units. Most of these traits displayed extensive natural variation within the aspen population and high broad-sense heritability, supporting their potential in genetic improvement of feedstocks towards improved bioconversion. Finally, a genome-wide association study (GWAS) revealed 13 genetic loci for saccharification yield (on a whole-tree-biomass basis), with six of them intersecting with associations for either height or stem diameter of the trees.

CONCLUSIONS

The simple growth traits of stem height and diameter were identified as good predictors of wood saccharification yield in aspen trees. GWAS elucidated the underlying genetics, revealing putative genetic markers for bioconversion of bioenergy tree feedstocks.

A Genomic Quantitative Study on the Contribution of the Ancestral-State Bases Relative to Derived Bases in the Divergence and Local Adaptation of Populus davidiana

Identifying alleles associated with adaptation to new environments will advance our understanding of evolution from the molecular level. Previous studies have found that the Populus davidiana southwest population in East Asia has differentiated from other populations in the range. We aimed to evaluate the contributions of the ancestral-state bases (ASBs) relative to derived bases (DBs) in the local adaptation of P. davidiana in the Yunnan-Guizhou Plateau from a quantitative perspective using whole-genome re-sequencing data from 90 P. davidiana samples from three regions across the species range. Our results showed that the uplift of the Qinghai-Tibet Plateau during the Neogene and associated climate fluctuations during the Middle Pleistocene were likely an important factor in the early divergence of P. davidiana. Highly differentiated genomic regions between populations were inferred to have undergone strong linked natural selection, and ASBs are the chief means by which populations of P. davidiana adapt to novel environmental conditions; however, when adapting to regions with high environmental differences relative to the ancestral range, the proportion of DBs was significantly higher than that of background regions, as ASBs are insufficient to cope with these environments. Finally, a number of genes were identified in the outlier region.

Differences in heat tolerance, water use efficiency and growth among Douglas-fir families and varieties evidenced by GWAS and common garden studies

Severe and frequent heat and drought events challenge the survival and development of long-generation trees. In this study, we investigated the genomic basis of heat tolerance, water use efficiency and growth by performing genome-wide association studies in coastal Douglas-fir (Pseudotsuga menziesii) and intervarietal (menziesii × glauca) hybrid seedlings. GWAS results identified 32 candidate genes involved in primary and secondary metabolism, abiotic stress and signaling, among other functions. Water use efficiency (inferred from carbon isotope discrimination), photosynthetic capacity (inferred from %N), height and heat tolerance (inferred from electrolyte leakage in a heat stress experiment) were significantly different among Douglas-fir families and varieties. High-elevation seed sources had increased water use efficiency, which could be a result of higher photosynthetic capacity. Similarly, families with greater heat tolerance also had higher water use efficiency and slower growth, suggesting a conservative growth strategy. Intervarietal hybrids showed increased heat tolerance (lower electrolyte leakage at 50 and 55 °C) and higher water use efficiency compared with coastal families, suggesting that hybridization might be a source of pre-adapted alleles to warming climates and should be considered for large-scale reforestation projects under increasingly arid conditions.

Genomic vulnerability to climate change in Quercus acutissima, a dominant tree species in East Asian deciduous forests

Understanding the evolutionary processes that shape the landscape of genetic variation and influence the response of species to future climate change is critical for biodiversity conservation. Here, we sampled 27 populations across the distribution range of a dominant forest tree, Quercus acutissima, in East Asia, and applied genome-wide analyses to track the evolutionary history and predict the fate of populations under future climate. We found two genetic groups (East and West) in Q. acutissima that diverged during Pliocene. We also found a heterogeneous landscape of genomic variation in this species, which may have been shaped by population demography and linked selections. Using genotype-environment association analyses, we identified climate-associated SNPs in a diverse set of genes and functional categories, indicating a model of polygenic adaptation in Q. acutissima. We further estimated three genetic offset metrics to quantify genomic vulnerability of this species to climate change due to the complex interplay between local adaptation and migration. We found that marginal populations are under higher risk of local extinction because of future climate change, and may not be able to track suitable habitats to maintain the gene-environment relationships observed under the current climate. We also detected higher reverse genetic offsets in northern China, indicating that genetic variation currently present in the whole range of Q. acutissima may not adapt to future climate conditions in this area. Overall, this study illustrates how evolutionary processes have shaped the landscape of genomic variation, and provides a comprehensive genome-wide view of climate maladaptation in Q. acutissima.

Evidence of hybrid breakdown among invasive hybrid cattails (Typha × glauca)

Interspecific hybridization has varied consequences for offspring fitness, with implications for the maintenance of species integrity. Hybrid vigour, when it occurs, can peak in first-generation (F1) hybrids and then decline in advanced-generation (F2+) hybrids. This hybrid breakdown, together with the processes affecting patterns of hybridization and hybrid fitness, determine the evolutionary stability of hybrid zones. An extensive hybrid zone in North America involving the cattails Typha latifolia, T. angustifolia, and their invasive hybrid T. × glauca is characterized by hybrid vigour among F1s, but the fitness of advanced-generation hybrids has not been studied. We compared seed germination and plant growth of T. latifolia (parental L), F1 T. × glauca (F1), hybrid backcrosses to T. angustifolia (bcA) and T. latifolia (bcL), and advanced-generation (F2) hybrids. Consistent with expectations under hybrid breakdown, we found reduced plant growth for F2 hybrids in comparison with F1s (plant height and above-ground biomass) and parental Ls (above-ground biomass). Backcrossed hybrids had intermediate measures of plant growth and bcLs were characterized by reduced seed germination in comparison with parental Ls. Hybrid breakdown could make the formation of F1s in North America finite because (1) hybridization among cattails is asymmetric, with T. angustifolia but not T. latifolia subject to genetic swamping, and (2) T. angustifolia is less common and subject to competitive displacement by F1s. Hybrid breakdown is therefore expected to reduce hybrid frequencies over time, contributing to the long-term maintenance of T. latifolia - the only native cattail in the study region.

Alternative modes of introgression-mediated selection shaped crop adaptation to novel climates

Recent plant genomic studies provide fine-grained details on the evolutionary consequences of adaptive introgression during crop domestication. Modern genomic approaches and analytical methods now make it possible to better separate the introgression signal from the demographic signal thus providing a more comprehensive and complex picture of the role of introgression in local adaptation. Adaptive introgression has been fundamental for crop expansion and has involved complex patterns of gene flow. In addition to providing new and more favorable alleles of large effect, introgression during the early stages of domestication also increased allelic diversity at adaptive loci. Previous studies have largely underestimated the effect of such increased diversity following introgression. Recent genomic studies in wheat, potato, maize, grapevine, and ryegrass show that introgression of multiple genes, of as yet unknown effect, increased the effectiveness of purifying selection, and promoted disruptive or fluctuating selection in early cultivars and landraces. Historical selection processes associated with introgression from crop wild relatives provide an instructive analog for adaptation to current climate change and offer new avenues for crop breeding research that are expected to be instrumental for strengthening food security in the coming years.

FLOWERING LOCUS T paralogs control the annual growth cycle in Populus trees

In temperate and boreal regions, perennials adapt their annual growth cycle to the change of seasons. These adaptations ensure survival in harsh environmental conditions, allowing growth at different latitudes and altitudes, and are therefore tightly regulated. Populus tree species cease growth and form terminal buds in autumn when photoperiod falls below a certain threshold.¹ This is followed by establishment of dormancy and cold hardiness over the winter. At the center of the photoperiodic pathway in Populus is the gene FLOWERING LOCUS T2 (FT2), which is expressed during summer and harbors significant SNPs in its locus associated with timing of bud set.^1-4 The paralogous gene FT1, on the other hand, is hyper-induced in chilling buds during winter.^3,5 Even though its function is so far unknown, it has been suggested to be involved in the regulation of flowering and the release of winter dormancy.^3,5 In this study, we employ CRISPR-Cas9-mediated gene editing to individually study the function of the FT-like genes in Populus trees. We show that while FT2 is required for vegetative growth during spring and summer and regulates the entry into dormancy, expression of FT1 is absolutely required for bud flush in spring. Gene expression profiling suggests that this function of FT1 is linked to the release of winter dormancy rather than to the regulation of bud flush per se. These data show how FT duplication and sub-functionalization have allowed Populus trees to regulate two completely different and major developmental control points during the yearly growth cycle.