The emerging importance of cross-ploidy hybridisation and introgression.

Natural hybridisation is now recognised as pervasive in its occurrence across the Tree of Life. Resurgent interest in natural hybridisation fuelled by developments in genomics has led to an improved understanding of the genetic factors that promote or prevent species cross-mating. Despite this body of work overturning many widely held assumptions about the genetic barriers to hybridisation, it is still widely thought that ploidy differences between species will be an absolute barrier to hybridisation and introgression. Here, we revisit this assumption, reviewing findings from surveys of polyploidy and hybridisation in the wild. In a case study in the British flora, 203 hybrids representing 35% of hybrids with suitable data have formed via cross-ploidy matings, while a wider literature search revealed 59 studies (56 in plants and 3 in animals) in which cross-ploidy hybridisation has been confirmed with genetic data. These results show cross-ploidy hybridisation is readily overlooked, and potentially common in some groups. General findings from these studies include strong directionality of hybridisation, with introgression usually towards the higher ploidy parent, and cross-ploidy hybridisation being more likely to involve allopolyploids than autopolyploids. Evidence for adaptive introgression across a ploidy barrier and cases of cross-ploidy hybrid speciation shows the potential for important evolutionary outcomes.

Evolutionary origin of genomic structural variations in domestic yaks.

Yak has been subject to natural selection, human domestication and interspecific introgression during its evolution. However, genetic variants favored by each of these processes have not been distinguished previously. We constructed a graph-genome for 47 genomes of 7 cross-fertile bovine species. This allowed detection of 57,432 high-resolution structural variants (SVs) within and across the species, which were genotyped in 386 individuals. We distinguished the evolutionary origins of diverse SVs in domestic yaks by phylogenetic analyses. We further identified 334 genes overlapping with SVs in domestic yaks that bore potential signals of selection from wild yaks, plus an additional 686 genes introgressed from cattle. Nearly 90% of the domestic yaks were introgressed by cattle. Introgression of an SV spanning the KIT gene triggered the breeding of white domestic yaks. We validated a significant association of the selected stratified SVs with gene expression, which contributes to phenotypic variations. Our results highlight that SVs of different origins contribute to the phenotypic diversity of domestic yaks.

Abiotic Niche Divergence of Hybrid Species from Their Progenitors.

AbstractAlthough more frequently discussed recently than previously, the role of ecology in homoploid hybrid and allopolyploid speciation has not been subjected to comparative analysis. We examined abiotic niche divergence of 22 assumed homoploid hybrid species and 60 allopolyploid species from that of their progenitors. Ecological niche modeling was employed in an analysis of each species' fundamental niche, and ordination methods were used in an analysis of realized niches. Both analyses utilized 100,000 georeferenced records. From estimates of niche overlap and niche breadth, we identified for both types of hybrid species four niche divergence patterns: niche novelty, niche contraction, niche intermediacy, and niche expansion. Niche shifts involving niche novelty were common and considered likely to play an important role in the establishment of both types of hybrid species, although more so for homoploid hybrid species than for allopolyploid species. Approximately 70% of homoploid hybrid species versus 37% of allopolyploid species showed shifts in the fundamental niche from their parents, and ∼86% versus ∼52%, respectively, exhibited shifts in the realized niche. Climate was shown to contribute more than soil and landform to niche shifts in both types of hybrid species. Overall, our results highlight the significance of abiotic niche divergence for hybrid speciation, especially without genome duplication.

Molecular signatures of parallel adaptive divergence causing reproductive isolation and speciation across two genera.

