Using genomic information for management planning of an endangered perennial, Viola uliginosa

Multiple introductions of divergent lineages and admixture conferred the high invasiveness in a widespread weed (Hypochaeris radicata)

Biological invasion consists of spatially and temporally varying stages, accompanied by ecological and evolutionary changes. Understanding the genomics underlying invasion dynamics provides critical insights into the geographic sources and genetic diversity, contributing to successful invasions across space and time. Here, we used genomic data and model-based approaches to characterize the invasion dynamics of Hypochaeris radicata L., a noxious weed in Korea. Genetic diversity and assignment patterns were investigated using 3563 SNPs of 283 individuals sampled from 22 populations. We employed a coalescent-based simulation method to estimate demographic changes for each population and inferred colonization history using both phylogenetic and population genetic model-based approaches. Our data suggest that H. radicata has been repeatedly been introduced to Korea from multiple genetic sources within the last 50 years, experiencing weak population bottlenecks followed by subsequent population expansions. These findings highlight the potential for further range expansion, particularly in the presence of human-mediated dispersal. Our study represents the first population-level genomic research documenting the invasion dynamics of the successful worldwide invader, H. radicata, outside of Europe.

ddRAD Sequencing and DNA Barcoding

Double-digest restriction site-associated DNA sequencing is a library preparation protocol that enables capturing variable sites across the genome including single-nucleotide polymorphisms (SNPs). These SNPs can be utilized to gain evolutionary insights into patterns observed in DNA barcodes, to infer population structure and phylogenies, to detect gene flow and introgression, and to perform species delimitation analyses. The protocol includes chemically shearing genomic DNA with restriction enzymes, unique tagging, size selection, and amplification of the resulting DNA fragments. Here we provide a detailed description of each step of the protocol, as well as information on essential equipment and common issues encountered during laboratory work.

RAD-sequencing improves the genetic characterization of a threatened tree peony (Paeonia ludlowii) endemic to China: Implications for conservation

Compared with traditional genetic markers, genomic approaches have proved valuable to the conservation of endangered species. Paeonia ludlowii having rarely and pure yellow flowers, is one of the world's most famous tree peonies. However, only several wild populations remain in the Yarlung Zangbo Valley (Nyingchi and Shannan regions, Xizang) in China due to increasing anthropogenic impact on the natural habitats. We used genome-wide single nucleotide polymorphisms to elucidate the spatial pattern of genetic variation, population structure and demographic history of P. ludlowii from the fragmented region comprising the entire range of this species, aiming to provide a basis for conserving the genetic resources of this species. Unlike genetic uniformity among populations revealed in previous studies, we found low but varied levels of intra-population genetic diversity, in which lower genetic diversity was detected in the population in Shannan region compared to those in Nyingzhi region. These spatial patterns may be likely associated with different population sizes caused by micro-environment differences in these two regions. Additionally, low genetic differentiation among populations (Fst = 0.0037) were detected at the species level. This line of evidence, combined with the result of significant genetic differentiation between the two closest populations and lack of isolation by distance, suggested that shared ancestry among now remnant populations rather than contemporary genetic connectivity resulted in subtle population structure. Demographic inference suggested that P. ludlowii probably experienced a temporal history of sharp population decline during the period of Last Glacial Maximum, and a subsequent bottleneck event resulting from prehistoric human activities on the Qinghai-Tibet Plateau. All these events, together with current habitat fragment and excavation might contribute to the endangered status of P. ludlowii. Our study improved the genetic characterization of the endangered tree peony (P. ludlowii) in China, and these genetic inferences should be considered when making different in situ and ex situ conservation actions for P. ludlowii in this evolutionary hotspot region.

Geographical Patterns of Genetic Variation in Locoto Chile (Capsicum pubescens) in the Americas Inferred by Genome-Wide Data Analysis

The locoto chile (Capsicum pubescens) is a regionally important food crop grown and marketed throughout the mid-highlands of South andCentral America, but little is known about its evolution and the diversity it harbours. An initial scan of genetic diversity and structure across its cultivation range was conducted, the first one using a genomic approach. The RAD-sequencing methodology was applied to a sampling of C. pubescens germplasm consisting of 67 accessions from different American countries, covering its range of distribution/cultivation on the continent. The RAD-seq SNP data obtained clustered the accessions into three major groups, with a high degree of admixture/reticulation among them. Moderate but significant differentiation and geographic structuration were found, depicting a south-north pattern in the distribution of genetic variation. The highest levels of diversity were found among central-western Bolivian individuals, while the lowest was found across Central America-Mexican germplasm. This study provides new genome-wide supported insights into the diversity and differentiation of C. pubescens, as well as a starting point for more efficient use of its genetic variation and germplasm conservation efforts. The findings also contribute to understanding the evolutionary history of C. pubescens, but further investigation is needed to disentangle its origin and diversification under domestication.

Phylogeographic structure of Heteroplexis (Asteraceae), an endangered endemic genus in the limestone karst regions of southern China

