Confirmation of natural hybrids between Gentiana straminea and G. siphonantha (Gentianaceae) based on molecular evidence

Insights from the first chromosome-level genome assembly of the alpine gentian Gentiana straminea Maxim

Gentiana straminea Maxim. is a perennial herb and mainly distributed in the Qinghai-Tibetan Plateau. To adapt to the extreme environment, it has developed particular morphological, physiological, and genetic structures. Also, rich in iridoids, it is one of the original plants of traditional Chinese herb 'Qinjiao'. Herein, we present its first chromosome-level genome sequence assembly and compare it with the genomes of other Gentiana species to facilitate the analysis of genomic characteristics. The assembled genome size of G. straminea was 1.25 Gb, with a contig N50 of 7.5 Mb. A total of 96.08% of the genome sequences was anchored on 13 pseudochromosomes, with a scaffold N50 of 92.70 Mb. A total of 54,310 protein-coding genes were predicted, 80.25% of which were functionally annotated. Comparative genomic analyses indicated that G. straminea experienced two whole-genome duplication events after the γ whole-genome triplication with other eudicots, and it diverged from other Gentiana species at ~3.2 Mya. A total of 142 enzyme-coding genes related to iridoid biosynthesis were identified in its genome. Additionally, we identified differences in the number and expression patterns of iridoid biosynthetic pathway genes in G. straminea compared with two other Gentiana species by integrating whole-genome sequence and transcriptomic analyses.

Comparison of plastid genomes and ITS of two sister species in Gentiana and a discussion on potential threats for the endangered species from hybridization

BACKGROUND

Gentiana rigescens Franchet is an endangered medicinal herb from the family Gentianaceae with medicinal values. Gentiana cephalantha Franchet is a sister species to G. rigescens possessing similar morphology and wider distribution. To explore the phylogeny of the two species and reveal potential hybridization, we adopted next-generation sequencing technology to acquire their complete chloroplast genomes from sympatric and allopatric distributions, as along with Sanger sequencing to produce the nrDNA ITS sequences.

RESULTS

The plastid genomes were highly similar between G. rigescens and G. cephalantha. The lengths of the genomes ranged from 146,795 to 147,001 bp in G. rigescens and from 146,856 to 147,016 bp in G. cephalantha. All genomes consisted of 116 genes, including 78 protein-coding genes, 30 tRNA genes, four rRNA genes and four pseudogenes. The total length of the ITS sequence was 626 bp, including six informative sites. Heterozygotes occurred intensively in individuals from sympatric distribution. Phylogenetic analysis was performed based on chloroplast genomes, coding sequences (CDS), hypervariable sequences (HVR), and nrDNA ITS. Analysis based on all the datasets showed that G. rigescens and G. cephalantha formed a monophyly. The two species were well separated in phylogenetic trees using ITS, except for potential hybrids, but were mixed based on plastid genomes. This study supports that G. rigescens and G. cephalantha are closely related, but independent species. However, hybridization was confirmed to occur frequently between G. rigescens and G. cephalantha in sympatric distribution owing to the lack of stable reproductive barriers. Asymmetric introgression, along with hybridization and backcrossing, may probably lead to genetic swamping and even extinction of G. rigescens.

CONCLUSION

G. rigescens and G. cephalantha are recently diverged species which might not have undergone stable post-zygotic isolation. Though plastid genome shows obvious advantage in exploring phylogenetic relationships of some complicated genera, the intrinsic phylogeny was not revealed because of matrilineal inheritance here; nuclear genomes or regions are hence crucial for uncovering the truth. As an endangered species, G. rigescens faces serious threats from both natural hybridization and human activities; therefore, a balance between conservation and utilization of the species is extremely critical in formulating conservation strategies.

Comparative analysis of mitochondrial genomes of two alpine medicinal plants of Gentiana (Gentianaceae)

