Next Article in Journal
Vitis rotundifolia Genes Introgressed with RUN1 and RPV1: Poor Recombination and Impact on V. vinifera Berry Transcriptome
Previous Article in Journal
Challenges and Opportunities in the Sustainable Improvement of Carrot Production
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 13
Issue 15
10.3390/plants13152093
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity Analysis and Prediction of Potential Suitable Areas for the Rare and Endangered Wild Plant Henckelia longisepala

by
Renfen Zhao
1,
Nian Huang
1,
Zhiyan Zhang
1,
Wei Luo
1,
Jianying Xiang
1,2,
Yuanjie Xu
3,* and
Yizhi Wang
1,2,*
1
College of Forestry, Southwest Forestry University, Kunming 650224, China
2
Yunnan Academy of Biodiversity, Southwest Forestry University, Kunming 650224, China
3
College of Soil and Water Conservation, Southwest Forestry University, Kunming 650224, China
*
Authors to whom correspondence should be addressed.
Plants 2024, 13(15), 2093; https://doi.org/10.3390/plants13152093
Submission received: 17 June 2024 / Revised: 27 July 2024 / Accepted: 27 July 2024 / Published: 29 July 2024
(This article belongs to the Section Plant Genetics, Genomics and Biotechnology)
Browse Figures
Review Reports Versions Notes

Abstract

:
Henckelia longisepala (H. W. Li) D. J. Middleton & Mich. Möller is a rare and endangered plant species found only in Southeastern Yunnan, China, and Northern Vietnam. It is listed as a threatened species in China and recognized as a plant species with extremely small populations (PSESP), while also having high ornamental value and utilization potential. This study used ddRAD-seq technology to quantify genetic diversity and structure for 32 samples from three extant populations of H. longisepala. The H. longisepala populations were found to have low levels of genetic diversity (Ho = 0.1216, He = 0.1302, Pi = 0.1731, FIS = 0.1456), with greater genetic differentiation observed among populations (FST = 0.3225). As indicated by genetic structure and phylogenetic analyses, samples clustered into three distinct genetic groups that corresponded to geographically separate populations. MaxEnt modeling was used to identify suitable areas for H. longisepala across three time periods and two climate scenarios (SSP1-2.6, SSP5-8.5). High-suitability areas were identified in Southeastern Yunnan Province, Northern Vietnam, and Eastern Laos. Future H. longisepala distribution was predicted to remain centered in these areas, but with a decrease in the total amount of suitable habitat. The present study provides key data on H. longisepala genetic diversity, as well as a theoretical basis for the conservation, development, and utilization of its germplasm resources.

1. Introduction

Gesneriaceae plants are highly appreciated by both domestic and international flower enthusiasts for their diverse morphology, vibrant and colorful flowers, and unique indoor ornamental characteristics [1]. From an early stage, certain developed countries, such as the United States, Canada, and Germany, have engaged in the development and research of Gesneriaceae plants. In addition to widely popular genera like Saintpaulia (commonly known as African violets) and Sinningia, numerous other genera and species have been cultivated, making them highly sought after in the current European and American markets [2]. Beyond their significant ornamental value, Gesneriaceae also hold important medicinal properties, particularly in the southern provinces of China such as Guangxi, Yunnan, and Guizhou. These plants are traditionally used in Chinese medicine for treating bruises, coughs, various sores, and poisons, with notable efficacy. For instance, species like Hemiboea cavaleriei and Hemiboea subcapitata are employed in their respective regions to treat conditions such as scabies and burns [3]. H. longisepala, a member of the Gesneriaceae family, is a subshrub with a stem height ranging from 25 to 80 cm and a stem diameter of 4–6 mm, and a short vertical rhizome in its underground part. Stems simple, apically densely brownish appressed when pubescent, basally glabrescent. Leaves opposite, equal to subequal in a pair; petioles 1.5–5.5 cm long, covered with long appressed brownish hairs when young, glabrescent; leaf blade ovate to oblanceolate, base attenuate to broadly cuneate, margins repand–crenulate to nearly entire, apex acuminate, glabrous to sparsely ciliate, especially abaxial on veins and especially when young, adaxially dark green, abaxially glaucous-green or sometimes dark purple; secondary veins 8–17 pairs; midvein and secondary veins prominent on both surfaces [4]. It has a beautiful plant shape, large and bright flowers, and high ornamental value.
In the wild, H. longisepala is found in the southeastern part of Yunnan Province, China (specifically in Hekou County, Jinping County, and Lüchun County), and Northern Vietnam, where it grows in valleys and shrublands at altitudes ranging from 250 to 800 m. It thrives within banana forests or in shaded areas along ditches. Its habitat requirements are highly specific, demanding stringent environmental conditions. Without a protective canopy from taller vegetation layers, H. longisepala individuals may fail to survive. The uniqueness of its habitat, coupled with its limited range and the ongoing destruction of valley rainforests, pose significant threats to the survival of existing H. longisepala populations. The threatened species list for China’s higher plants [5] lists H. longisepala as an endangered species; hence, its germplasm resources urgently need protection. Currently, one of the better-protected populations of H. longisepala is located in the Huanglianshan National Nature Reserve in Yunnan Province; the other populations are heavily anthropogenically disturbed, with their other habitats fragmented, resulting in lower numbers and poorer growth. The genetic diversity of a species, which is the product of historical evolution, supports its long-term survival and ability to adapt and change [6]. It is well known that preserving the genetic diversity of endangered species can benefit their long-term survival and evolution in changing environments [7,8]. Therefore, it is necessary for researchers to characterize the genetic diversity and population structure of endangered plant species, as a crucial step for their conservation and management [9].
Molecular marker technology is an important means of evaluating the genetic diversity of germplasm resources [10], whether of different populations of organisms or of different individuals from the same population [11,12]. Double digest restriction-site associated sequencing (ddRAD-seq) is a molecular marker technology developed by researchers to solve the problems of a relatively complicated library construction, high sequence dispersion, and high operation requirements [13,14,15,16]. ddRAD-seq utilizes multiplexed libraries obtained through the enzymatic digestion of whole genomic DNA followed by binding to specific adapters (creating reduced-representation libraries); it enables the acquisition of more SNP loci for a competitive cost [17]. ddRAD-seq technology not only improves sequencing efficiency, but also reduces the cost of experimentation [18]. Recently, ddRAD-seq technology has been widely applied in phylogenetics [19], molecular breeding [20], and population genetics [21], among other fields. Various enzymes suitable for ddRAD library construction in angiosperms have been described, enabling the rapid development of SNP markers and genome typing for multiple angiosperm species [22]. In addition, 18 QTLs (quantitative trait loci) closely related to yield traits in Asparagus officinalis were successfully identified using ddRAD, providing a solid theoretical foundation for the genetic analysis of yield traits and molecular marker-assisted breeding in Asparagus officinalis [23]. Furthermore, a study of Rhododendron meddianum with ddRAD technology revealed a high level of genetic diversity and moderate genetic differentiation among populations [6]. More recently, various studies have demonstrated that the ddRAD-seq technique can accurately determine SNP genotypes [24]. However, ddRAD-seq has not yet been utilized to study H. longisepala genetic diversity and genetic structure; to date, it has also not been used for any plant species in the same family (Gesneriaceae).
Owing to ongoing habitat destruction, H. longisepala populations are small and narrowly distributed. Understanding climate change impacts on H. longisepala distribution could provide valuable information for the conservation and management of this rare and endangered species. Climatic conditions are one of the most important factors affecting the natural geographical distributions of plant species, and climate change may lead to shifts in species’ distributions and also significantly impact their diversity [25]. Species distribution models are widely used in biogeography, conservation biology, and ecology because they can accurately predict the suitable habitat range for a given species [26]. The five most commonly used species distribution models are BIOCLIM, CLIMEX, ENFA, GARP, and MaxEnt [27]. The MaxEnt model is constructed based on known distribution points and associated environmental data, and then used to predict suitable areas based on relevant environmental variables [28]. Compared with other species distribution models, MaxEnt is simple to operate and can effectively predict suitable areas for narrowly distributed species with small sample sizes [29].
In this study, ddRAD-seq technology was used to study the genetic diversity and structure of H. longisepala, with the goal of better understanding its evolutionary potential and providing key reference data for its conservation. Using bioclimatic data and field distribution data for H. longisepala, MaxEnt modeling and ArcGIS software were used to predict the potential H. longisepala distribution under different climate scenarios. Key climatic factors affecting its distribution were identified to provide a theoretical basis for the conservation of H. longisepala.

2. Materials and Methods

