CPGAVAS2, an integrated plastome sequence annotator and analyzer

Insights into chloroplast genome structure, intraspecific variation, and phylogeny of Cyclamen species (Myrsinoideae)

Species from the flowering plant genus Cyclamen are popular amongst consumers. In particular Cyclamen persicum Mill. has been significantly used commercially, and certain small flowering species such as Cyclamen hederifolium and Cyclamen coum are gradually growing in popularity in the potted flower market. Here, the chloroplast genomes of nine Cyclamen samples including four Cyclamen species and five varieties of C. hederifolium were sequenced for genome structure comparison, White green septal striped leaves related gene screening and DNA molecular markers were developed for phylogenetic analysis. In comparing Cyclamen species' chloroplast genomes, gene content and gene order were found to be highly similar with the length of genomes ranging from 151,626 to 153,058 bp. The chloroplast genome of Cyclamen has 128 genes, including 84 protein-coding genes, 36 transfer RNA genes, and 8 ribosomal RNA genes. Based on intraspecific variation, seven hotspots, including three genes and four intergenic regions, were identified as variable markers for downstream species delimitation and interspecific relationship analyses. Moreover, a phylogenetic tree constructed with complete chloroplast genomes, revealed that Cyclamen are monophyletic with Lysimachia as the closest neighbor. Phylogenetic analyses of the 14 Cyclamen species with the seven variable regions showed five distinct clades within this genus. The highly supported topologies showed these seven regions may be used as candidate DNA barcode sequences to distinguish Cyclamen species. White green septal striped leaves is common in C. hederifolium, however the molecular mechanism of this has not yet been described. Here, we find that the intergenic region rps4-trnT-UGU seems related to white green septal striped leaves.

The complete chloroplast genome of a fast-growing tree Lophostemon confertus (Myrtaceae)

Lophostemon confertus (Myrtaceae), a fast-growing ornamental tree, is widely cultivated in tropical and subtropical regions. To determine its phylogenetic position within Myrtaceae, here we report its complete chloroplast (cp) genome, which is 160,297 bp long and contains two inverted repeats (IRs) of 26,490 bp each, separated by a small single-copy region of 18,826 bp and a large single-copy region of 88,491 bp. The cp genome contains 123 genes, including 73 unique protein-coding genes (six duplicated in the IR regions), 29 unique tRNA genes (seven duplicated in the IR regions), and four unique rRNA genes (all located in the IR regions). Phylogenetic analysis of 18 species of Myrtaceae showed that L. confertus is sister to Xanthostemon chrysanthus. The complete cp genome of L. confertus provides a valuable genetic resource for further phylogenetic studies.

The complete chloroplast genome sequence of Centaurea cyanus (Asteraceae)

Centaurea cyanus has been a weed in farmland for a long time. In this study, the chloroplast genome of C. cyanus was sequenced to establish the phylogenetic relationship between its genomic characteristics and other related species. The chloroplast gene structure of C. cyanus is a circular molecule with a length of 152,433 bp, including a large single-copy (LSC) region of 83,464 bp, a small single-copy (SSC) region of 18,545 bp, and a pair of inverted repeats sequences (IRs) region of 25,212 bp. The whole genome contains 130 genes, including 86 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenetic analysis showed that C. cyanus is close to Carthamus. tinctorius, C. tinctorius, C. diffusa, and C. maculosa, and all of them were in one clade. This study provides genetic resource information for the further study of Centaurea.

The characteristics of the complete chloroplast genome of Staurogyne concinnula (Hance) O. Kuntze (Acanthaceae)

Staurogyne concinnula (Hance) O. Kuntze (Acanthaceae) is an important ornamental herb mainly distributed in the southern region of China, including Fujian, Guangdong, Hainan, and Taiwan provinces. However, the complete chloroplast genome of S. concinnula, which could serve as a genetic resource for studies on its taxonomy and evolution, is poorly studied at present. In this study, we reported the complete chloroplast genome of S. concinnula that was assembled using high-throughput sequencing data. The chloroplast genome was 153,783 bp long, with a typical quadripartite structure containing a small single-copy region (SSC; 17,855 bp), a large single-copy region (LSC; 84,636 bp) and a pair of inverted repeats (IRs; each 25,646 bp). The overall GC content of the chloroplast genome was 38.04%. A total of 86 protein-coding genes (PCGs), 8 rRNA genes, and 37 tRNA genes were predicted. Phylogenetic analysis based on the combined sequences of 86 PCGs with the other 16 closely related species of Acanthaceae indicated that S. concinnula is closely related to Avicennia marina. The genomic data and finding from the phylogenetic studies of S. concinnula could provide useful information and give light to in-depth studies on the evolution pattern of the understudied species, as well as Staurogyne.

The complete chloroplast genome sequence of Zannichellia palustris (Potamogetonaceae)

Zannichellia palustris Linnaeus 1753 is a cosmopolitan submerged species capable of rapidly responding to environmental changes, with potential applications in the ecological treatment of heavy metal pollution in water bodies. This study aimed to characterize the complete chloroplast genome of Z. palustris, which has not been reported previously. The chloroplast genome of Z. palustris displays a quadripartite structure with a length of 155,262 base pairs (bp), comprising a large single copy (LSC) region of 85,397 bp, a small single copy (SSC) region of 18,057 bp, and a pair of inverted repeat (IR) regions of 25,904 bp. The GC content of the genome is 35.8%, with corresponding values of 33.4% for the LSC, 28.2% for the SSC, and 42.5% for the IR regions. The genome contained 130 genes, including 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Phylogenetic analysis within the order Alismatales revealed that Z. palustris clusters with the clade of Potamogeton perfoliatus, P. crispus and Stuckenia pectinata.

The complete chloroplast genome of Striga asiatica (L.) Kuntze 1891 (Orobanchaceae), a hemiparasitic weed from Guangxi China

Striga asiatica (L.) Kuntze 1891 is a hemiparasitic plant native to Asia and Africa. It is invasive and causes yield losses in crops such as corn, rice and sorghum. Lack of chloroplast genomic data has limited research into its obligate parasitic lifestyle. In this study, the complete chloroplast genome of Striga asiatica was sequenced and characterized. It is a quadripartite structure with a total length of 191,085 bp and a GC content of 37.86%. It has a large single copy region (LSC) of 51,406 bp, a small single copy region (SSC) of 273 bp, and two copies of the reverse repeat sequence (IRA and IRB) of 69,703 bp. A total of 122 protein-coding genes, 8 rRNA genes, and 44 tRNA genes were annotated in the chloroplast genome. There were a lot of ndh gene deletions and pseudogenizations in this chloroplast genome. For example, ndhA, D, E, H, I, and K were all pseudogenes because they were missing the 5' end start codon. ndhB, C, and J had shorter gene lengths than their homologs, and ndhF and ndhG were missing genes. The phylogenetic tree reveals that all Striga species form a clade, and a bootstrap value of 100 indicates that S. asiatica is closely related to Striga hermonthica and Striga sepera. The comprehensive chloroplast genomic resource of S. asiatica would assist researchers in comprehending hemiparasitic mechanisms, molecular markers, and evolutionary patterns of the genus Striga.