Parallel evolution of reproductive isolation (PERI) provides strong evidence for natural selection playing a fundamental role in the origin of species. However, PERI has been rarely demonstrated for well established species drawn from different genera. In particular, parallel molecular signatures for the same genes in response to similar habitat divergence in such different lineages is lacking. Here, based on whole-genome sequencing data, we first explore the speciation process in two sister species of Carpinus (Betulaceae) in response to divergence for temperature and soil-iron concentration in habitats they occupy in northern and southwestern China, respectively. We then determine whether parallel molecular mutations occur during speciation in this pair of species and also in another sister-species pair of the related genus, Ostryopsis, which occupy similarly divergent habitats in China. We show that gene flow occurred during the origin of both pairs of sister species since approximately 9.8 or approximately 2 million years ago, implying strong natural selection during divergence. Also, in both species pairs we detected concurrent positive selection in a gene (LHY) for flowering time and in two paralogous genes (FRO4 and FRO7) of a gene family known to be important for iron tolerance. These changes were in addition to changes in other major genes related to these two traits. The different alleles of these particular candidate genes possessed by the sister species of Carpinus were functionally tested and indicated likely to alter flowering time and iron tolerance as previously demonstrated in the pair of Ostryopsis sister species. Allelic changes in these genes may have effectively resulted in high levels of prezygotic reproductive isolation to evolve between sister species of each pair. Our results show that PERI can occur in different genera at different timescales and involve similar signatures of molecular evolution at genes or paralogues of the same gene family, causing reproductive isolation as a consequence of adaptation to similarly divergent habitats.

Occurrence and Prevention of Delayed Autonomous Selfing in Salvia umbratica (Lamiaceae).

Delayed autonomous selfing (DAS) provides reproductive assurance under conditions of pollinator and/or pollen-limitation. Few plant species have been investigated to determine if DAS is terminated when a flower is sufficiently pollinated by a pollen vector, thereby saving plant resources for other purposes. We examined this possibility in bumblebee-pollinated Salvia umbratica. We first showed that DAS resulting in high fruit set (100%) and seed set (>80%) per flower occurred in the absence of insect pollinators by means of style recurvature and was completed in 94% of flowers 72 h after they opened. In contrast, in flowers pollinated immediately after opening, DAS was prevented by corollas dropping away before styles recurve toward the upper thecae. We next showed that hand-pollination of flowers immediately after they opened resulted in high fruit set (100%) and seed set (>80%) when 5-10 pollen grains or more were deposited on their stigmas, whereas fruit set and seed set were reduced to 45.00 and 22.50%, respectively, when pollen loads were reduced to 1-3 pollen grains. Finally, we showed that on average single pollinator visits deposited 26 pollen grains on stigmas of flowers that had just opened, which is more than enough to ensure high fruit and seed set. Our results indicate that flower longevity is highly correlated with the pollinator environment and female fitness of S. umbratica, with extended flower longevity allowing DAS to occur being advantageous when pollinators are absent, while reduced floral longevity and prevention of DAS being favored when flowers are pollinated by pollinators. Thus, flower longevity in S. umbratica varies so as to optimize reproductive output and resource efforts, and is dependent on the availability and effectiveness of pollinators to pollinate flowers.

Hybrid speciation via inheritance of alternate alleles of parental isolating genes.

It is increasingly realized that homoploid hybrid speciation (HHS), which involves no change in chromosome number, is an important mechanism of speciation. HHS will likely increase in frequency as ecological and geographical barriers between species are continuing to be disrupted by human activities. HHS requires the establishment of reproductive isolation between a hybrid and its parents, but the underlying genes and genetic mechanisms remain largely unknown. In this study, we reveal by integrated approaches that reproductive isolation originates in one homoploid hybrid plant species through the inheritance of alternate alleles at genes that determine parental premating isolation. The parent species of this hybrid species are reproductively isolated by differences in flowering time and survivorship on soils containing high concentrations of iron. We found that the hybrid species inherits alleles of parental isolating major genes related to flowering time from one parent and alleles of major genes related to iron tolerance from the other parent. In this way, it became reproductively isolated from one parent by the difference in flowering time and from the other by habitat adaptation (iron tolerance). These findings and further modeling results suggest that HHS may occur relatively easily via the inheritance of alternate parental premating isolating genes and barriers.

Senecio as a model system for integrating studies of genotype, phenotype and fitness.

