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.

The genome sequence of common ivy, Hedera helix L., 1753.

We present a genome assembly from a specimen of Hedera helix (common ivy; Streptophyta; Magnoliopsida; Apiales; Araliaceae). The genome sequence is 1,199.4 megabases in span. Most of the assembly is scaffolded into 24 chromosomal pseudomolecules. The mitochondrial and plastid genomes have also been assembled and are 609.2 and 162.2 kilobases in length respectively.

Hybridization and divergent climatic preferences drive divergence of two allopatric Gentiana species on the Qinghai-Tibet Plateau.

BACKGROUND AND AIMS: Exploring how species diverge is vital for understanding the drivers of speciation. Factors such as geographical separation and ecological selection, hybridization, polyploidization and shifts in mating system are all major mechanisms of plant speciation, but their contributions to divergence are rarely well understood. Here we test these mechanisms in two plant species, Gentiana lhassica and G. hoae, with the goal of understanding recent allopatric species divergence on the Qinghai-Tibet Plateau (QTP). METHODS: We performed Bayesian clustering, phylogenetic analysis and estimates of hybridization using 561 302 nuclear genomic single nucleotide polymorphisms (SNPs). We performed redundancy analysis, and identified and annotated species-specific SNPs (ssSNPs) to explore the association between climatic preference and genetic divergence. We also estimated genome sizes using flow cytometry to test for overlooked polyploidy. KEY RESULTS: Genomic evidence confirms that G. lhassica and G. hoae are closely related but distinct species, while genome size estimates show divergence occurred without polyploidy. Gentiana hoae has significantly higher average FIS values than G. lhassica. Population clustering based on genomic SNPs shows no signature of recent hybridization, but each species is characterized by a distinct history of hybridization with congeners that has shaped genome-wide variation. Gentiana lhassica has captured the chloroplast and experienced introgression with a divergent gentian species, while G. hoae has experienced recurrent hybridization with related taxa. Species distribution modelling suggested range overlap in the Last Interglacial Period, while redundancy analysis showed that precipitation and temperature are the major climatic differences explaining the separation of the species. The species differ by 2993 ssSNPs, with genome annotation showing missense variants in genes involved in stress resistance. CONCLUSIONS: This study suggests that the distinctiveness of these species on the QTP is driven by a combination of hybridization, geographical isolation, mating system differences and evolution of divergent climatic preferences.

Temporal matches between monarch butterfly and milkweed population changes over the past 25,000 years.

In intimate ecological interactions, the interdependency of species may result in correlated demographic histories. For species of conservation concern, understanding the long-term dynamics of such interactions may shed light on the drivers of population decline. Here, we address the demographic history of the monarch butterfly, Danaus plexippus, and its dominant host plant, the common milkweed Asclepias syriaca (A. syriaca), using broad-scale sampling and genomic inference. Because genetic resources for milkweed have lagged behind those for monarchs, we first release a chromosome-level genome assembly and annotation for common milkweed. Next, we show that despite its enormous geographic range across eastern North America, A. syriaca is best characterized as a single, roughly panmictic population. Using approximate Bayesian computation with random forests (ABC-RF), a machine learning method for reconstructing demographic histories, we show that both monarchs and milkweed experienced population expansion during the most recent recession of North American glaciers 10,000-20,000 years ago. Our data also identify concurrent population expansions in both species during the large-scale clearing of eastern forests (∼200 years ago). Finally, we find no evidence that either species experienced a reduction in effective population size over the past 75 years. Thus, the well-documented decline of monarch abundance over the past 40 years is not visible in our genomic dataset, reflecting a possible mismatch of the overwintering census population to effective population size in this species.

Genetic factors predict hybrid formation in the British flora.

