Molecular identification of species in Prunus sect. Persica (Rosaceae), with emphasis on evaluation of candidate barcodes for plants

Abstract Species of Prunus L. sect. Persica are not only important fruit trees, but also popular ornamental and medicinal plants. Correct identification of seedlings, barks, or fruit kernels is sometimes required, but no reliable morphological characters are available. Nowadays, the technique of DNA barcoding has the potential to meet such requirements. In this study, we evaluated the suitability of 11 DNA loci (atpB‐rbcL, trnH‐psbA, trnL‐F, trnS‐G, atpF‐H, rbcL, matK, rpoB, rpoC1, nad1, and internal transcribed spacer [ITS]) as candidate DNA barcodes for peaches, using samples from 38 populations, covering all the species in sect. Persica. On the whole, the primers worked well in this group and sequencing difficulties were met only in the case of ITS locus. Five loci (rbcL, matK, rpoB, rpoC, and nad1) have very low variation rates, whereas atpB‐rbcL, atpF‐H, trnH‐psbA, trnL‐F and trnS‐G show more variability. The most variable loci, atpB‐rbcL and trnH‐psbA, can distinguish three of the five species. Two two‐locus combinations, atpB‐rbcL+trnL‐F and atpB‐rbcL+atpF‐H, can resolve all five species. We also find that identification powers of the loci are method‐dependent. The NeighborNet method shows higher species identification power than maximum parsimony, neighbor joining, and unweighted pair group method with arithmetic mean methods.     

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.     

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.     

Molecular phylogeny of Dipterocarpaceae in Southeast Asia using RFLP of PCR-amplified chloroplast genes

Dipterocarpaceae is the dominant family of Southeast Asia's climax tropical rain forest region, and it contains the region's most important commercial timber species. A molecular phylogeny of the Dipterocarpaceae subfamily Dipterocapoideae was constructed using restriction fragment length polymorphisms of polymerase chain reaction-amplified specific genes in chloroplast DNA. A total of 141 site changes were detected among ten genera and 30 species in 11 different genes: rbcL, psbA, psbD, rpoB, rpoC, petB, atpH, 16S, psaA, petA and trnK. Phylogenetic trees constructed by Wanger parsimony and neighbor-joining methods, using Upuna as the outgroup, displayed five monophytelic groups that included Upuna: HopeaShorea-Parashorea-Neobalanocarpus; Dryobalanops; Dipterocarpus; Anisoptera-Vatica-Cotylelobium; and Upuna. The phylogenetic trees clearly separate species with two different base chromosome numbers: the first group is x=7, and the other is x=11. The x=7 group is thought to be in a synapomorphic character state. Parashorea lucida is a sister to most Shorea species. Neobalanocarpus heimii and Hopea from a clade of a sister to two Shorea species, and Cotylelobium and Vatica are closely related species. Our conclusions agree with a phylogeny derived from wood anatomy data analysis, and with Symington's and Ashton's taxonomic classifications.The raw data of the PCR-RFLP analysis can be obtained from the authors     

RFLP analysis of cpDNA in the genus <Emphasis Type="Italic">Hypericum</Emphasis>

The chloroplast DNA of 43 species including 16 sections from the genus Hypericum was studied by PCR-RFLP analysis. The PCR-amplified products of four cpDNA regions, trnC-trnD, psbC-trnS, trnL-trnF and rbcL were digested with four restriction endonucleases. A high level of interspecific variation was detected while intraspecific diversity was not observed. The resulting parsimony analysis indicated the monophyletic assemblage of the sections Androsaemum, Olympia, Drosocarpium and Trigynobrathys. Monophyly of Hypericum is weakly supported, but close relationships of H. perforatum and H. maculatum are indicated. The members of Ascyreia are weakly resolved, but clustering of H. kouytchense and H. oblongifolium is well supported, however, H. reptans is nested with Olympia. CpDNA profiles and the positions on the parsimony tree indicate that the chloroplast donor among the putative parents of the hybrid species H. ×inodorum is H. androsaemum.     

PHYLOGENETIC ANALYSIS OF 32 STRAINS OF VAUCHERIA (XANTHOPHYCEAE) USING THE rbcL GENE AND ITS TWO FLANKING SPACER REGIONS1