The first complete chloroplast genome in Engelhardia sensu stricto, Engelhardia hainanensis Chen: genome characterization and its phylogenetic relationships within the family Juglandaceae

Trees of Engelhardia are important components of subtropical and tropical forests in South-eastern Asia with great ecological and economic values. However, phylogenetic relationships within Engelhardioideae (Juglandaceae) remains obscure. In this study, we report the first complete chloroplast genome sequences of Engelhardia sensu stricto, Engelhardia hainanensis Chen, a rare species endemic in southern China. Its complete chloroplast genome is 161,574 bp in length, with a typical quadripartite structure that includes a large single-copy region of 91,158 bp, a small single-copy region of 18,790 bp, and its GC content is 35.8%. A total of 128 genes were identified, including 83 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Furthermore, a phylogenetic tree of Juglandaceae was constructed based the complete chloroplast genome sequence, which strongly support the three-subfamily classification system in Juglandaceae, and E. hainanensis was resolved sister to two Alfaropsis species. This study provides valuable genomic information for the species identification and phylogenetic study of Juglandaceae.

Assembly of the Complete Mitochondrial Genome of Gelsemium elegans Revealed the Existence of Homologous Conformations Generated by a Repeat Mediated Recombination

Gelsemium elegans (G. elegans) is a Chinese medicinal plant with substantial economic and feeding values. There is a lack of detailed studies on the mitochondrial genome of G. elegans. In this study, the mitochondrial genome of G. elegans was sequenced and assembled, and its substructure was investigated. The mitochondrial genome of G. elegans is represented by two circular chromosomes of 406,009 bp in length with 33 annotated protein-coding genes, 15 tRNA genes, and three rRNA genes. We detected 145 pairs of repeats and found that four pairs of repeats could mediate the homologous recombination into one major conformation and five minor conformations, and the presence of conformations was verified by PCR amplification and Sanger sequencing. A total of 124 SSRs were identified in the G. elegans mitochondrial genome. The homologous segments between the chloroplast and mitochondrial genomes accounted for 5.85% of the mitochondrial genome. We also predicted 477 RNA potential editing sites and found that the nad4 gene was edited 38 times, which was the most frequent occurrence. Taken together, the mitochondrial genome of G. elegans was assembled and annotated. We gained a more comprehensive understanding on the genome of this medicinal plant, which is vital for its effective utilization and genetic improvement, especially for cytoplasmic male sterility breeding and evolution analysis in G. elegans.

Plastid Phylogenomics Provide Evidence to Accept Two New Members of Ligusticopsis (Apiaceae, Angiosperms)

Peucedanum nanum and P. violaceum are recognized as members of the genus Peucedanum because of their dorsally compressed mericarps with slightly prominent dorsal ribs and narrowly winged lateral ribs. However, these species are not similar to other Peucedanum taxa but resemble Ligusticopsis in overall morphology. To check the taxonomic positions of P. nanum and P. violaceum, we sequenced their complete plastid genome (plastome) sequences and, together with eleven previously published Ligusticopsis plastomes, performed comprehensively comparative analyses. The thirteen plastomes were highly conserved and similar in structure, size, GC content, gene content and order, IR borders, and the patterns of codon bias, RNA editing, and simple sequence repeats (SSRs). Nevertheless, twelve mutation hotspots (matK, ndhC, rps15, rps8, ycf2, ccsA-ndhD, petN-psbM, psbA-trnK, rps2-rpoC2, rps4-trnT, trnH-psbA, and ycf2-trnL) were selected. Moreover, both the phylogenetic analyses based on plastomes and on nuclear ribosomal DNA internal transcribed spacer (ITS) sequences robustly supported that P. nanum and P. violaceum nested in Ligusticopsis, and this was further confirmed by the morphological evidence. Hence, transferring P. nanum and P. violaceum into Ligusticopsis genus is reasonable and convincing, and two new combinations are presented.

Comparative phylogenetic analysis of oolong tea (Phoenix Dancong tea) using complete chloroplast genome sequences

Phoenix Dancong tea, a variety of oolong tea, is produced in Chaozhou, Guangdong Province, China, and is characterized by numerous hybridizations and polyploidization. To assess the genetic diversity and phylogenetic relationships among Phoenix Dancong tea and other oolong teas, an integrated circular chloroplast genome was constructed for thirty species of Phoenix Dancong tea from Chaozhou. The genome of Phoenix dancong tea is a circular molecule of 157,041-157,137 bp, with a pair of inverted repeats (26,072-26,610 bp each) separated by a large single copy (86,615-86,658 bp) and small single copy (18,264-18,284 bp). A total of 135 unique genes were encoded, including 90 protein coding genes, 37 tRNAs and 8 rRNAs. A comparative analysis with the other seven species in the oolong tea family that have been sequenced to date revealed similarities in structural organization, gene content and arrangement. Repeated sequence analysis identified 17-23 tandem repeats, 20-24 forward repeats and 25-27 palindromic repeats. Additionally, a total of 65-70 simple sequence repeats were detected, with mononucleotide repeats being the most common. Phylogenetic analyses showed that Phoenix Dancong tea and Fujian oolong tea were clustered with other cultivated Camellia sinensis in the genus Camellia of the family Theaceae, while the two oolong tea species were relatively independently cross-embedded in the genus, Camellia. Close genetic relationships were observed between Phoenix Dancong tea and other oolong teas, and the overall chloroplast genomes of oolong tea showed patterns with low variations and conserved evolution. The availability of Phoenix Dancong tea chloroplast genomes not only elucidated the relationship among oolong teas from different origins in Guangdong and Fujian but also provided valuable genetic resources to assist further molecular studies on the taxonomic and phylogenomic resolution of the genus Camellia.

Complete Chloroplast Genome Sequence of Endangered Species in the Genus Opisthopappus C. Shih: Characterization, Species Identification, and Phylogenetic Relationships

