Plastid genome data provide new insights into the phylogeny and evolution of the genus Epimedium

Comparative plastomes of eight subgenus Chamaesyce plants and system authentication of Euphorbiae Humifusae Herba

Euphorbiae Humifusae Herba (EHH) was provided with medicinal and edible uses, but frequently was adulterated with its closely related species. Hence, this study sought to identify EHH via an integrated approach comprising data from its morphological evaluation, HPLC analysis, comparative plastomes analysis and allele-specific PCR identification. First, the morphological characteristics of 8 subgenus Chamaesyce plants were summarized. Then, HPLC analysis showed that 18 batches of EHH were adulterated or unqualified. Furthermore, the plastomes of the 8 subg. Chamaesyce species were analyzed. Phylogenetic analysis revealed a sister relationship among the 8 subg. Chamaesyce species. The allele-specific PCR authentication was developed by the nucleotide polymorphisms (SNPs) and insertions or deletions (InDels) analysis. The results of allele-specific PCR showed that 27 batches of EHH were adulterated, indicating that the superior sensitivity of molecular authentication over the other methods used. This study provided a reference for rational use and phylogenetic research of EHH.

Complete mitochondrial genome of Melia azedarach L., reveals two conformations generated by the repeat sequence mediated recombination

Melia azedarach is a species of enormous value of pharmaceutical industries. Although the chloroplast genome of M. azedarach has been explored, the information of mitochondrial genome (Mt genome) remains surprisingly limited. In this study, we used a hybrid assembly strategy of BGI short-reads and Nanopore long-reads to assemble the Mt genome of M. azedarach. The Mt genome of M. azedarach is characterized by two circular chromosomes with 350,142 bp and 290,387 bp in length, respectively, which encodes 35 protein-coding genes (PCGs), 23 tRNA genes, and 3 rRNA genes. A pair of direct repeats (R1 and R2) were associated with genome recombination, resulting in two conformations based on the Sanger sequencing and Oxford Nanopore sequencing. Comparative analysis identified 19 homologous fragments between Mt and chloroplast genome, with the longest fragment of 12,142 bp. The phylogenetic analysis based on PCGs were consist with the latest classification of the Angiosperm Phylogeny Group. Notably, a total of 356 potential RNA editing sites were predicted based on 35 PCGs, and the editing events lead to the formation of the stop codon in the rps10 gene and the start codons in the nad4L and atp9 genes, which were verified by PCR amplification and Sanger sequencing. Taken together, the exploration of M. azedarach gap-free Mt genome provides a new insight into the evolution research and complex mitogenome architecture.

Chloroplast genomes of Caragana tibetica and Caragana turkestanica: structures and comparative analysis

BACKGROUND

The genus Caragana encompasses multiple plant species that possess medicinal and ecological value. However, some species of Caragana are quite similar in morphology, so identifying species in this genus based on their morphological characteristics is considerably complex. In our research, illumina paired-end sequencing was employed to investigate the genetic organization and structure of Caragana tibetica and Caragana turkestanica, including the previously published chloroplast genome sequence of 7 Caragana plants.

RESULTS

The lengths of C. tibetica and C. turkestanica chloroplast genomes were 128,433 bp and 129,453 bp, respectively. The absence of inverted repeat sequences in these two species categorizes them under the inverted repeat loss clade (IRLC). They encode 110 and 111 genes (4 /4 rRNA genes, 30 /31tRNA genes, and 76 /76 protein-coding genes), respectively. Comparison of the chloroplast genomes of C. tibetica and C. turkestanica with 7 other Caragana species revealed a high overall sequence similarity. However, some divergence was observed between certain intergenic regions (matK-rbcL, psbD-psbM, atpA-psbI, and etc.). Nucleotide diversity (π) analysis revealed the detection of five highly likely variable regions, namely rps2-atpI, accD-psaI-ycf4, cemA-petA, psbN-psbH and rpoA-rps11. Phylogenetic analysis revealed that C. tibetica's sister species is Caragana jubata, whereas C. turkestanica's closest relative is Caragana arborescens.

CONCLUSIONS

The present study provides worthwhile information about the chloroplast genomes of C. tibetica and C. turkestanica, which aids in the identification and classification of Caragana species.

A cornucopia of diversity-Ranunculales as a model lineage