The complete rbcL gene was sequenced for 21 species and 32 strains of Vaucheria and for five other Xanthophyceae (Asterosiphon dichotomus (Kützing) Rieth, Botrydium becharianum Vischer, B. cystosum Vischer, B. stoloniferum Mitra, Tribonema intermixtum Pascher). The psbA‐rbcL spacer, upstream of the rbcL gene, and the RUBISCO spacer between the rbcL and rbcS genes were also completely sequenced for the Vaucheria strains and Asterosiphon. The psbA‐rbcL spacer was the most variable region that was sequenced, and only the 3′ end of the spacer could be aligned. Phylogenetic analyses (maximum parsimony, neighbor joining, and maximum likelihood) were conducted using the DNA sequence and the amino acid sequence for the rbcL gene, and a second analysis was conducted using a portion of the psbA‐rbcL spacer +rbcL gene + RUBISCO spacer. All analyses showed that Vaucheria species formed monophyletic clades that corresponded with morphologically based subgeneric sections, including the section Racemosae. Species producing a gametophore (= fruiting branch, bearing both an antheridium and oogonium) formed a monophyletic clade in all analyses. The nongametophore species sometimes formed a monophyletic clade but other times formed a basal grade. Pair‐wise comparisons of nucleotides and amino acids showed that for some species, numerous nucleotide changes resulted in relatively few amino acid changes. Consequently, phylogenetic analysis of the amino acids produced numerous trees, which in a strict consensus tree resulted in numerous polychotomies. An original strain of V. terrestris that was deposited in two culture collections over 25 years ago had identical sequences, suggesting no rapid change was occurring in the sequenced regions. Two strains of V. prona, isolated from Europe and North America, had identical sequences. Other species, for which two or more strains were examined, had different sequences. These results suggest that cryptic species complexes exist within Vaucheria because the rbcL gene is a conservative gene that is identical in other protists.     

Species identification of the medicinal plant Tulipa edulis (Liliaceae) by DNA barcode marker

Tulipa edulis (Liliaceae) is the botanical origin of the traditional Chinese medicine (TCM) “Guangcigu”. Due to overexploitation that induced a decline in natural sources, many dried bulbs from other species of Tulipa have been used, adulterating the medicine in recent years. This practice may cause a series of inconsistent therapeutic effects and quality control problems in the herbal medicine industry. Hence, three DNA regions (matK, psbA-trnH and rbcL) were evaluated as barcodes for identifying T. edulis and its adulterants. All candidate DNA barcodes were successfully amplified from leaf samples. Based on the sequence divergences, rbcL and psbA-trnH can assign T. edulis and its adulterants to the correct genus, while matK can accurately differentiate T. edulis and its adulterants. Thus, at the DNA level, the matK intergenic region is a more suitable, accurate and applicable identification of T. edulis and its adulterants than rbcL and psbA-trnH.     

Phylogenetic relationships between disjunctly occurring groups of Tristicha trifaria (Podostemaceae)

Aim To reveal the phylogeographic relationship of disjunct specimens of Tristicha trifaria (Bory ex Willd.) Spreng., a member of the Podostemaceae river‐weed family, which is distributed exceptionally widely, but disjunctly, in Africa and the Americas. Location Brazil, Mexico, Ghana, Tanzania and Madagascar. Methods The chloroplast matK and rbcL genes, a trnK intron, the trnS‐trnG intergenic spacer (IGS), the two IGSs of trnT‐trnL‐trnF, a trnL intron, and nuclear ribosomal ITS regions were sequenced and analysed. Phylogenetic analyses were conducted using maximum likelihood and maximum parsimony methods. Results The T. trifaria samples analysed were separated into two groups in a rooted tree based on a combined matK/rbcL/ITS dataset; one contained the West African and all of the American samples, and the other contained the East African and Madagascan samples. An unrooted tree obtained from a combined analysis of all the chloroplast DNA and nuclear ITS data showed that a sample from West Africa was sister to an American T. trifaria group. Main conclusions The American and West African T. trifaria are closely related, despite the great distance between their locations. This observation, along with a tree of the whole Tristichoideae subfamily and estimated divergence times, suggests that an ancestor of T. trifaria migrated from Asia to Africa during the early Tertiary, and that this was followed by further westward migration to the Americas at the end of the Miocene or in the early Pliocene.     

Molecular phylogeny of mangroves VII. PCR-RFLP of trnS-psbC and rbcL gene regions in 24 mangrove and mangrove-associate species

Chloroplast DNA (cpDNA) regions, trnS-psbC and rbcL, from 120 individuals of 24 mangrove and mangrove associate species belonging to 11 orders, 13 families and 17 genera of Angiospermae were amplified by the polymerase chain reaction (PCR) and restriction-digested with HaeIII. Analysis of polymorphism in the restriction fragments (PCR-RFLP) revealed 18 classes of restriction banding pattern in trnS-psbC region. This has provided molecular evidence for diversity in the mangrove floral component at the above-species level. Intra-generic variations were observed in three genera, viz. Rhizophora, Avicennia and Suaeda. Species-specific restriction patterns were found in the genera Rhizophora and Suaeda. A natural hybrid belonging to the genus Rhizophora was also analysed, and its restriction pattern was the same as that of a putative parental species.PCR-RFLP analysis of rbcL gene region was less differentiating. However, it showed 13 different classes of restriction patterns and revealed the usefulness of these investigations for genome analysis at a higher taxonomic level. Intra-specific variation was not observed in any of the species in either of the cpDNA regions analysed. This is the first report which describes variations in the chloroplast genome of mangrove species. Received: 20 April 1999 / Accepted: 12 May 1999     

Phylogenetic relationships in Taxodiaceae and Cupressaceae sensu stricto based on matK gene, chlL gene, trnL-trnF IGS region, and trnL intron sequences

Nucleotide sequences from four chloroplast genes, the matK, chlL, intergenic spacer (IGS) region between trnL and trnF, and an intron of trnL, were determined from all species of Taxodiaceae and five species of Cupressaceae sensu stricto (s.s.). Phylogenetic trees were constructed using the maximum parsimony and the neighbor-joining methods with Cunninghamia as an outgroup. These analyses provided greater resolution of relationships among genera and higher bootstrap supports for clades compared to previous analyses. Results indicate that Taiwania diverged first, and then Athrotaxis diverged from the remaining genera. Metasequoia, Sequoia, and Sequoiadendron form a clade. Taxodium and Glyptostrobus form a clade, which is the sister to Cryptomeria. Cupressaceae s.s. are derived from within Taxodiaceae, being the most closely related to the Cryptomeria/Taxodium/Glyptostrobus clade. These relationships are consistent with previous morphological groupings and the analyses of molecular data. In addition, we found acceleration of evolutionary rates in Cupressaceae s.s. Possible causes for the acceleration are discussed.     

The evolution of actinorhizal symbioses: Evidence for multiple origins of the symbiotic association

According to morphologically based classification systems, actinorhizal plants, engaged in nitrogen-fixing symbioses with Frankia bacteria, are considered to be only distantly related. However, recent phylogenetic analyses of seed plants based on chloroplast rbcL gene sequences have suggested closer relationships among actinorhizal plants. A more thorough sampling of chloroplast rbcL gene sequences from actinorhizal plants and their nonsymbiotic close relatives was conducted in an effort to better understand the phylogenetic relationships of these plants, and ultimately, to assess the homology of the different actinorhizal symbioses. Sequence data from 70 taxa were analyzed using parsimony analysis. Strict consensus trees based on 24 equally parsimonious trees revealed evolutionary divergence between groups of actinorhizal species suggesting that not all symbioses are homologous. The arrangement of actinorhizal species, interspersed with nonactinorhizal taxa, is suggestive of multiple origins of the actinorhizal symbiosis. Morphological and anatomical characteristics of nodules from different actinorhizal hosts were mapped onto the rbclL-based consensus tree to further assess homology among rbcL-based actinorhizal groups. The morphological and anatomical features provide additional support for the rbcL-based groupings, and thus, together, suggest that actinorhizal symbioses have originated more than once in evolutionary history.     

Nutrient concentrations in litterfall from some western conifers with special reference to calcium

