Next Article in Journal
Effects of Prescribed Fire, Site Factors, and Seed Sources on the Spread of Invasive Triadica sebifera in a Fire-Managed Coastal Landscape in Southeastern Mississippi, USA
Previous Article in Journal
The Effects of DNA Methylation Inhibition on Flower Development in the Dioecious Plant Salix Viminalis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Torreya jackii (Taxaceae): A Special Species that is Genetically Admixed, Morphologically Distinct, and Geographically Sympatric with Parent Species

1
State Key Laboratory of Grassland Agro-Ecosystems, School of Life Sciences, Lanzhou University, Lanzhou 730000, China
2
Laboratory of Subtropical Biodiversity, School of Agricultural Sciences, Jiangxi Agricultural University, Nanchang 330045, China
*
Author to whom correspondence should be addressed.
Forests 2019, 10(2), 174; https://doi.org/10.3390/f10020174
Submission received: 14 December 2018 / Revised: 25 January 2019 / Accepted: 18 February 2019 / Published: 19 February 2019
(This article belongs to the Section Forest Ecology and Management)

Abstract

:
Torreya jackii Chun is an endangered species (Taxaceae) confined to a few localities in China. However, the species status of T. jackii within Torreya Arn. has not been clearly elucidated under a phylogenetic context. In this study, phylogenetic analyses based on the nuclear internal transcribed spacer (ITS) and amplified fragment length polymorphism (AFLP) indicated that T. jackii is closely related with a sympatric species T. grandis Fort. ex Lindl. that is present due to cultivation. However, analysis based on the concatenated sequences of seven chloroplast loci resolved T. jackii as the first branch within the genus. Given their overlapping distribution and synchronous blooming, we suggest that the plastid-nuclear incongruence was derived from the dilution of the nuclear genome of T. jackii by T. grandis via pollen-mediated introgression hybridization when the two species met due to cultivation. Introgressive hybridization is fairly common in plants but few cases have been recognized as independent species. Our study highlights the complexity of protecting endangered species and the need for caution to prevent the unreasonable expansion of economic crops into the distribution ranges of their wild relatives.

1. Introduction

Torreya Arn. is a genus in Taxaceae, which is distinguished by its completely enclosed seed within an aril and the production of two axillary female cones [1,2,3]. Fossils of this genus are found widely across the northern hemisphere and the earliest has been dated to a Jurassic deposit in Europe [4]. However, only six highly endangered members are extant and they have disjunct distributions in North America (NA) and eastern Asia (EA) [2]. Two species, T. taxifolia Arn. and T. californica Torr., are endemic in NA, while four species are present in EA comprising T. nucifera (L.) Sieb. et Zucc., T. fargesii Franch., T. jackii Chun, and T. grandis Fort. ex Lindl. The first three species in EA are confined to a few localities in Japan or China, and the latter is cultivated due to its edible seed and oil in some southern provinces of China, where it overlaps the range with T. jackii (Figure 1) [1].
The morphological divisions within Torreya are not congruent with the geographical distribution of the species. This genus is divided into two sections in the only proposed infrageneric system [5,6]. Sect. Ruminatae Hu is defined by deeply ruminated albumen and it contains two species from EA and one species from NA, whereas sect. Nuciferac Hu has only slightly ruminated albumen and it contains the remaining three species. It is unusual that species from the same continent are not morphologically approximated. Furthermore, despite their sympatry, T. jackii and T. grandis are morphologically distinct and ascribed to different sections due to differences in the albumen [5]. Moreover, their leaves differ in length, where those of T. jackii are at least 10 cm whereas those of T. grandis are less than 6 cm [1].
The first comprehensive molecular phylogeny of Torreya reconstructed based on the nuclear internal transcribed spacer (ITS) contradicts the morphologically infrageneric division, but it is highly congruent with the geographical distribution [7] where the genera were resolved into two clades distributed in NA and EA. In addition, T. jackii and T. grandis were found to be closely related [7]. However, a subsequent study based on chloroplast (cp) DNA obtained different phylogenies by resolving T. jackii as the most basal clade [8]. Moreover, a second analysis based on ITS sequences that were mostly identical to those employed in the previous study found a different phylogeny [7,8]. These discrepancies in terms of the geographical distribution, morphological division, genetic information from different genomes, and analyses based on different methods or samples suggest that the genus Torreya seems to have a complex history, especially the conflict of phylogenetic position of T. jackii in the genus.
In the present study, compared with previous studies [7,8] we increased the sample size and analyzed more genetic information from both the chloroplast and nuclear genomes in order to explore the phylogenetic relationships within Torreya, mainly to elucidate the underling mechanism responsible for the incongruent phylogenetic context of T. jackii according to different lines of evidence.