The Ranunculales are a hyperdiverse lineage in many aspects of their phenotype, including growth habit, floral and leaf morphology, reproductive mode, and specialized metabolism. Many Ranunculales species, such as opium poppy and goldenseal, have a high medicinal value. In addition, the order includes a large number of commercially important ornamental plants, such as columbines and larkspurs. The phylogenetic position of the order with respect to monocots and core eudicots and the diversity within this lineage make the Ranunculales an excellent group for studying evolutionary processes by comparative studies. Lately, the phylogeny of Ranunculales was revised, and genetic and genomic resources were developed for many species, allowing comparative analyses at the molecular scale. Here, we review the literature on the resources for genetic manipulation and genome sequencing, the recent phylogeny reconstruction of this order, and its fossil record. Further, we explain their habitat range and delve into the diversity in their floral morphology, focusing on perianth organ identity, floral symmetry, occurrences of spurs and nectaries, sexual and pollination systems, and fruit and dehiscence types. The Ranunculales order offers a wealth of opportunities for scientific exploration across various disciplines and scales, to gain novel insights into plant biology for researchers and plant enthusiasts alike.

Comparative analyses of complete chloroplast genomes reveal interspecific difference and intraspecific variation of Tripterygium genus

The genus Tripterygium was of great medicinal value and attracted much attention on the taxonomic study using morphological and molecular methods. In this study, we assembled 12 chloroplast genomes of Tripterygium to reveal interspecific difference and intraspecific variation. The sequence length (156,692-157,061 bp) and structure of Tripterygium were conserved. Comparative analyses presented abundant variable regions for further study. Meanwhile, we determined the ndhB gene under positive selection through adaptive evolution analysis. And the phylogenetic analyses based on 15 chloroplast genomes supported the monophyly of Tripterygium hypoglaucum and the potential sister relationship between Tripterygium wilfordii and Tripterygium regelii. Molecular dating analysis indicated that the divergence time within Tripterygium was approximately 5.99 Ma (95% HPD = 3.11-8.68 Ma). The results in our study provided new insights into the taxonomy, evolution process, and phylogenetic construction of Tripterygium using complete plastid genomes.

The complete chloroplast genome sequences of four Liparis species (Orchidaceae) and phylogenetic implications

Liparis Richard (Malaxideae, Epidendroideae) is a large and diverse genus of the family Orchidaceae, the taxonomy of which is complicated and controversial. In this study, we sequenced, assembled and analyzed four complete chloroplast genomes of Liparis species including L. kumokiri, L. makinoana, L. pauliana, and L. viridiflora, and evaluated their phylogenetic relationships with related species for the first time. These four chloroplast genomes (size range 153,095 to 158,239 bp) possess typical quadripartite structures that consist of a large single copy (LSC, 83,533-86,752 bp), a small single copy (SSC, 17,938-18,156 bp) and a pair of inverted repeats (IRs, 26,421-26,933 bp). The genomes contain 133 genes, including 87 protein coding genes, 38 tRNAs and 8 rRNA genes. The genome arrangements, gene contents, gene order, long repeats and simple sequence repeats were similar with small differences observed among these four chloroplast genomes. Five highly variable regions including ycf1, ndhA, ndhF, trnQ and trnK were identified from the comparative analysis with other nine related Liparis species, which had the potential to be used as DNA markers for species identification and phylogenetic studies of Liparis species. Phylogenetic analysis based on the complete chloroplast genome sequences strongly supported the polyphyly of Liparis and its further division into three branches. These results provided valuable information to illustrate the complicated taxonomy, phylogeny and evolution process of the Liparis genus.

DNA barcoding and comparative RNA-Seq analysis provide new insights into leaf formation using a novel resource of high-yielding Epimedium koreanum

Epimedium koreanum Nakai, a well-known traditional Chinese medicinal herb, has been widely used to treat osteoporosis and sexual dysfunction for thousands of years. However, due to the decreasing population of East Asian natural resources, yearly output of Epimedium crude herb has been in low supply year by year. In this study, an unusual variety of E. koreanum was discovered in Dunhua, Jilin Province, the northernmost area where this variety was found containing 6 individuals, with three branches that had 27 leaflets, which is much more than the typical leaflet number of 9. Firstly, the novel E. koreanum varety was identified using DNA barcodes. Then, 1171 differentially expressed genes (DEGs) were discovered through parallel RNA-seq analysis between the newly discovered variety and wild type (WT) E. koreanum plant. Furthermore, the results of bioinformatics investigation revealed that 914 positively and 619 negatively correlated genes associated with the number of leaflets. Additionally, based on RNA-Seq and qRT-PCR analysis, two homologous hub TCP genes, which were commonly implicated in plant leaf development, and shown to be up regulated and down regulated in the discovered newly variety, respectively. Thus, our study discovered a novel wild resource for leaf yield rewarding medicinal Epimedium plant breeding, provided insights into the relationship between plant compound leaf formation and gene expression of TCPs transcription factors and other gene candidates, providing bases for creating high yield cultivated Epimedium variety by using further molecular selection and breeding techniques in the future.