Natural hybridization can have a profound evolutionary impact, with consequences ranging from the extinction of rare taxa to the origin of new species. Natural hybridization is particularly common in plants; however, our understanding of the general factors that promote or prevent hybridization is hampered by the highly variable outcomes in different lineages. Here, we quantify the influence of different predictors on hybrid formation across species from an entire flora. We combine estimates of hybridization with ecological attributes and a new species-level phylogeny for over 1,100 UK flowering plant species. Our results show that genetic factors, particularly parental genetic distance, as well as phylogenetic position and ploidy, are key determinants of hybrid formation, whereas many other factors such as range overlap and genus size explain much less variation in hybrid formation. Overall, intrinsic genetic factors shape the evolutionary and ecological consequences of natural hybridization across species in a flora.

Is there hybridization between diploid and tetraploid Euphrasia in a secondary contact zone?

PREMISE: Strong postzygotic reproductive isolating barriers are usually expected to limit the extent of natural hybridization between species with contrasting ploidy. However, genomic sequencing has revealed previously overlooked examples of natural cross-ploidy hybridization in some flowering plant genera, suggesting that the phenomenon may be more common than once thought. We investigated potential cross-ploidy hybridization in British eyebrights (Euphrasia, Orobanchaceae), a group from which 13 putative cross-ploidy hybrid combinations have been reported based on morphology. METHODS: We analyzed a contact zone between diploid Euphrasia rostkoviana and tetraploid E. arctica in Wales. We sequenced part of the internal transcribed spacer (ITS) of nuclear ribosomal DNA and used genotyping by sequencing (GBS) to look for evidence of cross-ploidy hybridization and introgression. RESULTS: Common variant sites in the ITS region were fixed between diploids and tetraploids, indicating a strong barrier to hybridization. Clustering analyses of 356 single-nucleotide polymorphisms (SNPs) generated using GBS clearly separated samples by ploidy and revealed strong genetic structure (FST = 0.44). However, the FST distribution across all SNPs was bimodal, indicating potential differential selection on loci between diploids and tetraploids. Demographic inference suggested potential gene flow, limited to around one or fewer migrants per generation. CONCLUSIONS: Our results suggest that recent cross-ploidy hybridization is rare or absent in a site of secondary contact in Euphrasia. While a strong ploidy barrier prevents hybridization over ecological timescales, such hybrids may form in stable populations over evolutionary timescales, potentially allowing cross-ploidy introgression to take place.

Population genomics reveal deep divergence and strong geographical structure in gentians in the Hengduan Mountains.

Understanding the evolutionary and ecological processes driving population differentiation and speciation can provide critical insights into the formation of biodiversity. Here, we examine the link between population genetic processes and biogeographic history underlying the generation of diversity in the Hengduan Mountains (HM), a region harboring a rich and dynamic flora. We used restriction site-associated DNA sequencing to generate 1,907 single-nucleotide polymorphisms (SNPs) and four-kb of plastid sequence in species of the Gentiana hexaphylla complex (Gentianaceae). We performed genetic clustering with spatial and non-spatial models, phylogenetic reconstructions, and ancestral range estimation, with the aim of addressing the processes influencing diversification of G. hexaphylla in the HM. We find the G. hexaphylla complex is characterized by geographic genetic structure with clusters corresponding to the South, North and the central HM. Phylogenetic reconstruction and pairwise F ST analyses showed deep differentiation between Southern and Northern populations in the HM. The population in Mount Taibai exhibited the highest genetic similarity to the North HM. Ancestral range estimation indicated that the G. hexaphylla complex originated in the central HM and then diverged in the Pliocene and the Early Pleistocene, before dispersing widely, resulting in the current distinct lineages. Overall, we found deep genomic differentiation in the G. hexaphylla complex corresponds to geographic barriers to dispersal in the HM and highlights a critical role of the uplift of the Daxue Mountains and subsequent climatic fluctuations underlying diversification. The colonization of G. hexaphylla in the Mount Taibai region suggests directional dispersal between the alpine flora of the Qinling Mountains and the HM.

Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights (Euphrasia) Revealed by k-mers.

Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.

