Evaluating species nonmonophyly as a trait affecting genetic diversity: a case study of three endangered species of Antirrhinum L. (Scrophulariaceae)

Molecular markers are routinely used to assess levels of diversity within and among populations, particularly with regard to species of conservation concern. However, when interpreting the level and partitioning of diversity observed, an implicit assumption is often made that the populations of the species in question form a monophyletic group. We tested this assumption in three endemics of Antirrhinum (A. charidemi, A. subbaeticum, and A. valentinum) using 79 nuclear [internal transcribed spacer (ITS)] and 85 plastid (psbA-trnH, trnT-trnL, trnK-matK, trnS-trnG) sequences representing multiple accessions of each of 24 Antirrhinum species (single accession of A. cirrhigerum). These species share six life history traits implicated in levels of genetic diversity, and have been the subject of previous population genetic studies. Populations of all three species formed monophyletic groups on ITS analysis. In contrast, none of the three species formed monophyletic groups on plastid sequence analysis: populations of A. charidemi fall in a monophyletic group including one accession of A. mollissimum, populations of A. subbaeticum form a polyphyletic group with plastid sequences shared with A. pulverulentum, and populations of A. valentinum are unresolved within a clade containing six other species. Lack of monophyly using plastid sequences is interpreted as a combination of shared ancestral polymorphism and hybridization in a reticulate evolutionary history of these species. Monophyly in the ITS tree may reflect a more recent sequence homogenization. We draw attention to the evaluation of species monophyly alongside the contribution of other life history traits in the historical interpretation of the level and partitioning of genetic diversity, and its use in recommendations for species conservation programs.     

A molecular phylogeny of the genus Gagea (Liliaceae) in Germany inferred from non-coding chloroplast and nuclear DNA sequences

The systematics of the mainly yellow flowered Gagea species complex (Liliaceae) has long been considered difficult because only a few phenotypic features within this genus and as a result of hypothesized interspecific hybridisation. A molecular phylogenetic study of seven Gagea species (G. bohemica, G. lutea, G. minima, G. pomeranica, G. pratenis, G. spathacea and G. villosa) from Germany has been undertaken, based on plastid DNA sequences (trnL(UAA)-trnF(GAA), psbA-trnH) and on the nuclear ribosomal internal transcribed spacer (ITS). Sequence divergence within the Gagea species ranges up to 15.5% for psbA-trnH, 22.0% for trnL-trnF and 23.7% for ITS (ITS1 + 5.8S rRNA + ITS2). Two subspecies of Gagea bohemica: G. bohemica subsp. saxatilis and G. bohemica subsp. bohemica are characterized by trnL-trnF data and morphological features. Analysis of the ITS region shows that G. pomeranica represents a hybrid of G. pratensis and G. lutea. Lloydia serotina was initially used as an outgroup species, but it was placed within the investigated Gagea species in the psbA-trnH and the trnL-trnF phylogenetic tree.     

Potential use of DNA barcoding for the identification of Salvia based on cpDNA and nrDNA sequences

An effective DNA marker for authenticating the genus Salvia was screened using seven DNA regions (rbcL, matK, trnL–F, and psbA–trnH from the chloroplast genome, and ITS, ITS1, and ITS2 from the nuclear genome) and three combinations (rbcL + matK, psbA–trnH + ITS1, and trnL–F + ITS1). The present study collected 232 sequences from 27 Salvia species through DNA sequencing and 77 sequences within the same taxa from the GenBank. The discriminatory capabilities of these regions were evaluated in terms of PCR amplification success, intraspecific and interspecific divergence, DNA barcoding gaps, and identification efficiency via a tree-based method. ITS1 was superior to the other marker for discriminating between species, with an accuracy of 81.48%. The three combinations did not increase species discrimination. Finally, we found that ITS1 is a powerful barcode for identifying Salvia species, especially Salvia miltiorrhiza.     

Diversification and vicariance of desert plants: Evidence inferred from chloroplast DNA sequence variation of Lagochilus ilicifolius (Lamiaceae)

