Introgression from Populus balsamifera underlies adaptively significant variation and range boundaries in P. trichocarpa

In situ diversification and adaptive introgression in Taiwanese Scutellaria

Island habitats provide unique opportunities to study speciation. Recent work indicates that both ex situ origination and in situ speciation contribute to island species diversity. However, clear evidence of local adaptation of endemic plant species on islands requires in-depth studies, which are scarce. This study underscores the importance of local adaptation in maintaining species boundaries by examining how adaptive introgression, hybridization, and local adaptation contribute to genetic variation in island species. Multilocus genome scanning of 51 nuclear genes was used to investigate the evolutionary relationships of the Scutellaria species complex on Taiwan Island and assess the role of in situ diversification in generating high endemism and genetic diversity. Interspecies introgressions were detected by phylogenetic networks and ABBA-BABA-based analysis, suggesting ongoing or recent speciation processes. Coalescent-based simulation identified hybrid speciation in Scutellaria taiwanensis and Scutellaria hsiehii, with evidence of hybridization between more than two parental species. Genotype-environment association studies revealed that the influence of climate, particularly precipitation- and temperature-related factors, contributed to adaptive genetic divergence between species. Additionally, adaptive introgression related to environmental pressures that may have facilitated the colonization of new island habitats were identified. This research illustrates how hybridization, introgression, and adaptation shaped the evolutionary histories and divergence of this island-endemic plant species complex and sheds light on the multifaceted mechanisms of speciation on semi-isolated islands.

Insular environment-dependent introgression from an arid-grassland orchid to a wetland orchid on an oceanic island

Adaptive introgression plays a vital role in allowing recipient species to adapt and colonize new environments. However, our understanding of such environment-dependent introgressions is primarily limited to specific plant taxa in particular settings. In Japan, two related orchid species, the autonomously self-pollinating Pogonia minor and the outcrossing Pogonia japonica, typically inhabit dry grasslands and wetlands, respectively. Intriguingly, an island ecotype of P. japonica exists in arid, wind-swept, open sites on volcanic mountain slopes on Kozu Island, in the oceanic Izu Islands. To investigate potential introgression and its implications between P. japonica and P. minor on Kozu Island, we applied a comprehensive approach that included examining morphological traits, genome-wide SNP data, and plastid DNA sequences. We also examined the breeding systems of these species on Kozu Island through artificial pollination experiments to determine if introgression from P. minor has endowed the P. japonica ecotype with selfing capabilities. Extensive sampling on Kozu Island revealed that all P. japonica specimens exhibit signs of introgression from P. minor, suggesting the absence of pure P. japonica populations on the island. Furthermore, the chloroplast haplotypes of the insular P. japonica ecotype consistently match those of P. minor, indicating a predominantly asymmetrical initial hybridization with P. minor acting mainly as the maternal parent in the formation of F1 hybrids. Despite the advantages of self-fertilization in isolated environments, the insular P. japonica does not exhibit autogamy. Consequently, the scarcity of moist habitats, rather than selection pressure for selfing, likely contributes to the observed widespread introgression. Our study strongly suggests that the arid-environment-adapted P. minor has introgressed into the insular ecotype of P. japonica, enabling its successful colonization of arid volcanic mountain slopes of the oceanic island.

Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles

Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids, it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.

Comparative population genomics reveals convergent and divergent selection in the apricot-peach-plum-mei complex

The economically significant genus Prunus includes fruit and nut crops that have been domesticated for shared and specific agronomic traits; however, the genomic signals of convergent and divergent selection have not been elucidated. In this study, we aimed to detect genomic signatures of convergent and divergent selection by conducting comparative population genomic analyses of the apricot-peach-plum-mei (APPM) complex, utilizing a haplotype-resolved telomere-to-telomere (T2T) genome assembly and population resequencing data. The haplotype-resolved T2T reference genome for the plum cultivar was assembled through HiFi and Hi-C reads, resulting in two haplotypes 251.25 and 251.29 Mb in size, respectively. Comparative genomics reveals a chromosomal translocation of ~1.17 Mb in the apricot genomes compared with peach, plum, and mei. Notably, the translocation involves the D locus, significantly impacting titratable acidity (TA), pH, and sugar content. Population genetic analysis detected substantial gene flow between plum and apricot, with introgression regions enriched in post-embryonic development and pollen germination processes. Comparative population genetic analyses revealed convergent selection for stress tolerance, flower development, and fruit ripening, along with divergent selection shaping specific crop, such as somatic embryogenesis in plum, pollen germination in mei, and hormone regulation in peach. Notably, selective sweeps on chromosome 7 coincide with a chromosomal collinearity from the comparative genomics, impacting key fruit-softening genes such as PG, regulated by ERF and RMA1H1. Overall, this study provides insights into the genetic diversity, evolutionary history, and domestication of the APPM complex, offering valuable implications for genetic studies and breeding programs of Prunus crops.

