Development of Chloroplast and Nuclear DNA Markers for Chinese Oaks (Quercus Subgenus Quercus) and Assessment of Their Utility as DNA Barcodes

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

Comparison of chloroplast genomes and phylogenetic analysis of four species in Quercus section Cyclobalanopsis

The identification in Quercus L. species was considered to be difficult all the time. The fundamental phylogenies of Quercus have already been discussed by morphological and molecular means. However, the morphological characteristics of some Quercus groups may not be consistent with the molecular results (such as the group Helferiana), which may lead to blurring of species relationships and prevent further evolutionary researches. To understand the interspecific relationships and phylogenetic positions, we sequenced and assembled the CPGs (160,715 bp-160842 bp) of four Quercus section Cyclobalanopsis species by Illumina pair-end sequencing. The genomic structure, GC content, and IR/SC boundaries exhibited significant conservatism. Six highly variable hotspots were detected in comparison analysis, among which rpoC1, clpP and ycf1 could be used as molecular markers. Besides, two genes (petA, ycf2) were detected to be under positive selection pressure. The phylogenetic analysis showed: Trigonobalanus genus and Fagus genus located at the base of the phylogeny tree; The Quercus genus species were distincted to two clades, including five sections. All Compound Trichome Base species clustered into a single branch, which was in accordance with the results of the morphological studies. But neither of group Gilva nor group Helferiana had formed a monophyly. Six Compound Trichome Base species gathered together in pairs to form three branch respectively (Quercus kerrii and Quercus chungii; Quercus austrocochinchinensis with Quercus gilva; Quercus helferiana and Quercus rex). Due to a low support rate (0.338) in the phylogeny tree, the interspecies relationship between the two branches differentiated by this node remained unclear. We believe that Q. helferiana and Q. kerrii can exist as independent species due to their distance in the phylogeny tree. Our study provided genetic information in Quercus genus, which could be applied to further studies in taxonomy and phylogenetics.

Comparative and phylogenetic analysis of Asparagus meioclados Levl. and Asparagus munitus Wang et S. C. Chen plastomes and utility of plastomes mutational hotspots

Tiandong is a vital traditional Chinese herbal medicine. It is derived from the tuber root of the Asparagus cochinchinensis according to the Pharmacopoeia of the people's republic of China (2020 Edition). On account of the similar morphology, Asparagus meioclados and Asparagus munitus were used as Tian-Dong in southwest China. Chloroplast (cp) genomes are highly active genetic components of plants and play an extremely important role in improving the efficiency of the identification of plant species. To differentiate the medicinal plants belonging to the genus Asparagus, we sequenced and analyzed the complete plastomes (plastid genomes) of A. meioclados and A. munitus and obtained two plastomes whose length changed to 156,515 bp and 156,381 bp, respectively. A total of 111 unique genes have been detected in plastome, which included 78 protein-coding genes, 29 tRNA genes and 4 rRNA genes. In plastomes of A. meioclados and A. munitus, 14,685 and 14,987 codons were detected, among which 9942 and 10,207 had the relative synonymous codon usage (RSCU) values higher than 1, respectively. A. meioclados and A. munitus have 26 SSRs patterns, among which A. meioclados was 25 and A. munitus 21. The average Ka/Ks value was 0.36, and positive selection was detected in genes of the photosynthetic system (ndhF and rbcL) in Asparagus species. To perform the comparative analysis of plastomes, the two newly sequenced plastomes of the A. meioclados and A. munitus species were compared with that of A. cochinchinensis, and 12 hotspots, including 5 coding regions and 7 inter-genomic regions, were identified. Based on the whole plastome of Asparagus, 2 divergent hotspots (accD and rpl32-trnL-UAG) and 1 international barcode fragment (rbcL) were screened, which may be used as particular molecular markers for the identification of Asparagus species. In addition, we determined the phylogenetic relationship between A. meioclados and A. munitus in the genus Asparagus. This study enriches our knowledge of the molecular evolutionary relationships of the Asparagus genus and provides treasured data records for species identification, molecular breeding, and evolutionary analysis of this genus.

Comparative analysis of Diospyros (Ebenaceae) plastomes: Insights into genomic features, mutational hotspots, and adaptive evolution

Diospyros (Ebenaceae) is a widely distributed genus of trees and shrubs from pantropical to temperate regions, with numerous species valued for their fruits (persimmons), timber, and medicinal values. However, information regarding their plastomes and chloroplast evolution is scarce. The present study performed comparative genomic and evolutionary analyses on plastomes of 45 accepted Diospyros species, including three newly sequenced ones. Our study showed a highly conserved genomic structure across the Diospyros species, with 135-136 encoding genes, including 89 protein-coding genes, 1-2 pseudogenes (Ψycf1 for all, Ψrps19 for a few), 37 tRNA genes and 8 rRNA genes. Comparative analysis of Diospyros identified three intergenic regions (ccsA-ndhD, rps16-psbK and petA-psbJ) and five genes (rpl33, rpl22, petL, psaC and rps15) as the mutational hotspots in these species. Phylogenomic analysis identified the phylogenetic position of three newly sequenced ones and well supported a monophylogenetic (sub)temperate taxa and four clades in the pantropical taxa. The analysis codon usage identified 30 codons with relative synonymous codon usage (RSCU) values >1 and 29 codons ending with A and U bases. A total of three codons (UUA, GCU, and AGA) with highest RSCU values were identified as the optimal codons. Effective number of codons (ENC)-plot indicated the significant role of mutational pressure in shaping codon usage, while most protein-coding genes in Diospyros experienced relaxed purifying selection (d _N/d _S < 1). Additionally, the psbH gene showed positive selection (d _N/d _S > 1) in the (sub)temperate species. Thus, the results provide a meaningful foundation for further elaborating Diospyros's genetic architecture and taxonomy, enriching genetic diversity and conserving genetic resources.

Phylogenomics and plastome evolution of Indigofera (Fabaceae)

Introduction

Indigofera L. is the third largest genus in Fabaceae and includes economically important species that are used for indigo dye-producing, medicinal, ornamental, and soil and water conservation. The genus is taxonomically difficult due to the high level of overlap in morphological characters of interspecies, fewer reliability states for classification, and extensive adaptive evolution. Previous characteristic-based taxonomy and nuclear ITS-based phylogenies have contributed to our understanding of Indigofera taxonomy and evolution. However, the lack of chloroplast genomic resources limits our comprehensive understanding of the phylogenetic relationships and evolutionary processes of Indigofera.