Opisthopappus C. Shih is a rare genus of the Asteraceae family native to the Taihang Mountains in China. Due to the narrow distribution area, poor reproduction ability and human harvesting, Opisthopappus is threatened by extinction. However, the limited genetic information within Opisthopappus impede understanding of the conservation efforts and bioprospecting. Therefore, in this study, we reported the complete chloroplast (cp) genome sequences of two Opisthopappus species, including Opisthopappus taihangensis and Opisthopappus longilobus. The cp genomes of O. taihangensis and O. longilobus were 151,117 and 151,123 bp, which contained 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The repeat sequences, codon usage, RNA-editing sites, and comparative analyses revealed a high degree of conservation between the two species. The ycf1 gene was identified as a potential molecular marker. The phylogenetic tree demonstrated that O. longilobus was a separate species and not a synonym or variety of O. taihangensis. The molecular clock showed that two species diverge over a large time span, O. longilobus diverged at 15.24 Mya (Million years ago), whereas O. taihangensis diverged at 5.40 Mya We found that Opisthopappus and Ajania are closely related, which provides new ideas for the development of Opisthopappus. These results provide biological information and an essential basis to understand the evolutionary history of the Opisthopappus species, which will aid in the future the bioprospecting and conservation of endangered species.

The Complete Chloroplast Genome Sequence of Machilus chuanchienensis (Lauraceae): Genome Structure and Phylogenetic Analysis

Machilus chuanchienensis is an ecological tree distributed in southwestern China. It has a significant valuation with making Hawk tea using its leaves, an ethnic traditional tea-like beverage with a long history in Chinese tea culture. The whole chloroplast (cp) genome is an ideal model for the phylogenetic study of Lauraceae because of its simple structure and highly conserved features. There have been numerous reports of complete cp genome sequences in Lauraceae, but little is known about M. chuanchienensis. Here, the next-generation sequencing (NGS) was used to sequence the M. chuanchienensis cp genome. Then, a comprehensive comparative genome analysis was performed. The results revealed that the M. chuanchienensis's cp genome measured 152,748 base pairs (bp) with a GC content of 39.15% and coded 126 genes annotated, including comprising eight ribosomal RNA (rRNA), 36 transporter RNA (tRNA), and 82 protein-coding genes. In addition, the cp genome presented a typical quadripartite structure comprising a large single-copy (LSC; 93,811) region, a small single-copy (SSC; 18,803) region, and the inverted repeats (IRs; 20,067) region and contained 92 simple sequence repeat (SSR) locus in total. Phylogenetic relationships of 37 species indicated that M. chuanchienensis was a sister to M. balansae, M. melanophylla, and M. minutiflora. Further research on this crucial species may benefit significantly from these findings.

Complete chloroplast genomes of 11 Sabia samples: Genomic features, comparative analysis, and phylogenetic relationship

The genus Sabia is a woody climber belonging to the family Sabiaceae, order Proteales. Several species of this genus have been utilized as medicines for treating diseases, such as rheumatic arthritis, traumatism, hepatitis, etc. However, the lack of molecular data has prevented the accurate identification and refinement of taxonomic relationships in this genus. In this study, chloroplast genomes of 11 samples of the genus Sabia were assembled and analyzed. These chloroplast genomes showed a typical quadripartite structure and ranged in length from 160,956 to 162,209 bp. The structure of the genomes was found to be relatively conserved, with 130 genes annotated, including 85 coding genes, 37 tRNA genes, and eight rRNA genes. A total of 78-98 simple sequence repeats and 52-61 interspersed repeats were detected. Sequence alignment revealed 11 highly variable loci in chloroplast genomes. Among these loci, ndhF-ndhD achieved a remarkably higher resolution than the other regions. In addition, phylogenetic analysis indicated that Sect. Pachydiscus and Sect. Sabia of Sabia did not form two separate monophyletic groups. The divergence time calculated based on the Reltime method indicated that the evolutionary branches of Sabia and Meliosma started to form approximately 85.95 million years ago (Mya), and the species within Sabia began to diverge approximately 7.65 Mya. In conclusion, our study provides a basis for comprehensively exploring the phylogenetic relationships of Sabia. It also provides a methodological basis and data support for establishing a standardized and scientific identification system for this genus.

Structural Diversities and Phylogenetic Signals in Plastomes of the Early-Divergent Angiosperms: A Case Study in Saxifragales

As representative of the early-divergent groups of angiosperms, Saxifragales is extremely divergent in morphology, comprising 15 families. Within this order, our previous case studies observed significant structural diversities among the plastomes of several lineages, suggesting a possible role in elucidating their deep phylogenetic relationships. Here, we collected 208 available plastomes from 11 constituent families to explore the evolutionary patterns among Saxifragales. With thorough comparisons, the losses of two genes and three introns were found in several groups. Notably, 432 indel events have been observed from the introns of all 17 plastomic intron-containing genes, which could well play an important role in family barcoding. Moreover, numerous heterogeneities and strong intrafamilial phylogenetic implications were revealed in pttRNA (plastomic tRNA) structures, and the unique structural patterns were also determined for five families. Most importantly, based on the well-supported phylogenetic trees, evident phylogenetic signals were detected in combinations with the identified pttRNAs features and intron indels, demonstrating abundant lineage-specific characteristics for Saxifragales. Collectively, the results reported here could not only provide a deeper understanding into the evolutionary patterns of Saxifragales, but also provide a case study for exploring the plastome evolution at a high taxonomic level of angiosperms.

The complete chloroplast genome of Dendrobium moschatum (Buch.-Ham.) Sw. 1805 (Orchidaceae)

The morphological characteristic of Dendrobium moschatum (Buch.-Ham.) Sw. 1805 is very distinctive among Dendrobium Sw. 1799, and it has high medicinal and ornamental values. Here, we reported the first complete chloroplast genome of D. moschatum. The complete genome of D. moschatum was 159,701 bp in length with 130 genes, including 38 tRNA, 8 rRNA, and 84 protein-coding genes. Phylogenetic analysis showed that D. moschatum was strongly allied with D. denneanum Kerr. 1933.

Phylogenomic analyses based on the plastid genome and concatenated nrDNA sequence data reveal cytonuclear discordance in genus Atractylodes (Asteraceae: Carduoideae)