Conservation genomics provides insights into genetic resilience and adaptation of the endangered Chinese hazelnut, Corylus chinensis

Global climate change has increased concerns regarding biodiversity loss. However, many key conservation issues still required further research, including demographic history, deleterious mutation load, adaptive evolution, and putative introgression. Here we generated the first chromosome-level genome of the endangered Chinese hazelnut, Corylus chinensis, and compared the genomic signatures with its sympatric widespread C. kwechowensis-C. yunnanensis complex. We found large genome rearrangements across all Corylus species and identified species-specific expanded gene families that may be involved in adaptation. Population genomics revealed that both C. chinensis and the C. kwechowensis-C. yunnanensis complex had diverged into two genetic lineages, forming a consistent pattern of southwestern-northern differentiation. Population size of the narrow southwestern lineages of both species have decreased continuously since the late Miocene, whereas the widespread northern lineages have remained stable (C. chinensis) or have even recovered from population bottlenecks (C. kwechowensis-C. yunnanensis complex) during the Quaternary. Compared with C. kwechowensis-C. yunnanensis complex, C. chinensis showed significantly lower genomic diversity and higher inbreeding level. However, C. chinensis carried significantly fewer deleterious mutations than C. kwechowensis-C. yunnanensis complex, as more effective purging selection reduced the accumulation of homozygous variants. We also detected signals of positive selection and adaptive introgression in different lineages, which facilitated the accumulation of favorable variants and formation of local adaptation. Hence, both types of selection and exogenous introgression could have mitigated inbreeding and facilitated survival and persistence of C. chinensis. Overall, our study provides critical insights into lineage differentiation, local adaptation, and the potential for future recovery of endangered trees.

The genomes of Darwin's primroses reveal chromosome-scale adaptive introgression and differential permeability of species boundaries

Introgression is an important source of genetic variation that can determine species adaptation to environmental conditions. Yet, definitive evidence of the genomic and adaptive implications of introgression in nature remains scarce. The widespread hybrid zones of Darwin's primroses (Primula elatior, Primula veris, and Primula vulgaris) provide a unique natural laboratory for studying introgression in flowering plants and the varying permeability of species boundaries. Through analysis of 650 genomes, we provide evidence of an introgressed genomic region likely to confer adaptive advantage in conditions of soil toxicity. We also document unequivocal evidence of chloroplast introgression, an important precursor to species-wide chloroplast capture. Finally, we provide the first evidence that the S-locus supergene, which controls heterostyly in primroses, does not introgress in this clade. Our results contribute novel insights into the adaptive role of introgression and demonstrate the importance of extensive genomic and geographical sampling for illuminating the complex nature of species boundaries.

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.

Genomic divergence and introgression among three Populus species

Genomic divergence with gene flow is very common in both plants and animals. However, divergence and gene flow are two counteracting factors during speciation. Identifying the types of genes that are likely to be introgressed and what genetic factors restrict further effective reproduction of interspecific hybrids is of great interest to biologists. We aimed to address these issues using three related tree species, Populus alba (Pa), P. tremula (Pt), and P. tremuloides (Ps), and the interspecific hybrid of the former two species, P. × canescens (Pc). We collected 105 genomes for these four poplar lineages, including 28 Pa, 38Pt, 21 Ps, and 18 Pc individuals, to reconstruct their evolutionary histories. Our coalescence-based simulations indicated that Pa diverged earliest from Ps and Pt, and asymmetrical gene flow existed between any two lineages, with especially large ancient gene flow occurring between Pa and Pt. The genomic landscape of divergence between pairs of the three species are highly heterogeneous, which may have arisen through both divergent sorting of ancient polymorphisms and ongoing gene flow. We found that extant regions of the genome with introgressed ancestry reduced genetic divergence but elevated recombination rates and accounted for 5.76 % of the total genome. Introgressed genes were functionally associated with stress resistance, including innate immune response, anti-adversity response, and programmed cell death. However, candidate genes underlying postmating barriers of Pc were homozygous and resistant to introgression due to the incompatibility of alleles between loci after hybridization and were associated with endosperm and gamete formation and disease resistance. Our study revealed genomic dynamics during speciation with gene flow and identified regions of the genome that were likely introgressed and adaptive as well as candidate loci responsible for hybrid incompatibility that resulted in the formation of postmating barriers after hybridization.

