Automatic annotation of organellar genomes with DOGMA

Loss of the IR region in conifer plastomes: Changes in the selection pressure and substitution rate of protein-coding genes

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein-coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non-conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non-conifers. It is more obvious that in non-conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6- to 3.1-fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2- to 3.6-fold of non-conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein-coding genes due to gene function, plant habitat, or newly acquired IRs.

Complete Chloroplast Genome Analysis of Two Important Medicinal Alpinia Species: Alpinia galanga and Alpinia kwangsiensis

Most Alpinia species are valued as foods, ornamental plants, or plants with medicinal properties. However, morphological characteristics and commonly used DNA barcode fragments are not sufficient for accurately identifying Alpinia species. Difficulties in species identification have led to confusion in the sale and use of Alpinia for medicinal use. To mine resources and improve the molecular methods for distinguishing among Alpinia species, we report the complete chloroplast (CP) genomes of Alpinia galanga and Alpinia kwangsiensis species, obtained via high-throughput Illumina sequencing. The CP genomes of A. galanga and A. kwangsiensis exhibited a typical circular tetramerous structure, including a large single-copy region (87,565 and 87,732 bp, respectively), a small single-copy region (17,909 and 15,181 bp, respectively), and a pair of inverted repeats (27,313 and 29,705 bp, respectively). The guanine-cytosine content of the CP genomes is 36.26 and 36.15%, respectively. Furthermore, each CP genome contained 133 genes, including 87 protein-coding genes, 38 distinct tRNA genes, and 8 distinct rRNA genes. We identified 110 and 125 simple sequence repeats in the CP genomes of A. galanga and A. kwangsiensis, respectively. We then combined these data with publicly available CP genome data from four other Alpinia species (A. hainanensis, A. oxyphylla, A. pumila, and A. zerumbet) and analyzed their sequence characteristics. Nucleotide diversity was analyzed based on the alignment of the complete CP genome sequences, and five candidate highly variable site markers (trnS-trnG, trnC-petN, rpl32-trnL, psaC-ndhE, and ndhC-trnV) were found. Twenty-eight complete CP genome sequences belonging to Alpinieae species were used to construct phylogenetic trees. The results fully demonstrated the phylogenetic relationship among the genera of the Alpinieae, and further proved that Alpinia is a non-monophyletic group. The complete CP genomes of the two medicinal Alpinia species provides lays the foundation for the use of CP genomes in species identification and phylogenetic analyses of Alpinia species.

Chloroplast genome sequencing based on genome skimming for identification of Eriobotryae Folium

Background

Whole chloroplast genome (cpDNA) sequence is becoming widely used in the phylogenetic studies of plant and species identification, but in most cases the cpDNA were acquired from silica gel dried fresh leaves. So far few reports have been available to describe cpDNA acquisition from crude drugs derived from plant materials, the DNA of which usually was seriously damaged during their processing. In this study, we retrieved cpDNA from the commonly used crude drug Eriobotryae Folium (Pipaye in Chinese, which is the dried leaves of Eriobotrya japonica, PPY) using genome skimming technique.

Results

We successfully recovered cpDNA sequences and rDNA sequences from the crude drug PPY, and bioinformatics analysis showed a high overall consistency between the cpDNA obtained from the crude drugs and fresh samples. In the ML tree, each species formed distinct monophyletic clades based on cpDNA sequence data, while the phylogenetic relationships between Eriobotrya species were poorly resolved based on ITS and ITS2.

Conclusion

Our results demonstrate that both cpDNA and ITS/ITS2 are effective for identifying PPY and its counterfeits derived from distantly related species (i.e. Dillenia turbinata and Magnolia grandiflora), but cpDNA is more effective for distinguishing the counterfeits derived from the close relatives of Eriobotrya japonica, suggesting the potential of genome skimming for retrieving cpDNA from crude drugs used in Traditional Chinese Medicine for their identification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12896-021-00728-0.

Chloroplast genomes of two Pueraria DC. species: sequencing, comparative analysis and molecular marker development

Puerariae lobatae radix (Ge‐Gen in Chinese) and Puerariae thomsonii radix (Fen‐Ge) are widely used as medicine and health products, particularly in Chinese medicine. Puerarin and daidzein are the primary bioactive compounds in Puerariae radix. These isoflavones have been used to treat cardiovascular and cerebrovascular diseases, hypertension, diabetes, and osteoporosis. The content of puerarin in Ge‐Gen is about six times higher than that in Fen‐Ge, so its use has a higher pharmacological effect. It is therefore of great importance to effectively distinguish between these two species. However, because their basal plants, P. lobata (Willd.) Ohwi and P. thomsonii Benth., possess an extremely similar appearance, and detecting the level of chemical constituents is just a rough distinction, it is necessary to develop more efficient identification approaches. Here the complete chloroplast genomes of P. lobata and P. thomsonii were deciphered, including sequencing, assembly, comparative analysis, and molecular marker development. The results showed that they are 153,393 and 153,442 bp in length, respectively; both contain 124 annotated genes, including eight encoding rRNA, 29 encoding tRNA, and 87 encoding proteins. Phylogenetic analysis showed that they form a clade, indicating that they originate from the same ancestor. After obtaining 10 intergenic/intronic regions with a genetic distance greater than 0.5 cm, primers were designed to amplify regions of high variability in P. lobata and P. thomsonii. Finally, a 60‐bp differential base fragment, located in the intron of rpl16, was developed as a molecular marker to efficiently distinguish between these two species.

The complete mitochondrial genome of Pycanum ochraceum Distant 1893 (Hemiptera: Tessaratomidae)

In the present study, the complete mitochondrial genome of Pycanum ochraceum was identified for the first time. The entire genome is 17,198 bp in length with 73.6% A + T content. It contains 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), 13 protein-coding genes (PCGs) and 1 noncoding control region (D-loop). Phylogenetic analysis showed that Tessaratomidae bugs are monophyletic. This study can provide essential DNA molecular data for further phylogenetic and evolutionary analysis for Heteroptera.

Complete chloroplast genome of Baccaurea ramiflora (phyllanthaceae), a promising underutilized species

Baccaurea ramiflora is a high economic value for the exploitation and utilization of wild fruit tree species with edible and medicinal values in Southeast Asia. In this study, the complete chloroplast genome sequence of B. ramiflora was assembled and characterized. The chloroplast genome was 161,089 bp in length, consisting of a large single copy (LSC) of 89,515 bp and a small single copy (SSC) region of 18,826 bp, which were separated by a pair of 26,374 bp inverted repeat (IR) regions. The overall GC content was 36.71%. A total of 132 genes, including 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes were identified. A low intraspecies variation within B. ramiflora was found with 22 single nucleotide polymorphisms (SNPs) and 17 insertions and deletions (INDELs). Phylogenetic tree reconstructed by 14 chloroplast genomes revealed that B. ramiflora clusters together with species of Phyllanthus, Glochidion, and Flueggea.

