Comparative chloroplast-specific SNP and nSCoT markers analysis and population structure study in kiwifruit plants.

BACKGROUND: Kiwifruit (Actinidiaceae family) is an economically important fruit tree in China and New Zealand. It is a typical dioecious plant that has undergone frequent natural hybridization, along with chromosomal ploidy diversity within the genus Actinidia, resulting in higher genetic differences and horticultural diversity between interspecific and intraspecific traits. This diversity provides a rich genetic base for breeding. China is not only the original center of speciation for the Actinidia genus but also its distribution center, housing the most domesticated species: A. chinensis var. chinensis, A. chinensis var. deliciosa, A. arguta, and A. polygama. However, there have been relatively few studies on the application of DNA markers and the genetic basis of kiwifruit plants. By combining information from chloroplast-specific SNPs and nuclear SCoT (nSCoT) markers, we can uncover complementary aspects of genetic variation, population structure, and evolutionary relationships. In this study, one chloroplast DNA (cpDNA) marker pair was selected out of nine cpDNA candidate pairs. Twenty nSCoT markers were selected and used to assess the population structure and chloroplast-specific DNA haplotype diversity in 55 kiwifruit plants (Actinidia), including 20 samples of A. chinensis var. chinensis, 22 samples of A. chinensis var. deliciosa, 11 samples of A. arguta, and two samples of A. polygama, based on morphological observations collected from China. RESULTS: The average genetic distance among the 55 samples was 0.26 with chloroplast-specific SNP markers and 0.57 with nSCoT markers. The Mantel test revealed a very small correlation (r = 0.21). The 55 samples were categorized into different sub-populations using Bayesian analysis, the Unweighted Pair Group Method with the Arithmetic Mean (UPGMA), and the Principal Component Analysis (PCA) method, respectively. Based on the analysis of 205 variable sites, a total of 15 chloroplast-specific DNA haplotypes were observed, contributing to a higher level of polymorphism with an Hd of 0.78. Most of the chloroplast-specific DNA haplotype diversity was distributed among populations, but significant diversity was also observed within populations. H1 was shared by 24 samples, including 12 of A. chinensis var. chinensis and 12 of A. chinensis var. deliciosa, indicating that H1 is an ancient and dominant haplotype among the 55 chloroplast-specific sequences. H2 may not have evolved further.The remaining haplotypes were rare and unique, with some appearing to be exclusive to a particular variety and often detected in single individuals. For example, the H15 haplotype was found exclusively in A. polygama. CONCLUSION: The population genetic variation explained by chloroplast-specific SNP markers has greater power than that explained by nSCoTs, with chloroplast-specific DNA haplotypes being the most efficient. Gene flow appears to be more evident between A. chinensis var. chinensis and A. chinensis var. deliciosa, as they share chloroplast-specific DNA haplotypes, In contrast, A.arguta and A. polygama possess their own characteristic haplotypes, derived from the haplotype of A. chinensis var. chinensis. Compared with A. chinensis, the A.arguta and A. polygama showed better grouping. It also seems crucial to screen out, for each type of molecular marker, especially haplotypes, the core markers of the Actinidia genus.

Conserved DNA sequence analysis reveals the phylogeography and evolutionary events of Akebia trifoliata in the region across the eastern edge of the Tibetan Plateau and subtropical China.