To document arid Northern China as a diversification center for desert plants, and to better understand the mechanisms of desert taxa diversification, we used five cpDNA spacers (trnL-trnF, rps16, psbA-trnH, psbK-psbI and trnS-trnG) to investigate Lagochilus ilicifolius in all distribution areas, including Northern China, Mongolia and Russia. Phylogenetic analyses showed that L. ilicifolius comprises two distinctive lineages, one distributed in China, and another in Mongolia–Russia. Our data confirmed that arid Northern China, is a distinctive area with many endemic genera. Biogeographic inferences, based on a Bayesian uncorrelated lognormal model together with molecular dating, suggested that the main diversification within the species occurred in the Pleistocene (ca. 1.38–0.3 Ma), resulting from the transition of the climate of Eurasia to a dry-cold pattern as well as the desertification caused the vicariance of desert plants.     

Five molecular markers reveal extensive morphological homoplasy and reticulate evolution in the Malva alliance (Malvaceae)

The Malva alliance is a well-defined group with extensive morphological homoplasy. As a result, the relationships among the taxa as well as the evolution of morphological traits have remained elusive and the traditional classifications are highly artificial. Using five molecular markers (nuclear ITS, plastid matK plus trnK, ndhF, trnL-trnF, psbA-trnH), we arrived at a phylogenetic hypothesis of this group, the genera Alcea, Althaea and Malvalthaea being studied here for the first time with molecular data. Althaea and, in particular, Lavatera and Malva are highly polyphyletic as currently circumscribed, because their diagnostic characters, the number and degree of fusion of the epicalyx bracts, evolve in a highly homoplasious manner. In contrast, fruit morphology largely agrees with the molecularly delimited groups. Hybrid origins confirmed for the genus Malvalthaea and for Lavatera mauritanica and hybridization in the group of ruderal small-flowered mallows underline the importance of reticulate evolution in shaping the history of this group and complicating the interpretation of morphological evolution.     

Evaluation of six candidate DNA barcoding loci in Ficus (Moraceae) of China

Ficus, with about 755 species, diverse habits and complicated co‐evolutionary history with fig wasps, is a notoriously difficult group in taxonomy. DNA barcoding is expected to bring light to the identification of Ficus but needs evaluation of candidate loci. Based on five plastid loci (rbcL, matK, trnH‐psbA, psbK‐psbI, atpF‐atpH) and a nuclear locus [internal transcribed spacer (ITS)], we calculated genetic distances and DNA barcoding gaps individually and in combination and constructed phylogenetic trees to test their ability to distinguish the species of the genus. A total of 228 samples representing 63 putative species in Ficus (Moraceae) of China were included in this study. The results demonstrated that ITS has the most variable sites, greater intra‐ and inter‐specific divergences, the highest species discrimination rate (72%) and higher primer universality among the single loci. It is followed by psbK‐psbI and trnH‐psbA with moderate variation and considerably lower species discrimination rates (about 19%), whereas matK, rbcL and atpF‐atpH could not effectively separate the species. Among the possible combinations of loci, ITS + trnH‐psbA performed best but only marginally improved species resolution over ITS alone (75% vs. 72%). Therefore, we recommend using ITS as a single DNA barcoding locus in Ficus.     

The use of a non-coding region of chloroplast DNA in phylogenetic studies of the subtribeSonchinae (Asteraceae:Lactuceae)