Methods

Here, we newly assembled 18 chloroplast genomes of Indigofera. We performed a series of analyses of genome structure, nucleotide diversity, phylogenetic analysis, species pairwise Ka/Ks ratios, and positive selection analysis by combining with allied species in Papilionoideae.

Results and discussion

The chloroplast genomes of Indigofera exhibited highly conserved structures and ranged in size from 157,918 to 160,040 bp, containing 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Thirteen highly variable regions were identified, of which trnK-rbcL, ndhF-trnL, and ycf1 were considered as candidate DNA barcodes for species identification of Indigofera. Phylogenetic analysis using maximum likelihood (ML) and Bayesian inference (BI) methods based on complete chloroplast genome and protein-coding genes (PCGs) generated a well-resolved phylogeny of Indigofera and allied species. Indigofera monophyly was strongly supported, and four monophyletic lineages (i.e., the Pantropical, East Asian, Tethyan, and Palaeotropical clades) were resolved within the genus. The species pairwise Ka/Ks ratios showed values lower than 1, and 13 genes with significant posterior probabilities for codon sites were identified in the positive selection analysis using the branch-site model, eight of which were associated with photosynthesis. Positive selection of accD suggested that Indigofera species have experienced adaptive evolution to selection pressures imposed by their herbivores and pathogens. Our study provided insight into the structural variation of chloroplast genomes, phylogenetic relationships, and adaptive evolution in Indigofera. These results will facilitate future studies on species identification, interspecific and intraspecific delimitation, adaptive evolution, and the phylogenetic relationships of the genus Indigofera.

Lithocarpusdahuensis (Fagaceae), a new species from Fujian Province based on morphology and genomic data

Lithocarpusdahuensis, a new Fagaceae species from Fujian Province, China, is described and illustrated. The new species is morphologically similar to L.konishii, but its oblanceolate leaf blade has more pairs of acute teeth on the margin, denser lateral veins, smaller cupules enclosing up to 1/4-1/3 of the nut, and its nut is only half as long as those of L.konishii. The plastome of L.dahuensis was 161,303 bp in length and displayed the typical quadripartite structure. Phylogenetic analyses distinguished L.dahuensis from L.konishii with strong support based on whole plastome and nrITS, respectively.

Comparative chloroplast genome analyses of Amomum: insights into evolutionary history and species identification

BACKGROUND

Species in genus Amomum always have important medicinal and economic values. Classification of Amomum using morphological characters has long been a challenge because they exhibit high similarity. The main goals of this study were to mine genetic markers from cp genomes for Amomum species identification and discover their evolutionary history through comparative analysis.

RESULTS

Three species Amomum villosum, Amomum maximum and Amomum longipetiolatum were sequenced and annotated for the complete chloroplast (cp) genomes, and the cp genomes of A. longipetiolatum and A. maximum were the first reported. Three cp genomes exhibited typical quadripartite structures with 163,269-163,591 bp in length. Each genome encodes 130 functional genes including 79 protein-coding, 26 tRNAs and 3 rRNAs genes. 113-152 SSRs and 99 long repeats were identified in the three cp genomes. By designing specific primers, we amplified the highly variable loci and the mined genetic marker ccsA exhibited a relatively high species identification resolution in Amomum. The nonsynonymous and synonymous substitution ratios (Ka/Ks) in Amomum and Alpinia showed that most genes were subjected to a purifying selection. Phylogenetic analysis revealed the evolutionary relationships of Amomum and Alpinia species and proved that Amomum is paraphyletic. In addition, the sequenced sample of A. villosum was found to be a hybrid, becoming the first report of natural hybridization of this genus. Meanwhile, the high-throughput sequencing-based ITS2 analysis was proved to be an efficient tool for interspecific hybrid identification and with the help of the chloroplast genome, the hybrid parents can be also be determined.

CONCLUSION

The comparative analysis and mined genetic markers of cp genomes were conducive to species identification and evolutionary relationships of Amomum.

Comparative Genomic Analysis Uncovers the Chloroplast Genome Variation and Phylogenetic Relationships of Camellia Species

Camellia is the largest genus in the family Theaceae. Due to phenotypic diversity, frequent hybridization, and polyploidization, an understanding of the phylogenetic relationships between Camellia species remains challenging. Comparative chloroplast (cp) genomics provides an informative resource for phylogenetic analyses of Camellia. In this study, 12 chloroplast genome sequences from nine Camellia species were determined using Illumina sequencing technology via de novo assembly. The cp genome sizes ranged from 156,545 to 157,021 bp and were organized into quadripartite regions with the typical angiosperm cp genomes. Each genome harbored 87 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes in the same order and orientation. Differences in long and short sequence repeats, SNPs, and InDels were detected across the 12 cp genomes. Combining with the complete cp sequences of seven other species in the genus Camellia, a total of nine intergenic sequence divergent hotspots and 14 protein-coding genes with high sequence polymorphism were identified. These hotspots, especially the InDel (~400 bp) located in atpH-atpI region, had sufficient potential to be used as barcode markers for further phylogenetic analysis and species identification. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These cp genomes could facilitate the development of new molecular markers, accurate species identification, and investigations of the phylogenomic relationships of the genus Camellia.

Comparative analysis of chloroplast genomes of three medicinal Carpesium species: Genome structures and phylogenetic relationships

Carpesium (Asteraceae) is a genus that contains many plant species with important medicinal values. However, the lack of chloroplast genome research of this genus has greatly hindered the study of its molecular evolution and phylogenetic relationship. This study used the Illumina sequencing platform to sequence three medicinal plants of the Carpesium genus: Carpesium abrotanoides, Carpesium cernuum, and Carpesium faberi, obtaining three complete chloroplast genome sequences after assembly and annotation. It was revealed that the three chloroplast genomes were typical quadripartite structures with lengths of 151,389 bp (C. abrotanoides), 151,278 bp (C. cernuum), and 151,250 bp (C. faberi), respectively. A total of 114 different genes were annotated, including 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Abundant SSR loci were detected in all three chloroplast genomes, with most composed of A/T. The expansion and contraction of the IR region indicate that the boundary regions of IR/SC are relatively conserved for the three species. Using C. abrotanoides as a reference, most of the non-coding regions of the chloroplast genomes were significantly different among the three species. Five different mutation hot spots (trnC-GCA-petN, psaI, petA-psbJ, ndhF, ycf1) with high nucleotide variability (Pi) can serve as potential DNA barcodes of Carpesium species. Additionally, phylogenetic evolution analysis of the three species suggests that C. cernuum has a closer genetic relationship to C. faberi than C. abrotanoides. Simultaneously, Carpesium is a monophyletic group closely related to the genus Inula. Complete chloroplast genomes of Carpesium species can help study the evolutionary and phylogenetic relationships and are expected to provide genetic marker assistance to identify Carpesium species.