Genomic evidence of speciation by fusion in a recent radiation of grasshoppers

Postdivergence gene flow can trigger a number of creative evolutionary outcomes, ranging from the transfer of beneficial alleles across species boundaries (i.e., adaptive introgression) to the formation of new species (i.e., hybrid speciation). Although neutral and adaptive introgression has been broadly documented in nature, hybrid speciation is assumed to be rare and the evolutionary and ecological context facilitating this phenomenon still remains controversial. Through combining genomic and phenotypic data, we evaluate the hypothesis that the dual feeding regime (based on both scrub legumes and gramineous herbs) of the taxonomically controversial grasshopper Chorthippus saulcyi algoaldensis resulted from hybridization between the sister taxa C. binotatus (that exclusively feeds on scrub legumes) and C. saulcyi (that only feeds on gramineous herbs). Genetic clustering analyses and inferences from coalescent-based demographic simulations confirm that C. s. algoaldensis represents an independently evolving lineage and support the ancient hybrid origin of this taxon (about 1.4 Ma), which sheds light on its uncertain phylogenetic position and might explain its broader trophic niche. We propose a Pleistocene hybrid speciation model where range shifts resulting from climatic oscillations can promote the formation of hybrid swarms and facilitate their long-term persistence through geographic isolation from parental forms in topographically complex landscapes.

Genome-wide analyses of introgression between two sympatric Asian oak species

Introgression can be an important source of new alleles for adaption under rapidly changing environments, perhaps even more important than standing variation. Though introgression has been extensively studied in many plants and animals, key questions on the underlying mechanisms of introgression still remain unanswered. In particular, we are yet to determine the genomic distribution of introgressed regions along the genome; whether the extent and patterns of introgression are influenced by ecological factors; and when and how introgression contributes to adaptation. Here, we generated high-quality genomic resources for two sympatric widespread Asian oak species, Quercus acutissima and Q. variabilis, sampled in multiple forests to study introgression between them. We show that introgressed regions are broadly distributed across the genome. Introgression was affected by genetic divergence between pairs of populations and by the similarity of the environments in which they live-populations occupying similar ecological sites tended to share the same introgressed regions. Introgressed genomic footprints of adaptation were preferentially located in regions with suppressed recombination rate. Introgression probably confers adaptation in these oak populations by introducing allelic variation in cis-regulatory elements, in particular through transposable element insertions, thereby altering the regulation of genes related to stress. Our results provide new avenues of research for uncovering mechanisms of adaptation due to hybridization in sympatric species.

Expansion of the rare Eucalyptus risdonii under climate change through hybridization with a closely related species despite hybrid inferiority

BACKGROUND AND AIMS

Hybridization is increasingly recognized as an integral part of the dynamics of species range expansion and contraction. Thus, it is important to understand the reproductive barriers between co-occurring species. Extending previous studies that argued that the rare Eucalyptus risdonii was expanding into the range of the surrounding E. amygdalina by both seed and pollen dispersal, we here investigate the long-term fitness of both species and their hybrids and whether expansion is continuing.

METHODS

We assessed the survival of phenotypes representing a continuum between the two pure species in a natural hybrid swarm after 29 years, along with seedling recruitment. The performance of pure species as well as of artificial and natural hybrids was also assessed over 28 years in a common garden trial.

KEY RESULTS

In the hybrid zone, E. amygdalina adults showed greater mortality than E. risdonii, and the current seedling cohort is still dominated by E. risdonii phenotypes. Morphologically intermediate individuals appeared to be the least fit. Similar results were observed after growing artificial first-generation and natural hybrids alongside pure species families in a common garden trial. Here, the survival, reproduction, health and growth of the intermediate hybrids were significantly less than those of either pure species, consistent with hybrid inferiority, although this did not manifest until later reproductive ages. Among the variable progeny of natural intermediate hybrids, the most E. risdonii-like phenotypes were the most fit.

CONCLUSIONS

This study contributes to the increasing number of reports of hybrid inferiority in Eucalyptus, suggesting that post-zygotic barriers contribute to the maintenance of species integrity even between closely related species. However, with fitness rapidly recovered following backcrossing, it is argued that hybridization can still be an important evolutionary process, in the present case appearing to contribute to the range expansion of the rare E. risdonii in response to climate change.