Foliar litterfall nutrient concentrations were analysed for selected members of Taxodiaceae and Cupressaceae families andPseudotsuga menziesii for two arboreta in western Oregon and Washington. Nutrient results between arboreta show similar concentrations with the exception of magnesium, which may be the result of historical land use. Nutrient concentrations between species vary considerably.Pseudotsuga menziesii is particularly distinctive from the Cupressaceae and Taxodiaceae by retaining large amounts of phosphorus and potassium. Taxodiaceae is distinctive by high concentration of Mg while Cupressaceae retains calcium, especiallyChamaecyparis nootkatensis. Results suggest that all members of Taxodiaceae and Cupressaceae retain considerably more Ca than Pinaceae in foliar litter.     

Interspecific relationships and origins of Taxaceae and Cephalotaxaceae revealed by partitioned Bayesian analyses of chloroplast and nuclear DNA sequences

The precise delimitation of Taxaceae and Cephalotaxaceae is not totally resolved. Some contradicting taxonomic proposals have been published, which demonstrates the difficulties in establishing a natural classification of the families and especially in proposing a relevant treatment within the genera Taxus and Cephalotaxus. The aims of this study are to contribute to the phylogeny and specific delineation of the two conifer families on the basis of molecular data. A cladistic analysis of the sequences of five chloroplast (matK, rbcL, trnL, trnL-trnF spacer, and psbA-trnH spacer) and one nuclear (ITS) molecular markers was carried out, both individually and in combination, by distance, parsimony, likelihood, and Bayesian methods. The results confirm that the two families are monophyletic. In the genus Taxus, T. floridana is the first-branching taxon; T. brevifolia and T. globosa cluster together and are sister to T. baccata; the endemic T. yunnanensis clusters with T. wallichiana in subclade B and is only distantly related with the other four Taxus species in China (subclade A); T. fuana is closer to T. baccata than to other Taxus species. Torreya jackii and A. formosana are the first-branching species within Torreya and Amentotaxus, respectively. C. koreana and C. wilsoniana could be treated as two varieties of C. harringtonia. The ancestral distribution area of Taxaceae and Cephalotaxaceae is restricted either to Southwest China or Southeast China by DIVA analysis. The relaxed molecular clock indicates that the deepest divergences in Taxus go back to the late-Cretaceous. psbA-trnH, rbcL third codon position, and matK first codon position contributed most to the separation of taxa in Discriminant function analysis. Our results confirm, on a basis of multiple molecular markers and a complete sampling of basic species, the suggested monophyly of Taxaceae and Cephalotaxaceae and propose interspecific relationships within each group, with profound nomenclatural and taxonomic implications. Combination of partitioned Bayesian analysis and likelihood-based methods produced a more robust phylogenetic hypothesis for the two studied families.     

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.     

The architecture of the chloroplast psbA-trnH non-coding region in angiosperms

The psbA-trnH intergenic region is among the most variable regions in the angiosperm chloroplast genome. It is a popular tool for plant population genetics and species level phylogenetics and has been proposed as suitable for DNA barcoding studies. This region contains two parts differing in their evolutionary conservation: 1) the psbA 3′UTR (untranslated region) and 2) the psbA-trnH intergenic non-transcribed spacer. We compared the sequence and RNA secondary structure of the psbA 3′ UTR across angiosperms and found consensus motifs corresponding to the stem portions of the RNA stem-loop structures and a consensus TTAGTGTATA box. The psbA-trnH spacer exhibited patterns that can be explained by the independent evolution of large inversions in the psbA 3′UTR and mutational hot spots in the remaining portion of the psbA-trnH spacer. We conclude that a comparison of chloroplast UTRs across angiosperms offer clues to the identity of putative regulatory elements and information about selective constraints imposed on the chloroplast non-coding regions.     

中国葫芦科1新记录种——厚叶棒锤瓜

报道了葫芦科厚叶棒锤瓜[Neoalsomitra sarcophylla(Wall.)Hutch.]在中国的分布新记录。该物种形态上与N.balansae(Gagnep.)Hutch.近似,但蒴果长3~4cm;种子长6~7mm,宽3~4mm,两端呈短角状。本研究利用DNA条形码技术对该物种进行测序,获得matK、rbcL、psbA-trnH 3个基因序列。应用Blast法为该物种的分类处理提供佐证。凭证标本存放于广西药用植物园标本馆(GXMG)。     

Molecular phylogeny and evolution of host-plant use in conifer seed chalcids in the genus Megastigmus (Hymenoptera: Torymidae)