Parallel evolution, atavism, and extensive introgression explain the radiation of Epimedium sect. Diphyllon (Berberidaceae) in southern East Asia

East Asia is the richest region of plant biodiversity in the northern temperate zone, and its radiation provides key insights for understanding rapid speciation, including evolutionary patterns and processes. However, it is challenging to investigate the recent evolutionary radiation among plants because of the lack of genetic divergence, phenotypic convergence, and interspecific gene flow. Epimedium sect. Diphyllon is a rarely studied plant lineage endemic to East Asia, especially highly diversified in its southern part. In this study, we report a robust phylogenomic analysis based on genotyping-by-sequencing data of this lineage. The results revealed a clear biogeographic pattern for Epimedium sect. Diphyllon with recognition into two major clades corresponding to the Sino-Himalayan and Sino-Japanese subkingdoms of East Asian Flora and rapid diversification of the extant species dated to the Pleistocene. Evolutionary radiation of Epimedium sect. Diphyllon is characterized by recent and predominant parallel evolution and atavism between the two subkingdom regions, with extensive reticulating hybridization within each region during the course of diversification in southern East Asia. A parallel-atavism-introgression hypothesis is referred to in explaining the radiation of plant diversity in southern East Asia, which represents a potential model for the rapid diversification of plants under global climate cooling in the late Tertiary. Our study advances our understanding of the evolutionary processes of plant radiation in East Asia as well as in other biodiversity hotspot regions.

Assembly and analysis of the mitochondrial genome of Prunella vulgaris

Prunella vulgaris (Lamiaceae) is widely distributed in Eurasia. Former studies have demonstrated that P. vulgaris has a wide range of pharmacological effects. Nevertheless, no complete P. vulgaris mitochondrial genome has been reported, which limits further understanding of the biology of P. vulgaris. Here, we assembled the first complete mitochondrial genome of P. vulgaris using a hybrid assembly strategy based on sequencing data from both Nanopore and Illumina platforms. Then, the mitochondrial genome of P. vulgaris was analyzed comprehensively in terms of gene content, codon preference, intercellular gene transfer, phylogeny, and RNA editing. The mitochondrial genome of P. vulgaris has two circular structures. It has a total length of 297, 777 bp, a GC content of 43.92%, and 29 unique protein-coding genes (PCGs). There are 76 simple sequence repeats (SSRs) in the mitochondrial genome, of which tetrameric accounts for a large percentage (43.4%). A comparative analysis between the mitochondrial and chloroplast genomes revealed that 36 homologous fragments exist in them, with a total length of 28, 895 bp. The phylogenetic analysis showed that P. vulgaris belongs to the Lamiales family Lamiaceae and P. vulgaris is closely related to Salvia miltiorrhiza. In addition, the mitochondrial genome sequences of seven species of Lamiaceae are unconservative in their alignments and undergo frequent genome reorganization. This work reports for the first time the complete mitochondrial genome of P. vulgaris, which provides useful genetic information for further Prunella studies.

Mechanism of action of icaritin on uterine corpus endometrial carcinoma based on network pharmacology and experimental evaluation

Background

Uterine corpus endometrial carcinoma (UCEC) belongs to a group of epithelial malignant tumors. Icaritin is the main active compound of Epimedii Folium. Icaritin has been utilized to induce UCEC cells to death.

Methods

We wished to identify potential targets for icaritin in the treatment of UCEC, as well as to provide a groundwork for future studies into its pharmacologic mechanism of action. Network pharmacology was employed to conduct investigations on icaritin. Target proteins were chosen from the components of icaritin for UCEC treatment. A protein-protein interaction (PPI) network was established using overlapping genes. Analyses of enrichment of function and signaling pathways were undertaken using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively, to select "hub genes". Finally, experiments were carried out to ascertain the effect of icaritin on endometrial cancer (HEC-1-A) cells.