The systematic utility of sequences from a non-coding region of chloroplast DNA (cpDNA) betweenpsbA andtrnH^(GUG) was examined by assessing phylogenetic relationships in subtribeSonchinae (Asteraceae:Lactuceae). Primers constructed against highly conserved regions of tRNA genes were used for PCR amplification and sequencing. ThepsbA-trnH intergenic spacer contains several insertions and deletions (indels) inSonchinae with the length varying from 385 to 450 bp. Sequence divergence ranges from 0.00% to 7.54% withinSonchinae, with an average of 2.4%. Average sequence divergence inSonchus subg.Sonchus is 2.0%, while the mean for subg.Dendrosonchus and its close relatives in Macaronesia (the woodySonchus alliance) is 1.0%. Our results suggest that this region does not evolve rapidly enough to resolve relationships among closely related genera or insular endemics in theAsteraceae. The phylogenetic utility ofpsbA-trnH sequences of the non-coding cpDNA was compared to sequences from the ITS region of nuclear ribosomal DNA. The results suggest that ITS sequences evolve nearly four times faster thanpsbA-trnH intergenic spacer sequences. Furthermore, the ITS sequences provide more variable and phylogenetically informative sites and generate more highly resolved trees with more strongly supported clades, and thus are more suitable for phylogenetic comparisons at lower taxonomic levels than thepsbA-trnH intergenic chloroplast sequences.     

Molecular phylogenetic and biogeographical analysis of <Emphasis Type="Italic">Nitraria</Emphasis> based on nuclear and chloroplast DNA sequences

Based upon DNA sequences from six plastid regions (rbcL, psbB-psbH, trnL-trnF, rpS16, psbA-trnH, rpS16-trnK) and the internal transcribed spacer (ITS) region of nuclear ribosomal DNA, the phylogenetic relationships in the genus Nitraria and family Nitrariaceae are investigated by using methods of maximum parsimony, maximum likelihood, and Bayesian inference. Our study strongly supports the monophyly of Nitraria. Nitraria can be divided into four parts, namely, the N. sphaerocarpa group, N. retusa group, the N. roborowskii and N. tangutorum group, and a group consisting of N. schoberi, N. komarovii, N. sibirica, and N. billardieri. Ancestral area reconstruction using S-Diva shows that eastern Central Asia is most likely the place of origin, and then dispersals occurred to western Central Asia, Africa, and Australia.     

Barcoding Poplars (Populus L.) from Western China

Background

Populus is an ecologically and economically important genus of trees, but distinguishing between wild species is relatively difficult due to extensive interspecific hybridization and introgression, and the high level of intraspecific morphological variation. The DNA barcoding approach is a potential solution to this problem.

Methodology/Principal Findings

Here, we tested the discrimination power of five chloroplast barcodes and one nuclear barcode (ITS) among 95 trees that represent 21 Populus species from western China. Among all single barcode candidates, the discrimination power is highest for the nuclear ITS, progressively lower for chloroplast barcodes matK (M), trnG-psbK (G) and psbK-psbI (P), and trnH-psbA (H) and rbcL (R); the discrimination efficiency of the nuclear ITS (I) is also higher than any two-, three-, or even the five-locus combination of chloroplast barcodes. Among the five combinations of a single chloroplast barcode plus the nuclear ITS, H+I and P+I differentiated the highest and lowest portion of species, respectively. The highest discrimination rate for the barcodes or barcode combinations examined here is 55.0% (H+I), and usually discrimination failures occurred among species from sympatric or parapatric areas.

Conclusions/Significance

In this case study, we showed that when discriminating Populus species from western China, the nuclear ITS region represents a more promising barcode than any maternally inherited chloroplast region or combination of chloroplast regions. Meanwhile, combining the ITS region with chloroplast regions may improve the barcoding success rate and assist in detecting recent interspecific hybridizations. Failure to discriminate among several groups of Populus species from sympatric or parapatric areas may have been the result of incomplete lineage sorting, frequent interspecific hybridizations and introgressions. We agree with a previous proposal for constructing a tiered barcoding system in plants, especially for taxonomic groups that have complex evolutionary histories (e.g. Populus).     

Molecular and cytological evidences for the natural wheatgrass hybrids occurrence and origin in west China

Interspecies hybridization has been frequently observed in the tribe Triticeae. Natural hybridization between Kengyilia and Roegneria or Elymus species has not been reported as yet. Several sterile wheatgrass individuals exhibiting intermediately morphological traits between Kengyilia and Roegneria or Elymus species were identified in the meadow of Sichuan and Gansu provinces in China, suggesting their natural hybrid origin. The putative hybrids were analyzed by using the sequences of ITS and trnH-psbA together with cytological observation in order to assess the origin of hybrids. Both ITS and cytological data revealed the evidence of allopolyploid origin and confirmed the presence of StStYYP and StStYYHP genomes in the putative natural hybrids. The data suggest that the StStYYP hybrid originated from hybridization between Kengyilia and Roegneria and the hybrid with StStYYHP originated from hybridization between Kengyilia and Elymus. Chloroplast sequence data demonstrated that K. rigidula and K. melanthera were the likely maternal donors in the hybridization events.     

A range wide geographic pattern of genetic diversity and population structure of Hexinia polydichotoma (Asteraceae) in Tarim Basin and adjacent areas

In order to investigate the genetic diversity and influence of climate oscillations on evolutionary processes of organisms in northwest China, we selected Hexinia polydichotoma, a species endemic to China, and examined the geographic pattern of genetic variation in its entire cover range, Tarim Basin and adjacent areas. In the study, 22 chloroplast (cp) haplotypes were identified based on two cpDNA sequences (trnH–psbA and ycf6–psbM), and five ITS sequence variations were found. Shown in the cp haplotype network, the two common cp haplotypes mainly distributed along the northern and southern rims of the basin respectively and intersected in the center of the basin, whereas in the ITS network, ITS genotype 1 was widespread across the whole distribution area, and rare genotypes were concentrated in the western rim of the basin. Genetic variation primarily occurred among populations and SAMOVA groups. Fragmented desert habitat may have caused gene flow barrier among populations or groups far from each other, leading to significant genetic differentiation at these levels. It appears that expansion and contraction cycles of river systems and oases during the middle Pleistocene was the source of the observed fragmentation. Geographic range expansion in H. polydichotoma was supported by the significant value for Tajima's D, Fu's F_S, and by a unimodal mismatch distribution. It was possible that the enlargement of the Taklimakan Desert during the middle Pleistocene may have provided appropriate conditions for the range dispersal. We identified western rim of the basin as the center of genetic diversity of H. polydichotoma based on the present dataset.     

Phylogeny of the Neotropical sages (<Emphasis Type="Italic">Salvia</Emphasis> subg. <Emphasis Type="Italic">Calosphace</Emphasis>; Lamiaceae) and insights into pollinator and area shifts

Salvia subg. Calosphace (Lamiaceae, Lamiales) is a highly diverse clade endemic to the New World. The phylogenetic relationships of Calosphace have been previously investigated using DNA sequences of nuclear ITS region and plastid psbA–trnH intergenic spacer, but the resulting trees lack resolution and support for many clades. The present paper reassesses the phylogenetic relationships of subgenus Calosphace, including a broader taxon sampling, with a special focus on representing previously unsampled sections, and using an additional plastid marker (trnL–trnF region). Our results show increased resolution and overall patterns of support, recovering ten main clades. Within core Calosphace, the most inclusive of these main clades, 17 new subclades were identified. Of the 42 sections for which more than one species was analysed, only 12 are monophyletic. Our biogeographical analysis identified at least twelve migrations to South America from Mexican and Central American lineages, in agreement with previous suggestions of multiple origins of South American Calosphace diversity. This analysis also confirmed a colonization of the Antilles by Andean lineages. The reconstruction of ancestral states of pollination syndromes showed multiple shifts to ornithophily from melittophily and one reversal to the latter.     

Conservation genetics and geographic patterns of genetic variation of endangered shrub <Emphasis Type="Italic">Ammopiptanthus</Emphasis> (Fabaceae) in northwestern China