2. Materials and Methods

2.1. Sample Collections

Leaf samples were collected from at least two individuals of all six species of Torreya. T. jackii and T. grandis were the main focus of the present study, and we collected six individuals from four localities and three individuals from three localities for these two species, respectively. Amentotaxus argotaenia (Hance) Pilger and A. yunnanensis Li were selected as outgroup based on previous studies [9]. Details of the samples, their voucher information, and GenBank accession numbers are listed in Table 1.

2.2. DNA Extraction, PCR Amplification and Sequencing

DNA was extracted from approximately 20 mg of leaves dried in silica gel using the cetyltrimethylammonium bromide (CTAB) method [10]. Seven cp gene loci, rbcL, matK, trnL-trnF, psbB1-psbB2, rpl16, trnS-trnG, and trnD-trnT, and the nuclear internal transcribed spacer (ITS) were amplified with previously reported primers [11,12,13,14] according to the protocol described in our previous study [15]. PCR products from the ITS was purified using an Agarose Gel DNA Purification kit (TakaRa Inc., Dalian, China) and then cloned with the pMD19-T vector (TakaRa Inc., Dalian, China) according to the recommended protocol, before transformation into competent Escherichia coli JM109. The transformed bacteria were screened overnight on solid Luria–Bertani medium containing 100 mg/mL ampicillin at 37 °C. Five or more positive clones were amplified and sequenced using the universal primers M13-47 and RV-M. For all the cp loci, the primers used for amplification were employed for sequencing, which was conducted at the Beijing Genomics Institute. New sequences have been deposited in GenBank under the accession numbers listed in Table 1.
Amplified fragment length polymorphism (AFLP) analysis was conducted according to the protocol described in our previous study [16]. Selective amplification was performed using eight pairs of primer combinations, i.e., EcoRI-GTG/MseI-CAA, EcoRI-GAG/MseI-CTA, EcoRI-GTG/MseI-CTA, EcoRI-GTG/MseI-CTC, EcoRI-GAG/MseI-CTG, EcoRI-GAT/MseI-CTG, EcoRI-GTG/MseI-CTG, and EcoRI-GTG/MseI-CTT. Band data were collected using GeneScan 2.1 (Applied Biosystems, Foster City, CA, USA).

2.3. Data Analyses