Chloroplast Genome of Lithocarpus dealbatus (Hook.f. & Thomson ex Miq.) Rehder Establishes Monophyletic Origin of the Species and Reveals Mutational Hotspots with Taxon Delimitation Potential

There is phylogenetic ambiguity in the genus Lithocarpus and subfamily Quercoideae (Family: Fagaceae). Lithocarpus dealbatus, an ecologically important tree, is the dominant species among the Quercoideae in India. Although several studies have been conducted on the species' regeneration and ecological and economic significance, limited information is available on its phylo-genomics. To resolve the phylogeny in Quercoideae, we sequenced and assembled the 161,476 bp chloroplast genome of L. dealbatus, which has a large single-copy section of 90,732 bp and a small single-copy region of 18,987 bp, separated by a pair of inverted repeat regions of 25,879 bp. The chloroplast genome contained 133 genes, of which 86 were protein-coding genes, 39 were transfer RNAs, and eight were ribosomal RNAs. Analysis of repeat elements and RNA editing sites revealed interspecific similarities within the Lithocarpus genus. DNA diversity analysis identified five highly diverged coding and noncoding hotspot regions in the four genera, which can be used as polymorphic markers for species/taxon delimitation across the four genera of Quercoideae viz., Lithocarpus, Quercus, Castanea, and Castanopsis. The chloroplast-based phylogenetic analysis among the Quercoideae established a monophyletic origin of Lithocarpus, and a closer evolutionary lineage with a few Quercus species. Besides providing insights into the chloroplast genome architecture of L. dealbatus, the study identified five mutational hotspots having high taxon-delimitation potential across four genera of Quercoideae.

Fifteen complete chloroplast genomes of Trapa species (Trapaceae): insight into genome structure, comparative analysis and phylogenetic relationships

BACKGROUND

Trapa L. is a floating-leaved aquatic plant with important economic and ecological values. However, the species identification and phylogenetic relationship within Trapa are still controversial, which necessitates the need for plastid genome information of Trapa. In this study, complete chloroplast genomes of 13 Trapa species/taxa were sequenced and annotated. Combined with released sequences, comparative analyses of chloroplast genomes were performed on the 15 Trapa species/taxa for the first time.

RESULTS

The Trapa chloroplast genomes exhibited typical quadripartite structures with lengths from 155,453 to 155,559 bp. The gene orders and contents within Trapa were conservative, but several changes were found in the microstructure. The intron loss of rpl2, also detected in Lythraceae, was found in all Trapa species/taxa, suggesting close genetic relationship between Lythraceae and Trapaceae. Notably, two small-seed species (T. incisa and T. maximowiczii) showed the smallest genome size with 155,453 and 155,477 bp, respectively. Each cp genome contained the same 130 genes consisting of 85 protein-coding genes, 37 tRNA genes and 8 rRNA genes. Trapa species/taxa showed 37 (T. incisa and T. maximowiczii) to 41 (T. sibirica) long repeats, including forward, palindromic, reversed and complementary repeats. There were 110 (T. quadrispinosa) to 123 (T. incisa and T. maximowiczii) SSR (simple sequence repeat) loci in Trapa chloroplast genomes. Comparative analyses revealed that two hotspot regions (atpA-atpF and rps2-rpoC2) in Trapa chloroplast genomes could be served as potential molecular markers. Three phylogenetic analyses (ML, MP and BI) consistently showed that there were two clusters within Trapa, including large- and small-seed species/taxa, respectively; for the large-seed Trapa, they clustered according to their geographical origin and tubercle morphology on the surface of seeds.

CONCLUSION

In summary, we have acquired the sequences of 13 Trapa chloroplast genomes, and performed the comparative analyses within Trapa for the first time. The results have helped us better identify the Trapa species/taxa and deepen the understanding of genetic basis and phylogenetic relationship of Trapa, which will facilitate the effective management and utilization of the important genetic resources in the future.

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.-Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single-copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single-copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132-133 genes, including 85-87 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37-43 long repeats and 175-178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.

Chloroplast genome characteristics and phylogenetic analysis of the medicinal plant Blumea balsamifera (L.) DC

Blumea balsamifera (L.) DC., a medicinal plant with high economic value in the Asteraceae family, is widely distributed in China and Southeast Asia. However, studies on the population structure or phylogenetic relationships with other related species are rare owing to the lack of genome information. In this study, through high-throughput sequencing, we found that the chloroplast genome of B. balsamifera was 151,170 bp in length, with a pair of inverted repeat regions (IRa and IRb) comprising 24,982 bp, a large single-copy (LSC) region comprising 82,740 bp, and a small single-copy (SSC) region comprising 18,466 bp. A total of 130 genes were identified in the chloroplast genome of B. balsamifera, including 85 protein-coding, 37 transfer RNA, and 8 ribosomal RNA genes; furthermore, sequence analysis identified 53 simple sequence repeats. Whole chloroplast genome comparison indicated that the inverted regions (IR) were more conserved than large single-copy and SSC regions. Phylogenetic analysis showed that B. balsamifera is closely related to Pluchea indica. Conclusively, the chloroplast genome of B. balsamifera was helpful for species identification and analysis of the genetic diversity and evolution in the genus Blumea and family Asteraceae.

The complete chloroplast genome of Lithocarpus hancei (Benth.) Rehd (Fagaceae) from Zhejiang, China