BACKGROUND: The eastern edge of the Qinghai‒Tibet Plateau (QTP) and subtropical China have various regions where plant species originate and thrive, but these regions have been the focus of very few integrative studies. Here, we elucidated the phylogeographic structure of a continuous and widespread Akebia trifoliata population across these two regions. RESULTS: Sixty-one populations consisting of 391 genotypes were examined to assess population diversity and structure via network distribution analysis, maximum likelihood phylogenetic tree reconstruction, divergence time estimation, demographic history inference, and ancestral area reconstruction of both conserved internal transcribed spacer (ITS) and chloroplast (rps16) DNA sequences. The results showed that the ITS region was more variable than the rps16 region and could be suitable for studying intraspecific phylogeography. The A. trifoliata population displayed high genetic diversity, genetic differentiation and obvious phylogeographical structure, possibly originating on the eastern QTP, expanding during the last glacial-interglacial cycle, diverging in the early Pleistocene and middle Pleistocene, and extensively migrating thereafter. The migration route from west to east along rivers could be largely responsible for the long-distance dispersal of this species, while three main refuges (Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) with multiple ice shelters facilitated its wide distribution. CONCLUSIONS: Our results suggested that the from west to east long migration accompanying with the minor short reciprocal migration in the south-north direction, and the three main refuges (the Qinba Mountains, Nanling Mountains and Yunnan-Guizhou Plateau) contributed to the extant geographical distribution of A. trifoliata. In addition, this finding also strongly reduced the discrepancy between glacial contraction and postglacial expansion and the in situ survival hypothesis by simultaneously considering the existence of many similar climate-related ecological niches and migration influences.

Wasabi Gone Wild? Origin and Characterization of the Complete Plastomes of Ulleung Island Wasabi (Eutrema japonicum; Brassicaceae) and Other Cultivars in Korea.

Korean wasabi occurs naturally on the young oceanic, volcanic Ulleung Island off the east coast of the Korean Peninsula. Although the Ulleung Island wasabi is reported as Eutrema japonicum and has been suggested to be morphologically identical to cultivars in Korea, very little is known about its taxonomic identity and relationship with other cultivars. In this study, we sequenced the complete chloroplast DNA sequences of three naturally occurring Ulleung Island wasabi plants and six cultivars ('Daewang', 'Daruma', 'Micado', 'Orochi', 'Green Thumb', and 'Shogun') from continental Korea and determined the taxonomic identity of Korean wasabi on Ulleung Island. The size and organization of the complete chloroplast genomes of the nine accessions were nearly identical to those of previously reported wasabi cultivars. In addition, phylogenetic analysis based on the complete plastomes suggested that Ulleung Island wasabi most likely comprises various wasabi cultivars with three chlorotypes ('Shogun', 'Green Thumb', and a unique Chusan type). Based on the complete plastomes, we identified eight chlorotypes for the major wasabi cultivars and the Ulleung Island wasabi. Two major groups (1-'Mazuma' and 'Daruma', and 2-'Fujidaruma'/'Shimane No. 3'/Ulleung Island wasabi/five cultivars in Korea) were also identified based on mother line genealogical history. Furthermore, different types of variations (mutations, insertions/deletions (indels), mononucleotide repeats, and inversions) in plastomes were identified to distinguish different cultivar lines and five highly divergent hotspots. The nine newly obtained complete plastomes are valuable organelle genomic resources for species identification and infraspecific phylogeographic studies on wild and cultivated wasabi.

Extraction and analysis of high-quality chloroplast DNA with reduced nuclear DNA for medicinal plants.

BACKGROUND: Obtaining high-quality chloroplast genome sequences requires chloroplast DNA (cpDNA) samples that meet the sequencing requirements. The quality of extracted cpDNA directly impacts the efficiency and accuracy of sequencing analysis. Currently, there are no reported methods for extracting cpDNA from Erigeron breviscapus. Therefore, we developed a suitable method for extracting cpDNA from E. breviscapus and further verified its applicability to other medicinal plants. RESULTS: We conducted a comparative analysis of chloroplast isolation and cpDNA extraction using modified high-salt low-pH method, the high-salt method, and the NaOH low-salt method, respectively. Subsequently, the number of cpDNA copies relative to the nuclear DNA (nDNA ) was quantified via qPCR. As anticipated, chloroplasts isolated from E. breviscapus using the modified high-salt low-pH method exhibited intact structures with minimal cell debris. Moreover, the concentration, purity, and quality of E. breviscapus cpDNA extracted through this method surpassed those obtained from the other two methods. Furthermore, qPCR analysis confirmed that the modified high-salt low-pH method effectively minimized nDNA contamination in the extracted cpDNA. We then applied the developed modified high-salt low-pH method to other medicinal plant species, including Mentha haplocalyx, Taraxacum mongolicum, and Portulaca oleracea. The resultant effect on chloroplast isolation and cpDNA extraction further validated the generalizability and efficacy of this method across different plant species. CONCLUSIONS: The modified high-salt low-pH method represents a reliable approach for obtaining high-quality cpDNA from E. breviscapus. Its universal applicability establishes a solid foundation for chloroplast genome sequencing and analysis of this species. Moreover, it serves as a benchmark for developing similar methods to extract chloroplast genomes from other medicinal plants.