Three datasets were constructed where one comprised the nuclear ITS sequences, the second contained the concatenated sequences of the seven cpDNA regions, and the third comprised the AFLP bands. Multiple alignments of sequence data were obtained using ClustalX 1.83 [17] and refined manually with MEGA 5.0 [18]. MEGA was also used to calculate the genetic distances according to the Kimura two-parameter model. The phylogenetic relationship was constructed using the maximum parsimony (MP) and maximum likelihood (ML) models implemented in PAUP 4.0 [19], and Bayesian inference (BI) in MrBayes 3.1.2 [20]. MP analyses were conducted by employing heuristic search to obtain the starting tree with stepwise, tree-bisection-reconnection (TBR) branch swapping, and the options of steepest descent, MulTrees, and Collapse. Bootstrap supports (BPs) were calculated by bootstrap analysis with 1000 replicates and the same settings as above. The optimal substitution models for the ML and BI analyses were found with ModelTest [21] based on Akaike’s information criterion. The GTR+G+I model was selected for the seven combined cpDNA regions and the GTR+G model for the nuclear ITS gene. ML analyses were implemented via heuristic searches with 1000 replicates of random sequence addition, TBR branch swapping, and MulTrees. BPs were obtained using 1000 replicates by heuristic search with the same options. For BI, one cold and three hot Monte Carlo Markov chains were run twice for 2,000,000 generations, with sampling every 100 generations. Tracer v.1.5 (http://tree.bio.ed.ac.uk/software/tracer/) was used to choose a suitable burn-in period. PAUP program was used to calculate a consensus tree and posterior probabilities from the sampled trees after the burn-in period.
The multi-locus profiles of amplified fragment length polymorphism (AFLP) were scored for the presence (1) or absence (0) of fragments measuring between 50 bp and 300 bp. A neighbor-joining tree was constructed based on the Nei and Li distance [22] from the presence or absence matrix using the TREECON 1.3b program [23]. BPs obtained from 1000 pseudo-replicates were used as a measure of confidence in the reconstructed tree topology.

3. Results

The aligned lengths and variations in all the studied loci are listed in Table 2. The combined dataset of seven cp regions comprised 8241 bp, where 487 were parsimony-informative. The average genetic distances between species are listed in Table 3 (upper right). The overall genetic distance was 0.28% while those within T. grandis and T. jackii were both 0. One to two different ITS sequences were found in each species, with a total of 33 sequences and a combined length of 1152 bp, where 212 were parsimony-informative. The average genetic distances between species are listed in Table 3 (lower left). The overall genetic distance was 1.17% while those within T. grandis and T. jackii were 0.06% and 0.32%, respectively. The AFLP dataset comprised 813 bands and 623 were variable, where 507 were parsimony-informative.
The MP strict consensus tree reconstructed with cpDNA datasets was largely consistent with those obtained by ML and BI in terms of the topology, as shown in Figure 2, where the support values from all three methods are indicated above the branches. The species T. jackii from EA branched off first with strong support, followed by the two NA species, T. taxifolia and T. californica, and the remaining Torreya species were resolved into two clades, which one contained T. yunnanensis and T. fargesii with strong support, and the other contained T. grandis and T. nucifera with moderate support.
The MP strict consensus tree (Figure 3) based on the nuclear ITS dataset recovered the same topology as the ML and BI trees. Compared with the tree obtained using the cp gene loci, the major difference was the position of T. jackii, which appeared in the base position according to the cp phylogeny, whereas it was closely related to T. grandis in the tree based on the ITS phylogeny. The AFLP data recovered an almost identical topology (Figure 4) to that using the ITS sequences, except T. grandis was resolved into polyphyly.

4. Discussion

Our phylogenetic analyses based on seven cp loci showed that T. jackii is the most basal branch of Torreya. The remaining species were resolved into two clades, where one comprised two species from NA and the other contained three Asian species, T. grandis, T. fargesii, and T. nucifera (Figure 2). Both ITS and AFLP, which are mostly representative of the nuclear genome, recovered a largely similar phylogeny (Figure 3), except that T. jackii was resolved as closely related to T. grandis. Among the six T. jackii sequences, two were identical to that in T. grandis whereas the other four differed by a distance of only 0.32%, which was much smaller than the overall mean distance (1.17%) for ITS sequences (Table 3). Incongruent phylogenies were also reported in a previous study based on different samples and analyses [7,8], thereby suggesting that technical issues such as insufficient taxon sampling, long branch attraction, or sequencing errors might not explain the incongruence [24].
A number of biological factors could explain the incongruence between the phylogenies based on nuclear and cp DNA, including convergent evolution, lineage sorting, and introgressive hybridization [25]. Given our wide range of samples and the fact that we analyzed the cp and nuclear genomes, the first two explanations are unlikely to apply to T. jackii, but the last is difficult to reject [24]. Introgressive hybridization is a process where genes are transferred via the formation of an initial F1 hybrid that subsequently crosses with individuals from one or both of the parental species [26,27]. The occurrence of this process requires that the two species meet some requirements including contact in space, overlapping blossoming periods, the same number of chromosomes, and the maternal inherence of cp DNA [26]. We found that T. jackii and T. grandis may satisfy all these requirements. First, T. jackii inhabits southern China where T. grandis is cultivated (Figure 1). In most localities, the two species are well within the distance of pollen dispersal, which is mediated by the wind in both species. Second, the male plants from both species blossom from late March until middle April, and the females from late March until early April [28,29,30]. Third, both species possess 11 chromosomes including one sex chromosome [31]. Finally, maternal inheritance is documented in Torreya, although the genome is usually inherited by the paternal donor in gymnosperms [32].
Both parents could be ancestors of the backcross, but a pattern that is often seen in instances of introgression is the asymmetric transfer of genetic material [33]. Therefore, one of the hybridizing lineages acts mainly as a donor and the other taxon as a recipient. It has been suggested that the quantity of pollen may be important and that pollen might be preferentially dispersed from a species with more pollen to one with less [34,35]. T. grandis produce a much higher quantity of pollen than T. jackii due to cultivation. Therefore, the female individuals of T. jackii (actually, its ancestor) might have been fertilized by the pollen from T. grandis to yield the F1 hybrid, before multiple backcrosses with T. grandis. Finally, the genetic material from the mother donor was diluted over the generations to yield T. jackii.
Introgressive hybridization usually occurs between the ranges of two species and it leads either to speciation by adapting to a habitat different from those of the parents or extinction via genetic assimilation [26,27]. In each case, the offspring of hybridization are only a transient phase and they are seldom recognized as independent species. However, T. jackii is widely recognized as a species with little controversy, possibly due to its distinct morphology [1,2,36]. We suggest that this species is actually quite special and different from an ordinary species because it is not genetically or geographically isolated from one of the parents, and its cp donor or maternal ancestor might have become extinct, possibly due to the cultivation of T. grandis within the range of T. jackii. Therefore, we advise caution to prevent the unreasonable expansion of economic crops into the distribution ranges of their wild relatives, thereby highlighting the complex issues involved with the protection of endangered species.

5. Conclusions

In the present study, we explored the phylogenetic relationships of Torreya species based on the chloroplast and nuclear genome data, found the plastid-nuclear incongruent phylogenetic position of endangered species T. jackii, and elucidated the introgression hybridization responsible for the incongruence due to the cultivation of sympatric species T. grandis. This finding implicates the complexity of protecting endangered species and the need for caution to prevent the unreasonable expansion of economic crops into the distribution ranges of their wild relatives.

Author Contributions

Y.-J.W. and Y.-X.K. designed the study, analyzed the data, and wrote the paper. K.X. and Y.-X.K. performed the experiments. All authors reviewed the manuscript.

Funding

National Natural Science Foundation of China: 81274024; 41461008.

Acknowledgments

We thank Jian-Quan Liu for the samples and valuable comment on the draft. This research was supported by grants from the National Natural Science Foundation of China (Grant Nos. 81274024 and 41461008).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Fu, L.K.; Li, N.; Robert, R.M. Taxaceae. In Flora of China; Hong, D.Y., Raven, P.H., Eds.; Science Press: Beijing, China, 1999; pp. 89–96. (In Chinese) [Google Scholar]
  2. Farjon, A. A Handbook of the World’s Conifers; Brill Academic Publishers: Leiden, The Netherlands, 2010. [Google Scholar]
  3. Christenhusz, M.J.M.; Reveal, J.L.; Farjon, A.; Gardner, M.F.; Mill, R.R.; Chase, M.W. New classification and linear sequence of extant gymnosperms. Phytotaxa 2011, 19, 55–70. [Google Scholar] [CrossRef]
  4. Florin, R. The distribution of conifer and taxad genera in time and space. Acta Horti Bergiani 1963, 20, 122–312. [Google Scholar]
  5. Hu, H.H. Synoptical study of Chinese Torreya. Contr. Biol. Lab. Sci. Soc. China 1927, 3, 1–9. [Google Scholar]
  6. Kang, N.; Tang, Z.X. Studies on the taxonomy of the genus Torreya. Bullet. Bot. Res. 1995, 15, 349–362. (In Chinese) [Google Scholar]
  7. Li, J.H.; Davis, C.C.; Donoghue, M.J.; Kelley, S.; Tredici, P.D. Phylogenetic relationships of Torreya (Taxaceae) inferred from sequences of nuclear ribosomal DNA ITS region. Harv. Pap. Bot. 2001, 6, 275–281. [Google Scholar]
  8. Hao, D.C.; Xiao, P.G.; Huang, B.L.; Ge, G.B.; Yang, L. Interspecific relationships and origins of Taxaceae and Cephalotaxaceae revealed by partitioned Bayesian analyses of chloroplast and nuclear DNA sequences. Plant Syst. Evol. 2008, 276, 89–104. [Google Scholar] [CrossRef]
  9. Cheng, Y.C.; Nicolson, R.G.; Tripp, K.; Chaw, S.M. Phylogeny of Taxaceae and Cephalotaxaceae genera inferred from chloroplast matK gene and nuclear rDNA ITS region. Mol. Phylogenet. Evol. 2000, 14, 353–365. [Google Scholar] [CrossRef] [PubMed]
  10. Doyle, J.J.; Doyle, J.L. A rapid DNA isolation procedure for small amounts of fresh leaf tissue. Phytochem. Bull. 1987, 19, 11–15. [Google Scholar]
  11. Grivet, D.; Heinze, B.; Vendramin, G.G.; Petit, R.J. Genome walking with consensus primers: Application to the large single copy region of chloroplast DNA. Mol. Ecol. Notes 2001, 1, 345–349. [Google Scholar] [CrossRef]
  12. Little, D.P.; Schwarzbach, A.E.; Adams, R.P.; Hsieh, C.F. The circumscription and phylogenetic relationships of Callitropsis and the newly described genus Xanthocyparis (Cupressaceae). Am. J. Bot. 2004, 91, 1872–1881. [Google Scholar] [CrossRef]
  13. Kress, W.J.; Wurdack, K.J.; Zimmer, E.A.; Weigt, L.A.; Janzen, D.H. Use of DNA barcodes to identify flowering plants. Proc. Natl. Acad. Sci. USA 2005, 102, 8369–8374. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Shaw, J.; Lickey, E.B.; Beck, J.T.; Farmer, S.B.; Liu, W.; Miller, J.; Siripun, K.S.; Winder, C.T.; Schilling, E.E.; Small, R.L. The tortoise and the hare II: Relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. Am. J. Bot. 2005, 92, 142–166. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, Y.J.; Susanna, A.; Von Raab-Straube, E.; Milne, R.; Liu, J.Q. Island-like radiation of Saussurea (Asteraceae: Cardueae) triggered by uplifts of the Qinghai-Tibetan Plateau. Biol. J. Linn. Soc. 2009, 97, 893–903. [Google Scholar] [CrossRef]
  16. Kou, Y.X.; Zeng, J.; Liu, J.Q.; Zhao, C.M. Germplasm diversity and differentiation of Helianthus tuberosus L. revealed by AFLP marker and phenotypic traits. J. Agr. Sci. 2014, 152, 779–789. [Google Scholar] [CrossRef]
  17. Thompson, J.D.; Gibson, T.J.; Plewniak, F.; Jeanmougin, F.; Higgins, D.G. The CLUSTAL_X windows interface: Flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25, 4876–4882. [Google Scholar] [CrossRef] [PubMed]
  18. Tamura, K.; Peterson, D.; Peterson, N.; Stecher, G.; Nei, M.; Kumar, S. MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol. Biol. Evol. 2011, 28, 2731–2739. [Google Scholar] [CrossRef] [PubMed]
  19. Swofford, D.L. PAUP*: Phylogenetic Analysis Using Parsimony (*and Other Methods), Version 4.0b10; Sinauer Associates: Sunderland, MA, USA, 2002. [Google Scholar]
  20. Ronquist, F.; Teslenko, M.; Van der Mark, P.; Ayres, D.L.; Darling, A.; Höhna, S.; Larget, B.; Liu, L.; Suchard, M.A.; Huelsenbeck, J.P. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst. Biol. 2012, 61, 539–542. [Google Scholar] [CrossRef] [PubMed]
  21. Posada, D.; Crandall, K.A. Modeltest: Testing the model of DNA substitution. Bioinformatics 1998, 14, 817–818. [Google Scholar] [CrossRef] [PubMed]
  22. Nei, M.; Li, W.H. Mathematical model of studying genetic variation in terms of restriction endonucleases. Proc. Natl. Acad. Sci. USA 1979, 76, 5269–5273. [Google Scholar] [CrossRef] [PubMed]
  23. Van de Peer, Y.; De Wachter, R. TREECON for Windows: A software package for the construction and drawing of evolutionary trees for the Microsoft Windows environment. Comput. Appl. Biosci. 1994, 10, 569–570. [Google Scholar] [CrossRef] [PubMed]
  24. Som, A. Causes, consequences and solutions of phylogenetic incongruence. Brief. Bioinform. 2015, 16, 536–548. [Google Scholar] [CrossRef] [PubMed]
  25. McBreen, K.; Lockhart, P.J. Reconstructing reticulate evolutionary histories of plants. Trends Plant Sci. 2006, 11, 398–404. [Google Scholar] [CrossRef] [PubMed]
  26. Abbott, R.; Albach, D.; Ansell, S.; Arntzen, J.W.; Baird, S.J.E.; Bierne, N.; Boughman, J.; Brelsford, A.; Buerkle, C.A.; Buggs, R.; et al. Hybridization and speciation. J. Evol. Biol. 2013, 26, 229–246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Todesco, M.; Pascual, M.A.; Owens, G.L.; Ostevik, K.L.; Moyers, B.T.; Hübner, S.; Heredia, S.M.; Hahn, K.L.; Caseys, C.; Bock, D.G.; et al. Hybridization and extinction. Evol. Appl. 2016, 9, 892–908. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Gao, Z.W. Study on the biological character and conservation of endemic species Torreya jackii. Chin. Biodivers. 1997, 5, 206–209. (In Chinese) [Google Scholar]
  29. Wang, X.M.; Zhou, B.Z.; Cao, Y.H.; Wang, K.; Lin, G.X. A GIS-based study on the habitat characteristics of Chinese Torreya in Kuaiji mountains. Acta Agric. Uni. Jiangxiensis 2010, 32, 523–527. (In Chinese) [Google Scholar]
  30. Dong, L.M.; Shen, D.F.; Yu, W.W.; Zeng, Y.R.; Wu, Z.M.; Dai, W.S. Trait variations of male Torreya grandis trees. J. Zhejiang A F Uni. 2012, 29, 715–721. (In Chinese) [Google Scholar]
  31. Guan, Q.; Yuan, M.; Yu, Z. Early discrimination of karyotype and sexuality for Chinese Torreya. Sci. Silvae Sinicae 1993, 29, 389–392. (In Chinese) [Google Scholar]
  32. Buchholz, J.T. The embryogeny of Torreya, with a note on Austrotaxus. Bullet. Torrey Bot. Club 1940, 67, 731–754. [Google Scholar] [CrossRef]
  33. Vallejo-Marín, M.; Cooley, A.M.; Lee, M.Y.; Folmer, M.; McKain, M.R.; Puzey, J.R. Strongly asymmetric hybridization barriers shape the origin of a new polyploid species and its hybrid ancestor. Am. J. Bot. 2016, 103, 1–17. [Google Scholar] [CrossRef]
  34. Mallet, J. Hybridization as an invasion of the genome. Trends Ecol. Evol. 2005, 20, 229–237. [Google Scholar] [CrossRef] [PubMed]
  35. Currat, M.; Ruedi, M.; Petit, R.J.; Excoffier, L. The hidden side of invasions: Massive introgression by local genes. Evolution 2008, 62, 1908–1920. [Google Scholar] [CrossRef] [PubMed]
  36. Cheng, W.C.; Fu, L.K. Torreya Arn. In Flora Reipublicae Popularis Sinica; Cheng, W.C., Fu, L.K., Eds.; Science Press: Beijing, China, 1978; pp. 457–467. (In Chinese) [Google Scholar]
Figure 1. Geographical distribution of the extant Torreya species. (A), worldwide distribution (orange regions); (B), distribution in Eastern Asia; (C), distribution in North America.
Figure 1. Geographical distribution of the extant Torreya species. (A), worldwide distribution (orange regions); (B), distribution in Eastern Asia; (C), distribution in North America.
Forests 10 00174 g001
Figure 2. Strict consensus of the most parsimonious tree obtained for Torreya based on the seven combined chloroplast regions. Numbers above the branches denote the bootstrap values based on 1000 replicates from the maximum parsimony (MP) and maximum likelihood (ML) analyses, and the posterior probabilities from Bayesian analyses. EA, Eastern Asia; NA, North America.
Figure 2. Strict consensus of the most parsimonious tree obtained for Torreya based on the seven combined chloroplast regions. Numbers above the branches denote the bootstrap values based on 1000 replicates from the maximum parsimony (MP) and maximum likelihood (ML) analyses, and the posterior probabilities from Bayesian analyses. EA, Eastern Asia; NA, North America.
Forests 10 00174 g002
Figure 3. Strict consensus of the most parsimonious tree obtained for Torreya based on the nuclear internal transcribed spacer (ITS). The other details are the same as those in Figure 2.
Figure 3. Strict consensus of the most parsimonious tree obtained for Torreya based on the nuclear internal transcribed spacer (ITS). The other details are the same as those in Figure 2.
Forests 10 00174 g003
Figure 4. Neighbor-joining tree based on the amplified fragment length polymorphism (AFLP) data. Numbers above the branches denote the bootstrap values based on 1000 replicates from MP.
Figure 4. Neighbor-joining tree based on the amplified fragment length polymorphism (AFLP) data. Numbers above the branches denote the bootstrap values based on 1000 replicates from MP.
Forests 10 00174 g004
Table 1. Species, geographic origins, vouchers, and GenBank numbers used in this study.
Table 1. Species, geographic origins, vouchers, and GenBank numbers used in this study.
TaxonGeographic Origin/Voucher No.GenBank No. (matK, rbcL, trnL-trnF, trnD-trnT, rpl16, psbB1-psbB2, trnS-trnG, ITS)
Torreya
T. jackii -1Xianju, Zhejiang/Liu and Noshiro 005KJ588955KJ589001KJ589070KJ589047KJ589024KJ588978KJ589093KJ588928
T. jackii -2Shaowu, Fujian/Kou 0104KJ588956KJ589002KJ589071KJ589048KJ589025KJ588979KJ589094KJ588929
T. jackii -3Jinyun, Zhejiang/Kou 1202KJ588957KJ589003KJ589072KJ589049KJ589026KJ588980KJ589095KJ588930-31
T. jackii -4Xianju, Zhejiang/Kou 1506KJ588958KJ589004KJ589073KJ589050KJ589027KJ588981KJ589096KJ588932-32
T. jackii -5Lishui, Zhejiang/Kou 1101KJ588959KJ589005KJ589074KJ589051KJ589028KJ588982KJ589097KJ588934-35
T. jackii -6Lishui, Zhejiang/Kou 1102KJ588960KJ589006KJ589075KJ589052KJ589029KJ588983KJ589098KJ588936-37
T. grandis -1Lin’an, Zhejiang/Liu and Noshiro 010KJ588940KJ588986KJ589055KJ589032KJ589009KJ588963KJ589078KJ588907
T. grandis -2Zhuji, Zhejiang/Liu and Noshiro 013KJ588941KJ588987KJ589056KJ589033KJ589010KJ588964KJ589079KJ588908-09
T. grandis -3Jinyun, Zhejiang/Kou 1208KJ588942KJ588988KJ589057KJ589034KJ589011KJ588965KJ589080KJ588910-11
T. fargesii -1Jinfoshan, Chongqing/Liu and Noshiro 014KJ588947KJ588993KJ589062KJ589039KJ589016KJ588970KJ589085KJ588917-18
T. fargesii -2Jinfoshan, Chongqing/Liu and Noshiro 016KJ588948KJ588994KJ589063KJ589040KJ589017KJ588971KJ589086KJ588919-20
T. fargesii -3Shenlongjia, Hubei/Liu 200807-2KJ588949KJ588995KJ589064KJ589041KJ589018KJ588972KJ589087KJ588921-22
T. fargesii -4Lijiang, Yunnan/Liu and Noshiro 018KJ588945KJ588991KJ589060KJ589037KJ589014KJ588968KJ589083KJ588914
T. fargesii -5Weixi, Yunnan/Liu and Noshiro 019KJ588946KJ588992KJ589061KJ589038KJ589015KJ588969KJ589084KJ588915-16
T. nucifera -1Forestry and Forest Products Research Institute, Japan/Nakais 200203KJ588943KJ588989KJ589058KJ589035KJ589012KJ588966KJ589081KJ588912
T. nucifera -2Kunming Bot. Gard. (from Japan)/8KJ588944KJ588990KJ589059KJ589036KJ589013KJ588967KJ589082KJ588913
T. taxifolia -1Kunming Bot. Gard./1062-89-c (Arnold Arborefum)KJ588950KJ588996KJ589065KJ589042KJ589019KJ588973KJ589088KJ588923
T. taxifolia -2Kunming Bot. Gard./857-87*aKJ588951KJ588997KJ589066KJ589043KJ589020KJ588974KJ589089KJ588924
T. taxifolia -3Kunming Bot. Gard./1054-89-IKJ588952KJ588998KJ589067KJ589044KJ589021KJ588975KJ589090KJ588925
T. californica -1Philipps-Universität Marburg Bot. Gard./7KJ588953KJ588999KJ589068KJ589045KJ589022KJ588976KJ589091KJ588926
T. californica -2San Francisco, California/Bruce Bartholomew, s.n. SFBG# XY-226KJ588954KJ589000KJ589069KJ589046KJ589023KJ588977KJ589092KJ588927
Amentotaxus
A. argotaeniaJinfoshan, Chongqing/3025KJ588962KJ589008KJ589077KJ589054KJ589031KJ588985KJ589100KJ588939
A. yunnanensisHeilongtan, Kunming/3007KJ588961KJ589007KJ589076KJ589053KJ589030KJ588984KJ589099KJ588938
Table 2. Statistic for the sequences used in this study. PI: parsimony-informative sites; CI: consistency index; RI: retention index.
Table 2. Statistic for the sequences used in this study. PI: parsimony-informative sites; CI: consistency index; RI: retention index.
LociLengthVariable SitesPICIRI
ITS128152460.8390.960
matK143414131.0001.000
psbB1-psbB21434431.0001.000
rbcL129912110.9440.969
rpl1696213111.0001.000
trnD-trnT151921140.9800.986
trnL-trnF91415130.9840.987
trnS-trnG679840.9050.918
Table 3. Pairwise distances (%) among ITS (lower left) and combined plastid (upper right) sequences from six Torreya species.
Table 3. Pairwise distances (%) among ITS (lower left) and combined plastid (upper right) sequences from six Torreya species.
Plastid/ITST. fargesiiT. nuciferaT. taxifoliaT. californicaT. grandisT. jackii
T. fargesii 0.200.420.400.270.33
T. nucifera0.92 0.280.260.120.21
T. taxifolia2.292.10 0.390.370.39
T. californica1.901.720.68 0.350.36
T. grandis1.350.742.542.16 0.33
T. jackii1.200.602.392.010.57

Share and Cite

MDPI and ACS Style

Wang, Y.-J.; Xiao, K.; Kou, Y.-X. Torreya jackii (Taxaceae): A Special Species that is Genetically Admixed, Morphologically Distinct, and Geographically Sympatric with Parent Species. Forests 2019, 10, 174. https://doi.org/10.3390/f10020174

AMA Style

Wang Y-J, Xiao K, Kou Y-X. Torreya jackii (Taxaceae): A Special Species that is Genetically Admixed, Morphologically Distinct, and Geographically Sympatric with Parent Species. Forests. 2019; 10(2):174. https://doi.org/10.3390/f10020174

Chicago/Turabian Style

Wang, Yu-Jin, Kun Xiao, and Yi-Xuan Kou. 2019. "Torreya jackii (Taxaceae): A Special Species that is Genetically Admixed, Morphologically Distinct, and Geographically Sympatric with Parent Species" Forests 10, no. 2: 174. https://doi.org/10.3390/f10020174

APA Style

Wang, Y. -J., Xiao, K., & Kou, Y. -X. (2019). Torreya jackii (Taxaceae): A Special Species that is Genetically Admixed, Morphologically Distinct, and Geographically Sympatric with Parent Species. Forests, 10(2), 174. https://doi.org/10.3390/f10020174

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