The complete plastome of a cultivar of Lannea coromandelica

Lannea coromandelica (Houtt.) Merr. is a deciduous tree in the family Anacardiaceae, which grows in lowland and hill forests; 100–1800 m. SW Guangdong, S Guangxi, S Yunnan [Bhutan, India, Myanmar, Nepal, Sri Lanka; cultivated elsewhere in continental SE Asia, such as in Cambodia, Laos, Malaysia, Thailand, Vietnam, where it is probably naturalized]. The length of the complete plastome is 162,460 bp, including 130 genes consisting of 85 protein-coding genes, 37 tRNA genes and 8 rRNA genes. The assembled plastome has the typical structure and gene content of angiosperms plastome, which includes two inverted repeats (IRs) regions of 26,877 bp, a large single copy (LSC) region of 89,599 bp and a small single-copy (SSC) region of 19,107 bp. The total G/C content in the plastome of L. coromandelica is 37.7%. The complete plastome sequence of L. coromandelica will provide contributions to the conservation genetics of this species as well as to phylogenetic studies in Anacardiaceae.

Plastid phylogenomics provides novel insights into the infrafamilial relationship of Polypodiaceae

The polygrammoids (Polypodiaceae) are the most species-rich and diversified epiphytic fern lineages, and hold an important role to understand the deep diverging events and rapid adaptation to changing environments in the plant tree of life. Despite progress in the phylogeny of this group of ferns in previous multilocus phylogenetic studies, uncertainty remains especially in backbone relationships among closely related clades, and the phylogenetic placement of recalcitrant species or lineages. Here, we investigated the deep phylogenetic relationships within Polypodiaceae by sampling all major lineages and using 81 plastid genomes (plastomes), of which 70 plastomes were newly sequenced with high-throughput sequencing technology. Based on parsimony, maximum-likelihood, Bayesian and multispecies coalescent analyses of genome skimming data, we achieved a better resolution of the backbone phylogeny of Polypodiaceae. Using simulated data matrices, we detected that potential phylogenetic artefacts, such as long-branch attraction and insufficient taxonomic sampling, may have a confounding impact on the incongruence of phylogenetic inferences. Furthermore, our phylogenetic analyses offer greater resolution than previous multilocus studies, providing a robust framework for future phylogenetic implications on the subfamilial taxonomy of Polypodiaceae. Our phylogenomic study not only demonstrates the advantage of a character-rich plastome dataset for resolving the recalcitrant lineages that have undergone rapid radiation, but also sheds new light on integrative explorations understanding the evolutionary history of large fern groups in the genomic era.

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

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

Molecular evolution of chloroplast genomes in subfamily Zingiberoideae (Zingiberaceae)

BACKGROUND

Zingiberoideae is a large and diverse subfamily of the family Zingiberaceae. Four genera in subfamily Zingiberoideae each possess 50 or more species, including Globba (100), Hedychium (> 80), Kaempferia (50) and Zingiber (150). Despite the agricultural, medicinal and horticultural importance of these species, genomic resources and suitable molecular markers for them are currently sparse.

RESULTS

Here, we have sequenced, assembled and analyzed ten complete chloroplast genomes from nine species of subfamily Zingiberoideae: Globba lancangensis, Globba marantina, Globba multiflora, Globba schomburgkii, Globba schomburgkii var. angustata, Hedychium coccineum, Hedychium neocarneum, Kaempferia rotunda 'Red Leaf', Kaempferia rotunda 'Silver Diamonds' and Zingiber recurvatum. These ten chloroplast genomes (size range 162,630-163,968 bp) possess typical quadripartite structures that consist of a large single copy (LSC, 87,172-88,632 bp), a small single copy (SSC, 15,393-15,917 bp) and a pair of inverted repeats (IRs, 29,673-29,833 bp). The genomes contain 111-113 different genes, including 79 protein coding genes, 28-30 tRNAs and 4 rRNA genes. The dynamics of the genome structures, gene contents, amino acid frequencies, codon usage patterns, RNA editing sites, simple sequence repeats and long repeats exhibit similarities, with slight differences observed among the ten genomes. Further comparative analysis of seventeen related Zingiberoideae species, 12 divergent hotspots are identified. Positive selection is observed in 14 protein coding genes, including accD, ccsA, ndhA, ndhB, psbJ, rbcL, rpl20, rpoC1, rpoC2, rps12, rps18, ycf1, ycf2 and ycf4. Phylogenetic analyses, based on the complete chloroplast-derived single-nucleotide polymorphism data, strongly support that Globba, Hedychium, and Curcuma I + "the Kaempferia clade" consisting of Curcuma II, Kaempferia and Zingiber, form a nested evolutionary relationship in subfamily Zingiberoideae.

CONCLUSIONS

Our study provides detailed information on ten complete Zingiberoideae chloroplast genomes, representing a valuable resource for future studies that seek to understand the molecular evolutionary dynamics in family Zingiberaceae. The identified divergent hotspots can be used for development of molecular markers for phylogenetic inference and species identification among closely related species within four genera of Globba, Hedychium, Kaempferia and Zingiber in subfamily Zingiberoideae.

Authentication of Zingiber Species Based on Analysis of Metabolite Profiles

Zingiber corallinum and Zingiber montanum, which belong to the Zingiberaceae family, are traditional Chinese folk medicinal herbs in Guizhou and Yunnan Province of China. They share great similarities in morphology, chemical constituent, and DNA barcoding sequence. The taxonomy of the two Zingiber species is controversial and discrimination of traditional Chinese medicines directly affects the pharmacological and clinical effects. In the present study, we performed a systemic analysis of "super-barcode" and untargeted metabolomics between Z. corallinum and Z. montanum using chloroplast (cp) genome sequencing and gas chromatography-mass spectrometry (GC-MS) analysis. Comparison and phylogenetic analysis of cp genomes of the two Zingiber species showed that the cp genome could not guarantee the accuracy of identification. An untargeted metabolomics strategy combining GC-MS with chemometric methods was proposed to distinguish the Zingiber samples of known variety. A total of 51 volatile compounds extracted from Z. corallinum and Z. montanum were identified, and nine compounds were selected as candidate metabolic markers to reveal the significant difference between Z. corallinum and Z. montanum. The performance of the untargeted metabolomic approach was verified with unknown Zingiber samples. Although the cp genomes could not be used to identify Zingiber species in this study, it will still provide a valuable genomics resource for population studies in the Zingiberaceae family, and the GC-MS based metabolic fingerprint is more promising for species identification and safe application of Z. corallinum and Z. montanum.