The size diversity of the Pteridaceae family chloroplast genome is caused by overlong intergenic spacers.

BACKGROUND: While the size of chloroplast genomes (cpDNAs) is often influenced by the expansion and contraction of inverted repeat regions and the enrichment of repeats, it is the intergenic spacers (IGSs) that appear to play a pivotal role in determining the size of Pteridaceae cpDNAs. This provides an opportunity to delve into the evolution of chloroplast genomic structures of the Pteridaceae family. This study added five Pteridaceae species, comparing them with 36 published counterparts. RESULTS: Poor alignment in the non-coding regions of the Pteridaceae family was observed, and this was attributed to the widespread presence of overlong IGSs in Pteridaceae cpDNAs. These overlong IGSs were identified as a major factor influencing variations in cpDNA size. In comparison to non-expanded IGSs, overlong IGSs exhibited significantly higher GC content and were rich in repetitive sequences. Species divergence time estimations suggest that these overlong IGSs may have already existed during the early radiation of the Pteridaceae family. CONCLUSIONS: This study reveals new insights into the genetic variation, evolutionary history, and dynamic changes in the cpDNA structure of the Pteridaceae family, providing a fundamental resource for further exploring its evolutionary research.

Genetic diversities in wild and cultivated populations of the two closely-related medical plants species, Tripterygium Wilfordii and T. Hypoglaucum (Celastraceae).

BACKGROUND: The sustainable supply of medicinal plants is important, and cultivating and domesticating them has been suggested as an optimal strategy. However, this can lead to a loss of genetic diversity. Tripterygium wilfordii Hook. f. is a medicinal plant commonly used in traditional Chinese medicine, but its wild populations are dwindling due to excessive harvesting. To protect the species and meet the increasing demand, it is urgent to cultivate it on a large scale. However, distinguishing between T. wilfordii and T. hypoglaucum, two similar species with different medicinal properties, is challenging. Therefore, it is crucial to understand the genetic diversity and population structure of these species for their sustainable utilization. RESULTS: In this study, we investigated the genetic diversity and population structure of the two traditional medicinal semiwoody vines plant species, Tripterygium wilfordii and T. hypoglaucum, including wild and cultivated populations using chloroplast DNA (cpDNA) sequences and microsatellite loci. Our results indicated that the two species maintain a high level of genetic divergence, indicating possible genetic bases for the different contents of bioactive compounds of the two species. T. wilfordii showed lower genetic diversity and less subdivided population structures of both markers than T. hypoglaucum. The potential factors in shaping these interesting differences might be differentiated pollen-to-seed migration rates, interbreeding, and history of population divergence. Analyses of cpDNA and microsatellite loci supported that the two species are genetically distinct entities. In addition, a significant reduction of genetic diversity was observed for cultivated populations of the two species, which mainly resulted from the small initial population size and propagated vegetative practice during their cultivation. CONCLUSION: Our findings indicate significant genetic divergence between T. wilfordii and T. hypoglaucum. The genetic diversity and population structure analyses provide important insights into the sustainable cultivation and utilization of these medicinal plants. Accurate identification and conservation efforts are necessary for both species to ensure the safety and effectiveness of crude drug use. Our study also highlighted the importance of combined analyses of different DNA markers in addressing population genetics of medicinal plants because of the contrasts of inheritance and rates of gene flow. Large-scale cultivation programs should consider preserving genetic diversity to enhance the long-term sustainability of T. wilfordii and T. hypoglaucum. Our study proposed that some populations showed higher genetic diversity and distinctness, which can be considered with priority for conservation and as the sources for future breeding and genetic improvement.