Results

We demonstrated that icaritin has bioactive components and putative targets that are therapeutically important. Icaritin treatment induced sustained activation of the phosphoinositide 3-kinase/protein kinase B (PI3K/Akt pathway) and inhibited growth of HEC-1-A cells.

Conclusion

Our data provide a rationale for preclinical and clinical evaluations of icaritin for UCEC therapy.

Comparative plastomes of Pueraria montana var. lobata (Leguminosae: Phaseoleae) and closely related taxa: insights into phylogenomic implications and evolutionary divergence

BACKGROUND

Pueraria montana var. lobata (kudzu) is an important food and medicinal crop in Asia. However, the phylogenetic relationships between Pueraria montana var. lobata and the other two varieties (P. montana var. thomsonii and P. montana var. montana) remain debated. Although there is increasing evidence showing that P. montana var. lobata adapts to various environments and is an invasive species in America, few studies have systematically investigated the role of the phylogenetic relationships and evolutionary patterns of plastomes between P. montana var. lobata and its closely related taxa.

RESULTS

26 newly sequenced chloroplast genomes of Pueraria accessions resulted in assembled plastomes with sizes ranging from 153,360 bp to 153,551 bp. Each chloroplast genome contained 130 genes, including eight rRNA genes, 37 tRNA genes, and 85 protein-coding genes. For 24 newly sequenced accessions of these three varieties of P. montana, we detected three genes and ten noncoding regions with higher nucleotide diversity (π). After incorporated publically available chloroplast genomes of Pueraria and other legumes, 47 chloroplast genomes were used to construct phylogenetic trees, including seven P. montana var. lobata, 14 P. montana var. thomsonii and six P. montana var. montana. Phylogenetic analysis revealed that P. montana var. lobata and P. montana var. thomsonii formed a clade, while all sampled P. montana var. montana formed another cluster based on cp genomes, LSC, SSC and protein-coding genes. Twenty-six amino acid residues were identified under positive selection with the site model. We also detected six genes (accD, ndhB, ndhC, rpl2, rpoC2, and rps2) that account for among-site variation in selective constraint under the clade model between accessions of the Pueraria montana var. lobata clade and the Pueraria montana var. montana clade.

CONCLUSION

Our data provide novel comparative plastid genomic insights into conservative gene content and structure of cp genomes pertaining to P. montana var. lobata and the other two varieties, and reveal an important phylogenetic clue and plastid divergence among related taxa of P. montana come from loci that own moderate variation and underwent modest selection.

Identification of three cultivated varieties of Scutellaria baicalensis using the complete chloroplast genome as a super-barcode

Scutellaria baicalensis has been one of the most commonly used traditional Chinese medicinal plants in China for more than 2000 years. The three new varieties cultivated could not be distinguished by morphology before flowering. It will hinder the promotion of later varieties. Chloroplast DNA has been widely used in species identification. Moreover, previous studies have shown that complete chloroplast genome sequences have been suggested as super barcodes for identifying plants. Therefore, we sequenced and annotated the complete chloroplast genomes of three cultivated varieties. The chloroplast genomes of SBW, SBR, and SBP were 151,702 bp, 151,799 bp, and 151,876 bp, which contained 85 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. The analysis of the repeat sequences, codon usage, and comparison of chloroplast genomes shared a high degree of conservation. However, the sliding window results show significant differences among the three cultivated varieties in matK-rps16 and petA-psbJ. And we found that the matK-rps16 sequence can be used as a barcode for the identification of three varieties. In addition, the complete chloroplast genome contains more variations and can be used as a super-barcode to identify these three cultivated varieties. Based on the protein-coding genes, the phylogenetic tree demonstrated that SBP was more closely related to SBW, in the three cultivated varieties. Interestingly, we found that S. baicalensis and S. rehderiana are closely related, which provides new ideas for the development of S. baicalensis. The divergence time analysis showed that the three cultivated varieties diverged at about 0.10 Mya. Overall, this study showed that the complete chloroplast genome could be used as a super-barcode to identify three cultivated varieties of S. baicalensis and provide biological information, and it also contributes to bioprospecting.

Systematic Quality Evaluation of Epimedium wushanense T. S. Ying Based on Two Quality Control Standards: Total Flavonoid Glycosides and Epimedin C