Atractylodes species are widely distributed across East Asia and are cultivated as medicinal herbs in China, Japan, and Korea. Their unclear morphological characteristics and low levels of genetic divergence obscure the taxonomic relationships among these species. In this study, 24 plant samples were collected representing five species of Atractylodes located in China; of these, 23 belonged to members of the A. lancea complex. High-throughput sequencing was used to obtain the concatenated nrDNA sequences (18S-ITS1-5.8S-ITS2-28S) and plastid genomes. The concatenated nrDNA sequence lengths for all the Atractylodes species were 5,849 bp, and the GC content was 55%. The lengths of the whole plastid genome sequences ranged from 152,138 bp (A. chinensis) to 153,268 bp (A. lancea), while their insertion/deletion sites were mainly distributed in the intergenic regions. Furthermore, 33, 34, 36, 31, and 32 tandem repeat sequences, as well as 30, 30, 29, 30, and 30 SSR loci, were detected in A. chinensis, A. koreana, A. lancea, A. japonica, and A. macrocephala, respectively. In addition to these findings, a considerable number of heteroplasmic variations were detected in the plastid genomes, implying a complicated phylogenetic history for Atractylodes. The results of the phylogenetic analysis involving concatenated nrDNA sequences showed that A. lancea and A. japonica formed two separate clades, with A. chinensis and A. koreana constituting their sister clade, while A. lancea, A. koreana, A. chinensis, and A. japonica were found based on plastid datasets to represent a mixed clade on the phylogenetic tree. Phylogenetic network analysis suggested that A. lancea may have hybridized with the common ancestor of A. chinensis and A. japonica, while ABBA-BABA tests of SNPs in the plastid genomes showed that A. chinensis was more closely related to A. japonica than to A. lancea. This study reveals the extensive discordance and complexity of the relationships across the members of the A. lancea complex (A. lancea, A. chinensis, A. koreana, and A. japonica) according to cytonuclear genomic data; this may be caused by interspecific hybridization or gene introgression.

New insights into the phylogeny of the complex thalloid liverworts (Marchantiopsida) based on chloroplast genomes

Marchantiopsida (complex thalloid liverworts) are one of the earliest lineages of embryophytes (land plants), and well-known for their air pores and chambers, pegged rhizoids, and absence of organellular RNA editing sites. Despite their importance to an understanding of early embryophyte evolution, many key nodes within this class remain poorly resolved, owing to the paucity of genetic loci previously available for phylogenetic analyses. Here, we sequenced 54 plastomes, representing 28 genera, nearly all families, and all orders of Marchantiopsida. Based on these plastomes, we present a hypothesis of deep relationships within the class, and make the first investigations of gene contents and synteny. Overall, the Marchantiopsida plastomes were well-conserved, with the exception of the genus Cyathodium that has plastomes with higher GC content, fewer single sequence repeats (SSRs), and more structural variations, implying that this genus might possess RNA editing sites. Abundant repetitive elements and six highly divergent regions were identified as suitable for future infrafamilial taxonomic studies. The phylogenetic topology of Sphaerocarpales, Neohodgsoniales and Blasiales within Marchantiopsida was essentially congruent with previous studies but generally we obtained higher support values. Based on molecular evidence and previous morphological studies, we include Lunulariales in Marchantiales and suggest the retention of narrowed delimitation of monotypic families. The phylogenetic relationships within Marchantiales were better resolved, and 13 monophyletic families were recovered. Our analyses confirmed that the loss of intron 2 of ycf3 is a synapomorphy of Marchantiidae. Finally, we propose a new genus, Asterellopsis (Aytoniaceae), and present an updated classification of Marchantiopsida. The highly supported phylogenetic backbone provided here establishes a framework for future comparative and evolutionary studies of the complex thalloid liverworts.

The Complete Chloroplast Genome Sequence of Laportea bulbifera (Sieb. et Zucc.) Wedd. and Comparative Analysis with Its Congeneric Species

Laportea bulbifera (L. bulbifera) is an important medicinal plant of Chinese ethnic minorities, with high economic and medicinal value. However, the medicinal materials of the genus Laportea are prone to be misidentified due to the similar morphological characteristics of the original plants. Thus, it is crucial to discover their molecular marker points and to precisely identify these species for their exploitation and conservation. Here, this study reports detailed information on the complete chloroplast (cp) of L. bulbifera. The result indicates that the cp genome of L. bulbifera of 150,005 bp contains 126 genes, among them, 37 tRNA genes and 81 protein-coding genes. The analysis of repetition demonstrated that palindromic repeats are more frequent. In the meantime, 39 SSRs were also identified, the majority of which were mononucleotides Adenine-Thymine (A-T). Furthermore, we compared L. bulbifera with eight published Laportea plastomes, to explore highly polymorphic molecular markers. The analysis identified four hypervariable regions, including rps16, ycf1, trnC-GCA and trnG-GCC. According to the phylogenetic analysis, L. bulbifera was most closely related to Laportea canadensis (L. canadensis), and the molecular clock analysis speculated that the species originated from 1.8216 Mya. Overall, this study provides a more comprehensive analysis of the evolution of L. bulbifera from the perspective of phylogenetic and intrageneric molecular variation in the genus Laportea, which is useful for providing a scientific basis for further identification, taxonomic, and evolutionary studies of the genus.

Comparative and phylogenetic analyses of the chloroplast genome reveal the taxonomy of the Morus genus

Mulberry (genus Morus) is an economically important woody plant with an altered ploidy level. The variable number of Morus species recognized by different studies indicates that the genus is in need of revision. In this study, the chloroplast (CP) genomes of 123 Morus varieties were de novo assembled and systematically analyzed. The 123 varieties represented six Morus species, namely, Morus alba, Morus nigra, Morus notabilis, Morus rubra, Morus celtidifolia, and Morus serrata. The Morus CP genome was found to be 158,969~159,548 bp in size with 125 genes, including 81 protein coding, 36 tRNA, and 8 rRNA genes. The 87 out of 123 mulberry accessions were assigned to 14 diverse groups with identical CP genome, which indicated that they are maternally inherited and share 14 common ancestors. Then 50 diverse CP genomes occurred in 123 mulberry accessions for further study. The CP genomes of the Morus genus with a quadripartite structure have two inverted repeat (IR) regions (25,654~25,702 bp) dividing the circular genome into a large single-copy (LSC) region (87,873~88,243 bp) and small single-copy (SSC) region (19,740~19,994 bp). Analysis of the phylogenetic tree constructed using the complete CP genome sequences of Morus revealed a monophyletic genus and that M. alba consisted of two clades, M. alba var. alba and M. alba var. multicaulis. The Japanese cultivated germplasms were derived from M. alba var. multicaulis. We propose that the Morus genus be classified into six species, M. nigra, M. notabilis, M. serrata, M. celtidifolia, M. rubra, and M. alba with two subspecies, M. alba var. alba and M. alba var. multicaulis. Our findings provide a valuable resource for the classification, domestication, and breeding improvement of mulberry.

Evolutionary history of genus Coptis and its dynamic changes in the potential suitable distribution area

The genus Coptis belongs to the Ranunculaceae family, containing 15 recognized species highly diverse in morphology. It is a conspicuous taxon with special evolutionary position, distribution pattern and medicinal value, which makes it to be of great research and conservation significance. In order to better understand the evolutionary dynamics of Coptis and promote more practical conservation measures, we performed plastome sequencing and used the sequencing data in combination with worldwide occurrence data of Coptis to estimate genetic diversity and divergence times, rebuild biogeographic history and predict its potential suitable distribution area. The average nucleotide diversity of Coptis was 0.0067 and the hotspot regions with the highest hypermutation levels were located in the ycf1 gene. Coptis is most likely to have originated in North America and Japanese archipelago and has a typical Eastern Asian and North American disjunct distribution pattern, while the species diversity center is located in Mid-West China and Japan. The crown age of the genus is estimated at around 8.49 Mya. The most suitable climatic conditions for Coptis were as follows: precipitation of driest quarter > 25.5 mm, annual precipitation > 844.9 mm and annual mean temperature -3.1 to 19 °C. The global and China suitable area shows an upward trend in the future when emission of greenhouse gases is well controlled, but the area, especially in China, decreases significantly without greenhouse gas policy interventions. The results of this study provide a comprehensive insight into the Coptis evolutionary dynamics and will facilitate future conservation efforts.

Plastid phylogenomic insights into relationships, divergence, and evolution of Apiales

Members of Apiales are monophyletic and radiated in the Late Cretaceous. Fruit morphologies are critical for Apiales evolution and negative selection and mutation pressure play important roles in environmental adaptation. Apiales include many foods, spices, medicinal, and ornamental plants, but the phylogenetic relationships, origin and divergence, and adaptive evolution remain poorly understood. Here, we reconstructed Apiales phylogeny based on 72 plastid genes from 280 species plastid genomes representing six of seven families of this order. Highly supported phylogenetic relationships were detected, which revealed that each family of Apiales is monophyletic and confirmed that Pennanticeae is a member of Apiales. Genera Centella and Dickinsia are members of Apiaceae, and the genus Hydrocotyle previously classified into Apiaceae is confirmed to belong to Araliaceae. Besides, coalescent phylogenetic analysis and gene trees cluster revealed ten genes that can be used for distinguishing species among families of Apiales. Molecular dating suggested that the Apiales originated during the mid-Cretaceous (109.51 Ma), with the families' radiation occurring in the Late Cretaceous. Apiaceae species exhibit higher differentiation compared to other families. Ancestral trait reconstruction suggested that fruit morphological evolution may be related to shifts in plant types (herbaceous or woody), which in turn is related to the distribution areas and species numbers. Codon bias and positive selection analyses suggest that negative selection and mutation pressure may play important roles in environmental adaptation of Apiales members. Our results improve the phylogenetic framework of Apiales and provide insights into the origin, divergence, and adaptive evolution of this order and its members.

Phylogenomic analysis and development of molecular markers for the determination of twelve plum cultivars (Prunus, Rosaceae)

Background

Plums are one of the most important economic crops of the Rosaceae family and are produced all over the world. China has many local varieties, but the genomic information is limited for genetic studies. Here, we first sequenced, assembled, and analyzed the plastomes of twelve plum cultivars and developed molecular markers to distinguish them.

Results

The twelve plastomes of plum cultivars have a circular structure of 157,863–157,952 bp containing a large single-copy region (LSC) of 86,109–86,287 bp, a small copy region (SSC) of 18,927–19,031 bp, and two inverted repeats (IR) of 26,353–26,387 bp each. The plastomes of plum cultivars encode 131 genes, including 86 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. We detected 50, 54, 54, 53, 53, 50, 54, 54, 54, 49, 50, 54 SSRs in the twelve analyzed varieties, respectively. For repeat sequences, we identified 553 tandem repeats, 204 direct repeats, and 270 palindromic repeats. We also analyzed the expansion/contraction of IR regions. The genes rpl22, rps19, rpl2, ycf1, ndhF, and the trnH span on or near the boundary of IR and single-copy regions. Phylogenetic analysis showed that the twelve cultivars were clustered with the P. salicina and P. domestica. We developed eight markers LZ01 to LZ08 based on whole plastomes and nuclear genes and validated them successfully with six repetitions.

Conclusions

The results obtained here could fill in the blanks of the plastomes of these twelve plum cultivars and provide a wider perspective based on the basis of the plastomes of Prunus to the molecular identification and phylogenetic construction accurately. The analysis from this study provides an important and valuable resource for studying the genetic basis for agronomic and adaptive differentiation of the Prunus species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08965-z.

Danxiaorchismangdangshanensis (Orchidaceae, Epidendroideae), a new species from central Fujian Province based on morphological and genomic data

Danxiaorchismangdangshanensis, a new mycoheterotrophic species from Fujian Province, China, is described and illustrated. The new species is morphologically similar to D.singchiana, but its callus of labellum is a less distinctive Y-shape with three auricles on the apex, four pollinia that are narrowly elliptic in shape and equal in size, and it lacks fine roots. The plastome of D.mangdangshanensis is highly degraded. Phylogenetic analyses distinguished D.mangdangshanensis from its congeners, D.singchiana and D.yangii, with strong support based on nrITS + matK and plastomes, respectively.

Complete chloroplast genome sequences of Phlomis fruticosa and Phlomoides strigosa and comparative analysis of the genus Phlomis sensu lato (Lamiaceae)

The taxonomic terms "Phlomis" and "Phlomoides" had been used to describe two sections within the genus Phlomis belonging to the Lamiaceae family. Recently, phylogenetic analyses using molecular markers showed that Phlomis and Phlomoides formed two monophyletic clades, and thus they are generally accepted as separate genera. In this study, we assembled the complete chloroplast genome of Phlomis fruticosa, which is the first reported chloroplast genome belonging to Phlomis genus, as well as the complete chloroplast genome of Phlomoides strigosa belonging to Phlomoides genus. The results showed that the length of chloroplast genome was 151,639 bp (Phlomis fruticosa) and 152,432 bp (Phlomoides strigosa), with conserved large single copy regions, small single copy regions, and inverted repeat regions. 121 genes in Phlomis fruticosa and 120 genes in Phlomoides strigosa were annotated. The chloroplast genomes of Phlomis fruticosa, Phlomoides strigosa, and three reported Phlomoides species, as well as those of 51 species from the Lamiaceae family, which covered 12 subfamilies, were subjected to phylogenetic analyses. The Phlomis and Phlomoides species were split into two groups, which were well supported by both maximum likelihood and Bayesian inference tree analyses. Our study provided further evidence to recognize Phlomis and Phlomoides as independent genera.