Evolutionary history of an Irano-Turanian cushion-forming legume (Onobrychis cornuta).

The Irano-Turanian region is one of the largest floristic regions in the world and harbors a high percentage of endemics, including cushion-like and dwarf-shrubby taxa. Onobrychis cornuta is an important cushion-forming element of the subalpine/alpine flora of the Irano-Turanian floristic region. To specify the genetic diversity among the populations of this species (including individuals of O. elymaitica), we employed nrDNA ITS and two noncoding regions of plastid DNA (rpl32-trnL(UAG) and trnT(UGU)-trnL(UAA)). The most striking feature of O. cornuta assemblages was the unexpectedly high nucleotide diversity in both the nDNA and cpDNA dataset. In the analyses of nuclear and plastid regions, 25 ribotypes and 42 haplotypes were found among 77 and 59 accessions, respectively, from Iran, Turkey, and Afghanistan. Network analysis of the datasets demonstrated geographic differentiation within the species. Phylogenetic analyses of all dataset retrieved O. cornuta as a non-monophyletic species due to the inclusion of O. elymaitica, comprising four distinct lineages. In addition, our analyses showed cytonuclear discordance between both nuclear and plastid topologies regarding the position of some O. cornuta individuals. The underlying causes of this inconsistency remain unclear. However, we speculate that chloroplast capture, incomplete lineage sorting, and introgression were the main reasons for this event. Furthermore, molecular dating analysis indicated that O. cornuta originated in the early Pliocene (around 4.8 Mya) and started to diversify throughout the Pliocene and in particular the Pleistocene. Moreover, O. elymaitica was reduced to a subspecific rank within the species.

Population Structure and Genetic Diversity of the Toona ciliata (Meliaceae) Complex Assayed with Chloroplast DNA Markers.

Toona ciliata is a deciduous or semi-deciduous tree species and belongs to the Toona genus of the Meliaceae family. Owing to low natural regeneration and over-exploitation, the species is listed as an endangered species at level II in China and its conservation has received increasing concern. Here, we sampled 447 individuals from 29 populations across the range-wide distribution of the T. ciliata complex in China and assessed their genetic variation using two chloroplast DNA markers. The results showed that the overall haplotype diversity and nucleotide diversity per site were high at h = 0.9767 and π = 0.0303 for the psbA-trnH fragment and h= 0.8999 and π = 0.0189 for the trnL-trnL fragment. Phylogenetic analysis supported the division of the natural distribution of T. ciliata complex into western and eastern regions. The genetic diversity was higher in the western region than in the eastern region, showing significant phylogeographic structure. Genetic differentiation among populations was moderate (Φst=42.87%), and the effects of isolation by distance (IBD) were significant. A neutrality test and mismatch distribution analysis indicated that the distribution of the T. ciliata complex generally did not expand, although a few local populations could likely expand after bottleneck effects. The overall results were complementary to and consolidated previous studies using mitochondrial and nuclear DNA markers. We finally discussed strategies for the genetic conservation of the T. ciliata complex.

Phylogeographical and population genetics of Polyspora sweet in China provides insights into its phylogenetic evolution and subtropical dispersal.