2.1. Plant Materials and Data Collection

A total of 32 individual samples of H. longisepala were collected in January 2024 from three wild populations in the Hekou, Jinping, and Lüchun counties of Honghe Prefecture, Yunnan Province (Table 1). Several individuals were sampled per population, and one or two healthy, young leaves were selected per individual, and the distance between sampling points of the same populations was 50 m (between populations, it was more than 70 km). After collection, the leaves were quickly dried and preserved with silica gel. Further information on the sampling sites, number of samples, and population codes is provided in Table 1.
Distribution data were obtained from the Chinese Virtual Herbarium (https://www.cvh.ac.cn/), the Global Biodiversity Information Facility (https://www.gbif.org/zh/), the published literature, local flora guides, and field investigations. A total of 13 distribution records were obtained for H. longisepala after removing repetitive and unclear datapoints. Climatic data (bio1-19) were downloaded from the WorldClim database (https://worldclim.org/data/index.html, accessed on 13 April 2024); the data included 19 climate variables with a resolution of 2.5′ (about 5 km geographical resolution) for the current period (1970–2000), the 2030s (2021–2040), and the 2050s (2041–2060). Climate data for both future periods were predicted using the climate system model BCCCSM2-MR of The Beijing Climate Center (BCC) in the Global Climate Modeling (GCM) data block, a reliable source of future climate distribution data. For example, correlations between the model predictions and observed values were as high as 0.86 [28]. For the future climate data, there were four greenhouse gas emission scenarios: SSP1-2.6 (low-concentration CO2 emission scenario), SSP2-4.5 (low- and medium-concentration CO2 emission scenario), SSP3-7.0 (medium- and high-concentration CO2 emission scenario), and SSP5-8.5 (high-concentration CO2 emission scenario). Two greenhouse gas emission scenarios (SSP1-2.6, SSP5-8.5) were utilized in this study.

2.2. Genetic Diversity Assessment

2.2.1. DNA Extraction, Library Construction, and Sequencing

Total genomic DNA was extracted from silica-dried leaf tissues using a modified CTAB method [30]. The quality of the extracted DNA was assessed via 0.8% agarose gel electrophoresis, and DNA samples were quantified using an ultraviolet spectrophotometer [31], thereby ensuring compliance with subsequent library-building requirements. A combination of DpnII and MspI endonucleases was used for the double digestion of the quality-inspected DNA samples. Based on the selected endonuclease, the most suitable reaction conditions were determined, typically setting the digestion temperature at 37 °C for 5 h. After enzyme digestion, 5 µL of each sample was taken for gel electrophoresis to verify digestion. The DNA enzyme digestion was required to be complete, meaning no main band should be visible post-digestion. Instead, there should be clearly dispersed bands between 250 and 1000 bp. VAHTSTM DNA Clean Beads (Vazyme, Nanjing, China) were then used to remove over- and under-sized fragments after digestion. Subsequently, the sheared fragments were ligated to the P2 adapter. A PCR was performed to enrich the sequencing library template using the DNA fragments connected to the adapters. Finally, the PCR products were purified using VAHTSTM DNA Clean Beads. Final fragment selection and purification of the DNA library was performed using 2% agarose gel electrophoresis (220~450 bp). After passing quality control, the libraries were paired-end sequenced on the Illumina NovaSeq platform by Shanghai Personalbio Co. (Shanghai, China).

2.2.2. Sequencing Data Quality Control and Filtering

After sequencing, low-quality reads were removed using Fastp (v0.20.0) [32] to obtain high-quality sequencing data. The filtering criteria included the following: (1) The removal of adapter sequences, specifically at the 3′ end. (2) Quality filtering, using a sliding window approach, with a window size of 5 bp. Windows moved from the 3′ end to the 5′ end, and the average Q-value was calculated for each position; bases in the window were deleted if the Q-value was less than 20. (3) Length filtering: reads ≤50 bp in length were removed. (4) Fuzzy base N-filtering: reads with more than five N bases were removed.

2.2.3. SNP Detection and Screening

As there is no published genome for the Gesneriaceae plant family, no reference genome was available for this study. Therefore, reads in each sample were clustered according to sequence similarity using ustacks (http://catchenlab.life.illinois.edu/stacks/, accessed on 26 April 2024) [33], and loci were combined across all samples using cstacks to obtain consensus sequences for each loci. SNPs were obtained using the populations program in Stacks. Finally, filtering in PLINK (v1.90b6) [34] was implemented to obtain high-quality SNPs (parameters: --geno 0.05 --mind 0.1 --maf 0.05).

2.2.4. Data Analysis

The SNP data were analyzed with the following software. GCTA (v1.94.1) [35] was used to perform a principal component analysis based on the dataset of high-quality SNPs. The populations module of Stacks (v2.58) [33] was utilized to calculate genetic diversity parameters, including the expected heterozygosity (He), observed heterozygosity (Ho), the inbreeding coefficient (Fis), and nucleotide diversity (Pi). Admixture (v1.3.0) [36] was used to assess the population genetic structure with values of K ranging from two to ten (i.e., assuming the existence of 2–10 ancestral populations). A mixed model approach was used in Admixture, and the CV error was obtained for each value of K. The maximum likelihood algorithm in FastTree (2.1.11) [37] was used to construct a phylogeny for the H. longisepala samples (FastTree -gtr -nt alignment file --boot 1000 > tree file).

2.3. Suitable Area Prediction for H. longisepala

2.3.1. Data Preprocessing

Distribution data for H. longisepala were processed in preparation for subsequent analyses. The longitude and latitude values for the distribution sites were converted into decimal degree values, input into Excel tables, and saved in ASC format for later use.
After downloading the 19 climatic factors from Worldclim (Table 2), the prediction area was set to include Yunnan and the surrounding areas of Guangxi, Vietnam, and Laos. Climate data for the prediction area were obtained via mask extraction and saved in ASC format.

2.3.2. Correlation Analysis of Climatic Factors

As many climate variables are strongly correlated, their joint inclusion in statistical models can affect model outcomes, such as assessments of variable contributions and observed relationships among variables, thereby affecting the accuracy of the model fitting [38]. The bioclimatic data for each period included 19 climatic factors related to temperature and precipitation (Table 2). Therefore, an initial model was run in MaxEnt to evaluate the contributions of each climatic factor to the model predictions. Climatic factors with a contribution rate < 0% were removed in order of their rank. Next, a correlation analysis was performed in ENMTools (https://rdocumentation.org/packages/ENMTools/, accessed on 2 May 2024) for the remaining climatic factors. Highly correlated (|r| > 0.8) factors were grouped into pairs, and the variable with the lowest contribution rate was eliminated from each pair.

2.3.3. Parameter Selection

The H. longisepala climate and distribution data were imported into MaxEnt, where response curves were drawn for each climatic factor, predicted distribution maps were created, and a variable importance analysis was conducted using the jackknife method. The model was run ten times, with 25% of the data points included in the validation set. The accuracy of the model predictions was verified by comparing the actual distribution points in the validation set with the predicted values. The receiver operating characteristic (ROC) curve generated by the MaxEnt (https://biodiversityinformatics.amnh.org/open_source/maxent/, accessed on 1 May 2024) model was used to evaluate the model accuracy, with the area under the curve (AUC) being used to measure the impact of the climatic variables on H. longisepala distribution predictions. The larger the AUC value, the greater the influence of a given variable, indicating the greater accuracy of the prediction results. The AUC ranges from zero to one, with AUC values from 0.5 to 0.6 indicating a failed prediction, 0.6 to 0.7 indicating poor prediction results, 0.7 to 0.9 good prediction results, and AUC values greater than 0.9 indicating a high level of confidence in the predictions [39].

2.3.4. Classification and Description of Suitable Areas

The ArcGIS (https://www.esri.com/en-us/arcgis/products/arcgis-desktop/resources/, versions 10.7, accessed on 26 May 2024) reclassification function was used to assess suitability for H. longisepala based on the MaxEnt simulation results. Four suitability classes were established based on probability values, and the range of suitability values was determined via the equal distance method [40]. Non-suitable areas ranged from 0 to 0.2142, low-suitability areas from 0.2142 to 0.4284, moderate-suitability areas from 0.4284 to 0.6426, and high-suitability areas from 0.6426 to 1. The proportion of each study site suitable for H. longisepala was determined using a regional statistical function.

3. Results

3.1. ddRAD Sequencing and Data Processing

Double digest restriction-site associated (ddRAD) sequencing was performed for 32 H. longisepala individuals, and a total of 607,138,062 high-quality reads (or 18,973,064 per sample) were obtained after quality control (Table 3). High-quality reads accounted for more than 96% of all reads per sample, with an average of 97.96%. The average GC content per sample was 40.33%, and the average Q30 was 97.29%. The low GC content and high Q30 suggest that the library construction was successful, and the amount of data obtained met the requirements for subsequent data analyses.

3.2. SNP Statistics

SNPs were screened and filtered using PLINK [31], and a total of 1,613,563 SNPs were obtained, with an average of 50,423 valid SNPs per sample. More SNPs were heterozygous than homozygous: there were 822,245 heterozygous SNPs (or 52.87% heterozygous sites on average) and 791,318 homozygous SNPs (47.13%). Transitions or transversions occurred at 1,613,563 SNPs, with transitions being more common than transversions. The ratio of transitions to transversions (Ts/Tv) ranged from 1.31 to 1.48, with an average of 1.37. There were 934,743 SNPs with purine-to-purine or pyrimidine-to-pyrimidine transitions, with a transition ratio of 57.93%, and 678,820 SNPs with purine-to-pyrimidine substitutions, with a transversion ratio of 42.07% (Table 4).

3.3. Genetic Diversity Analysis

Heterozygosity is an important index of genetic variation in natural populations, which can reveal the population structure and even historical dynamics [41]. The observed heterozygosity (Ho) of the three H. longisepala populations ranged from 0.1182 to 0.1256 (mean Ho = 0.1216), and the expected heterozygosity (He) ranged from 0.1151 to 0.1883 (mean He = 0.1302) (Table 5). In the HSL and TMC populations, the observed heterozygosity was lower than the expected heterozygosity, indicating a history of self-fertilization in these populations. In the NXC population, the observed heterozygosity was greater than the expected heterozygosity, suggesting a heterozygote advantage. The inbreeding coefficient (FIS) ranged from 0.0420 to 0.207, with an average value of 0.1456 (FIS > 0.78%); FIS was highest in HLS, followed by TMC and NXC. This suggests a high level of inbreeding in some populations. Nucleotide diversity (Pi) is a measure of within-population genetic diversity and reflects the average variation between samples [42]. The nucleotide diversity of the three populations ranged from 0.1408 to 0.1990, with an average value of 0.1731. Among populations, Pi was highest in HLS; in other words, differences among HLS samples were more pronounced than for samples from other populations.
The population differentiation index (FST) quantifies changes in population structure, and the magnitude of FST is correlated with the level of genetic differentiation, which may be affected by mutation, drift, and natural selection, among other factors [43,44]. The FST was calculated for each population pair (Table 6), with values ranging from 0.1924 to 0.4071 and a mean of 0.3225. While the HLS and TMC populations were well-differentiated (FST of 0.1924), genetic differentiation was highest between NXC and the other two populations (FST values greater than 0.25). The genetic distance was calculated between populations of H. longisepala (Table 7), with values ranging from 0.4524 to 0.8276 and a mean value of 0.6824. While the genetic distance between the HLS and TMC populations was the highest (0.8276), the genetic distance between the HLS and TMC population was the smallest (0.4524). The genetic distance between different samples ranged from 0.0072 to 0.1869 (S1), with an average value of 0.0795. The genetic distance between the TM-B6 and TM-C2 samples was the smallest, while the genetic distance between the NX-3 and TM-C6 samples was the highest (Table S1).

3.4. Genetic Structure Analysis

Phylogenetic trees can be used to show the evolutionary relationships among populations within a species and to measure how closely related the populations are [45]. As can be seen from the phylogenetic tree for H. longisepala (Figure 1), the 32 samples were classified into three clusters, one for each population, and there was no evidence of intermingling between populations.
In population genetic studies, SNP data can be used to classify individuals into subgroups via principal component analysis [42]. In this study, a principal component analysis (PCA) was performed on all 32 samples (Figure 2), with the first principal component (PC1) and second principal component (PC2) accounting for 27% and 16.7%, respectively, of the variation. Together, PC1 and PC2 categorized the samples into three groups, one for each population, which were widely spaced in the PCA diagram and genetically distinct.
Characterizing the population genetic structure can reveal ongoing evolutionary processes within a species [46]. In population structure analyses, CV error values are compared for different values of K; the smaller the CV error, the more accurate the K value. As shown in Figure 3, the CV error was minimized for K = 3, indicating that the optimal number of clusters for the study samples is three. When K = 2, NXC and HLS were grouped together, separately from TMC, which formed its own group. When K = 3, the three populations each represented their own group (Figure 4), although there was evidence of gene exchange between some HLS and TMC samples. The genetic structure results are consistent with those of the phylogenetic tree and principal component analysis.

3.5. MaxEnt Model Prediction and Major Climate Factors

Using MaxEnt, AUC values were obtained to monitor the simulation progress over ten repetitions. Under current climate conditions, the training set AUC for the potential distribution model was 0.974. Under the SSP1-2.6 scenario, the average AUC of the training set was 0.922 for the 2030s and 0.903 for the 2050s. Under the SSP5-8.5 scenario, the average AUC of the training set measured 0.925 in the 2030s and 0.919 in the 2050s. According to standard AUC evaluation criteria, the prediction model therefore achieved a high level of accuracy, effectively predicting the H. longisepala potential distribution with accurate and reliable simulation results.
Through the initial modeling and correlation analysis, this study ultimately identified seven bioclimatic factors suitable for constructing a MaxEnt model (Table 8). Among these seven factors, bio18 (precipitation of warmest quarter), bio7 (annual temperature range), bio16 (precipitation of wettest quarter), and bio9 (mean temperature of driest quarter) made significant contributions, with a cumulative contribution rate of 91.5%. Precipitation in the warmest quarter (bio18) had the greatest impact on the plant, with a contribution of 36.7%.

3.6. Current Potential Distribution of H. longisepala

Species distribution data generated by the MaxEnt model were imported into ArcGIS to create a geographical distribution map of suitable areas for H. longisepala under current climatic conditions (Figure 5). Predicted values were classified into four habitat suitability groups: non-suitable, low suitability, moderate suitability, and high suitability. The amount of habitat encompassed by each group was calculated. The potential H. longisepala habitat covered about 32.42 × 104 km2, primarily in Southeastern Yunnan Province, Western Guangxi Province, Northern Vietnam, and Eastern Laos. Within this area, low-suitability habitat covered 15.67 × 104 km2 (or 48.33% of the total suitable area), moderate-suitability habitat covered 11.25 × 104 km2, and high-suitability habitat covered 5.50 × 104 km2.

3.7. Potential Future Distribution of H. longisepala

Climate values for the 2030s and 2050s, under either the SSP1-2.6 or SSP5-8.5 scenario, were used to predict the future distribution of H. longisepala (Figure 6). Under the SSP1-2.6 scenario, the total suitable area covered 26.63 × 104 km2 from 2021 to 2040, which is 17.86% less than the current habitat size (Table 9). Low-suitability habitat covered 14.12 × 104 km2, moderate-suitability habitat 7.26 × 104 km2, and high-suitability habitat 5.24 × 104 km2; these values represent decreases of 9.89%, 35.47%, and 4.73%, respectively, from current values. From 2041 to 2060, suitable H. longisepala habitat covered 27.79 × 104 km2, a decrease of 14.28% compared with the current value. Low-suitability habitat covered 14.93 × 104 km2 and moderate-suitability habitat covered 7.76 × 104 km2, a decrease of 31.02%. High-suitability habitat covered 5.10 × 104 km2, a decrease of 7.27%.
Under the SSP5-8.5 scenario, the total H. longisepala habitat covered 19.78 × 104 km2 from 2021 to 2040, which is 38.99% less than the current habitat size. Low-suitability habitat covered 10.76 × 104 km2, moderate-suitability habitat 4.99 × 104 km2, high-suitability habitat 4.03 × 104 km2; these represent reductions of 31.33%, 55.64%, and 26.73%, respectively, from current values. From 2041 to 2060, the total habitat area was 25.27 × 104 km2, which is 22.05% less than the current habitat size. Low-suitability habitat covered 12.63 × 104 km2, moderate-suitability habitat 7.66 × 104 km2, and high-suitability habitat 4.98 × 104 km2, representing decreases of 19.40%, 31.91%, and 9.45%, respectively, from current values. For both future emission scenarios, the total amount of suitable habitat for H. longisepala decreased, with suitable areas primarily located in Southeastern Yunnan Province, Northern and Western Vietnam, Eastern Laos, and Western Guangxi Province. The amount of high-suitability habitat also decreased, indicating that the future climate is predicted to have a negative impact on H. longisepala growth.

4. Discussion

4.1. Genetic Diversity

Breeding and improvement work depend upon adequate germplasm resources, but the effective utilization of germplasm resources relies on a characterization of their genetic diversity and other background studies [47]. In this study, genetic variation within and among populations of H. longisepala was investigated for the first time using ddRAD-seq. The genetic diversity within a species is influenced by a number of factors, including its breeding system, evolutionary history, geographic range, life history, and mode of seed dispersal [7]. Within populations, the observed heterozygosity averaged 0.1216, which was lower than the expected heterozygosity (0.1302), indicating the occurrence of self-fertilization within populations. Genetic diversity in H. longisepala (0.1731) was lower than that of Ramonda myconi [48], as assessed with RAPD (H = 0.2591, I = 0.3883). To date, only older molecular markers such as AFLPs and ISSRs have been utilized for the study of genetic diversity in the Gesneriaceae. For example, ISSRs were used to quantify genetic diversity in Oreocharis benthamii (H = 0.1014, I = 0.1528) [49], and AFLPs were utilized in Primulina tabacum (H = 0.220, I = 0.32) [50]. In general, the genetic diversity of endangered and very small populations of plants is low, as reported for Vatica xishuangbannaensis (0.169; RAPD markers) [51], Bretschneidera sinensis (0.141; ISSR markers) [52], and Horsfieldia pandurifolia (0.152; AFLPs) [53]. The low level of genetic diversity observed here for H. longisepala may be attributed to its highly specific habitat requirements, restricted pollen movement, and limited seed dispersal. According to the field investigation, the natural distribution of H. longisepala is narrow, with plants primarily located alongside highways, either under plantain forests or in the shade of a ditch. Seeds are oval, about 0.25 mm long, and frequently do not germinate after dispersal, such as when they fall into thick litter. Therefore, H. longisepala seedlings are few and mostly distributed along roads. These findings may also be due to the size of the study sample and its geographical attributes.

4.2. Genetic Structure

In this study, the principal component analysis, population genetic structure analysis, and phylogenetic analysis revealed that the three H. longisepala populations (HLS, TMC, NXC) were genetically separated concordant with their geographical distribution; little genetic exchange occurred among the sample populations during their evolutionary history. This pattern may be due to the limited number of individuals sampled and the small sampling range; for example, only four samples were collected in NXC. It may also be due to the narrow natural distribution of H. longisepala and obvious geographical isolation among populations. In addition, the FST among populations ranged from 0.1924 to 0.4071, with a mean value of 0.3225, indicating a high level of genetic differentiation among populations. The geographical distance between the study populations was more than 50 km, making it difficult for gene flow (as seeds or pollen) to occur between populations. Furthermore, according to the field surveys, the H. longisepala habitat was found to be severely disturbed by human activities. When the gene flow is low, it can produce genetic differentiation among populations [54]. The analyses in this study were complementary, and together indicated that the observed population structure was a reliable finding.

4.3. Prediction of Potential Suitable Areas

The geographic distribution of a given species is determined by a variety of factors, including anthropological, biological, and environmental factors, as well as the species’ own ability to evolve, spread, and adapt [55]. Climate is considered to be the most important factor shaping species’ distributions [56,57]. As documented here, the H. longisepala distribution is concentrated in the Southeastern Yunnan Province and Northern Vietnam. The predicted distribution indicated a wider habitat range for H. longisepala than that documented, and in addition to the existing records, there are potentially suitable areas for this species in Western Guangxi, as well as in Laos. According to the modeling results, precipitation in the warmest quarter, annual temperature range, precipitation in the wettest quarter, and the mean temperature of the driest quarter were the most important factors shaping H. longisepala distribution, indicating that rainfall and temperature were limiting factors in determining habitat suitability. H. longisepala distribution falls within a subtropical monsoon climate with high rainfall and temperatures, consistent with the observed H. longisepala distribution described here.
The prediction accuracy of the MaxEnt modeling depends on both the actual distribution data and environmental factors [58]. As a rare and endangered species, there are few reliably documented H. longisepala sites in the wild. Human activity also substantial affects H. longisepala survival, but no metric of human disturbance was included in this study, which may have affected the accuracy of the simulation results. At present, H. longisepala only exists in fragmented habitat patches in the wild, the availability of suitable habitat is continuously declining, and there is no current management scheme to introduce seedlings to new areas. Therefore, using study predictions, H. longisepala could be introduced to high-suitability areas as part of a conservation scheme. Future studies should examine aspects of H. longisepala physiology and biochemistry, tissue culture techniques, and relevant ecological characteristics.

4.4. Conservation Implications

The ultimate goal of species’ conservation is to ensure the continuous survival of populations and to maintain their evolutionary potential by preserving natural levels of genetic diversity [59,60]. The conservation status of H. longisepala in China has been determined in the present study by understanding the population structure and genetic variation between and within different populations. An in situ conservation plan should be implemented, with a prime focus on protecting the core areas with vulnerable populations, to preserve and maintain the prevailing genetic diversity. The ex situ approach of conservation should involve the establishment of a germplasm or seed bank for H. longisepala. However, priority should be given to the populations harboring higher genetic diversity for seed bank development. Fortunately, the HLS population is located inside a national nature reserve, which is less disturbed by anthropogenic activities, has better growth conditions, and has larger populations. From the field survey, the natural distribution of H. longisepala is narrow, and its habitat mostly consists of highly disturbed areas along highways. Hence, conservation of the H. longisepala habitat is key to preventing further losses of genetic variation. Conservation strategies should consider genetic aspects, in view of the need to reduce genetic drift and inbreeding [61]. Conservation measures should also aim to reduce anthropogenic impacts on H. longisepala survival, as well as on other Gesneriaceae species; species in this family are typically characterized by narrow distributions, small populations, and specialized growth conditions. Considering the low genetic diversity and significant genetic differentiation in H. longisepala, we suggest that in situ conservation should be established in all existing populations (especially the NXC population) to protect the natural habitats and individuals of H. longisepala. However, the genetic diversity and numbers of individuals of H. longisepala cannot be increased by short-term in situ conservation. We should implement artificial breeding, to transplant large numbers of seedlings to the natural populations to increase the size of the population. Then, introductions might be carried out in suitable areas, as identified here, with future studies of plant growth and the development needed to support population recovery. Apart from the conservation strategies mentioned above, the sensitization of local people about the alarming genetic situation of this important endangered H. longisepala, and involving local bodies in framing policies and programs for proper plant management, are also highly recommended.

5. Conclusions

In this study, the genetic diversity and genetic structure of 32 individual H. longisepala samples were investigated for the first time using ddRAD-seq technology. Low levels of genetic diversity were found in the H. longisepala populations, with the greatest variation among samples occurring in the HLS population; more genetic differentiation existed among populations than within populations. Using population structure analyses and PCA, samples were clustered into three distinct genetic groups. MaxEnt modeling was implemented to determine the potential H. longisepala distribution under two climate scenarios (SSP1-2.6 and SSP5-8.5) during the current period, 2021–2040, and 2041–2060. Bio18 (precipitation of warmest quarter), bio7 (annual temperature range), bio16 (precipitation of wettest quarter), and bio9 (mean temperature of driest quarter) were the main factors affecting H. longisepala distribution, with a cumulative contribution rate of 91.5%. Precipitation in the warmest quarter (bio18) had the greatest impact on the plant, with a contribution of 36.7%. At present, high-suitability areas for H. longisepala are mostly found in Southeastern Yunnan Province, Northern and Western Vietnam, Eastern Laos, and Western Guangxi Province. In the future, the potential H. longisepala distribution is predicted to remain centered in these areas, but the amount of suitable habitat will be reduced, especially for high-suitability habitat.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants13152093/s1, Table S1: The table of genetic distances between the samples.

Author Contributions

Conceptualization, R.Z. and Y.W.; methodology, R.Z.; investigation, N.H., Z.Z., W.L. and Y.W.; writing—original draft preparation, R.Z. and J.X.; writing—review and editing, Y.X. and Y.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a grant from the Yunnan Key Laboratory of Forest Plant Cultivation and Utilization/State Forestry Administration Key Laboratory of Yunnan Rare and Endangered Species Conservation and Propagation (KFJJ21-02); the Agricultural Joint Special Project of Yunnan Province (202301BD070001-122); and the Ten Thousand Talents Program of Yunnan Province (Top-notch Young Talents, YNWR-QNBJ-2019-244).

Data Availability Statement

The original contributions presented in the study are included in the article and Supplementary Material; further inquiries can be directed to the corresponding author.

Acknowledgments

We thank Bo Xiao, Jinguo Zhang, Zhengming Zhu, and Zhiyong Yu for their assistance in the field investigation.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Tatsuzawa, F.; Hosokawa, M. Flower Colors and Their Anthocyanins in Saintpaulia Cultivars (Gesneriaceae). Hortic. J. 2016, 85, 63–69. [Google Scholar] [CrossRef]
  2. Jiang, N.; Ning, S.J.; Pan, B.; Shen, X.L.; Huang, S.J. The Resources of Ornamental Plants in Gesneriaceae in Guangxi. Guihaia 2012, 32, 494–500. [Google Scholar]
  3. Wang, N.; Hu, C.Y.; Xiao, P.G.; Liu, Y.; Bai, Z.F. A Survey of the Traditional Pharmacognosy of Domestic Plants of the Family Gesneriaceae. J. Chin. Med. Mater. 2023, 46, 2975–2979. [Google Scholar]
  4. Bui, H.Q.; Nuraliev, M.S.; Möller, M.; Kuznetsov, A.N.; Kuznetsova, S.P.; Middleton, D.J.; Wen, F. Henckelia Longisepala (Gesneriaceae), a New Record for Vietnam. Rheedea 2020, 30, 176–186. [Google Scholar] [CrossRef]
  5. Qin, H.N.; Yang, Y.; Dong, S.Y.; He, Q.; Jia, Y.; Zhao, L.; Yu, S.X.; Liu, H.Y.; Liu, B.; Yan, Y.H.; et al. Threatened Species List of China’s Higher Plants. Biodivers. Sci. 2017, 25, 696–744. [Google Scholar] [CrossRef]
  6. Zhang, X.J.; Liu, X.F.; Liu, D.T.; Cao, Y.R.; Li, Z.H.; Ma, Y.P.; Ma, H. Genetic Diversity and Structure of Rhododendron Meddianum, a Plant Species with Extremely Small Populations. Plant Divers. 2021, 43, 472–479. [Google Scholar] [CrossRef] [PubMed]
  7. Wu, F.Q.; Shen, S.K.; Zhang, X.J.; Wang, Y.H.; Sun, W.B. Genetic Diversity and Population Structure of an Extremely Endangered Species: The World’s Largest Rhododendron. AoB Plants 2015, 7, plu082. [Google Scholar] [CrossRef] [PubMed]
  8. Nielsen, N.H.; Backes, G.; Stougaard, J.; Andersen, S.U.; Jahoor, A. Genetic Diversity and Population Structure Analysis of European Hexaploid Bread Wheat (Triticum aestivum L.) Varieties. PLoS ONE 2014, 9, e94000. [Google Scholar] [CrossRef] [PubMed]
  9. Guizado, S.J.V.; Nadeem, M.A.; Ali, F.; Barut, M.; Habyarimana, E.; Gómez, T.P.; Santillan, J.A.V.; Canales, E.T.; Gómez, J.C.C.; Chung, G.; et al. Genetic Diversity and Population Structure of Endangered Rosewood from the Peruvian Amazon Using ISSR Markers. Acta Amaz. 2020, 50, 204–212. [Google Scholar] [CrossRef]
  10. Cao, Y.Y.; Diao, Q.N.; Chen, Y.Y.; Zhang, Y.P. Analysis of Genetic Diversity of Cucumis Melo Based on 2b-RAD Simplified Genome Sequencing. Acta Bot. Boreali-Occident. Sin. 2021, 41, 96–106. [Google Scholar]
  11. Lavretsky, P.; DaCosta, J.M.; Sorenson, M.D.; McCracken, K.G.; Peters, J.L. ddRAD-seq Data Reveal Significant Genome-wide Population Structure and Divergent Genomic Regions That Distinguish the Mallard and Close Relatives in North America. Mol. Ecol. 2019, 28, 2594–2609. [Google Scholar] [CrossRef] [PubMed]
  12. Garcia, C.B.; Silva, A.V.D.; Carvalho, I.A.S.D.; Nascimento, W.F.D.; Ramos, S.L.F.; Rodrigues, D.P.; Zucchi, M.I.; Costa, F.M.; Alves-Pereira, A.; Batista, C.E.D.A.; et al. Low Diversity and High Genetic Structure for Platonia Insignis Mart., an Endangered Fruit Tree Species. Plants 2024, 13, 1033. [Google Scholar] [CrossRef] [PubMed]
  13. Yangkun, W.; Yan, H.; Tianzhen, Z. Current status and perspective of RAD-seq in genomic research. Hereditas 2014, 36, 41–49. [Google Scholar] [CrossRef] [PubMed]
  14. Mu, X.Y.; Tong, L.; Sun, M.; Zhu, Y.X.; Wen, J.; Lin, Q.W.; Liu, B. Phylogeny and Divergence Time Estimation of the Walnut Family (Juglandaceae) Based on Nuclear RAD-Seq and Chloroplast Genome Data. Mol. Phylogenet. Evol. 2020, 147, 106802. [Google Scholar] [CrossRef]
  15. Kajiya-Kanegae, H.; Takanashi, H.; Fujimoto, M.; Ishimori, M.; Ohnishi, N.; Wacera, W.F.; Omollo, E.A.; Kobayashi, M.; Yano, K.; Nakano, M.; et al. RAD-Seq-Based High-Density Linkage Map Construction and QTL Mapping of Biomass-Related Traits in Sorghum Using the Japanese Landrace Takakibi NOG. Plant Cell Physiol. 2020, 61, 1262–1272. [Google Scholar] [CrossRef] [PubMed]
  16. Peterson, B.K.; Weber, J.N.; Kay, E.H.; Fisher, H.S.; Hoekstra, H.E. Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species. PLoS ONE 2012, 7, e37135. [Google Scholar] [CrossRef] [PubMed]
  17. Aballay, M.M.; Aguirre, N.C.; Filippi, C.V.; Valentini, G.H.; Sánchez, G. Fine-Tuning the Performance of ddRAD-Seq in the Peach Genome. Sci. Rep. 2021, 11, 6298. [Google Scholar] [CrossRef] [PubMed]
  18. Li, L.L.; Huang, S.; Hou, J.; Gao, L.Y.; Li, Q.; Li, Y.; Zhu, H.Y.; Yang, L.M.; Hu, J.B. Construction of a High-Density Genetic Map for Melon Using ddRAD-Seq Technology from a Population Derived from Flexuosus and Reticulatus Botanical Groups. Sci. Hortic. 2020, 272, 10953. [Google Scholar] [CrossRef]
  19. Kobayashi, H.; Haino, Y.; Iwasaki, T.; Tezuka, A.; Nagano, A.J.; Shimada, S. ddRAD-Seq Based Phylogeographic Study of Sargassum Thunbergii (Phaeophyceae, Heterokonta) around Japanese Coast. Mar. Environ. Res. 2018, 140, 104–113. [Google Scholar] [CrossRef] [PubMed]
  20. Kaviriri, D.K.; Zhang, Q.; Zhang, X.; Jiang, L.; Zhang, J.; Wang, J.; Khasa, D.P.; You, X.; Zhao, X. Phenotypic Variability and Genetic Diversity in a Pinus Koraiensis Clonal Trial in Northeastern China. Genes 2020, 11, 673. [Google Scholar] [CrossRef] [PubMed]
  21. Liu, J.X.; Zhou, M.Y.; Yang, G.Q.; Zhang, Y.X.; Ma, P.F.; Guo, C.; Vorontsova, M.S.; Li, D.Z. ddRAD Analyses Reveal a Credible Phylogenetic Relationship of the Four Main Genera of Bambusa-Dendrocalamus-Gigantochloa Complex (Poaceae: Bambusoideae). Mol. Phylogenet. Evol. 2020, 146, 106758. [Google Scholar] [CrossRef] [PubMed]
  22. Yang, G.Q.; Chen, Y.M.; Wang, J.P.; Guo, C.; Zhao, L.; Wang, X.-Y.; Guo, Y.; Li, L.; Li, D.-Z.; Guo, Z.-H. Development of a Universal and Simplified ddRAD Library Preparation Approach for SNP Discovery and Genotyping in Angiosperm Plants. Plant Methods 2016, 12, 39. [Google Scholar] [CrossRef] [PubMed]
  23. García, V.; Castro, P.; Die, J.V.; Millán, T.; Gil, J.; Moreno, R. QTL Analysis of Morpho-Agronomic Traits in Garden Asparagus (Asparagus Officinalis L.). Horticulturae 2022, 9, 41. [Google Scholar] [CrossRef]
  24. Feng, J.Y.; Zhao, S.; Li, M.; Zhang, C.; Qu, H.J.; Li, Q.; Li, J.W.; Lin, Y.; Pu, Z.G. Genome-Wide Genetic Diversity Detection and Population Structure Analysis in Sweetpotato (Ipomoea Batatas) Using RAD-Seq. Genomics 2020, 112, 1978–1987. [Google Scholar] [CrossRef] [PubMed]
  25. Shi, X.D.; Yin, Q.; Sang, Z.Y.; Zhu, Z.L.; Jia, Z.K.; Ma, L.Y. Prediction of Potentially Suitable Areas for the Introduction of Magnolia Wufengensis under Climate Change. Ecol. Indic. 2021, 127, 107762. [Google Scholar] [CrossRef]
  26. Elith, J.; Leathwick, J.R. Species Distribution Models: Ecological Explanation and Prediction Across Space and Time. Annu. Rev. Ecol. Evol. Syst. 2009, 40, 677–697. [Google Scholar] [CrossRef]
  27. Liu, Y.; Tian, B.; Ou, G.L. Displacement Distribution and Climate Explanation on Humid and Semi-Humid Ever green Broadleaved Forests Using Niche Model of Cyclobalanopsis Glauca and C. Glaucoides in China. Guihaia 2022, 42, 460–469. [Google Scholar] [CrossRef]
  28. Merow, C.; Smith, M.J.; Silander, J.A. A Practical Guide to MaxEnt for Modeling Species’ Distributions: What It Does, and Why Inputs and Settings Matter. Ecography 2013, 36, 1058–1069. [Google Scholar] [CrossRef]
  29. Hernandez, P.A.; Graham, C.H.; Master, L.L.; Albert, D.L. The Effect of Sample Size and Species Characteristics on Performance of Different Species Distribution Modeling Methods. Ecography 2006, 29, 773–785. [Google Scholar] [CrossRef]
  30. Yu, Y.L.; Wang, H.C.; Yu, Z.X.; Schinnerl, J.; Tang, R.; Geng, Y.P.; Chen, G. Genetic Diversity and Structure of the Endemic and Endangered Species Aristolochia Delavayi Growing along the Jinsha River. Plant Divers. 2021, 43, 225–233. [Google Scholar] [CrossRef] [PubMed]
  31. Yang, F.M.; Cai, L.; Dao, Z.L.; Sun, W.B. Genomic Data Reveals Population Genetic and Demographic History of Magnolia Fistulosa (Magnoliaceae), a Plant Species With Extremely Small Populations in Yunnan Province, China. Front. Plant Sci. 2022, 13, 811312. [Google Scholar] [CrossRef] [PubMed]
  32. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. Fastp: An Ultra-Fast All-in-One FASTQ Preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef] [PubMed]
  33. Catchen, J.; Bassham, S.; Wilson, T.; Currey, M.; O’Brien, C.; Yeates, Q.; Cresko, W.A. The Population Structure and Recent Colonization History of Oregon Threespine Stickleback Determined Using RAD-Seq. Mol. Ecol. 2013, 22, 2864–2883. [Google Scholar] [CrossRef] [PubMed]
  34. Chang, C.C.; Chow, C.C.; Tellier, L.C.; Vattikuti, S.; Purcell, S.M.; Lee, J.J. Second-Generation PLINK: Rising to the Challenge of Larger and Richer Datasets. GigaScience 2015, 4, s13742-015-0047-8. [Google Scholar] [CrossRef]
  35. Yang, J.; Lee, S.H.; Goddard, M.E.; Visscher, P.M. GCTA: A Tool for Genome-Wide Complex Trait Analysis. Am. J. Hum. Genet. 2011, 88, 76–82. [Google Scholar] [CrossRef] [PubMed]
  36. Pfaff, C.L.; Parra, E.J.; Bonilla, C.; Hiester, K.; McKeigue, P.M.; Kamboh, M.I.; Hutchinson, R.G.; Ferrell, R.E.; Boerwinkle, E.; Shriver, M.D. Population Structure in Admixed Populations: Effect of Admixture Dynamics on the Pattern of Linkage Disequilibrium. Am. J. Hum. Genet. 2001, 68, 198–207. [Google Scholar] [CrossRef] [PubMed]
  37. Rao, V.R.; Hodgkin, T. Genetic Diversity and Conservation and Utilization of Plant Genetic Resources. Plant Cell Tissue Organ. Cult. 2002, 68, 1–19. [Google Scholar] [CrossRef]
  38. Gao, R.; Liu, L.; Zhao, L.; Cui, S. Potentially Suitable Geographical Area for Monochamus Alternatus under Current and Future Climatic Scenarios Based on Optimized MaxEnt Model. Insects 2023, 14, 182. [Google Scholar] [CrossRef] [PubMed]
  39. Chen, Y.N.; Liu, Z.G.; Yu, T.; Liao, H.; Zhou, J.Y. Prediction of Potential Distribution of Prunus Mume Based on MaxEnt Model. Chin. Wild Plant Resour. 2024, 43, 107–113+126. [Google Scholar]
  40. Wu, T.; Lu, Y.; Fang, Y.; Xin, X.; Li, L.; Li, W.; Jie, W.; Zhang, J.; Liu, Y.; Zhang, L.; et al. The Beijing Climate Center Climate System Model (BCC-CSM): The Main Progress from CMIP5 to CMIP6. Geosci. Model Dev. 2019, 12, 1573–1600. [Google Scholar] [CrossRef]
  41. Zhang, R.L.; Liu, W.Y.; Zhang, Y.X.; Tu, D.D.; Zhu, L.Y.; Zhang, W.J. Genetic Diversity Analysis of Melocalamus Arrectus Based on Reduced-Representation Genome Sequencing. Mol. Plant Breed. 2021. [Google Scholar]
  42. Willoughby, J.R.; Sundaram, M.; Wijayawardena, B.K.; Kimble, S.J.A.; Ji, Y.; Fernandez, N.B.; Antonides, J.D.; Lamb, M.C.; Marra, N.J.; DeWoody, J.A. The Reduction of Genetic Diversity in Threatened Vertebrates and New Recommendations Regarding IUCN Conservation Rankings. Biol. Conserv. 2015, 191, 495–503. [Google Scholar] [CrossRef]
  43. Jorde, L.B. Linkage Disequilibrium and the Search for Complex Disease Genes. Genome Res. 2000, 10, 1435–1444. [Google Scholar] [CrossRef] [PubMed]
  44. Tikendra, L.; Potshangbam, A.M.; Amom, T.; Dey, A.; Nongdam, P. Understanding the Genetic Diversity and Population Structure of Dendrobium Chrysotoxum Lindl.-An Endangered Medicinal Orchid and Implication for Its Conservation. S. Afr. J. Bot. 2021, 138, 364–376. [Google Scholar] [CrossRef]
  45. Miller, M.R.; Dunham, J.P.; Amores, A.; Cresko, W.A.; Johnson, E.A. Rapid and Cost-Effective Polymorphism Identification and Genotyping Using Restriction Site Associated DNA (RAD) Markers. Genome Res. 2007, 17, 240–248. [Google Scholar] [CrossRef] [PubMed]
  46. Berihulay, H.; Li, Y.; Liu, X.; Gebreselassie, G.; Islam, R.; Liu, W.; Jiang, L.; Ma, Y. Genetic Diversity and Population Structure in Multiple Chinese Goat Populations Using a SNP Panel. Anim. Genet. 2019, 50, 242–249. [Google Scholar] [CrossRef] [PubMed]
  47. Marrano, A.; Birolo, G.; Prazzoli, M.L.; Lorenzi, S.; Valle, G.; Grando, M.S. SNP-Discovery by RAD-Sequencing in a Germplasm Collection of Wild and Cultivated Grapevines (V. Vinifera L.). PLoS ONE 2017, 12, e0170655. [Google Scholar] [CrossRef] [PubMed]
  48. Dubreuil, M.; Riba, M.; Mayol, M. Genetic Structure and Diversity in Ramonda Myconi (Gesneriaceae): Effects of Historical Climate Change on a Preglacial Relict Species. Am. J. Bot. 2008, 95, 577–587. [Google Scholar] [CrossRef] [PubMed]
  49. Fu, Q.; Lu, G.H.; Fu, Y.H.; Wang, Y.Q. Genetic Differentiation between Two Varieties of Oreocharis Benthamii (Gesneriaceae) in Sympatric and Allopatric Regions. Ecol. Evol. 2020, 10, 7792–7805. [Google Scholar] [CrossRef] [PubMed]
  50. Ni, X.W.; Huang, Y.L.; Wu, L.; Zhou, R.C.; Deng, S.; Wu, D.R.; Wang, B.S.; Su, G.H.; Tang, T.; Shi, S.H. Genetic Diversity of the Endangered Chinese Endemic Herb Primulina Tabacum (Gesneriaceae) Revealed by Amplified Fragment Length Polymorphism (AFLP). Genetica 2006, 127, 177–183. [Google Scholar] [CrossRef] [PubMed]
  51. Zhu, X.; Zou, R.; Tang, J.; Deng, L.; Wei, X. Genetic Diversity Variation during the Natural Regeneration of Vatica Guangxiensis, an Endangered Tree Species with Extremely Small Populations. Glob. Ecol. Conserv. 2023, 42, e02400. [Google Scholar] [CrossRef]
  52. Xu, G.; Liang, Y.; Jiang, Y.; Liu, X.S.; Hu, S.L.; Xiao, Y.F.; Hao, B.B. Genetic diversity and population structure of Bretschneidera sinensis, an endangered species: Genetic diversity and population structure of Bretschneidera sinensis, an endangered species. Biodivers. Sci. 2014, 21, 723–731. [Google Scholar]
  53. Mao, C.L.; Zhang, F.L.; Li, X.Q.; Yang, T.; Zhao, Q.; Hu, Y.H.; Wu, Y. Genetic Diversity of Horsfieldia Pandurifolia Based on AFLP Markers. J. Trop. Subtrop. Bot. 2020, 28, 271–276. [Google Scholar] [CrossRef]
  54. Pan, Y.Z.; Wang, X.Q.; Sun, G.L.; Li, F.S.; Gong, X. Application of RAD Sequencing for Evaluating the Genetic Diversity of Domesticated Panax Notoginseng (Araliaceae). PLoS ONE 2016, 11, e0166419. [Google Scholar] [CrossRef] [PubMed]
  55. Zhu, G.P.; Liu, G.Q.; Bu, W.J.; Gao, Y. Ecological niche modeling and its applications in biodiversity conservation: Ecological niche modeling and its applications in biodiversity conservation. Biodivers. Sci. 2013, 21, 90–98. [Google Scholar]
  56. Garzón, M.B.; Blazek, R.; Neteler, M.; Dios, R.S.D.; Ollero, H.S.; Furlanello, C. Predicting Habitat Suitability with Machine Learning Models: The Potential Area of Pinus Sylvestris L. in the Iberian Peninsula. Ecol. Model. 2006, 197, 383–393. [Google Scholar] [CrossRef]
  57. Manthey, M.; Box, E.O. Realized Climatic Niches of Deciduous Trees: Comparing Western Eurasia and Eastern North America. J. Biogeogr. 2007, 34, 1028–1040. [Google Scholar] [CrossRef]
  58. Kramer-Schadt, S.; Niedballa, J.; Pilgrim, J.D.; Schröder, B.; Lindenborn, J.; Reinfelder, V.; Stillfried, M.; Heckmann, I.; Scharf, A.K.; Augeri, D.M.; et al. The Importance of Correcting for Sampling Bias in MaxEnt Species Distribution Models. Divers. Distrib. 2013, 19, 1366–1379. [Google Scholar] [CrossRef]
  59. Coelho, N.; Gonçalves, S.; Romano, A. Endemic Plant Species Conservation: Biotechnological Approaches. Plants 2020, 9, 345. [Google Scholar] [CrossRef] [PubMed]
  60. Desai, P.; Gajera, B.; Mankad, M.; Shah, S.; Patel, A.; Patil, G.; Narayanan, S.; Kumar, N. Comparative Assessment of Genetic Diversity among Indian Bamboo Genotypes Using RAPD and ISSR Markers. Mol. Biol. Rep. 2015, 42, 1265–1273. [Google Scholar] [CrossRef] [PubMed]
  61. Gutiérrez-Lozano, M.; Sánchez-González, A.; Octavio-Aguilar, P.; Galván-Hernández, D.M.; Vázquez-García, J.A. Genetic Diversity and Structure of Magnolia Mexicana (Magnoliacea): A Threatened Species in Eastern Mexico. Silvae Genet. 2023, 72, 132–142. [Google Scholar] [CrossRef]
Plants 13 02093 g001
Figure 1. Phylogenetic tree for Henckelia longisepala samples based on SNP data, with 14 samples from the HLS population in blue; 4 samples from the NXC population in red; and 14 samples from the TMC population in green.
Figure 1. Phylogenetic tree for Henckelia longisepala samples based on SNP data, with 14 samples from the HLS population in blue; 4 samples from the NXC population in red; and 14 samples from the TMC population in green.
Plants 13 02093 g001
Plants 13 02093 g002
Figure 2. Principal components analysis (PCA) of all 32 Henckelia longisepala samples, with the proportion of the variance explained being 27% for PC1 and 16.7% for PC2.
Figure 2. Principal components analysis (PCA) of all 32 Henckelia longisepala samples, with the proportion of the variance explained being 27% for PC1 and 16.7% for PC2.
Plants 13 02093 g002
Plants 13 02093 g003
Figure 3. Distribution of CV error values corresponding to different K values; the CV error was minimized for K = 3, indicating that the optimal number of clusters for the study samples is three.
Figure 3. Distribution of CV error values corresponding to different K values; the CV error was minimized for K = 3, indicating that the optimal number of clusters for the study samples is three.
Plants 13 02093 g003
Plants 13 02093 g004
Figure 4. Genetic structure of Henckelia longisepala based on analysis; the results are K = 2 and K = 3.
Figure 4. Genetic structure of Henckelia longisepala based on analysis; the results are K = 2 and K = 3.
Plants 13 02093 g004
Plants 13 02093 g005
Figure 5. Potential geographic distribution of Henckelia longisepala in the current period.
Figure 5. Potential geographic distribution of Henckelia longisepala in the current period.
Plants 13 02093 g005
Plants 13 02093 g006
Figure 6. Potential geographic distribution of Henckelia longisepala for different time periods under two future climate change scenarios (SSP1-2.6, SSP5-8.5). Forecast maps are shown for the SSP1-2.6 climate change scenario in the 2030s (a) and 2050s (b) and for the SSP5-8.5 climate change scenario in the 2030s (c) and 2050s (d).
Figure 6. Potential geographic distribution of Henckelia longisepala for different time periods under two future climate change scenarios (SSP1-2.6, SSP5-8.5). Forecast maps are shown for the SSP1-2.6 climate change scenario in the 2030s (a) and 2050s (b) and for the SSP5-8.5 climate change scenario in the 2030s (c) and 2050s (d).
Plants 13 02093 g006
Table 1. Collection site details for Henckelia longisepala.
Table 1. Collection site details for Henckelia longisepala.
PopLongitude (° E)Latitude (° N)Altitude (m)LocalityNumber of Samples
HLS102.349422.64571009Yunnan Lüchun14
NXC103.889722.6095462Yunnan Hekou4
TMC103.026422.6250947Yunnan Jinping14
Table 2. Overview of the 19 climatic variables.
Table 2. Overview of the 19 climatic variables.
Climate VariableDescriptionClimate VariableDescription
bio1Annual mean temperaturebio11Mean temperature of coldest quarter
bio2Mean diurnal rangebio12Annual precipitation
bio3Isothermalitybio13Precipitation of wettest month
bio4Temperature seasonalitybio14Precipitation of driest month
bio5Max temperature of warmest monthbio15Precipitation seasonality
bio6Min temperature of coldest monthbio16Precipitation of wettest quarter
bio7Temperature annual rangebio17Precipitation of driest quarter
bio8Mean temperature of wettest quarterbio18Precipitation of warmest quarter
bio9Mean temperature of driest quarterbio19Precipitation of coldest quarter
bio10Mean temperature of warmest quarter
Table 3. ddRAD sequencing data of Henckelia longisepala.
Table 3. ddRAD sequencing data of Henckelia longisepala.
Sample_NameHQ ReadsTotal_Base (bp)HQ Bases (bp)HQ Bases %GC (%)Q30 (%)
HLS_A112,508,536 1,847,065,824 1,781,995,529 96.48 40.97 94.94
HLS_A214,332,744 2,120,032,512 2,041,843,119 96.31 43.12 94.92
HLS_A336,818,088 5,410,729,440 5,257,449,583 97.17 40.22 96.41
HLS_A434,107,414 5,009,117,184 4,872,148,968 97.27 39.26 96.19
HLS_A530,046,774 4,413,880,800 4,292,836,294 97.26 38.49 96.10
HLS_B132,220,506 4,733,449,056 4,602,590,836 97.24 39.03 96.29
HLS_B232,318,314 4,827,110,569 4,675,400,724 96.86 38.35 96.19
HLS_B332,792,914 4,906,151,916 4,745,839,024 96.73 38.04 96.00
HLS_B414,576,362 2,136,863,520 2,087,739,170 97.70 38.95 97.89
HLS_B515,292,504 2,244,107,808 2,189,839,488 97.58 39.19 97.82
HLS_C115,307,248 2,253,815,136 2,191,488,136 97.23 42.68 97.67
HLS_C218,084,322 2,656,833,696 2,589,175,360 97.45 42.18 97.71
HLS_D116,292,676 2,398,851,360 2,331,478,275 97.19 41.94 97.71
HLS_D216,406,160 2,414,157,408 2,347,565,196 97.24 43.57 97.66
NX_115,163,360 2,221,348,320 2,170,983,575 97.73 38.85 97.84
NX_217,780,208 2,614,287,456 2,545,830,403 97.38 39.43 97.73
NX_315,545,240 2,282,696,928 2,225,675,817 97.50 39.43 97.73
NX_414,758,408 2,167,298,496 2,112,253,626 97.46 38.15 97.72
TM_A114,494,752 2,130,803,712 2,076,078,685 97.43 43.20 97.82
TM_A216,837,792 2,481,160,320 2,410,517,148 97.15 43.16 97.59
TM_B114,951,876 2,195,641,152 2,140,724,678 97.50 39.91 97.74
TM_B214,226,138 2,093,088,960 2,037,690,222 97.35 42.23 97.68
TM_B314,467,912 2,133,347,904 2,071,285,801 97.09 41.71 97.65
TM_B414,766,900 2,175,391,872 2,115,346,007 97.24 41.65 97.70
TM_B513,956,478 2,053,735,200 1,999,096,791 97.34 40.38 97.66
TM_B615,690,434 2,297,819,808 2,247,479,371 97.81 39.11 97.97
TM_C115,504,368 2,272,298,112 2,218,737,176 97.64 38.23 97.86
TM_C217,817,304 2,612,644,992 2,550,693,042 97.63 39.69 97.88
TM_C316,814,138 2,466,775,008 2,407,381,170 97.59 40.25 97.84
TM_C417,449,764 2,558,971,008 2,495,422,887 97.52 41.10 97.80
TM_C518,474,068 2,713,303,584 2,644,282,676 97.46 39.40 97.75
TM_C617,334,360 2,539,707,264 2,482,303,852 97.74 38.78 97.87
Table 4. SNP information for Henckelia longisepala samples from ddRAD sequencing.
Table 4. SNP information for Henckelia longisepala samples from ddRAD sequencing.
SampleSNP NumberTransitionsTransversionsNumber of Heterozygous SNPsHeterozygosity (%)Number of Homozygous SNPsHomozygosity (%)Ts/Tv
HLS_A137,02821,25115,77719,18551.8117,843 48.191.35
HLS_A237,91521,85716,05821,56956.8916,346 43.111.36
HLS_A350,95729,20221,75529,49957.8921,458 42.111.34
HLS_A449,24628,24521,00128,97458.8420,272 41.161.34
HLS_A546,99226,99619,99626,90357.2520,089 42.751.35
HLS_B148,14227,60320,53927,44257.0020,700 43.001.34
HLS_B252,10429,69622,40829,23556.1122,869 43.891.33
HLS_B342,69924,25018,44923,37254.7419,327 45.261.31
HLS_B448,88427,79921,08527,86757.0121,017 42.991.32
HLS_B548,64427,82120,82323,86249.0524,782 50.951.34
HLS_C144,79625,95418,84224,13753.8820,659 46.121.38
HLS_C250,28128,95221,32929,74559.1620,536 40.841.36
HLS_D149,45328,73020,72331,82864.3617,625 35.641.39
HLS_D243,26325,06918,19420,52947.4522,734 52.551.38
NX_178,99646,84232,15427,57134.9051,425 65.101.46
NX_265,40339,06426,33917,77027.1747,633 72.831.48
NX_383,26849,29833,97029,08034.9254,188 65.081.45
NX_479,20846,79432,41427,65834.9251,550 65.081.44
TM_A140,82723,64717,18022,51255.1418,315 44.861.38
TM_A242,95324,84318,11024,40756.8218,546 43.181.37
TM_B137,41721,61315,80422,60360.4114,814 39.591.37
TM_B240,85423,63717,21722,33254.6618,522 45.341.37
TM_B340,33323,21117,12223,88359.2116,450 40.791.36
TM_B443,11724,86318,25425,11458.2518,003 41.751.36
TM_B538,97422,39616,57821,08654.1017,888 45.901.35
TM_B647,32727,38819,93927,72458.5819,603 41.421.37
TM_C140,34123,22617,11518,16645.0322,17554.971.36
TM_C248,28627,92320,36329,16560.4019,12139.601.37
TM_C348,60127,85920,74230,90263.5817,69936.421.34
TM_C443,14224,79618,34624,09655.8519,04644.151.35
TM_C548,57727,99120,58629,01259.7219,56540.281.36
TM_C695,53555,92739,60835,01736.6560,51863.351.41
Table 5. Genetic diversity statistics for each Henckelia longisepala population.
Table 5. Genetic diversity statistics for each Henckelia longisepala population.
Pop Obs HetObs HomExp HetExp HomPiFIS
HLS0.12560.87440.18830.81170.19900.2079
TMC0.11820.88180.15720.84280.17950.1868
NXC0.12100.87900.11510.88490.14080.0420
Table 6. Inter-population FST statistics.
Table 6. Inter-population FST statistics.
PopHLSTMCNXC
HLS 0.19240.3680
TMC 0.4071
Table 7. Genetic distance between populations of Henckelia longisepala.
Table 7. Genetic distance between populations of Henckelia longisepala.
PopHLSTMC
TMC0.4524
NXC0.82760.7672
Table 8. Contributions of the bioclimatic variables to the MaxEnt model simulations.
Table 8. Contributions of the bioclimatic variables to the MaxEnt model simulations.
Environment VariablePercent ContributionPermutation Importance
bio1836.743.9
bio725.433.9
bio1617.516.3
bio911.92.5
bio173.22.3
bio123.10.6
bio32.20.5
Table 9. Amount of suitable Henckelia longisepala habitat predicted for different time periods (/104 km2).
Table 9. Amount of suitable Henckelia longisepala habitat predicted for different time periods (/104 km2).
ScenarioPeriodHigh SuitabilityRatio %Medium SuitabilityRatio
%		Low SuitabilityRatio
%		TotalRatio
%
Current —5.50 —11.25 —15.67 —32.42 —
ssp1-2.62030s5.24 −4.73 7.26 −35.47 14.12 −9.89 26.63 −17.86
2050s5.10 −7.27 7.76 −31.02 14.93 −4.72 27.79 −14.28
ssp5-8.52030s4.03 −26.73 4.99 −55.64 10.76 −31.33 19.78 −38.99
2050s4.98 −9.45 7.66 −31.91 12.63 −19.40 25.27 −22.05
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhao, R.; Huang, N.; Zhang, Z.; Luo, W.; Xiang, J.; Xu, Y.; Wang, Y. Genetic Diversity Analysis and Prediction of Potential Suitable Areas for the Rare and Endangered Wild Plant Henckelia longisepala. Plants 2024, 13, 2093. https://doi.org/10.3390/plants13152093

AMA Style

Zhao R, Huang N, Zhang Z, Luo W, Xiang J, Xu Y, Wang Y. Genetic Diversity Analysis and Prediction of Potential Suitable Areas for the Rare and Endangered Wild Plant Henckelia longisepala. Plants. 2024; 13(15):2093. https://doi.org/10.3390/plants13152093

Chicago/Turabian Style

Zhao, Renfen, Nian Huang, Zhiyan Zhang, Wei Luo, Jianying Xiang, Yuanjie Xu, and Yizhi Wang. 2024. "Genetic Diversity Analysis and Prediction of Potential Suitable Areas for the Rare and Endangered Wild Plant Henckelia longisepala" Plants 13, no. 15: 2093. https://doi.org/10.3390/plants13152093

APA Style

Zhao, R., Huang, N., Zhang, Z., Luo, W., Xiang, J., Xu, Y., & Wang, Y. (2024). Genetic Diversity Analysis and Prediction of Potential Suitable Areas for the Rare and Endangered Wild Plant Henckelia longisepala. Plants, 13(15), 2093. https://doi.org/10.3390/plants13152093

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop