Characterization of the Chloroplast Genome Sequence of Acer miaotaiense: Comparative and Phylogenetic Analyses

Comparative complete chloroplast genome of Geum japonicum: evolution and phylogenetic analysis

Geum japonicum (Rosaceae) has been widely used in China as a traditional herbal medicine due to its high economic and medicinal value. However, the appearance of Geum species is relatively similar, making identification difficult by conventional phenotypic methods, and the studies of genomics and species evolution are lacking. To better distinguish the medicinal varieties and fill this gap, we carried out relevant research on the chloroplast genome of G. japonicum. Results show a typical quadripartite structure of the chloroplast genome of G. japonicum with a length of 156,042 bp. There are totally 131 unique genes in the genome, including 87 protein-coding genes, 36 tRNA genes, and 8 rRNA genes, and there were also 87 SSRs identified and mostly mononucleotide Adenine-Thymine. We next compared the plastid genomes among four Geum species and obtained 14 hypervariable regions, including ndhF, psbE, trnG-UCC, ccsA, trnQ-UUG, rps16, psbK, trnL-UAA, ycf1, ndhD, atpA, petN, rps14, and trnK-UUU. Phylogenetic analysis revealed that G. japonicum is most closely related to Geum aleppicum, and possibly has some evolutionary relatedness with an ancient relic plant Taihangia rupestris. This research enriched the genome resources and provided fundamental insights for evolutionary studies and the phylogeny of Geum.

Comparative analysis of the complete chloroplast genomes of six threatened subgenus Gynopodium (Magnolia) species

BACKGROUND

The subgenus Gynopodium belonging to genus Magnolia have high ornamental, economic, and ecological value. Subgenus Gynopodium contains eight species, but six of these species are threatened. No studies to date have characterized the characteristics of the chloroplast genomes (CPGs) within subgenus Gynopodium species. In this study, we compared the structure of CPGs, identified the mutational hotspots and resolved the phylogenetic relationship of subgenus Gynopodium.

RESULTS

The CPGs of six subgenus Gynopodium species ranged in size from 160,027 bp to 160,114 bp. A total of 131 genes were identified, including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. We detected neither major expansions or contractions in the inverted repeat region, nor rearrangements or insertions in the CPGs of six subgenus Gynopodium species. A total of 300 large repeat sequences (forward, reverse, and palindrome repeats), 847 simple sequence repeats, and five highly variable regions were identified. One gene (ycf1) and four intergenic regions (psbA-trnH-GUG, petA-psbJ, rpl32-trnL-UAG, and ccsA-ndhD) were identified as mutational hotspots by their high nucleotide diversity (Pi) values (≥ 0.004), which were useful for species discrimination. Maximum likelihood and Bayesian inference trees were concordant and indicated that Magnoliaceae consisted of two genera Liriodendron and Magnolia. Six species of subgenus Gynopodium clustered as a monophyletic clade, forming a sister clade with subgenus Yulania (BS = 100%, PP = 1.00). Due to the non-monophyly of subgenus Magnolia, subgenus Gynopodium should be treated as a section of Magnolia. Within section Gynopodium, M. sinica diverged first (posterior probability = 1, bootstrap = 100), followed by M. nitida, M. kachirachirai and M. lotungensis. M. omeiensis was sister to M. yunnanensis (posterior probability = 0.97, bootstrap = 50).

CONCLUSION

The CPGs and characteristics information provided by our study could be useful in species identification, conservation genetics and resolving phylogenetic relationships of Magnoliaceae species.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

Background

