- Research
- Open access
- Published:
Identification and characterization of a new species of Taxus — Taxus qinlingensis by multiple taxonomic methods
BMC Plant Biology volume 24, Article number: 658 (2024)
Abstract
Background
The taxonomy of Taxus Linn. remains controversial due to its continuous phenotypic variation and unstable topology, thus adversely affecting the formulation of scientific conservation strategies for this genus. Recently, a new ecotype, known as Qinling type, is mainly distributed in the Qinling Mountains and belongs to a monophyletic group. Here, we employed multiple methods including leaf phenotype comparison (leaf shapes and microstructure), DNA barcoding identification (ITS + trnL-trnF + rbcL), and niche analysis to ascertain the taxonomic status of the Qinling type.
Results
Multiple comparisons revealed significant differences in the morphological characters (length, width, and length/width ratio) among the Qinling type and other Taxus species. Leaf anatomical analysis indicated that only the Qinling type and T. cuspidata had no papilla under the midvein or tannins in the epicuticle. Phylogenetic analysis of Taxus indicated that the Qinling type belonged to a monophyletic group. Moreover, the Qinling type had formed a relatively independent niche, it was mainly distributed around the Qinling Mountains, Ta-pa Mountains, and Taihang Mountains, situated at an elevation below 1500 m.
Conclusions
Four characters, namely leaf curvature, margin taper, papillation on midvein, and edges were put forward as primary indexes for distinguishing Taxus species. The ecotype Qingling type represented an independent evolutionary lineage and formed a unique ecological niche. Therefore, we suggested that the Qingling type should be treated as a novel species and named it Taxus qinlingensis Y. F. Wen & X. T. Wu, sp. nov.
Background
Species are the basic units of biodiversity, yet distinguishing closely related species can be challenging [1]. The fundamental goals of species delimitation are not only to identify species, but also to clarify the phylogenetic relationships among and within species [2]. In addition, species delimitation also plays a vital role in the formulation of conservation policies [3,4,5,6]. Incorrect species delimitation can lead to erroneous or misleading conclusions about the evolutionary history and adaptive mechanisms of the species, consequently resulting in inappropriate conservation management decisions and adverse effects on management strategies [7, 8].
An independently evolved species (namely, the so-called good species) can be effectively delineated from different perspectives based on the classical taxonomy. However, most genealogies are always in divergence state in nature, making it challenging to discriminate species solely from one perspective [9]. More than 30 species concepts have been proposed [10], with the morphological species concept is the most commonly used [9,10,11]. With the development of DNA sequencing, the phylogenetic species concept has become indispensable in diagnosing species [9]. In addition, difference in the micro-structure of vegetative organs can be used as a basis for plant taxonomy [12, 13]. Furthermore, speciation is usually accompanied by the niche differentiation, and thus ecological concepts are also introduced to discriminate species [11, 14,15,16,17,18]. Integrating these multiple concepts enhances the rationality of species delimitation and improves the accuracy of biodiversity research results [19]. Now “the integrative species concept” is widely used in species classification [1, 9, 20].
Taxus is the most widespread genus with the largest number of species in Taxaceae [21, 22]. However, the classification of Taxus remains controversial due to its high phenotypic plasticity, niche overlapping across different lineages, and limited morphological diagnostic characters, particularly within Asian taxa [23,24,25,26]. There are three species (Taxus wallichiana, Taxus fuana and Taxus cuspidata) and three varieties (Taxus wallichiana var. wallichiana, Taxus wallichiana var. chinensis, Taxus wallichiana var. mairei) among the Taxus according to Flora of China (FOC) [27]. However, Farjon [22] categorized Taxus into five species in China (T. wallichiana, T. contorta (synonym of T. fuana), T. cuspidata, Taxus chinensis, and Taxus mairei). Möller et al. [26] divided Taxus into eight species (T. cuspidata, T. chinensis, T. mairei, T. wallichiana, T. contorta, Taxus calcicola, Taxus phytonii, Taxus florinii) and three lineages (Qinling, Emei, and Huangshan type) in Asia based on the morphological, cpDNA, and combined nuclear evidence (ITS and NEEDLY). Among these, the Qinling type mainly distributed from the northern parts of the Qinling Mountains to the southern parts of Taihang Mountains. Emei type is mostly distributed around Mount Emei and Qiaojia in the northeast of Yunnan province [28]. Huangshan type is mainly located in the east of China such as Huangshan, Sanqingshan, and Jiuhua Mountains [26]. Thus, Qinling type has a wide distribution range and obvious geographical boundary, compared with the other two lineages. Qinling type was thought to be closely related to Taxus mairei based on the trnL-trnF sequences in the phylogeographic study of T. wallichiana (Taxaceae) [23]. However, molecular data (rbcL, matK, trnH - psbA, trnL-F) revealed Qinling type forming a distinct lineage and being treated as an ecotype [24, 29]. Complete plastomes analysis suggested the Qinling type had formed a monophyletic clade and it was closely related to T. contorta [25]. Möller et al. [26] conjectured that Qinling type was the hybrid of T. contorta (♀) and T. Huangshan type (♂) based on 13 cpDNA markers and two nuclear regions, suggesting that its morphological variation was not sufficient for species classification. Moreover, prior research has indicated phenotypic variation within the Qinling type, warranting further classification of its taxonomic status.
In this study, we employed “the integrative species concept” to determine the taxonomic status of the Qinling type based on a comprehensive analysis integrating morphological observations, microstructure of leaves, phylogenetic analysis, and niche analysis. Moreover, we established a retrieval list of Taxus taxonomy. Our findings aim to enhance the understanding of Taxus taxonomy, facilitating the formulation of effective conservation strategies for the species and improving our awareness of biodiversity.
Results
Comparison of morphological characters and microstructure among different Taxus species
Leaf morphology comparison of six Taxus species was shown in Fig. 1. It shown the differences of leaf curvature, such as falcate (T. wallichiana), sword-shaped (T. contorta), and linear (T. chinensis). Also, the leaf margin taper were different, particularly in tip and curvature. Leaf scanning electron microscopy (SEM) shown that stomatal apparatus and cuticular papillae are only distributed on the lower epidermis, with two stomatal bands on either side of the midvein (Fig. 2B-C). Leaf anatomy analysis indicated that Taxus leaves were a typical bifacial leaf, presenting the characteristic of shady plants (Fig. 3). The leaf blade is composed of three parts, epidermis, mesophyll, and veins. The epidermal is mono-layer and covered with cuticle (Cu). The mesophyll comprise palisade tissue (Pa) and spongy tissue (Sp).Within the palisade tissue, cells are irregularly shaped and arranged in single-layer, double-layer, or triple-layer, exhiting a clear hierarchal structure (Table 1). Meanwhile, spongy tissue cells are irregularly shaped, differing in cell sizes, with large gaps, and a loose arrangement (Fig. 3). Leaf veins, located centrally within the leaf, comprise xylem (Xy) and phloem (Ph). Notably, among the six Taxus species, only the Qinling type and T. cuspidata lack cuticular papillae under the midvein and tannins in the epicuticle (Fig. 3). However, the hard texture of leaves, the high rise of midvein on the adaxial side, and the involution of leaf edges were typical characters to distinguish T. cuspidata from Qinling type. Thus, Taxus species can be discriminated based on the characters of leaf curvature, margin taper, papillation on midvein, and edges.
Phenotype of the six Taxus species (Left: the front of leaves; Right: the back of leaves) (Qinling type from Lushi, Henan; T. wallichiana from Kunming, Yunnan; T. contorta from Jilong, Xizang; T. chinensis from Jiangyou, Sichuan; T. mairei from Changsha, Hunan; T. cuspidata from Harbin, Heilongjiang)
Leaf anatomy of the six Taxus species. Notes Xy: Xylem; Ph: Phloem; Pa: Palisade tissue; Sp: Spongy tissue; Cu: Cuticle; Ue: Epicuticle; Le: Lower epidermis; Pi: Papillate. Red rectangle shown whether there were cuticular papillae under the midvein; Red circle shown the rise degree of midvein on the adaxial side. The black substance shown in the epicuticle means the leaves contain tannins (Species of T. wallichiana, T. contota, T. chinensis and T. mairei)
For leaf morphology, we compared the leaf length, width, and length/width ratio (Table S1) among the six Taxus species based on non-parametric statistics. Significant differences were observed in the these morphological characters in most pairwise species comparison with the following exceptions (Fig. 4, Table S1). As for the leaf length (Fig. 4A), no significant difference was found between Qinling type and T. contorta, between Qinling type and T. mairei, between T. contorta and T. mairei. Leaf width multiple comparison revealed no significant difference between T. wallichiana and T. contorta, between T. chinensis and T. cuspidata (Fig. 4B). No significant difference was observed in length/width ratio between Qinling type and T. chinensis, between T. wallichiana and T. cuspidata (Fig. 4C).
Multiple leaf comparison of the six Taxus species based on non-parametric test. (A) Leaf length, Kruskal-Wallis, χ2(5) = 308.86; (B) Leaf width, Kruskal-Wallis, χ2(5) = 532.62; (C) Leaf length/width ratio. Kruskal-Wallis, χ2(5) = 480.91. Comparisons marked “ns” did not reach significance, the other all reached significance
Species delimitation
The length of individual sequence for ITS, trnL-trnF, and rbcL, as well as the concatenated sequences (the sum of the length of the three sequences), were 1104, 947, 746, and 2797 bp, respectively. The incongruence length difference test (ILD) indicated there was no significant difference (0.05 < p = 1.00) in the concatenated sequence length among species. Thus, the taxonomic discrimination of Taxus species was conducted based on concatenated sequence. Substitution saturation test (Iss < Iss.c and p < 0.05) showed that the combination of the three sequences was suitable for phylogenetic analysis. Bayesian estimation (BI) and maximum likelihood estimation (ML) yielded congruent topologies (Fig. 5 and 6). 18 populations from the Qinling Mountains formed a monophyletic group and were identified as the Qinling type, while the other Taxus species also formed their respective monophyletic groups. Consequently, these results suggest that six Taxus species had a stable systematic position.
Niche overlap analyses
Niche overlap analyses showed that T. chinensis had the largest niche breadth (0.310), followed by T. cuspidata (0.176), while the Qinling type (0.063) exhibited the smallest niche breadth (Fig. 7A-E). Schoener’s D metric revealed varying levels of niche overlap, ranging from very limited overlap (D = 0.052, between Qinling type and T. mairei) to moderate overlap (D = 0.440, between T. qinlingensis and T. chinensis) (Fig. 7F-J). The similarity score for ecological niche models (ENMs) constructed for the actual occurrences of the two species was lower than expected under the null hypothesis, indicating that the two species are more divergent than expected based on the available habitat for the Qinling type (Fig. 7F-J). However, pairwise comparisons did not yield significant difference (p > 0.05). Niche similarity tests showed no significant difference between the two species in pairwise comparison (p > 0.05), except between Qinling type and T. chinensis (p < 0.05) (Fig. S1 A-E).
Principal component analysis (PCA) of the driving factors (14 variables) revealed that PC1 and PC2 explained 38% and 18.5% of the total variance, respectively. Rooting conditions, elevation, and oxygen availability to roots were the dominant factors that led to the niche differentiation of the Taxus species in the current study (Fig. S1 F).
Niche quantification and overlap based on principal component analysis (PCA-env). Panel A-E: Each color (red or green) represents one species, purple represents the degree of overlap between the two species. Histograms (F) - (J) show the observed niche overlap D between the two ranges (bars with a diamond) and simulated niche overlaps (gray bars) on the tests of niche equivalency. D represents the degree of niche overlap, the significance of test is shown (Non-significant, p > 0.05; ***, p < 0.001)
Discussion
Unique morphological characters and anatomical structure of Qinling type
A previous study suggested using one bud scale and 26 leaf characters to discriminate five Taxus species [30], discovering that T. contorta showed the most significant morphological difference compared to other Taxus species. These 27 morphological characters have been used to identify new Taxus species and all other species and types in Eurasia [28, 30], and lay the foundation for identifying Taxus species [31]. However, considering the leaf morphological characters and anatomical structure, we proposed four key morphological characters to discriminate Taxus species: curvature, margin taper, papillation on midvein, and edges (Fig. 2A). These characters can effectively identify Taxus species and simplify the species identification procedure. Additionally, significant differences were observed in multiple comparisons of leaf phenotype (length, width, and length/width ratio). Therefore, these quantitative traits can serve as indictors for categorizing Taxus species.
Tannins can serve as the complementary indicator for phytochemistry classification [32]. Among the six Taxus species, only the Qinling type and T. cuspidata lacked tannins in the epidermal cells. Therefore, the presence or absence of tannins, along with the cuticular papillae under the midvein, are two important characters to distinguish the Qinling type from the T. chinensis and T. mairei in their coexisting region in the warm temperate and subtropical China. Furthermore, altitude difference can also be used as one of the criterias to discriminate Qinling type (low altitude) from T. chinensis and T. wallichiana (high altitude) in their sympatric distribution areas, as evidenced by field investigations and prior findings [22, 33]. In conclusion, the analysis of leaf phenotype and micro-structure analysis indicated that the Qinling type was significantly different from other Taxus species, exhibiting unique morphological and anatomical characteristics.
Stable phylogenetic position and distribution pattern of Qinling type
Qinling type was once thought to be closely associated with T. wallichiana var. mairei based on the RFLP markers [23]. Then, Liu et al. [24] proposed Qinling type as an independent lineage by the DNA barcoding. Fu et al. [25] pointed out that the Qinling type might be most closely related to T. contorta based on the complete plastid genomes. Möller et al. [26] proposed that the Qinling type is the hybrid of T. contorta (♀) × T. Huangshan type (♂) based on 13 cpDNA markers and two nuclear regions, and it was an independent lineage. In current study, Qinling type formed a monophyletic lineage based on one nuclear regions (ITS) and two cpDNA makers (trnL-trnF + rbcL), and its divergence time was prior to T. chinensis and T. mairei (Fig. 5 and 6). Taken together, these results confirmed that Qinling type is an independent lineage. However, the position of Qinling type in the phylogenetic tree is constantly changing. The inconsistent findings of Qinling type may be attributed to the conflicts between nuclear genomes and chloroplast genomes, the ancient hybridization, and the reciprocal chloroplast capture [34, 35]. From the phylogenetic perspective, the Qinling type represent a steady lineage. This may be attributed to geographic isolation and short migration distances, leading to reduced gene exchange between different species over time [36,37,38].
Previous studies have revealed that Rhinopithecus brelichi [39] and Juglans regia [40] all originated from ancient hybridization events, with their respective hybridization periods being 1.8 Mya and 3.45 Mya. In Taxus species, Qin et al. revealed the divergence between T. wallichiana and T. florinii occurred in 5.9 Mya (± 0.39) based on population-genomic analysis from dd-RAD data [41]. This time approximately aligns with the predication made by Möller et al. [26] (4.9–6.8 Mya). Furthermore, Möller et al. [26] also proposed that the hybridization event of Qinling type occurred in the late Miocene, around 6.8 Mya (95% HPD = 1.9–12.7 Mya) years ago. Hence, Qinling type belongs to the ancient hybridization event rather than a recent one.
Möller et al. [28] investigated the geographical distribution of the Qinling type and found it mainly distributed around the western Sichuan basin, spanning from the northern Qinling Mountains to the southern Taihang Mountains, indicating a relatively independent niche of the Qinling type in terms of geographical distribution and elevation. Furthermore, our previous field investigation showed that Qinling type predominantly inhabited the lower elevations of the Qinling Mountains and Taihang Mountains. When Qinling type was sympatric with T. chinensis and T. wallichiana, its altitude was lower than 1500 m [28, 42]. Nevertheless, our niche similarity tests demonstrated that Qinling type had niche overlap with T. chinensis, possibly because niche overlap analysis fails to capture the vertical spatial niche differences between the two species [16]. Notably, elevation was the main driving factor in the differentiation of Taxus species from the PCA-env analysis and species distribution models (SDMs) [33].
Suggestions for future research priorities on Taxus
DNA barcoding is a powerful taxonomic tool to identify and discover species [43,44,45]. Liu et al. [29] showed that trnL - trnF is an ideal barcode for Taxus, effectively distinguishing all Taxus species and attaching to ITS for hybird identification. The purpose of this study was to determine the taxonomic status of the Qinling type. Therefore, we used the combination of ITS + trnL-trnF + rbcL to identify Taxus species. Möller et al. [26] reconstructed the phylogeny of the extant Taxus lineages, yet only one sample was included in each taxon. In our study, a total of 66 individuals of Qinling type were collected from 18 populations in the Qinling mountains to confirm its taxonomic status. We found Qinling type showed obvious phenotypic and niche differenation compared with the other Taxus species, and it formed an independent lineage. Furthermore, we have completed Double Digest Restriction Associated DNA (dd-RAD) Sequencing for Qinling type and its related species in their sympatric distribution. Phylogenetic analysis, admixture analysis, and demographic analysis indicated that Qinling type was an independent lineage not originated from recent hybridization (unpublished data). In addition, Qinling type is also found in Sichuan and Yunnan provinces [29]. However, due to the sporadic distribution of these samples, they were not included in this study (morphological and phylogenetic analysis). Next, we will further to explore the origin and diffusion corridors of Taxus qinlingensis based on a more comprehensive sampling strategy.
According to the International Union for Conservation of Nature (IUCN), three Taxus species (T. wallichiana, T. contorta, and T. cuspidata) are classified as endangered, while two Taxus species are considered vulnerable (T. chinensis and T. mairei) (IUCN, 2022a) [46]. In China, all the Taxus species are under the first-grade state protection (National Forestry and Grassland Administration, 2021). In general, Taxus species are dioecious [22]. However, our previous investigation found that some individuals of T. qinlingensis around the Qinling Mountains were monoecious (unpublished data). Up to now, the sex transformation mechanisms of Taxus remain largely unknown [47, 48]. Generally speaking, the mechanisms of sex determination systems including genetic sex determination (GSD), environment sex determination (ESD), and the combination of GSD and ESD [49]. These mechanisms reveals that plants have to adopt different strategies to cope with diverse biotic and abiotic stresses for survival in adverse environmental conditions [50, 51]. In the future, with the application of whole-genome sequencing [52,53,54], the conservation strategies for T. qinlingensis can be explored from its habitat water scarcity, forest regeneration mechanism, seed after-ripening, and sex determination mechanism [55, 56].
Conclusion
This study employed a combination of leaf phenotype comparison, phylogenetic analysis, and niche analysis to reveal the classification status of the Qinling type, and we named it Taxus qinlingensis Y. F. Wen & X. T. Wu sp. nov (with the detailed information provided in Taxonomic Key). Soil conditions and elevation played a vital role in the speciation of Taxus. These results enhance the importance of using “the integrative species concept” in species categorization, and reasonable species assessments will help to raise awareness of the importance of biodiversity. Future studies should focus on the origin and diffusion path of Taxus qinlingensis, exploring its sex determination locus.
Methods
Sample collection and DNA extraction
Several field investigations of Qinling type were performed from the year 2019 to 2022, covering the northern part of Qinling Mountains to the southern part of Taihang Mountains. Habitat, seed branches, and cones of Qinling type were shown in Fig. 8. A total of 66 samples were collected from 18 Qinling type populations, with 2–6 samples per population. Detailed information of species and populations were presented in Table S2. The formal identification of all voucher specimens was identified by Professor Yafeng Wen according to Flora of China (1999) [27] and previous studies [42], and subsequently deposited in the herbarium of Central South University of Forestry and Technology. Accession numbers of the specimens are shown in Table S3. DNA was extracted from fresh samples by a modified 2×CTAB (Cetyltrimethylammonium Bromide) protocol [42]. The quality and integrity of DNA were evaluated by 1.0% agarose gels.
Morphological characters comparison among different species
Considering the wide applicability and popularity of Farjon’s classification system [22], we gathered five Taxus species to analyze their morphological characteristics, namely T. wallichiana, T. contorta, T. cuspidata, T. chinensis and T. mairei. Among them, T. contorta and T. mairei were most closely related to Qinling type, while T. chinensis and Qinling type are sympatrically distributed across the Qinling mountains [23, 25, 26]. These Taxus species were gathered from June to August in 2022, and their distribution is shown in Fig. S2. Six characters were employed, with their respective character states were defined in Fig. 2 and 3, and Table 2. For each species, ten mature leaves were selected for measuring leaf length (LL) and leaf width (LW, leaf blade at widest point). Ten measurements were taken for each mature leaf, resulting in a total of 100 measurements for each species. Length/width (L/W) is the ratio of leaf length (LL) to leaf width (LW). Considering the data did not conform to the normal distribution and homogeneity test of variance, significant differences among different Taxus species were analyzed using Kruskal-Wallis’s multiple range tests based on non-parametric statistics in “ggpubr” and “rstatix” package [57, 58]. Statistically significant differences between two Taxus species (p ≤ 5.0e-3) were indicated as “*” .
Microstructure observation of the leaf and lower epidermis
To observe the differences in leaf microstructure among the Taxus species, healthy and mature leaves were collected from the middle part of the tree and washed with distilled water. Square leaf samples with 3 mm×3 mm size were cut off at 1/3 distance from the base of the leaf, fixed by FAA fixative (containing 37% formaldehyde, glacial acetic acid, 70% ethyl alcohol; 1:1:18), labeled, and stored at 4℃ in the refrigerator. By the conventional paraffin sectioning method [59], the slices with a thickness of 10 μm were prepared as follows: leaf samples were dehydrated, transparentized, wax dipped, embedded, sectioned, dewaxed, rehydrated, stained (safranin O - fast green), dehydrated, transparentized (dimethylbenzene), and sealed (neutral balsam). Photographs of the sections were taken by Leica inverted microscope equipped with a 10× and 20× objective lens, respectively.
Scanning electron microscope (SEM) observation was used to make the stomatal band and papilla on the abaxial side more clearly. First, leaves were put into the stationary liquid of glutaric dialdehyde (2.5%). Then, the leaves were washed with PBS buffer (biosharp) and distilled water, dehydrated in different concentrations of alcohol (30%, 50%, 70%, 80%, 90%, 95%, 100%), transitioned to tert-butanol, and then freeze-dried. After drying, the samples were plated in a GVC-1000 ion sputterer for 20 min and then subjected to scanning electron microscopy (SEM) on a JSM-6380LV for sub-microscopic observation. Representative fields of view were selected to observe the leaf morphology and photographed for analysis.
Barcoding sequencing and sequence assembly
One nuclear internal transcribed spacer sequence (ITS) and two chloroplast DNA fragments (trnL-trnF and rbcL) were chosen as DNA barcoding to identify Taxus species by the results reported by Liu et al. [29](Table S4, [60,61,62]). PCR amplification was carried out in 20 µL reaction system containing 15 µL 2 × Taq PCR Master mix, 1 µL forward primer, 1 µL reverse primer, and 3 µL DNA template. PCR amplification procedures for each sequence are shown in Table S3. PCR products were inspected on 1.0% agarose gels, and purified using the SanPrep Column PCR Product Purification Kit, and the specific DNA fragments were recovered using the SanPrep Column DNA Gel Extraction Kit. Subsequently, the purified PCR products were sequenced in both directions on ABI3770XL DNA Sequencer (Applied Biosystems, Foster City, CA, USA). Furthermore, one sequence of Pseudotaxus chienii, three sequences of Qinling type, and thirty-two sequences of Taxus were downloaded from NCBI for phylogenetic analysis (Table S3). Consequently, this study incorporated Qinling type samples from the provinces of Shanxi, Henan, and Shaanxi.
Firstly, nucleotide sequences were aligned using MAFFT v7.471 with default parameters [63]. Secondly, the sequences were manually adjusted using MEGA X [64]. Thirdly, sequence fragments were made into the combined matrix by phylotools package (https://github.com/helixcn/phylotools). Afterwards, the incongruence length difference test (ILD) was carried out using PAUP 4.0 to determine whether the supermatrix was suitable for subsequent analyses [65, 66]. Finally, the saturation of the supermatrix was determined through substitution saturation test using DAMBE [67].
Phylogenetic relationship
The phylogenetic relationship was assessed based on the combination of the three sequences (ITS + trnL-trnF + rbcL). Bayesian inference (BI) analyses were conducted by Mrbayes 3.2.6 [68] with the GTRGAMMA nucleotide substitution model determined by PartitionFinder2 [69] in PhyloSuite software [70] according to the principle of minimum Akaike information criterion (AIC). Then, the maximum likelihood (ML) analyses were conducted using GTRGAMMA model in RAxML v 8.2.X with 1000 bootstrapping replicates [71]. Pseudotaxus chienii were rooted as the out-group for all the analysis. The cladogram tree was edited by the Evolview software [72].
Niche equivalency and similarity tests
19 bioclimatic variables and elevation were downloaded from WorldClim 2.1 (https://worldclim.org) with a spatial resolution of 2.5 min from 1970 to 2000. To avoid biased estimates of model coefficients and spurious significance levels resulting from multicollinearity, we excluded highly correlated climate variables based on Pearson’s correlation coefficient (|r|>0.75), finally, six climatic variables (Bio2, Bio3, Bio11, Bio15, Bio18 and elevation) were retained. Slope and aspect were obtained from Digital Elevation Model using the 3D analyst tools in the software ArcGIS 10.4. In addition, five soil variables were downloaded from the Harmonized World Soil Database (HWSD,https://www.fao.org/soils-portal/data-hub/soil-maps-anddatabases/harmonized-world-soil-database-v12/en/) and one human interference index (HII) from Socioeconomic Data and Applications Center (SEDAC, http://sedac.ciesin.columbia.edu) based on previous studies [33, 73, 74]. Finally, a total of 14 variables were retained for subsequent analysis (Fig. S3).
Ecological niche models and multivariate analysis in the “ecospat” package were used to quantify and compare the niches in geographic space and environmental space between two species (including subspecies and varieties) [16, 75, 76]. The degree of niche overlap between two species was estimated using the Schoener’s D index, ranging from 0 (no niche overlap) to 1 (perfect niche overlap) [16, 77]. Rödder et al. [78] proposed the following grading: 0-0.2 = no or very limited overlap, 0.2–0.4 = low overlap, 0.4–0.6 = moderate overlap, 0.6–0.8 = high overlap, 0.8-1.0 = very high overlap. Niche equivalency test evaluated the congruence of the two species’ geographic niches, whereas the niche similarity test was conducted to investigate the similarity between actual environmental niches and the expected ones [16]. The null hypothesis of niche equivalency is rejected when empirical values are significantly less than the critical values, indicating the two niches are not equivalent to each other [19, 79]. Hypotheses of niche similarity (i.e., the hypothesis of niche conservatism) are accepted when there is significant difference (P < 0.05) between the empirically observed Schoener’s D values (i.e., niche overlap) and the expected simulated overlap from the 100 pseudo-replicated datasets [18]. Furthermore, the contribution rate of driving factors to niche differentiation was calculated by the “factoextra” package [80].
Taxonomic treatment
Taxus qinlingensis
Y. F. Wen & X. T. Wu sp. nov. (秦岭红豆杉 is the Chinese name of Taxus qinlingensis.)
Diagnosis
This newly species resembles Taxus mairei in possessing similar leaf phenotype and elevation distribution, yet it is distinct due to the absence of papillate on midvein and tannis in the epidermis.
Type
CHINA. Shanxi Province: Xi’an City, Zhouzhi County, 108.23 N, 33.87 E, elev. 1100 m, July 23, 2019. Yafeng Wen WZZ201901 (holotype: CSFI).
Description
Evergreen tree to 25 m, trunk to 1.3 m DBH (diameter at breast height). Diaecious or monoecious. Flowering from February to April, fruiting from September to December. Bark thin, light red, purple-brown, or gray, split into strips or irregular flakes. Branches are numerous, ascending to erect, then spreading or drooping, forming a spreading, rounded, or pyramidal crown. Foliage branchlets irregularly alternate, terete with fine grooves alongside decurrent leaf bases, yellowish green turning green to brown. Winter buds are ovoid, scales persistent at base of shoots. Leaves distichously arranged, relatively loose on branchlets with leaf base twisted and spirally inserted, nearly sessile, thinned leathery, linear or falcate, 2.0–3.0 cm long, 2.5–3.5 mm wide, apex acute and aristulate, with revolute margins. Leaf adaxial side color was yellowish green to dark green with glossiness. Leaf abaxial side has two yellowish stomatal bands with an irregular arrangement, no papillae on the midvein, and both the midvein and the leaf edges color are shiny. Midvein length was 300–400 μm. Male cones axillary, solitary, ovoid with bracts at the base, with short stems under cones, and form rows on either side along the fertile shoots, yellowish green to yellow, 8–14 peltate microsporophylls, each with 4–6(-8) pollensacs. Female cones axillary, solitary, ovoid, and subsessile. Aril green at first, covering the lower part of the seed, swelling to succulent red (usually more translucent), and cover seeds, leaving its apex free, cup-like, 8–10 mm long, 7–10 mm wide. Seeds ovoid or obovoid, slightly flattened, with two obtuse ridges and a protruding apex, 5.0–8.0 mm long, 3.5–5.0 mm in diameter, green to brown or black.
Distribution and habitat
Taxus qinlingensis is mainly grown in deciduous broad-leaved forests on mountain slopes at altitudes between 500 and 1600 m in Shaanxi, Henan, Shanxi, and Hubei provinces. However, it can also found in the theropencedrymion of Yunnan and northwestern Sichuan (Danba, Songpan), where elevations ranging from 1700–2000 m and 2300–2500 m, respectively.
Specimens examined
CHINA. Shaanxi province: Zhouzhi county, evergreen mixed angiosperm/deciduous broad-leaved forest, 990–1150 m, 23 July 2019, M. Q. Wang & Y. F. Wen (Herbarium, Central South University of Forestry and Technology, CSFI); Weinan city, mixed angiosperm/conifer forest hill, 1450–1500 m, 1 August 2019, M. Q. Wang & X. T. Wu (CSFI); Zhashui county, mixed angiosperm/conifer forest hill, 900–1000 m, 18 April 2021, Z. P. Liao (CSFI); Foping county, mixed angiosperm/conifer forest hill, 900–1000 m, 3 March 2022, Z. P. Liao (CSFI); Lantian county, mixed angiosperm/conifer forest hill, 1000–1200 m, 31 March 2022, Z. P. Liao (CSFI). Gansu province: Liangdang county, mixed angiosperm/conifer forest on hill, 1100–1200 m, 30 July 2019, M. Q. Wang & Y. F. Wen (CSFI). Henan province: Lushi county, evergreen mixed angiosperm/deciduous broad-leaved forest, 1000–1300 m, 12 December 2022, X. T. Wu (CSFI).
Conservation status
Endangered (EN A2cd; C1). The entire area of occupancy (AOO) of this species in nature is around 164 km2. Some are distributed around the village, known as the fengshui trees, while others are distributed around the mountains that are not easy to access. The total number is less than 30,000. Taxus exhibits weak regeneration ability and poor habitat quality. Consequently, we proposed that T. qinlingensis should be classified as Endangered (EN) by IUCN standards [81,82,83].
Key to the Taxus species
1 midvein without papillae | 2 |
---|---|
1 midvein with densely and evenly distributed papillae | 3 |
2 Leaves thinned leathery, linear, 2.0–3.0 cm × 2.5–3.5 mm, leaf margins parallel, apex acute, midvein color different from stomatal band color, leaf edges shiny on abaxial leaf surface, with relatively low leaf length/width ratio (6.50–8.06, mean > 7). | T. qinlingensis |
2 Leaves leathery, linear and straight, 2.0–4.0 cm × 2.0–3.5 mm, leaf edges involution, midvein convex, apex mucronate, midvein color different from stomatal band color, leaf edges shiny on abaxial leaf surface, with relatively high leaf length/width ratio (9.00-10.25, mean > 9). | T. cuspidata |
3 Leaves leathery, linear, 1.5–2.5 cm × 2.5–3.5 mm, leaf margins parallel, apex acute, midvein color same as stomatal band color, leaf edges not shiny on abaxial leaf surface, with relatively low leaf length/width ratio (5.87–7.20, mean > 6) | T. chinensis |
3 Leaves texture relatively hard, falcate, 1.0–2.5 cm × 1.5–3.0 mm, leaf margins parallel, apex acuminate, midvein color different from stomatal band color, leaf edges shiny on abaxial leaf surface, with relatively high leaf length/width ratio (5.00–13.00, mean > 8) | T. wallichiana |
3 Leaves thinned leathery, sword-shaped, 1.2–2.5 cm × 1.0–3.5 mm, leaf margins parallel, apex acute, aristate, midvein color same as stomatal band color, leaf edges not shiny on abaxial leaf surface, with relatively high leaf length/width ratio (8.00–25.00, mean > 10) | T. contorta |
3 Leaves thinned leathery, slightly falcate, 2.0–3.5 cm × 2.5–4.0 mm, leaf margins parallel, apex acuminate, midvein color different from stomatal band color, leaf edges shiny on abaxial leaf surface, with relatively low leaf length/width ratio (5.87–7.15, mean > 6) | T. mairei |
Data availability
The original phenotypic datasets in current study are available in Table S5. All sequences (ITS, trnL - trnF and rbcL sequences) used in this study have been submitted to the National Center for Biotechnology Information (NCBI, https://www.ncbi.nlm.nih.gov/) with accession numbers (MW788517 - MW788520, MZ208838 - MZ208861, and OQ891381 - OQ891399 for ITS sequence; MW792129 - MW792132, MZ220779 - MZ220802, and OQ978853 - OQ978871 for trnL - trnF sequences; MW792092 - MW792095, MW893342 - MW893365, and OQ913085 - OQ913103 for rbcL sequence) (Table S3). All the sequences have been available online. Information for other samples used for phylogenetic analysis download from NCBI can be found in Table S3.
References
Morigengaowa, Shang H, Liu BD, Kang M, Yan YH. One or more species? GBS sequencing and morphological traits evidence reveal species diversification of Sphaeropteris brunoniana in China. Biodivers Sci. 2019;27(11):1196–204. https://doi.org/10.17520/biods.2019146.
Orr RJS, Haugen MN, Berning B, Bock P, Cumming RL, Florence WK, et al. A genome-skimmed phylogeny of a widespread bryozoan family, Adeonidae. BMC Evol Biol. 2019;19(1):235. https://doi.org/10.1186/s12862-019-1563-4.
Mace GM. The role of taxonomy in species conservation. Philos Trans R Soc Lond B Biol Sci. 2004;359(1444):711–9. https://doi.org/10.1098/rstb.2003.1454.
Gutiérrez EE, Helgen KM. Mammalogy: outdated taxonomy blocks conservation. Nature. 2013;495(7441):314. https://doi.org/10.1038/495314e.
Cunha RL, Forsman ZH, Belderok R, Knapp ISS, Castilho R, Toonen RJ. Rare coral under the genomic microscope: timing and relationships among Hawaiian Montipora. BMC Evol Biology. 2019;19(1):153. https://doi.org/10.1186/s12862-019-1476-2.
Conix S. Taxonomy and conservation science: interdependent and value-laden. Hist. Philos Life Sci. 2019;41(2):15. https://doi.org/10.1007/s40656-019-0252-3.
Zachos F. Species splitting puts conservation at risk. Nature. 2013;494:35. https://doi.org/10.1038/494035c.
Hu Y, Thapa A, Fan H, Ma T, Wu Q, Ma S. Genomic evidence for two phylogenetic species and long-term population bottlenecks in red pandas. Sci Adv. 2020;6(9):eaax5751. https://doi.org/10.1126/sciadv.aax5751.
Liu JQ. The integrative species concept and species on the speciation way. Biodivers Sci. 2016;24(9):1004–8. https://doi.org/10.17520/biods.2016222.
Wilkins JS, Species. A history of the idea. Berkeley: University of California Press; 2009. p. 320.
Lin X, Shih C, Hou Y, Shu X, Zhang M, Hu J, et al. Climatic-niche evolution with key morphological innovations across clades within Scutiger boulengeri (Anura: Megophryidae). Ecol Evol. 2021;11(15):10353–68. https://doi.org/10.1002/ece3.7838.
Metcalfe CR, Chalk L. Anatomy of the dicotyledons: wood structure and conclusion of the general introduction. Volume II. New York: Oxford University Press; 1983.
Kolb KL, Gomes SMA, Lombardi JA. Leaf anatomy as a taxonomy tool for the identification of Brazilian native species of Chionanthus (Oleaceae). Flora. 2020;266:151590. https://doi.org/10.1016/j.flora.2020.151590.
Wiens JJ, Graham CH. Niche conservatism: integrating evolution, ecology, and conservation biology. Annu Rev Ecol Evol Syst. 2005;36:519–39. https://doi.org/10.1146/annurev.ecolsys.36.102803.095431.
De Queiroz K. Species concepts and species delimitation. Syst Biol. 2007;56:879–86. https://doi.org/10.1080/10635150701701083.
Broennimann O, Fitzpatrick MC, Pearman PB, Petitpierre B, Pellissier L, Yoccoz NG, et al. Measuring ecological niche overlap from occurrence and spatial environmental data. Glob Ecol Biogeogr. 2012;21(4):481–97. https://doi.org/10.1111/j.1466-8238.2011.00698.x.
Sun Y, Surget-Groba Y, Gao SX. Divergence maintained by climatic selection despite recurrent gene flow: a case study of Castanopsis carlesii (Fagaceae). Mol Ecol. 2016;25(18):4580–92. https://doi.org/10.1111/mec.13764.
Moreno-Contreras I, Sánchez-González LA, Arizmendi MC, Prieto-Torres DA, Navarro-Sigüenza AG. Climatic niche evolution in the Arremon brunneinucha Complex (Aves: Passerellidae) in a Mesoamerican Landscape. Evol Biol. 2020;47:123–32. https://doi.org/10.1007/s11692-020-09498-7.
Hong DY. Biodiversity pursuits need a scientific and operative species concept. Biodivers Sci. 2016;24(9):979–99. https://doi.org/10.17520/biods.2016203.
Wang MZ, Fan XK, Zhang YH, Wu J, Mao LM, Zhang SL, et al. Phylogenomics and integrative taxonomy reveal two new species of Amana (Liliaceae). Plant Divers. 2023;45(1):54–68. https://doi.org/10.1016/j.pld.2022.03.001.
Spjut RW. Taxonomy and nomenclature of Taxus (Taxaceae). J Bot Res Inst Tex. 2007;1(1):203–89. https://biostor.org/reference/157469.
Farjon A. A handbook of the world’s conifers (volume 1, second revised edition). Boston; Brill Academic; 2017.
Gao LM, Möller M, Zhang XM, Hollingsworth ML, Liu J, Mill RR, et al. High variation and strong phylogeographic pattern among cpDNA haplotypes in Taxus Wallichiana (Taxaceae) in China and North Vietnam. Mol Ecol. 2007;16(22):4684–98. https://doi.org/10.1111/j.1365-294X.2007.03537.x.
Liu J, Möller M, Gao LM, Zhang DQ, Li DZ. DNA barcoding for the discrimination of eurasian yews (Taxus L., Taxaceae) and the discovery of cryptic species. Mol Ecol Resour. 2011;11(1):89–100. https://doi.org/10.1111/j.1755-0998.2010.02907.x.
Fu CN, Wu CS, Ye LJ, Mo ZQ, Liu J, Chang YW, et al. Prevalence of isomeric plastomes and effectiveness of plastome super-barcodes in yews (Taxus) worldwide. Sci Rep. 2019;9(1):2773–83. https://doi.org/10.1038/s41598-019-39161-x.
Möller M, Liu J, Li Y, Li JH, Ye LJ, Mill R, et al. Repeated intercontinental migrations and recurring hybridizations characterize the evolutionary history of yew (Taxus L). Mol Phylogenet Evol. 2020;153:106952. https://doi.org/10.1016/j.ympev.2020.106952.
Fu LG, Li N, Mill R, Taxaceae In: Flora of China. Science Press, Beijing and, Garden MB. St. Louis, 1999; 4: 89–96.
Möller M, Gao LM, Mill RR, Liu L, Zhang DQ, Poudel RC, et al. A multidisciplinary approach reveals hidden taxonomic diversity in the morphologically challenging Taxus wallichiana complex. Taxon. 2013;62(6):1161–77. https://doi.org/10.12705/626.9.
Liu J, Milne RI, Möller M, Zhu GF, Ye LJ, Luo YH, et al. Integrating a comprehensive DNA barcode reference library with a global map of yews (Taxus L.) for forensic identification. Mol Ecol Resour. 2018;18:1115–31. https://doi.org/10.1111/1755-0998.12903.
Poudel RC, Möller M, Gao LM, Ahrends A, Baral SR, Liu J, et al. Using morphological, molecular and climatic data to delimitate yews along the Hindu Kush-Himalaya and adjacent regions. PLoS ONE. 2012;7(10):e46873. https://doi.org/10.1371/journal.pone.0046873.
Möller M, Gao LM, Mill RR, Li DZ, Hollingsworth ML, Gibby M. Morphometric analysis of the Taxus Wallichiana complex (Taxaceae) based on herbarium material. Bot J Linn Soc. 2007;155(3):307–35. https://doi.org/10.1111/j.1095-8339.2007.00697.x.
Hu ZH. Plant anatomy. Beijing: Higher Education Press; 2010.
Wu XT, Wang MQ, Li XY, Yan YD, Dai MJ, Xie WY, et al. Response of distribution patterns of two closely related species in Taxus genus to climate change since last inter-glacial. Ecol Evol. 2022;12(9):e9302. https://doi.org/10.1002/ece3.9302.
Font M, Garcia-Jacas N, Vilatersana R, Roquet C, Susanna A. Evolution and biogeography of Centaurea section acrocentron inferred from nuclear and plastid DNA sequence analyses. Ann Botany. 2009;103(6):985–97. https://doi.org/10.1093/aob/mcp022.
Renoult JP, Kjellberg F, Grout C, Santoni S, Khadari B. Cyto-nuclear discordance in the phylogeny of Ficus section Galoglychia and host shifts in plant-pollinator associations. BMC Evol Biol. 2009;9:248–67. https://doi.org/10.1186/1471-2148-9-248.
Nelson TC, Crandall JG, Ituarte CM, Catchen JM, Cresko WA. Selection, linkage, and population structure interact to shape genetic variation among three spine stickleback genomes. Genetics. 2019;212(4):1367–82. https://doi.org/10.1534/genetics.119.302261.
Ducret V, Schaub M, Goudet J, Roulin A. Female-biased dispersal and non-random gene flow of MC1R variants do not result in a migration load in barn owls. Heredity. 2019;122(3):305–14. https://doi.org/10.1038/s41437-018-0115-9.
Liu J, Möller M, Provan J, Gao LM, Poudel CR, Li DZ. Geological and ecological factors drive cryptic speciation of Yews in a Biodiversity Hotspot. New Phytol. 2013;199:1093–108. https://doi.org/10.1111/nph.12336.
Wu H, Wang ZF, Zhang YX, Frantz L, Roos C, Irwin DM, et al. Hybrid origin of a primate, the gray snub-nosed monkey. Science. 2023;380:eabl4997. https://doi.org/10.1126/science.abl4997.
Zhang BW, Xu LL, Li N, Yan PC, Jiang XH, Woeste KE, et al. Phylogenomics reveals an ancient hybrid origin of the Persian Walnut. Mol Biol Evol. 2019;36(11):2451–61. https://doi.org/10.1093/molbev/msz112.
Qin HT, Mӧller M, Milne R, Luo YH, Zhu GF, Li DZ, et al. Multiple paternally inherited chloroplast capture events associated with Taxus speciation in the Hengduan Mountains. Mol Phylogenet Evol. 2023;189:107915. https://doi.org/10.1016/j.ympev.2023.107915.
Wu XT. Genetic divergence and local adaption for Taxus Chinensis and Taxus mairei. Changsha: Central South University of Forestry and Technology; 2021.
Yao H, Song J, Liu C, et al. Use of ITS2 region as the universal DNA barcode for plants and animals. PLoS ONE. 2010;5:e13102. https://doi.org/10.1371/journal.pone.0013102.
Soledispa P, Santos-Ordóñez E, Miranda M, et al. Molecular barcode and morphological analysis of Smilax Purhampuy Ruiz. Ecuador. Peer J. 2001;9:e11028. https://doi.org/10.7717/peerj.11028.
Antil S, Abraham JS, Sripoorna S, Maurya S, Dagar J, Makhija S, et al. DNA barcoding, an effective tool for species identification: a review. Mol Biol Rep. 2023;50:761–75. https://doi.org/10.1007/s11033-022-08015-7.
IUCN. The IUCN Red List of Threatened Species. 2022a. https://www.iucnredlist.org.
Natri HM, Merilä J, Shikano T. The evolution of sex determination associated with a chromosomal inversion. Nat Communication. 2019;10:145. https://doi.org/10.1038/s41467-018-08014-y.
Chen JH, Huang Y, Brachi B, Yun QZ, Zhang W, Lu W, et al. Genome-wide analysis of Cushion willow provides insights into alpine plant divergence in a biodiversity hotspot. Nat Communication. 2019;10(1):5230–42. https://doi.org/10.1038/s41467-019-13128-y.
Fu Y, Xu Z, Wen B, Gao J, Chen Z. Gonad-specific transcriptomes reveal Differential expression of Gene and miRNA between male and female of the Discus Fish (Symphysodon Aequifasciatus). Front Physiol. 2020;11:754. https://doi.org/10.3389/fphys.2020.00754.
Jiang Y, Duan Y, Yin J, Ye S, Zhu J, Zhang F, et al. Genome-wide identification and characterization of the Populus WRKY transcription factor family and analysis of their expression in response to biotic and abiotic stresses. J Exp Bot. 2014;65(22):6629–44. https://doi.org/10.1093/jxb/eru381.
Liu C, Xie T, Chen C, Luan A, Long J, Li C, Ding Y, He Y. Genome-wide organization and expression profiling of the R2R3-MYB transcription factor family in pineapple (Ananas comosus). BMC Genomics. 2017;18(1):503. https://doi.org/10.1186/s12864-017-3896-y.
Cheng J, Wang X, Liu X, Zhu X, Li Z, Chu H, et al. Chromosome-level genome of Himalayan yew provides insights into the origin and evolution of the paclitaxel biosynthetic pathway. Mol Plant. 2021;14(7):1199–209. https://doi.org/10.1016/j.molp.2021.04.015.
Song C, Fu F, Yang L, Niu Y, Tian Z, He X, et al. Taxus yunnanensis genome offers insights into gymnosperm phylogeny and taxol production. Communication Biology. 2021;4:1203. https://doi.org/10.1038/s42003-021-02697-8.
Xiong X, Gou J, Liao Q, Li Y, Zhou Q, Bi G, et al. The Taxus genome provides insights into paclitaxel biosynthesis. Nat Plants. 2021;7:1026–36. https://doi.org/10.1038/s41477-021-00963-5.
Harrison S, Damschen E, Fernandez-Going B, Eskelinen A, Copeland S. Plant communities on infertile soils are less sensitive to climate change. Ann Botany. 2015;116(6):1017–22. https://doi.org/10.1093/aob/mcu230.
Wen YF, Uchiyama K, Ueno S, Han WJ, Xie WD, Tsumura Y. Assessment of the genetic diversity and population structure of Maire yew (Taxus chinensis var. Mairei) for conservation purposes. Can J for Res. 2018;48(5):589–98. https://doi.org/10.1139/cjfr-2017-0429.
Kassambara K, rstatix. Pipe-Friendly Framework for Basic Statistical Tests; 2022a, https://CRAN.R-project.org/package=rstatix.
Kassambara K. ggpubr: ‘ggplot2’ Based Publication Ready Plots; 2022b. https://CRAN.R-project.org/package=ggpubr.
Li ZL. Plant slicing technology (Second Edition). Beijing: Science Press;1978.
Fay MF, Swensen SM, Chase MW. Taxonomic affinities of Medusagyne oppositifolia (Medusagynaceae). Kew Bull. 1997;52(1):111–20. https://doi.org/10.2307/4117844.
Taberlet P, Gielly L, Pautou G, Bouvet J. Universal primers for amplification of three non-coding regions of chloroplast DNA. Plant Mol Biol. 1991;17:1105–9. https://doi.org/10.1007/BF00037152.
Yin WX, Zhang MZ, Lou CJ, Wang JZ, Yu WT, Hu ML, et al. Identification of species of Taxus by DNA Barcoding based on ITS sequence. J Zhejiang Forestry Sci Technol. 2020;40(2):111–5. https://doi.org/10.3969/j.issn.1001-3776.2020.02.017.
Yamada N, Katoh T, Tomii K, Katoh K. Parallelization of MAFFT for large-scale multiple sequence alignments. Bioinformatics. 2018;34(14):2490–2. https://doi.org/10.1093/bioinformatics/bty121.
Hall BG. Building phylogenetic trees from molecular data with MEGA. Mol Biol Evol. 2013;30(5):1229–35. https://doi.org/10.1093/molbev/mst012.
Dolphin K, Belshaw R, Orme CDL, Quicke DLJ. Noise and incongruence: interpreting results of the incongruence length difference test. Mol Phylogenet Evol. 2000;17(3):401–6. https://doi.org/10.1006/mpev.2000.0845.
Swofford DL. PAUP*: phylogenetic analysis using parsimony. Sunderland (MA). Sinauer Associates; 2003.
Xia X. DAMBE7: New and improved tools for data analysis in molecular biology and evolution. Mol Biol Evol. 2018;35(6):1550–2. https://doi.org/10.1093/molbev/msy073.
Ronquist F, Teslenko M, van der Mark P, Ayres DL, Höhna S, Larget B, et al. MrBayes 3.2: efficient bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61(3):539–42. https://doi.org/10.1093/sysbio/sys029.
Lanfear R, Frandsen PB, Wright AM, Senfeld T, Calcott B. PartitionFinder 2: new methods for selecting partitioned models of evolution for molecular and morphological phylogenetic analyses. Mol Biol Evol. 2016;34(3):772–3. https://doi.org/10.1093/molbev/msw260.
Zhang D, Gao FL, Li WX, Jakovlić I, Zou H, Zhang J, et al. PhyloSuite: an integrated and scalable desktop platform for streamlined molecular sequence data management and evolutionary phylogenetics studies. Mol Ecol Resour. 2020;20(1):348–55. https://doi.org/10.1111/1755-0998.
Stamatakis, A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312-13. https://doi.org/10.1093/bioinformatics/btu033.
Subramanian B, Gao SH, Lercher MJ, Hu SN, Chen WH. Evolview v3: a webserver for visualization, annotation, and management of phylogenetic trees. Nucleic Acids Res. 2019;47(W1):W270–5. https://doi.org/10.1093/nar/gkz357.
Gallardo B, Zieritz A, Aldridge DC. The importance of the human footprint in shaping the global distribution of terrestrial, freshwater and marine invaders. PLoS ONE. 2015;10(5):e0125801. https://doi.org/10.1371/journal.pone.0125801.
Morrison KD, Hammer E, Boles O, Madella M, Whitehouse N, Gaillard MJ, et al. Mapping past human land use using archaeological data: a newclassification for global land use synthesis and data harmonization. PLoS ONE. 2021;16(4):e0246662. https://doi.org/10.1371/journal.pone.0246662.
Di Cola V, Broennimann O, Petitpierre B, Breiner FT, D’Amen M, Randin C, et al. Ecospat: an R package to support spatial analyses and modeling of species niches and distributions. Ecography. 2017;40(6):774–87. https://doi.org/10.1111/ecog.02671.
Dai X, Li X, Huang Y, Liu X. The speciation and adaptation of the polyploids: a case study of the Chinese Isoetes L. diploid-polyploid complex. BMC Evol Biol. 2020;20(1):118. https://doi.org/10.1186/s12862-020-01687-4.
Schoener TW. Nonsynchronous Spatial Overlap of Lizards in patchy habitats. Ecology. 1970;51:408–18. https://doi.org/10.2307/1935376.
Rödder D, Engler JO. Quantitative metrics of overlaps in Grinnellian niches: advances and possible drawbacks. Glob Ecol Biogeogr. 2011;20(6):915–27. https://doi.org/10.1111/j.1466-8238.2011.00659.x.
Warren DL, Glor RE, Turelli M. ENMTools: a toolbox for comparative studies of environmental niche models. Ecography. 2010;33:607–11. https://doi.org/10.1111/j.1600-0587.2009.06142.x.
Kassambara A, Mundt F, factoextra. Extract and visualize the results of multivariate data analyses. R package version 1.0.7. 2020. Retrieved from https://CRAN.R-project.org/package=factoextra.
IUCN. IUCN Red List Categories and Criteria, Version 3.1, Gland, Switzerland and Cambridge, UK. 2012b.
Xie Y. Progress and application of IUCN Red list of threatened species. Biodivers Sci. 2022;30(10):22445. https://doi.org/10.17520/biods.2022445.
Yang Y. An updated red list assessment of gymnosperms from China (Version 2021). Biodivers Sci. 2021;29:1599–606. https://doi.org/10.17520/biods.2021342.
Acknowledgements
We are grateful to Xiaohang Min for providing samples of Taxus contorta, which was planted in Tibet Academy of Forest Sciences. The provenance of T. contorta was Jilong, Xizang.
Funding
This study was supported by the National Natural Science Foundation of China under Grant (Grant No. 31470666), Hunan Provincial Innovation Foundation for Postgraduate (Grant No. CX2018B435), and Scientific Innovation Fund for Graduate of Central South University of Forestry and Technology (Grant No. 20181009). Discipline of China Forestry Bureau (Grant No. 2016 21), Double First-class Initiative Cultivation Discipline in Hunan Province (Grant No. 2018469). National Research Base for Landscape Architecture of Qingxiushan, Nanning, Guangxi (Grant No. 202196).
Author information
Authors and Affiliations
Contributions
WYF, WXT contributed to the conception of the study. WMQ, CY, LZP collected the plant material. WXT, WMQ, LXY performed the experiment. WXT, LXY contributed significantly to the analysis and manuscript preparation. WXT performed the data analyses and wrote the manuscript. WYF, WS, ZDL helped revise the paper with constructive discussions. All authors read and approved the manuscript.
Corresponding authors
Ethics declarations
Ethics approval and consent to participate
The collection of Taxus had been permitted by natural reserves or botanical garden. Collection of plant material complied with institutional or national guidelines and was conducted in accordance with local legislation. No materials from animal or human were used in this research.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Wu, X., Wang, M., Li, X. et al. Identification and characterization of a new species of Taxus — Taxus qinlingensis by multiple taxonomic methods. BMC Plant Biol 24, 658 (2024). https://doi.org/10.1186/s12870-024-05338-4
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12870-024-05338-4