Demographic history and identification of threats revealed by population genomic analysis provide insights into conservation for an endangered maple

Plant species with extremely small populations conservation program: achieving Kunming-Montreal global biodiversity targets

Conservation programs for plant species with extremely small populations (PSESP) have been successfully implemented for several decades in China. Here we highlight how their inclusion in several national conservation policies helps meet targets of the Kunming-Montreal Global Biodiversity Framework (KMGBF) and show how lessons from these programs can be applied more widely.

Conserving Evolutionary Potential: Combining Landscape Genomics with Established Methods to Inform Plant Conservation

Biodiversity conservation requires conserving evolutionary potential-the capacity for wild populations to adapt. Understanding genetic diversity and evolutionary dynamics is critical for informing conservation decisions that enhance adaptability and persistence under environmental change. We review how emerging landscape genomic methods provide plant conservation programs with insights into evolutionary dynamics, including local adaptation and its environmental drivers. Landscape genomic approaches that explore relationships between genomic variation and environments complement rather than replace established population genomic and common garden approaches for assessing adaptive phenotypic variation, population structure, gene flow, and demography. Collectively, these approaches inform conservation actions, including genetic rescue, maladaptation prediction, and assisted gene flow. The greatest on-the-ground impacts from such studies will be realized when conservation practitioners are actively engaged in research and monitoring. Understanding the evolutionary dynamics shaping the genetic diversity of wild plant populations will inform plant conservation decisions that enhance the adaptability and persistence of species in an uncertain future.

Chromosome-scale genomes of Quercus sichourensis and Quercus rex provide insights into the evolution and adaptation of Fagaceae

The Fagaceae, a plant family with a wide distribution and diverse adaptability, has garnered significant interest as a subject of study in plant speciation and adaptation. Meanwhile, certain Fagaceae species are regarded as highly valuable wood resources due to the exceptional quality of their wood. In this study, we present two high-quality, chromosome-scale genome sequences for Quercus sichourensis (848.75 Mb) and Quercus rex (883.46 Mb). Comparative genomics analysis reveals that the difference in the number of plant disease resistance genes and the nonsynonymous and synonymous substitution ratio (Ka/Ks) of protein-coding genes among Fagaceae species are related to different environmental adaptations. Interestingly, most genes related to starch synthesis in the investigated Quercoideae species are located on a single chromosome, as compared to the outgroup species, Fagus sylvatica. Furthermore, resequencing and population analysis of Q. sichourensis and Q. rex reveal that Q. sichourensis has lower genetic diversity and higher deleterious mutations compared to Q. rex. The high-quality, chromosome-level genomes and the population genomic analysis of the critically endangered Q. sichourensis and Q. rex will provide an invaluable resource as well as insights for future study in these two species, even the genus Quercus, to facilitate their conservation.

Conservation genomics provides insights into genetic resilience and adaptation of the endangered Chinese hazelnut, Corylus chinensis

Global climate change has increased concerns regarding biodiversity loss. However, many key conservation issues still required further research, including demographic history, deleterious mutation load, adaptive evolution, and putative introgression. Here we generated the first chromosome-level genome of the endangered Chinese hazelnut, Corylus chinensis, and compared the genomic signatures with its sympatric widespread C. kwechowensis-C. yunnanensis complex. We found large genome rearrangements across all Corylus species and identified species-specific expanded gene families that may be involved in adaptation. Population genomics revealed that both C. chinensis and the C. kwechowensis-C. yunnanensis complex had diverged into two genetic lineages, forming a consistent pattern of southwestern-northern differentiation. Population size of the narrow southwestern lineages of both species have decreased continuously since the late Miocene, whereas the widespread northern lineages have remained stable (C. chinensis) or have even recovered from population bottlenecks (C. kwechowensis-C. yunnanensis complex) during the Quaternary. Compared with C. kwechowensis-C. yunnanensis complex, C. chinensis showed significantly lower genomic diversity and higher inbreeding level. However, C. chinensis carried significantly fewer deleterious mutations than C. kwechowensis-C. yunnanensis complex, as more effective purging selection reduced the accumulation of homozygous variants. We also detected signals of positive selection and adaptive introgression in different lineages, which facilitated the accumulation of favorable variants and formation of local adaptation. Hence, both types of selection and exogenous introgression could have mitigated inbreeding and facilitated survival and persistence of C. chinensis. Overall, our study provides critical insights into lineage differentiation, local adaptation, and the potential for future recovery of endangered trees.