Gentiana crassicaulis and G. straminea are alpine plants of Gentiana with important medicinal value and complex genetic backgrounds. In this study, the mitochondrial genomes (mtDNAs) of these two species were sequenced. The mtDNAs of G. crassicaulis and G. straminea are 368,808 and 410,086 bp long, respectively, 52 and 49 unique genes are annotated in the two species, and the gene arrangement varies widely. Compared to G. crassicaulis, G. straminea loses three effective genes, namely atp6, trnG-GCC and trnV-GAC. As a pseudogene, the atp6 gene of G. straminea is incomplete, which is rare in higher plants. We detected 1696 and 1858 pairs of long repeats and 213 SSRs and 250 SSs in the mtDNAs of G. crassicaulis and G. straminea, respectively. There are 392 SNPs and 18 InDels between the two genomes, and syntenic sequence and structural variation analysis show low collinearity between the two genomes. Chloroplast DNA transferring to mtDNA is observed in both species, and 46,511 and 55,043 bp transferred segments containing three tRNA genes are identified, respectively. Comparative analysis of mtDNAs of G. crassicaulis, G. straminea and four species of Gentianales determined 18 core genes, and there is no specific gene in G. crassicaulis and G. straminea. The phylogenetic tree based on mtDNAs places Gentianaceae in a branch of Gentianales. This study is the first to analyze the mtDNAs of Gentianaceae, which could provide information for analysis of the structure of mtDNAs of higher plants and phylogenetic research of Gentianaceae and Gentianales.

Population Transcriptomics Reveals Gene Flow and Introgression Between Two Non-sister Alpine Gentians

Distributional shifts driven by Quaternary climatic oscillations have been suggested to cause interspecific hybridization and introgression. In this study, we aimed to test this hypothesis by using population transcriptomes and coalescent modeling of two alpine none-sister gentians. Previous studies suggested that historical hybridizations occurred between Gentiana siphonantha and G. straminea in the high-altitude Qinghai-Tibet Plateau although both species are not sister to each other with the most recent divergence. In the present study, we sequenced transcriptomes of 33 individuals from multiple populations of G. siphonantha and G. straminea. The two species are well delimited by nuclear genomic SNPs while phylogenetic analyses of plastomes clustered one G. straminea individual into the G. siphonantha group. Further population structure analyses of the nuclear SNPs suggested that two populations of G. siphonantha were admixed with around 15% ancestry from G. straminea. These analyses suggested genetic introgressions from G. straminea to G. siphonantha. In addition, our coalescent-based modeling results revealed that gene flow occurred between the two species since Last Glacier Maximum after their initial divergence, which might have leaded to the observed introgressions. Our results underscore the significance of transcriptome population data in determining timescale of interspecific gene flow and direction of the resulting introgression.

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.

Population subdivision and hybridization in a species complex of Gentiana in the Qinghai-Tibetan Plateau

BACKGROUND AND AIMS

Hosting several global biodiversity hotspots, the region of the Qinghai-Tibetan Plateau (QTP) is exceptionally species-rich and harbours a remarkable level of endemism. Yet, despite a growing number of studies, factors fostering divergence, speciation and ultimately diversity remain poorly understood for QTP alpine plants. This is particularly the case for the role of hybridization. Here, we explored the evolutionary history of three closely related Gentiana endemic species, and tested whether our results supported the mountain geo-biodiversity hypothesis (MGH).

METHODS

We genotyped 69 populations across the QTP with one chloroplast marker and 12 nuclear microsatellite loci. We performed phylogeographical analysis, Bayesian clustering, approximate Bayesian computation and principal components analysis to explore their genetic relationship and evolutionary history. In addition, we modelled their distribution under different climates.

KEY RESULTS

Each species was composed of two geographically distinct clades, corresponding to the south-eastern and north-western parts of their distribution. Thus Gentiana veitchiorum and G. lawrencei var. farreri, which diverged recently, appear to have shared at least refugia in the past, from which their range expanded later on. Indeed, climatic niche modelling showed that both species went through continuous expansion from the Last Interglacial Maximum to the present day. Moreover, we have evidence of hybridization in the northwest clade of G. lawrencei var. farreri, which probably occurred in the refugium located on the plateau platform. Furthermore, phylogenetic and population genetic analyses suggested that G. dolichocalyx should be a geographically limited distinct species with low genetic differentiation from G. lawrencei var. farreri.

CONCLUSIONS

Climatic fluctuations in the region of the QTP have played an important role in shaping the current genetic structure of G. lawrencei var. farreri and G. veitchiorum. We argue that a species pump effect did occur prior to the Last Interglacial Maximum, thus lending support to the MGH. However, our results do depart from expectations as suggested in the MGH for more recent distribution range and hybridization dynamics.