Estimating the time since admixture from phased and unphased molecular data

After admixture, recombination breaks down genomic blocks of contiguous ancestry. The breakdown of these blocks forms a new “molecular clock” that ticks at a much faster rate than the mutation clock, enabling accurate dating of admixture events in the recent past. However, existing theory on the breakdown of these blocks, or the accumulation of delineations between blocks, so‐called “junctions”, has mostly been limited to using regularly spaced markers on phased data. Here, we present an extension to the theory of junctions using the ancestral recombination graph that describes the expected number of junctions for any distribution of markers along the genome. Furthermore, we provide a new framework to infer the time since admixture using unphased data. We demonstrate both the phased and unphased methods on simulated data and show that our new extensions have improved accuracy with respect to previous methods, especially for smaller population sizes and more ancient admixture times. Lastly, we demonstrate the applicability of our method on three empirical data sets, including labcrosses of yeast (Saccharomyces cerevisae) and two case studies of hybridization in swordtail fish and Populus trees.

Finding common ground: Toward comparable indicators of adaptive capacity of tree species to a changing climate

Adaptive capacity, one of the three determinants of vulnerability to climate change, is defined as the capacity of species to persist in their current location by coping with novel environmental conditions through acclimation and/or evolution. Although studies have identified indicators of adaptive capacity, few have assessed this capacity in a quantitative way that is comparable across tree species. Yet, such multispecies assessments are needed by forest management and conservation programs to refine vulnerability assessments and to guide the choice of adaptation measures. In this paper, we propose a framework to quantitatively evaluate five key components of tree adaptive capacity to climate change: individual adaptation through phenotypic plasticity, population phenotypic diversity as influenced by genetic diversity, genetic exchange within populations, genetic exchange between populations, and genetic exchange between species. For each component, we define the main mechanisms that underlie adaptive capacity and present associated metrics that can be used as indices. To illustrate the use of this framework, we evaluate the relative adaptive capacity of 26 northeastern North American tree species using values reported in the literature. Our results show adaptive capacity to be highly variable among species and between components of adaptive capacity, such that no one species ranks consistently across all components. On average, the conifer Picea glauca and the broadleaves Acer rubrum and A. saccharinum show the greatest adaptive capacity among the 26 species we documented, whereas the conifers Picea rubens and Thuja occidentalis, and the broadleaf Ostrya virginiana possess the lowest. We discuss limitations that arise when comparing adaptive capacity among species, including poor data availability and comparability issues in metrics derived from different methods or studies. The breadth of data required for such an assessment exemplifies the multidisciplinary nature of adaptive capacity and the necessity of continued cross-collaboration to better anticipate the impacts of a changing climate.

Range and niche expansion through multiple interspecific hybridization: a genotyping by sequencing analysis of Cherleria (Caryophyllaceae)

BACKGROUND

Cherleria (Caryophyllaceae) is a circumboreal genus that also occurs in the high mountains of the northern hemisphere. In this study, we focus on a clade that diversified in the European High Mountains, which was identified using nuclear ribosomal (nrDNA) sequence data in a previous study. With the nrDNA data, all but one species was monophyletic, with little sequence variation within most species. Here, we use genotyping by sequencing (GBS) data to determine whether the nrDNA data showed the full picture of the evolution in the genomes of these species.

RESULTS

The overall relationships found with the GBS data were congruent with those from the nrDNA study. Most of the species were still monophyletic and many of the same subclades were recovered, including a clade of three narrow endemic species from Greece and a clade of largely calcifuge species. The GBS data provided additional resolution within the two species with the best sampling, C. langii and C. laricifolia, with structure that was congruent with geography. In addition, the GBS data showed significant hybridization between several species, including species whose ranges did not currently overlap.

CONCLUSIONS

The hybridization led us to hypothesize that lineages came in contact on the Balkan Peninsula after they diverged, even when those lineages are no longer present on the Balkan Peninsula. Hybridization may also have helped lineages expand their niches to colonize new substrates and different areas. Not only do genome-wide data provide increased phylogenetic resolution of difficult nodes, they also give evidence for a more complex evolutionary history than what can be depicted by a simple, branching phylogeny.

Transcriptome-wide analysis of introgression-resistant regions reveals genetic divergence genes under positive selection in Populus trichocarpa

