DNA barcoding ofRhododendron(Ericaceae), the largest Chinese plant genus in biodiversity hotspots of the Himalaya-Hengduan Mountains

The first complete mitochondrial genome of the genus Laelaps with novel gene arrangement reveals extensive rearrangement and phylogenetics in the superfamily Dermanyssoidea

We collected 56 specimens of Laelaps chini from the endemic Hengduan Mountain rat species (Eothenomys miletus) and obtained the first complete mitochondrial genome of L. chini by next-generation sequencing (NGS). The L. chini mitogenome is 16,507 bp in size and contains 37 genes and a control region of 2380 bp in length. The L. chini mitogenome has a high AT content and a compact arrangement with four overlapping regions ranging from 1 to 2 bp and 16 spacer regions ranging from 1 to 48 bp. We analyzed 13 protein-coding genes of L. chini mitogenome  and found that protein-coding genes in the L. chini mitogenome preferred codons ending in A/U and codon usage pattern was mainly influenced by natural selection. Cox1 has the slowest evolution rate and cox3 has the fastest evolution rate. We combined the mitochondrial genome of eight species of gamasid mites in the superfamily Dermanyssoidea from Genbank and the L. chini mitochondrial genome to analyze its rearrangement patterns and breakpoint numbers. We found that the L. chini mitogenome showed a novel arrangement pattern and nine species of gamasid mites in the superfamily Dermanyssoidea, which have been sequenced complete mitochondrial genomes to date, all showed different degrees of rearrangement. Laelaps chini, Echinolaelaps echidninus and Echinolaelaps fukinenensis were closely related species based on genetic distance and phylogenetic analyses. Notably they are clustered with Varroa destructor of the family Varroidae, suggesting that the family Varroidae is more closely related to the family Laelapidae, but more data are needed to test whether Varroa can be classified under the family Laelapidae. The L. chini mitogenome is the first complete mitochondrial genome for the genus Laelaps, and contributes to further exploration of the mitochondrial gene rearrangements and phylogeny for the superfamily Dermanyssoidea.

Unlocking the Genetic Identity of Endangered Paphiopedilum Orchids: A DNA Barcoding Approach

Orchids of the genus Paphiopedilum, also called slippers, are among the most valued representatives of the Orchidaceae family due to their aesthetic qualities. Due to overexploitation, deforestation, and illegal trade in these plants, especially in the vegetative phase, Paphiopedilum requires special protection. This genus is listed in Appendix I of the Convention on International Trade in Endangered Species of Wild Fauna and Flora. Their precise identification is of great importance for the preservation of genetic resources and biodiversity of the orchid family (Orchidaceae). Therefore, the main objective of the study was to investigate the usefulness of the DNA barcoding technique for the identification of endangered orchids of the genus Paphiopedilum and to determine the effectiveness of five loci: matK, rbcL, ITS2, atpF-atpH and trnH-psbA as potential molecular markers for species of this genus. Among single locus barcodes, matK was the most effective at identifying species (64%). Furthermore, matK, ITS2, matK + rbcL, and matK + trnH-psbA barcodes can be successfully used as a complementary tool to identify Paphiopedilum orchids while supporting morphological data provided by taxonomists.

DNA Barcoding of Invasive Terrestrial Plant Species in India

Invasive plants are known to cause biodiversity loss and pose a major risk to human health and environment. Identification of invasive plants and distinguishing them from native species has been relied on morphological examination. Stringent requirement of floral characters and decreasing number of expert taxonomists are making conventional morphology-based identification system tedious and resource-intensive. DNA barcoding may help in quick identification of invasive species if distinct sequence divergence pattern at various taxonomic levels is observed. The present work evaluates the utility of four molecular markers; rbcL, matK, their combination (rbcL + matK), and psbA-trnH for identification of 37 invasive plant species from India and also in distinguishing them from 97 native species. A psbA-trnH locus was found to be of restricted utility in this work as it was represented by the members of a single family. A hierarchical increase in K2P mean divergence across different taxonomic levels was found to be the maximum for matK alone followed by rbcL + matK and rbcL alone, respectively. NJ clustering analysis, however, confirmed the suitability of combined locus (rbcL + matK) over individual rbcL and matK as the DNA barcode. RbcL showed the lowest resolution power among the three markers studied. MatK exhibited much better performance compared to rbcL alone in identifying most of the species accurately although it failed to show monophyly of genus Dinebra. Two families; Asteraceae and Poaceae, remained polyphyletic in the trees constructed by all three markers. Combined locus (rbcL + matK) was found to be the most suitable marker as it raised the resolution power of both the markers and could identify more than 90% of genera correctly. Phylogenetic tree constructed by Maximum-Parsimony method using combined locus as a molecular marker exhibited the best resolution, thus, supporting the significance of two-locus combination of rbcL + matK for barcoding invasive plant species from India. Present study contributes to the global barcode data of invasive plant species by adding fifty-one new sequences to it. Effective barcoding of additional number of native as well as invasive plant species from India is possible using this dual locus if it is combined with one or more new molecular plastid markers. Expansion of barcode database with a focus on barcode performance optimisation to improve discrimination ability at species level can be undertaken in future.

Exploring DNA barcode for accurate identification of threatened Aconitum L. species from Western Himalaya

BACKGROUND

Aconitum species, belonging to Ranunculaceae, have high medicinal importance but due to their overexploitation come under IUCN (International Union for Conservation of Nature) red list. The precise identification of the Aconitum species is equally important because they are used in herbal formulations. The present study aimed to develop an efficient DNA barcode system for the authentic identification of Aconitum species.

METHODS AND RESULTS

A set of 92 barcode gene sequences (including 12 developed during the present study and 80 retrieved from NCBI) of 5 Aconitum species (A. heterophyllum, A. vialoceum, A. japonicum, A. napellus, and A. stapfianum) were analyzed using three methods (tree-based, distance-based, and similarity-based) for species discrimination. The PWG-distance method was found most effective for species discrimination. The discrimination rate of PWG- distance ranged from 33.3% (rbcL + trnH-psbA) to 100% (ITS, rbcL + ITS, ITS + trnH-psbA and rbcL + ITS + trnH-psbA). Among DNA barcodes and their combinations, the ITS marker had the highest degree of species discrimination (NJ-40%, PWG-100% and BLAST-40%), followed by trnH-psbA (NJ-20%, PWG-60% and BLAST-20%). ITS also had higher barcoding gap as compared to other individual barcodes and their combinations. Further, we also analyzed six Aconitum species (A. balfourii, A. ferox, A. heterophyllum, A. rotundifolium, A. soongaricum and A. violaceum) existing in Western Himalaya. These species were distinguished clearly through tree-based method using the ITS barcode gene with 100% species resolution.

CONCLUSION

ITS showed the best species discrimination power and was used to develop species-specific barcodes for Aconitum species. DNA barcodes developed during the present study can be used to identify Aconitum species.

Plastome phylogenomics and morphological traits analyses provide new insights into the phylogenetic position, species delimitation and speciation of Triplostegia (Caprifoliaceae)

BACKGROUND

The genus Triplostegia contains two recognized species, T. glandulifera and T. grandiflora, but its phylogenetic position and species delimitation remain controversial. In this study, we assembled plastid genomes and nuclear ribosomal DNA (nrDNA) cistrons sampled from 22 wild Triplostegia individuals, each from a separate population, and examined these with 11 recently published Triplostegia plastomes. Morphological traits were measured from herbarium specimens and wild material, and ecological niche models were constructed.

RESULTS