Pervasive Phylogenomic Incongruence Underlies Evolutionary Relationships in Eyebrights (Euphrasia, Orobanchaceae).

Disentangling the phylogenetic relationships of taxonomically complex plant groups is often mired by challenges associated with recent speciation, hybridization, complex mating systems, and polyploidy. Here, we perform a phylogenomic analysis of eyebrights (Euphrasia), a group renowned for taxonomic complexity, with the aim of documenting the extent of phylogenetic discordance at both deep and at shallow phylogenetic scales. We generate whole-genome sequencing data and integrate this with prior genomic data to perform a comprehensive analysis of nuclear genomic, nuclear ribosomal (nrDNA), and complete plastid genomes from 57 individuals representing 36 Euphrasia species. The species tree analysis of 3,454 conserved nuclear scaffolds (46 Mb) reveals that at shallow phylogenetic scales postglacial colonization of North Western Europe occurred in multiple waves from discrete source populations, with most species not being monophyletic, and instead combining genomic variants from across clades. At a deeper phylogenetic scale, the Euphrasia phylogeny is structured by geography and ploidy, and partially by taxonomy. Comparative analyses show Southern Hemisphere tetraploids include a distinct subgenome indicative of independent polyploidy events from Northern Hemisphere taxa. In contrast to the nuclear genome analyses, the plastid genome phylogeny reveals limited geographic structure, while the nrDNA phylogeny is informative of some geographic and taxonomic affinities but more thorough phylogenetic inference is impeded by the retention of ancestral polymorphisms in the polyploids. Overall our results reveal extensive phylogenetic discordance at both deeper and shallower nodes, with broad-scale geographic structure of genomic variation but a lack of definitive taxonomic signal. This suggests that Euphrasia species either have polytopic origins or are maintained by narrow genomic regions in the face of extensive homogenizing gene flow. Moreover, these results suggest genome skimming will not be an effective extended barcode to identify species in groups such as Euphrasia, or many other postglacial species groups.

The first genome for the Cape Primrose Streptocarpus rexii (Gesneriaceae), a model plant for studying meristem-driven shoot diversity.

Cape Primroses (Streptocarpus, Gesneriaceae) are an ideal study system for investigating the genetics underlying species diversity in angiosperms. Streptocarpus rexii has served as a model species for plant developmental research for over five decades due to its unusual extended meristem activity present in the leaves. In this study, we sequenced and assembled the complete nuclear, chloroplast, and mitochondrial genomes of S. rexii using Oxford Nanopore Technologies long read sequencing. Two flow cells of PromethION sequencing resulted in 32 billion reads and were sufficient to generate a draft assembly including the chloroplast, mitochondrial and nuclear genomes, spanning 776 Mbp. The final nuclear genome assembly contained 5,855 contigs, spanning 766 Mbp of the 929-Mbp haploid genome with an N50 of 3.7 Mbp and an L50 of 57 contigs. Over 70% of the draft genome was identified as repeats. A genome repeat library of Gesneriaceae was generated and used for genome annotation, with a total of 45,045 genes annotated in the S. rexii genome. Ks plots of the paranomes suggested a recent whole genome duplication event, shared between S. rexii and Primulina huaijiensis. A new chloroplast and mitochondrial genome assembly method, based on contig coverage and identification, was developed, and successfully used to assemble both organellar genomes of S. rexii. This method was developed into a pipeline and proved widely applicable. The nuclear genome of S. rexii and other datasets generated and reported here will be invaluable resources for further research to aid in the identification of genes involved in morphological variation underpinning plant diversification.

The symmetry spectrum in a hybridising, tropical group of rhododendrons.