The complete chloroplast genome sequence of the medicinal plant Cercis chinensis and phylogenetic analysis

Cercis chinensis is a deciduous shrub or small landscape tree with exceptional ornamental characteristics. Here, we report the complete chloroplast genome of C. chinensis to provide a foundation for further phylogenetic studies of the Leguminosae. The chloroplast (cp) genome was 158,999 bp in size containing a large single-copy (LSC) region 88,176 bp in length, a small single-copy (SSC) region 19,583 bp in length, and two inverted repeat (IR) regions that were 25,620 bp each. The total GC content of the cp genome was 42.9% with the LSC, SSC, and IR regions 36.2, 33.8, and 29.4%, respectively. The cp genome contains 129 genes, including 85 protein-coding, 36 tRNA, and 8 rRNA genes. The phylogenetic analysis revealed that C. chinensis is closely related to Cercis glabra. These results provide valuable information about the evolutionary processes and conservation genetics of C. chinensis.

Evolutionary characteristics and phylogeny of cotton chloroplast tRNAs

The novel structural variations were identified in cotton chloroplast tRNAs and gene loss events were more obvious than duplications in chloroplast tRNAs. Transfer RNAs (tRNA) have long been believed an evolutionary-conserved molecular family, which play the key roles in the process of protein biosynthesis in plant life activities. In this study, we detected the evolutionary characteristics and phylogeny of chloroplast tRNAs in cotton plants, an economic and fibered important taxon in the world. We firstly annotated the chloroplast tRNAs of 27 Gossypium species to analyze their genetic composition, structural characteristics and evolution. Compared with the traditional view of evolutionary conservation of tRNA, some novel tRNA structural variations were identified in cotton plants. I.g., tRNA^Val-UAC and tRNA^Ile-GAU only contained one intron in the anti-condon loop region of tRNA secondary structure, respectively. In the variable region, some tRNAs contained a circle structure with a few nucleotides. Interestingly, the calculation result of free energy indicated that the variation of novel tRNAs contributed to the stability of tRNA structure. Phylogenetic analysis suggested that chloroplast tRNAs have evolved from multiple common ancestors, and the tRNA^Met seemed to be an ancestral tRNA, which can be duplicated and diversified to produce other tRNAs. The chloroplast tRNAs contained a group I intron in cotton plants, and the evolutionary analysis of introns indicated that group I intron of chloroplast tRNA originated from cyanobacteria. Analysis of gene duplication and loss events showed that gene loss events were more obvious than duplications in Gossypium chloroplast tRNAs. Additionally, we found that the rate of transition was higher than the ones of transversion in cotton chloroplast tRNAs. This study provided new insights into the structural characteristics and evolution of chloroplast tRNAs in cotton plants.

Sequencing of the Complete Mitochondrial Genome of Pingus sinensis (Spirurina: Quimperiidae): Gene Arrangements and Phylogenetic Implications

Despite several decades of intensive research on spirurine nematodes, molecular data on some of the main lineages are still absent, which makes taxonomic classification insufficiently resolved. In the present study, we sequenced the first complete mitogenome for the family Quimperiidae, belonging to P. sinensis (Spirurina: Quimperiidae), a parasite living in the intestines of snakehead (Ophiocephalus argus). The circular mitogenome is 13,874 bp long, and it contains the standard nematode gene set: 22 transfer RNAs, 2 ribosomal RNAs and 12 protein-coding genes. There are also two long non-coding regions (NCR), in addition to only 8 other intergenic regions, ranging in size from 1 to 58 bp. To investigate its phylogenetic position and study the relationships among other available Spirurina, we performed the phylogenetic analysis using Bayesian inference and maximum likelihood approaches by concatenating the nucleotide sequences of all 36 genes on a dataset containing all available mitogenomes of the suborder Spirurina from NCBI and compared with gene order phylogenies using the MLGO program. Both supported the closer relationship of Ascaridoidea to Seuratoidea than to Spiruroidea. Pingus formed a sister-group with the Cucullanus genus. The results provide a new insights into the relationships within Spirurina.

The complete chloroplast genome of Clerodendrum lindleyi Decne. ex Planch. (Tubiflorae: Verbenaceae) from nanjing, China

Clerodendrum lindleyi Decne. ex Planch. is a Chinese medicinal plant in the Lingnan region of China. In this study, the complete chloroplast genome sequence of C. lindleyi was assembled and characterized from high-throughput sequencing data. The chloroplast genome is 151,678 bp in length, consisting of a large single-copy (LSC) and a small single-copy (SSC) regions of 83,043 bp and 17,311 bp, respectively, which are separated by a pair of 25,662 bp inverted repeat (IR) regions. The overall GC content of the genome is 38.18%. The genome contains 133 genes, including 88 protein-coding, 37 tRNA, and 8 rRNA genes. A phylogenetic tree reconstructed by using 16 chloroplast genomes reveals that C. lindleyi is most closely related to C. trichotomum which together forms a group that is a sister to genus Caryopteris. The work reported here is the first complete chloroplast genome of C. lindleyi which will provide useful information to the evolutionary studies on the genus of Clerodendrum.

A complete sequence of mitochondrial genome of Neolamarckia cadamba and its use for systematic analysis

Neolamarckia cadamba is an important tropical and subtropical tree for timber industry in southern China and is also a medicinal plant because of the secondary product cadambine. N. cadamba belongs to Rubiaceae family and its taxonomic relationships with other species are not fully evaluated based on genome sequences. Here, we report the complete sequences of mitochondrial genome of N. cadamba, which is 414,980 bp in length and successfully assembled in two genome circles (109,836 bp and 305,144 bp). The mtDNA harbors 83 genes in total, including 40 protein-coding genes (PCGs), 31 transfer RNA genes, 6 ribosomal RNA genes, and 6 other genes. The base composition of the whole genome is estimated as 27.26% for base A, 22.63% for C, 22.53% for G, and 27.56% for T, with the A + T content of 54.82% (54.45% in the small circle and 54.79% in the large circle). Repetitive sequences account for ~ 0.14% of the whole genome. A maximum likelihood (ML) tree based on DNA sequences of 24 PCGs supports that N. cadamba belongs to order Gentianales. A ML tree based on rps3 gene of 60 species in family Rubiaceae shows that N. cadamba is more related to Cephalanthus accidentalis and Hymenodictyon parvifolium and belongs to the Cinchonoideae subfamily. The result indicates that N. cadamba is genetically distant from the species and genera of Rubiaceae in systematic position. As the first sequence of mitochondrial genome of N. cadamba, it will provide a useful resource to investigate genetic variation and develop molecular markers for genetic breeding in the future.