Transcriptome data analysis provides insights into the conservation of Michelia lacei, a plant species with extremely small populations distributed in Yunnan province, China

BACKGROUND

Michelia lacei W.W.Smith (Magnoliaceae), was classified as a Plant Species with Extremely Small Populations (PSESP) by the Yunnan Provincial Government in both action plans of 2012 and 2021. This evergreen tree is known for its high ornamental and scientific value, but it faces significant threats due to its extremely small population size and narrow geographical distribution. The study aims to understand the genetic structure, diversity, and demographic history of this species to inform its conservation strategies.

RESULTS

The analysis of transcriptome data from 64 individuals across seven populations of M. lacei identified three distinct genetic clusters and generated 104,616 single-nucleotide polymorphisms (SNPs). The KM ex-situ population, originating from Longling County, exhibited unique genetic features, suggesting limited gene flow. The genetic diversity was substantial, with significant differences between populations, particularly between the KM lineage and the OTHER lineage. Demographic history inferred from the data indicated population experienced three significant population declines during glaciations, followed by periods of recovery. We estimated the effective population size (Ne) of the KM and OTHER lineages 1,000 years ago were 85,851 and 416,622, respectively. Gene flow analysis suggested past gene flow between populations, but the KM ex-situ population showed no recent gene flow. A total of 805 outlier SNPs, associated with four environmental factors, suggest potential local adaptation and showcase the species' adaptive potential. Particularly, the BZ displayed 515 adaptive loci, highlighting its strong potential for adaptation within this group.

CONCLUSIONS

The comprehensive genomic analysis of M. lacei provides valuable insights into its genetic background and highlights the urgent need for conservation efforts. The study underscores the importance of ex-situ conservation methods, such as seed collection and vegetative propagation, to safeguard genetic diversity and promote population restoration. The preservation of populations like MC and BZ is crucial for maintaining the species' genetic diversity. In-situ conservation measures, including the establishment of in-situ conservation sites and community engagement, are essential to enhance protection awareness and ensure the long-term survival of this threatened plant species.

Anthropogenic disturbance driving population decline of a dominant tree in East Asia evergreen broadleaved forests over the last 11,000 years

Current biodiversity loss is generally considered to have been caused by anthropogenic disturbance, but it is unclear when anthropogenic activities began to affect biodiversity loss. One hypothesis suggests it began with the Industrial Revolution, whereas others propose that anthropogenic disturbance has been associated with biodiversity decline since the early Holocene. To test these hypotheses, we examined the unique vegetation of evergreen broadleaved forests (EBLFs) in East Asia, where humans have affected landscapes since the early Holocene. We adopted a genomic approach to infer the demographic history of a dominant tree (Litsea elongata) of EBLFs. We used Holocene temperature and anthropogenic disturbance factors to calculate the correlation between these variables and the historical effective population size of L. elongata with Spearman statistics and integrated the maximum-entropy niche model to determine the impact of climate change and anthropogenic disturbance on fluctuation in its effective population size. We identified 9 well-defined geographic clades for the populations of L. elongata. Based on the estimated historical population sizes of these clades, all the populations contracted, indicating persistent population decline over the last 11,000 years. Demographic history of L. elongata and human population change, change in cropland use, and change in irrigated rice area were significantly negatively correlated, whereas climate change in the Holocene was not correlated with demographic history. Our results support the early human impact hypothesis and provide comprehensive evidence that early anthropogenic disturbance may contribute to the current biodiversity crisis in East Asia.

Transcriptome data reveals the conservation genetics of Cypripedium forrestii, a plant species with extremely small populations endemic to Yunnan, China

The Cypripedium forrestii is an orchid species with extremely small populations (PSESP) in Yunnan, China. C. forrestii is range-restricted and less-studied than many orchid species, and it is exposed to various threats to its survival. We investigated its potential habitats and collected 52 samples from eight locations, as well as two outgroup species for reference. We developed genetic markers (SNPs) for C. forrestii based on transcriptome sequencing (RNA-seq) data, and analyzed the genetic diversity, population structure, gene flow and demographic history of C. forrestii in detail. C. forrestii is a taxonomically independent species to protect. We found that the genetic diversity of C. forrestii was very low (1.7e-4) compared with other endangered species. We identified three genetic clusters, and several populations with distinct genetic backgrounds. Most genetic diversity was found within sampling sites (87.87%) and genetic clusters (91.39%). Gene flow has been greatly limited over the most recent generations, probably due to geographical distance, historical climate change and habitat fragmentation. We also detected a severe bottleneck event brought about by the recent population constraints. These factors, together with its reproductive characteristics, contribute to the population fragmentation and low genetic diversity of C. forrestii. Based on our findings, we suggest an integrative conservation strategy to protect and recover the genetic diversity of C. forrestii and a further comprehensive study of its ecological traits in the future.

The chromosome-scale genome of Magnolia sinica (Magnoliaceae) provides insights into the conservation of plant species with extremely small populations (PSESP)

Magnolia sinica (Magnoliaceae) is a highly threatened tree endemic to southeast Yunnan, China. In this study, we generated for the first time a high-quality chromosome-scale genome sequence from M. sinica, by combining Illumina and ONT data with Hi-C scaffolding methods. The final assembled genome size of M. sinica was 1.84 Gb, with a contig N50 of ca. 45 Mb and scaffold N50 of 92 Mb. Identified repeats constituted approximately 57% of the genome, and 43,473 protein-coding genes were predicted. Phylogenetic analysis shows that the magnolias form a sister clade with the eudicots and the order Ceratophyllales, while the monocots are sister to the other core angiosperms. In our study, a total of 21 individuals from the 5 remnant populations of M. sinica, as well as 22 specimens belonging to 8 related Magnoliaceae species, were resequenced. The results showed that M. sinica had higher genetic diversity (θw = 0.01126 and θπ = 0.01158) than other related species in the Magnoliaceae. However, population structure analysis suggested that the genetic differentiation among the 5 M. sinica populations was very low. Analyses of the demographic history of the species using different models consistently revealed that 2 bottleneck events occurred. The contemporary effective population size of M. sinica was estimated to be 10.9. The different patterns of genetic loads (inbreeding and numbers of deleterious mutations) suggested constructive strategies for the conservation of these 5 different populations of M. sinica. Overall, this high-quality genome will be a valuable genomic resource for conservation of M. sinica.

Genome-wide diversity evaluation and core germplasm extraction in ex situ conservation: A case of golden Camellia tunghinensis

Whether ex situ populations constructed in the limited nursery resources of botanical gardens can preserve enough genetic diversity of endangered plants in the wild remains uncertain. Here, a case study was conducted with Camellia tunghinensis, which is one of the species with the lowest natural distribution area in the sect. Chrysantha (golden camellia) of the family Theaceae. We investigated the genetic diversity and population structure of 229 samples from wild and ex situ populations using genotyping by sequencing (GBS). Core germplasm was constructed from these samples. The results showed that wild C. tunghinensis exhibited high genetic diversity, with observed heterozygosity of 0.257-0.293 and expected heterozygosity of 0.247-0.262. Compared with wild populations, the genetic diversity of ex situ populations established by transplanting wild seedlings was close to or even higher. However, the genetic diversity of those established by seed or cuttings of a few superior trees was lower. The Admixture analysis revealed that the structure of the ex situ populations derived from seeds and cuttings was relatively simple compared with the ex situ populations derived from transplanted wild seedlings and wild populations. These results suggested that direct transplanting of wild seedlings was more conducive to preserving the genetic diversity of endangered plants in the wild. In addition, wild populations demonstrated a small differentiation (mean F ST = 0.044) among themselves, possibly due to long-term and frequent gene flow between the wild populations. In contrast, moderate differentiation (mean F ST > 0.05) was detected among ex situ populations and between ex situ and wild populations. This may be the combined result of the absence of gene flow pathways and strong selection pressure in various ex situ environments. Finally, 77 core germplasms were extracted from 229, likely representing the genetic diversity of C. tunghinensis. This study provides future strategies for the ex situ conservation and management of the golden camellia species and other rare and endangered plants.

Seed Germination and Seedling Growth Influenced by Genetic Features and Drought Tolerance in a Critically Endangered Maple

Understanding the adaptation of plant species will help us develop effective breeding programs, guide the collection of germplasm, and improve the success of population restoration projects for threatened species. Genetic features correlate with species adaptation. Acer yangbiense is a critically endangered plant species with extremely small populations (PSESP). However, no information was available on its seed germination and seedling growth in populations with different genetic characteristics. In this study, we investigated seed germination and compared the performance of 566 seedlings in 10 maternal half-sib families cultivated in Kunming Botanical Garden. The results showed that A. yangbiense seeds required an average of 44 days to start germinating, with a 50% germination rate estimated to take about 47-76 days, indicating slow and irregular germination. There is a trade-off between the growth and survival in A. yangbiense seedlings, with fast growth coming at the cost of low survival. Groups that were able to recover from a recent bottleneck consistently had higher relative growth rates. High genetic diversity and low levels of inbreeding are likely to be responsible for their improved survival during drought conditions and rapid growth under optimal environmental conditions. Our results suggest that maternal genetic traits might be used as indicators for conservation and population restoration. These findings provide us with new information that could be applied to support ex situ conservation and reintroduction of threatened species.

The genetic structure and demographic history revealed by whole-genome resequencing provide insights into conservation of critically endangered Artocarpus nanchuanensis

Introduction

Whole-genome resequencing technology covers almost all nucleotide variations in the genome, which makes it possible to carry out conservation genomics research on endangered species at the whole-genome level.

Methods

In this study, based on the whole-genome resequencing data of 101 critically endangered Artocarpus nanchuanensis individuals, we evaluated the genetic diversity and population structure, inferred the demographic history and genetic load, predicted the potential distributions in the past, present and future, and classified conservation units to propose targeted suggestions for the conservation of this critically endangered species.

Results

Whole-genome resequencing for A. nanchuanensis generated approximately 2 Tb of data. Based on abundant mutation sites (25,312,571 single nucleotide polymorphisms sites), we revealed that the average genetic diversity (nucleotide diversity, π) of different populations of A. nanchuanensis was relatively low compared with other trees that have been studied. And we also revealed that the NHZ and QJT populations harboured unique genetic backgrounds and were significantly separated from the other five populations. In addition, positive genetic selective signals, significantly enriched in biological processes related to terpene synthesis, were identified in the NHZ population. The analysis of demographic history of A. nanchuanensis revealed the existence of three genetic bottleneck events. Moreover, abundant genetic loads (48.56% protein-coding genes) were identified in Artocarpus nanchuanensis, especially in genes related to early development and immune function of plants. The predication analysis of suitable habitat areas indicated that the past suitable habitat areas shifted from the north to the south due to global temperature decline. However, in the future, the actual distribution area of A. nanchuanensis will still maintain high suitability.

Discussion

Based on total analyses, we divided the populations of A. nanchuanensis into four conservation units and proposed a number of practical management suggestions for each conservation unit. Overall, our study provides meaningful guidance for the protection of A. nanchuanensis and important insight into conservation genomics research.

Resequencing of Rosa rugosa accessions revealed the history of population dynamics, breed origin, and domestication pathways

BACKGROUND

Rosa rugosa is a shrub that originated in China and has economic and ecological value. However, during the development of R. rugosa, the genetic background was chaotic, and the genetic structure among different wild populations was unclear, as well as wild and cultivated accessions. Here, we report whole-genome resequencing of wild and cultivated R. rugosa accessions.

RESULTS

A total of 19,041,284 SNPs were identified in 188 R. rugosa accessions and 3 R. chinensis accessions by resequencing. Population genetic analysis revealed that cultivated and wild groups were separated very early. All R. rugosa accessions were divided into 8 categories based on genetic structure: (1) Weihai, Yantai, and Liaoning category, (2) Jilin category, and (3) Hammonasset category (above three are wild); (4) traditional varieties, (5) hybrids between R. rugosa and R. chinensis, (6) Zizhi Rose, (7) Kushui Rose, (8) hybrids between R. rugosa and R. multiflora. We found that the heterozygosity and genetic diversity of wild accessions were generally lower than those of cultivated individuals. The genes that were selected during cultivation were identified, and it was found that these genes were mainly related to environmental adaptation and growth.

CONCLUSIONS

The Jilin population was the oldest population and later migrated to Liaoning and then migrated to Yantai and Weihai by sea regression in the Bohai Basin. The Hammonasset naturalized population probably originated from the Jilin population and then experienced separate differentiation. The long-term asexual reproduction pattern of R. rugosa decreased genetic diversity in the wild population. During R. rugosa cultivation, the ancestors of the Jilin population were involved in breeding traditional varieties, after which almost no wild individuals were engaged in breeding. However, in recent decades, cross breeding of R. rugosa started the utilization of wild germplasms. In comparison, some other species play important roles in variety formation. Few genes related to economic traits were selected, suggesting no directional domestication in the R. rugosa cultivation process.

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.