Two major developments have made it possible to use examples of ecological radiations as model systems to understand evolution and ecology. First, the integration of quantitative genetics with ecological experiments allows detailed connections to be made between genotype, phenotype, and fitness in the field. Second, dramatic advances in molecular genetics have created new possibilities for integrating field and laboratory experiments with detailed genetic sequencing. Combining these approaches allows evolutionary biologists to better study the interplay between genotype, phenotype, and fitness to explore a wide range of evolutionary processes. Here, we present the genus Senecio (Asteraceae) as an excellent system to integrate these developments, and to address fundamental questions in ecology and evolution. Senecio is one of the largest and most phenotypically diverse genera of flowering plants, containing species ranging from woody perennials to herbaceous annuals. These Senecio species exhibit many growth habits, life histories, and morphologies, and they occupy a multitude of environments. Common within the genus are species that have hybridized naturally, undergone polyploidization, and colonized diverse environments, often through rapid phenotypic divergence and adaptive radiation. These diverse experimental attributes make Senecio an attractive model system in which to address a broad range of questions in evolution and ecology.

Genome-wide analysis of Cushion willow provides insights into alpine plant divergence in a biodiversity hotspot.

The Hengduan Mountains (HDM) biodiversity hotspot exhibits exceptional alpine plant diversity. Here, we investigate factors driving intraspecific divergence within a HDM alpine species Salix brachista (Cushion willow), a common component of subnival assemblages. We produce a high-quality genome assembly for this species and characterize its genetic diversity, population structure and pattern of evolution by resequencing individuals collected across its distribution. We detect population divergence that has been shaped by a landscape of isolated sky island-like habitats displaying strong environmental heterogeneity across elevational gradients, combined with population size fluctuations that have occurred since approximately the late Miocene. These factors are likely important drivers of intraspecific divergence within Cushion willow and possibly other alpine plants with a similar distribution. Since intraspecific divergence is often the first step toward speciation, the same factors can be important contributors to the high alpine species diversity in the HDM.

Late Pleistocene speciation of three closely related tree peonies endemic to the Qinling-Daba Mountains, a major glacial refugium in Central China.

Determining the factors promoting speciation is a major task in ecological and evolutionary research and can be aided by phylogeographic analysis. The Qinling-Daba Mountains (QDM) located in central China form an important geographic barrier between southern subtropical and northern temperate regions, and exhibit complex topography, climatic, and ecological diversity. Surprisingly, few phylogeographic analyses and studies of plant speciation in this region have been conducted. To address this issue, we investigated the genetic divergence and evolutionary histories of three closely related tree peony species (Paeonia qiui, P. jishanensis, and P. rockii) endemic to the QDM. Forty populations of the three tree peony species were genotyped using 22 nuclear simple sequence repeat markers (nSSRs) and three chloroplast DNA sequences to assess genetic structure and phylogenetic relationships, supplemented by morphological characterization and ecological niche modeling (ENM). Morphological and molecular genetic analyses showed the three species to be clearly differentiated from each other. In addition, coalescent analyses using DIYABC conducted on nSSR variation indicated that the species diverged from each other in the late Pleistocene, while ecological niche modeling (ENM) suggested they occupied a larger area during the Last Glacial Maximum (LGM) than at present. The combined genetic evidence from nuclear and chloroplast DNA and the results of ENM indicate that each species persisted through the late Pleistocene in multiple refugia in the Qinling, Daba, and Taihang Mountains with divergence favored by restricted gene flow caused by geographic isolation, ecological divergence, and limited pollen and seed dispersal. Our study contributes to a growing understanding of the origin and population structure of tree peonies and provides insights into the high level of plant endemism present in the Qinling-Daba Mountains of Central China.

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.

Completing the hybridization triangle: the inheritance of genetic incompatibilities during homoploid hybrid speciation in ragworts (Senecio).