Triplostegia is a monophyletic genus within the subfamily Dipsacoideae comprising three monophyletic species, T. glandulifera, T. grandiflora, and an unrecognized species Triplostegia sp. A, which occupies much higher altitude than the other two. The new species had previously been misidentified as T. glandulifera, but differs in taproot, leaf, and other characters. Triplotegia is an old genus, with stem age 39.96 Ma, and within it T. glandulifera diverged 7.94 Ma. Triplostegia grandiflora and sp. A diverged 1.05 Ma, perhaps in response to Quaternary climate fluctuations. Niche overlap between Triplostegia species was positively correlated with their phylogenetic relatedness.

CONCLUSIONS

Our results provide new insights into the species delimitation of Triplostegia, and indicate that a taxonomic revision of Triplostegia is needed. We also identified that either rpoB-trnC or ycf1 could serve as a DNA barcode for Triplostegia.

Genome-wide identification and bioinformatics analysis of the WD40 transcription factor family and candidate gene screening for anthocyanin biosynthesis in Rhododendron simsii

WD40 transcription factors (TFs) constitute a large gene family in eukaryotes, playing diverse roles in cellular processes. However, their functions in the major ornamental plant, Rhododendron simsii, remain poorly understood. In this study, we identified 258 WD40 proteins in the R. simsii genome, which exhibited an uneven distribution across chromosomes. Based on domain compositions and phylogenetic analysis, we classified these 258 RsWD40 proteins into 42 subfamilies and 47 clusters. Comparative genomic analysis suggested that the expansion of the WD40 gene family predates the divergence of green algae and higher plants, indicating an ancient origin. Furthermore, by analyzing the duplication patterns of RsWD40 genes, we found that transposed duplication played a major role in their expansion. Notably, the majority of RsWD40 gene duplication pairs underwent purifying selection during evolution. Synteny analysis identified significant orthologous gene pairs between R. simsii and Arabidopsis thaliana, Oryza sativa, Vitis vinifera, and Malus domestica. We also investigated potential candidate genes involved in anthocyanin biosynthesis during different flower development stages in R. simsii using RNA-seq data. Specifically, we identified 10 candidate genes during the bud stage and 7 candidate genes during the full bloom stage. GO enrichment analysis of these candidate genes revealed the potential involvement of the ubiquitination process in anthocyanin biosynthesis. Overall, our findings provide a valuable foundation for further investigation and functional analysis of WD40 genes, as well as research on the molecular mechanisms underlying anthocyanin biosynthesis in Rhododendron species.

A haplotype-resolved genome assembly of Rhododendron vialii based on PacBio HiFi reads and Hi-C data

Rhododendron vialii (subgen. Azaleastrum) is an evergreen shrub with high ornamental value. This species has been listed as a plant species with extremely small populations (PSESP) for urgent protection by China's Yunnan provincial government in 2021, due to anthropogenic habitat fragmentation. However, limited genomic resources hinder scientifically understanding of genetic threats that the species is currently facing. In this study, we assembled a high-quality haplotype-resolved genome of R. vialii based on PacBio HiFi long reads and Hi-C reads. The assembly contains two haploid genomes with sizes 532.73 Mb and 521.98 Mb, with contig N50 length of 35.67 Mb and 34.70 Mb, respectively. About 99.92% of the assembled sequences could be anchored to 26 pseudochromosomes, and 14 gapless assembled chromosomes were included in this assembly. Additionally, 60,926 protein-coding genes were identified, of which 93.82% were functionally annotated. This is the first reported genome of R. vialii, and hopefully it will lay the foundations for further research into the conservation genomics and horticultural domestication of this ornamentally important species.

Ploidy variation in Rhododendron subsection Maddenia and its implications for conservation

Polyploidy, which is common in plants, can confound taxon recognition and hence conservation assessments. In the taxonomically complex genus Rhododendron, 25 % of the over 1,300 taxa are considered under threat and 27 % Near Threatened or Data Deficient, with their taxonomy needing to be resolved urgently. Although ploidy levels of Rhododendron taxa range from diploid (2x) to dodecaploid (12x) according to previous reports, the extent of polyploidy across the genus has not been examined. We first summarized the taxonomic distribution of polyploids in the genus based on the literature. Then as a case study, we estimated ploidy levels of 47 taxa in subsection Maddenia (subgenus Rhododendron, section Rhododendron) using flow cytometry, together with verification of meiotic chromosome counts for representative taxa. The summary of reported ploidy in Rhododendron indicates that polyploidy is most common in subgenera Pentanthera and Rhododendron. In subsection Maddenia, all examined taxa are diploids except for the R. maddenii complex that shows a high ploidy variation (2-8x, 12x). We investigated ploidy level of 12 taxa in subsection Maddenia for the first time, and estimated genome sizes of two Rhododendron species. Knowledge of ploidy levels will inform phylogenetic analysis of unresolved species complexes. Overall, our study of subsection Maddenia provides a model for examining multiple issues including taxonomic complexity, ploidy variation and geographic distribution in relation to biodiversity conservation.

Flora diversity survey and establishment of a plant DNA barcode database of Lomas ecosystems in Peru

Lomas formations or "fog oases" are islands of vegetation in the desert belt of the west coast of South America, with a unique vegetation composition among the world's deserts. However, plant diversity and conservation studies have long been neglected, and there exists a severe gap in plant DNA sequence information. To address the lack of DNA information, we conducted field collections and laboratory DNA sequencing to establish a DNA barcode reference library of Lomas plants from Peru. This database provides 1,207 plant specimens and 3,129 DNA barcodes data corresponding with collections from 16 Lomas locations in Peru, during 2017 and 2018. This database will facilitate both rapid species identification and basic studies on plant diversity, thereby enhancing our understanding of Lomas flora's composition and temporal variation, and providing valuable resources for conserving plant diversity and maintaining the stability of the fragile Lomas ecosystems.

DNA barcoding in herbal medicine: Retrospective and prospective

DNA barcoding has been widely used for herb identification in recent decades, enabling safety and innovation in the field of herbal medicine. In this article, we summarize recent progress in DNA barcoding for herbal medicine to provide ideas for the further development and application of this technology. Most importantly, the standard DNA barcode has been extended in two ways. First, while conventional DNA barcodes have been widely promoted for their versatility in the identification of fresh or well-preserved samples, super-barcodes based on plastid genomes have rapidly developed and have shown advantages in species identification at low taxonomic levels. Second, mini-barcodes are attractive because they perform better in cases of degraded DNA from herbal materials. In addition, some molecular techniques, such as high-throughput sequencing and isothermal amplification, are combined with DNA barcodes for species identification, which has expanded the applications of herb identification based on DNA barcoding and brought about the post-DNA-barcoding era. Furthermore, standard and high-species coverage DNA barcode reference libraries have been constructed to provide reference sequences for species identification, which increases the accuracy and credibility of species discrimination based on DNA barcodes. In summary, DNA barcoding should play a key role in the quality control of traditional herbal medicine and in the international herb trade.

Variation in gene expression along an elevation gradient of Rhododendron sanguineum var. haemaleum assessed in a comparative transcriptomic analysis

Selection along environmental gradients may play a vital role in driving adaptive evolution. Nevertheless, genomic variation and genetic adaptation along environmental clines remains largely unknown in plants in alpine ecosystems. To close this knowledge gap, we assayed transcriptomic profiles of late flower bud and early leaf bud of Rhododendron sanguineum var. haemaleum from four different elevational belts between 3,000 m and 3,800 m in the Gaoligong Mountains. By comparing differences in gene expression of these samples, a gene co-expression network (WGCNA) was constructed to identify candidate genes related to elevation. We found that the overall gene expression patterns are organ-specific for the flower and leaf. Differentially expressed unigenes were identified in these organs. In flowers, these were mainly related to terpenoid metabolism (RsHMGR, RsTPS), while in leaves mainly related to anthocyanin biosynthesis (RsCHS, RsF3'5'H). Terpenoids are the main components of flower scent (fragrance) likely attracting insects for pollination. In response to fewer pollinators at higher elevation zone, it seems relatively less scent is produced in flower organs to reduce energy consumption. Secondary metabolites in leaves such as anthocyanins determine the plants' alternative adaptive strategy to extreme environments, such as selective pressures of insect herbivory from environmental changes and substrate competition in biosynthesis pathways at high elevations. Our findings indicated that the gene expression profiles generated from flower and leaf organs showed parallel expression shifts but with different functionality, suggesting the existence of flexibility in response strategies of plants exposed to heterogeneous environments across elevational gradients. The genes identified here are likely to be involved in the adaptation of the plants to these varying mountainous environments. This study thus contributes to our understanding of the molecular mechanisms of adaptation in response to environmental change.