Two quality control standards, total flavonoid glycosides of Epimedii Folium and epimedin C of Epimedii Wushanensis Folium, were used to systematically evaluate the quality of Epimedium wushanense T. S. Ying, so as to provide reference for its germplasm screening and resource utilization. Seven representative populations of E. wushanense covering its main distribution areas were uniformly sampled during the flowering period. There were significant quality differences among the populations of E. wushanense. According to the quality standard of total flavonoid glycosides, all populations were superior to the quality standard of the Chinese Pharmacopoeia for Epimedii Folium, with more than 1.5 % total flavonoid glycosides. The variation ranges of epimedin A, epimedin B, epimedin C, icariin and total flavonoid glycosides were 0.40-0.76 %, 0.51-0.83 %, 1.70-9.31 %, 0.40-1.23 % and 3.05-10.61 %, respectively. According to the quality standard of epimedin C, all populations were better than the quality standard of the Chinese Pharmacopoeia for Epimedii Wushanensis Folium, with more than 1.0 % epimedin C. The variation range of epimedin C was 2.22-10.06 %. When comparing the results of the two methods, a trend of slightly lower mean values was found for total flavonoid glycosides, except for the HBXW population. The quality of E. wushanense was superior to both the quality standard of Epimedii Folium and Epimedii Wushanense Folium in the Chinese Pharmacopoeia. Epimedin C was the most abundant component. Among the investigated populations, HBXW and HBGK exhibited the highest quality, and may provide excellent genetic resources for standardized cultivation. In addition, the habitat of these populations can also serve a reference for cultivation conditions.

Comparison of plastid genomes and ITS of two sister species in Gentiana and a discussion on potential threats for the endangered species from hybridization

BACKGROUND

Gentiana rigescens Franchet is an endangered medicinal herb from the family Gentianaceae with medicinal values. Gentiana cephalantha Franchet is a sister species to G. rigescens possessing similar morphology and wider distribution. To explore the phylogeny of the two species and reveal potential hybridization, we adopted next-generation sequencing technology to acquire their complete chloroplast genomes from sympatric and allopatric distributions, as along with Sanger sequencing to produce the nrDNA ITS sequences.

RESULTS

The plastid genomes were highly similar between G. rigescens and G. cephalantha. The lengths of the genomes ranged from 146,795 to 147,001 bp in G. rigescens and from 146,856 to 147,016 bp in G. cephalantha. All genomes consisted of 116 genes, including 78 protein-coding genes, 30 tRNA genes, four rRNA genes and four pseudogenes. The total length of the ITS sequence was 626 bp, including six informative sites. Heterozygotes occurred intensively in individuals from sympatric distribution. Phylogenetic analysis was performed based on chloroplast genomes, coding sequences (CDS), hypervariable sequences (HVR), and nrDNA ITS. Analysis based on all the datasets showed that G. rigescens and G. cephalantha formed a monophyly. The two species were well separated in phylogenetic trees using ITS, except for potential hybrids, but were mixed based on plastid genomes. This study supports that G. rigescens and G. cephalantha are closely related, but independent species. However, hybridization was confirmed to occur frequently between G. rigescens and G. cephalantha in sympatric distribution owing to the lack of stable reproductive barriers. Asymmetric introgression, along with hybridization and backcrossing, may probably lead to genetic swamping and even extinction of G. rigescens.

CONCLUSION

G. rigescens and G. cephalantha are recently diverged species which might not have undergone stable post-zygotic isolation. Though plastid genome shows obvious advantage in exploring phylogenetic relationships of some complicated genera, the intrinsic phylogeny was not revealed because of matrilineal inheritance here; nuclear genomes or regions are hence crucial for uncovering the truth. As an endangered species, G. rigescens faces serious threats from both natural hybridization and human activities; therefore, a balance between conservation and utilization of the species is extremely critical in formulating conservation strategies.

Insights into taxonomy and phylogenetic relationships of eleven Aristolochia species based on chloroplast genome

INTRODUCTION

The Aristolochia, as an important genus comprised of over 400 species, has attracted much interest because of its unique chemical and pharmacological properties. However, the intrageneric taxonomy and species identification within Aristolochia have long been difficult because of the complexity of their morphological variations and lack of high-resolution molecular markers.

METHODS

In this study, we sampled 11 species of Aristolochia collected from distinct habitats in China, and sequenced their complete chloroplast (cp) genomes.

RESULTS