Comparative Chloroplast Genome Analyses of Species in Gentiana section Cruciata (Gentianaceae) and the Development of Authentication Markers

Gentiana section Cruciata is widely distributed across Eurasia at high altitudes, and some species in this section are used as traditional Chinese medicine. Accurate identification of these species is important for their utilization and conservation. Due to similar morphological and chemical characteristics, correct discrimination of these species still remains problematic. Here, we sequenced three complete chloroplast (cp) genomes (G. dahurica, G. siphonantha and G. officinalis). We further compared them with the previously published plastomes from sect. Cruciata and developed highly polymorphic molecular markers for species authentication. The eight cp genomes shared the highly conserved structure and contained 112 unique genes arranged in the same order, including 78 protein-coding genes, 30 tRNAs, and 4 rRNAs. We analyzed the repeats and nucleotide substitutions in these plastomes and detected several highly variable regions. We found that four genes (accD, clpP, matK and ycf1) were subject to positive selection, and sixteen InDel-variable loci with high discriminatory powers were selected as candidate barcodes. Our phylogenetic analyses based on plastomes further confirmed the monophyly of sect. Cruciata and primarily elucidated the phylogeny of Gentianales. This study indicated that cp genomes can provide more integrated information for better elucidating the phylogenetic pattern and improving discriminatory power during species authentication.

The Complete Plastome Sequences of Seven Species in Gentiana sect. Kudoa (Gentianaceae): Insights Into Plastid Gene Loss and Molecular Evolution

The chloroplast (cp) genome is useful in the study of phylogenomics, molecular dating, and molecular evolution. Gentiana sect. Kudoa is a predominantly alpine flowering plant that is valued for its contributions to medicine, ecology, and horticulture. Previous evolutionary studies showed that the plastid gene loss pattern and intra-sectional phylogenetics in sect. Kudoa are still unclear. In this study, we compared 11 Gentiana plastomes, including 7 newly sequenced plastomes from sect. Kudoa, to represent its three serious: ser. Ornatae, ser. Verticillatae, and ser. Monanthae. The cp genome sizes of the seven species ranged from 137,278 to 147,156 bp. The plastome size variation mainly occurred in the small single-copy and long single-copy regions rather than the inverted repeat regions. Compared with sect. Cruciata, the plastomes in ser. Ornatae and ser. Verticillatae had lost approximately 11 kb of sequences containing 11 ndh genes. Conversely, far fewer losses were observed in ser. Monanthae. The phylogenetic tree revealed that sect. Kudoa was not monophyletic and that ser. Monanthae was more closely related to other sections rather than sect. Kudoa. The molecular dating analysis indicated that ser. Monanthae and sect. Kudoa diverged around 8.23 Ma. In ser. Ornatae and ser. Verticillatae, the divergence occurred at around 0.07-1.78 Ma. The nucleotide diversity analysis indicated that the intergenic regions trnH-psbA, trnK-trnQ, ycf3-trnS and rpl32-trnL constituted divergence hotspots in both sect. Kudoa and Gentiana, and would be useful for future phylogenetic and population genetic studies.

Genetic homogenization of the nuclear ITS loci across two morphologically distinct gentians in their overlapping distributions in the Qinghai-Tibet Plateau

Interspecific hybridization and introgression can lead to partial genetic homogenization at certain neutral loci between morphologically distinct species and may obscure the species delimitation using nuclear genes. In this study, we investigated this phenomenon through population genetic survey of two alpine plants (Gentiana siphonantha and G. straminea) in the Qinghai-Tibet Plateau, where the distributions of two species are partly overlapped. We identified two clusters of chloroplast DNA haplotypes which correspond to the two species, and three clusters of ITS ribotypes. In addition to clusters specific to each species, the third ITS cluster, which was most likely derived from hybridization between the other two clusters and subsequent recombination and concerted evolution, was widely shared by two species in their adjacent areas. In contrast to the morphological distinctiveness of the two species, interspecific gene flow possibly led to genetic homogenization at their ITS loci. The new ITS lineage recovered for species in adjacent areas is distinctly different from original lineages found in allopatric areas. These findings may have general implications for our understanding of cryptic changes at some genetic loci caused by interspecific gene flow in the history, and they indicate that species delimitation should be based on a combination of both nuclear and chloroplast DNA sequence variations.