Phylogeographic patterns of Ammopiptanthus in northwestern China were examined with internal transcribed spacer (ITS) and three chloroplast intergenic spacers (trnH–psbA, trnL–trnF, and trnS–trnG). Two ITS genotypes (a–b) and 8 chloroplast haplotypes (A–H) were detected. Both ITS genotypes and chloroplast lineages were split in two geographic regions: western Xinjiang and the Alxa Desert. This lineage split was also supported by AMOVA analysis and the Mantel test. AMOVA showed that 89.81 % of variance in Ammopiptanthus occurred between the two geographic regions, and correlation between genetic distances and geographical distances was significant (r = 0.757, p < 0.0001). All populations in western Xinjiang shared haplotype A with high frequency, and range expansion was strongly supported by negative Fu’s F_S value, and mismatch distribution analysis, whereas populations in the Alxa Desert had higher genetic diversity and structure. We speculate that the cold and dry climate during the early Quaternary fragmented habitats of the species, limiting gene flow between regions, and interglacial periods most likely led to the range expansion in western Xinjiang. The low genetic diversity of Ammopiptanthus indicate a significant extinction risk, and protective measures should be taken immediately.     

Evolution and biogeography of Centaurea section Acrocentron inferred from nuclear and plastid DNA sequence analyses

Background and Aims

Section Acrocentron of the genus Centaurea is one of the largest sections of Centaurea with approx. 100 species. The geographic distribution, centred in the Mediterranean, makes it an excellent example for studies of the biogeographic history of this biodiversity-rich region.

Methods

Plastid (trnH-psbA) and nuclear (ITS and ETS) DNA sequence analysis was used for phylogenetic reconstruction. Ancestral biogeographic patterns were inferred by dispersal-vicariance analysis (DIVA).

Key Results

The resulting phylogeny has implications for the sectional classification of Acrocentron and confirms merging sect. Chamaecyanus into Acrocentron as a subsection. Previous suggestions of an eastern Mediterranean origin of the group are confirmed. The main centres of diversification established in previous studies are now strongly supported. Expansion of the group in two different radiations that followed patently diverse paths is inferred.

Conclusions

Radiation followed two waves, widely separated in time scale. The oldest one, from Turkey to Greece and the northern Balkans and then to North Africa and Iberia, should be dated at the end of the Miocene in the Messinian period. It reached the Iberian Peninsula from the south, following a route that is landmarked by several relictic taxa in Sicily and North Africa. A later radiation during the Holocene interglacial periods followed, involving species from the north of the Balkan Peninsula, along a Eurasian pathway running from Central Iberia to the steppes of Kazakhstan. A generalized pattern of reticulation is also evident from the results, indicating past contacts between presently separated species. Molecular data also confirmed the extent of hybridization within Acrocentron and were successful in reconstructing the paleogeography of the section.Key words: Centaurea sect. Acrocentron, Cardueae, dispersal-vicariance, ETS, ITS, Mediterranean, phylogeny, psbA-trnH, radiation, reticulation     

Application of DNA Barcodes in Asian Tropical Trees – A Case Study from Xishuangbanna Nature Reserve,Southwest China

Background

Within a regional floristic context, DNA barcoding is more useful to manage plant diversity inventories on a large scale and develop valuable conservation strategies. However, there are no DNA barcode studies from tropical areas of China, which represents one of the biodiversity hotspots around the world.

Methodology and Principal Findings

A DNA barcoding database of an Asian tropical trees with high diversity was established at Xishuangbanna Nature Reserve, Yunnan, southwest China using rbcL and matK as standard barcodes, as well as trnH–psbA and ITS as supplementary barcodes. The performance of tree species identification success was assessed using 2,052 accessions from four plots belonging to two vegetation types in the region by three methods: Neighbor-Joining, Maximum-Likelihood and BLAST. We corrected morphological field identification errors (9.6%) for the three plots using rbcL and matK based on Neighbor-Joining tree. The best barcode region for PCR and sequencing was rbcL (97.6%, 90.8%), followed by trnH–psbA (93.6%, 85.6%), while matK and ITS obtained relative low PCR and sequencing success rates. However, ITS performed best for both species (44.6–58.1%) and genus (72.8–76.2%) identification. With trnH–psbA slightly less effective for species identification. The two standard barcode rbcL and matK gave poor results for species identification (24.7–28.5% and 31.6–35.3%). Compared with other studies from comparable tropical forests (e.g. Cameroon, the Amazon and India), the overall performance of the four barcodes for species identification was lower for the Xishuangbanna Nature Reserve, possibly because of species/genus ratios and species composition between these tropical areas.