Calanthe (Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of the Calanthe alliance (Calanthe, Cephalantheropsis, and Phaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed a de novo assembly, and did a comparative analysis of 8 Calanthe group species' plastomes: 6 Calanthe and 2 Phaius species. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

Results

The complete plastomes of the Calanthe group species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost in C. delavayi (P. delavayi), including ndhC, ndhF, and ndhK genes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in the Calanthe group, and confirmed the position of Phaius delavayi in the genus Calanthe as opposed to its previous placement in Phaius.

Conclusion

This study provides a report on the complete plastomes of 6 Calanthe and 2 Phaius species and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

BACKGROUND

Calanthe (Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of the Calanthe alliance (Calanthe, Cephalantheropsis, and Phaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed a de novo assembly, and did a comparative analysis of 8 Calanthe group species' plastomes: 6 Calanthe and 2 Phaius species. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

RESULTS

The complete plastomes of the Calanthe group species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost in C. delavayi (P. delavayi), including ndhC, ndhF, and ndhK genes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in the Calanthe group, and confirmed the position of Phaius delavayi in the genus Calanthe as opposed to its previous placement in Phaius.

CONCLUSION

This study provides a report on the complete plastomes of 6 Calanthe and 2 Phaius species and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

Structural Characterization of the Acer ukurunduense Chloroplast Genome Relative to Related Species in the Acer Genus

Acer ukurunduense refers to a deciduous tree distributed in Northeast Asia and is a widely used landscaping tree species. Although several studies have been conducted on the species’ ecological and economic significance, limited information is available on its phylo-genomics. Our study newly constitutes the complete chloroplast genome of A. ukurunduense into a 156,645-bp circular DNA, which displayed a typical quadripartite structure. In addition, 133 genes were identified, containing 88 protein-coding genes, 37 tRNA genes, and eight rRNA genes. In total, 107 simple sequence repeats and 49 repetitive sequences were observed. Thirty-two codons indicated that biased usages were estimated across 20 protein-coding genes (CDS) in A. ukurunduense. Four hotspot regions (trnK-UUU/rps16, ndhF/rpl32, rpl32/trnL-UAG, and ycf1) were detected among the five analyzed Acer species. Those hotspot regions may be useful molecular markers and contribute to future population genetics studies. The phylogenetic analysis demonstrated that A. ukurunduense is most closely associated with the species of Sect. Palmata. A. ukurunduense and A. pubipetiolatum var. pingpienense diverged in 22.11 Mya. We selected one of the hypervariable regions (trnK-UUU/rps16) to develop a new molecular marker and designed primers and confirmed that the molecular markers could accurately discriminate five Acer species through Sanger sequencing. By sequencing the cp genome of A. ukurunduense and comparing it with the relative species of Acer, we can effectively address the phylogenetic problems of Acer at the species level and provide insights into future research on population genetics and genetic diversity.

Comparative Genomic and Phylogenetic Analysis of Chloroplast Genomes of Hawthorn (Crataegus spp.) in Southwest China

The hawthorns (Crataegus spp.) are widely distributed and famous for their edible and medicinal values. There are ∼18 species and seven varieties of hawthorn in China distributed throughout the country. We now report the chloroplast genome sequences from C. scabrifolia, C. chungtienensis and C. oresbia, from the southwest of China and compare them with the previously released six species in Crataegus and four species in Rosaceae. The chloroplast genome structure of Crataegus is typical and can be divided into four parts. The genome sizes are between 159,654 and 159,898bp. The three newly sequenced chloroplast genomes encode 132 genes, including 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Comparative analysis of the chloroplast genomes revealed six divergent hotspot regions, including ndhA, rps16-trnQ-UUG, ndhF-rpl32, rps16-psbK, trnR-UCU-atpA and rpl32-trnL-UAG. According to the correlation and co-occurrence analysis of repeats with indels and SNPs, the relationship between them cannot be ignored. The phylogenetic tree constructed based on the complete chloroplast genome and intergenic region sequences indicated that C. scabrifolia has a different origin from C. chungtienensis and C. oresbia. We support the placement of C. hupehensis, C. cuneata, C. scabrifolia in C. subg. Crataegus and C. kansuensis, C. oresbia, C. kansuensis in C. subg. Sanguineae. In addition, based on the morphology, geographic distribution and phylogenetic relationships of C. chungtienensis and C. oresbia, we speculate that these two species may be the same species. In conclusion, this study has enriched the chloroplast genome resources of Crataegus and provided valuable information for the phylogeny and species identification of this genus.

The Complete Chloroplast Genome of Carya cathayensis and Phylogenetic Analysis

Carya cathayensis, an important economic nut tree, is narrowly endemic to eastern China in the wild. The complete cp genome of C. cathayensis was sequenced with NGS using an Illumina HiSeq2500, analyzed, and compared to its closely related species. The cp genome is 160,825 bp in length with an overall GC content of 36.13%, presenting a quadripartite structure comprising a large single copy (LSC; 90,115 bp), a small single copy (SSC; 18,760 bp), and a pair of inverted repeats (IRs; 25,975 bp). The genome contains 129 genes, including 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. A total of 252 simple sequence repeats (SSRs) and 55 long repeats were identified. Gene selective pressure analysis showed that seven genes (rps15, rpoA, rpoB, petD, ccsA, atpI, and ycf1-2) were possibly under positive selection compared with the other Juglandaceae species. Phylogenetic relationships of 46 species inferred that Juglandaceae is monophyletic, and that C. cathayensis is sister to Carya kweichowensis and Carya illinoinensis. The genome comparison revealed that there is a wide variability of the junction sites, and there is higher divergence in the noncoding regions than in coding regions. These results suggest a great potential in phylogenetic research. The newly characterized cp genome of C. cathayensis provides valuable information for further studies of this economically important species.

Comparative Analysis of the Chloroplast Genome for Four Pennisetum Species: Molecular Structure and Phylogenetic Relationships

The genus Pennisetum (Poaceae) is both a forage crop and staple food crop in the tropics. In this study, we obtained chloroplast genome sequences of four species of Pennisetum (P. alopecuroides, P. clandestinum, P. glaucum, and P. polystachion) using Illumina sequencing. These chloroplast genomes have circular structures of 136,346–138,119 bp, including a large single-copy region (LSC, 79,380–81,186 bp), a small single-copy region (SSC, 12,212–12,409 bp), and a pair of inverted repeat regions (IRs, 22,284–22,372 bp). The overall GC content of these chloroplast genomes was 38.6–38.7%. The complete chloroplast genomes contained 110 different genes, including 76 protein-coding genes, 30 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes. Comparative analysis of nucleotide variability identified nine intergenic spacer regions (psbA-matK, matK-rps16, trnN-trnT, trnY-trnD-psbM, petN-trnC, rbcL-psaI, petA-psbJ, psbE-petL, and rpl32-trnL), which may be used as potential DNA barcodes in future species identification and evolutionary analysis of Pennisetum. The phylogenetic analysis revealed a close relationship between P. polystachion and P. glaucum, followed by P. clandestinum and P. alopecuroides. The completed genomes of this study will help facilitate future research on the phylogenetic relationships and evolution of Pennisetum species.

The complete chloroplast genome and characteristics analysis of Callistemon rigidus R.Br

Callistemon rigidus R.Br. one of the traditional Chinese medicinal plants, is acrid-flavored and mild-natured, with the prominent effects reducing swelling, resolving phlegm, and dispelling rheumatism. Clinically, it has been commonly used to treat cold, cough and asthma, pain and swelling from impact injuries, eczema, rheumatic arthralgia. The chloroplast genome study on Callistemon rigidus R.Br. is a few seen. This study demonstrates the data collected from the assembly and annotation of the chloroplast (cp) genome of Callistemon rigidus R.Br., followed by furthers comparative analysis with the cp genomes of closely related species. C. rigidus R.Br. showed a cp genome in the size of 158, 961 bp long with 36.78% GC content, among which a pair of inverted repeats (IRs) of 26, 671 bp separated a large single-copy (LSC) region of 87, 162 bp and a small single-copy (SSC) region of 18, 457 bp. Altogether 131 genes were hosted, including 37 transfer RNAs, 8 ribosomal RNAs, and 86 protein-coding genes. 284 simple sequence repeats (SSRs) were also marked out. A comparative analysis of the genome structure and the sequence data of closely related species unveiled the conserved gene order in the IR and LSC/SSC regions, a quite constructive finding for future phylogenetic research. Overall, this study providing C. rigidus R.Br. genomic resources could positively contribute to the evolutionary study and the phylogenetic reconstruction of Myrtaceae.

Comparative Plastome Analyses and Phylogenetic Applications of the Acer Section Platanoidea

The Acer L. (Sapindaceae) is one of the most diverse and widespread genera in the Northern Hemisphere. Section Platanoidea harbours high genetic and morphological diversity and shows the phylogenetic conflict between A. catalpifolium and A. amplum. Chloroplast (cp) genome sequencing is efficient for the enhancement of the understanding of phylogenetic relationships and taxonomic revision. Here, we report complete cp genomes of five species of Acer sect. Platanoidea. The length of Acer sect. Platanoidea cp genomes ranged from 156,262 bp to 157,349 bp and detected the structural variation in the inverted repeats (IRs) boundaries. By conducting a sliding window analysis, we found that five relatively high variable regions (trnH-psbA, psbN-trnD, psaA-ycf3, petA-psbJ and ndhA intron) had a high potential for developing effective genetic markers. Moreover, with an addition of eight plastomes collected from GenBank, we displayed a robust phylogenetic tree of the Acer sect. Platanoidea, with high resolutions for nearly all identified nodes, suggests a promising opportunity to resolve infrasectional relationships of the most species-rich section Platanoidea of Acer.

The complete chloroplast genome of Myriophyllum spicatum reveals a 4-kb inversion and new insights regarding plastome evolution in Haloragaceae

Myriophyllum, among the most species-rich genera of aquatic angiosperms with ca. 68 species, is an extensively distributed hydrophyte lineage in the cosmopolitan family Haloragaceae. The chloroplast (cp) genome is useful in the study of genetic evolution, phylogenetic analysis, and molecular dating of controversial taxa. Here, we sequenced and assembled the whole chloroplast genome of Myriophyllum spicatum L. and compared it to other species in the order Saxifragales. The complete chloroplast genome sequence of M. spicatum is 158,858 bp long and displays a quadripartite structure with two inverted repeats (IR) separating the large single copy (LSC) region from the small single copy (SSC) region. Based on sequence identification and the phylogenetic analysis, a 4-kb phylogenetically informative inversion between trnE-trnC in Myriophyllum was determined, and we have placed this inversion on a lineage specific to Myriophyllum and its close relatives. The divergence time estimation suggested that the trnE-trnC inversion possibly occurred between the upper Cretaceous (72.54 MYA) and middle Eocene (47.28 MYA) before the divergence of Myriophyllum from its most recent common ancestor. The unique 4-kb inversion might be caused by an occurrence of nonrandom recombination associated with climate changes around the K-Pg boundary, making it interesting for future evolutionary investigations.

Characterization of the complete chloroplast genome of sugar maple (Acer saccharum)

Acer saccharum is one ecologically and economically important tree species cultivated widely across the world. In this study we generated the complete chloroplast (cp) genome of A. saccharum via genome-skimming method. The assembled genome is 155,684 base-pairs (bp) in size, with one large single copy region of 85,393 bp and one small single copy region of 18,033 bp separated by two inverted repeats of 26,129 bp. The genome contains a total of 133 genes, including 85 protein-coding genes, 8 rRNAs and 40 tRNAs. Furthermore, phylogenomic estimation strongly supported A. saccharum as a distinct lineage within the monophyletic Acer.

Species Identification of Dracaena Using the Complete Chloroplast Genome as a Super-Barcode

The taxonomy and nomenclature of Dracaena plants are much disputed, particularly for several Dracaena species in Asia. However, neither morphological features nor common DNA regions are ideal for identification of Dracaena spp. Meanwhile, although multiple Dracaena spp. are sources of the rare traditional medicine dragon's blood, the Pharmacopoeia of the People's Republic of China has defined Dracaena cochinchinensis as the only source plant. The inaccurate identification of Dracaena spp. will inevitably affect the clinical efficacy of dragon's blood. It is therefore important to find a better method to distinguish these species. Here, we report the complete chloroplast (CP) genomes of six Dracaena spp., D. cochinchinensis, D. cambodiana, D. angustifolia, D. terniflora, D. hokouensis, and D. elliptica, obtained through high-throughput Illumina sequencing. These CP genomes exhibited typical circular tetramerous structure, and their sizes ranged from 155,055 (D. elliptica) to 155,449 bp (D. cochinchinensis). The GC content of each CP genome was 37.5%. Furthermore, each CP genome contained 130 genes, including 84 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. There were no potential coding or non-coding regions to distinguish these six species, but the maximum likelihood tree of the six Dracaena spp. and other related species revealed that the whole CP genome can be used as a super-barcode to identify these Dracaena spp. This study provides not only invaluable data for species identification and safe medical application of Dracaena but also an important reference and foundation for species identification and phylogeny of Liliaceae plants.

Comparative analysis of the complete chloroplast genome sequences of six species of Pulsatilla Miller, Ranunculaceae

BACKGROUND

Baitouweng is a traditional Chinese medicine with a long history of different applications. Although referred to as a single medicine, Baitouweng is actually comprised of many closely related species. It is therefore critically important to identify the different species that are utilized in these medicinal applications. Knowledge about their phylogenetic relationships can be derived from their chloroplast genomes and may provide additional insights into development of molecular markers.

METHODS

Genomic DNA was extracted from six species of Pulsatilla and then sequenced on an Illumina HiSeq 4000. Sequences were assembled into contigs by SOAPdenovo 2.04, aligned to the reference genome using BLAST, and then manually corrected. Genome annotation was performed by the online DOGMA tool. General characteristics of the cp genomes of the six species were analyzed and compared with closely related species. Additionally, phylogenetic trees were constructed, based on single nucleotide polymorphisms (SNPs) and 51 shared protein-coding gene sequences in the cp genome among all 31 species via maximum likelihood.

RESULTS

The size of cp genomes of P. chinensis (Bge.) Regel, P. chinensis (Bge.) Regel var. kissii (Mandl) S. H. Li et Y. H. Huang, P. cernua (Thunb.) Bercht. et Opiz f. plumbea J. X. Ji et Y. T. zhao, P. dahurica (Fisch.) Spreng, P. turczaninovii Kryl. et Serg, and P. cernua (Thunb.) Bercht. et Opiz. were 163,851 bp, 163,756 bp, 162,481 bp, 162,450 bp, 162,795 bp, and 162,924 bp, respectively. Each species included two inverted repeat regions, a small single-copy region, and a large single-copy region. A total of 134 genes were annotated, including 90 protein-coding genes, 36 tRNAs, and eight rRNAs across all species. In simple sequence repeat analysis, only P. dahurica was found to contain hexanucleotide repeats. A total of 26, 39, 32, 37, 32 and 43 large repeat sequences were identified in the genic regions of the six Pulsatilla species. Nucleotide diversity analysis revealed that the rpl36 gene and ccsA-ndhD region have the highest Pi value. In addition, two phylogenetic trees of the cp genomes were constructed, which laced all Pulsatilla species into one branch within Ranunculaceae.

CONCLUSIONS

We identified and analyzed the cp genome features of six species of P. Miller, with implications for species identification and phylogenetic analysis.

Characterization of the Complete Chloroplast Genome of Acer truncatum Bunge (Sapindales: Aceraceae): A New Woody Oil Tree Species Producing Nervonic Acid

Acer truncatum, which is a new woody oil tree species, is an important ornamental and medicinal plant in China. To assess the genetic diversity and relationships of A. truncatum, we analyzed its complete chloroplast (cp) genome sequence. The A. truncatum cp genome comprises 156,492 bp, with the large single-copy, small single-copy, and inverted repeat (IR) regions consisting of 86,010, 18,050, and 26,216 bp, respectively. The A. truncatum cp genome contains 112 unique functional genes (i.e., 4 rRNA, 30 tRNA, and 78 protein-coding genes) as well as 78 simple sequence repeats, 9 forward repeats, 1 reverse repeat, 5 palindromic repeats, and 7 tandem repeats. We analyzed the expansion/contraction of the IR regions in the cp genomes of six Acer species. A comparison of these cp genomes indicated the noncoding regions were more diverse than the coding regions. A phylogenetic analysis revealed that A. truncatum is closely related to A. miaotaiense. Moreover, a novel ycf4-cemA indel marker was developed for distinguishing several Acer species (i.e., A. buergerianum, A. truncatum, A. henryi, A. negundo, A. ginnala, and A. tonkinense). The results of the current study provide valuable information for future evolutionary studies and the molecular barcoding of Acer species.

The Complete Plastid Genome of Magnolia zenii and Genetic Comparison to Magnoliaceae species

Magnolia zenii is a critically endangered species known from only 18 trees that survive on Baohua Mountain in Jiangsu province, China. Little information is available regarding its molecular biology, with no genomic study performed on M. zenii until now. We determined the complete plastid genome of M. zenii and identified microsatellites. Whole sequence alignment and phylogenetic analysis using BI and ML methods were also conducted. The plastome of M. zenii was 160,048 bp long with 39.2% GC content and included a pair of inverted repeats (IRs) of 26,596 bp that separated a large single-copy (LSC) region of 88,098 bp and a small single-copy (SSC) region of 18,757 bp. One hundred thirty genes were identified, of which 79 were protein-coding genes, 37 were transfer RNAs, and eight were ribosomal RNAs. Thirty seven simple sequence repeats (SSRs) were also identified. Comparative analyses of genome structure and sequence data of closely-related species revealed five mutation hotspots, useful for future phylogenetic research. Magnolia zenii was placed as sister to M. biondii with strong support in all analyses. Overall, this study providing M. zenii genomic resources will be beneficial for the evolutionary study and phylogenetic reconstruction of Magnoliaceae.