BACKGROUND: Geological movements and climatic fluctuations stand as pivotal catalysts driving speciation and phylogenetic evolution. The genus Polyspora Sweet (Theaceae), prominently found across the Malay Archipelagos and Indochina Peninsula in tropical Asia, exhibits its northernmost distribution in China. In this study, we investigated the evolutionary and biogeographical history of the genus Polyspora in China, shedding light on the mechanisms by which these species respond to ancient geological and climatic fluctuations. METHODS: Phylogenetic relationships of 32 representative species of Theaceae were reconstructed based on the chloroplast genome and ribosome 18-26 S rRNA datasets. Species divergence time was estimated using molecular clock and five fossil calibration. The phylogeography and population genetics in 379 individuals from 32 populations of eight species were analyzed using chloroplast gene sequences (trnH-psbA, rpoB-trnC and petN-psbM), revealing the glacial refugia of each species, and exploring the causes of the phylogeographic patterns. RESULTS: We found that Chinese Polyspora species diverged in the middle Miocene, showing a tropical-subtropical divergence order. A total of 52 haplotypes were identified by the combined chloroplast sequences. Chinese Polyspora exhibited a distinct phylogeographical structure, which could be divided into two clades and eight genealogical subdivisions. The divergence between the two clades occurred approximately 20.67 Ma. Analysis of molecular variance revealed that the genetic variation mainly occurred between species (77.91%). At the species level, Polyspora axillaris consists of three lineages, while P. speciosa had two lineages. The major lineages of Chinese Polyspora diverged between 12 and 15 Ma during the middle to late Miocene. The peak period of haplotype differentiation in each species occurred around the transition from the last interglacial to the last glacial period, approximately 6 Ma ago. CONCLUSION: The primary geographical distribution pattern of Chinese Polyspora was established prior to the last glacial maximum, and the population historical dynamics were relatively stable. The geological and climatic turbulence during the Quaternary glacial period had minimal impact on the distribution pattern of the genus. The genus coped with Quaternary climate turbulence by glacial in situ survival in multiple refuges. The Sino-Vietnam border and Nanling corridor might be the genetic mixing center of Polyspora.

Two lineages of Lemna aequinoctialis (Araceae, Lemnoideae) based on physiology, morphology, and phylogeny.

Lemna aequinoctialis Welw. is a widely spread species that has diverse physiological and molecular properties. Flower characteristics are important factors in deducing taxonomical status; however, owing to the rarity of flowering observations in Lemna, studying them has been a prolonged challenge. In this study, physiological and morphological analyses were conducted by inducing flowering, and molecular analysis was done based on the two chloroplast DNA loci (matK, atpF-atpH intergeneric spacer) of L. aequinoctialis sensu Landolt (1986) from 70 strains found in 70 localities in Japan, Korea, Thailand, and the US. In total, 752 flowering fronds from 13 strains were observed based on axenic conditions. Two different trends in flower organ development-protogyny and adichogamy-were detected in these strains. Their physiological traits were divided into two groups, showing different morphological features based on frond thickness, root cap, and anther sizes. Molecular analysis showed two lineages corresponding to two physiological groups. These were identified as L. aequinoctialis sensu Beppu et al. (1985) and L. aoukikusa Beppu et Murata based on the description of the nomenclature of L. aoukikusa. These were concluded as independent taxa and can be treated as different species. Furthermore, the distribution of L. aoukikusa is not only limited to Japan.

From forest to savanna and back to forest: Evolutionary history of the genus Dimorphandra (Fabaceae).

The tree genus Dimorphandra (Fabaceae), which contains 26 species divided into three subgenera, was studied using DNA sequence data from six chloroplast genome regions (cpDNA) and the nuclear internal transcribed spacer (ITS). The analyses, which included Bayesian phylogenies and haplotype networks, ancestral area reconstructions, and ecological niche modeling, allowed for exploring the evolutionary history of Dimorphandra. Within the subgenus Phaneropsia, the cpDNA sequence data were more closely-related to species from the genus Mora, while the ITS sequence data displayed a closer phylogenetic relationship with the subgenus Pocillum. This incongruence may be due to incomplete lineage sorting associated with ancient polymorphisms. The Amazonian Dimophandra lineages were highly polymorphic and divergent, while those from the Cerrado and the Atlantic Forest had low levels of polymorphisms. The Amazon likely gave rise to the Dimophandra lineage that produced the Cerrado species, while a Cerrado lineage likely gave rise to the Atlantic Forest species. Habitat shifts were identified as a key factor in shaping the late evolutionary history of Dimorphandra.

Complete mitochondrial and chloroplast DNA sequences of the freshwater green microalga Medakamo hakoo.