Lithocarpus hancei (Benth.) Rehd is a widely distributed evergreen tree with broad-leaves that dominates the lower stories of the forest in China. Here, we sequenced and assembled the complete chloroplast genome of L. hancei. The genome is 161,304 bp with one large single copy (LSC: 90,585 bp), one small single copy (SSC: 18,959 bp), and two inverted repeat (IR) regions (IRa and IRb, each 25,880 bp). It contains 117 genes, including 80 protein-coding genes, 33 tRNA genes, and four rRNA genes. Phylogenetic analysis of 21 representative cp genomes of the Fagaceae suggests Lithocarpus is monophyletic with strong bootstrap support and also that L. hancei is closely related to L. polystachyus. The cp genome is important for constructing a robust phylogeny of Lithocarpus and Fagaceae for future study.

Comparative Survey of Morphological Variations and Plastid Genome Sequencing Reveals Phylogenetic Divergence between Four Endemic Ilex Species

Holly (Ilex L.), from the monogeneric Aquifoliaceae, is a woody dioecious genus cultivated as pharmaceutical and culinary plants, ornamentals, and industrial materials. With distinctive leaf morphology and growth habitats, but uniform reproductive organs (flowers and fruits), the evolutionary relationships of Ilex remain an enigma. To date, few contrast analyses have been conducted on morphology and molecular patterns in Ilex. Here, the different phenotypic traits of four endemic Ilex species (I. latifolia, I. suaveolens, I. viridis, and I. micrococca) on Mount Huangshan, China, were surveyed through an anatomic assay and DNA image cytometry, showing the unspecified link between the examined morphology and the estimated nuclear genome size. Concurrently, the newly-assembled plastid genomes in four Ilex have lengths ranging from 157,601 bp to 157,857 bp, containing a large single-copy (LSC, 87,020–87,255 bp), a small single-copy (SSC, 18,394–18,434 bp), and a pair of inverted repeats (IRs, 26,065–26,102 bp) regions. The plastid genome annotation suggested the presence of numerable protein-encoding genes (89–95), transfer RNA (tRNA) genes (37–40), and ribosomal RNA (rRNA) genes (8). A comprehensive comparison of plastomes within eight Ilex implicated the conserved features in coding regions, but variability in the junctions of IRs/SSC and the divergent hotspot regions potentially used as the DNA marker. The Ilex topology of phylogenies revealed the incongruence with the traditional taxonomy, whereas it informed a strong association between clades and geographic distribution. Our work herein provided novel insight into the variations in the morphology and phylogeography in Aquifoliaceae. These data contribute to the understanding of genetic diversity and conservation in the medicinal Ilex of Mount Huangshan.

Framework Phylogeny, Evolution and Complex Diversification of Chinese Oaks

Oaks (Quercus L.) are ideal models to assess patterns of plant diversity. We integrated the sequence data of five chloroplast and two nuclear loci from 50 Chinese oaks to explore the phylogenetic framework, evolution and diversification patterns of the Chinese oak's lineage. The framework phylogeny strongly supports two subgenera Quercus and Cerris comprising four infrageneric sections Quercus, Cerris, Ilex and Cyclobalanopsis for the Chinese oaks. An evolutionary analysis suggests that the two subgenera probably split during the mid-Eocene, followed by intergroup divergence within the subgenus Cerris around the late Eocene. The initial diversification of sections in the subgenus Cerris was dated between the mid-Oligocene and the Oligocene-Miocene boundary, while a rapid species radiation in section Quercus started in the late Miocene. Diversification simulations indicate a potential evolutionary shift on section Quercus, while several phenotypic shifts likely occur among all sections. We found significant negative correlations between rates of the lineage diversification and phenotypic turnover, suggesting a complex interaction between the species evolution and morphological divergence in Chinese oaks. Our infrageneric phylogeny of Chinese oaks accords with the recently proposed classification of the genus Quercus. The results point to tectonic activity and climatic change during the Tertiary as possible drivers of evolution and diversification in the Chinese oak's lineage.

Comparative analyses of chloroplast genomes from 13 Lagerstroemia (Lythraceae) species: identification of highly divergent regions and inference of phylogenetic relationships

Seven divergence hotspots as plastid markers for DNA barcoding was selected, and the phylogeny of 13 Lagerstroemia species based on the cp genome data was reconstructed within Myrtales. The Lagerstroemia species used in this study originated in China and have high economic and ecological value. The shared interspecific morphological characteristics and intraspecific morphological variation resulting from hybridization among Lagerstroemia taxa have made resolving their classification problems and phylogenetic relationships difficult. Systematic comparative genomic analysis has been shown to resolve phylogenetic relationships. We sequenced and annotated 6 Lagerstroemia cp genomes (Lagerstroemia excelsa, Lagerstroemia limii, Lagerstroemia siamica, Lagerstroemia tomentosa, Lagerstroemia venusta, and Lagerstroemia calyculata) for the first time and combined them with previously published genomes for Lagerstroemia species. Bioinformatics was used to analyse the 13 cp genomes in terms of gene structure and organization, codon usage, contraction and expansion of inverted repeat regions, repeat structure, divergence hotspots, species pairwise Ka/Ks ratios and phylogenetic relationships. The length varied between 152,049 bp in Lagerstroemia subcostata and 152,521 bp in L. venusta. We selected seven divergence hotspots in the cp genomes that had the potential to act as plastid markers to distinguish Lagerstroemia species. The phylogenetic relationships within Myrtales inferred from the cp genomes of 13 Lagerstroemia species and 27 other Myrtales species were highly supported, which illustrated several novel relationships within Myrtales. Taken together, our results provide comprehensive chloroplast genomic resources, which can be used further for species identification and molecular breeding of Lagerstroemia species.

A high level of chloroplast genome sequence variability in the Sawtooth Oak Quercus acutissima

The Sawtooth Oak, Quercus acutissima Carruth., is an economically and ecologically important tree species in the family Fagaceae with a wide distribution in China. Here, we examined its intraspecific chloroplast (cp) genome variability using available and a newly sequenced genome. The new cp genome comes from a Q. acutissima individual collected from Shenyang (Northeast China; "Q. acutissima Shenyang" in the following), and then is compared with two recently published cp genomes from Tongchuan (Northwest China) and Nanjing (East China). The cp genome of Q. acutissima Shenyang exhibits a slightly larger genome size than the other two individuals, although each encodes 86 protein-coding genes, 40 tRNA genes and eight rRNA genes. We also found the length difference for the IR/SC boundary region among the three cp genomes. Sequence comparison revealed a high intraspecific genetic divergence: the three cp genomes differ by 332 sequence patterns including 77 single nucleotide polymorphisms, and 255 indels (each gap considered) scattering across 67 regions. Phylogenetic analyses based on the cp genome recovered the split between the subgenus Cerris and the subgenus Quercus, but revealed that three Q. acutissima individuals did not cluster together, indicating that even complete cp genome data fail to reproduce species boundaries in Asian members of section Cerris. Our results show that more complete plastomes covering remote ranges needs to be sequenced to provide a solid backbone for future population-scale in-depth studies and phylogenetic analysis of section Cerris.

DNA barcodes for the identification of Stephania (Menispermaceae) species

Stephania is a medicinal plants-rich genus of Menispermaceae. However, the identification of morphologically-similar species in Stephania is difficult using the currently reported methods. The indiscriminate overexploitation of Stephania plants has resulted in clinical misuse and endangerment of many species, which necessitates the development of an efficient and reliable method for species authentication. Therefore, six candidate DNA barcode sequences (ITS, ITS2, psbA-trnH, matK, rbcL, and trnL-F) were tested for their capacity to identify Stephania species. The barcodes were analyzed either as a single region or in combination by tree-based [neighbor-joining (NJ) and Bayesian inference (BI)], distance-based (PWG-distance), and sequence similarity-based (TaxonDNA) methods. Amplification and sequencing success rates were 100% for all six candidate barcodes. A comparison of six barcode regions showed that ITS exhibited the highest number of variable and informative sites (182/179), followed by psbA-trnH (173/162). DNA barcoding gap assessment showed that interspecific distances of the six barcodes were greater than intraspecific distances. The identification results showed that species discrimination rates of combination barcodes were higher than those of single-region barcodes. Based on best match and best close match methods, the ITS+psbA-trnH combination exhibited the highest discrimination power (93.93%). Further, all Stephania species could be resolved in the phylogenetic trees based on ITS+psbA-trnH (NJ, BI). This study demonstrates that DNA barcoding is an efficient method to identify Stephania species and recommends that the ITS+psbA-trnH combination is the best DNA barcode for the identification of Stephania species.

Species Identification of Oaks (Quercus L., Fagaceae) from Gene to Genome

Species identification of oaks (Quercus) is always a challenge because many species exhibit variable phenotypes that overlap with other species. Oaks are notorious for interspecific hybridization and introgression, and complex speciation patterns involving incomplete lineage sorting. Therefore, accurately identifying Quercus species barcodes has been unsuccessful. In this study, we used chloroplast genome sequence data to identify molecular markers for oak species identification. Using next generation sequencing methods, we sequenced 14 chloroplast genomes of Quercus species in this study and added 10 additional chloroplast genome sequences from GenBank to develop a DNA barcode for oaks. Chloroplast genome sequence divergence was low. We identified four mutation hotspots as candidate Quercus DNA barcodes; two intergenic regions (matK-trnK-rps16 and trnR-atpA) were located in the large single copy region, and two coding regions (ndhF and ycf1b) were located in the small single copy region. The standard plant DNA barcode (rbcL and matK) had lower variability than that of the newly identified markers. Our data provide complete chloroplast genome sequences that improve the phylogenetic resolution and species level discrimination of Quercus. This study demonstrates that the complete chloroplast genome can substantially increase species discriminatory power and resolve phylogenetic relationships in plants.

Complete Chloroplast Genomes of Vachellia nilotica and Senegalia senegal: Comparative Genomics and Phylogenomic Placement in a New Generic System

Vachellia and Senegalia are the most important genera in the subfamily Mimosoideae (Fabaceae). Recently, species from both genera were separated from the long-characterized Acacia due to their macro-morphological characteristics. However, this morpho-taxonomic differentiation struggles to discriminate some species, for example, Vachellia nilotica and Senegalia senegal. Therefore, sequencing the chloroplast (cp) genomes of these species and determining their phylogenetic placement via conserved genes may help to validate the taxonomy. Hence, we sequenced the cp genomes of V. nilotica and S. senegal, and the results showed that the sizes of the genomes are 165.3 and 162.7 kb, respectively. The cp genomes of both species comprised large single-copy regions (93,849~91,791 bp) and pairs of inverted repeats (IR; 26,093~26,008 bp). The total numbers of genes found in the V. nilotica and S. senegal cp genomes were 135 and 132, respectively. Approximately 123:130 repeats and 290:281 simple sequence repeats were found in the S. senegal and V. nilotica cp genomes, respectively. Genomic characterization was undertaken by comparing these genomes with those of 17 species belonging to related genera in Fabaceae. A phylogenetic analysis of the whole genome dataset and 56 shared genes was undertaken by generating cladograms with the same topologies and placing both species in a new generic system. These results support the likelihood of identifying segregate genera from Acacia with phylogenomic disposition of both V. nilotica and S. senegal in the subfamily Mimosoideae. The current study is the first to obtain complete genomic information on both species and may help to elucidate the genome architecture of these species and evaluate the genetic diversity among species.

Genome Comparison Reveals Mutation Hotspots in the Chloroplast Genome and Phylogenetic Relationships of Ormosia Species

The papilionoid legume genus Ormosia comprises approximately 130 species, which are distributed mostly in the Neotropics, with some species in eastern Asia and northeastern Australia. The taxonomy and evolutionary history remain unclear due to the lack of a robust species-level phylogeny. Chloroplast genomes can provide important information for phylogenetic and population genetic studies. In this study, we determined the complete chloroplast genome sequences of five Ormosia species by Illumina sequencing. The Ormosia chloroplast genomes displayed the typical quadripartite structure of angiosperms, which consisted of a pair of inverted regions separated by a large single-copy region and a small single-copy region. The location and distribution of repeat sequences and microsatellites were determined. Comparative analyses highlighted a wide spectrum of variation, with trnK-rbcL, atpE-trnS-rps4, trnC-petN, trnS-psbZ-trnG, trnP-psaJ-rpl33, and clpP intron being the most variable regions. Phylogenetic analysis revealed that Ormosia is in the Papilionoideae clade and is sister to the Lupinus clade. Overall, this study, which provides Ormosia chloroplast genomic resources and a comparative analysis of Ormosia chloroplast genomes, will be beneficial for the evolutionary study and phylogenetic reconstruction of the genus Ormosia and molecular barcoding in population genetics and will provide insight into the chloroplast genome evolution of legumes.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Quercus bawanglingensis Huang, Li et Xing, a Vulnerable Oak Tree in China