Comparing gene expression patterns and genetic polymorphisms between populations is of central importance for understanding the origin and maintenance of biodiversity. Based on population-specific gene expression levels and allele frequency differences, we sought to identify population divergence (PD) genes across the introgression-resistant genomic regions of Populus trichocarpa. Genes containing highly diverged loci [i.e., genetic divergence (GD)] or showing expression divergence (ED) between populations were widely distributed in the genome and substantially enriched in functional categories related to stress responses, disease resistance, timing of flowering, cell cycle regulation, plant growth, and development. Nine genomic regions showing evidence of strong positive selection were overlapped with GD genes, which had significant differences between Oregon (a southernmost peripheral deme) and the other demes. However, we did not find evidence that genes under positive selection show an enrichment for ED. PD genes and genes under selection pertained to the same gene classes, such as SERINE/CYSTEINE PROTEASE, ABC TRANSPORTER, GLYCOSYLTRANSFERASE and other transferases. Our analysis also revealed that GD genes were polymorphic within the species (41.9 ± 3.66 biallelic variants per gene), as previously reported in herbaceous plants. By contrast, ED genes contained less genetic variants (10.73 ± 1.14) and were likely highly expressed. In addition, we found that trans- rather than cis-acting variants considerably contribute to the evolution of >90% PD genes. Overall, this study elucidates that cohorts of PD genes agree with the general attributes of known speciation genes and GD genes will provide substrates for positive selection to operate on.

Interspecific introgression and natural selection in the evolution of Japanese apricot (Prunus mume)

Domestication and population differentiation in crops involve considerable phenotypic changes. The logs of these evolutionary paths, including natural/artificial selection, can be found in the genomes of the current populations. However, these profiles have been little studied in tree crops, which have specific characters, such as long generation time and clonal propagation, maintaining high levels of heterozygosity. We conducted exon-targeted resequencing of 129 genomes in the genus Prunus, mainly Japanese apricot (Prunus mume), and apricot (Prunus armeniaca), plum (Prunus salicina), and peach (Prunus persica). Based on their genome-wide single-nucleotide polymorphisms merged with published resequencing data of 79 Chinese P. mume cultivars, we inferred complete and ongoing population differentiation in P. mume. Sliding window characterization of the indexes for genetic differentiation identified interspecific fragment introgressions between P. mume and related species (plum and apricot). These regions often exhibited strong selective sweeps formed in the paths of establishment or formation of substructures of P. mume, suggesting that P. mume has frequently imported advantageous genes from other species in the subgenus Prunus as adaptive evolution. These findings shed light on the complicated nature of adaptive evolution in a tree crop that has undergone interspecific exchange of genome fragments with natural/artificial selections.

Genomics of Evolutionary Novelty in Hybrids and Polyploids

It has long been recognized that hybridization and polyploidy are prominent processes in plant evolution. Although classically recognized as significant in speciation and adaptation, recognition of the importance of interspecific gene flow has dramatically increased during the genomics era, concomitant with an unending flood of empirical examples, with or without genome doubling. Interspecific gene flow is thus increasingly thought to lead to evolutionary innovation and diversification, via adaptive introgression, homoploid hybrid speciation and allopolyploid speciation. Less well understood, however, are the suite of genetic and genomic mechanisms set in motion by the merger of differentiated genomes, and the temporal scale over which recombinational complexity mediated by gene flow might be expressed and exposed to natural selection. We focus on these issues here, considering the types of molecular genetic and genomic processes that might be set in motion by the saltational event of genome merger between two diverged species, either with or without genome doubling, and how these various processes can contribute to novel phenotypes. Genetic mechanisms include the infusion of new alleles and the genesis of novel structural variation including translocations and inversions, homoeologous exchanges, transposable element mobilization and novel insertional effects, presence-absence variation and copy number variation. Polyploidy generates massive transcriptomic and regulatory alteration, presumably set in motion by disrupted stoichiometries of regulatory factors, small RNAs and other genome interactions that cascade from single-gene expression change up through entire networks of transformed regulatory modules. We highlight both these novel combinatorial possibilities and the range of temporal scales over which such complexity might be generated, and thus exposed to natural selection and drift.

Environment-dependent introgression from Quercus dentata to a coastal ecotype of Quercus mongolica var. crispula in northern Japan