Many diverse plant clades possess bilaterally symmetrical flowers and specialised pollination syndromes, suggesting that these traits may promote diversification. We examined the evolution of diverse floral morphologies in a species-rich tropical radiation of Rhododendron. We used restriction-site associated DNA sequencing on 114 taxa from Rhododendron sect. Schistanthe to reconstruct phylogenetic relationships and examine hybridisation. We then captured and quantified floral variation using geometric morphometric analyses, which we interpreted in a phylogenetic context. We uncovered phylogenetic conflict and uncertainty caused by introgression within and between clades. Morphometric analyses revealed flower symmetry to be a morphological continuum without clear transitions between radial and bilateral symmetry. Tropical Rhododendron species that began diversifying into New Guinea c. 6 million years ago expanded into novel floral morphological space. Our results showed that the evolution of tropical Rhododendron is characterised by recent speciation, recurrent hybridisation and the origin of floral novelty. Floral variation evolved via changes to multiple components of the corolla that are only recognised in geometric morphometrics with both front and side views of flowers.

A taxonomic, genetic and ecological data resource for the vascular plants of Britain and Ireland.

The vascular flora of Britain and Ireland is among the most extensively studied in the world, but the current knowledge base is fragmentary, with taxonomic, ecological and genetic information scattered across different resources. Here we present the first comprehensive data repository of native and alien species optimized for fast and easy online access for ecological, evolutionary and conservation analyses. The inventory is based on the most recent reference flora of Britain and Ireland, with taxon names linked to unique Kew taxon identifiers and DNA barcode data. Our data resource for 3,227 species and 26 traits includes existing and unpublished genome sizes, chromosome numbers and life strategy and life-form assessments, along with existing data on functional traits, species distribution metrics, hybrid propensity, associated biomes, realized niche description, native status and geographic origin of alien species. This resource will facilitate both fundamental and applied research and enhance our understanding of the flora's composition and temporal changes to inform conservation efforts in the face of ongoing climate change and biodiversity loss.

A phylogeny of Antirrhinum reveals parallel evolution of alpine morphology.

Parallel evolution of similar morphologies in closely related lineages provides insight into the repeatability and predictability of evolution. In the genus Antirrhinum (snapdragons), as in other plants, a suite of morphological characters are associated with adaptation to alpine environments. We tested for parallel trait evolution in Antirrhinum by investigating phylogenetic relationships using restriction-site associated DNA (RAD) sequencing. We then associated phenotypic information to our phylogeny to reconstruct the patterns of morphological evolution and related this to evidence for hybridisation between emergent lineages. Phylogenetic analyses showed that the alpine character syndrome is present in multiple groups, suggesting that Antirrhinum has repeatedly colonised alpine habitats. Dispersal to novel environments happened in the presence of intraspecific and interspecific gene flow. We found support for a model of parallel evolution in Antirrhinum. Hybridisation in natural populations, and a complex genetic architecture underlying the alpine morphology syndrome, support an important role of natural selection in maintaining species divergence in the face of gene flow.

A Draft Genome of the Ginger Species Alpinia nigra and New Insights into the Genetic Basis of Flexistyly.

Angiosperms possess various strategies to ensure reproductive success, such as stylar polymorphisms that encourage outcrossing. Here, we investigate the genetic basis of one such dimorphism that combines both temporal and spatial separation of sexual function, termed flexistyly. It is a floral strategy characterised by the presence of two morphs that differ in the timing of stylar movement. We performed a de novo assembly of the genome of Alpinia nigra using high-depth genomic sequencing. We then used Pool-seq to identify candidate regions for flexistyly based on allele frequency or coverage differences between pools of anaflexistylous and cataflexistylous morphs. The final genome assembly size was 2 Gb, and showed no evidence of recent polyploidy. The Pool-seq did not reveal large regions with high FST values, suggesting large structural chromosomal polymorphisms are unlikely to underlie differences between morphs. Similarly, no region had a 1:2 mapping depth ratio which would be indicative of hemizygosity. We propose that flexistyly is governed by a small genomic region that might be difficult to detect with Pool-seq, or a complex genomic region that proved difficult to assemble. Our genome will be a valuable resource for future studies of gingers, and provides the first steps towards characterising this complex floral phenotype.