The 11 cp genomes of Aristolochia ranged in size from 159,375bp (A. tagala) to 160,626 bp (A. tubiflora), each containing a large single-copy (LSC) region (88,914-90,251 bp), a small single-copy (SSC) region (19,311-19,917 bp), and a pair of inverted repeats (IR) (25,175-25,698 bp). These cp genomes contained 130-131 genes each, including 85 protein-coding genes (CDS), 8 ribosomal RNA genes, and 37-38 transfer RNA genes. In addition, the four types of repeats (forward, palindromic, reverse, and complement repeats) were examined in Aristolochia species. A. littoralis had the highest number of repeats (168), while A. tagala had the lowest number (42). The total number of simple sequence repeats (SSRs) is at least 99 in A. kwangsiensis, and, at most, 161 in A. gigantea. Interestingly, we detected eleven highly mutational hotspot regions, including six gene regions (clpP, matK, ndhF, psbT, rps16, trnK-UUU) and five intergenic spacer regions (ccsA-ndhD, psbZ-trnG-GCC, rpl33-rps18, rps16-trnQ-UUG, trnS-GCU-trnG-UCC). The phylogenetic analysis based on the 72 protein-coding genes showed that 11 Aristolochia species were divided into two clades which strongly supported the generic segregates of the subgenus Aristolochia and Siphisia.

DISCUSSION

This research will provide the basis for the classification, identification, and phylogeny of medicinal plants of Aristolochiaceae.

Epimediumlongnanense (Berberidaceae), a new species from Gansu, China

Epimediumlongnanense, a new species of Epimedium (Berberidaceae) from Longnan Prefecture, Gansu Province, China, is here proposed and illustrated. E.longnanense has large flowers with petals possessing long spur and obvious basal lamina, and thus should be grouped into series Davidianae. The species closely resembles E.flavum of ser. Davidianae in morphology. However, it can be easily distinguished by its elongated rhizome (vs. compact), trifoliolate leaves (vs. five leaflets, sometimes trifoliolate), pale pink or purplish-red inner sepals with 6-8 × 2-3 mm (vs. pale sulphur-yellow, ca. 11 × 4 mm).

Plastome phylogenomics provide new perspective into the phylogeny and evolution of Betulaceae (Fagales)

BACKGROUND

Betulaceae is a relatively small but morphologically diverse family, with many species having important economic and ecological values. Although plastome structure of Betulaceae has been reported sporadically, a comprehensive exploration for plastome evolution is still lacking. Besides, previous phylogenies had been constructed based on limited gene fragments, generating unrobust phylogenetic framework and hindering further studies on divergence ages, biogeography and character evolution. Here, 109 plastomes (sixteen newly assembled and 93 previously published) were subject to comparative genomic and phylogenomic analyses to reconstruct a robust phylogeny and trace the diversification history of Betulaceae.

RESULTS

All Betulaceae plastomes were highly conserved in genome size, gene order, and structure, although specific variations such as gene loss and IR boundary shifts were revealed. Ten divergent hotspots, including five coding regions (P_i > 0.02) and five noncoding regions (P_i > 0.035), were identified as candidate DNA barcodes for phylogenetic analysis and species delimitation. Phylogenomic analyses yielded high-resolution topology that supported reciprocal monophyly between Betula and Alnus within Betuloideae, and successive divergence of Corylus, Ostryopsis, and Carpinus-Ostrya within Coryloideae. Incomplete lineage sorting and hybridization may be responsible for the mutual paraphyly between Ostrya and Carpinus. Betulaceae ancestors originated from East Asia during the upper Cretaceous; dispersals and subsequent vicariance accompanied by historical environment changes contributed to its diversification and intercontinental disjunction. Ancestral state reconstruction indicated the acquisition of many taxonomic characters was actually the results of parallel or reversal evolution.

CONCLUSIONS

Our research represents the most comprehensive taxon-sampled and plastome-level phylogenetic inference for Betulaceae to date. The results clearly document global patterns of plastome structural evolution, and established a well-supported phylogeny of Betulaceae. The robust phylogenetic framework not only provides new insights into the intergeneric relationships, but also contributes to a perspective on the diversification history and evolution of the family.

Comprehensive Comparative Analysis and Development of Molecular Markers for Dianthus Species Based on Complete Chloroplast Genome Sequences