Plastid phylogenomics of the Gynoxoid group (Senecioneae, Asteraceae) highlights the importance of motif-based sequence alignment amid low genetic distances

PREMISE

The genus Gynoxys and relatives form a species-rich lineage of Andean shrubs and trees with low genetic distances within the sunflower subtribe Tussilaginineae. Previous molecular phylogenetic investigations of the Tussilaginineae have included few, if any, representatives of this Gynoxoid group or reconstructed ambiguous patterns of relationships for it.

METHODS

We sequenced complete plastid genomes of 21 species of the Gynoxoid group and related Tussilaginineae and conducted detailed comparisons of the phylogenetic relationships supported by the gene, intron, and intergenic spacer partitions of these genomes. We also evaluated the impact of manual, motif-based adjustments of automatic DNA sequence alignments on phylogenetic tree inference.

RESULTS

Our results indicate that the inclusion of all plastid genome partitions is needed to infer well-supported phylogenetic trees of the Gynoxoid group. Whole plastome-based tree inference suggests that the genera Gynoxys and Nordenstamia are polyphyletic and form the core clade of the Gynoxoid group. This clade is sister to a clade of Aequatorium and Paragynoxys and also includes some but not all representatives of Paracalia.

CONCLUSIONS

The concatenation and combined analysis of all plastid genome partitions and the construction of manually-curated, motif-based DNA sequence alignments are found to be instrumental in the recovery of well-supported relationships of the Gynoxoid group. We demonstrate that the correct assessment of homology in genome-level plastid sequence data sets is crucial for subsequent phylogeny reconstruction and that the manual post-processing of multiple sequence alignments improves the reliability of such reconstructions amid low genetic distances between taxa.

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

Phylogenetic relationships in the Sorghum genus based on sequencing of the chloroplast and nuclear genes

Sorghum [Sorghum bicolor (L.) Moench] is an important food crop with a diverse gene pool residing in its wild relatives. A total of 15 sorghum accessions from the unexploited wild gene pool of the Sorghum genus, representing the five subgenera, were sequenced, and the complete chloroplast genomes and 99 common single-copy concatenated nuclear genes were assembled. Annotation of the chloroplast genomes identified a total of 81 protein-coding genes, 38 tRNA, and four rRNA genes. The gene content and gene order among the species was identical. A total of 153 nonsynonymous amino acid changes in 40 genes were identified across the species. Phylogenetic analysis of both the whole chloroplast genome and nuclear genes revealed a similar topology with two distinct clades within the genus. The species within the subgenera Eusorghum, Chaetosorghum, and Heterosorghum clustered in one clade, whereas the species within the subgenera Parasorghum and Stiposorghum clustered in a second clade. However, the subgenera Parasorghum and Stiposorghum were not monophyletic, suggesting the need for further research to resolve the relationships within this group. The close relationship between the two monotypic subgenera Chaetosorghum and Heterosorghum suggests that species within these subgenera could be considered as one group. This analysis provides an improved understanding of the genetic relationships within the Sorghum genus and defines diversity in wild sorghum species that may be useful for crop improvement.

Sequencing and phylogenetic analysis of the complete chloroplast genome of Arisaema heterophyllum Blume

Arisaema heterophyllum Blume is a perennial medicinal herb widely distributed in China, Korea and Japan. In this study, the complete chloroplast genome sequence of A. heterophyllum was assembled and characterized based on high-throughput sequencing data. The whole chloroplast genome is 170,610 bp in length and contains 95,485 bp large single-copy (LSC) and 22,605 bp small single-copy (SSC) regions separated by a pair of 26,260 bp inverted repeat (IR) regions. It contained a total of 129 genes, including 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes, with an overall GC content of 34.5%. A phylogenetic tree reconstructed by 30 chloroplast genomes reveals that A. heterophyllum is mostly related to the same genus A. ringens, A. franchetianum and A. erubescens. The complete chloroplast genome of A. heterophyllum was the firstly reported and deposited at GenBank under accession number MZ424448.

The complete chloroplast genome of Agave angustifolia

Agave angustifolia is commonly used for the production of bacanora, a kind of fermented and distilled beverage in Mexico. In the present study, we have successfully assembled its chloroplast genome. The full length of the genome is 157,274 bp with the GC content of 37.84%. There is a large single copy region (LSC) of 85,895 bp, a pair of inverted repeat regions (IR) of 26,575 bp and a small single copy region (SSC) of 18,229 bp in the genome. A total of 132 genes are annotated in the cp genome. Among these, there are 86 protein-coding genes, 38 tRNAs and 8 rRNAs. Phylogenetic analysis reveals that A. angustifolia is closely related with A. H11648.

Whole genome based insights into the phylogeny and evolution of the Juglandaceae

BACKGROUND

The walnut family (Juglandaceae) contains commercially important woody trees commonly called walnut, wingnut, pecan and hickory. Phylogenetic relationships and diversification within the Juglandaceae are classic and hot scientific topics that have been elucidated by recent fossil, morphological, molecular, and (paleo) environmental data. Further resolution of relationships among and within genera is still needed and can be achieved by analysis of the variation of chloroplast, mtDNA, and nuclear genomes.

RESULTS

We reconstructed the backbone phylogenetic relationships of Juglandaceae using organelle and nuclear genome data from 27 species. The divergence time of Juglandaceae was estimated to be 78.7 Mya. The major lineages diversified in warm and dry habitats during the mid-Paleocene and early Eocene. The plastid, mitochondrial, and nuclear phylogenetic analyses all revealed three subfamilies, i.e., Juglandoideae, Engelhardioideae, Rhoipteleoideae. Five genera of Juglandoideae were strongly supported. Juglandaceae were estimated to have originated during the late Cretaceous, while Juglandoideae were estimated to have originated during the Paleocene, with evidence for rapid diversification events during several glacial and geological periods. The phylogenetic analyses of organelle sequences and nuclear genome yielded highly supported incongruence positions for J. cinerea, J. hopeiensis, and Platycarya strobilacea. Winged fruit were the ancestral condition in the Juglandoideae, but adaptation to novel dispersal and regeneration regimes after the Cretaceous-Paleogene boundary led to the independent evolution of zoochory among several genera of the Juglandaceae.

CONCLUSIONS

A fully resolved, strongly supported, time-calibrated phylogenetic tree of Juglandaceae can provide an important framework for studying classification, diversification, biogeography, and comparative genomics of plant lineages. Our addition of new, annotated whole chloroplast genomic sequences and identification of their variability informs the study of their evolution in walnuts (Juglandaceae).

Whole genome based insights into the phylogeny and evolution of the Juglandaceae

BACKGROUND

The walnut family (Juglandaceae) contains commercially important woody trees commonly called walnut, wingnut, pecan and hickory. Phylogenetic relationships and diversification within the Juglandaceae are classic and hot scientific topics that have been elucidated by recent fossil, morphological, molecular, and (paleo) environmental data. Further resolution of relationships among and within genera is still needed and can be achieved by analysis of the variation of chloroplast, mtDNA, and nuclear genomes.

RESULTS

We reconstructed the backbone phylogenetic relationships of Juglandaceae using organelle and nuclear genome data from 27 species. The divergence time of Juglandaceae was estimated to be 78.7 Mya. The major lineages diversified in warm and dry habitats during the mid-Paleocene and early Eocene. The plastid, mitochondrial, and nuclear phylogenetic analyses all revealed three subfamilies, i.e., Juglandoideae, Engelhardioideae, Rhoipteleoideae. Five genera of Juglandoideae were strongly supported. Juglandaceae were estimated to have originated during the late Cretaceous, while Juglandoideae were estimated to have originated during the Paleocene, with evidence for rapid diversification events during several glacial and geological periods. The phylogenetic analyses of organelle sequences and nuclear genome yielded highly supported incongruence positions for J. cinerea, J. hopeiensis, and Platycarya strobilacea. Winged fruit were the ancestral condition in the Juglandoideae, but adaptation to novel dispersal and regeneration regimes after the Cretaceous-Paleogene boundary led to the independent evolution of zoochory among several genera of the Juglandaceae.

CONCLUSIONS

A fully resolved, strongly supported, time-calibrated phylogenetic tree of Juglandaceae can provide an important framework for studying classification, diversification, biogeography, and comparative genomics of plant lineages. Our addition of new, annotated whole chloroplast genomic sequences and identification of their variability informs the study of their evolution in walnuts (Juglandaceae).

Complete chloroplast genome of Chondria tumulosa (Ceramiales, Rhodophyta), a recently described cryptogenic species with invasive traits from Papahānaumokuākea Marine National Monument, Hawai'i

The complete chloroplast genome of Chondria tumulosa, a red alga from Manawai (Pearl and Hermes Atoll), Hawai‘i, was determined and analyzed using next-generation sequencing and de novo assembly approaches. The chloroplast genome sequence of C. tumulosa was 172,617 bp and contained 231 genes, consisting of 197 protein-coding genes, 29 transfer RNA genes, three ribosomal RNA genes, one transfer-messenger RNA gene, one non-coding RNA gene, and one intron inserted into the trnM gene. The number of genes and genome structure was largely similar to other members of the family Rhodomelaceae. The phylogenomic analysis of 32 complete cpDNA from the red algal order Ceramiales showed that C. tumulosa is a distinct species within the Chondrieae tribe, and is a diverging early relative to the other three available Chondria chloroplast genomes.

Complete chloroplast genome sequence of Mahonia duclouxiana (Berberidaceae), a medicinal plant in China

Mahonia duclouxiana is a member of the genus Mahonia of Berberidaceae and is distributed in South Asia. Here, the complete chloroplast genome sequence of M. duclouxiana was reported. The complete chloroplast genome is 165,384 bp in length, which has a large single-copy (LSC) region of 73,477 bp, a small single-copy (SSC) region of 18,563 bp, and two inverted repeat (IR) regions of 36,672 bp. The G/C content in the chloroplast genome is 38.1%. The whole chloroplast genome contains 151 genes, including 38 unique tRNA genes, 105 unique protein-coding genes, and 8 unique rRNA genes. The phylogenetic analysis supported that this species should be included in Maddenia. The complete chloroplast genome sequence of M. duclouxiana will provide extremely important information in tracing the evolutionary history of the genus Mahonia and the development of the medicinal value.

The complete chloroplast genome sequence of Arundo formosana Hack. (Poaceae)

Arundo formosana Hack. belongs to the Arundionideae subfamily of Poaceae. In this study, we sequenced and assembled the complete chloroplast genome of A. formosana. The complete chloroplast genome was 136,919 bp in size, including a large single copy region of 82,039 bp, a small single-copy region of 12,108 bp and a pair of reverse repeats of 21,386 bp in size. The annotation of A. formosana indicates that it contained 81 protein-coding genes, 47 tRNA and 8 rRNA. Our phylogenetic analysis of the 36 grass complete chroloplast genomes of protein-coding genes using Cyperus rotundus as outgroup showed that A. formosana is closely related to Crinipes species to form the Arundionideae lineage of the grass family.

The complete chloroplast genome of Prinsepia uniflora (Rosaceae), a medicinal plant found in northwest China

Prinsepia uniflora Batalin 1892 is a medicinal plant widely distributed in northwest China. In this study, we report and characterize the complete chloroplast (cp) genome sequence of P. uniflora. The entire sequence is 159,179 bp in length, consisting of the large single-copy region (LSC) and small single copy region (SSC) (87,239 and 19,180 bp, respectively); these two regions are separated by a pair of 26,380-bp inverted repeat (IR) regions. The genome contains 131 genes, including 86 protein-coding genes, 37 tRNA genes, and eight rRNA genes. The overall GC content of the genome is 36.7%. A phylogenetic tree constructed from 18 chloroplast genomes revealed that P. uniflora was clustered with Prinsepia sinensis and Prinsepia utilis, all of which belong to the genus Prinsepia, which is supported as a sister group by a moderate bootstrap support value of 55% with the Malus and Pyrus genera.

Mitochondrial genome of Theileria uilenbergi endemic in sheep and goats in China

Mitochondrial genomes provide new insights that help elucidating biological features, genetic evolution, and classification of protozoans. Theileria uilenbergi (T. uilenbergi), transmitted by Haemaphysalis qinghaiensis and H. longicornis, is considered as highly pathogenic to sheep and goats in China. This study reports and outlines features of its mitochondrial genome. The T. uilenbergi mitochondrial genome is a linear monomeric molecule of 6.0 kb length, which encodes three protein-coding genes named cytochrome c oxidase I (cox1), cytochrome b (cob), and cytochrome c oxidase III (cox3), as well as six large subunit (LSU) rRNA gene fragments, and ends in terminal inverted repeats (TIRs). The array structure and organization of the mitochondrial genome of T. uilenbergi is identical to that of T. parva. Phylogenetic analysis based on the amino acid sequences of cox1, cob, and cox3 genes suggests that T. uilenbergi is distantly related to the group of transforming Theileria species such as T. parva. This study contributes to a comprehensive understanding of the phylogeny and evolution of the mitochondrial genome of piroplasms and provides useful information of diagnostic marker for T. uilenbergi.