We report the complete organellar genome sequences of an ultrasmall green alga, Medakamo hakoo strain M-hakoo 311, which has the smallest known nuclear genome in freshwater green algae. Medakamo hakoo has 90.8-kb chloroplast and 36.5-kb mitochondrial genomes containing 80 and 33 putative protein-coding genes, respectively. The mitochondrial genome is the smallest in the Trebouxiophyceae algae studied so far. The GC content of the nuclear genome is 73%, but those of chloroplast and mitochondrial genomes are 41% and 35%, respectively. Codon usages in the organellar genomes have a different tendency from that in the nuclear genome. The organellar genomes have unique characteristics, such as the biased encoding of mitochondrial genes on a single strand and the absence of operon structures in chloroplast ribosomal genes. Medakamo hakoo will be helpful for understanding the evolution of the organellar genome and the regulation of gene expression in chloroplasts and mitochondria.

Demographic dynamics and molecular evolution of the rare and endangered subsect. Gerardianae of Pinus: insights from chloroplast genomes and mitochondrial DNA markers.

MAIN CONCLUSION: The divergence of subsect. Gerardianae was likely triggered by the uplift of the Qinghai-Tibetan Plateau and adjacent mountains. Pinus bungeana might have probably experienced expansion since Last Interglacial period. Historical geological and climatic oscillations have profoundly affected patterns of nucleotide variability, evolutionary history, and species divergence in numerous plants of the Northern Hemisphere. However, how long-lived conifers responded to geological and climatic fluctuations in East Asia remain poorly understood. Here, based on paternally inherited chloroplast genomes and maternally inherited mitochondrial DNA markers, we investigated the population demographic history and molecular evolution of subsect. Gerardianae (only including three species, Pinus bungeana, P. gerardiana, and P. squamata) of Pinus. A low level of nucleotide diversity was found in P. bungeana (π was 0.00016 in chloroplast DNA sequences, and 0.00304 in mitochondrial DNAs). The haplotype-based phylogenetic topology and unimodal distributions of demographic analysis suggested that P. bungeana probably originated in the southern Qinling Mountains and experienced rapid population expansion since Last Interglacial period. Phylogenetic analysis revealed that P. gerardiana and P. squamata had closer genetic relationship. The species divergence of subsect. Gerardianae occurred about 27.18 million years ago (Mya) during the middle to late Oligocene, which was significantly associated with the uplift of the Qinghai-Tibetan Plateau and adjacent mountains from the Eocene to the mid-Pliocene. The molecular evolutionary analysis showed that two chloroplast genes (psaI and ycf1) were under positive selection, the genetic lineages of P. bungeana exhibited higher transition and nonsynonymous mutations, which were involved with the strongly environmental adaptation. These findings shed light on the population evolutionary history of white pine species and provide striking insights for comprehension of their species divergence and molecular evolution.

Peripatric speciation within Torreya fargesii (Taxaceae) in the Hengduan Mountains inferred from multi-loci phylogeography.

BACKGROUND: The Hengduan Mountains (HDM) are one of the major global biodiversity hotspots in the world. Several evolutionary scenarios, especially in-situ diversification, have been proposed to account for the high species richness of temperate plants. However, peripatric speciation, an important mode of allopatric speciation, has seldom been reported in this region. RESULTS: Here, two chloroplast DNA regions and 14 nuclear loci were sequenced for 112 individuals from 10 populations of Torreya fargesii var. fargesii and 63 individuals from 6 populations of T. fargesii var. yunnanensis. Population genetic analyses revealed that the two varieties are well differentiated genetically (FST, 0.5765) and have uneven genetic diversity (π, 0.00221 vs. 0.00073 on an average of nuclear loci). The gene genealogical relationship showed that T. fargesii var. yunnanensis is inferred as derived from T. fargesii var. fargesii, which was further supported by the coalescent simulations (DIYABC, fastsimcoal2 and IMa2). By the coalescent simulations, the divergence time (~ 2.50-3.65 Ma) and the weak gene flow between the two varieties were detected. The gene flow was asymmetrical and only occurred in later stages of divergence, which is caused by second contact due to the population expansion (~ 0.61 Ma) in T. fargesii var. fargesii. In addition, niche modeling indicated that the two varieties are differentiated geographically and ecologically and have unbalanced distribution range. CONCLUSIONS: Overall, T. fargesii var. fargesii is always parapatric with respect to T. fargesii var. yunnanensis, and the latter derived from the former in peripatry of the HDM following a colonization from central China during the late Pliocene. Our findings demonstrate that peripatric speciation following dispersal events may be an important evolutionary scenario for the formation of biodiversity hotspot of the HDM.

