Phylogeography and allopatric divergence of cypress species (Cupressus L.) in the Qinghai-Tibetan Plateau and adjacent regions
© Xu et al; licensee BioMed Central Ltd. 2010
Received: 5 August 2009
Accepted: 22 June 2010
Published: 22 June 2010
Although allopatric speciation is viewed as the most common way in which species originate, allopatric divergence among a group of closely related species has rarely been examined at the population level through phylogeographic analysis. Here we report such a case study on eight putative cypress (Cupressus) species, which each have a mainly allopatric distribution in the Qinghai-Tibetan Plateau (QTP) and adjacent regions. The analysis involved sequencing three plastid DNA fragments (trnD-trnT, trnS-trnG and trnL-trnF) in 371 individuals sampled from populations at 66 localities.
Both phylogenetic and network analyses showed that most DNA haplotypes recovered or haplotype-clustered lineages resolved were largely species-specific. Across all species, significant phylogeographic structure (NST > GST, P < 0.05) implied a high correlation between haplotypes/lineages and geographic distribution. Two species, C. duclouxiana and C. chengiana, which are distributed in the eastern QTP region, contained more haplotypes and higher diversity than five species with restricted distributions in the western highlands of the QTP. The remaining species, C. funebris, is widely cultivated and contained very little cpDNA diversity.
It is concluded that the formation of high mountain barriers separating deep valleys in the QTP and adjacent regions caused by various uplifts of the plateau since the early Miocene most likely promoted allopatric divergence in Cupressus by restricting gene flow and fixing local, species-specific haplotypes in geographically isolated populations. The low levels of intraspecific diversity present in most species might stem from population bottlenecks brought about by recurrent periods of unfavorable climate and more recently by the negative impacts of human activities on species' distributions. Our findings shed new light on the importance of geographical isolation caused by the uplift of the QTP on the development of high plant species diversity in the QTP biodiversity hotspot.
Geographic isolation between populations within species is a major precursor of allopatric speciation [1, 2]. In the absence of gene flow, different alleles become fixed in different populations through the action of selection and/or genetic drift potentially leading to the origin of new species . It is expected, therefore, that in regions where geographic isolation is common, for example, in areas where high mountains separate populations occurring in deep valleys, allopatric divergence will be promoted and accelerated due to the formation of locally adapted gene pools and possibly the origin of intrinsic postzygotic isolation according to the Bateson-Dobzhansky-Muller model [3–5]. A phylogeographic analysis of closely related species that occur in such areas can be instructive on the importance of geographic isolation and allopatric divergence in their origin and subsequent evolution .
We used plastid (mainly chloroplast) DNA sequence data to examine interspecific divergence and to assess intraspecific diversity of these species. Plastid DNA is paternally inherited and dispersed mainly by pollen in all species of the Cupressaceae investigated [20–23]. Recent studies suggest that paternally inherited plastid DNA sequences in conifers may be species-specific and therefore effective in tracing lineage sorting between species [24–27]. Plastid DNA sequence variation has been used previously in phylogenetic studies of cypresses, but only single sample from most species were examined . In addition, surveys of RAPD variation  and nuclear ITS variation  have been conducted within and among cypress species including those examined in the current study. However, the results of these studies were inconsistent with regard to interspecific relationships . In this study, we aimed to address the following questions based on surveys of sequence variation for three plastid DNA fragments within and between populations at the species level: (1) Do plastid DNA haplotypes or clusters of such haplotypes correspond with morphological differentiation and the allopatric distributions of the Cupressus species examined? (2) Does intraspecific DNA diversity correlate with the natural population size of a species and, in turn, provide guidance for conservation management?
Haplotype distribution and relationships
Because the number of mutations distinguishing any two lineages was always greater than six, we constructed separate haplotype networks for the three lineages that contained most haplotypes (i.e. lineages A, C and E). Within lineage E, two haplotypes H10 and H3 occurred in a small sub-lineage that was separated by four mutations from other haplotypes (Fig. 2). However, H10 co-occurred with the most frequent haplotype in this lineage, H4, in most sampled populations. The single haplotype (H5) recovered from C. jiangeensis was also present in lineage E and is derived from H4. In lineage A, H17 was the most distantly related haplotype relative to other haplotypes (Fig. 2).
Diversity among and within species
Estimates of average genetic diversity within populations (HS), total gene diversity (HT), interpopulation haplotype differentiation (GST), and interpopulation haplotype differentiation taking into account sequence difference (NST) (mean ± SE in parentheses)
NST - GST
0.558 (0.056) *
C. chengiana (lineage E,
-0.007 (0.075) ns
C. chengiana (lineage A,
-0.021 (0.055) ns
C. cashmeriana &
C. torulosa s.l. & C. gigantea
Analysis of molecular variance (AMOVA) for cpDNA haplotypes of cypress species
Source of variation
C. chengiana (lineage E,
FST = 0.068
C. chengiana (lineage A,
FST = -0.002
FCT = 0.439*
Among populations within groups
FSC = 0.032
FST = 0. 457*
FST = -0.035
FST = 0.233*
FST = 0
FST = 1.000
C. cashmeriana &
FST = 0
C. torulosa s.l. and
FCT = 0.976*
Among populations within groups
FSC = 1.000
FST = 1.000*
FCT = 0.635*
Among populations within species
FSC = 0.210 *
FST = 0.712 *
Discussion and Conclusions
Two broad findings emerge from our population-level phygeographic study of the eight Cupressus species that occur in the QTP and adjacent regions. First, there was generally good correspondence between the major plastid DNA lineages/haplotypes identified and the morphological taxonomy of species in the group . Such divergence in plastid DNA sequence between species in this group probably results from their allopatric origins and the maintenance of broadly allopatric distributions following extensive uplifts of the QTP since the Miocene. Second, the amounts and pattern of intraspecific plastid DNA variation recorded differed greatly among the species examined. Species with restricted natural distributions contained low levels of intraspecific diversity relative to those having wide natural distributions. These two findings are important to an improved understanding of the importance of geographical isolation on the development of high plant species diversity in the QTP biodiversity hotspot and also the distribution of intraspecific diversity within this region.
Phylogenetic relationships among the major plastid DNA lineages detected were largely consistent with those resolved for species in a previous plastid DNA phylogeny constructed by Little  based on other plastid DNA sequences and mainly single sample per species. Interestingly, in both phylogenies the Mediterranean C. sempervirens and C. atlantica were nested within the Asian lineage indicating that they probably originated from this lineage. Phylogenetic relationships resolved in both plastid DNA phylogenies were a little different from those constructed from RAPD or nuclear ITS sequence data [12, 28, 29]. Because these latter studies sampled only one or very few individuals from each species, it remains unknown whether the nuclear variation reported was species-specific or not. In contrast to the other studies, our investigation is the only one to have examined molecular variation in cypresses at the population level. Our results suggest that two species (C. austrotibetica and C. cashmeriana) share the same haplotype, which is very similar to the single haplotype detected in C. torulosa. However, all three of these species differ in morphology . The fact that two very distinct cpDNA lineages were recovered from northern and southern populations of C. chengiana, is of further interest in that plants from these two regions have sometimes been treated as separate subspecies [13, 14]. Thus, our cpDNA data support this division; however, morphological as well as further molecular studies are necessary before a taxonomic revision of this and other Asian cypress species is contemplated.
Our data show that plastid DNA lineage sorting is nearly complete between recognized morphological units of Asian Cupressus (Fig. 2). Furthermore, across all populations and species examined, a strong signal for phylogeographic structure is present: N ST (0.962) was significantly larger than G ST (0.791, P < 0.05) (Table 1), demonstrating that gene-flow is very low between these broadly allopatric species (Fig. 1). According to preliminary calibrations, C. funebris diverged from the other lineage/species around 26 million years ago (Mya), while the two major east-west divergences between C. duclouxiana and C. gigantea, and between C. chengiana and the C. torulosa/C. austrotibetica/C. cashmeriana lineage, were dated to around 10 Mya and 13 Mya, respectively. Although these estimates of dates of origin should be treated with caution, they correspond well with geological evidence that the QTP was extensively uplifted from the early Miocene to the Pliocene [30–34]. It is likely that these extensive uplifts created fragmentation and isolation of habitats, which promoted repeated allopatric speciation of Cupressus in the QTP and adjacent regions. The formation of high mountains separated by deep valleys will have reduced gene flow between Cupressus populations that were diverging into different species in different valleys and hastened the fixation of locally specific haplotypes. However, allopatric differentiation usually results in the origin of 'paraphyletic' species because of incomplete lineage sorting . Thus, the formation of almost completely mutual 'monophyletic lineages' correlated with morphological units, which we recorded here, is likely to have been brought about by repeated extinctions and/or genetic drift in the isolated habitats [3–5]. A hypothesis of the allopatric origin of each Cupressus species in a different valley or valley system is supported by the low level of inter-population differentiation and minimal population substructure observed in most species/lineages except C. duclouxiana and C. chengiana (Table 2). In these latter two species, a high level of population differentiation was recorded among geographically isolated regions, but very little population divergence within such regions (Fig. 1, Table 2). Two lineages of C. chengiana may partly reflect a long history of divergence with the species fragmented into two separate areas in the north and south of western China approximately nine Mya.
The presence of strong geographical isolating barriers in the form of high mountains might also have contributed to rapid species diversification in other plant groups present in the QTP [9, 35–38]. For example, in the genus Ligularia, most species with distinct morphologies and restricted distributions are thought to have originated through allopatric isolation from a common ancestor within a very short time correlated with the period of extensive uplift of the QTP . Our phylogeographic analysis at the population level of closely related cypress species in the QTP, thus provides further evidence for geographic isolation acting as a driving force of high rates of speciation in this and adjacent regions.
Whereas high levels of plastid DNA diversity were recorded in C. chengiana and C. duclouxiana, only one or two haplotypes were present in each of C. gigantea, C. austrotibetica, C. cashmeriana, C. torulosa and C. funebris. This difference in intraspecific diversity is highly correlated with the natural distribution ranges of the sampled species [10–19]. Cupressus gigantea, C. austrotibetica, C. cashmeriana and C. torulosa occur sparsely at high altitude in the QTP and the west Himalayas, which have been highly vulnerable to the effects of global climatic change . These species might have been more widespread during periods of more favorable climate in their history. Alternatively or in addition, their restricted distributions might be a consequence of the negative impacts of human activities. Either way, they may have experienced severe population bottlenecks in the relatively recent past which reduced significantly their genetic diversity. In contrast to the other four species, Cupressus funebris is distributed in central and eastern China and adjacent regions where adverse climatic effects should have been less intense. It seems, however, that most natural populations of this species have become extinct in the wild and only a few, comprising trees of large size, survived in south China and north Vietnam [10–19]. However, according to our own field observations, even these populations (for example, the Ruyuan population in Guangdong) (Additional File 1) are likely to have been artificially cultivated because these large trees grow around one ancient temple. The presence of only two haplotypes in this species may be a legacy of its widespread cultivation and eradication from the wild as seems also to be the case for C. sempervirens in Europe . In contrast to these Asian cypress species showing low levels of plastid DNA diversity, C. duclouxiana, and C. chengiana contain high levels of diversity and occur at lower altitude at the eastern margin of the QTP, where they might have experienced fewer bottlenecks and/or founder effects caused by historical climatic change . It is also to be noted that most extant populations of C. chengiana are natural, as are populations of C. duclouxiana except those cultivated in North Yunnan, and therefore are likely to have been less affected by human activities and subsequent bottleneck effects.
All of the species examined, except for C. duclouxiana and C. funebris, are listed as endangered . Genetic diversity is critically important for a species' evolutionary potential in the face of environmental change, and loss of genetic diversity is often associated with reduced fitness . The amounts and distribution of genetic diversity within these cypress species are thus important in regard to formulating appropriate plans for conservation purposes. Based on the amounts and distribution of plastid DNA sequence variation within these species, it would seem that for most species, conservation of only a few populations will be sufficient to maintain plastid DNA diversity at current levels. However, the high population differentiation found in C. duclouxiana (Table 1) suggests a need to conserve more populations for maintenance of plastid DNA diversity in this species. Similarly, in C. chengiana our results make clear that conservation of multiple populations in both Gansu and Sichuan provinces will be necessary to ensure that plastid DNA diversity is maintained in this species.
Species and samples
Leaf samples were collected from eight Asian Cupressus species that were identified according to the recent taxonomic revision of Little : C. chengiana S. Y. Hu, C. jiangeensis N. Zhao, C. cashmeriana Royle ex Carrière, C. duclouxiana Hickel ex Camus, C. funebris Endl., C. torulosa D. Don, C. austrotibetica Silba and C. gigantea W. C. Cheng & L. K. Fu. The one Asian species not included was C. tonkinensis Silba, which occurs in Vietnam and across the Vietnam-China border. In total, samples were collected from 371 trees over 66 populations (Additional File 1), which covered most of the natural or cultivated distribution of the species examined. All individuals of C. austrotibetica sampled appeared to have been artificially cultivated because they grew around ancient temples, while three populations of C. duclouxiana (29, 30 and 31) entirely comprised cultivated trees. Most populations of C. chengiana and C. gigantea were natural, although a few individuals in the Lixian seemed to have been cultivated. Only one individual of C. jiangeensis was available for analysis. When more than one sample was collected from a population, samples were taken from trees at least 50 m apart. The latitude, longitude and altitude of localities for most populations sampled were recorded using an Etrex GIS monitor (Garmin, Taiwan). Also surveyed were three cultivated individuals of C. sempervirens sampled from England, Italy and China, respectively, and one individual of C. atlantica. We sequenced three individuals of the New World cypress species as outgroups because recent phylogenetic analyses suggest that all Old World cypress species comprise a monophyletic lineage that is distinct from the New World cypress lineage and juniper lineage .
DNA extraction, amplification and sequencing
Genomic DNA was isolated from approximately 50 mg of silica-gel dried, leaf-needle material using the CTAB method . Three to five individuals from each of 25 populations sampled across most species were used in an initial screen of plastid DNA polymorphisms. Three different plastid DNA primer pairs were selected for use in this scan, and all revealed polymorphism. Subsequently, we used these three primer pairs to survey all individuals sampled for the trnD-trnT , trnS-trnG  and trnL-trnF  regions, respectively. Polymerase chain reaction (PCR) was performed following the procedure of Zhang et al.  in a 25-μL volume, containing 10-40 ng plant DNA, 50 mM Tris-HCl (pH 8.0), 1.5 mM MgCl2, 0.5 mM dNTPs, 2 μM of each primer, and 0.75 unit of Taq polymerase.
PCR products were purified using a TIANquick Midi Purification Kit following the recommended protocol (TIANGEN). Sequencing reactions were performed with the PCR primers described above, to cover the whole PCR segment using ABI Prism BigdyeTM Terminator Cycle Sequencing Ready Reaction Kit. Sequencing was performed with BigDye Terminator 3.1 (Applied Biosystems) on an ABI PRISM 3130×l sequencer using the same primers as used for PCR amplifications. Sequences were aligned with CLUSTAL X  and double-checked manually. A matrix of combined sequences was constructed for the 371 individuals examined, and showed that 30 different cpDNA sequences (haplotypes) were present. All sequences have been deposited in the GenBank databases under accession numbers HM008329-HM008379.
We used PAUP 4.0  to perform maximum parsimony and neighbour joining analyses to examine phylogenetic relationships among cpDNA haplotypes (Additional File 3). All indels were coded as binary characters (0 or 1) using the GAPCODER program  and those with further indels were sequentially assumed as multiple evolutionary events. Maximum parsimony analyses (equally weighted characters and nucleotide transformations) involved a heuristic search strategy with 100 replicates of random addition of sequences, in combination with ACCTRAN character optimization and MULPARS+TBR branch-swapping and STEEPEST DESCENT options on. Bootstrap values for MP analyses (1000 replicates) were used to assess the relative support for monophyletic groups .
We used BEAST (under the common model GTR)  to estimate the mutation rate of Cupressus and the possible divergence timescales of the main cpDNA lineages across the eight Asian cypress species. We edited two cpDNA (trnD-trnT, trnS-trnG and trnL-trnF) datasets respectively. The first comprised three Cupressus species, three Juniperus species, three Calitropsis species, one species for Platycladus, Calocedrus, Chamaecyparis and Thuja. The other included the most common haplotypes of the Cupressus species and three Calitropsis species. For analyses of both datasets, we excluded all indels and only nucleotide mutations were considered. A strict molecular clock assumption was rejected for both datasets because of significant rate heterogeneity among lineages. We used BEAST version 1.5.3 employing a Bayesian Markov chain Monte Carlo (MCMC) chain and co-estimated topology to estimate the possible mutation rate of Cupressus and divergence timescales across the Asian cypress species. Under the GTR model of nucleotide substitution with a gamma distribution and four rate categories, three tree prior models (Coalescent: constant size, exponential growth; and Speciation: Yule process) were implemented with rate variation across branches assumed to be uncorrelated and log-normally distributed. The model that yielded the highest posterior probability estimates was employed to perform a final analysis. For all analyses, posterior distributions of parameters were approximated using two independent MCMC analyses of 100,000,000 generations with 20% burn-in. The program Tracer v1.5  was used to check Effective Sample Size (ESS) and the program TreeAnnotator 1.5.3 (part of the BEAST 1.5.3 package) was used to combine all samples and to converge and/or summarize output results. A tree with ages for each node and their 95% HPD were displayed and modified in FigTree v1.3.1 . According to the fossil calibrations of the first dataset, a mutation rate of cpDNA (3.2 × 10-10 substitutions per site per year) was obtained (Additional File 4). This mutation rate was used to calibrate the interspecific or interlineage divergences across the cypress species (Additional File 4).
An examination of interrelationships between haplotypes within lineages recovered in phylogenetic trees was conducted using TCS 1.21 . TCS was run with the default parsimony connection limit of 95%. We assumed that all coded indels and mutations evolved at an equal rate. Estimates of unbiased genetic diversity (HE) [equivalent to expected heterozygosity for diploid data] were calculated as ś: ś = n/(n-1) (1- p i 2), in which n is the number of samples in the population, k is the number of haplotypes, and p is the population frequency of the ith haplotype . Where there were insufficient numbers of individuals in a population, multiple nearby populations were combined into a larger one for analysis (for example, populations 43-46, each comprising a single tree, were treated together as one population). Average gene diversity within populations (HS), total gene diversity (HT) and two measures of population differentiation, GST (coefficient of genetic variation over all populations ) and NST (equivalent coefficient taking into account sequence similarities between haplotypes) were calculated using the program PERMUT (available at http://www.pierroton.inra.fr/genetics/labo/Software/) . Significant phylogeographic structure was inferred by testing whether NST was significantly greater than GST using PERMUT with 1000 permutations. Genetic structure was further examined by analysis of molecular variance (AMOVA, ) as implemented in Arlequin version 3.01 . The significance of isolation by distance between populations was tested using a Mantel test with 1000 random permutations on matrices of pairwise population FST values and the natural logarithm of geographical distances . Pairwise FST values between populations were estimated using Arlequin while geographical distances between populations were calculated with the help of the program at http://www.indo.com/distance/ (e.g. ).
This research was supported by grants from the National Natural Science Foundation of China (30725004, 30430560) and Key Innovation Project of Ministry of Education of China to JQL, and a Royal Society-NSA China International Joint Project award 2010/R4 to RJA and JQL.
- Mayr E: Systematics and the Origin of Species. 1942, New York: Columbia University Press
- Sobel JM, Chen GF, Watt LR, Schemske DW: The biology of speciation. Evolution. 2010, 64: 295-315. 10.1111/j.1558-5646.2009.00877.x.View ArticlePubMed
- Coyne JA, Orr HA: Speciation. 2004, Massachusetts: Sinauer Associates Inc
- Orr HA: Dobzhansky, Bateson, and the genetics of speciation. Genetics. 1996, 144: 1331-1335.PubMed CentralPubMed
- Jaramillo-Correa JP, Aguirre-Planter E, Khasa DP, Eguiarte LE, Piñero D, Furnier GR, Bousquet J: Ancestry and divergence of subtropical montane forest isolates: molecular biogeography of the genus Abies (Pinaceae) in southern México and Guatemala. Mol Ecol. 2008, 17: 2476-2490. 10.1111/j.1365-294X.2008.03762.x.View ArticlePubMed
- Avise JC: Phylogeography: The History and Formation of Species. 2000, Cambridge, Massachusetts: Harvard University Press
- Myers N, Mittermeier RA, Mittermeiri CG, da Fonseca GAB, Kent J: Biodiversity hotspots for conservation priorities. Nature. 2000, 403: 853-858. 10.1038/35002501.View ArticlePubMed
- Yang FS, Li YF, Ding X, Wang XQ: Extensive population expansion of Pedicularis longiflora (Orobanchaceae) on the Qinghai-Tibetan Plateau and its correlation with the Quaternary climate change. Mol Ecol. 2008, 17: 5135-5145. 10.1111/j.1365-294X.2008.03976.x.View ArticlePubMed
- Wang LY, Abbott RJ, Zheng W, Chen P, Wang YJ, Liu JQ: History and evolution of alpine plants endemic to the Qinghai-Tibetan Plateau: Aconitum gymnandrum (Ranunculaceae). Mol Ecol. 2009, 18: 709-721. 10.1111/j.1365-294X.2008.04055.x.View ArticlePubMed
- Farjon A: A Monograph of Cupressaceae and Sciadopitys. Kew: Royal Botanic Gardens. 2005
- Fu LG, Jin JM: China Plant Red Data Book--Rare and Endangered Plants. 1992, Beijing, China: Science Press
- Little DP: Evolution and circumscription of the True Cypresses (Cupressaceae: Cupressus). Syst Bot. 2006, 31: 461-480.View Article
- Silba J: The Trans-Pacific relationship of Cupressus in India and North America. J Int Conifer Preserv Soc. 1994, 1: 1-28.
- Silba J: A monograph of the genus Cupressus L. J Int Conifer Preserv Soc. 1998, 5: 1-100.
- Silba J: A monograph of the genus Cupressus L. in the twenty-first century. J Int Conifer Preserv Soc. 2005, 12: 30-103.
- Fu LG, Yu YF, Farjon A: Cupressaceae. Flora of China. 1999, St. Louis, USA and Beijing, China: Science Press and Missouri Botanical Garden Press, 65-67.
- Zhao N: Species nova generis Cupressi. Acta Phytotax Sin. 1980, 18: 210-
- Silba J: Revised generic concepts of Cupressus L. (Cupressaceae). Phytologia. 1981, 49: 394-
- Silba J: A new species of Cupressus L. from Tibet (Cupressaceae). Phytologia. 1988, 65: 333-336.
- Neale DB, Marshall KA, Harry DE: Inheritance of chloroplast and mitochondrial DNA in incense-cedar (Calocedrus decurrens). Can J For Res. 1991, 21: 717-720. 10.1139/x91-100.View Article
- Mogensen HL: The hows and whys of cytoplasmic inheritance in seed plants. Am J Bot. 1996, 83: 383-404. 10.2307/2446172.View Article
- Kondo T, Sumura Y, Kawahara T, Okamura M: Inheritance of chloroplast and mitochondrial DNA in interspecific hybrids of Chamaecyparis spp. Breed Sci. 1998, 78: 177-179.
- Ennos RA, Sinclair WT, Hu X-S, Langdon A: Using organelle markers to elucidate the history, ecology and evolution of plant populations. Molecular Sytematics and Plant Evolution. Edited by: Hollingsworth PM, Bateman RM, Gornall RJ. 1999, London: Taylor and Francis, 1-19.View Article
- Petit RJ, Duminil J, Fineschi S, Hampe A, Salvini D, Vendramin GG: Comparative organization of chloroplast, mitochondrial and nuclear diversity in plant populations. Mol Ecol. 2005, 14: 689-701. 10.1111/j.1365-294X.2004.02410.x.View ArticlePubMed
- Currat M, Ruedi M, Petit RJ, Excoffier L: The hidden side of invasions: massive introgression by local genes. Evolution. 2008, 62: 1908-1920.PubMed
- Du KF, Petit RJ, Liu JQ: More introgression with less gene flow: Chloroplast vs. Mitochondrial DNA in the Picea asperata complex in China, and comparison with other conifers. Mol Ecol. 2009, 18: 1396-1407. 10.1111/j.1365-294X.2009.04107.x.View ArticlePubMed
- Little DP, Stevenson D: A comparison of algorithms for the identification of specimens using DNA barcodes: examples from gymnosperms. Cladistics. 2007, 23: 1-21. 10.1111/j.1096-0031.2006.00126.x.View Article
- Rushforth K, Adams RP, Zhong M, XQiang M, Pandley RN: Variation among Cupressus species from the eastern hemisphere based on Random Amplified Polymorphic DNAs (RAPDs). Biochem Syst Ecol. 2003, 31: 17-24. 10.1016/S0305-1978(02)00073-X.View Article
- Xiang Q, Li J: Derivation of Xanthocyparis and Juniperus from within Cupressus: evidence from sequences of nrDNA internal transcribed spacer region. Harvard Pap Bot. 2005, 9: 375-382.
- Shi YF, Li JJ, Li BY: Uplift and Environmental Changes of Qinghai-Tibetan Plateau in the Late Cenozoic. 1998, Guangzhou, China: Guangdong Science and Technology Press
- Harrison TM, Copeland P, Kidd WSF, Yin A: Raising Tibet. Science. 1992, 255: 1663-1670. 10.1126/science.255.5052.1663.View ArticlePubMed
- Spicer R, Harris N, Widdowson M, Herman A, Guo S, Valdes P: Constant elevation of southern Tibet over the past 15 million years. Nature. 2003, 421: 622-624. 10.1038/nature01356.View ArticlePubMed
- Guo ZT, Ruddiman WF, Hao QZ, Wu HB, Qiao YS, Zhu RX, Peng SZ, Wei JJ, Yuan BY, Liu TS: Onset of Asian desertification by 22 Myr ago inferred from loess deposits in China. Nature. 2002, 416: 159-163. 10.1038/416159a.View ArticlePubMed
- An ZS, Kutzbach JE, Prell WL, Port SC: Evolution of Asian monsoons and phased uplift of the Himalayan-Tibetan Plateau since late Miocene times. Nature. 2001, 411: 62-66. 10.1038/35075035.View Article
- Wang LY, Ikeda H, Liu TL, Wang YJ, Liu JQ: Repeated range Expansion and Glacial Endurance of Potentilla glabra (Rosaceae) in the Qinghai-Tibetan Plateau. J Integr Plant Biol. 2009, 51: 698-706. 10.1111/j.1744-7909.2009.00818.x.View ArticlePubMed
- Liu JQ, Wang YJ, Wang AL, Ohba H, Abbott R: Radiation and diversification of Ligularia-Cremanthodium-Parasenecio triggered by uplifts of the Qinghai-Tibetan Plateau. Mol Phylogenet Evol. 2006, 38: 31-49. 10.1016/j.ympev.2005.09.010.View ArticlePubMed
- Wang YJ, Susanna A, von Raab-Straube E, Milne R, Liu JQ: Island-like radiation of Saussurea (Asteraceae: Cardueae) trigged by uplifts of the Qinghai-Tibetan Plateau. Biol J Linn Soc. 2009, 97: 893-903. 10.1111/j.1095-8312.2009.01225.x.View Article
- Liu JQ, Gao TG, Chen ZD, Lu AM: Molecular phylogeny and biogeography of the Qinghai-Tibet Plateau endemic Nannoglottis (Asteraceae). Mol Phylogenet Evol. 2002, 23: 307-325. 10.1016/S1055-7903(02)00039-8.View ArticlePubMed
- Bagnoli F, Vendramin GG, Buonamici A, Doulis AG, Gonzálezmartínez SC, Laporta N, Magri D, Raddi P, Sebastiani F, Fineschi S: Is Cupressus sempervirens native in Italy? An answer from genetic and palaeobotanical data. Mol Ecol. 2009, 18: 2276-2286. 10.1111/j.1365-294X.2009.04182.x.View ArticlePubMed
- Guerrant EO: Genetic and demographic considerations in the sampling and reintroduction of rare plants. Conservation Biology: The Theory and Practice of Nature Conservation, Preservation and Management. Edited by: Fiedler PL, Jain SK. New York: Chapman & Hall, 1992: 321-344.
- Mao KS, Hao G, Liu JQ, Adams RP, Milne RI: Diversification and biogeography of Juniperus (Cupressaceae): variable diversification rates and multiple intercontinental dispersals. New Phytol. 2010
- Doyle JJ, Doyle JL: A rapid DNA isolation procedure for small quantities of fresh leaf material. Phytochem Bull. 1987, 19: 11-15.
- Demesure B, Sodzi N, Petit RJ: A set of universal primers for amplification of polymorphic non-coding regions of mitochondrial and chloroplast DNA in plants. Mol Ecol. 1995, 4: 129-131. 10.1111/j.1365-294X.1995.tb00201.x.View ArticlePubMed
- Hamilton MB: Four primer pairs for the amplification of chloroplast intergenic regions with intraspecific variation. Mol Ecol. 1999, 8: 521-523.PubMed
- 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-1109. 10.1007/BF00037152.View ArticlePubMed
- Zhang Q, Chiang TY, George M, Liu JQ, Abbott RJ: Phylogeography of the Qinghai-Tibetan Plateau endemic Juniperus przewalskii (Cupressaceae) inferred from chloroplast DNA sequence variation. Mol Ecol. 2005, 14: 3513-3524. 10.1111/j.1365-294X.2005.02677.x.View ArticlePubMed
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMed
- Swofford DL: PAUP* Phylogenetic Analysis Using Parsimony (* and Other Methods), Version 4. 2000, Sunderland, Massachusetts: Massachusetts Sinauer Associates
- Young ND, Healy J: GapCoder automates the use of indel characters in phylogenetic analysis. BMC Bioinformatics. 2003, 4: 1-6. 10.1186/1471-2105-4-6.View Article
- Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.View Article
- Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.PubMed CentralView ArticlePubMed
- Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1660. 10.1046/j.1365-294x.2000.01020.x.View ArticlePubMed
- Nei M: Molecular Evolutionary Genetics. 1987, New York: Columbia University Press
- Nei M: Analysis of gene diversity in subdivided populations. PNAS. 1973, 70: 3321-3323. 10.1073/pnas.70.12.3321.PubMed CentralView ArticlePubMed
- Pons O, Petit RJ: Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics. 1996, 144: 1237-1245.PubMed CentralPubMed
- Excoffier L, Smouse PE, Quattro JM: Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics. 1992, 131: 479-491.PubMed CentralPubMed
- Excoffier L, Laval G, Schneider S: ARLEQUIN (version 3.01): an integrated software package for population genetic data analysis. Evol Bioinform Online. 2005, 1: 47-50.PubMed Central
- Rousset F: Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics. 1997, 145: 1219-1228.PubMed CentralPubMed
- Chen KM, Abbott RJ, Milne RI, Tian XM, Liu JQ: Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China. Mol Ecol. 2008, 17: 4276-4288. 10.1111/j.1365-294X.2008.03911.x.View ArticlePubMed
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.