Molecular characterization of a chalcone synthase gene RhCHS from Rhododendron × hybridum Hort

Chalcone synthase (CHS) plays a vital role in anthocyanin biosynthesis pathway, which is associated with petal color of flower. To date, lots of CHS genes have been obtained from plants, while few were from Rhododendron genus. In this study we got a new CHS gene named RhCHS (MW358095) from Rhododendron × hybridum Hort. It had a 2040 bp coding region consisting of two exons and one intron. By using the deduced RhCHS protein as a query sequence, 15 CHS homologous family genes with sequence similarity from 60% to 98%, designated as RgCHS-D(x), were retrieved from the genome assembly of Rhododendron griersonianum (RGv1.1) by TBlastN. 12 CHS family genes were found locating in No.9 chromosome arranged in clusters, while only 3 of them exhibited in No.1, 2, and 8 chromosomes, respectively. The results revealed gene duplication of CHS in evolutionary process. Multiple alignment of the deduced amino acid sequence of RhCHS showed high similarity of the active site, the catalytic residue, and the signature motif, the conserved characteristics of which were also exhibited in the tertiary structure prediction of the RhCHS, as well as the phylogenetic tree, all these demonstrated the RhCHS belonging to the type III PKS superfamily. HPLC-MS/MS of flower petals detected the total concentration of CC, DC, and PelC. These anthocyanidins showed an overall increasing trend during the flowering period and reached the peak in the full-blooming stage, which was consistence with the changeable rule of RhCHS expression level. The promoter, which was 1507 bp exhibiting high β-glucuronidase (GUS) staining activity, was predicted containing many cis-acting elements, especially light and transcription factor such as bHLH, MYB, WRKY, Dof, and ERF. In short, this study may provide the help to Rhododendron × hybridum Hort. not only in the mechanism research of petals color exhibition, but also in molecular breeding of CHS practice value.

Gapless genome assembly of azalea and multi-omics investigation into divergence between two species with distinct flower color

The genus Rhododendron (Ericaceae), with more than 1000 species highly diverse in flower color, is providing distinct ornamental values and a model system for flower color studies. Here, we investigated the divergence between two parental species with different flower color widely used for azalea breeding. Gapless genome assembly was generated for the yellow-flowered azalea, Rhododendron molle. Comparative genomics found recent proliferation of long terminal repeat retrotransposons (LTR-RTs), especially Gypsy, has resulted in a 125 Mb (19%) genome size increase in species-specific regions, and a significant amount of dispersed gene duplicates (13 402) and pseudogenes (17 437). Metabolomic assessment revealed that yellow flower coloration is attributed to the dynamic changes of carotenoids/flavonols biosynthesis and chlorophyll degradation. Time-ordered gene co-expression networks (TO-GCNs) and the comparison confirmed the metabolome and uncovered the specific gene regulatory changes underpinning the distinct flower pigmentation. B3 and ERF TFs were found dominating the gene regulation of carotenoids/flavonols characterized pigmentation in R. molle, while WRKY, ERF, WD40, C2H2, and NAC TFs collectively regulated the anthocyanins characterized pigmentation in the red-flowered R simsii. This study employed a multi-omics strategy in disentangling the complex divergence between two important azaleas and provided references for further functional genetics and molecular breeding.

Investigation of insecticidal activity of two Rhododendron species on stored-product insects

This study was designed to investigate the insecticidal activity of the essential oils (EOs) and extracts from Rhododendron rufum and Rhododendron przewalskii. The EOs were extracted from the leaves of R. Rufum and R. przewalskii by hydro-distillation and their chemical components were analyzed by gas chromatography-mass spectrometry (GC-MS). The repellency, contact toxicity and antifeedant activity of the EOs and extracts were evaluated against Sitophilus oryzae and Tribolium castaneum along with those of their main components. A total of nine compounds were identified from the EO of R. Rufum, and the most abundant component was myristicin (79.72%). The EO of R. Rufum exhibited repellent activities at different levels and its main compound myristicin showed contact toxicity and repellent effects against S. oryzae and T. castaneum. Meanwhile, by bioassay-guided fractionation, four compounds with strong antifeedant activities against T. castaneum, 24-methylenecycloartanyl-2'E, 4'Z-tetradecadienoate (1), methyl thyrsiflorin B acetate (2), friedelin (3) and Excoecarin R1 methyl ester (4) were separated and identified from the ethanol extract of R. przewalskii for the first time. Considering the significant anti-insect activities, the EOs and extracts of R. Rufum and R. przewalskii might be used in integrated pest strategies, establishing a good perspective for the comprehensive use of natural plant resources of Rhododendron genus.

An integrated study of Violae Herba (Viola philippica) and five adulterants by morphology, chemical compositions and chloroplast genomes: insights into its certified plant origin

BACKGROUND

Viola philippica Cav. is the only original plant for Violae Herba, as described in the Chinese Pharmacopoeia. The quality of this crude drug is affected by several adulterants from congeneric Viola species, and the authentic plant origin of Violae Herba is still controversial. Genome-based identification offers abundant genetic information and potential molecular markers that can be used for the authentication of closely related species. This study aims to investigate the certified origin of Violae Herba and to develop more effective markers for these easily confused species at the genetic level.

METHODS

We compared the morphology and chemical composition of 18 batches of commercial samples and six widespread medicinal Viola plants used as Violae Herba or its substitutes by TLC and HPLC-Triple-TOF-MS/MS analyses. The complete chloroplast genomes of these species were sequenced and analyzed, including the general features, repeat sequences, mutational hotspots and phylogeny. The complete chloroplast genomes used as superbarcodes and some specific barcodes screened from mutational hotspots were tested for their ability to distinguish Viola species.

RESULTS

A comparative study showed that Violae Herba is a multi-origin traditional Chinese medicine. Commercial decoction pieces and the standard reference drug were mainly derived from V. prionantha, clashing with the record in the Chinese Pharmacopoeia. Chloroplast genome analyses of V. philippica and five adulterants indicated that sequence divergence was relatively low within Viola species. By tree-based approaches, the complete chloroplast genomes showed a better discrimination ability and phylogenetic resolution for each Viola species. These results indicate that the whole chloroplast genomes can be used as superbarcodes to differentiate Viola medicinal plants. More specific DNA barcodes could be further developed from the Viola chloroplast genomes for more efficient and rapid identification of commercial Violae Herba and its adulterants.

CONCLUSIONS

This study has implications for chloroplast genome-based phylogenetic analysis and the authentication of multiple Viola species used as Violae Herba. The legal origin recorded in the Chinese Pharmacopoeia should be further revised to V. prionantha, in line with the commercial Violae Herba in the TCM markets.

DNA Barcoding and Phylogenomic Analysis of the Genus Fritillaria in China Based on Complete Chloroplast Genomes

The Fritillaria is an extremely complicated genus in taxonomy and phylogeny, which contains numerous medicinal species in China. Both traditional characteristic-based taxonomy and universal DNA barcodes (ITS, trnH-psbA, and rbcL) are difficult to effectively identify the species. Here, we generated a large dataset of chloroplast genomes from multiple accessions per species of Fritillaria to evaluate their effectiveness in species discrimination. Moreover, phylogeny of species in China was explored based on the complete chloroplast genomes, and then divergence times of each node were estimated. The results showed that all 21 species in Fritillaria here (including two suspicious species) could be correctly discriminated using cpDNA genomes except F. cirrhosa, which suggested that DNA super-barcode could greatly enhance species discriminatory resolution for complicated genera. Furthermore, four regions (ycf1, matK-trnG-GCC, rpoC1, and matK) gained remarkably higher resolution than that of other plastid regions, but only matK might be suitable to identify Fritillaria species in consideration of its lengths. Phylogenomic analysis showed that the subgenus Fritillaria in China was divided into four major clades with obvious geographic structure. Among them, Clade I, mainly distributed in southwest China, was a young and complicated group. Moreover, according to the analysis, taxonomic treatments of the two suspicious species, namely "F. omeiensis" and "F. hupehensis" in Flora of China (2000) are questionable and might need further revision. Molecular dating revealed that both origin and divergence of subgenus Fritillaria, as well as its four major clades, were significantly associated with geological and climatic fluctuations during the Middle to Late Miocene. This study would enrich case studies of DNA super-barcode and provide new insights on speciation, lineage diversification, and biogeography of the Fritillaria in China.

Testing the Complete Plastome for Species Discrimination, Cryptic Species Discovery and Phylogenetic Resolution in Cephalotaxus (Cephalotaxaceae)

Species of Cephalotaxus have great economic and ecological values. However, the taxonomy and interspecific phylogenetic relationships within the genus have been controversial and remained not fully resolved until now. To date, no study examined the efficiency of the complete plastome as super-barcode across Cephalotaxus species with multiple samples per taxon. In this study, we have evaluated the complete plastome in species discrimination and phylogenetic resolution in Cephalotaxus by including 32 individuals of all eight recognized species and five varieties following Farjon's classification (2010) with multiple samples per taxon. Our results indicated that not all species recognized in recent taxonomic revisions of Cephalotaxus could be distinguished and not all were monophyletic. Based on the plastome phylogeny, a new taxonomic classification for the genus comprising nine species and two varieties, including a cryptic species, was proposed. The phylogeny also resolved all interspecific relationships. Compared to the plastome based classification, standard DNA barcodes, alone or in combination, only recognized a maximum of seven out of the nine species. Moreover, two highly variable single loci, ycf1 and rps16, each alone achieved full species discrimination. With the moderate length of 1079 bp, rps16 is proposed as a specific barcode to discriminate Cephalotaxus species. The super-barcodes and specific barcode candidates will aid in the identification of endangered Cephalotaxus species, and to help focus conservation measures.

DNA barcoding of medicinal orchids in Asia

Growing popularity of herbal medicine has increased the demand of medicinal orchids in the global markets leading to their overharvesting from natural habitats for illegal trade. To stop such illegal trade, the correct identification of orchid species from their traded products is a foremost requirement. Different species of medicinal orchids are traded as their dried or fresh parts (tubers, pseudobulbs, stems), which look similar to each other making it almost impossible to identify them merely based on morphological observation. To overcome this problem, DNA barcoding could be an important method for accurate identification of medicinal orchids. Therefore, this research evaluated DNA barcoding of medicinal orchids in Asia where illegal trade of medicinal orchids has long existed. Based on genetic distance, similarity-based and tree-based methods with sampling nearly 7,000 sequences from five single barcodes (ITS, ITS2, matK, rbcL, trnH-psbA and their seven combinations), this study revealed that DNA barcoding is effective for identifying medicinal orchids. Among single locus, ITS performed the best barcode, whereas ITS + matK exhibited the most efficient barcode among multi-loci. A barcode library as a resource for identifying medicinal orchids has been established which contains about 7,000 sequences of 380 species (i.e. 90%) of medicinal orchids in Asia.

Natural hybridization among three Rhododendron species (Ericaceae) revealed by morphological and genomic evidence

BACKGROUND

Natural hybridization can influence the adaptive response to selection and accelerate species diversification. Understanding the composition and structure of hybrid zones may elucidate patterns of hybridization processes that are important to the formation and maintenance of species, especially for taxa that have experienced rapidly adaptive radiation. Here, we used morphological traits, ddRAD-seq and plastid DNA sequence data to investigate the structure of a Rhododendron hybrid zone and uncover the hybridization patterns among three sympatric and closely related species.

RESULTS

Our results show that the hybrid zone is complex, where bi-directional hybridization takes place among the three sympatric parental species: R. spinuliferum, R. scabrifolium, and R. spiciferum. Hybrids between R. spinuliferum and R. spiciferum (R. ×duclouxii) comprise multiple hybrid classes and a high proportion of F₁ generation hybrids, while a novel hybrid taxon between R. spinuliferum and R. scabrifolium dominated the F₂ generation, but no backcross individuals were detected. The hybrid zone showed basically coincident patterns of population structure between genomic and morphological data.

CONCLUSIONS

Natural hybridization exists among the three Rhododendron species in the hybrid zone, although patterns of hybrid formation vary between hybrid taxa, which may result in different evolutionary outcomes. This study represents a unique opportunity to dissect the ecological and evolutionary mechanisms associated with adaptive radiation of Rhododendron species in a biodiversity hotspot.

Spatiotemporal evolution of the global species diversity of Rhododendron

Evolutionary radiation is a widely recognized mode of species diversification, but its underlying mechanisms have not been unambiguously resolved for species-rich cosmopolitan plant genera. In particular, it remains largely unknown how biological and environmental factors have jointly driven its occurrence in specific regions. Here, we use Rhododendron, the largest genus of woody plants in the Northern Hemisphere, to investigate how geographic and climatic factors, as well as functional traits, worked together to trigger plant evolutionary radiations and shape the global patterns of species richness based on a solid species phylogeny. Using 3,437 orthologous nuclear genes, we reconstructed the first highly supported and dated backbone phylogeny of Rhododendron comprising 200 species that represent all subgenera, sections, and nearly all multispecies subsections, and found that most extant species originated by evolutionary radiations when the genus migrated southward from circumboreal areas to tropical/subtropical mountains, showing rapid increases of both net diversification rate and evolutionary rate of environmental factors in the Miocene. We also found that the geographically uneven diversification of Rhododendron led to a much higher diversity in Asia than in other continents, which was mainly driven by two environmental variables, that is, elevation range and annual precipitation, and were further strengthened by the adaptation of leaf functional traits. Our study provides a good example of integrating phylogenomic and ecological analyses in deciphering the mechanisms of plant evolutionary radiations, and sheds new light on how the intensification of the Asian monsoon has driven evolutionary radiations in large plant genera of the Himalaya-Hengduan Mountains.

Developing an efficient DNA barcoding system to differentiate between Lilium species

BACKGROUND

Lilium is an important ornamental bulb, possesses medicinal properties, and is also edible. Species within the Lilium genus share very similar morphology and macroscopic characteristics, thus they cannot be easily and clearly distinguished from one another. To date, no efficient species-specific markers have been developed for classifying wild lily species, which poses an issue with further characterizing its medicinal properties.

RESULTS

To develop a simple and reliable identification system for Lilium, 45 representative species from 6 sections were used to develop a DNA barcoding system, which was based on DNA sequence polymorphisms. In this study, we assessed five commonly used DNA barcode candidates (ITS, rbcL, ycf1b, matK and psbA-trnH) and five novel barcode candidates obtained from highly variable chloroplast genomic regions (trnL-trnF, trnS-trnG, trnF-ndhJ, trnP-psaJ-rpI33 and psbB-psbH). We showed that a set of three novel DNA barcodes (ITS + trnP-psaJ-rpI33 + psbB-psbH) could be efficiently used as a genetic marker to distinguish between lily species, as assessed by methods including DNAsp, BI and ML tree, and Pair Wise Group (PWG).

CONCLUSIONS

A rapid and reliable DNA barcoding method was developed for all 45 wild Lilium species by using ITS, trnP-psaJ-rpI33, and psbB-psbH as DNA barcoding markers. The method can be used in the classification of wild Lilium species, especially endangered species, and also provides an effective method for selective lily breeding.

Incongruent phylogenies and their implications for the study of diversification, taxonomy, and genome size evolution of Rhododendron

PREMISE

Classification of taxa depends on the quality of inferred phylogenies. Rhododendron, a highly species-rich genus (>1156 species) of woody plants, has a highly debated infrageneric classification, due to its huge diversity, homoplasy in key characters, and incongruence among data sets. We provide a broad coverage of representative species to resolve Rhododendron infrageneric phylogeny and highlight the areas of incongruence. We further investigate the effect of polyploidy and genome size evolution on diversification of Rhododendron.

METHODS

We generated two plastid and two nuclear loci for 260 Rhododendron species. We analyzed the loci separately as well as concatenated, utilizing both likelihood and Bayesian methods. We tested incongruence both among the data sets and with previous studies. We estimated genome sizes for 125 species through flow cytometry.

RESULTS

Our results suggest stronger support for larger subgenera; however, the smaller subgenera pose several problems; for example, R. tomentosum (former genus Ledum) occupies incongruent positions based on different DNA regions. The main shift to higher diversification in the genus occurs in the Himalayan/Southeast Asian clade of R. subg. Hymenanthes. We found that polyploidy occurs in almost all subgenera but most frequently within R. subg. Rhododendron sections Rhododendron and Schistanthe.

CONCLUSIONS

We endorse the recognition of five major clades at the subgeneric level, but a number of species cannot be confidently assigned to these clades due to incongruency. With regard to genome size evolution, results support previous reports that genome sizes of tropical plants are lower than those of colder and temperate regions and that genome downsizing promotes diversification.

Rhododendron kuomeianum (Ericaceae), a new species from northeastern Yunnan (China), based on morphological and genomic data

Rhododendron kuomeianum Y.H. Chang, J. Nielsen & Y.P. Ma, a new species of Rhododendron (Ericaceae) within subsect. Maddenia in sect. Rhododendron from Yiliang County, NE Yunnan, China, is described and illustrated. The new species is similar to R. valentinianum, but it can be easily distinguished by its sparse scales on the abaxial surface of the leaf blade, fewer flowers per inflorescence and white corolla with pale red margins. There are also differences in the widths of calyx lobes, leaf blade shape and indumentum characteristics of the petiole between the new species and Rhododendron linearilobum. We confirmed that R. kuomeianum is a new species closely related to R. valentinianum and R. changii with phylogenomic studies of 10 species within this subsection based on restriction site-associated DNA sequencing (RAD-seq) data. These phylogenomic analyses also clarified additional taxonomic problems in this subsection previously raised by morphological analysis. Our findings make a strong case for using next-generation sequencing to explore phylogenetic relationships and identify new species, especially in plants groups with complicated taxonomic problems.

Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron

Standard plant DNA barcodes based on 2-3 plastid regions, and nrDNA ITS show variable levels of resolution, and fail to discriminate among species in many plant groups. Genome skimming to recover complete plastid genome sequences and nrDNA arrays has been proposed as a solution to address these resolution limitations. However, few studies have empirically tested what gains are achieved in practice. Of particular interest is whether adding substantially more plastid and nrDNA characters will lead to an increase in discriminatory power, or whether the resolution limitations of standard plant barcodes are fundamentally due to plastid genomes and nrDNA not tracking species boundaries. To address this, we used genome skimming to recover near-complete plastid genomes and nuclear ribosomal DNA from Rhododendron species and compared discrimination success with standard plant barcodes. We sampled 218 individuals representing 145 species of this species-rich and taxonomically difficult genus, focusing on the global biodiversity hotspots of the Himalaya-Hengduan Mountains. Only 33% of species were distinguished using ITS+matK+rbcL+trnH-psbA. In contrast, 55% of species were distinguished using plastid genome and nrDNA sequences. The vast majority of this increase is due to the additional plastid characters. Thus, despite previous studies showing an asymptote in discrimination success beyond 3-4 plastid regions, these results show that a demonstrable increase in discriminatory power is possible with extensive plastid genome data. However, despite these gains, many species remain unresolved, and these results also reinforce the need to access multiple unlinked nuclear loci to obtain transformative gains in species discrimination in plants.

The Rhododendron Plant Genome Database (RPGD): a comprehensive online omics database for Rhododendron

Background

The genus Rhododendron L. has been widely cultivated for hundreds of years around the world. Members of this genus are known for great ornamental and medicinal value. Owing to advances in sequencing technology, genomes and transcriptomes of members of the Rhododendron genus have been sequenced and published by various laboratories. With increasing amounts of omics data available, a centralized platform is necessary for effective storage, analysis, and integration of these large-scale datasets to ensure consistency, independence, and maintainability.

Results

Here, we report our development of the Rhododendron Plant Genome Database (RPGD; http://bioinfor.kib.ac.cn/RPGD/), which represents the first comprehensive database of Rhododendron genomics information. It includes large amounts of omics data, including genome sequence assemblies for R. delavayi, R. williamsianum, and R. simsii, gene expression profiles derived from public RNA-Seq data, functional annotations, gene families, transcription factor identification, gene homology, simple sequence repeats, and chloroplast genome. Additionally, many useful tools, including BLAST, JBrowse, Orthologous Groups, Genome Synteny Browser, Flanking Sequence Finder, Expression Heatmap, and Batch Download were integrated into the platform.

Conclusions

RPGD is designed to be a comprehensive and helpful platform for all Rhododendron researchers. Believe that RPGD will be an indispensable hub for Rhododendron studies.

Differential expressions of anthocyanin synthesis genes underlie flower color divergence in a sympatric Rhododendron sanguineum complex

BACKGROUND

The Rhododendron sanguineum complex is endemic to alpine mountains of northwest Yunnan and southeast Tibet of China. Varieties in this complex exhibit distinct flower colors even at the bud stage. However, the underlying molecular regulations for the flower color variation have not been well characterized. Here, we investigated this via measuring flower reflectance profiles and comparative transcriptome analyses on three coexisting varieties of the R. sanguineum complex, with yellow flush pink, bright crimson, and deep blackish crimson flowers respectively. We compared the expression levels of differentially-expressed-genes (DEGs) of the anthocyanin / flavonoid biosynthesis pathway using RNA-seq and qRT-PCR data. We performed clustering analysis based on transcriptome-derived Single Nucleotide Polymorphisms (SNPs) data, and finally analyzed the promoter architecture of DEGs.

RESULTS

Reflectance spectra of the three color morphs varied distinctively in the range between 400 and 700 nm, with distinct differences in saturation, brightness, hue, and saturation/hue ratio, an indirect measurement of anthocyanin content. We identified 15,164 orthogroups that were shared among the three varieties. The SNP clustering analysis indicated that the varieties were not monophyletic. A total of 40 paralogous genes encoding 12 enzymes contributed to the flower color polymorphism. These anthocyanin biosynthesis-related genes were associated with synthesis, modification and transportation properties (RsCHS, RsCHI, RsF3H, RsF3'H, RsFLS, RsANS, RsAT, RsOMT, RsGST), as well as genes involved in catabolism and degradation (RsBGLU, RsPER, RsCAD). Variations in sequence and cis-acting elements of these genes might correlate with the anthocyanin accumulation, thus may contribute to the divergence of flower color in the R. sanguineum complex.

CONCLUSIONS

Our results suggested that the varieties are very closely related and flower color variations in the R. sanguineum complex correlate tightly with the differential expression levels of genes involved in the anabolic and catabolic synthesis network of anthocyanin. Our study provides a scenario involving intricate relationships between genetic mechanisms for floral coloration accompanied by gene flow among the varieties that may represent an early case of pollinator-mediated incipient sympatric speciation.

The complete mitochondrial and plastid genomes of Rhododendron simsii, an important parent of widely cultivated azaleas

The genus Rhododendron of the heather family (Ericaceae) is well known and widely cultivated for their highly ornamental value. The most widely cultivated Rhododendron species is Rhododendron simsii (Indoor azalea). In this study, we assembled the complete linear mitochondrial genome (GenBank accession number MW030508) and quadripartite plastid genome (GenBank accession number MW030509). The mitochondrial genome is 802,707 bp in length with containing 53 unique genes (33 protein-coding, 17 tRNA, and 3 rRNA genes), while the 152,214 bp long plastid genome is smaller and containing 105 unique genes (4 rRNA, 26 tRNA, and 75 protein-coding genes). Phylogenetic analysis showed that the same species relationship with APG system as well as the low supports of branches which is the common characteristic of resolved Ericales phylogenetics.

Identification of medicinal plants within the Apocynaceae family using ITS2 and psbA-trnH barcodes

To ensure the safety of medications, it is vital to accurately authenticate species of the Apocynaceae family, which is rich in poisonous medicinal plants. We identified Apocynaceae species by using nuclear internal transcribed spacer 2 (ITS2) and psbA-trnH based on experimental data. The identification ability of ITS2 and psbA-trnH was assessed using specific genetic divergence, BLAST1, and neighbor-joining trees. For DNA barcoding, ITS2 and psbA-trnH regions of 122 plant samples of 31 species from 19 genera in the Apocynaceae family were amplified. The PCR amplification for ITS2 and psbA-trnH sequences was 100%. The sequencing success rates for ITS2 and psbA-trnH sequences were 81% and 61%, respectively. Additional data involved 53 sequences of the ITS2 region and 38 sequences of the psbA-trnH region were downloaded from GenBank. Moreover, the analysis showed that the inter-specific divergence of Apocynaceae species was greater than its intra-specific variations. The results indicated that, using the BLAST1 method, ITS2 showed a high identification efficiency of 97% and 100% of the samples at the species and genus levels, respectively, via BLAST1, and psbA-trnH successfully identified 95% and 100% of the samples at the species and genus levels, respectively. The barcode combination of ITS2/psbA-trnH successfully identified 98% and 100% of samples at the species and genus levels, respectively. Subsequently, the neighbor joining tree method also showed that barcode ITS2 and psbA-trnH could distinguish among the species within the Apocynaceae family. ITS2 is a core barcode and psbA-trnH is a supplementary barcode for identifying species in the Apocynaceae family. These results will help to improve DNA barcoding reference databases for herbal drugs and other herbal raw materials.

Species Identification in the Rhododendron vernicosum-R. decorum Species Complex (Ericaceae)

Delimitating species boundaries is the primary aim of biological classification and could be critical for evaluating the evolving process of species and conserving biodiversity. Rhododendron is an iconic group with an extraordinary diversity in southwest China. However, it remains unknown whether the recorded species therein comprise independently evolving lineages or artificially delimitated morphological entities. In this study, we carried out species delimitation of four Rhododendron species in the R. vernicosum-R. decorum species complex based on morphological analyses and population genetic data from nuclear simple sequence repeats (SSR) markers. We randomly selected a total of 105 specimens of different individuals identified as four species across their distributional ranges to examine the statistically distinct phenotypic clusters based on 19 morphological traits. Similarly, we genotyped 55 individuals of four species from 21 populations using 15 SSR markers. The morphological analyses sorted R. decorum and the other three species into two different phenotypic clusters. The genetic clusters were consistent with the morphological clusters. However, we also recovered the third genetic cluster, comprising six R. vernicosum populations and containing the admixed genetic compositions of the other two distinct genetic clusters. This hybrid group was morphologically similar to the typical R. vernicosum (including the samples from its type specimen locality and both R. verruciferum and R. gonggashanense) but with more genetic ancestry from R. decorum. Based on our findings, we identify two distinct species and one putative hybrid group due to introgression in the R. vernicosum-R. decorum species complex. We propose to merge R. verruciferum and R. gonggashanense into R. vernicosum based on genetic compositions and our morphological analyses. The hybrid group inferred from our findings, however, needs further investigations.

Species Quantification in Complex Herbal Formulas-Vector Control Quantitative Analysis as a New Method

Product mislabeling and/or species fraud in Traditional Chinese Medicine (TCM) not only decrease TCM quality, but also pose a potential health issue to the end user. Up to now, methods to control TCM quality have been developed to detect specific metabolites or identify the original species. However, species quantification in complex herbal formulas is rarely concerned. Here, we reported a simple Vector Control Quantitative Analysis (VCQA) method for flexible and accurate multiplex species quantification in traditional Chinese herbal formulas. We developed PCR-based strategy to quickly generate the integrated DNA fragments from multiple targeted species, which can be assembled into the quantitative vector in one round of cloning by Golden Gate ligation and Gateway recombination technique. With this method, we recruited the nuclear ribosomal DNA Internal Transcribed Spacer (ITS) region for the quantification of Ligusticum sinense “Chuanxiong,” Angelica dahurica (Hoffm.) Benth. & Hook.f. ex Franch. & Sav., Notopterygium incisum K. C. Ting ex H. T. Chang, Asarum sieboldii Miq., Saposhnikovia divaricata (Turcz.) Schischk., Nepeta cataria L., Mentha canadensis L., and Glycyrrhiza uralensis Fisch. ex DC. in ChuanXiong ChaTiao Wan, a classic Chinese herbal formula with very long historical background. We found that, firstly, VCQA method could eliminate the factors affecting such as the variations in DNA extracts when in combination with the use of universal and species-specific primers. Secondly, this method detected the limit of quantification of A. sieboldii Miq. in formula products down to 1%. Thirdly, the stability of quality of ChuanXiong ChaTiao Wan formula varies significantly among different manufacturers. In conclusion, VCQA method has the potential power and can be used as an alternative method for species quantification of complex TCM formulas.

Comparative population genetic analyses suggest hybrid origin of R hododendron pubicostatum, an endangered plant species with extremely small populations endemic to Yunnan, China

Gene flow between sympatric congeneric plants is thought to be very common and may pose serious threats to endangered species. In the present study, we evaluate the genetic diversity and divergence of three sympatric Rhododendron species in Jiaozi Mountain using newly developed microsatellites through the Illumina MiSeq sequencing approach. Genetic diversity of all three Rhododendron species studied was moderate in comparison to genetic parameters previously reported from species of this genus. Interestingly, genetic structure analysis of the three species identified a possible hybrid origin of the threatened Rh. pubicostatum. This sympatry should be considered a unimodal hybrid zone, since R h. pubicostatum is predominant here. Unimodal hybrid zones are uncommon in Rhododendron, despite the fact that hybridization frequently occurs in the genus. Issues pertaining to the conservation of R h. pubicostatum resulting from admixture of genetic material from its parental species are discussed.

The utility of DNA barcodes to confirm the identification of palm collections in botanical gardens

The palm family (Arecaceae) is of high ecological and economic value, yet identification in the family remains a challenge for both taxonomists and horticulturalists. The family consists of approximately 2600 species across 181 genera and DNA barcoding may be a useful tool for species identification within the group. However, there have been few systematic evaluations of DNA barcodes for the palm family. In the present study, five DNA barcodes (rbcL, matK, trnH-psbA, ITS, ITS2) were evaluated for species identification ability across 669 samples representing 314 species and 100 genera in the Arecaceae, employing four analytical methods. The ITS gene region was found to not be a suitable barcode for the palm family, due in part, to low recovery rates and paralogous gene copies. Among the four analyses used, species resolution for ITS2 was much higher than that achieved with the plastid barcodes alone (rbcL, matK, trnH-psbA), and the barcode combination ITS2 + matK + rbcL gave the highest resolution among all single barcodes and their combinations, followed by ITS2 + matK. Among 669 palm samples analyzed, 110 samples (16.3%) were found to be misidentified. The 2992 DNA barcode sequences generated in this study greatly enriches the existing identification toolbox available to plant taxonomists that are interested in researching genetic relationships among palm taxa as well as for horticulturalists that need to confirm palm collections for botanical garden curation and horticultural applications. Our results indicate that the use of the ITS2 DNA barcode gene region provides a useful and cost-effective tool to confirm the identity of taxa in the Palm family.

Machine learning approaches outperform distance- and tree-based methods for DNA barcoding of Pterocarpus wood

Machine-learning approaches (MLAs) for DNA barcoding outperform distance- and tree-based methods on identification accuracy and cost-effectiveness to arrive at species-level identification of wood. DNA barcoding is a promising tool to combat illegal logging and associated trade, and the development of reliable and efficient analytical methods is essential for its extensive application in the trade of wood and in the forensics of natural materials more broadly. In this study, 120 DNA sequences of four barcodes (ITS2, matK, ndhF-rpl32, and rbcL) generated in our previous study and 85 downloaded from National Center for Biotechnology Information (NCBI) were collected to establish a reference data set for six commercial Pterocarpus woods. MLAs (BLOG, BP-neural network, SMO and J48) were compared with distance- (TaxonDNA) and tree-based (NJ tree) methods based on identification accuracy and cost-effectiveness across these six species, and also were applied to discriminate the CITES-listed species Pterocarpus santalinus from its anatomically similar species P. tinctorius for forensic identification. MLAs provided higher identification accuracy (30.8-100%) than distance- (15.1-97.4%) and tree-based methods (11.1-87.5%), with SMO performing the best among the machine learning classifiers. The two-locus combination ITS2 + matK when using SMO classifier exhibited the highest resolution (100%) with the fewest barcodes for discriminating the six Pterocarpus species. The CITES-listed species P. santalinus was discriminated successfully from P. tinctorius using MLAs with a single barcode, ndhF-rpl32. This study shows that MLAs provided higher identification accuracy and cost-effectiveness for forensic application over other analytical methods in DNA barcoding of Pterocarpus wood.

Prevalence of isomeric plastomes and effectiveness of plastome super-barcodes in yews (Taxus) worldwide

Taxus (yew) is both the most species-rich and taxonomically difficult genus in Taxaceae. To date, no study has elucidated the complexities of the plastid genome (plastome) or examined the possibility of whole plastomes as super-barcodes across yew species worldwide. In this study, we sequenced plastomes from two to three individuals for each of the 16 recognized yew species (including three potential cryptics) and Pseudotaxus chienii. Our comparative analyses uncovered several gene loss events that independently occurred in yews, resulting in a lower plastid gene number than other Taxaceous genera. In Pseudotaxus and Taxus, we found two isomeric arrangements that differ by the orientation of a 35 kb fragment flanked by "trnQ-IRs". These two arrangements exist in different ratios within each sampled individual, and intraspecific shifts in major isomeric arrangements are first reported here in Taxus. Moreover, we demonstrate that entire plastomes can be used to successfully discriminate all Taxus species with 100% support, suggesting that they are useful as super-barcodes for species identification. We also propose that accD and rrn16-rrn23 are promising special barcodes to discriminate yew species. Our newly developed Taxus plastomic sequences provide a resource for super-barcodes and conservation genetics of several endangered yews and serve as comprehensive data to improve models of plastome complexity in Taxaceae as a whole and authenticate Taxus species.

DNA barcoding of Corydalis, the most taxonomically complicated genus of Papaveraceae

The genus Corydalis is recognized as one of the most taxonomically challenging plant taxa. It is mainly distributed in the Himalaya-Hengduan Mountains, a global biodiversity hotspot. To date, no effective solution for species discrimination and taxonomic assignment in Corydalis has been developed. In this study, five nuclear and chloroplast DNA regions, ITS, ITS2, matK, rbcL, and psbA-trnH, were preliminarily assessed based on their ability to discriminate Corydalis to eliminate inefficient regions, and the three regions showing good performance (ITS, ITS2 and matK) were then evaluated in 131 samples representing 28 species of 11 sections of four subgenera in Corydalis using three analytical methods (NJ, ML, MP tree; K2P-distance and BLAST). The results showed that the various approaches exhibit different species identification power and that BLAST shows the best performance among the tested approaches. A comparison of different barcodes indicated that among the single barcodes, ITS (65.2%) exhibited the highest identification success rate and that the combination of ITS + matK (69.6%) provided the highest species resolution among all single barcodes and their combinations. Three Pharmacopoeia-recorded medicinal plants and their materia medica were identified successfully based on the ITS and ITS2 regions. In the phylogenetic analysis, the sections Thalictrifoliae, Sophorocapnos, Racemosae, Aulacostigma, and Corydalis formed well-supported separate lineages. We thus hypothesize that the five sections should be classified as an independent subgenus and that the genus should be divided into three subgenera. In this study, DNA barcoding provided relatively high species discrimination power, indicating that it can be used for species discrimination in this taxonomically complicated genus and as a potential tool for the authentication of materia medica belonging to Corydalis.

Species delimitation and interspecific relationships of the endangered herb genus Notopterygium inferred from multilocus variations

Species identification and discrimination is the basis of biodiversity research. In general, it is considered that numerous nucleotide variations (e.g., whole chloroplast genomes) can identify species with higher resolution than a few loci, e.g., partial chloroplast or nuclear gene fragments. In this study, we tested this hypothesis by sampling population genetics samples of the endangered herb genus Notopterygium. We sequenced the complete plastomes, five nuclear gene regions, three chloroplast DNA fragments, and a nuclear internal transcribed spacer (nrITS) region for 18 populations sampled throughout most of the geographic ranges of all six Notopterygium species. Species identification analysis showed that four DNA barcodes (matK, rbcL, trnS-trnG, and nrITS) and/or combinations of these markers achieved Notopterygium species discrimination at higher resolution than the general plastomes and nuclear gene sequences. In particular, nrITS had the highest discriminatory power among all of the individual markers. Molecular data sets and morphological evidence indicated that all six Notopterygium species could be reclassified unambiguously to four putative species clades. N. oviforme and N. franchetii had the closest relationship. Molecular dating showed that the origin and divergence of Notopterygium species was significantly associated with geological and climatic fluctuations during the middle of the Pliocene. In conclusion, our results suggest that a few nucleotide variations can achieve species discrimination with higher resolution than numerous plastomes and general nuclear gene fragments when discerning related Notopterygium species.

Efficient Identification of Pulsatilla (Ranunculaceae) Using DNA Barcodes and Micro-Morphological Characters

Pulsatilla (Ranunculaceae) comprises about 40 species, many of which have horticultural and/or medicinal importance. However, the recognition and identification of wild Pulsatilla species is difficult due to the presence of complex morphological characters. DNA barcoding is a powerful molecular tool capable of rapidly and accurately distinguishing between species. Here, we assessed the effectiveness of four commonly used DNA barcoding loci-rbcL (R), trnH-psbA ( T ), matK (M), and ITS (I)-to identify species of Pulsatilla from a comprehensive sampling group. Among the four barcoding single loci, the nuclear ITS marker showed the highest interspecific distances and the highest rate of correct identification. Among the eleven combinations, the chloroplast multi-locus R+T and R+M+T combinations were found to have the best species discrimination rate, followed by R+M. Overall, we propose that the R+M+T combination and the ITS marker on its own are, respectively, the best multi- and single-locus barcodes for discriminating among species of Pulsatilla. The phylogenetic analysis was able to distinguish species of Pulsatilla to the subgenus level, but the analysis also showed relatively low species resolution. This may be caused by incomplete lineage sorting and/or hybridization events in the evolutionary history of the genus, or by the resolution limit of the candidate barcodes. We also investigated the leaf epidermis of eight representative species using scanning electronic microscopy. The resulting micro-morphological characters were valuable for identification of related species. Using additional genome fragments, or even whole chloroplast genomes combined with micro-morphological data may permit even higher resolution of species in Pulsatilla.

Prediction of Lime Tolerance in Rhododendron Based on Herbarium Specimen and Geochemical Data

Rhododendrons are typically known to be calcifuges that cannot grow well in lime soils. Data on lime tolerance of different taxa in Rhododendron are scarce. Habitats of naturally distributed specimens of genus Rhododendron were compiled as Chinese text-based locations from the Chinese Virtual Herbarium. The locations were then geocoded into latitude/longitude pairs and subsequently connected to soil characteristics including pH and CaCO₃ from the Harmonized World Soil Database (HWSD). Using the upper quartile values of pH > 7.2 and CaCO₃ > 2% weight in topsoil as threshold, we predicted the lime tolerant taxa. A dataset of 31,146 Rhododendron specimens including the information on taxonomy, GPS locations and soil parameters for both top- and subsoil was built. The majority of the specimens were distributed in soils with moderately acidic pH and without presence of CaCO₃. 76 taxa with potential lime tolerance were predicted out of 525 taxa. The large scale data analysis based on combined data of geocoded herbarium specimens and HWSD allows identification of valuable Rhododendron species, subspecies or botanical varieties with potential tolerance to lime soils with higher pH. The predicted tolerant taxa are valuable resources for an in-depth evaluation of lime tolerance or for further use in horticulture and breeding.

DNA barcoding herbaceous and woody plant species at a subalpine forest dynamics plot in Southwest China

Although DNA barcoding has been widely used to identify plant species composition in temperate and tropical ecosystems, relatively few studies have used DNA barcodes to document both herbaceous and woody components of forest plot. A total of 201 species (72 woody species and 129 herbaceous species) representing 135 genera distributed across 64 families of seed plants were collected in a 25 ha CForBio subalpine forest dynamics plot. In total, 491 specimens were screened for three DNA regions of the chloroplast genome (rbcL, matK, and trnH-psbA) as well as the internal transcribed spacers (ITS) of nuclear ribosomal DNA. We quantified species resolution for each barcode separately or in combination using a ML tree-based method. Amplification and sequencing success were highest for rbcL, followed by trnH-psbA, which performed better than ITS and matK. The rbcL + ITS barcode had slightly higher species resolution rates (88.60%) compared with rbcL + matK (86.60%) and rbcL + trnH-psbA (86.01%). The addition of trnH-psbA or ITS to the rbcL + matK barcode only marginally increased species resolution rates, although in combination the four barcodes had the highest discriminatory power (90.21%). The situations where DNA barcodes did not discriminate among species were typically associated with higher numbers of co-occurring con-generic species. In addition, herbaceous species were much better resolved than woody species. Our study represents one of the first applications of DNA barcodes in a subalpine forest dynamics plot and contributes to our understanding of patterns of genetic divergence among woody and herbaceous plant species.

DNA barcoding a taxonomically complex hemiparasitic genus reveals deep divergence between ploidy levels but lack of species-level resolution

DNA barcoding is emerging as a useful tool not only for species identification but also for studying evolutionary and ecological processes. Although plant DNA barcodes do not always provide species-level resolution, the generation of large DNA barcode data sets can provide insights into the mechanisms underlying the generation of species diversity. Here, we study evolutionary processes in taxonomically complex British Euphrasia (Orobanchaceae), a group with multiple ploidy levels, frequent self-fertilization, young species divergence and widespread hybridization. We use a phylogenetic approach to investigate the colonization history of British Euphrasia, followed by a DNA barcoding survey and population genetic analyses to reveal the causes of shared sequence variation. Phylogenetic analysis shows Euphrasia have colonized Britain from mainland Europe on multiple occasions. DNA barcoding reveals that no British Euphrasia species has a consistent diagnostic sequence profile, and instead, plastid haplotypes are either widespread across species, or are population specific. The partitioning of nuclear genetic variation suggests differences in ploidy act as a barrier to gene exchange, while the divergence between diploid and tetraploid ITS sequences supports the polyploids being allotetraploid in origin. Overall, these results show that even when lacking species-level resolution, analyses of DNA barcoding data can reveal evolutionary patterns in taxonomically complex genera.

Integrating a comprehensive DNA barcode reference library with a global map of yews (Taxus L.) for forensic identification

Rapid and accurate identification of endangered species is a critical component of biosurveillance and conservation management, and potentially policing illegal trades. However, this is often not possible using traditional taxonomy, especially where only small or preprocessed parts of plants are available. Reliable identification can be achieved via a comprehensive DNA barcode reference library, accompanied by precise distribution data. However, these require extensive sampling at spatial and taxonomic scales, which has rarely been achieved for cosmopolitan taxa. Here, we construct a comprehensive DNA barcode reference library and generate distribution maps using species distribution modelling (SDM), for all 15 Taxus species worldwide. We find that trnL-trnF is the ideal barcode for Taxus: It can distinguish all Taxus species and in combination with ITS identify hybrids. Among five analysis methods tested, NJ was the most effective. Among 4,151 individuals screened for trnL-trnF, 73 haplotypes were detected, all species-specific and some population private. Taxonomical, geographical and genetic dimensions of sampling strategy were all found to affect the comprehensiveness of the resulting DNA barcode library. Maps from SDM showed that most species had allopatric distributions, except T. mairei in the Sino-Himalayan region. Using the barcode library and distribution map data, two unknown forensic samples were identified to species (and in one case, population) level and another was determined as a putative interspecific hybrid. This integrated species identification system for Taxus can be used for biosurveillance, conservation management and to monitor and prosecute illegal trade. Similar identification systems are recommended for other IUCN- and CITES-listed taxa.

The Application and Limitation of Universal Chloroplast Markers in Discriminating East Asian Evergreen Oaks

The East Asian subtropics mostly occupied by evergreen broad-leaved forests (EBLFs), is one of the global diversity centers for evergreen oaks. Evergreen oaks are keystone canopy trees in EBLFs with important ecosystem function and crucial significance for regional biodiversity conservation. However, the species composition and diversity of Asian evergreen oaks are poorly understood. Here, we test whether the four chloroplast markers atpI-atpH, matK, psbA-trnH, and ycf1, can discriminate the two evergreen oak sections in Asia - Cyclobalanopsis and Ilex. Two hundred and seventy-two individuals representing 57 species were scanned and 17 species from other oaks sections were included for phylogenetic reconstruction. The genetic diversity of the Quercus sections was also compared. Overall, we found that universal chloroplast DNA (cpDNA) barcoding markers could resolve two clades in Quercus, i.e., subgenus Cerris (Old World Clade) and subgenus Quercus (New World Clade). The chloroplast markers distinguished the main sections, with few exceptions. Each cpDNA region showed no barcoding gap and none of them provided good resolution at the species level. The best species resolution (27.78%) was obtained when three or four markers were combined and analyzed using BLAST. The high conservation of the cpDNA and complicated evolutionary patterns, due to incomplete lineage sorting, interspecific hybridization and introgressions may hinder the ability of cpDNA markers to discriminate different species. When comparing diversification pattern across Quercus sections (Cyclobalanopsis, Ilex, Cerris, Quercus, and Protobalanus), we found that section Ilex was the most genetically diverse, and section Cyclobalanopsis was lower genetically diverse. This diversification pattern may have resulted from the interplay of the Eurasia Cenozoic tectonic movements, climate changes and different niches of their ancestral lineages.