The first complete chloroplast genome sequences in Resedaceae: Genome structure and comparative analysis

Caylusea hexagyna and Ochradenus baccatus are two species in the Resedaceae family. In this study, we analysed the complete plastid genomes of these two species using high-throughput sequencing technology and compared their genomic data. The length of the plastid genome of C. hexagyna was 154,390 bp while that of O. baccatus was 153,380 bp. The lengths of the inverted repeats (IR) regions were 26,526 bp and 26,558 bp, those of the large single copy (LSC) regions were 83,870 bp and 83,023 bp; and those of the small single copy (SSC) regions were 17,468 bp and 17,241 bp in C. hexagyna and O. baccatus, respectively. Both genomes consisted of 113 genes: 79 protein-coding genes, 30 tRNA genes and 4 rRNA genes. Repeat analysis showed that the plastid genome included all types of repeats, with more frequent occurrences of palindromic sequences. Comparative studies of SSR markers showed that there were 256 markers in C. hexagyna and 255 in O. baccatus; the majority of the SSRs in these plastid genomes were mononucleotide repeats (A/T). All the clusters in the phylogenetic tree had high support. This study reported the first complete plastid genomes of the genera Caylusea and Ochradenus and the first for the Resedaceae family.

The complete chloroplast genome sequence of Asplenium komarovii Akasawa, a rare fern in South Korea

We sequenced the complete chloroplast genome of Asplenium komarovii Akasawa (syn: Asplenium scolopendrium L. subsp. japonicum (Komarov) Rasbach, Reichstein & Viane), which is designated as a rare species in South Korea. The complete chloroplast genome is 149,393 bp in total length and comprised of the following regions: large single copy (82,464 bp), small single copy (21,345 bp), and a pair of inverted repeats (22,792 bp). The overall GC content is 40.9% and the genome encoded a total of 115 genes, including 84 protein-coding, 27 transfer RNA, and 4 ribosomal RNA genes. Phylogenetic analysis based on 21 representative chloroplast genomes of the suborder Aspleniineae (and one outgroup) indicates that Aspleniaceae is monophyletic and sister to Diplaziopsidaceae, with Rhadidosoraceae as the basal group in this three family clade. Asplenium komarovii is sister to A. nidus and A. prolongatum with strong bootstrap support. The chloroplast genome of A. komarovii will be useful in establishing its relationships within the A. scolopendrium complex, which is currently unresolved.

Comparative Chloroplast Genomes of Zosteraceae Species Provide Adaptive Evolution Insights Into Seagrass

Seagrasses are marine flowering plants found in tropical and sub-tropical areas that live in coastal regions between the sea and land. All seagrass species evolved from terrestrial monocotyledons, providing the opportunity to study plant adaptation to sea environments. Here, we sequenced the chloroplast genomes (cpGenomes) of three Zostera species, then analyzed and compared their cpGenome structures and sequence variations. We also performed a phylogenetic analysis using published seagrass chloroplasts and calculated the selection pressure of 17 species within seagrasses and nine terrestrial monocotyledons, as well as estimated the number of shared genes of eight seagrasses. The cpGenomes of Zosteraceae species ranged in size from 143,877 bp (Zostera marina) to 152,726 bp (Phyllospadix iwatensis), which were conserved and displayed similar structures and gene orders. Additionally, we found 17 variable hotspot regions as candidate DNA barcodes for Zosteraceae species, which will be helpful for studying the phylogenetic relationships and interspecies differences between seagrass species. Interestingly, nine genes had positive selection sites, including two ATP subunit genes (atpA and atpF), two ribosome subunit genes (rps4 and rpl20), two DNA-dependent RNA polymerase genes (rpoC1 and rpoC2), as well as accD, clpP, and ycf2. These gene regions may have played key roles in the seagrass adaptation to diverse environments. The Branch model analysis showed that seagrasses had a higher rate of evolution than terrestrial monocotyledons, suggesting that seagrasses experienced greater environmental pressure. Moreover, a branch-site model identified positively selected sites (PSSs) in ccsA, suggesting their involvement in the adaptation to sea environments. These findings are valuable for further investigations on Zosteraceae cpGenomes and will serve as an excellent resource for future studies on seagrass adaptation to sea environments.

The complete chloroplast genome and characteristics analysis of Musa basjoo Siebold

BACKGROUND

An ornamental plant often seen in gardens and farmhouses, Musa basjoo Siebold can also be used as Chinese herbal medicine. Its pseudostem and leaves are diuretic; its root can be decocted together with ginger and licorice to cure gonorrhea and diabetes; the decoct soup of its pseudostem can help relieve heat, and the decoct soup of its dried flower can treat cerebral hemorrhage. There have not been many chloroplast genome studies on M. basjoo Siebold.

METHODS AND RESULTS

We characterized its complete chloroplast genome using Novaseq 6000 sequencing. This paper shows that the length of the chloroplast genome M. basjoo Siebold is 172,322 bp, with 36.45% GC content. M. basjoo Siebold includes a large single-copy region of 90,160 bp, a small single-copy region of 11,668 bp, and a pair of inverted repeats of 35,247 bp. Comparing the genomic structure and sequence data of closely related species, we have revealed the conserved gene order of the IR and LSC/SSC regions, which has provided a very inspiring discovery for future phylogenetic research.

CONCLUSIONS

Overall, this study has constructed an evolutionary tree of the genus Musa species with the complete chloroplast genome sequence for the first time. As can be seen, there is no obvious multi-branching in the genus, and M. basjoo Siebold and Musa itinerans are the closest relatives.

tRNAscan-SE 2.0: improved detection and functional classification of transfer RNA genes