Abstract. Phylogenetic relationships amongst Megastigmus species (Chalcidoidea: Torymidae) associated with conifer seeds were inferred from DNA sequence data. Twenty‐nine species of seed chalcids were analysed using two different genes, cytochrome b (mitochondrial DNA) and the D2 domain of the 28S ribosomal DNA. Maximum‐parsimony and maximum‐likelihood analyses showed that taxa formed two monophyletic groups, one clade comprising all species associated with Cupressaceae and Taxodiaceae hosts with the exception of Chamaecyparis, and the other clade composed of species associated with Pinaceae. Species infesting Cupressaceae and Taxodiaceae seemed to be specialized to particular host genera or even to be species specific, which was consistent with a taxonomic radiation following initial host adaptation. By contrast, Megastigmus species associated with Pinaceae appeared capable of shifting onto different congeneric species or even onto a new host genus, with their evolution apparently less constrained by plant association. We hypothesized that the Megastigmus group associated with Pinaceae may have a much higher invasive potential than that related to Cupressaceae. The study also confirmed the presence of invasive Nearctic species in the Palaearctic, and demonstrated the existence of a cryptic species complex.     

Systematic Evaluation of Isoëtes asiatica Makino (Isoëtaceae) based on AFLP, nrITS, and Chloroplast DNA Sequences

On the basis of amplified fragment length polymorphism (AFLP) and nucleotide sequence data from nuclear ribosomal internal transcribed spacer (nrITS) and three chloroplast DNA regions (rbcL, cpITS, and trnS-psbC spacer), we investigated the species delimitation and the evolutionary lineage of Isoëtes asiatica from Hokkaido, Japan. The neighbor-joining (NJ) dendrogram based on AFLP markers revealed the well-defined clusters (bootstrap value?=?100%) of I. asiatica. Results from the principal component analysis are largely congruent with those obtained in the NJ dendrogram. The maximum parsimony analysis, based on data from nrITS and three chloroplast DNA sequences, supported a monophyly of three species, I. asiatica, Isoëtes echinospora, and Isoëtes maritima from Hokkaido, Kamchatka, and Alaska regions. The distinct species status of I. asiatica was also well supported in the combined chloroplast DNA phylogeny. Therefore, I. asiatica appear to represent example of gradual speciation due to spatial isolation of ancestral populations followed by genetic divergence. Our results also suggest that I. asiatica is probably not the ancestral diploid of the polyploids occurring in East Asia.     

Phylogenetic relationships in gnetophyta deduced fromrbcL gene sequences

Part of the large subunit of the ribulose-1,5-bisphosphate carboxylase gene (rbcL) was sequenced (1333 base pairs) from three species of gymnosperms:Ephedra sinica Gnetum parvifolium, Welwitschia mirabilis. Phylogenetic trees inferred from the neighbor joining, Wagner parsimony and maximum likelihood methods showed thatGnetum andWelwitschia were more closely related to each other than either is toEphedra within Gnetophyta, and the result supports previous cladistical analysis of morphological data.     

Molecular phylogeny of Dipterocarpaceae in Indonesia based on chloroplast DNA

In order to construct a molecular phylogeny of Indonesian Dipterocarpoideae (Dipterocarpaceae), PCR-RFLP of the chloroplast regions rbcL, petB, psbA, psaA, and trnL-F was performed with seven restriction enzymes in 129 samples including 58 species from nine genera. In the strict consensus tree with Monotes kerstingii as outgroup Indonesian Dipterocarpaceae were divided into two major clades. One clade (bootstrap value=71) consisted of Upuna, Cotylelobium, Anisoptera, Vatica, Dipterocarpus (tribe Dipterocarpeae, bootstrap value=83) and Dryobalanops (tribe Shoreae, bootstrap value=99) in a basal position. The second clade consisted of Hopea, Parashorea, and Shorea (tribe Shoreae) with 95% bootstrap support. Tribe Dipterocarpeae is monophyletic, tribe Shoreae is polyphyletic since Dryobalanops is sister to tribe Dipterocarpeae. In the neighbour-joining tree the sister group position of Dryobalanops to tribe Dipterocarpeae is not supported by the bootstrap analysis. Alternatively, we used Upuna borneensis as outgroup. The effect of outgroup selection on tree topology, taxonomic classification and the interpretation of character evolution is discussed.