Performance of generalist hemiparasitic Euphrasia across a phylogenetically diverse host spectrum.

Generalist hemiparasites may attach to many different host species and experience complex parasite-host interactions. How these parasite-host interactions impact on the fitness of hemiparasitic plants remain largely unknown. We used experimentally tractable eyebrights (Euphrasia, Orobanchaceae) to understand parasite-host species interactions affecting the performance of a generalist hemiparasitic plant. Common garden experiments were carried out measuring Euphrasia performance across 45 diverse hosts and in different parasite-host combinations. We show that variation in hemiparasite performance can be attributed mainly to host species and host phylogenetic relationships (λ = 0.82; 0.17-1.00 CI). When variation in performance is considered temporally, annual host species cause earlier flowering, and lead to poorer performance late in the season. While Euphrasia species typically perform similarly on a given host species, some eyebrights show more specialized parasite-host species interactions. Our results show that generalist hemiparasites only benefit from attaching to a limited, but phylogenetically divergent, subset of hosts. The conserved responses of divergent Euphrasia species suggest hemiparasite performance is affected by common host attributes. However, evidence for more complex parasite-host species interactions show that a generalist hemiparasite can potentially respond to individual host selection pressures and may adapt to local host communities.

The nature of intraspecific and interspecific genome size variation in taxonomically complex eyebrights.

BACKGROUND AND AIMS: Genome size varies considerably across the diversity of plant life. Although genome size is, by definition, affected by genetic presence/absence variants, which are ubiquitous in population sequencing studies, genome size is often treated as an intrinsic property of a species. Here, we studied intra- and interspecific genome size variation in taxonomically complex British eyebrights (Euphrasia, Orobanchaceae). Our aim is to document genome size diversity and investigate underlying evolutionary processes shaping variation between individuals, populations and species. METHODS: We generated genome size data for 192 individuals of diploid and tetraploid Euphrasia and analysed genome size variation in relation to ploidy, taxonomy, population affiliation and geography. We further compared the genomic repeat content of 30 samples. KEY RESULTS: We found considerable intraspecific genome size variation, and observed isolation-by-distance for genome size in outcrossing diploids. Tetraploid Euphrasia showed contrasting patterns, with genome size increasing with latitude in outcrossing Euphrasia arctica, but with little genome size variation in the highly selfing Euphrasia micrantha. Interspecific differences in genome size and the genomic proportions of repeat sequences were small. CONCLUSIONS: We show the utility of treating genome size as the outcome of polygenic variation. Like other types of genetic variation, such as single nucleotide polymorphisms, genome size variation may be affected by ongoing hybridization and the extent of population subdivision. In addition to selection on associated traits, genome size is predicted to be affected indirectly by selection due to pleiotropy of the underlying presence/absence variants.

Lineage-specific plastid degradation in subtribe Gentianinae (Gentianaceae).

The structure and sequence of plastid genomes is highly conserved across most land plants, except for a minority of lineages that show gene loss and genome degradation. Understanding the early stages of plastome degradation may provide crucial insights into the repeatability and predictability of genomic evolutionary trends. We investigated these trends in subtribe Gentianinae of the Gentianaceae, which encompasses ca. 450 species distributed around the world, particularly in alpine and subalpine environments. We sequenced, assembled, and annotated the plastomes of 41 species, representing all six genera in subtribe Gentianinae and all main sections of the species-rich genus Gentiana L. We reconstructed the phylogeny, estimated divergence times, investigated the phylogenetic distribution of putative gene losses, and related these to substitution rate shifts and species' habitats. We obtained a strongly supported topology consistent with earlier studies, with all six genera in Gentianinae recovered as monophyletic and all main sections of Gentiana having full support. While closely related species have very similar plastomes in terms of size and structure, independent gene losses, particularly of the ndh complex, have occurred in multiple clades across the phylogeny. Gene loss was usually associated with a shift in the boundaries of the small single-copy and inverted repeat regions. Substitution rates were variable between clades, with evidence for both elevated and decelerated rate shifts. Independent lineage-specific loss of ndh genes occurred at a wide range of times, from Eocene to Pliocene. Our study illustrates that diverse degradation patterns shape the evolution of the plastid in this species-rich plant group.