Insights into the phylogenetic relationships and species boundaries of the Myricaria squamosa complex (Tamaricaceae) based on the complete chloroplast genome.

Myricaria plants are widely distributed in Eurasia and are helpful for windbreak and embankment protection. Current molecular evidence has led to controversy regarding species boundaries within the Myricaria genus and interspecific phylogenetic relationships between three specific species-M. bracteata, M. paniculata and M. squamosa-which have remained unresolved. This study treated these three unresolved taxa as a species complex, named the M. squamosa complex. The genome skimming approach was used to determine 35 complete plastome sequences and nuclear ribosomal DNA sequences for the said complex and other closely related species, followed by de novo assembly. Comparative analyses were conducted across Myricaria to identify the genome size, gene content, repeat type and number, SSR (simple sequence repeat) abundance, and codon usage bias of chloroplast genomes. Tree-based species delimitation results indicated that M. bracteata, M. paniculata and M. squamosa could not be distinguished and formed two monophyletic lineages (P1 and P2) that were clustered together. Compared to plastome-based species delimitation, the standard nuclear DNA barcode had the lowest species resolution, and the standard chloroplast DNA barcode and group-specific barcodes delimitated a maximum of four out of the five species. Plastid phylogenomics analyses indicated that the monophyletic M. squamosa complex is comprised of two evolutionarily significant units: one in the western Tarim Basin and the other in the eastern Qinghai-Tibet Plateau. This finding contradicts previous species discrimination and promotes the urgent need for taxonomic revision of the threatened genus Myricaria. Dense sampling and plastid genomes will be essential in this effort. The super-barcodes and specific barcode candidates outlined in this study will aid in further studies of evolutionary history.

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.

Molecular phylogeography and historical demography of a widespread herbaceous species from eastern North America, Podophyllum peltatum.

PREMISE: Glacial/interglacial cycles and topographic complexity are both considered to have shaped today's diverse phylogeographic patterns of taxa from unglaciated eastern North America (ENA). However, few studies have focused on the phylogeography and population dynamics of wide-ranging ENA herbaceous species occurring in forest understory habitat. We examined the phylogeographic pattern and evolutionary history of Podophyllum peltatum L., a widely distributed herb inhabiting deciduous forests of ENA. METHODS: Using chloroplast DNA (cpDNA) sequences and nuclear microsatellite loci, we investigated the population structure and genetic diversity of the species. Molecular dating, demographic history analyses, and ecological niche modeling were also performed to illustrate the phylogeographic patterns. RESULTS: Our cpDNA results identified three main groups that are largely congruent with boundaries along the Appalachian Mountains and the Mississippi River, two major geographic barriers in ENA. Populations located to the east of the Appalachians and along the central Appalachians exhibited relatively higher levels of genetic diversity. Extant lineages may have diverged during the late Miocene, and range expansions of different groups may have happened during the Pleistocene glacial/interglacial cycles. CONCLUSIONS: Our findings indicate that geographic barriers may have started to facilitate the population divergence in P. peltatum before the Pleistocene. Persistence in multiple refugia, including areas around the central Appalachians during the Quaternary glacial period, and subsequent expansions under hospitable climatic condition, especially westward expansion, are likely responsible for the species' contemporary genetic structure and phylogeographic pattern.

Nuclear and chloroplast DNA phylogeography reveals high genetic diversity and postglacial range expansion in Quercus mexicana.