tRNAscan-SE has been widely used for transfer RNA (tRNA) gene prediction for over twenty years, developed just as the first genomes were decoded. With the massive increase in quantity and phylogenetic diversity of genomes, the accurate detection and functional prediction of tRNAs has become more challenging. Utilizing a vastly larger training set, we created nearly one hundred specialized isotype- and clade-specific models, greatly improving tRNAscan-SE’s ability to identify and classify both typical and atypical tRNAs. We employ a new comparative multi-model strategy where predicted tRNAs are scored against a full set of isotype-specific covariance models, allowing functional prediction based on both the anticodon and the highest-scoring isotype model. Comparative model scoring has also enhanced the program's ability to detect tRNA-derived SINEs and other likely pseudogenes. For the first time, tRNAscan-SE also includes fast and highly accurate detection of mitochondrial tRNAs using newly developed models. Overall, tRNA detection sensitivity and specificity is improved for all isotypes, particularly those utilizing specialized models for selenocysteine and the three subtypes of tRNA genes encoding a CAU anticodon. These enhancements will provide researchers with more accurate and detailed tRNA annotation for a wider variety of tRNAs, and may direct attention to tRNAs with novel traits.

The complete chloroplast genome sequence of Bromus catharticus Vahl. (Poaceae)

Bromus catharticus Vahl. belongs to the Pooideae subfamily of Poaceae. In this study, we sequenced and assembled the complete chloroplast genome of B. catharticus. The complete chloroplast genome was 134,718 bp in size, including a large single-copy region of 80,540 bp, a small single-copy region of 11,806 bp and a pair of reverse repeats of 21,186 bp in size. The annotation of B. catharticus indicates that it contained 89 protein-coding genes, 47 tRNA genes and eight rRNA genes. Our phylogenetic analysis of all protein-coding genes of the 36 grass complete chroloplast genomes using Cyperus rotundus as outgroup showed that B. catharticus is closely related to the Koeleria and Avena species to form the Pooideae lineage of the grass family.

Characterization of the chloroplast genome sequence of Bonia amplexicaulis (L.C.Chia, H.L.Fung & Y.L.Yang) N.H.Xia (Poaceae)

Bonia amplexicaulis (L.C.Chia, H.L.Fung & Y.L.Yang) N.H.Xia is a member of the Bambusoideae subfamily in Poaceae. In this study, we sequenced, assembled and characterized the complete chloroplast genome of B. amplexicaulis. The complete chloroplast genome was 139,935 bp in size, including a large single copy region of 83,453 bp, a small single-copy region of 12,860 bp and a pair of reverse repeats of 21,811 bp in size. The annotation of the B. amplexicaulis chloroplast genome indicates that it contained 83 protein-coding genes, 36 tRNA genes and 8 rRNA genes. Our phylogenetic analysis of all protein-coding genes from the 36 complete chroloplast grass genomes using Cyperus rotundus as outgroup showed that B. amplexicaulis is closely related to Otatea glauca and Pariana campestris to form the Bambusoideae lineage of the grass family.

The chloroplast genome of Amygdalus L. (Rosaceae) reveals the phylogenetic relationship and divergence time

Background

Limited access to genetic information has greatly hindered our understanding of the molecular evolution, phylogeny, and differentiation time of subg. Amygdalus. This study reported complete chloroplast (cp) genome sequences of subg. Amygdalus, which further enriched the available valuable resources of complete cp genomes of higher plants and deepened our understanding of the divergence time and phylogenetic relationships of subg. Amygdalus.

Results

The results showed that subg. Amygdalus species exhibited a tetrad structure with sizes ranging from 157,736 bp (P. kansuensis) to 158,971 bp (P. davidiana), a pair of inverted repeat regions (IRa/IRb) that ranged from 26,137–26,467 bp, a large single-copy region that ranged from 85,757–86,608 bp, and a small single-copy region that ranged from 19,020–19,133 bp. The average GC content of the complete cp genomes in the 12 species was 36.80%. We found that the structure of the subg. Amygdalus complete cp genomes was highly conserved, and the 12 subg. Amygdalus species had an rps19 pseudogene. There was not rearrangement of the complete cp genome in the 12 subg. Amygdalus species. All 12 subg. Amygdalus species clustered into one clade based on both Bayesian inference and maximum likelihood. The divergence time analyses based on the complete cp genome sequences showed that subg. Amygdalus species diverged approximately 15.65 Mya.

Conclusion

Our results provide data on the genomic structure of subg. Amygdalus and elucidates their phylogenetic relationships and divergence time.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07968-6.

Proposal to recognise the tribes Adinobotryeae and Glycyrrhizeae (Leguminosae subfamily Papilionoideae) based on chloroplast phylogenomic evidence

Within the legume family, the taxonomic status of subtribe Glycyrrhizinae of tribe Galegeae and of the genus Adinobotrys has been re-assessed. Based on genome skimming data, we conducted phylogenomic analyses of the inverted repeat-lacking clade within subfamily Papilionoideae. The results support the sister relationship between Glycyrrhizeae and Adinobotrys. Glycyrrhizeae is resurrected based on Glycyrrhiza and Glycyrrhizopsis, and a new tribe, Adinobotryeae, is proposed to accommodate Adinobotrys.

The complete chloroplast genome of Keteleeria evelyniana

Here, we report the complete chloroplast genome of Keteleeria evelyniana. The genome is 116,940 bp in size, which is comprised of a large single-copy (LSC) region of 74,075 bp, a small single-copy (SSC) region of 40,425 bp, and two short inverted repeat (IR) regions of 1,220 bp. The overall GC content of the plastome was 38.5%. The new sequence comprised 103 unique genes, including 74 protein-coding genes, 4 rRNA genes, and 25 tRNA genes. Phylogenetic analysis showed that K.evelyniana was close to Keteleeria hainanensis and Keteleeria davidiana.

Characterization and Comparative Analysis of Complete Chloroplast Genomes of Three Species From the Genus Astragalus (Leguminosae)

Astragalus is the largest genus in Leguminosae. Several molecular studies have investigated the potential adulterants of the species within this genus; nonetheless, the evolutionary relationships among these species remain unclear. Herein, we sequenced and annotated the complete chloroplast genomes of three Astragalus species—Astragalus adsurgens, Astragalus mongholicus var. dahuricus, and Astragalus melilotoides using next-generation sequencing technology and plastid genome annotator (PGA) tool. All species belonged to the inverted repeat lacking clade (IRLC) and had similar sequences concerning gene contents and characteristics. Abundant simple sequence repeat (SSR) loci were detected, with single-nucleotide repeats accounting for the highest proportion of SSRs, most of which were A/T homopolymers. Using Astragalus membranaceus var. membranaceus as reference, the divergence was evident in most non-coding regions of the complete chloroplast genomes of these species. Seven genes (atpB, psbD, rpoB, rpoC1, trnV, rrn16, and rrn23) showed high nucleotide variability (Pi), and could be used as DNA barcodes for Astragalus sp. cemA and rpl33 were found undergoing positive selection by the section patterns in the coded protein. Phylogenetic analysis showed that Astragalus is a monophyletic group closely related to the genus Oxytropis within the tribe Galegeae. The newly sequenced chloroplast genomes provide insight into the unresolved evolutionary relationships within Astragalus spp. and are expected to contribute to species identification.