A new homoploid hybrid lineage needs to establish a degree of reproductive isolation from its parent species if it is to persist as an independent entity, but the role hybridization plays in this process is known in only a handful of cases. The homoploid hybrid ragwort species, Senecio squalidus (Oxford ragwort), originated following the introduction of hybrid plants to the UK approximately 320 years ago. The source of the hybrid plants was from a naturally occurring hybrid zone between S. aethnensis and S. chrysanthemifolius on Mount Etna, Sicily. Previous studies of the parent species found evidence for multiple incompatibility loci causing transmission ratio distortion of genetic markers in their hybrid progeny. This study closes the hybridization triangle by reporting a genetic mapping analysis of the remaining two paired cross combinations between S. squalidus and its parents. Genetic maps produced from F2 mapping families were generally collinear but with half of the linkage groups showing evidence of genomic reorganization between genetic maps. The new maps produced from crosses between S. squalidus and each parent showed multiple incompatibility loci distributed across the genome, some of which co-locate with previously reported incompatibility loci between the parents. These findings suggest that this young homoploid hybrid species has inherited a unique combination of genomic rearrangements and incompatibilities from its parents that contribute to its reproductive isolation.

Reticulate evolution within a spruce (Picea) species complex revealed by population genomic analysis.

The role of reticulation in the rapid diversification of organisms is attracting greater attention in evolutionary biology. Evidence of genetic exchange between diverging taxa is reported frequently, although most studies fail to show how hybridization and introgression contribute to the adaptation and differentiation of introgressed taxa. Here, we report a population genomics approach to test the role of hybridization and introgression in the evolution of the Picea likiangensis species complex, which comprises four taxa occurring in the biodiversity hotspot of the Hengduan-Himalayan mountains. Based on 84,793 SNPs detected in transcriptomes of 82 trees collected from 35 localities, we identified 18 hybrids (including backcrosses) distributed within the range boundaries of the four taxa. Coalescent simulations, for each pair of taxa and for all taxa taken together, rejected several tree-like divergence models and supported instead a reticulate evolution model with secondary contacts occurring during Pleistocene glacial cycles after initial divergence in the late Pliocene. Significant gene flow occurred among some taxa after secondary contact according to an analysis based on modified ABBA-BABA statistics that accommodated a rapid diversification scenario. A novel finding was that introgression between certain taxa can contribute to increasing divergence (and possibly reproductive isolation) between those taxa and other taxa within a complex at some loci. These results illuminate the reticulate nature of evolution within the P. likiangensis complex and highlight the value of population genomic data in detecting the effects of introgression in the rapid diversification of related taxa.

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.

Demographic expansion and genetic load of the halophyte model plant Eutrema salsugineum.

The halophyte model plant Eutrema salsugineum (Brassicaceae) disjunctly occurs in temperate to subarctic Asia and North America. This vast, yet extremely discontinuous distribution constitutes an ideal system to examine long-distance dispersal and the ensuing accumulation of deleterious mutations as expected in expanding populations of selfing plants. In this study, we resequenced individuals from 23 populations across the range of E. salsugineum. Our population genomic data indicate that E. salsugineum migrated "out of the Altai region" at least three times to colonize northern China, northeast Russia and western China. It then expanded its distribution into North America independently from northeast Russia and northern China, respectively. The species colonized northern China around 33.7 thousand years ago (kya) and underwent a considerable expansion in range size approximately 7-8 kya. The western China lineage is likely a hybrid derivative of the northern China and Altai lineages, originating approximately 25-30 kya. Deleterious alleles accumulated in a stepwise manner from (a) Altai to northern China and North America and (b) Altai to northeast Russia and North America. In summary, E. salsugineum dispersed from Asia to North America and deleterious mutations accumulated in a stepwise manner during the expansion of the species' distribution.

Ancient polymorphisms and divergence hitchhiking contribute to genomic islands of divergence within a poplar species complex.