Population genomic and historical analysis suggests a global invasion by bridgehead processes in Mimulus guttatus.

Imperfect historical records and complex demographic histories present challenges for reconstructing the history of biological invasions. Here, we combine historical records, extensive worldwide and genome-wide sampling, and demographic analyses to investigate the global invasion of Mimulus guttatus from North America to Europe and the Southwest Pacific. By sampling 521 plants from 158 native and introduced populations genotyped at >44,000 loci, we determined that invasive M. guttatus was first likely introduced to the British Isles from the Aleutian Islands (Alaska), followed by admixture from multiple parts of the native range. We hypothesise that populations in the British Isles then served as a bridgehead for vanguard invasions worldwide. Our results emphasise the highly admixed nature of introduced M. guttatus and demonstrate the potential of introduced populations to serve as sources of secondary admixture, producing novel hybrids. Unravelling the history of biological invasions provides a starting point to understand how invasive populations adapt to novel environments.

Barcode UK: A complete DNA barcoding resource for the flowering plants and conifers of the United Kingdom.

DNA barcoding and metabarcoding provide new avenues for investigating biological systems. These techniques require well-curated reference libraries with extensive coverage. Generating an exhaustive national DNA barcode reference library can open up new avenues of research in ecology, evolution and conservation, yet few studies to date have created such a resource. In plant DNA barcoding, herbarium collections provide taxonomically robust material but also pose challenges in lab processing. Here, we present a national DNA barcoding resource covering all of the native flowering plants and conifers of the United Kingdom. This represents 1,482 plant species, with the majority of specimens (81%) sourced from herbaria. Using Sanger sequencing of the plant DNA barcode markers, rbcL, matK, and ITS2, at least one DNA barcode was retrieved from 98% of the UK flora. We sampled from multiple individuals, resulting in a species coverage for rbcL of 96% (4,477 sequences), 90% for matK (3,259 sequences) and 75% for ITS2 (2,585 sequences). Sequence recovery was lower for herbarium material compared to fresh collections, with the age of the specimen having a significant effect on the success of sequence recovery. Species level discrimination was highest with ITS2, however, the ability to successfully retrieve a sequence was lowest for this region. Analyses of the genetic distinctiveness of species across a complete flora showed DNA barcoding to be informative for all but the most taxonomically complex groups. The UK flora DNA barcode reference library provides an important resource for many applications that require plant identification from DNA.

Recurrent hybridization underlies the evolution of novelty in Gentiana (Gentianaceae) in the Qinghai-Tibetan Plateau.

The Qinghai-Tibetan Plateau (QTP) and adjacent areas are centres of diversity for several alpine groups. Although it is known that the QTP acted as a source area for diversification of the alpine genus Gentiana, the evolutionary processes underlying diversity in this genus, especially the formation of narrow endemics, are still poorly understood. Hybridization has been proposed as a driver of plant endemism in the QTP but few cases have been documented with genetic data. Here, we describe a new endemic species in Gentiana section Cruciata as G. hoae sp. nov., and explore its evolutionary history with complete plastid genomes and nuclear ribosomal internal transcribed spacer sequence data. Genetic divergence within G. hoae ~3 million years ago was followed by postglacial expansion on the QTP, suggesting Pleistocene glaciations as a key factor shaping the population history of G. hoae. Furthermore, a mismatch between plastid and nuclear data suggest that G. hoae participated in historical hybridization, while population sequencing show this species continues to hybridize with the co-occurring congener G. straminea in three locations. Our results indicate that hybridization may be a common process in the evolution of Gentiana and may be widespread among recently diverged taxa of the QTP.