Characterization of the DNA molecular sequence of complete plastid genome of Paeonia rockii subsp. taibaishanica, an endemic species in China

Paeonia rockii subsp. taibaishanica (Paeoniaceae), one of the tree peony species, is endemic to the Qinling Mountains in central China. In this study, we characterized its whole plastid genome sequence using the Illumina sequencing platform. The complete plastid genome size of P. rockii subsp. taibaishanica is 153,368 bp in length, including a large single copy (LSC) region of 85,030 bp, a small single copy (SSC) region of 17,042 bp, and a pair of inverted repeats (IRs) of 25,648 bp. The genome contains 131 genes, including 83 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The GC contents in chloroplast genome, LSC region, SSC region, and IR region were 38.3%, 36.6%, 32.6%, and 43.1%, respectively. A total of 16 species are used to construct the phylogenetic tree of Paeoniaceae, the results showed that P. rockii subsp. taibaishanica is more closely related with congeneric Paeonia suffruticosa and Paeonia ostii, these species were clustered into a clade with high bootstrap support.

The complete chloroplast genome of Rhodobryum laxelimbatum (Hampe ex Ochi) Z. Iwatsuki and T. J. Koponen

Rhodobryum laxelimbatum (Bryophyta, Bryaceae) is one of the folk medicine resources in Southwest China, which has excellent potential for application in treating cardiovascular diseases. In this study, R. laxelimbatum was sequenced by high-throughput sequencing technology. The complete chloroplast genome is 124,632 bp in length with a quadripartite structure. Two inverted repeat regions are 9837 bp, separated by a large single copy region of 86,444 bp and a small single copy region of 18,514 bp. It encodes 118 unique genes, including 82 protein-coding genes, 32 tRNA genes, and four rRNA genes. The phylogenetic tree was constructed based on the complete chloroplast genome sequences of 18 bryophytes, downloaded from GenBank and acquired in this study. The phylogenetic analysis strongly indicated that R. laxelimbatum was the sister group of the clade which consists of Mnium marginatum, Pohlia cruda and Pohlia nutans. The R. laxelimbatum chloroplast genome sequence provides new genomic resources, which will improve its research, conservation, and application in the future.

The complete chloroplast genome sequence of Oxystophyllum changjiangense (Orchidaceae)

Oxystophyllum changjiangense has high economic value due to its wide applications in horticultural and medicinal fields. Here, the first chloroplast genome of O. changjiangense was sequenced and reported. The chloroplast genome displayed the typical quadripartite structure containing a pair of inverted repeats (IR), a long single-copy region (LSC), and a short single-copy region (SSC). Total 110 genes were found including 76 protein-coding genes, 30 tRNA genes, and four rRNA genes. Phylogeny analysis showed O. changjiangense has a close relationship with Phaius species.

The complete chloroplast genome of Rhododendron molle and its phylogenetic position within Ericaceae

Rhododendron mole (Blume) G. Don is an attractive ornamental and valuable medicinal plant which widely distributed in the southern regions of China. In order to promote the studies on the genetic diversity of this species, we assembled the complete chloroplast (cp) genome of R. molle by using the genome skimming approach. The results showed that the cp genome of R. molle exhibited a quadripartite cycle with 197,877 bp, comprising of two inverted repeats (IRs) of 43,831 bp separated by a large single copy (LSC) region of 110,189 bp and a quite small single copy (SSC) region of 26 bp. It encodes 146 genes, including 92 protein-coding, 46 tRNA, and eight rRNA genes. The overall GC content of the cp genome was 36.0%. The phylogenetic analysis indicated that R. molle is closely related to R. delavayi. Thus, the cp genome sequence of R. molle provides a rich source of genetic information for studies on Rhododendron taxonomy, phylogeny, and evolution, as well as lays the foundation for further development and utilization of R. molle.

Characterization of the complete chloroplast genome of an annual herb, Chenopodium album (Amaranthaceae)

Chenopodium album is an annual herb from Amaranthaceae with worldwide distribution. It is a leafy vegetable as well as an important subsidiary grain crop with high nutritional value and medicinal value. In this study, we reported the complete chloroplast genome of C. album. The total chloroplast genome was 152,167 bp in length, containing a large single-copy region (LSC, 83,676 bp), a small single-copy region (SSC, 18,105 bp), and a pair of inverted repeat regions (IRs, 25,193 bp). The complete chloroplast genome contains 110 genes, including 78 protein-coding genes, 28 transfer RNA (tRNA) genes, and 4 ribosomal RNA (rRNA) genes with an overall GC content of 37.3%. Phylogenetic analysis showed that C. album was sister to C. acuminatum within Chenopodioideae. The complete chloroplast genome of C. album will provide useful resources for the development and utilization of this species and the phylogenetic study of Amaranthaceae.

Comparative analysis of plastomes in Oxalidaceae: Phylogenetic relationships and potential molecular markers

The wood sorrel family, Oxalidaceae, is mainly composed of annual or perennial herbs, a few shrubs, and trees distributed from temperate to tropical zones. Members of Oxalidaceae are of high medicinal, ornamental, and economic value. Despite the rich diversity and value of Oxalidaceae, few molecular markers or plastomes are available for phylogenetic analysis of the family. Here, we reported four new whole plastomes of Oxalidaceae and compared them with plastomes of three species in the family, as well as the plastome of Rourea microphylla in the closely related family Connaraceae. The eight plastomes ranged in length from 150,673 bp (Biophytum sensitivum) to 156,609 bp (R. microphylla). Genome annotations revealed a total of 129-131 genes, including 83-84 protein-coding genes, eight rRNA genes, 37 tRNA genes, and two to three pseudogenes. Comparative analyses showed that the plastomes of these species have minor variations at the gene level. The smaller plastomes of herbs B. sensitivum and three Oxalis species are associated with variations in IR region sizes, intergenic region variation, and gene or intron loss. We identified sequences with high variation that may serve as molecular markers in taxonomic studies of Oxalidaceae. The phylogenetic trees of selected superrosid representatives based on 76 protein-coding genes corroborated the Oxalidaceae position in Oxalidales and supported it as a sister to Connaraceae. Our research also supported the monophyly of the COM (Celastrales, Oxalidales, and Malpighiales) clade.