Original Research ARTICLE
DNA Barcoding Evaluation and Its Taxonomic Implications in the Recently Evolved Genus Oberonia Lindl. (Orchidaceae) in China
- 1Key Laboratory of Plant Resources Conservation and Sustainable Utilization, Guangdong Provincial Key Laboratory of Applied Botany, South China Botanical Garden, Chinese Academy of Sciences, Guangzhou, China
- 2College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
- 3Shanghai Chenshan Plant Science Research Center, Chinese Academy of Sciences, Shanghai, China
- 4Shanghai Key Laboratory of Plant Functional Genomics and Resources, Shanghai Chenshan Botanical Garden, Shanghai, China
The orchid genus Oberonia Lindl., is a taxonomically complex genus characterized by recent species radiations and many closely related species. All Oberonia species are under conservation as listed in the CITES and the IUCN Red List Categories and Criteria. Given its difficulties in taxonomy and conservation status, Oberonia is an excellent model for developing DNA barcodes. Three analytical methods and five DNA barcoding regions (rbcL, matK, trnH-psbA, ITS, and ITS2) were evaluated on 127 individuals representing 40 species and 1 variety of Oberonia from China. All the three plastid candidates tested (rbcL, matK, and trnH-psbA) have a lower discriminatory power than the nuclear regions (ITS and ITS2), and ITS had the highest resolution rate (82.14%). Two to four combinations of these gene sets were not better than the ITS alone, but when considering modes of inheritance, rbcL+ITS and matK+ITS were the best barcodes for identifying Oberonia species. Furthermore, the present barcoding system has many new insights in the current Oberonia taxonomy, such as correcting species identification, resolving taxonomic uncertainties, and the underlying presence of new or cryptic species in a genus with a complex speciation history.
Oberonia Lindl. is a monophyletic genus (Tang G. D. et al., 2015) with its infrageneric classification still unclear with involving many closely related and recently radiated species. It comprises about 150–200 species centered in tropical Asia and extending to tropical Africa, NE Australia and the Pacific islands. There are 44 species and 2 varieties distributed in China (Su, 2000; Chen et al., 2009; Lin and Lin, 2009; Ormerod, 2010; Xu et al., 2010; Tian et al., 2013), including some taxonomically complex groups lacking clear species delimitations, such as O. acaulis Griff., O. acaulis var. luchunensis S. C. Chen and O. gongshanensis Ormerod; O. arisanensis Hayata, O. delicata Z. H. Tsi and S. C. Chen, O. japonica (Maxim.) Makino and O. menghaiensis S. C. Chen; O. austro-yunnanensis S. C. Chen and Z. H. Tsi and O. jenkinsiana Griff. ex Lindl.; O. insularis Hayata and O. sinica (S. C. Chen and K. Y. Lang) Ormerod, etc. All of these taxonomically complex groups show slightly differences in the morphology of their leaves or flowers, while in some cases, characters are overlaps between species. For example, O. japonica and O. delicata are morphologically similar in most aspects except for some differences in flower character. The former is characterized by sepals broader than petals (sepals 0.9 mm, petals 0.7 mm), while the latter has equal sepals and petals in width. However, flowers with sepals broader than petals (sepals 0.8 mm, petals 0.6 mm) were also found in the population of O. delicata from type location in Yunnan Province (LYL037). The overlapping characters among species make the discrimination and delineation of this genus more challenging. In addition, Oberonia species are vegetative consensus and have diminutive flowers that are not easily discernable by the naked eye, making their identification and taxonomy exceedingly difficult even for professional taxonomists. Thus, a tool such as DNA barcoding, to aid rapid and accurate identification of these species is critically needed.
Oberonia species in China are mostly distributed in Yunnan, Guangxi, Guangdong, and Hainan provinces which belong to the Indo-Burma hot spots. All of them are endangered and listed in CITES (Conventions on International Trade of Endangered Species of Fauna and Flora, http://www.cites.org/eng/disc/species.shtml) and the IUCN Red List categories and criteria (http://www.zhb.gov.cn/gkml/hbb/bgg/201309/t20130912_260061.htm.) Habitat loss and illegal collection in this region poses a great threat to species survival, particularly in the case of narrow endemic Oberonia species that confined to one location in Yunnan Province (Cardinale et al., 2012). Rapidly and correctly identifying Oberonia species in China could promote the monitoring of endangered taxa.
DNA barcoding is a relatively new tool for species identification (Hebert et al., 2003, 2004) and has been applied in many areas, such as taxonomy (Meier et al., 2006; Huang et al., 2013; Mutanen et al., 2015), the discovery of new or cryptic species (e.g., Burns et al., 2008; Karanovic, 2015; Saitoh et al., 2015), biodiversity assessment and conservation (e.g., Taberlet et al., 2012; Ji et al., 2013; Liu et al., 2014), and monitoring the illegal wildlife trades (e.g., Baker, 2008; Gathier et al., 2013). After more than 10 years development, many plastid regions and nuclear regions have been suggested as universal barcodes for land plants, such as rbcL+matK, trnH-psbA, ITS, ITS2, etc. (CBOL Plant Working Group, 2009; Chen et al., 2010; Ren et al., 2010; Li et al., 2011). Although, significant progress has been made in the DNA barcoding, the discrimination of closely related species in recently evolved genera such as Oberonia Lindl., Willows (Salix L.) and Curcuma L. (Twyford, 2014; Chen et al., 2015), still poses a great challenge. Hence, testing DNA barcodes in such taxonomically difficult Oberonia genus could help to further understand the potential of these barcodes. The establishment of an available barcoding system for Oberonia could also facilitate the taxonomy and conservation of these taxa.
Five DNA barcode regions were assessed (rbcL, matK, psbA-trnH, ITS, and ITS2) in 40 species and 1 variety of Oberonia obtained from China. The objectives of this study are to: test the effectiveness of suggested core DNA barcodes (rbcL+matK) in Oberonia; evaluate the resolution of these five barcodes and in 2- to 4-region combinations to correctly identify individuals and discriminate among closely related species; and explore some of the taxonomic implications in Oberonia.
Materials and Methods
Taxon Sampling, DNA Extraction, Amplification, Sequencing, and Sequence Downloads
In this study, 127 sequences from 40 species and 1 variety were collected for DNA sequencing, in which 123 sequences representing 39 species and 1 variety were from 6 provinces (Guangdong, Guangxi, Fujian, Hainan, Yunnan, and Tibet) in China. Four sequences from 1 species in Genbank were also included. In order to cover the morphological variability and geographical ranges, more individuals (>7) of widespread species such as O. caulescens Lindl. and O. ensiformis (Sm.) Lindl. were sampled. To test the potential of DNA barcodes, individuals of the taxonomically complex groups (e.g., O. acaulis, O. acaulis var. luchunensis, and O. gongshanensis; O. jenkinsiana and O. austro-yunnanensis; O. arisanensis, O. delicata, O. japonica, and O. menghaiensis; O. insularis and O. sinica) were sampled as much as possible. Detailed information about the samples is included in Table S1. The voucher specimens are deposited in the herbarium of the South China Botanical Garden, Chinese Academy of Sciences, Guangzhou (IBSC).
Total DNA was isolated from fresh or silica-dried leaves using a modified CTAB method (Doyle and Doyle, 1987). Primers are listed in Table S2. The ITS2 sequence was derived from the ITS (ITS, including ITS1, 5.8 s and ITS2) data directly. Two primer pairs of trnH-psbA were used for amplification and sequencing. Polymerase chain reaction (PCR) was conducted in a reaction mix (30 μl) each containing 10–20 ng (1–2 μl) of template DNA, 15 μl of 2 × PCR mix (0.005 units/μl Taq DNA polymerase, 4 mM MgCl2, and 0.4 mM dNTPs, TIANGEN), 1 μl of each primer (10 μM) and 11–12 μl of ddH2O. The PCR program started with a 2 min pre-melt stage at 98°C, followed by 36 cycles of 10 s at 98°C, annealing at 51–54°C (51°C for rbcL and ITS, 52°C for matK, 54°C for both two primer pairs of trnH-psbA) for 30 s, followed by 50 s at 68°C, and a final 8 min extension at 68°C. The PCR products were run on 1% agarose gels to check the quality of the amplified DNA. Then, PCR products with high quality were sent to Invitrogen (Shanghai) for purification and sequencing from both directions to reduce sequencing error.
Sequences for each marker were edited and assembled using SEQUENCHER 4.14 (GenCodes, Corp. Ann Arbor) and then manually adjusted using Bioedit v22.214.171.124 (Hall, 1999). To assess the barcoding resolution for all barcodes (rbcL, matK, ITS, ITS2, trnH-psbA, rbcL+matK, rbcL+ITS, rbcL+ITS2, rbcL+trnH-psbA, matK+ITS, matK+ITS2, matK+trnH-psbA, rbcL+matK+ITS, rbcL+matK+ITS2, rbcL+matK+trnH-psbA, rbcL+matK+ITS+trnH-psbA), three analytical methods were employed, i.e., the pair-wise genetic distance method (PWG-distance), the sequence similarity method (TAXONDNA) and a phylogenetic-based method (Neighbor-Joining trees and Bayesian inference trees).
For the pair-wise genetic-based method, six parameters [average inter-specific distance, average theta (Θ) prime and smallest inter-specific distance; average intra-specific distance, average theta (Θ) and largest intra-specific distance] were calculated in Mega7 using the Kimura two-parameter distance model (K2P), to explore the intra- and inter-species variation (Chen et al., 2010; Pang et al., 2013; Kumar et al., 2016) We considered discrimination to be successful if the minimum inter-species distance involving a species, represented by more than one individual was larger than its maximum intra-species distance.
The sequence similarity method used the proportion of correct identifications identified with TAXONDNA/ Species Identifier 1.8 program, to assess the potential of all markers for accurate species identification. The “Best Match” (BM) and “Best Close Match” (BCM) tests in TAXONDNA were run for species that were represented by more than one individual (Meier et al., 2006).
For the phylogenetic-based method, Neighbor Joining (NJ) trees of all markers were conducted in Mega7 with K2P model and Bayesian inference (BI) trees were conducted in MrBayes v. 3.1 (Huelsenbeck and Ronquist, 2001; Ronquist and Huelsenbeck, 2003). If all the individuals of one species were clustered into a monophyletic group with the support nodes above 70% (NJ) or 95% (BI), then the species was considered as successfully identified. Species with a single specimen were included, but lacked depth to calculate significance. Dendrobium strongylanthum Rchb. f. (KF177656, KF143723, KF177553, and GU339107) was the outgroup for the tree-based analysis (Tang G. D. et al., 2015; Xu et al., 2015).
PCR Amplification, Sequencing, and Genetic Divergence
The characteristics of five DNA barcoding regions and their combinations are shown in Table 1. When evaluated separately, rbcL, matK, and ITS had very high success rates (100%) for PCR amplification and sequencing using a single primer pair. Although, trnH-psbA also exhibited relatively high amplification success of 95.12% with two primer pairs, only 7 samples (5.69%) were successfully amplified and sequenced using the commonly used primer pair trnH2/psbAF (Sang et al., 1997; Tate and Simpson, 2003). The remaining samples were, however, generated for trnH-psbA using another primer trnH(GUG)/psbA (Hamilton, 1999). Overall, the aligned length of the five markers ranged from 269 bp (ITS2) to 1187 bp (rbcL). A total of 486 new sequences were submitted to NCBI, which included 123 sequences of each rbcL, matK, and ITS, separately, and 117 sequences of trnH-psbA, were submitted to NCBI (Table S1). In addition, we downloaded 12 sequences of O. mucronata (D. Don) Ormerod and Seidenf. from NCBI, including 4 sequences each of rbcL, matK, and ITS, respectively (JN005593, JN005592, JN005591, JN005590; JN004531, JN004530, JN004529, JN004528; JN114637, JN114636, JN114635, and JN114634). In total, 127 accessions of rbcL, matK, and ITS representing 40 species and 1 variety and 117 accessions of trnH-psbA representing 39 species were obtained for further analysis.
Among the five markers, ITS2 had the highest proportion of variable sites (52.42%) and parsimony-informative sites (41.26%), and rbcL had the lowest. RbcL showed the lowest intra-specific and inter-specific divergence as well, whilst trnH-psbA showed the highest intra-specific divergence (0.82%), followed by ITS2 (0.57%). However, the greatest inter-specific distance was found in ITS2 (9.92%), followed by ITS (6.74%).
The average intra- and inter-specific divergence in the 11 combinations varied from 0.14 to 0.45% and 1.00 to 4.21%, respectively (Table 1). The combinations of matK+trnH-psbA and matK+ITS showed the highest intra- and inter-specific genetic divergence, respectively. The core barcode rbcL+matK and the combination of rbcL+trnH-psbA had the lowest intra-specific genetic divergence, respectively (Table 1).
DNA Barcoding Gap Assessment
The relative distribution of K2P distances based on single barcodes and their combinations indicated that ITS and ITS2 had relatively clear barcoding gaps, while the remaining three tested candidate barcodes and their combinations had overlaps between their inter- and intraspecific distances (Figure S1).
Species Resolution of Candidate Barcodes
For the PWG-distance method, we used the local barcoding gap to reveal the species resolution power of candidate barcodes. That is, when a minimum inter-specific distance larger than the maximum intra-specific distance of a species, we considered it successfully identified. The proportion of the local barcoding gap varied among the markers tested (Figure S1 and Figure 1). Among single regions, the ITS exhibited the best species resolution (82.14%), followed by ITS2 (71.43%). In contrast, rbcL had the lowest species resolution (25.81%). Of the 11 combinations, rbcL+ITS showed the highest species resolution (82.14%), while rbcL+trnH-psbA showed the lowest species resolution (71.86%).
Figure 1. Species discrimination rates of multiple candidate barcodes in Oberonia. R, rbcL; M, matK; I, ITS; I2, ITS2; T, trnH-psbA.
The TAXONDNA analysis based on BM and BCM methods exhibited similar discrimination success (Table 2). The ITS had the highest success rate for the correct identification of species (BM and BCM: 73.22%) among the five single barcodes, followed by the ITS2 (BM and BCM: 64.56%), whilst rbcL had the lowest resolution rate (BM and BCM: 33.85%). Among the 11 combinations, rbcL+ITS and matK+ITS performed the best (BM and BCM: 74.59%), followed by rbcL+matK+ITS (BM and BCM: 74.80%). However, the core barcode rbcL+matK demonstrated the lowest species resolution rate (BM and BCM: 55.11%).
The Neighbor Joining (NJ) and Bayesian inference (BI) tree-building method showed similar discriminatory results (Figure 1). Among the five single barcodes, the ITS demonstrated the best discrimination power (NJ tree and BI tree: 75%), followed by ITS2, while rbcL and trnH-psbA had the lowest discrimination power. In the 11 combination barcodes, matK+ITS reached the highest resolution power (NJ tree: 78.57%; BI tree: 75%), followed by rbcL+ITS, rbcL+matK+ITS, and rbcL+matK+ITS+trnH-psbA. The core barcode rbcL+matK recommended by COBL had poor resolution (NJ tree and BI tree: 64.29%), and only slightly better than all combinations with trnH-psbA (rbcL+trnH-psbA, matK+ trnH-psbA, and rbcL+matK+trnH-psbA).
Universality of DNA Barcodes
Primer universality is an important criterion for an ideal DNA barcode (Kress and Erickson, 2007). In this regard, the rbcL, matK, and ITS had the best performance in PCR amplification and sequencing among the four regions (successfully amplifying and sequencing 100% samples), consistent with many previous studies (Xu et al., 2015; Yan et al., 2015). However, compared to the above three barcodes, trnH-psbA had a relatively low sequencing success rate of 95.12% when two primer pairs were used. This was due largely to poly (T) tracts at about 100 bp in the forward direction when sequencing. In addition, a 230 bp indel in three sequences of two species, i.e., O. intermedia King and Pantl and O. delicata resulted in alignment difficulties. As for the insertion events, small inversions associated with palindromes, and sequencing problems related to mononucleotide repeats within this non-coding chloroplast region will complicate its use as a barcode (Whitlock et al., 2010). Thus, sequence alignment of this region must be careful to avoid overestimates of the substitution events.
The Resolution of Tested Candidate Barcodes in Oberonia
When evaluated alone, the three plastid regions studied (rbcL, matK, trnH-psbA) had a resolution ranging from 21.43 to 57.14% (based on tree-building analysis), which is much lower than the discriminatory rate of nuclear region (Figure 1). Thus, all the single plastid regions are not recommended as DNA barcodes for the genus Oberonia. The low resolution of chloroplast regions has been previously reported in other plants including Paphiopedilum (25.74%), Curcuma (21.66%), and Quercus (0%) (Piredda et al., 2011; Chen et al., 2015; Guo et al., 2015). This could be due to the lower substitution rates that are found in plastid genomes, relative to their nuclear regions. Consequently, this highlights a need to search for nuclear DNA barcodes.
For the nuclear genome, the ITS and ITS2 generally provided better identification rates than the chloroplast sequences (Chen et al., 2010; Yao et al., 2010; Li et al., 2011). Likewise, the ITS and ITS2 have more parsimony informative sites, larger inter-specific distances and more discriminatory power than chloroplast regions in this study. In comparison the ITS distinguished 75% of monophyletic species, which was the best discrimination performance among the five loci (Figure 1). Meanwhile, any combinations with ITS produced better results than those combinations without (Figure 1 and Table 1). According to some previous results, it is difficult to amplify and directly sequence the region in some taxa because of incomplete concerted evolution of this multiple-copy region caused by hybridization or other factors (Alvarez and Wendel, 2003). However, it is not a problem in Oberonia and it has been widely used to generate phylogenies in many orchid taxa (Koehler et al., 2002; Zhai et al., 2014; Tang Y. et al., 2015). The amplification and sequencing rates of the ITS in this study were nearly 100%. Overall, for a single barcode, the ITS is the best candidate for Oberonia. Meanwhile, our results indicate that the ITS2 alone or in combinations with plastid markers did not have a higher discriminatory power than the ITS and/or its corresponding combinations (Figure 1 and Table 1). However, considering the ease of amplification, we suggest that the ITS2 may be an ideal supplementary barcode when ITS amplification has failed.
Multi-locus barcodes have been suggested as DNA barcodes for land plants and can often improve the resolution rate of species identification (CBOL Plant Working Group, 2009; Li et al., 2011). In this study, the two-locus barcode rbcL+matK, recommended by the CBOL Plant Working Group (2009), had a low discrimination rate of 64.29% based on the Tree-building method (Figure 1), which was lower than the identification rate of 72% proposed by the CBOL Plant Working Group (2009). One of the most plausible explanations for this discrepancy is that the CBOL Plant Working Group focused on assessing the relative, rather than the absolute discriminatory power of the tested barcodes. We sampled more closely related species within the Oberonia genus and while rbcL and matK discriminates among genera well, the resolution rates of these two markers, alone and in combination, decreased at infrageneric levels, especially within recently evolved genera. Of the 2- to 4- combinations of the five loci tested, rbcL+ITS and matK+ITS exhibited the best discriminatory performance, almost the same as the single ITS barcode (Figure 2, Figures S2, S3). In plant DNA barcoding studies, the use of markers from different genomes with different modes of inheritance has been suggested, because such combinations of DNA markers will further our understanding of species delimitation and the evolutionary processes of speciation (Li et al., 2011). Although the resolution rates of rbcL+ITS and matK+ITS were not better than the single marker ITS itself, we suggest that the combination marker either rbcL+ITS or matK+ITS, should be the first choice to barcode Oberonia species.
Figure 2. Bayesian inference (BI) tree based on the combination of matK+ITS sequences in Oberonia. Numbers on branches represent BI and NJ support values, respectively.
Implications of DNA Barcoding for the Current Taxonomy of Oberonia
The results obtained in this study shed some light on the identification and taxonomy of the genus Oberonia. In previous DNA barcoding studies on a single genus or between closely related groups, species misidentification can be corrected with DNA barcoding (Pryer et al., 2010; Zhang et al., 2012; Yan et al., 2015). In the present study, the sample O. kanburiensis 1 was initially identified as O. acaulis Griff., based on leaf morphology (lacked flowers). DNA barcoding showed that the sequences of this sample were different from the other samples of this species, and it continually clustered with O. kanburiensis Seidenf. sequences. Confusion often occurs between O. acaulis and O. kanburiensis, because their morphology is similar prior to flowering. After blooming we rechecked the specimen and confirmed that this accession was misidentified and was in-fact O. kanburiensis. Such misidentifications were also found in O. delacourii 1 and O. delacourii 2 which were initially identified as O. ensiformis (Smith) Lindl. based on leaf morphology (lacked flowers). This finding indicates that DNA barcoding can differentiate species, with small sample input, great speed, and higher reliability, than previous methods. Thus, a more robust method has been demonstrated for endangered species monitoring in Oberonia genus.
In addition to the effectiveness of correcting misidentified specimen, DNA barcoding could also help in resolving taxonomic uncertainties. For example, O. austro-yunnanensis S. C. Chen et Z. H. Tsi and O. jenkinsiana Griff. ex Lindl. were considered as two separate species, discriminated by a joint at the base of the leaf (O. austro-yunnanensis with a basal joint vs. O. jenkinsiana without a basal joint) and subtle differences in the lip of flowers. After carefully checking the specimen and the original literature, we found that neither O. austro-yunnanensis nor O. jenkinsiana have joints at the base of the leaf, and are morphologically similar in most respects. We suspect that these are the same species. In the phylogenetic tree (Figure 2) O. austro-yunnanensis continually clustered with O. jenkinsiana. Besides, O. jenkinsiana is a widespread species and the geographical distribution ranges of these two species overlap. Therefore, our analysis implies that species with identical morphology, distribution, and sequences should be treated as one species. It is also probable that O. sinica and O. insularis reflect a similar situation. In another situation, Oberonia acaulis var. luchunensis are treated as a variety of O. acaulis, but in the light of these DNA barcoding results and the large difference in morphology, O. acaulis var. luchunensis could be raised to an independent species rather than a variety of O. acaulis, although further work is necessary (Figure 2).
Discovery of new or cryptic species is an important application of DNA barcoding within taxonomy (DeSalle et al., 2005). Many studies have employed the DNA barcoding to discover new or cryptic species in a broad range of animals and plants (Burns et al., 2008; Liu et al., 2011; Saitoh et al., 2015). In our study, some taxa such as samples Oberonia sp. 1, 2, 3 are similar to Oberonia caulescens Lindl. with subtle morphological differences in flowers and leaves, not withholding geographical divergences. However, they did not cluster with O. caulescens in the phylogenetic trees, indicating the existence of a new or cryptic species, yet additional study is necessary.
Morphological characterization associated with geographical, ecological, reproductive, and molecular data will facilitate the construction of a robust taxonomic system for any particular taxa. However, taxonomy and species delimitation within a single genus, especially genera with closely related species, is more difficult. In our study, despite the excellent performance of DNA barcoding in the Oberonia species from China, DNA barcoding does have difficulties in discriminating closely related species. For example, taxonomic complex group O. arisanensis, O. delicata, O. japonica, and O. menghaiensis are morphological similar with subtle differences in the diagnostic characteristic of their flowers. DNA barcoding did not discriminate these four species and they all clustered together in the rbcL+ITS tree and matK+ITS tree (Figure 2, Figures S2, S3). Another case for DNA barcoding failure occurred in O. anthropophora, which displayed incongruent signals of nuclear and plastid gene regions (Figures S2, S3). The slow rate of molecular evolution, paralogy, incidence of hybridization, introgression, and incomplete sorting of ancestral polymorphisms are the most likely sources of DNA barcoding failure for closely related Oberonia species (Funk and Omland, 2003; Hollingsworth et al., 2011; León-Romero et al., 2012). The exploration of more molecular markers, such as SNP and SSR, are needed to develop DNA barcodes to assist in species identification (Liu et al., 2008; Yuan et al., 2012; Zeng et al., 2012). This data will promote progress in DNA barcoding, while also facilitate the identification of endangered species.
Designed the study: YL, YT, and FX. Performed the experiment and analyzed the data: YL and YT. Wrote the paper: YL, YT, and FX. All authors read and approved the final manuscript.
This study was financially supported by the National Natural Science Foundation of China (Grant no. 31370231).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are grateful to Prof. Hong-Feng Chen, Prof. Qi-Fei Yi, Dr. Huai-Zhen Tian, Prof. Yan-Liu, De-Ping Ye, Jian-Wu Li, Qiang Liu, Ming-Zhong Huang, and Ya-Ping Chen for their assistance in the sample collection. We also thank Dr. Hai-Fei Yan and Yu Zhang (South China Botanical Garden, Chinese Academy of Sciences) for their help in experiments, data analysis and revising the manuscripts. We thank Dr. Lei Duan and Dr. Juan Chen for revising the manuscripts.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.01791/full#supplementary-material
Burns, J. M., Janzen, D. H., Hajibabaei, M., Hallwachs, W., and Hebert, P. D. N. (2008). DNA barcodes and cryptic species of skipper butterflies in the genus Perichares in Area de Conservación Guanacaste, Costa Rica. Proc. Natl. Acad. Sci. U.S.A. 105, 6350–6355. doi: 10.1073/pnas.0712181105
Chen, J., Zhao, J. T., Erickson, D. L., Xia, N., and Kress, W. J. (2015). Testing DNA barcodes in closely related species of Curcuma (Zingiberaceae) from Myanmar and China. Mol. Ecol. Resour. 15, 337–348. doi: 10.1111/1755-0998.12319
Chen, S. C., Ormerod, P., and Wood, J. J. (2009). “Oberonia Lindl.,” in Flora of China, Vol. 25, eds Z. Y. Wu, P. H. Raven, and D. Y. Hong (Beijing; St. Louis, MO: Science Press and Missouri Botanical Garden Press), 236–245.
Chen, S., Yao, H., Han, J., Liu, C., Song, J., Shi, L., et al. (2010). Validation of the ITS2 region as a novel DNA barcode for identifying medicinal plant species. PLoS ONE 5:e8613. doi: 10.1371/journal.pone.0008613
Cuénoud, P., Savolainen, V., Chatrou, L. W., Powell, M., Grayer, R. J., and Chase, M. W. (2002). Molecular phylogenetics of Caryophyllales based on nuclear 18S rDNA and plastid rbcL, atpB, and matK DNA sequences. Am. J. Bot. 89, 132–144. doi: 10.3732/ajb.89.1.132
Funk, D. J., and Omland, K. E. (2003). Species-level paraphyly and polyphyly: frequency, causes, and consequences, with insights from animal mitochondrial DNA. Annu. Rev. Ecol. Evol. Syst. 34, 397–423. doi: 10.1146/annurev.ecolsys.34.011802.132421
Gathier, G., van der Niet, T., Peelen, T., van Vugt, R. R., Eurlings, M. C., and Gravendeel, B. (2013). Forensic identification of CITES protected slimming cactus (Hoodia) using DNA barcoding. J. Forensic Sci. 58, 1467–1471. doi: 10.1111/1556-4029.12184
Hebert, P. D. N., Penton, E. H., Burns, J. M., Janzen, D. H., and Hallwachs, W. (2004). Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proc. Natl. Acad. Sci. U.S.A. 101, 14812–14817. doi: 10.1073/pnas.0406166101
Huang, J., Zhang, A., Mao, S., and Huang, Y. (2013). DNA Barcoding and species boundary delimitation of selected species of Chinese Acridoidea (Orthoptera: Caelifera). PLoS ONE 8:e82400. doi: 10.1371/journal.pone.0082400
Ji, Y., Ashton, L., Pedley, S. M., Edwards, D. P., Tang, Y., Nakamura, A., et al. (2013). Reliable, verifiable and efficient monitoring of biodiversity via metabarcoding. Ecol. Lett. 16, 1245–1257. doi: 10.1111/ele.12162
Koehler, S., Williams, N. H., Whitten, W. M., and do Amaral, M. C. E. (2002). Phylogeny of the Bifrenaria (Orchidaceae) complex based on morphology and sequence data from nuclear rDNA internal transcribed spacers (ITS) and chloroplast trnL-trnF region. Int. J. Plant Sci. 163, 1055–1066. doi: 10.1086/342035
Kress, W. J., and Erickson, D. L. (2007). A two-locus global DNA barcode for land plants, the coding rbcL gene complements the non-coding trnH-psbA spacer region. PLoS ONE 2:e508. doi: 10.1371/journal.pone.0000508
León-Romero, Y., Mejía, O., and Soto-Galera, E. (2012). DNA barcoding reveals taxonomic conflicts in the Herichthys bartoni species group (Pisces: Cichlidae). Mol. Ecol. Resour. 12, 1021–1026. doi: 10.1111/1755-0998.12018
Li, D. Z., Gao, L. M., Li, H. T., Wang, H., Ge, X. J., and Liu, J. Q. (2011). Comparative analysis of a large dataset indicates that internal transcribed spacer (ITS) should be incorporated into the core barcode for seed plants. Proc. Natl. Acad. Sci. U.S.A. 108, 19641–19646. doi: 10.1073/pnas.1104551108
Liu, J., Möller, M., Gao, L. M., Zhang, D. Q., and Li, D. Z. (2011). DNA barcoding for the discrimination of Eurasian yews (Taxus L., Taxaceae) and the discovery of cryptic species. Mol. Ecol. Resour. 11, 89–100. doi: 10.1111/j.1755-0998.2010.02907.x
Liu, J., Yan, H. F., Newmaster, S. G., Pei, N., Ragupathy, S., and Ge, X. J. (2014). The use of DNA barcoding as a tool for the conservation biogeography of subtropical forests in China. Divers. Distrib. 21, 188–199. doi: 10.1111/ddi.12276
Liu, Z., Meng, J., and Sun, X. (2008). A novel feature-based method for whole genome phylogenetic analysis without alignment: application to HEV genotyping and subtyping. Biochem. Biophys. Res. Commun. 368, 223–230. doi: 10.1016/j.bbrc.2008.01.070
Meier, R., Shiyang, K., Vaidya, G., and Ng, P. K. (2006). DNA barcoding and taxonomy in Diptera: a tale of high intraspecific variability and low identifications success. Syst. Biol. 55, 715–728. doi: 10.1080/10635150600969864
Mutanen, M., Kekkonen, M., Prosser, S. W., Hebert, P. D., and Kaila, L. (2015). One species in eight: DNA barcodes from type specimens resolve a taxonomic quagmire. Mol. Ecol. Resour. 15, 967–984. doi: 10.1111/1755-0998.12361
Piredda, R., Simeone, M. C., Attimonelli, M., Bellarosa, R., and Schirone, B. (2011). Prospects of barcoding the Italian wild dendroflora: oaks reveal severe limitations to tracking species identity. Mol. Ecol. Resour. 11, 72–83. doi: 10.1111/j.1755-0998.2010.02900.x
Pryer, K. M., Schuettpelz, E., Huiet, L., Grusz, A. L., Rothfels, C. J., Avent, T., et al. (2010). DNA barcoding exposes a case of mistaken identity in the fern horticultural trade. Mol. Ecol. Resour. 10, 979–985. doi: 10.1111/j.1755-0998.2010.02858.x
Reeves, G., Chase, M. W., Goldblatt, P., Rudall, P., Fay, M. F., Cox, A. V., et al. (2001). Molecular systematics of Iridaceae: evidence from four plastid DNA regions. Am. J. Bot. 88, 2074–2087. doi: 10.2307/3558433
Ren, B. Q., Xiang, X. G., and Chen, Z. D. (2010). Species identification of Alnus (Betulaceae) using nrDNA and cpDNA genetic markers. Mol. Ecol. Resour. 10, 594–605. doi: 10.1111/j.1755-0998.2009.02815.x
Saitoh, T., Sugita, N., Someya, S., Iwami, Y., Kobayashi, S., Kamigaichi, H., et al. (2015). DNA barcoding reveals 24 distinct lineages as cryptic bird species candidates in and around the Japanese Archipelago. Mol. Ecol. Resour. 15, 177–186. doi: 10.1111/1755-0998.12282
Sun, Y., Skinner, D. Z., Liang, G. H., and Hulbert, S. H. (1994). Phylogenetic analysis of Sorghum and related taxa using internal transcribed spacers of nuclear ribosomal DNA. Theor. Appl. Genet. 89, 26–32. doi: 10.1007/BF00226978
Taberlet, P., Coissac, E., Pompanon, F., Brochmann, C., and Willerslev, E. (2012). Towards next-generation biodiversity assessment using DNA metabarcoding. Mol. Ecol. 21, 2045–2050. doi: 10.1111/j.1365-294X.2012.05470.x
Tang, G. D., Zhang, G. Q., Hong, W. J., Liu, Z. J., and Zhuang, X. Y. (2015). Phylogenetic analysis of Malaxideae (Orchidaceae: Epidendroideae): two new species based on the combined nrDNA ITS and chloroplast matK sequences. Guihaia 35, 447–463. doi: 10.11931/guihaia.gxzw201506015
Tang, Y., Yukawa, T., Bateman, R. M., Jiang, H., and Peng, H. (2015). Phylogeny and classification of the East Asian Amitostigma alliance (Orchidaceae: Orchideae) based on six DNA markers. BMC Evol. Biol. 15:96. doi: 10.1186/s12862-015-0376-3
Xu, S., Li, D., Li, J. W., Xiang, X., Jin, W., Huang, W., et al. (2015). Evaluation of the DNA barcodes in Dendrobium (Orchidaceae) from mainland Asia. PLoS ONE 10:e0115168. doi: 10.1371/journal.pone.0115168
Yan, L. J., Liu, J., Möller, M., Zhang, L., Zhang, X. M., Li, D. Z., et al. (2015). DNA barcoding of Rhododendron (Ericaceae), the largest Chinese plant genus in biodiversity hotspots of the Himalaya–Hengduan Mountains. Mol. Ecol. Resour. 15, 932–944. doi: 10.1111/1755-0998.12353
Yuan, Z., Chu, G., Dan, Y., Li, J., Zhang, L., Gao, X., et al. (2012). BRCA1: a new candidate gene for bovine mastitis and its association analysis between single nucleotide polymorphisms and milk somatic cell score. Mol. Biol. Rep. 39, 6625–6631. doi: 10.1007/s11033-012-1467-5
Zeng, X., Yuan, Z., Tong, X., Li, Q., Gao, W., Qin, M., et al. (2012). Phylogenetic study of Oryzoideae species and related taxa of the Poaceae based on atpB-rbcL and ndhF DNA sequences. Mol. Biol. Rep. 39, 5737–5744. doi: 10.1007/s11033-011-1383-0
Zhai, J. W., Zhang, G. Q., Li, L., Wang, M., Chen, L. J., Chung, S. W., et al. (2014). A new phylogenetic analysis sheds new light on the relationships in the Calanthe alliance (Orchidaceae) in China. Mol. Phylogenet. Evol. 77, 216–222. doi: 10.1016/j.ympev.2014.04.005
Keywords: DNA barcoding, Oberonia, taxonomy, radiation, species identification, ITS
Citation: Li YL, Tong Y and Xing FW (2016) DNA Barcoding Evaluation and Its Taxonomic Implications in the Recently Evolved Genus Oberonia Lindl. (Orchidaceae) in China. Front. Plant Sci. 7:1791. doi: 10.3389/fpls.2016.01791
Received: 23 September 2016; Accepted: 14 November 2016;
Published: 05 December 2016.
Edited by:Renchao Zhou, Sun Yat-sen University, China
Reviewed by:János Taller, University of Pannonia, Hungary
Łukasz Kajtoch, Institute of Systematics and Evolution of Animals (PAN), Poland
Junwen Zhai, College of Landscape Architecture, China
Copyright © 2016 Li, Tong and Xing. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Fuwu Xing, email@example.com