Introgression from one species in a specific environment to another may facilitate colonization of the environment by the recipient species. However, such environment-dependent introgression has been clarified in limited plant taxa. In northern Japan, there are two interfertile oak species: Quercus dentata (Qd) in coastal areas and Q. mongolica var. crispula (Qc) in inland areas. However, at higher latitudes where Qd is rare, a coastal Qc ecotype with Qd-like traits is distributed in the coastal areas. We distinguished inland Qc, coastal Qc, and coastal Qd populations based on genome-wide genotypes and multitrait phenotypes and verified introgression from coastal Qd to coastal Qc using reduced library sequencing. Genotypes and phenotypes differed among the populations, and coastal Qc was intermediate between inland Qc and coastal Qd. The ABBA-BABA test showed introgression from coastal Qd to coastal Qc. In coastal Qc, we found various stages of introgression after the first generation of backcross but detected no genomic regions where introgression was enhanced. Overall, we show evidence for introgression from a coastal species to an ecotype of an inland species, which has colonized the coastal environment. It remains unclear whether introgressed alleles are selected in the coastal environment.

Phylogenomics of the genus Populus reveals extensive interspecific gene flow and balancing selection

Phylogenetic analysis is complicated by interspecific gene flow and the presence of shared ancestral polymorphisms, particularly those maintained by balancing selection. In this study, we aimed to examine the prevalence of these factors during the diversification of Populus, a model tree genus in the Northern Hemisphere. We constructed phylogenetic trees of 29 Populus taxa using 80 individuals based on re-sequenced genomes. Our species tree analyses recovered four main clades in the genus based on consensus nuclear phylogenies, but in conflict with the plastome phylogeny. A few interspecific relationships remained unresolved within the multiple-species clade because of inconsistent gene trees. Our results indicated that gene flow has been widespread within each clade and also occurred among the four clades during their early divergence. We identified 45 candidate genes with ancient polymorphisms maintained by balancing selection. These genes were mainly associated with mating compatibility, growth and stress resistance. Both gene flow and selection-mediated ancient polymorphisms are prevalent in the genus Populus. These are potentially important contributors to adaptive variation. Our results provide a framework for the diversification of model tree genus that will facilitate future comparative studies.

Forest genomics: Advancing climate adaptation, forest health, productivity, and conservation

Forest ecosystems provide important ecological services and resources, from habitat for biodiversity to the production of environmentally friendly products, and play a key role in the global carbon cycle. Humanity is counting on forests to sequester and store a substantial portion of the anthropogenic carbon dioxide produced globally. However, the unprecedented rate of climate change, deforestation, and accidental importation of invasive insects and diseases are threatening the health and productivity of forests, and their capacity to provide these services. Knowledge of genetic diversity, local adaptation, and genetic control of key traits is required to predict the adaptive capacity of tree populations, inform forest management and conservation decisions, and improve breeding for productive trees that will withstand the challenges of the 21st century. Genomic approaches have well accelerated the generation of knowledge of the genetic and evolutionary underpinnings of nonmodel tree species, and advanced their applications to address these challenges. This special issue of Evolutionary Applications features 14 papers that demonstrate the value of a wide range of genomic approaches that can be used to better understand the biology of forest trees, including species that are widespread and managed for timber production, and others that are threatened or endangered, or serve important ecological roles. We highlight some of the major advances, ranging from understanding the evolution of genomes since the period when gymnosperms separated from angiosperms 300 million years ago to using genomic selection to accelerate breeding for tree health and productivity. We also discuss some of the challenges and future directions for applying genomic tools to address long-standing questions about forest trees.

Admixture mapping in interspecific Populus hybrids identifies classes of genomic architectures for phytochemical, morphological and growth traits

The genomic architecture of functionally important traits is key to understanding the maintenance of reproductive barriers and trait differences when divergent populations or species hybridize. We conducted a genome-wide association study (GWAS) to study trait architecture in natural hybrids of two ecologically divergent Populus species. We genotyped 472 seedlings from a natural hybrid zone of Populus alba and Populus tremula for genome-wide markers from reduced representation sequencing, phenotyped the plants in common gardens for 46 phytochemical (phenylpropanoid), morphological and growth traits, and used a Bayesian polygenic model for mapping. We detected three classes of genomic architectures: traits with finite, detectable associations of genetic loci with phenotypic variation in addition to highly polygenic heritability; traits with indications for polygenic heritability only; and traits with no detectable heritability. For the first class, we identified genome regions with plausible candidate genes for phenylpropanoid biosynthesis or its regulation, including MYB transcription factors and glycosyl transferases. GWAS in natural, recombinant hybrids represent a promising step towards resolving the genomic architecture of phenotypic traits in long-lived species. This facilitates the fine-mapping and subsequent functional characterization of genes and networks causing differences in hybrid performance and fitness.