PREMISE: Phylogeographical studies are fundamental for understanding factors that influence the spatial distribution of genetic lineages within species. Population expansions and contractions, distribution shifts, and climate changes are among the most important factors shaping the genetic compositions of populations. METHODS: We investigated the phylogeography of an endemic oak, Quercus mexicana (Fagaceae), which has a restricted distribution in northeastern Mexico along the Sierra Madre Oriental and adjacent areas. Nuclear and chloroplast DNA microsatellite markers were used to describe the genetic diversity and structure of 39 populations of Q. mexicana along its entire distribution area. We tested whether population expansion or contraction events influenced the genetic diversity and structure of the species. We also modeled the historical distributional range of Q. mexicana (for the Mid Holocene, the Last Glacial Maximum, and the Last Interglacial) to estimate the extent to which climate fluctuations have impacted the distribution of this oak species. RESULTS: Our results revealed high genetic diversity and low genetic structure in Q. mexicana populations. Ecological niche models suggested historical fluctuations in the distributional range of Q. mexicana. Historical range changes, gene flow, and physical barriers seem to have played an important role in shaping the phylogeographic structure of Q. mexicana. CONCLUSIONS: Our study indicates that the genetic structure of Q. mexicana may have been the result of responses of oak trees not only to heterogeneous environments present in the Sierra Madre Oriental and adjacent areas, but also to elevational and latitudinal shifts in response to climate changes in the past.

Adaptive evolution and co-evolution of chloroplast genomes in Pteridaceae species occupying different habitats: overlapping residues are always highly mutated.

BACKGROUND: The evolution of protein residues depends on the mutation rates of their encoding nucleotides, but it may also be affected by co-evolution with other residues. Chloroplasts function as environmental sensors, transforming fluctuating environmental signals into different physiological responses. We reasoned that habitat diversity may affect their rate and mode of evolution, which might be evidenced in the chloroplast genome. The Pteridaceae family of ferns occupy an unusually broad range of ecological niches, which provides an ideal system for analysis. RESULTS: We conducted adaptive evolution and intra-molecular co-evolution analyses of Pteridaceae chloroplast DNAs (cpDNAs). The results indicate that the residues undergoing adaptive evolution and co-evolution were mostly independent, with only a few residues being simultaneously involved in both processes, and these overlapping residues tend to exhibit high mutations. Additionally, our data showed that Pteridaceae chloroplast genes are under purifying selection. Regardless of whether we grouped species by lineage (which corresponded with ecological niches), we determined that positively selected residues mainly target photosynthetic genes. CONCLUSIONS: Our work provides evidence for the adaptive evolution of Pteridaceae cpDNAs, especially photosynthetic genes, to different habitats and sheds light on the adaptive evolution and co-evolution of proteins.

Population Genetic Analysis of Paris polyphylla var. yunnanensis Based on cpDNA Fragments.

Paris polyphylla var. yunnanensis is a well-known medicinal plant that is mainly distributed in Southwest China; however, its genetic diversity and biodiversity processes are poorly understood. In this study, the sequences of cpDNA trnL-trnF fragments of 15 wild populations and 17 cultivated populations of P. polyphylla var. yunnanensis were amplified, sequenced, and aligned to study the population genetics of this species. Genetic diversity was analyzed based on nucleotide diversity, haplotype diversity, Watterson diversity, population-level diversity, and species-level genetic diversity. Genetic structure and genetic differentiation were explored using haplotype distribution maps and genetic distance matrices. A total of 15 haplotypes were identified in the 32 populations of P. polyphylla var. yunnanensis. Five unique haplotypes were identified from the fourteen haplotypes of the cultivated populations, while only one unique haplotype was identified from the ten haplotypes of the wild populations. The haplotype richness and genetic diversity of the cultivated populations were higher than those of the wild populations (HT = 0.900 vs. 0.861). In addition, there were no statistically significant correlations between geographic distance and genetic distance in the cultivated populations (r = 0.16, p > 0.05), whereas there was a significant correlation between geographical distance and genetic structure in the wild populations (r = 0.32, p > 0.05), indicating that there was a geographical and genetic connection between the wild populations. There was only 2.5% genetic variation between the wild populations and cultivated populations, indicating no obvious genetic differentiation between the wild and cultivated populations. Overall, the genetic background of the cultivated populations was complex, and it was hypothesized that the unique haplotypes and higher diversity of the cultivated populations were caused by the mixed provenance of the cultivated populations.