De novo assembly and characterization of the complete chloroplast genome of Elymus magellanicus (É.Desv.) Á.Löve, 1984 (Poaceae, Pooideae)

Elymus magellanicus (É.Desv.) Á.Löve is a foliage accent plant that originated in South America. In this study, the complete chloroplast genome of E. magellanicus is reported. It was found to have a total size of 133,249 bp. The chloroplast genome was found to consist of two inverted repeats (IRA and IRB) of 21,421 bp each, a small single-copy region of 12,709 bp, and a large single-copy region (77,698 bp). The annotation results show the GC content of the chloroplast genome to be 38.47%, including 40 tRNA genes, 82 protein-coding genes, and 8 rRNA genes. Phylogenetic analysis of 29 species revealed that E. magellanicus is closely related to E. arenarius.

The chloroplast genome of Spiraea thunbergii (Rosaceae)

Spiraea thunbergii (S. thunbergia) is a very common ornamental shrub in the garden, with important horticultural and economic value. In this study, we assembled the complete chloroplast (cp) genome of S. thunbergia into a typical quadripartite structure. The genome size and GC content of the S. thunbergia cp genome are 155,922 bp and 36.76%, respectively. It contains 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The phylogenetic tree supported that S. thunbergii is closely related to Spiraea mongolica in the Rosaceae family. The study will provide significant genomic resource for elucidating the phylogenetic relationship of Spiraea.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

Complete chloroplast genome of Leptopus chinensis (Bunge) Pojark (Euphorbiaceae), a traditional Chinese herbal medicine

Leptopus chinensis (Bunge) Pojark is a medicinal herb endemic to China. In this study, its complete chloroplast (cp) genome was characterized with a discussion of its phylogenetic placement. The cp genome of L. chinensis was 154,600 bp long, with a double-stranded circular tetrad structure containing a small single-copy (SSC) region (17,717 bp), a large single-copy (LSC) region (83,159 bp), and a pair of inverted repeat (IR) regions (26,862 bp each). The overall GC content of the genome was 36.8% (LSC, 34.7%; SSC, 29.8%; IR, 42.3%). Phylogenetic analysis indicated that L. chinensis is a sister species to L. cordifolius.

The complete chloroplast genome of Salix matsudana f. tortuosa

In this study, the complete chloroplast genome of Salix matsudana f. tortuosa was sequenced and analyzed. The genome of Salix matsudana f. tortuosa was 155,673 bp in length and was quadripartite in structure, containing a large single-copy region with a length of 84,447 bp, a small single-copy region with a length of 16,320 bp, and two inverted repeats of 27,453 bp in length. The chloroplast genome contains 130 genes, including 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The GC content is 36.64%. The phylogenetic tree shows that Salix matsudana f. tortuosa, Salix matsudana, and Salix babylonica are closely related and located on the same branch. The chloroplast genome of Salix matsudana f. tortuosa will provide important data for further systematic study of Salicaceae and the genus Salix.

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.

Chloroplast genome characterization of Rubus arcticus L

Rubus arcticus Linnaeus (1753) is a medicinal and edible plant in the Rosaceae with wide distribution in northeast China. The total length of the genome was 156,668 bp with a GC content of 37.1%, including a large single-copy (LSC, 85,958 bp) region, a small single-copy region (SSC, 18,756 bp), and inverted repeat (IR, 51,954 bp) regions. A total of 129 genes were identified. The numbers of protein genes tRNAs and rRNAs were 85, 36, and 8, respectively. Phylogenetic analysis indicated that R. arcticus belongs to the Rubus genus. Published R. arcticus chloroplast genomes have yielded insights into the closely related species identification, phylogenetic position and Rubus evolution.

The chloroplast genome of the Iris japonica Thunberg (Butterfly flower) reveals the genomic and evolutionary characteristics of Iris species

Iris japonica Thunberg is one of the horticultural species belonging to the Iris genus and Iridaceae family. Previous studies have revealed its hepatoprotective activity and ornamental values. However, little genetic and genomic information about this species is available. Here, to decipher the chloroplast genome and reveal its evolutionary characteristics, we sequenced, de novo assembled, and comprehensively analyzed the chloroplast genome of I. japonica. The genome was 152,453 bp in length and displayed a circular structure with a large single-copy region, a small single-copy region, and two inverted repeat regions. It contained 131 genes, including 85 protein-coding genes, eight ribosomal RNA genes, and 38 transfer RNA genes. We also identified 23 microsatellite repeat sequences, 34 tandem repeat sequences, and 60 dispersed repeat sequences in the chloroplast genome of I. japonica. Sequence divergence analyses of the chloroplast genomes of 20 Iris species revealed that the top four most highly variable regions were ndhC-trnV-UAC, rpl22-rps19, rps16-trnQ-UUG, and trnG-UCC-trnR-UCU. Phylogenetic analysis showed that I. japonica was most closely related to I. tectorum. This study reported a new chloroplast genome of I. japonica and performed comparative analyses of 20 Iris chloroplast genomes. The results would facilitate the evolutionary research and development of molecular markers for Iris species.

Comparative Genomics and Phylogenetic Analysis of the Chloroplast Genomes in Three Medicinal Salvia Species for Bioexploration

To systematically determine their phylogenetic relationships and develop molecular markers for species discrimination of Salvia bowleyana, S. splendens, and S. officinalis, we sequenced their chloroplast genomes using the Illumina Hiseq 2500 platform. The chloroplast genomes length of S. bowleyana, S. splendens, and S. officinalis were 151,387 bp, 150,604 bp, and 151,163 bp, respectively. The six genes ndhB, rpl2, rpl23, rps7, rps12, and ycf2 were present in the IR regions. The chloroplast genomes of S. bowleyana, S. splendens, and S. officinalis contain 29 tandem repeats; 35, 29, 24 simple-sequence repeats, and 47, 49, 40 interspersed repeats, respectively. The three specific intergenic sequences (IGS) of rps16-trnQ-UUG, trnL-UAA-trnF-GAA, and trnM-CAU-atpE were found to discriminate the 23 Salvia species. A total of 91 intergenic spacer sequences were identified through genetic distance analysis. The two specific IGS regions (trnG-GCC-trnM-CAU and ycf3-trnS-GGA) have the highest K2p value identified in the three studied Salvia species. Furthermore, the phylogenetic tree showed that the 23 Salvia species formed a monophyletic group. Two pairs of genus-specific DNA barcode primers were found. The results will provide a solid foundation to understand the phylogenetic classification of the three Salvia species. Moreover, the specific intergenic regions can provide the probability to discriminate the Salvia species between the phenotype and the distinction of gene fragments.

Characterization of the complete chloroplast genome of Clematis orientalis (Ranunculaceae)

The medicinal plant Clematis orientalis L. belongs to the family Ranunculaceae. In this study, we determined the complete chloroplast genome sequence of C. orientalis and its phylogenetic relationships with other species. The complete chloroplast genome of C. glauca is 159,518 bp in length, circular in structure, and has four regions including a large single-copy (LSC) region of 79,453 bp; a small single-copy (SSC) region of 17,997 bp; and two inverted repeat (IR) regions of 31,034 bp. The GC content of the genome is 38%, and those of LSC, SSC, and IR regions are 36.2, 31.4, and 42%, respectively. The genome encodes 129 unique genes, including 85 protein-coding genes, 36 tRNA genes, and eight rRNA genes. Phylogenomic analysis reveals that C. orientalis is most closely related to C. aethusifolia. This study contributes to better understanding of phylogenetic relationships of Ranunculaceae.

Characterization of the complete Salix viminalis var. gmelinii Turcz 1854 chloroplast genome from the northeast of China

The Salix viminalis var. gmelinii Turcz 1854 is a variant of the Salix genus, from the Salicaceae family, and possesses an extremely high economic value. In this study the complete chloroplast genome of the woody plant S. viminalis var. gmelinii was characterized for the first time using a high-throughput approach in conjunction with de novo assembly technology. The S.viminalis var. gmelinii chloroplast genome is 155,405 base pairs (bp) in length and contains 36.71% GC content. It incorporates a large single-copy region (LSC, 84,287bp) alongside one small-copy region (SSC, 16,198bp), and two inverted repeat regions (IRA and IRB, 27,460bp). Moreover, this chloroplast genome encodes 128 genes, which comprises 83 protein-coding genes, 37 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. Furthermore, the phylogenetic analysis revealed that S.viminalis var. gmelinii is closely related to S. cupularis and S.gordejevii.

Characterization and Comparative Analysis of Chloroplast Genomes in Five Uncaria Species Endemic to China

Uncaria, a perennial vine from the Rubiaceae family, is a typical Chinese traditional medicine. Currently, uncertainty exists over the Uncaria genus’ evolutionary relationships and germplasm identification. The complete chloroplast genomes of four Uncaria species mentioned in the Chinese Pharmacopoeia and Uncaria scandens (an easily confused counterfeit) were sequenced and annotated. The findings demonstrated that the whole chloroplast genome of Uncaria genus is 153,780–155,138 bp in full length, encoding a total of 128–131 genes, containing 83–86 protein-coding genes, eight rRNAs and 37 tRNAs. These regions, which include eleven highly variable loci and 31–49 SSRs, can be used to create significant molecular markers for the Uncaria genus. The phylogenetic tree was constructed according to protein-coding genes and the whole chloroplast genome sequences of five Uncaria species using four methods. The topology of the two phylogenetic trees showed no difference. The sequences of U. rhynchophylla and U. scandens are clustered in one group, while the U. hirsuta and U. macrophylla are clustered in another group. U. sessilifructus is clustered together with the above two small clades. New insights on the relationship were revealed via phylogenetic research in five Uncaria species. This study will provide a theoretical basis for identifying U. rhynchophylla and its counterfeits, as well as the species of the Uncaria genus. This research provides the initial chloroplast genome report of Uncaria, contributes to elucidating the chloroplast genome evolution of Uncaria in China.

Comparative analysis of the chloroplast genomes of eight Piper species and insights into the utilization of structural variation in phylogenetic analysis

With more than 2000 species, Piper is regarded as having high medicinal, cosmetic, and edible value. There also remain some taxonomic and evolutionary uncertainties about the genus. This study performed chloroplast genome sequencing of eight poorly studied Piper species and a comparative analysis with black pepper (Piper nigrum). All examined species were highly similar in gene content, with 79 protein-coding genes, 24 tRNAs, and four rRNAs. They also harbored significant structural differences: The number of SSRs ranged from 63 to 87, over 10,000 SNPs were detected, and over 1,000 indels were found. The spatial distribution of structural differences was uneven, with the IR and LSC being relatively more conserved and the SSC region highly variable. Such structural variations of the chloroplast genome can help in evaluating the phylogenetic relationships between species, deciding some hard-to-distinguish evolutionary relationships, or eliminating improper markers. The SSC region may be evolving at high speed, and some species showed a high degree of sequence variation in the SSC region, which seriously affected marker sequence detection. Conversely, CDS sequences tended to lack variation, and some CDSs can serve as ideal markers for phylogenetic reconstruction. All told, this study provides an effective strategy for selecting chloroplast markers, analyzing difficult-to-distinguish phylogenetic relationships and avoiding the taxonomic errors caused by high degree of sequence variations.

Chloroplast Genome Annotation Tools: Prolegomena to the Identification of Inverted Repeats

The development of next-generation sequencing technology and the increasing amount of sequencing data have brought the bioinformatic tools used in genome assembly into focus. The final step of the process is genome annotation, which works on assembled genome sequences to identify the location of genome features. In the case of organelle genomes, specialized annotation tools are used to identify organelle genes and structural features. Numerous annotation tools target chloroplast sequences. Most chloroplast DNA genomes have a quadripartite structure caused by two copies of a large inverted repeat. We investigated the strategies of six annotation tools (Chloë, Chloroplot, GeSeq, ORG.Annotate, PGA, Plann) for identifying inverted repeats and analyzed their success using publicly available complete chloroplast sequences of taxa belonging to the asterid and rosid clades. The annotation tools use two different approaches to identify inverted repeats, using existing general search tools or implementing stand-alone solutions. The chloroplast sequences studied show that there are different types of imperfections in the assembled data and that each tool performs better on some sequences than the others.

Complete chloroplast genomes of two medicinal Swertia species: the comparative evolutionary analysis of Swertia genus in the Gentianaceae family