Though the karst regions in south and southwest China are plant diversity hotspots, our understanding of the phylogeography and evolutionary history of the plants there remains limited. The genus Heteroplexis (Asteraceae) is one of the typical representative plants isolated by karst habitat islands, and is also an endangered and endemic plant to China. In this study, species-level phylogeographic analysis of the genus Heteroplexis was conducted using restriction site-associated DNA sequencing (RADseq). The genetic structure showed a clear phylogeographic structure consistent with the current species boundaries in the H. microcephala, H. incana, H. vernonioides, H. sericophylla, and H. impressinervia. The significant global (R = 0.37, P < 0.01) and regional (R = 0.650.95, P < 0.05) isolation by distance (IBD) signals among species indicate strong geographic isolation in the karst mountains, which may result in chronically restricted gene flow and increased genetic drift and differentiation. Furthermore, the phylogeographic structure of Heteroplexis suggested a southward migration since the last glacial period. Demographic analysis revealed the karst mountains as a refuge for Heteroplexis species. Finally, both Treemix and ABBA-BABA statistic detected significant historical gene flow between species. Significant historical gene flow and long-term stability of effective population size (Ne) together explain the high genome-wide genetic diversity among species (π = 0.05370.0838). However, the recent collapse of Ne, widespread inbreeding within populations, and restricted contemporary gene flow suggest that Heteroplexis species are probably facing a high risk of genetic diversity loss. Our results help to understand the evolutionary history of karst plants and guide conservation.

Genomic Data Reveals Population Genetic and Demographic History of Magnolia fistulosa (Magnoliaceae), a Plant Species With Extremely Small Populations in Yunnan Province, China

Elucidating the genetic background of threatened species is fundamental to their management and conservation, and investigating the demographic history of these species is helpful in the determination of the threats facing them. The woody species of the genus Magnolia (Magnoliaceae) have high economic, scientific and ecological values. Although nearly half of all Magnolia species have been evaluated as threatened, to date there has been no population genetic study employing Next Generation Sequencing (NGS) technology in this genus. In the present study, we investigate the conservation genomics of Magnolia fistulosa, a threatened species endemic to the limestone area along the Sino-Vietnamese border, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. To increase the reliability of our statistical inferences, we employed two approaches, Stacks and ipyrad, for SNP calling. A total of 15,272 and 18,960, respectively, putatively neutral SNPs were generated by Stacks and ipyrad. Relatively high genetic diversity and large population divergence were detected in M. fistulosa. Although higher absolute values were calculated using the ipyrad data set, the two data sets showed the same trends in genetic diversity (π, H e), population differentiation (F ST) and inbreeding coefficients (F IS). A change in the effective population size of M. fistulosa within the last 1 Ma was detected, including a population decline about 0.5-0.8 Ma ago, a bottleneck event about 0.2-0.3 Ma ago, population fluctuations during the last glacial stage, and the recovery of effective population size after the last glacial maximum. Our findings not only lay the foundation for the future conservation of this species, but also provide new insights into the evolutionary history of the genus Magnolia in southeastern Yunnan, China.

Genetic diversity and structure of Rhododendron meddianum, a plant species with extremely small populations

Rhododendron meddianum is a critically endangered species with important ornamental value and is also a plant species with extremely small populations. In this study, we used double digest restriction-site-associated DNA sequencing (ddRAD) technology to assess the genetic diversity, genetic structure and demographic history of the three extant populations of R. meddianum. Analysis of SNPs indicated that R. meddianum populations have a high genetic diversity (π = 0.0772 ± 0.0024, H _E  = 0.0742 ± 0.002). Both F _ST values (0.1582-0.2388) and AMOVA showed a moderate genetic differentiation among the R. meddianum populations. Meanwhile, STRUCTURE, PCoA and NJ trees indicated that the R. meddianum samples were clustered into three distinct genetic groups. Using the stairway plot, we found that R. meddianum underwent a population bottleneck about 70,000 years ago. Furthermore, demographic models of R. meddianum and its relative, Rhododendron cyanocarpum, revealed that these species diverged about 3.05 (2.21-5.03) million years ago. This divergence may have been caused by environmental changes that occurred after the late Pliocene, e.g., the Asian winter monsoon intensified, leading to a drier climate. Based on these findings, we recommend that R. meddianum be conserved through in situ, ex situ approaches and that its seeds be collected for germplasm.

Conservation Genomics of a Threatened Rhododendron: Contrasting Patterns of Population Structure Revealed From Neutral and Selected SNPs

Though it is well-acknowledged that next generation sequencing (NGS) technologies can provide further insights into plant conservation management than traditional molecular markers, studies employing NGS to address conservation genomics and subsequent conservation strategies for threatened plants are still rare. Rhododendron is the largest genus of woody plants in China, and many species are threatened, however, to date there has been no conservation genetic research using NGS in this genus. In the present study, we investigated the conservation genetics of R. cyanocarpum, a threatened species endemic to the Cangshan Mountains in Yunnan, China, using a double digest restriction-site-associated DNA-sequencing (ddRAD-seq) approach. Due to the availability of sufficient SNPs, we were able to distinguish between neutral and putatively selected SNPs and were able to further investigate the genetic diversity, population structure, and differentiation in R. cyanocarpum, as well as make an estimation of its demographic history. A total of 6,584 SNPs were obtained, of which 5,729 were neutral (detected using Tajima’s D). In terms of the 5,729 neutral SNPs, R. cyanocarpum had a higher genetic diversity (π = 0.0702 ± 0.0017, H_e = 0.0675 ± 0.0016) than other plant species assessed using Rad-seq methods, while population differentiation (F_st from 0.0314 to 0.0452) was weak. Interestingly, contrasting patterns of population structure were revealed from all neutral and selected SNPs, with distinct genetic clusters forming for all SNPs and neutral SNPs, but no distinct subgroups for selected ones. Moreover, we were able to detect changes in effective population size (N_e) of R. cyanocarpum from 150,000 years ago, including a bottleneck event ca. 60,000 years ago, followed by recovery of N_e over a short period, and a subsequent gradual decline in N_e to date. Implications for conserving R. cyanocarpum based on these main results are then discussed.