Quercus bawanglingensis Huang, Li et Xing, an endemic evergreen oak of the genus Quercus (Fagaceae) in China, is currently listed in the Red List of Chinese Plants as a vulnerable (VU) plant. No chloroplast (cp) genome information is currently available for Q. bawanglingensis, which would be essential for the establishment of guidelines for its conservation and breeding. In the present study, the cp genome of Q. bawanglingensis was sequenced and assembled into double-stranded circular DNA with a length of 161,394 bp. Two inverted repeats (IRs) with a total of 51,730 bp were identified, and the rest of the sequence was separated into two single-copy regions, namely, a large single-copy (LSC) region (90,628 bp) and a small single-copy (SSC) region (19,036 bp). The genome of Q. bawanglingensis contains 134 genes (86 protein-coding genes, 40 tRNAs and eight rRNAs). More forward (29) than inverted long repeats (21) are distributed in the cp genome. A simple sequence repeat (SSR) analysis showed that the genome contains 82 SSR loci, involving 84.15% A/T mononucleotides. Sequence comparisons among the nine complete cp genomes, including the genomes of Q. bawanglingensis, Q. tarokoensis Hayata (NC036370), Q. aliena var. acutiserrata Maxim. ex Wenz. (KU240009), Q. baronii Skan (KT963087), Q. aquifolioides Rehd. et Wils. (KX911971), Q. variabilis Bl. (KU240009), Fagus engleriana Seem. (KX852398), Lithocarpus balansae (Drake) A. Camus (KP299291) and Castanea mollissima Bl. (HQ336406), demonstrated that the diversity of SC regions was higher than that of IR regions, which might facilitate identification of the relationships within this extremely complex family. A phylogenetic analysis showed that Fagus engleriana and Trigonobalanus doichangensis form the basis of the produced evolutionary tree. Q. bawanglingensis and Q. tarokoensis, which belong to the group Ilex, share the closest relationship. The analysis of the cp genome of Q. bawanglingensis provides crucial genetic information for further studies of this vulnerable species and the taxonomy, phylogenetics and evolution of Quercus.

Complete chloroplast genome of Camellia japonica genome structures, comparative and phylogenetic analysis

Camellia is an economically, ecologically and phylogenetically valuable genus in the family Theaceae. The frequent interspecific hybridization and polyploidization makes this genus phylogenetically and taxonomically under controversial and require detailed investigation. Chloroplast (cp) genome sequences have been used for cpDNA marker development and genetic diversity evaluation. Our research newly sequenced the chloroplast genome of Camellia japonica using Illumina HiSeq X Ten platform, and retrieved five other chloroplast genomes of Camellia previously published for comparative analyses, thereby shedding lights on a deeper understanding of the applicability of chloroplast information. The chloroplast genome sizes ranged in length from 156,607 to 157,166 bp, and their gene structure resembled those of other higher plants. There were four categories of SSRs detected in six Camellia cpDNA sequences, with the lengths ranging from 10 to 17bp. The Camellia species exhibited different evolutionary routes that lhbA and orf188, followed by orf42 and psbZ, were readily lost during evolution. Obvious codon preferences were also shown in almost all protein-coding cpDNA and amino acid sequences. Selection pressure analysis revealed the influence of different environmental pressures on different Camellia chloroplast genomes during long-term evolution. All Camellia species, except C. crapnelliana, presented the identical rate of amplification in the IR region. The datasets obtained from the chloroplast genomes are highly supportive in inferring the phylogenetic relationships of the Camellia taxa, indicating that chloroplast genome can be used for classifying interspecific relationships in this genus.

The complete chloroplast genomes of three Betulaceae species: implications for molecular phylogeny and historical biogeography

Background

Previous phylogenetic conclusions on the family Betulaceae were based on either morphological characters or traditional single loci, which may indicate some limitations. The chloroplast genome contains rich polymorphism information, which is very suitable for phylogenetic studies. Thus, we sequenced the chloroplast genome sequences of three Betulaceae species and performed multiple analyses to investigate the genome variation, resolve the phylogenetic relationships, and clarify the divergence history.

Methods

Chloroplast genomes were sequenced using the high-throughput sequencing. A comparative genomic analysis was conducted to examine the global genome variation and screen the hotspots. Three chloroplast partitions were used to reconstruct the phylogenetic relationships using Maximum Likelihood and Bayesian Inference approaches. Then, molecular dating and biogeographic inferences were conducted based on the whole chloroplast genome data.

Results

Betulaceae chloroplast genomes consisted of a small single-copy region and a large single copy region, and two copies of inverted repeat regions. Nine hotspots can be used as potential DNA barcodes for species delimitation. Phylogenies strongly supported the division of Betulaceae into two subfamilies: Coryloideae and Betuloideae. The phylogenetic position of Ostryopsis davidiana was controversial among different datasets. The divergence time between subfamily Coryloideae and Betuloideae was about 70.49 Mya, and all six extant genera were inferred to have diverged fully by the middle Oligocene. Betulaceae ancestors were probably originated from the ancient Laurasia.

Discussions

This research elucidates the potential of chloroplast genome sequences in the application of developing molecular markers, studying evolutionary relationships and historical dynamic of Betulaceae.It also reveals the advantages of using chloroplast genome data to illuminate those phylogenies that have not been well solved yet by traditional approaches in other plants.

Comparative systematics and phylogeography of Quercus Section Cerris in western Eurasia: inferences from plastid and nuclear DNA variation