Ancient introgression drives adaptation to cooler and drier mountain habitats in a cypress species complex

Introgression may act as an important source of new genetic variation to facilitate the adaptation of organisms to new environments, yet how introgression might enable tree species to adapt to higher latitudes and elevations remains unclear. Applying whole-transcriptome sequencing and population genetic analyses, we present an example of ancient introgression from a cypress species (Cupressus gigantea) that occurs at higher latitude and elevation on the Qinghai-Tibet Plateau into a related species (C. duclouxiana), which has likely aided the latter species to extend its range by colonizing cooler and drier mountain habitats during postglacial periods. We show that 16 introgressed candidate adaptive loci could have played pivotal roles in response to diverse stresses experienced in a high-elevation environment. Our findings provide new insights into the evolutionary history of Qinghai-Tibet Plateau plants and the importance of introgression in the adaptation of species to climate change.

The role of interspecific hybridization in adaptive potential at range margins

Is interspecific hybridization an ordinary part of species biology? And if so, how evolutionarily important is it? These questions have been discussed in the botanical literature, in one form or another, at least since J.P. Lotsy early in the last century. He coined the term syngameon, now defined as "a group of otherwise distinct species interconnected by limited gene exchange, i.e. the most inclusive interbreeding evolutionary unit" (Suarez-Gonzalez, Lexer, & Cronk, Biology Letters, 14, 20170688, ). North American poplars (Populus, Salicaceae) form one such syngameon. In this issue of Molecular Ecology, a new study (Chhatre, Evans, DiFazio, & Keller, Molecular Ecology, 27, ) uses three species from the North American poplar syngameon to tackle the twin issues of (a) the extent of gene exchange and (b) the significance of this gene exchange to the biology of these trees. They demonstrate that a hybrid zone exists where the ranges of Populus angustifolia and Populus balsamifera overlap in the Rocky Mountains, and postulate that this hybridization may facilitate population survival at the range edges. Indeed, the authors show that a remarkable number of loci are introgressing under selection. Very remarkably, they detect additional hybridity (making a trihybrid zone) with Populus trichocarpa (a species that does not occur in the area). Intriguingly, there is some genomic evidence of ancient introgression events. This suggests a model of episodic species divergence and hybridization, in which the syngameon is dynamic and behaving as a supraspecific metapopulation over geological time.

Adaptive introgression: a plant perspective

Introgression is emerging as an important source of novel genetic variation, alongside standing variation and mutation. It is adaptive when such introgressed alleles are maintained by natural selection. Recently, there has been an explosion in the number of studies on adaptive introgression. In this review, we take a plant perspective centred on four lines of evidence: (i) introgression, (ii) selection, (iii) phenotype and (iv) fitness. While advances in genomics have contributed to our understanding of introgression and porous species boundaries (task 1), and the detection of signatures of selection in introgression (task 2), the investigation of adaptive introgression critically requires links to phenotypic variation and fitness (tasks 3 and 4). We also discuss the conservation implications of adaptive introgression in the face of climate change. Adaptive introgression is particularly important in rapidly changing environments, when standing genetic variation and mutation alone may only offer limited potential for adaptation. We conclude that clarifying the magnitude and fitness effects of introgression with improved statistical techniques, coupled with phenotypic evidence, has great potential for conservation and management efforts.

From transects to transcripts: Teasing apart the architecture of reproductive isolation

Understanding the processes underlying speciation has long been a challenge to evolutionary biologists. This spurs from difficulties teasing apart the various mechanisms that contribute to the evolution of barriers to reproduction. The study by Rafati et al. () in this issue of Molecular Ecology combines spatially explicit whole-genome resequencing with evaluation of differential gene expression across individuals with mixed ancestry to associate the genomic architecture of reproductive barriers with expression of reproductive incompatibilities. In a natural hybrid zone between rabbit subspecies, Oryctolagus cuniculus cuniculus and O. c. algirus (Figure ), Rafati et al. () use landscape-level patterns of allele frequency variation to identify potential candidate regions of the genome associated with reproductive isolation. These candidate regions are used to test predictions associated with the genomic architecture of reproductive barriers, including the role of structural rearrangements, enrichment of functional categories associated with incompatibilities, and the contribution of protein-coding versus regulatory changes. A lack of structural rearrangements and limited protein-coding changes in candidate regions point towards the importance of regulatory variation as major contributors to genetic incompatibilities, while functional enrichments indicate overrepresentation of genes associated with male infertility. To quantify phenotypic expression of proposed incompatibilities, the authors assess gene expression of experimental crosses. Extensive misregulation of gene expression within the testes of backcross hybrids relative to F1 and parental individuals provides an important link between genotype and phenotype, validating hypotheses developed from assessment of genomic architectures. Together, this work shows how pairing natural hybrid zones with experimental crosses can be used to link observations in nature to mechanistic underpinnings that may be tested experimentally.