Dianthus spp. is a genus with high economic and ornamental value in the Caryophyllaceae, which include the famous fresh-cut carnation and the traditional Chinese herbal medicine, D. superbus. Despite the Dianthus species being seen everywhere in our daily lives, its genome information and phylogenetic relationships remain elusive. Thus, we performed the assembly and annotation of chloroplast genomes for 12 individuals from seven Dianthus species. On this basis, we carried out the first comprehensive and systematic analysis of the chloroplast genome sequence characteristics and the phylogenetic evolution of Dianthus. The chloroplast genome of 12 Dianthus individuals ranged from 149,192 bp to 149,800 bp, containing 124 to 126 functional genes. Sequence repetition analysis showed the number of simple sequence repeats (SSRs) ranged from 75 to 80, tandem repeats ranged from 23 to 41, and pair-dispersed repeats ranged from 28 to 43. Next, we calculated the synonymous nucleotide substitution rates (Ks) of all 76 protein coding genes to obtain the evolution rate of these coding genes in Dianthus species; rpl22 showed the highest Ks (0.0471), which suggested that it evolved the swiftest. By reconstructing the phylogenetic relationships within Dianthus and other species of Caryophyllales, 16 Dianthus individuals (12 individuals reported in this study and four individuals downloaded from NCBI) were divided into two strongly supported sister clades (Clade A and Clade B). The Clade A contained five species, namely D. caryophyllus, D. barbatus, D. gratianopolitanus, and two cultivars (‘HY’ and ‘WC’). The Clade B included four species, in which D. superbus was a sister branch with D. chinensis, D. longicalyx, and F1 ‘87M’ (the hybrid offspring F1 from D. chinensis and ‘HY’). Further, based on sequence divergence analysis and hypervariable region analysis, we selected several regions that had more divergent sequences, to develop DNA markers. Additionally, we found that one DNA marker can be used to differentiate Clade A and Clade B in Dianthus. Taken together, our results provide useful information for our understanding of Dianthus classification and chloroplast genome evolution.

The complete chloroplast genome of Epimedium muhuangense (Berberidaceae)

Epimedium muhuangense S. Z. He & Y. Y. Wang 2017, one of the rare unifoliolate species in the Epimedium genus of Berberidaceae, is distributed in the Guizhou province of China. In present research, we sequenced the complete chloroplast genome of E. muhuangense with Illumina sequencing technology. The whole genome was 157,264 bp in length, which consisted of a large single-copy region (LSC, 88,588 bp), a small single-copy region (SSC, 17,036 bp), and a pair of inverted repeat regions (IRa and IRb, 25,820 bp). A total of 112 unique genes were successfully annotated, consisting of 78 protein-encoding genes, 30 rRNA, and four tRNA. Phylogenetic analysis demonstrated that E. muhuangense is closely related to E. elachyphyllum.

The complete chloroplast genome of Epimedium jinchengshanense Y. J. Zhang & J. Q. Li (Berberidaceae), an ornamental and medicinal species

Epimedium jinchengshanense Y. J. Zhang & J. Q. Li 2014 is an important ornamental and medicinal herb, but of unclear taxonomy. In this study, the complete chloroplast genome of E. jinchengshanense was sequenced. The genome was 157,169 bp in length, with a large single-copy region of 88, 520 bp, a small single-copy region of 17,075 bp and 2 inverted repeat regions of 25, 787 bp. The genome consisted of 113 unique genes, including 79 protein-coding genes, 30 tRNA genes and 4 rRNA genes. The GC contents were 38.78%. Phylogenetic analysis showed a close relationship between E. jinchengshanense and E. ilicifolium, which was explained by the morphological similarity of flowers and leaves of the two species.

The complete chloroplast genome of Epimedium platypetalum K. Mey. (Berberidaceae), a rare plant species from China

Epimedium L. is an important medicinal herbaceous genus in the family Berberidaceae. Epimedium platypetalum K. Mey. is a plant species only narrowly distributed in the western part of China. Here, the complete chloroplast genome of Epimedium platypetalum was assembled. The chloroplast genome of E. platypetalum was 159,088 bp in length, with a total GC content of 38.79%. A total of 112 unique genes were identified, among which 78 are protein-coding genes, 30 tRNA genes, and four rRNA genes. Phylogenetic results revealed that E. platypetalum formed a sister relationship with E. membranaceum K. Mey. Our findings provided valuable data for future research on phylogenetic relationship and germplasm exploration within the genus Epimedium.