Oaks (Quercus) comprise more than 400 species worldwide and centres of diversity for most sections lie in the Americas and East/Southeast Asia. The only exception is the Eurasian sect. Cerris that comprises about 15 species, most of which are confined to western Eurasia. This section has not been comprehensively studied using molecular tools. Here, we assess species diversity and provide a first comprehensive taxonomic and phylogeographic scheme of western Eurasian members of sect. Cerris using plastid (trnH-psbA) and nuclear (5S-IGS) DNA variation with a dense intra-specific and geographic sampling. Chloroplast haplotypes primarily reflected phylogeographic patterns originating from interspecific cytoplasmic gene flow within sect. Cerris and its sister section Ilex. We identified two widespread and ancestral haplotypes, and locally restricted derived variants. Signatures shared with Mediterranean species of sect. Ilex, but not with the East Asian Cerris oaks, suggest that the western Eurasian lineage came into contact with Ilex only after the first (early Oligocene) members of sect. Cerris in Northeast Asia had begun to radiate and move westwards. Nuclear 5S-IGS diversification patterns were more useful for establishing a molecular-taxonomic framework and to reveal hybridization and reticulation. Four main evolutionary lineages were identified. The first lineage is comprised of Q. libani, Q. trojana and Q. afares and appears to be closest to the root of sect. Cerris. These taxa are morphologically most similar to the East Asian species of Cerris, and to both Oligocene and Miocene fossils of East Asia and Miocene fossils of western Eurasia. The second lineage is mainly composed of the widespread Q. cerris and the narrow endemic species Q. castaneifolia, Q. look, and Q. euboica. The third lineage comprises three Near East species (Q. brantii, Q. ithaburensis and Q. macrolepis), well adapted to continental climates with cold winters. The forth lineage appears to be the most derived and comprises Q. suber and Q. crenata. Q. cerris and Q.  trojana displayed high levels of variation; Q. macrolepis and Q. euboica, previously treated as subspecies of Q. ithaburensis and Q. trojana, likely deserve independent species status. A trend towards inter-specific crosses was detected in several taxa; however, we found no clear evidence of a hybrid origin of Q. afares and Q. crenata, as currently assumed.

Interspecific chloroplast genome sequence diversity and genomic resources in Diospyros

BACKGROUND

Fruits of persimmon plants are traditional healthy food in China, Korea and Japan. However, due to the shortage of morphological and DNA markers, the development of persimmon industry has been heavily inhibited.

RESULTS

Chloroplast genomes of Diospyros cathayensis, D. virginiana, D. rhombifolia and D. deyangensis were newly sequenced. Comparative analyses of ten chloroplast genomes including six previously published chloroplast genomes of Diospyros provided new insights into the genome sequence diversity and genomic resources of the genus. Eight hyper-variable regions, trnH-psbA, rps16-trnQ, rpoB-trnC, rps4-trnT-trnL, ndhF, ndhF-rpl32-trnL, ycf1a, and ycf1b, were discovered and can be used as chloroplast DNA markers at/above species levels. The complete chloroplast genome sequences provided the best resolution at inter-specific level in comparison with different chloroplast DNA sequence datasets.

CONCLUSION

Diospyros oleifera, D. deyangensis, D. virginiana, D. glaucifolia, D. lotus and D. jinzaoshi are important wild species closely related to the cultivated persimmon D. kaki. The hyper-variable regions can be used as DNA markers for global genetic diversity detection of Diospyros. Deeper study on these taxa would be helpful for elucidating the origin of D. kaki.

A nuclear DNA barcode for eastern North American oaks and application to a study of hybridization in an Arboretum setting

DNA barcoding has proved difficult in a number of woody plant genera, including the ecologically important oak genus Quercus. In this study, we utilized restrictionsite-associated DNA sequencing (RAD-seq) to develop an economical single nucleotide polymorphism (SNP) DNA barcoding system that suffices to distinguish eight common, sympatric eastern North American white oak species. Two de novo clustering pipelines, PyRAD and Stacks, were used in combination with postclustering bioinformatic tools to generate a list of 291 potential SNPs, 80 of which were included in a barcoding toolkit that is easily implemented using MassARRAY mass spectrometry technology. As a proof-of-concept, we used the genotyping toolkit to infer potential hybridization between North American white oaks transplanted outside of their native range (Q. michauxii, Q. montana, Q muehlenbergii/Q. prinoides, and Q. stellata) into a horticultural collection surrounded by natural forests of locally native trees (Q. alba and Q. macrocarpa) in the living collection at The Morton Arboretum (Lisle, IL, USA). Phylogenetic and clustering analyses suggested low rates of hybridization between cultivated and native species, with the exception of one Q. michauxii mother tree, the acorns of which exhibited high admixture from either Q. alba or Q. stellata and Q. macrocarpa, and a hybrid between Q. stellata that appears to have backcrossed almost exclusively to Q. alba. Together, RAD-seq and MassARRAY technologies allow for efficient development and implementation of a multispecies barcode for one of the more challenging forest tree genera.

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.

Genetic variation and structure of Ubame oak, Quercus phillyraeoides, in Japan revealed by chloroplast DNA and nuclear microsatellite markers

Genetic variation and structure of Ubame oak (Quercus phillyraeoides A. Gray), a shrub tree on the Pacific coast of Japan, were examined to elucidate historical population dynamics using five chloroplast DNA (cpDNA) marker regions and 11 nuclear microsatellite loci. Three cpDNA haplotypes (A, B and C) were identified in a screen of 41 populations across the entire distribution range in Japan. Haplotype A was the most prevalent and was found in the entire range. Haplotype B was locally restricted to the Kii Peninsula, while haplotype C was restricted to the southwestern part of Japan. These haplotypes corresponded with the nuclear genetic constitution revealed by microsatellite markers. Two genetically differentiated major groups were identified by STRUCTURE analysis applied to 536 individuals from 28 populations, and they mostly corresponded with the two major cpDNA haplotypes, A and C. These populations were further divided into three geographically identified groups: group 1 in the area including Kanto and Tokai regions, the Kii Peninsula and the Muroto-misaki Cape on the Pacific coast; group 2 in the Bungo-suido Channel area; and group 3 in southern Kyushu and Okinawa. Populations in the Seto Inland Sea were divided into two groups: one was included in group 1 and appears to have originated from last glacial maximum (LGM) refugia located in the Kii Peninsula, while the other was included in group 2 and appears to have originated from LGM refugia located in southern Kyushu. These groups can be considered as conservation units for the preservation of unique seashore ecosystems, or as a seed source to foster coastal protection forests and next-generation production forests. Considerable care should be taken to protect isolated populations that may be specialized to unique local environments, such as those on the islands of Koshikijima and Izenajima.