How genome divergence eventually leads to speciation is a topic of prime evolutionary interest. Genomic islands of elevated divergence are frequently reported between diverging lineages, and their size is expected to increase with time and gene flow under the speciation-with-gene-flow model. However, such islands can also result from divergent sorting of ancient polymorphisms, recent ecological selection regardless of gene flow, and/or recurrent background selection and selective sweeps in low-recombination regions. It is challenging to disentangle these nonexclusive alternatives, but here we attempt to do this in an analysis of what drove genomic divergence between four lineages comprising a species complex of desert poplar trees. Within this complex we found that two morphologically delimited species, Populus euphratica and Populus pruinosa, were paraphyletic while the four lineages exhibited contrasting levels of gene flow and divergence times, providing a good system for testing hypotheses on the origin of divergence islands. We show that the size and number of genomic islands that distinguish lineages are not associated with either rate of recent gene flow or time of divergence. Instead, they are most likely derived from divergent sorting of ancient polymorphisms and divergence hitchhiking. We found that highly diverged genes under lineage-specific selection and putatively involved in ecological and morphological divergence occur both within and outside these islands. Our results highlight the need to incorporate demography, absolute divergence measurement, and gene flow rate to explain the formation of genomic islands and to identify potential genomic regions involved in speciation.

Chasing ghosts: allopolyploid origin of Oxyria sinensis (Polygonaceae) from its only diploid congener and an unknown ancestor.

Reconstructing the origin of a polyploid species is particularly challenging when an ancestor has become extinct. Under such circumstances, the extinct donor of a genome found in the polyploid may be treated as a 'ghost' species in that its prior existence is recognized through the presence of its genome in the polyploid. In this study, we aimed to determine the polyploid origin of Oxyria sinensis (2n = 40) for which only one congeneric species is known, that is diploid O. digyna (2n = 14). Genomic in situ hybridization (GISH), transcriptome, phylogenetic and demographic analyses, and ecological niche modelling were conducted for this purpose. GISH revealed that O. sinensis comprised 14 chromosomes from O. digyna and 26 chromosomes from an unknown ancestor. Transcriptome analysis indicated that following divergence from O. digyna, involving genome duplication around 12 million years ago (Ma), a second genome duplication occurred approximately 6 Ma to give rise to O. sinensis. Oxyria sinensis was shown to contain homologous gene sequences divergent from those present in O. digyna in addition to a set that clustered with those in O. digyna. Coalescent simulations indicated that O. sinensis expanded its distribution approximately 6-7 Ma, possibly following the second polyploidization event, whereas O. digyna expanded its range much later. It was also indicated that the distributions of both species contracted and re-expanded during the Pleistocene climatic oscillations. Ecological niche modelling similarly suggested that both species experienced changes in their distributional ranges in response to Quaternary climatic changes. The extinction of the unknown 'ghost' tetraploid species implicated in the origin of O. sinensis could have resulted from superior adaptation of O. sinensis to repeated climatic changes in the region where it now occurs.

Morphological Convergence Between an Allopolyploid and One of its Parental Species Correlates with Biased Gene Expression and DNA Loss.

The contribution of gene expression modulation to phenotypic evolution is of major importance to an understanding of the origin of divergent or convergent phenotypes during and following polyploid speciation. Here, we analyzed genome-wide gene expression in 2 subspecies of the allotetraploid species, Senecio mohavensis A. Gray, and its diploid parents S. flavus (Decne.) Sch. Bip. and S. glaucus L. The tetraploid is morphologically much more similar to S. flavus, leading to earlier confusion over its taxonomic status. By means of an analysis of transcriptomes of all 3 species, we show that gene expression divergence between the parent species is relatively low (ca. 14% of loci), whereas there is significant unequal expression between ca. 20-25% of the parental homoeologues (gene copies) in the tetraploid. The majority of the expression bias in the tetraploid is in favor of S. flavus homoeologues (ca. 65% of the differentially expressed loci), and overall expression of this parental species subgenome is higher than that of the S. glaucus subgenome. To determine whether absence of expression of a particular S. glaucus homoeologue in the allotetraploid could be due to loss of DNA, we carried out a PCR-based assay and confirmed that in 3 out of 10 loci the S. glaucus homoeologue appeared absent. Our results suggest that biased gene expression is one cause of the allotetraploid S. mohavensis being more similar in morphology to one of its parent, S. flavus, and that such bias could result, in part, from loss of S. glaucus homoeologues at some loci in the allotetraploid.