Functional and structural profiles of GST gene family from three Populus species reveal the sequence-function decoupling of orthologous genes

A common assumption in comparative genomics is that orthologous genes are functionally more similar than paralogous genes. However, the validity of this assumption needs to be assessed using robust experimental data. We conducted tissue-specific gene expression and protein function analyses of orthologous groups within the glutathione S-transferase (GST) gene family in three closely related Populus species: Populus trichocarpa, Populus euphratica and Populus yatungensis. This study identified 21 GST orthologous groups in the three Populus species. Although the sequences of the GST orthologous groups were highly conserved, the divergence in enzymatic functions was prevalent. Through site-directed mutagenesis of orthologous proteins, this study revealed that nonsynonymous substitutions at key amino acid sites played an important role in the divergence of enzymatic functions. In particular, a single amino acid mutation (Arg39→Trp39) contributed to P. euphratica PeGSTU30 possessing high enzymatic activity via increasing the hydrophobicity of the active cavity. This study provided experimental evidence showing that orthologues belonging to the gene family have functional divergences. The nonsynonymous substitutions at a few amino acid sites resulted in functional divergence of the orthologous genes. Our findings provide new insights into the evolution of orthologous genes in closely related species.

Population genomic analysis reveals that homoploid hybrid speciation can be a lengthy process

An increasing number of species are thought to have originated by homoploid hybrid speciation (HHS), but in only a handful of cases are details of the process known. A previous study indicated that Picea purpurea, a conifer in the Qinghai-Tibet Plateau (QTP), originated through HHS from P. likiangensis and P. wilsonii. To investigate this origin in more detail, we analysed transcriptome data for 114 individuals collected from 34 populations of the three Picea species from their core distributions in the QTP. Phylogenetic, principal component and admixture analyses of nuclear SNPs showed the species to be delimited genetically and that P. purpurea was admixed with approximately 60% of its ancestry derived from P. wilsonii and 40% from P. likiangensis. Coalescent simulations revealed the best-fitting model of origin involved formation of an intermediate hybrid lineage between P. likiangensis and P. wilsonii approximately 6 million years ago (mya), which backcrossed to P. wilsonii to form P. purpurea approximately one mya. The intermediate hybrid lineage no longer exists and is referred to as a "ghost" lineage. Our study emphasizes the power of population genomic analysis combined with coalescent analysis for reconstructing the stages involved in the origin of a homoploid hybrid species over an extended period. In contrast to other studies, we show that these stages can in some instances span a relatively long period of evolutionary time.

The role of diversification in community assembly of the oaks (Quercus L.) across the continental U.S

PREMISE OF THE STUDY

Evolutionary and biogeographic history, including past environmental change and diversification processes, are likely to have influenced the expansion, migration, and extinction of populations, creating evolutionary legacy effects that influence regional species pools and the composition of communities. We consider the consequences of the diversification process in shaping trait evolution and assembly of oak-dominated communities throughout the continental United States (U.S.).

METHODS

Within the U.S. oaks, we tested for phylogenetic and functional trait patterns at different spatial scales, taking advantage of a dated phylogenomic analysis of American oaks and the U.S. Forest Service (USFS) Forest Inventory and Analysis (FIA).

KEY RESULTS

We find (1) phylogenetic overdispersion at small grain sizes throughout the U.S. across all spatial extents and (2) a shift from overdispersion to clustering with increasing grain sizes. Leaf traits have evolved in a convergent manner, and these traits are clustered in communities at all spatial scales, except in the far west, where species with contrasting leaf types co-occur.

CONCLUSIONS

Our results support the hypotheses that (1) interspecific interactions were important in parallel adaptive radiation of the genus into a range of habitats across the continent and (2) that the diversification process is a critical driver of community assembly. Functional convergence of complementary species from distinct clades adapted to the same local habitats is a likely mechanism that allows distantly related species to coexist. Our findings contribute to an explanation of the long-term maintenance of high oak diversity and the dominance of the oak genus in North America.