Comparative Genomics and Phylogenetic Analysis Revealed the Chloroplast Genome Variation and Interspecific Relationships of Corylus (Betulaceae) Species

Corylus L. is an economically and phylogenetically important genus in the family Betulaceae. Taxonomic and phylogenetic relationships of Corylus species have long been controversial for lack of effective molecular markers. In this study, the complete chloroplast (cp) genomes of six Corylus species were assembled and characterized using next-generation sequencing. We compared the genome features, repeat sequences, sequence divergence, and constructed the phylogenetic relationships of the six Corylus species. The results indicated that Corylus cp genomes were typical of the standard double-stranded DNA molecule, ranging from 160,445 base pairs (bp) (C. ferox var. thibetca) to 161,621 bp (C. yunnanensis) in length. Each genome contained a pair of inverted repeats (IRs), a large single-copy (LSC) region and a small single-copy (SSC) region. Each of the six cp genomes possessed 113 unique genes arranged in the same order, including 80 protein-coding, 29 tRNA, and 4 rRNA genes. C. yunnanensis contained the highest number of repeat sequences, and the richest SSRs in six cp genomes were A/T mononucleotides. Comparative analyses of six Corylus cp genomes revealed four hotspot regions (trnH-psbA, rpoB-trnC, trnF-ndhJ, and rpl32-trnL) that could be used as potential molecular markers. Phylogenetic analyses of the complete chloroplast genomes and 80 protein-coding genes exhibited nearly identical topologies that strongly supported the monophyly of Corylus and simultaneously revealed the generic relationships among Betulaceae. The availability of these genomes can offer valuable genetic information for further taxonomy, phylogeny, and species delimitation in Corylus or even Betulaceae plants.

DNA barcoding of vouchered xylarium wood specimens of nine endangered Dalbergia species

ITS2+ trnH - psbA was the best combination of DNA barcode to resolve the Dalbergia wood species studied. We demonstrate the feasibility of building a DNA barcode reference database using xylarium wood specimens. The increase in illegal logging and timber trade of CITES-listed tropical species necessitates the development of unambiguous identification methods at the species level. For these methods to be fully functional and deployable for law enforcement, they must work using wood or wood products. DNA barcoding of wood has been promoted as a promising tool for species identification; however, the main barrier to extensive application of DNA barcoding to wood is the lack of a comprehensive and reliable DNA reference library of barcodes from wood. In this study, xylarium wood specimens of nine Dalbergia species were selected from the Wood Collection of the Chinese Academy of Forestry and DNA was then extracted from them for further PCR amplification of eight potential DNA barcode sequences (ITS2, matK, trnL, trnH-psbA, trnV-trnM1, trnV-trnM2, trnC-petN, and trnS-trnG). The barcodes were tested singly and in combination for species-level discrimination ability by tree-based [neighbor-joining (NJ)] and distance-based (TaxonDNA) methods. We found that the discrimination ability of DNA barcodes in combination was higher than any single DNA marker among the Dalbergia species studied, with the best two-marker combination of ITS2+trnH-psbA analyzed with NJ trees performing the best (100% accuracy). These barcodes are relatively short regions (<350 bp) and amplification reactions were performed with high success (≥90%) using wood as the source material, a necessary factor to apply DNA barcoding to timber trade. The present results demonstrate the feasibility of using vouchered xylarium specimens to build DNA barcoding reference databases.

Development of Chloroplast Genomic Resources for Oryza Species Discrimination

Rice is the most important crop in the world as the staple food for over half of the population. The wild species of Oryza represent an enormous gene pool for genetic improvement of rice cultivars. Accurate and rapid identification of these species is critical for effective utilization of the wild rice germplasm. In this study, we developed valuable chloroplast molecular markers by comparing the chloroplast genomes for species identification. Four chloroplast genomes of Oryza were newly sequenced on the Illumina HiSeq platform and other 14 Oryza species chloroplast genomes from Genbank were simultaneously taken into consideration for comparative analyses. Among 18 Oryza chloroplast genomes, five variable regions (rps16-trnQ, trnTEYD, psbE-petL, rpoC2 and rbcL-accD) were detected for DNA barcodes, in addition to differences in simple sequence repeats (SSR) and repeat sequences. The highest species resolution (72.22%) was provided by rpoC2 and rbcL-accD with distance-based methods. Three-marker combinations (rps16-trnQ + trnTEYD + rbcL-accD, rps16-trnQ + trnTEYD + rpoC2 and rpoC2 + trnTEYD + psbE-petL) showed the best species resolution (100%). Phylogenetic analysis based on the chloroplast genome provided the best resolution of Oryza. In the comparison of chloroplast genomes in this study, identification of the most variable regions and assessment of the focal regions of divergence were efficient in developing species-specific DNA barcodes. Based on evaluation of the chloroplast genomic resources, we conclude that chloroplast genome sequences are a reliable and valuable molecular marker for exploring the wild rice genetic resource in rice improvement.

The Chloroplast Genome Sequence of Scutellaria baicalensis Provides Insight into Intraspecific and Interspecific Chloroplast Genome Diversity in Scutellaria

Scutellaria baicalensis Georgi (Lamiaceae) is the source of the well-known traditional Chinese medicine “HuangQin” (Radix Scutellariae). Natural sources of S. baicalensis are rapidly declining due to high market demand and overexploitation. Moreover, the commercial products of Radix Scutellariae have often been found to contain adulterants in recent years, which may give rise to issues regarding drug efficacy and safety. In this study, we developed valuable chloroplast molecular resources by comparing intraspecific and interspecific chloroplast genome. The S. baicalensis chloroplast genome is a circular molecule consisting of two single-copy regions separated by a pair of inverted repeats. Comparative analyses of three Scutellaria chloroplast genomes revealed six variable regions (trnH-psbA, trnK-rps16, petN-psbM, trnT-trnL, petA-psbJ, and ycf1) that could be used as DNA barcodes. There were 25 single nucleotide polymorphisms(SNPs) and 29 indels between the two S. baicalensis genotypes. All of the indels occurred within non-coding regions. Phylogenetic analysis suggested that Scutellarioideae is a sister taxon to Lamioideae. These resources could be used to explore the variation present in Scutellaria populations and for further evolutionary, phylogenetic, barcoding and genetic engineering studies, in addition to effective exploration and conservation of S. baicalensis.