Conclusions/Significance

Although the core barcodes rbcL and matK were not suitable for species identification of tropical trees from Xishuangbanna Nature Reserve, they could still help with identification at the family and genus level. Considering the relative sequence recovery and the species identification performance, we recommend the use of trnH–psbA and ITS in combination as the preferred barcodes for tropical tree species identification in China.     

Phylogenetic analysis of Lappula Moench (Boraginaceae) based on molecular and morphological data

Lappula Moench includes ca. 70 species and exhibits a wide range of nutlet variation. Currently, the evolutionary relationships among species of Lappula have not been examined; therefore, to elucidate phylogenetic relationships and morphological evolution within Lappula and related genera, we conducted phylogenetic analyses with matrices that include 48 species as well as four DNA regions (ITS, trnL-trnF, rpS16 and psbA-trnH) and 18 morphological characters. These matrices were analyzed using maximum parsimony, maximum likelihood, and Bayesian inference methods. Analyses of the combined molecular and morphological data result in a phylogeny that is better resolved than that based solely on molecular sequence data. Phylogenetic results suggest that the current infrageneric classification of Lappula, at least at the subsectional and series level, is artificial. The evolutionary patterns of 18 morphological characters are investigated in a phylogenetic context. In Lappula, nutlet homomorphism and small corollas are resolved as ancestral, while nutlet heteromorphism and larger corollas are derived.     

Phylogeny of the Southeast Asian endemic genus Neocinnamomum H. Liu (Lauraceae)

A phylogenetic analysis of Neocinnamomum H. Liu and related genera was conducted using psbA–trnH, trnK cpDNA regions, and the ITS nrDNA segment. Neocinnamomum was confirmed to be monophyletic, and an evolutionary series of inflorescence development within the genus was recognized. The compound thyrse seen in N. caudatum is reduced to the few- to many-flowered condensed inflorescences with a poorly defined branching system seen in most species and ultimately to the 1-flowered inflorescence seen in N. atjehense. Consensus network analysis (CNA) suggested that long-branch attraction is responsible for the observed close relationship between Neocinnamomum and Cassytha L. in a combined analysis of the complete data. In contrast, the sister relationship of Neocinnamomum and Caryodaphnopsis seen in the Bayesian analyses of the partial combined matrix was supported by CNA and is also supported by morphology and wood and bark anatomy. The close similarity of the compound thyrse of less derived Neocinnamomum species to the thyrsoid inflorescences of some Caryodaphnopsis species is also seen as strong support for their affinity.     

Phylogeny and biogeography of Valerianaceae (Dipsacales) with special reference to the South American valerians

Species of Valerianaceae are a common component of the alpine flora throughout the Northern Hemisphere as well as the Andes of South America. Sequence data from three chloroplast markers (psbA-trnH intron, trnK-matK intron, and the trnL-F region) along with the internal transcribed spacer region (ITS) of nuclear ribosomal DNA were used to infer relationships within Valerianaceae. Both genomes, as well as a combined data set, provide support for the major clades within the group and do not support a monophyletic Valeriana. In addition, these data indicate that Plectritis is nested within South American Valeriana, as opposed to being sister to Centhranthus as previously hypothesized. Valerianaceae appear to have originated in Asia, probably in the Himalayas, and subsequently to have dispersed several times to Europe and to the New World. Our results imply that Valerianaceae colonized South America on multiple occasions from the north. In one of these cases there appears to have been a substantial and rapid radiation, primarily in the high elevation paramo habitat. A variety of methods were used to estimate divergence times to determine when Valerianaceae might have colonized South America. Regardless of the method and fossil constraints applied, our estimates suggest that Valerianaceae colonized South America prior to the formation of the Isthmus of Panama.