The complete chloroplast genome of Swertia kouitchensis has been sequenced and assembled, compared with that of S. bimaculata to determine the evolutionary relationships among species of the Swertia in the Gentianaceae family. Swertia kouitchensis and S. bimaculata are from the Gentianaceae family. The complete chloroplast genome of S. kouitchensis was newly assembled, annotated, and analyzed by Illumina Hiseq 2500 platform. The chloroplast genomes of the two species encoded a total of 133, 134 genes, which included 88-89 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA genes. One intron was contained in each of the eight protein-coding genes and eight tRNA-coding genes, whereas two introns were found in two genes (ycf3 and clpP). The most abundant codon of the two species was for isoleucine, and the least abundant codon was for cysteine. The number of microsatellite repeat sequences was twenty-eight and thirty-two identified in the chloroplast genomes of S. kouitchensis and S. bimaculata, respectively. A total of 1127 repeat sequences were identified in all the 23 Swertia chloroplast genomes, and they fell into four categories. Furthermore, five divergence hotspot regions can be applied to discriminate these 23 Swertia species through genomes comparison. One pair of genus-specific DNA barcodes primer has been accurately identified. Therefore, the diverse regions cloned by a specific primer may become an effective and powerful molecular marker for the identification of Swertia genus. Moreover, four genes (ccsA, ndhK, rpoC1, and rps12) were positive selective pressure. The phylogenetic tree showed that the 23 Swertia species were clustered into a large clade including four evident subbranches, whereas the two species of S. kouitchensis and S. bimaculata were separately clustered into the diverse but correlated species group.

Characterization of the complete chloroplast genome of Hevea pauciflora (Euphorbiaceae), an important wild relative of the rubber tree

Hevea pauciflora belongs to the Euphorbiaceae family, an important wild relative of the rubber tree. This study sequenced, assembled, and annotated the complete chloroplast genome of H. pauciflora. The complete chloroplast genome is 161,123 bp with a canonical quadripartite structure containing a large single-copy (LSC) region (89,109 bp), a small single-copy (SSC) region (18,376 bp), and two inverted repeat regions (IRa and IRb) (26,819 bp, each). A total of 134 genes were annotated, including 86 protein-coding genes, four pseudogenes, 36 tRNA genes, and eight rRNA genes. The 134 genes include four major groups: ‘self-replication’, ‘photosynthesis’, ‘unknown function’, and ‘others’. A phylogenetic analysis clustered H. pauciflora, H. brasiliensis, H. camargoana, and H. benthamiana into one clade, consistent with traditional taxonomy. This study provides useful data for further studies of Hevea genus and the phylogenetic relationships of Euphorbiaceae species.

Characterization of the complete chloroplast genome of the Helleborus atrorubens Waldst. & Kit. (Ranunculaceae)

Helleborus atrorubens is an economically important perennial garden plant with medicinal value. Here, we sequenced the complete chloroplast genome of H. atrorubens. The results revealed the chloroplast genome to be 166,695 bp in length. It possesses a typical quadripartite structure containing one large single copy (LSC) region (84994 bp), one small single copy (SSC) region (17,825 bp), and a pair of inverted repeat (IR) regions (31938 bp). This chloroplast genome encoded 130 genes, out of which 85 code for proteins, 37 for transfer RNAs, and 8 for ribosomal RNAs. Simple sequence repeat (SSR) markers and the top variable coding regions were identified.

Our study lays a foundation for further research, such as species differentiation and phylogenetic reconstruction of the Helleborus genus.

The complete chloroplast genome of Zanthoxylum stenophyllum Hemsl. (Rutaceae), a traditional Chinese medicinal plant

Zanthoxylum stenophyllum Hemsl., a type species for the genus Zanthoxylum (Rutaceae), is a traditional medicinal plant. We studied the complete chloroplast genome of this species using BGISEQ-500 platform. The chloroplast genome was 158,314 bp in size with a GC content of 38.45%. The genome contained two short inverted repeat (IRa and IRb) regions of 27,052 bp, a large single-copy region (LSC, 86,029 bp) and a small single-copy region (SSC, 18,181 bp). The annotated complete chloroplast genome contains 133 distinct genes, including 88 protein-coding genes, 37 transfer RNAs (tRNAs), and 8 ribosomal RNAs (rRNAs). Phylogenetic analysis indicated that Z. stenophyllum is clustered with Z. schiniflium and Z. pinnatum in the same branch with 100% bootstrap support. This complete chloroplast genome provides valuable genomic information for the molecular phylogeny and sustainable utilization of Zanthoxylum.

The complete chloroplast genome sequence of Litsea coreana var. lanuginosa (Lauraceae): genome structure and phylogenetic analysis

Litsea coreana var. lanuginose is a perennial, indeciduous, and broad-leaved tree used as an essential medicinal and edible plant. In addition, this species is well-known for its leaves are rich in aromatic oil. In this study, we firstly assembled and characterized the complete chloroplast genome of L. coreana var. lanuginose using Illumina pair-end sequencing and performed a phylogenetic analysis with other 13 species in Lauraceae. The results revealed that its chloroplast genome was 152,859 bp in total length with 39% of GC content, containing a pair of inverted repeats of 20,084 bp (IRA and IRB), separated by a large single-copy (LSC) region of 93,795 bp and a small single-copy (SSC) region of 18,896 bp. The plastid genome of L. coreana var. lanuginose encoded 125 genes, including 81 protein-coding genes, 36 transfer RNA (tRNA), and eight ribosomal RNA (rRNA) genes. The phylogenetic analysis suggested that L. coreana var. lanuginose was closely related to the clade of Litsea monopetala, Litsea garrettii, and Litsea elongate in Lauraceae family.

Phylogenetic analysis and complete chloroplast genome of Salix wilhelmsiana M.B

Salix wilhelmsiana M.B. Bieberstein 1819 is a perennial woody plant with high economic and ecological value. In this study, we annotated the chloroplast (cp) genome of Salix wilhelmsiana M.B. The results showed that the length of the complete cp genome is 155,577 bp, which is typically composed of two single-copy regions (large single-copy (LSC) of 84,439 bp and small single-copy (SSC) of 16,221 bp) and a pair of IR regions of 27,457 bp with a quadripartite structure. The genome contains 129 genes, including 84 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The GC content was 36.70%. Phylogenetic analysis based on cp genome sequences of 19 species from the Salicaceae family revealed that S. wilhelmsiana M.B. is closely related to S. viminalis var. gmelinii.