Methods for obtaining and analyzing whole chloroplast genome sequences

Comparative chloroplast genomics and phylogenetic analysis of Oreomecon nudicaulis (Papaveraceae)

Oreomecon nudicaulis, commonly known as mountain poppy, is a significant perennial herb. In 2022, the species O. nudicaulis, which was previously classified under the genus Papaver, was reclassified within the genus Oreomecon. Nevertheless, the phylogenetic status and chloroplast genome within the genus Oreomecon have not yet been reported. This study elucidates the chloroplast genome sequence and structural features of O. nudicaulis and explores its evolutionary relationships within Papaveraceae. Using Illumina sequencing technology, the chloroplast genome of O. nudicaulis was sequenced, assembled, and annotated. The results indicate that the chloroplast genome of O. nudicaulis exhibits a typical circular quadripartite structure. The chloroplast genome is 153,903 bp in length, with a GC content of 38.87%, containing 84 protein-coding genes, 8 rRNA genes, 38 tRNA genes, and 2 pseudogenes. The genome encodes 25,815 codons, with leucine (Leu) being the most abundant codon, and the most frequently used codon is AUU. Additionally, 129 microsatellite markers were identified, with mononucleotide repeats being the most abundant (53.49%). Our phylogenetic analysis revealed that O. nudicaulis has a relatively close relationship with the genus Meconopsis within the Papaveraceae family. The phylogenetic analysis supported the taxonomic status of O. nudicaulis, as it did not form a clade with other Papaver species, consistent with the revised taxonomy of Papaveraceae. This is the first report of a phylogenomic study of the complete chloroplast genome in the genus Oreomecon, which is a significant genus worldwide. This analysis of the O. nudicaulis chloroplast genome provides a theoretical basis for research on genetic diversity, molecular marker development, and species identification, enriching genetic information and supporting the evolutionary relationships among Papaveraceae.

The complete chloroplast genome of Oxytropis ochrocephala Bunge 1874 (Fabaceae) and its phylogenetic analysis

This study presents the first-ever complete chloroplast (cp) genome sequence of Oxytropis ochrocephala Bunge 1874, a member of the Fabaceae family. The cp genome spans 126,996 base pairs and includes 109 genes, comprising 76 protein-coding genes, 29 tRNA genes, and four rRNA genes. Notably, the genome lacks an inverted repeat (IR) region. Additionally, we constructed phylogenetic trees for 34 species within Trib. Galegeae, employing both maximum likelihood (ML) and Bayesian inference (BI) methods. These analyses robustly support the monophyly of the Oxytropis species, evidenced by high bootstrap values (BP = 100) and posterior probabilities (PP = 1). Within this clade, O. ochrocephala exhibits a sister relationship with other Oxytropis species. Our findings provide valuable insights into the genetic makeup and evolutionary relationships of O. ochrocephala within the Galegeae tribe.

The complete chloroplast genome of Mussaenda pubescens and phylogenetic analysis

The chloroplast (cp) genome sequence of Mussaenda pubescens, a promising resource that is used as a traditional medicine and drink, is important for understanding the phylogenetic relationships among the Mussaenda family and genetic improvement and reservation. This research represented the first comprehensive description of the morphological characteristics of M. pubescens, as well as an analysis of the complete cp genome and phylogenetic relationship. The results indicated a close relationship between M. pubescens and M. hirsutula based on the morphological characteristics of the flower and leaves. The cp was sequenced using the Illumina NovaSeq 6000 platform. The results indicated the cp genome of M. pubescens spanned a total length of 155,122 bp, including a pair of inverted repeats (IRA and IRB) with a length of 25,871 bp for each region, as well as a large single-copy (LSC) region and a small single-copy (SSC) region with lengths of 85,370 bp and 18,010 bp, respectively. The results of phylogenetic analyses demonstrated that species within the same genus displayed a tendency to group closely together. It was suggested that Antirhea, Cinchona, Mitragyna, Neolamarckia, and Uncaria might have experienced an early divergence. Furthermore, M. hirsutula showed a close genetic connection to M. pubescens, with the two species having partially overlapping distributions in China. This study presents crucial findings regarding the identification, evolution, and phylogenetic research on Mussaenda plants, specifically targeting M. pubescens.

Extraction and analysis of high-quality chloroplast DNA with reduced nuclear DNA for medicinal plants

BACKGROUND

Obtaining high-quality chloroplast genome sequences requires chloroplast DNA (cpDNA) samples that meet the sequencing requirements. The quality of extracted cpDNA directly impacts the efficiency and accuracy of sequencing analysis. Currently, there are no reported methods for extracting cpDNA from Erigeron breviscapus. Therefore, we developed a suitable method for extracting cpDNA from E. breviscapus and further verified its applicability to other medicinal plants.

RESULTS

We conducted a comparative analysis of chloroplast isolation and cpDNA extraction using modified high-salt low-pH method, the high-salt method, and the NaOH low-salt method, respectively. Subsequently, the number of cpDNA copies relative to the nuclear DNA (nDNA ) was quantified via qPCR. As anticipated, chloroplasts isolated from E. breviscapus using the modified high-salt low-pH method exhibited intact structures with minimal cell debris. Moreover, the concentration, purity, and quality of E. breviscapus cpDNA extracted through this method surpassed those obtained from the other two methods. Furthermore, qPCR analysis confirmed that the modified high-salt low-pH method effectively minimized nDNA contamination in the extracted cpDNA. We then applied the developed modified high-salt low-pH method to other medicinal plant species, including Mentha haplocalyx, Taraxacum mongolicum, and Portulaca oleracea. The resultant effect on chloroplast isolation and cpDNA extraction further validated the generalizability and efficacy of this method across different plant species.

CONCLUSIONS

The modified high-salt low-pH method represents a reliable approach for obtaining high-quality cpDNA from E. breviscapus. Its universal applicability establishes a solid foundation for chloroplast genome sequencing and analysis of this species. Moreover, it serves as a benchmark for developing similar methods to extract chloroplast genomes from other medicinal plants.

Comparative analysis of chloroplast genomes of Pulsatilla species reveals evolutionary and taxonomic status of newly discovered endangered species Pulsatilla saxatilis

BACKGROUND

Pulsatilla saxatilis, a new species of the genus Pulsatilla has been discovered. The morphological information of this species has been well described, but its chloroplast genome characteristics and comparison with species of the same genus remain to be reported.

RESULTS

Our results showed that the total length of chloroplast (cp.) genome of P. saxatilis is 162,659 bp, with a GC content of 37.5%. The cp. genome contains 134 genes, including 90 known protein-coding genes, 36 tRNA genes, and 8 rRNA genes. P. saxatilis demonstrated similar characteristics to other species of genus Pulsatilla. Herein, we compared cp. genomes of 10 species, including P. saxatilis, and found that the cp. genomes of the genus Pulsatilla are extremely similar, with a length of 162,322-163,851 bp. Furthermore, The SSRs of Pulsatilla ranged from 10 to 22 bp in length. Among the four structural regions of the cp. genome, most long repeats and SSRs were detected in the LSC region, followed by that in the SSC region, and least in IRA/ IRB regions. The most common types of long repeats were forward and palindromic repeats, followed by reverse repeats, and only a few complementary repeats were found in 10 cp. genomes. We also analyzed nucleotide diversity and identified ccsA_ndhD, rps16_trnK-UUU, ccsA, and rbcL, which could be used as potential molecular markers for identification of Pulsatilla species. The results of the phylogenetic tree constructed by connecting the sequences of high variation regions were consistent with those of the cp. gene phylogenetic tree, and the species more closely related to P. saxatilis was identified as the P. campanella.

CONCLUSION

It was determined that the closest species to P. saxatilis is P. campanella, which is the same as the conclusion based on pollen grain characteristics, but different from the P. chinensis determined based on morphological characteristics. By revealing information on the chloroplast characteristics, development, and evolution of the cp. genome and the potential molecular markers, this study provides effective molecular data regarding the evolution, genetic diversity, and species identification of the genus Pulsatilla.

Design and assembly of the 117-kb Phaeodactylum tricornutum chloroplast genome

There is growing impetus to expand the repertoire of chassis available to synthetic biologists. Chloroplast genomes present an interesting alternative for engineering photosynthetic eukaryotes; however, development of the chloroplast as a synthetic biology chassis has been limited by a lack of efficient techniques for whole-genome cloning and engineering. Here, we demonstrate two approaches for cloning the 117-kb Phaeodactylum tricornutum chloroplast genome that have 90% to 100% efficiency when screening as few as 10 yeast (Saccharomyces cerevisiae) colonies following yeast assembly. The first method reconstitutes the genome from PCR-amplified fragments, whereas the second method involves precloning these fragments into individual plasmids from which they can later be released. In both cases, overlapping fragments of the chloroplast genome and a cloning vector are homologously recombined into a singular contig through yeast assembly. The cloned chloroplast genome can be stably maintained and propagated within Escherichia coli, which provides an exciting opportunity for engineering a delivery mechanism for bringing DNA directly to the algal chloroplast. Also, one of the cloned genomes was designed to contain a single SapI site within the yeast URA3 (coding for orotidine-5'-phosphate decarboxylase) open-reading frame, which can be used to linearize the genome and integrate designer cassettes via golden-gate cloning or further iterations of yeast assembly. The methods presented here could be extrapolated to other species-particularly those with a similar chloroplast genome size and architecture (e.g. Thalassiosira pseudonana).

The complete plastome genome sequence of Cynanchum otophyllum (Asclepiadaceae), a unique medicinal species in China

Cynanchum otophyllum Schneid is an important medicinal plant in China. In this paper, the chloroplast genome of C. otophyllum was sequenced based on high-throughput technology, and the chloroplast genome structure characteristics and phylogenetic relationship of C. otophyllum were analyzed. The results showed the complete plastome genome size of C. otophyllumis 160,874bp, including one small single copy (SSC, 19,851bp) and one large single copy (LSC, 92,009bp) regions isolated by a pair of inverted repeat regions (IRs, 24,507bp). The whole plastome genome including 84 protein encoding genes, 8 rRNA and 37 tRNA. Based on the phylogenetic topologies, C. otophyllum shows close association with additional Gomphocarpus and Asclepias genus. This study contributes to an enhanced understanding of the genetic information of C. otophyllum and provides a theoretical basis for the development of molecular markers and phylogeographic of the species, as well as for constructing the phylogenetic tree of Asclepiadaceae.

A Scalable and Robust Chloroplast Genotyping Solution: Development and Application of SNP and InDel Markers in the Maize Chloroplast Genome

Maize(Zea mays. L) is a globally important crop, and understanding its genetic diversity is crucial for plant breeding phylogenetic analyses and comparative genetics. While nuclear markers have been extensively used for mapping agriculturally important genes, they are limited in recognizing characteristics, such as cytoplasmic male sterility and reciprocal cross hybrids. In this study, we performed next-generation sequencing of 176samples, and the maize cultivars represented five distinct groups. A total of 89 single nucleotide polymorphisms (SNPs) and 11 insertion/deletion polymorphisms (InDels) were identified. To enable high-throughput detection, we successfully amplified and confirmed 49 SNP and InDel markers, which were defined as a Varietal Chloroplast Panel (VCP) using the Kompetitive Allele Specific PCR (KASP). The specific markers provided a valuable tool for identifying chloroplast groups. The verification experiment, focusing on the identification of reciprocal cross hybrids and cytoplasmic male sterility hybrids, demonstrated the significant advantages of VCP markers in maternal inheritance characterization. Furthermore, only a small subset of these markers is needed to provide useful information, showcasing the effectiveness of these markers in elucidating the artificial selection process of elite maize lines.

Comparative Analysis of Chloroplast Pan-Genomes and Transcriptomics Reveals Cold Adaptation in Medicago sativa

Alfalfa (Medicago sativa) is a perennial forage legume that is widely distributed all over the world; therefore, it has an extremely complex genetic background. Though population structure and phylogenetic studies have been conducted on a large group of alfalfa nuclear genomes, information about the chloroplast genomes is still lacking. Chloroplast genomes are generally considered to be conservative and play an important role in population diversity analysis and species adaptation in plants. Here, 231 complete alfalfa chloroplast genomes were successfully assembled from 359 alfalfa resequencing data, on the basis of which the alfalfa chloroplast pan-genome was constructed. We investigated the genetic variations of the alfalfa chloroplast genome through comparative genomic, genetic diversity, phylogenetic, population genetic structure, and haplotype analysis. Meanwhile, the expression of alfalfa chloroplast genes under cold stress was explored through transcriptome analysis. As a result, chloroplast genomes of 231 alfalfa lack an IR region, and the size of the chloroplast genome ranges from 125,192 bp to 126,105 bp. Using population structure, haplotypes, and construction of a phylogenetic tree, it was found that alfalfa populations could be divided into four groups, and multiple highly variable regions were found in the alfalfa chloroplast genome. Transcriptome analysis showed that tRNA genes were significantly up-regulated in the cold-sensitive varieties, while rps7, rpl32, and ndhB were down-regulated, and the editing efficiency of ycf1, ycf2, and ndhF was decreased in the cold-tolerant varieties, which may be due to the fact that chloroplasts store nutrients through photosynthesis to resist cold. The huge number of genetic variants in this study provide powerful resources for molecular markers.

Comparative Phylogenetic Analysis of Ancient Korean Tea "Hadong Cheon-Nyeon Cha (Camellia sinensis var. sinensis)" Using Complete Chloroplast Genome Sequences

Wild teas are valuable genetic resources for studying evolution and breeding. Here, we report the complete chloroplast genome of the ancient Korean tea 'Hadong Cheon-nyeon Cha' (C. sinensis var. sinensis), which is known as the oldest tea tree in Korea. This study determined seven Camellia sinensis var. sinenesis, including Hadong Cheon-nyeon Cha (HCNC) chloroplast genome sequences, using Illumina sequencing technology via de novo assembly. The chloroplast genome sizes ranged from 157,019 to 157,114 bp and were organized into quadripartite regions with the typical chloroplast genomes. Further, differences in SNPs and InDels were detected across the seven chloroplast genomes through variance analysis. Principal component and phylogenetic analysis suggested that regional constraints, rather than functional constraints, strongly affected the sequence evolution of the cp genomes in this study. These genomic resources provide evolutionary insight into Korean tea plant cultivars and lay the foundation for a better understanding of the ancient Korean tea plant HCNC.

A novel in-situ-process technique constructs whole circular cpDNA library

BACKGROUND

The chloroplast genome (cp genome) is directly related to the study and analysis of molecular phylogeny and evolution of plants in the phylogenomics era. The cp genome, whereas, is highly plastic and exists as a heterogeneous mixture of sizes and physical conformations. It is advantageous to purify/enrich the circular chloroplast DNA (cpDNA) to reduce sequence complexity in cp genome research. Large-insert, ordered DNA libraries are more practical for genomics research than conventional, unordered ones. From this, a technique of constructing the ordered BAC library with the goal-insert cpDNA fragment is developed in this paper.

RESULTS

This novel in-situ-process technique will efficiently extract circular cpDNA from crops and construct a high-quality cpDNA library. The protocol combines the in-situ chloroplast lysis for the high-purity circular cpDNA with the in-situ substitute/ligation for the high-quality cpDNA library. Individually, a series of original buffers/solutions and optimized procedures for chloroplast lysis in-situ is different than bacterial lysis in-situ; the in-situ substitute/ligation that reacts on the MCE membrane is suitable for constructing the goal-insert, ordered cpDNA library while preventing the large-insert cpDNA fragment breakage. The goal-insert, ordered cpDNA library is arrayed on the microtiter plate by three colonies with the definite cpDNA fragment that are homologous-corresponds to the whole circular cpDNA of the chloroplast.

CONCLUSION

The novel in-situ-process technique amply furtherance of research in genome-wide functional analysis and characterization of chloroplasts, such as genome sequencing, bioinformatics analysis, cloning, physical mapping, molecular phylogeny and evolution.

The complete chloroplast genome of Cicer reticulatum and comparative analysis against relative Cicer species

The chloroplast (cp) genome is an adequate genomic resource to investigate evolutionary relationships among plant species and it carries marker genes available for species identification. The Cicer reticulatum is one of perennial species as the progenitor of cultivated chickpeas. Although a large part of the land plants has a quadruple chloroplast genome organization, the cp genome of C. reticulatum consists of one LSC (Large Single Copy Region), one SSC (Small Single Copy Region), and one IR (Inverted Repeat) region, which indicates that it has an untypical and unique structure. This type of chloroplast genome belongs to the IR-lacking clade. Chloroplast DNA (cpDNA) was extracted from fresh leaves using a high salt-based protocol and sequencing was performed using DNA Nanoball Sequencing technology. The comparative analysis employed between the species to examine genomic differences and gene homology. The study also included codon usage frequency analysis, hotspot divergence analysis, and phylogenetic analysis using various bioinformatics tools. The cp genome of C. reticulatum was found 125,794 bp in length, with an overall GC content of 33.9%. With a total of 79 protein-coding genes, 34 tRNA genes, and 4 rRNA genes. Comparative genomic analysis revealed 99.93% similarity between C. reticulatum and C. arietinum. Phylogenetic analysis further indicated that the closest evolutionary relative to C. arietinum was C. reticulatum, whereas the previously sequenced wild Cicer species displayed slight distinctions across their entire coding regions. Several genomic regions, such as clpP and ycf1, were found to exhibit high nucleotide diversity, suggesting their potential utility as markers for investigating the evolutionary relationships within the Cicer genus. The first complete cp genome sequence of C. reticulatum will provide novel insights for future genetic research on Cicer crops.

Comparative analysis of the complete chloroplast genome of Papaveraceae to identify rearrangements within the Corydalis chloroplast genome

Chloroplast genomes are valuable for inferring evolutionary relationships. We report the complete chloroplast genomes of 36 Corydalis spp. and one Fumaria species. We compared these genomes with 22 other taxa and investigated the genome structure, gene content, and evolutionary dynamics of the chloroplast genomes of 58 species, explored the structure, size, repeat sequences, and divergent hotspots of these genomes, conducted phylogenetic analysis, and identified nine types of chloroplast genome structures among Corydalis spp. The ndh gene family suffered inversion and rearrangement or was lost or pseudogenized throughout the chloroplast genomes of various Corydalis species. Analysis of five protein-coding genes revealed simple sequence repeats and repetitive sequences that can be potential molecular markers for species identification. Phylogenetic analysis revealed three subgenera in Corydalis. Subgenera Cremnocapnos and Sophorocapnos represented the Type 2 and 3 genome structures, respectively. Subgenus Corydalis included all types except type 3, suggesting that chloroplast genome structural diversity increased during its differentiation. Despite the explosive diversification of this subgenus, most endemic species collected from the Korean Peninsula shared only one type of genome structure, suggesting recent divergence. These findings will greatly improve our understanding of the chloroplast genome of Corydalis and may help develop effective molecular markers.

Comparative Analysis of the Chloroplast Genome of Sicyos angulatus with Other Seven Species of Cucurbitaceae Family

Sicyos angulatus (SA) is an annual plant from the Cucurbitaceae family that is native to the eastern part of North America. This study aims to assemble and annotate the chloroplast genome of S. angulatus, and then compare it with plastomes of the other species representing the Cucurbitaceae family. The chloroplast genome size of S. angulatus is 154,986 bp, including a pair of inverted repeats (IR) of 26,276 bp, and small single-copy region (SSC) of 18,079 bp and large single-copy region (LSC) of 84,355 bp. Compared to other Cucurbitaceae species, the chloroplast genome of S. angulatus is almost 4222 bp smaller than the plastome Gynostemma pentaphyllum. All other seven species have an identical set of tRNA (37), except Citrullus laevigata, which contains 36 tRNA. The IRa/LSC junction in all eight species is located upstream of rpl2 and downstream of trnH gene. Moreover, variation in the size of the gene and the presence of pseudogene ycf1 has been seen because of the IR contraction and expansion. The highest number of tandem repeats was seen in G. pentaphyllum, and then Corynocarpus leavigata. The sequence divergence analysis and topology of the phylogenetic tree indicate that S. angulatus is more similar to genus Citrullus as compared to genus Gynostemma. These findings contribute to developing the genomic marker for the purpose of future genetic studies.

The First Complete Chloroplast Genome of Campanula carpatica: Genome Characterization and Phylogenetic Diversity

Campanula carpatica is an ornamental flowering plant belonging to the family Campanulaceae. The complete chloroplast genome of C. carpatica was obtained using Illumina HiSeq X and Oxford Nanopore (Nanopore GridION) platforms. The chloroplast genome exhibited a typical circular structure with a total length of 169,341 bp, comprising a large single-copy region of 102,323 bp, a small single-copy region of 7744 bp, and a pair of inverted repeats (IRa/IRb) of 29,637 bp each. Out of a total 120 genes, 76 were protein-coding genes, 36 were transfer RNA genes, and eight were ribosomal RNA genes. The genomic characteristics of C. carpatica are similar to those of other Campanula species in terms of repetitive sequences, sequence divergence, and contraction/expansion events in the inverted repeat regions. A phylogenetic analysis of 63 shared genes in 16 plant species revealed that Campanula zangezura is the closest relative of C. carpatica. Phylogenetic analysis indicated that C. carpatica was within the Campanula clade, and C. pallida occupied the outermost position of that clade.

Comparative analysis of dioecious Amaranthus plastomes and phylogenomic implications within Amaranthaceae s.s

BACKGROUND

The genus Amaranthus L. consists of 70-80 species distributed across temperate and tropical regions of the world. Nine species are dioecious and native to North America; two of which are agronomically important weeds of row crops. The genus has been described as taxonomically challenging and relationships among species including the dioecious ones are poorly understood. In this study, we investigated the phylogenetic relationships among the dioecious amaranths and sought to gain insights into plastid tree incongruence. A total of 19 Amaranthus species' complete plastomes were analyzed. Among these, seven dioecious Amaranthus plastomes were newly sequenced and assembled, an additional two were assembled from previously published short reads sequences and 10 other plastomes were obtained from a public repository (GenBank).

RESULTS

Comparative analysis of the dioecious Amaranthus species' plastomes revealed sizes ranged from 150,011 to 150,735 bp and consisted of 112 unique genes (78 protein-coding genes, 30 transfer RNAs and 4 ribosomal RNAs). Maximum likelihood trees, Bayesian inference trees and splits graphs support the monophyly of subgenera Acnida (7 dioecious species) and Amaranthus; however, the relationship of A. australis and A. cannabinus to the other dioecious species in Acnida could not be established, as it appears a chloroplast capture occurred from the lineage leading to the Acnida + Amaranthus clades. Our results also revealed intraplastome conflict at some tree branches that were in some cases alleviated with the use of whole chloroplast genome alignment, indicating non-coding regions contribute valuable phylogenetic signals toward shallow relationship resolution. Furthermore, we report a very low evolutionary distance between A. palmeri and A. watsonii, indicating that these two species are more genetically related than previously reported.

CONCLUSIONS

Our study provides valuable plastome resources as well as a framework for further evolutionary analyses of the entire Amaranthus genus as more species are sequenced.

Co-evolution of large inverted repeats and G-quadruplex DNA in fungal mitochondria may facilitate mitogenome stability: the case of Malassezia

Mitogenomes are essential due to their contribution to cell respiration. Recently they have also been implicated in fungal pathogenicity mechanisms. Members of the basidiomycetous yeast genus Malassezia are an important fungal component of the human skin microbiome, linked to various skin diseases, bloodstream infections, and they are increasingly implicated in gut diseases and certain cancers. In this study, the comparative analysis of Malassezia mitogenomes contributed to phylogenetic tree construction for all species. The mitogenomes presented significant size and gene order diversity which correlates to their phylogeny. Most importantly, they showed the inclusion of large inverted repeats (LIRs) and G-quadruplex (G4) DNA elements, rendering Malassezia mitogenomes a valuable test case for elucidating the evolutionary mechanisms responsible for this genome diversity. Both LIRs and G4s coexist and convergently evolved to provide genome stability through recombination. This mechanism is common in chloroplasts but, hitherto, rarely found in mitogenomes.

Plastid Phylogenomic Analyses Reveal a Cryptic Species of Ligusticopsis (Apiaceae, Angiosperms)

Ligusticopsis litangensis is identified and described as a cryptic species from Sichuan Province, China. Although the distribution of this cryptic species overlaps with that of Ligusticopsis capillacea and Ligusticopsis dielsiana, the morphological boundaries between them are explicit and have obviously distinguishable characters. The main distinguishing features of the cryptic species are as follows: long conical multi-branched roots, very short pedicels in compound umbels, unequal rays, oblong-globose fruits, 1-2 vittae per furrow and 3-4 vittae on the commissure. The above-mentioned features differ somewhat from other species within the genus Ligusticopsis, but generally coincide with the morphological boundaries defined for the genus Ligusticopsis. To determine the taxonomic position of L. litangensis, we sequenced and assembled the plastomes of L. litangensis and compared them with the plastomes of 11 other species of the genus Ligusticopsis. Notably, both phylogenetic analyses based on ITS sequences and the complete chloroplast genome robustly supported that three accessions of L. litangensis are monophyletic clade and then nested in Ligusticopsis genus. Moreover, the plastid genomes of 12 Ligusticopsis species, including the new species, were highly conserved in terms of gene order, gene content, codon bias, IR boundaries and SSR content. Overall, the integration of morphological, comparative genomic and phylogenetic evidence indicates that Ligusticopsis litangensis actually represents a new species.

Characteristics, Comparative Analysis, and Phylogenetic Relationships of Chloroplast Genomes of Cultivars and Wild Relatives of Eggplant (Solanum melongena)

The eggplant (Solanum melongena) is a popular vegetable around the world. However, the origin and evolution of eggplant has long been considered complex and unclear, which has become the barrier to improvements in eggplant breeding. Sequencing and comparative analyses of 13 complete chloroplast (cp) genomes of seven Solanum species were performed. Genome sizes were between 154,942 and 156,004 bp, the smallest genome was from S. torvum and the largest from S. macrocapon. Thirteen cp genomes showed highly conserved sequences and GC contents, particularly at the subgenus level. All genes in the 13 genomes were annotated. The cp genomes in this study comprised 130 genes (i.e., 80 protein-coding genes, 8 rRNA genes, and 42 tRNA genes), apart from S. sisymbriifolium, which had 129 (79 protein-coding genes, 8 rRNA genes, and 42 tRNA genes.). The rps16 was absent from the cp genome of S. sisymbriifolium, resulting in a nonsense mutation. Twelve hotspot regions of the cp genome were identified, which showed a series of sequence variations and differed significantly in the inverted repeat/single-copy boundary regions. Furthermore, phylogenetic analysis was conducted using 46 cp genomic sequences to determine interspecific genetic and phylogenetic relationships in Solanum species. All species formed two branches, one of which contained all cultivars of the subgenus Leptostemonum. The cp genome data and phylogenetic analysis provides molecular evidence revealing the origin and evolutionary relationships of S. melongena and its wild relatives. Our findings suggest precise intra- and interspecies relatedness within the subgenus Leptostemonum, which has positive implications for work on improvements in eggplant breeding, particularly in producing heterosis, expanding the source of species variation, and breeding new varieties.

Genetic relationship of Pleione based on the chloroplast genome

Pleione (Orchidaceae) is a very famous horticultural plant with a high international reputation for its unique flower shape and abundant color. The small difference in morphological characteristics among Pleione species caused by weak reproductive isolation and easy hybridization makes the taxonomic status of individual species very confusing. Chloroplast (cp) genome information is of great significance for the study of plant phylogeny and taxonomy. In this study, the cp genomes of Pleione were sequenced and the complete cp structure and sequence characteristics of 19 species were compared and analyzed. The cp genome of Pleione species exhibited a conserved tetrad structure and each species encoded 135 protein-coding genes, 38 tRNA and 8 RNA genes. The cp genome sizes of 19 Pleione were 157964-159269 bp and the length of LSC, SSC and IR were 85953-87172 bp, 18499-18712 bp, 26459-26756 bp, respectively. Palindromic and forward repeats accounted for a high proportion and the SSRs were mainly mononucleotide repeats in Pleione. Analysis of chloroplast sequence differences indicated that the differences in coding regions were smaller than those in non-coding regions, and the variation in LSC and SSC regions was greater than that in IR regions. Phylogenetic analysis showed that all Pleione species inferred from the cp genome were clustered together and received high support. However, the genetic relationship of Pleione plants is different from the current update system of this genus. Therefore, the demarcation of Pleione interspecific relationships still needs further investigation due to the lack of sufficient evidence. The cp genome serves as valuable information for the identification of Pleione species and the study of phylogenetic relationships.

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

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

Poaceae Chloroplast Genome Sequencing: Great Leap Forward in Recent Ten Years

The first complete chloroplast genome of rice (Oryza sativa) was published in 1989, ushering in a new era of studies of chloroplast genomics in Poaceae. Progresses in Next-Generation Sequencing (NGS) and Third-Generation Sequencing (TGS) technologiesand in the development of genome assembly software, have significantly advanced chloroplast genomics research. Poaceae is one of the most targeted families in chloroplast genome research because of its agricultural, ecological, and economic importance. Over the last 30 years, 2,050 complete chloroplast genome sequences from 40 tribes and 282 genera have been generated, most (97%) of them in the recent ten years. The wealth of data provides the groundwork for studies on species evolution, phylogeny, genetic transformation, and other aspects of Poaceae chloroplast genomes. As a result, we have gained a deeper understanding of the properties of Poaceae chloroplast genomes. Here, we summarize the achievements of the studies of the Poaceae chloroplast genomes and envision the challenges for moving the area ahead.

New Insights into Plastid and Mitochondria Evolution in Wild Peas (Pisum L.)

Plastids and mitochondria are organelles of plant cells with small genomes, which may exhibit discordant microevolution as we earlier revealed in pea crop wild relatives. We sequenced 22 plastid and mitochondrial genomes of Pisum sativum subsp. elatius and Pisum fulvum using Illumina platform, so that the updated sample comprised 64 accessions. Most wild peas from continental southern Europe and a single specimen from Morocco were found to share the same organellar genome constitution; four others, presumably hybrid constitutions, were revealed in Mediterranean islands and Athos Peninsula. A mitochondrial genome closely related to that of Pisum abyssinicum, from Yemen and Ethiopia, was unexpectedly found in an accession of P. sativum subsp. elatius from Israel, their plastid genomes being unrelated. Phylogenetic reconstructions based on plastid and mitochondrial genomes revealed different sets of wild peas to be most related to cultivated P. sativum subsp. sativum, making its wild progenitor and its origin area enigmatic. An accession of P. fulvum representing ‘fulvum-b’ branch, according to a nuclear marker, appeared in the same branch as other fulvum accessions in organellar trees. The results stress the complicated evolution and structure of genetic diversity of pea crop wild relatives.

Comparative analysis of 84 chloroplast genomes of Tylosema esculentum reveals two distinct cytotypes

Tylosema esculentum (marama bean) is an important orphan legume from southern Africa that has long been considered to have the potential to be domesticated as a crop. The chloroplast genomes of 84 marama samples collected from various geographical locations in Namibia and Pretoria were compared in this study. The cp genomes were analyzed for diversity, including SNPs, indels, structural alterations, and heteroplasmy. The marama cp genomes ranged in length from 161,537 bp to 161,580 bp and contained the same sets of genes, including 84 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The genes rpoC2 and rpoB, and the intergenic spacers trnT-trnL and ndhG-ndhI were found to be more diverse than other regions of the marama plastome. 15 haplotypes were found to be divided into two groups, differing at 122 loci and at a 230 bp inversion. One type appears to have greater variability within the major genome present, and variations amongst individuals with this type of chloroplast genome seems to be distributed within specific geographic regions but with very limited sampling for some regions. However, deep sequencing has identified that within most of the individuals, both types of chloroplast genomes are present, albeit one is generally at a very low frequency. The inheritance of this complex of chloroplast genomes appears to be fairly constant, providing a conundrum of how the two genomes co-exist and are propagated through generations. The possible consequences for adaptation to the harsh environment in which T. esculentum survives are considered. The results pave the way for marama variety identification, as well as for understanding the origin and evolution of the bean.

Assembly, annotation, and comparative analysis of Ipomoea chloroplast genomes provide insights into the parasitic characteristics of Cuscuta species

In the Convolvulaceae family, around 1650 species belonging to 60 genera are widely distributed globally, mainly in the tropical and subtropical regions of America and Asia. Although a series of chloroplast genomes in Convolvulaceae were reported and investigated, the evolutionary and genetic relationships among the chloroplast genomes of the Convolvulaceae family have not been extensively elucidated till now. In this study, we first reported the complete chloroplast genome sequence of Ipomoea pes-caprae, a widely distributed coastal plant with medical values. The chloroplast genome of I. pes-caprae is 161667 bp in length, and the GC content is 37.56%. The chloroplastic DNA molecule of I. pes-caprae is a circular structure composed of LSC (large-single-copy), SSC (small-single-copy), and IR (inverted repeat) regions, with the size of the three regions being 88210 bp, 12117 bp, and 30670 bp, respectively. The chloroplast genome of I. pes-caprae contains 141 genes, and 35 SSRs are identified in the chloroplast genome. Our research results provide important genomic information for the molecular phylogeny of I. pes-caprae. The Phylogenetic analysis of 28 Convolvulaceae chloroplast genomes showed that the relationship of I. pes-caprae with I. involucrata or I. obscura was much closer than that with other Convolvulaccae species. Further comparative analyses between the Ipomoea species and Cuscuta species revealed the mechanism underlying the formation of parasitic characteristics of Cuscuta species from the perspective of the chloroplast genome.

CPJSdraw: analysis and visualization of junction sites of chloroplast genomes

Background

Chloroplast genomes are usually circular molecules, and most of them are tetrad structures with two inverted repeat (IR) regions, a large single-copy region, and a small single-copy region. IR contraction and expansion are among the genetic diversities during the evolution of plant chloroplast genomes. The only previously released tool for the visualization of junction sites of the regions does not consider the diversity of the starting point of genomes, which leads to incorrect results or even no results for the examination of IR contraction and expansion.

Results

In this work, a new tool named CPJSdraw was developed for visualizing the junction sites of chloroplast genomes. CPJSdraw can format the starting point of the irregular linearized genome, correct the junction sites of IR and single-copy regions, display the tetrad structure, visualize the junction sites of any number (≥1) of chloroplast genomes, show the transcription direction of genes adjacent to junction sites, and indicate the IR expansion or contraction of chloroplast genomes.

Conclusions

CPJSdraw is a software that is universal and reliable in analysis and visualization of IR expansion or contraction of chloroplast genomes. CPJSdraw has more accurate analysis and more complete functions when compared with previously released tool. CPJSdraw as a perl package and tested data are available at http://dx.doi.org/10.5281/zenodo.7669480 for English users. In addition, an online version with a Chinese interface is available at http://cloud.genepioneer.com:9929/#/tool/alltool/detail/335.

Characterization of the Complete Chloroplast Genome and Phylogenetic Implications of Euonymus microcarpus (Oliv.) Sprague

Euonymus microcarpus (Oliv.) Sprague, is a species of evergreen shrub of the genus Euonymus, family Celastraceae. Here, we extracted the genomic DNA from the leaves of E. microcarpus and constructed a paired-end library. The chloroplast genome of E. microcarpus was generated with the high-throughput sequencing by the illumina Hiseq X Ten platform and de novo assembly. The chloroplast genome had a quadripartite structure, containing a long single copy region with a size of 85,386 bp and a short single copy region with a size of 18,456 bp, separated by two inverted repeat regions of 26,850 bp. The chloroplast genome contained 133 genes identified in total, including 87 potential protein-coding genes, 38 transfer RNA genes, and eight ribosomal RNA genes. A total of 282 simple sequence repeats and 63 long repeats were found. Furthermore, the phylogenetic relationships inferred that E. microcarpus is sister to E. japonicus and E. schensianus. A comparison of the structure of the chloroplast genomes of eight Euonymus species suggests a nucleotide variability of the junction sites and a higher divergence of non-coding regions, compared to the coding regions. The original findings of the study serves as a good reference for chloroplast genome assembly and a valuable foundation for the genetic diversity and evolution of E. microcarpus.

Comparative chloroplast genomics provides insights into the genealogical relationships of endangered Tetraena mongolica and the chloroplast genome evolution of related Zygophyllaceae species

A comprehensive understanding of genetic background for rare species will provide an important theoretical basis for the future species management, monitoring and conservation. Tetraena mongolica is restrictedly distributed in the western Ordos plateau of China and has been listed as a national protected plant. We generated 13 chloroplast (cp) genomes of T. mongolica (size range of 106,062-106,230 bp) and conducted a series of comparative analyses of six Zygophyllaceae cp genomes. T. mongolica cp genome exhibited a quadripartite structure with drastically reduced inverted repeats (IRs, 4,315 bp) and undergone the loss of a suit of ndh genes and a copy of rRNAs. Furthermore, all the T. mongolica populations were divided into two genetic groups based on complete cp phylogenomics. In addition, notably variable genome size, gene order and structural changes had been observed among the six Zygophyllaceae cp genomes. Overall, our findings provide insights into the cp genome evolution mode and intraspecific relationships of T. mongolica, and provide a molecular basis for scientific conservation of this endangered plant.

Comparative and phylogenetic analysis of complete chloroplast genomes from five Artemisia species

Artemisia Linn. is a large genus within the family Asteraceae that includes several important medicinal plants. Because of their similar morphology and chemical composition, traditional identification methods often fail to distinguish them. Therefore, developing an effective identification method for Artemisia species is an urgent requirement. In this study, we analyzed 15 chloroplast (cp) genomes, including 12 newly sequenced genomes, from 5 Artemisia species. The cp genomes from the five Artemisia species had a typical quadripartite structure and were highly conserved across species. They had varying lengths of 151,132-151,178 bp, and their gene content and codon preferences were similar. Mutation hotspot analysis identified four highly variable regions, which can potentially be used as molecular markers to identify Artemisia species. Phylogenetic analysis showed that the five Artemisia species investigated in this study were sister branches to each other, and individuals of each species formed a monophyletic clade. This study shows that the cp genome can provide distinguishing features to help identify closely related Artemisia species and has the potential to serve as a universal super barcode for plant identification.

Assembly and comparative analysis of the complete mitochondrial genome of Salix wilsonii using PacBio HiFi sequencing

Salix L. (willows) is one of the most taxonomically complex genera of flowering plants, including shrubs, tall trees, bushes, and prostrate plants. Despite the high species diversity, only five mitochondrial genomes (mitogenomes) have been released in this genus. Salix wilsonii is an important ornamental and economic willow tree in section Wilsonia of the genus Salix. In this study, the S. wilsonii mitogenome was assembled into a typical circular structure with a size of 711,456 bp using PacBio HiFi sequencing. A total of 58 genes were annotated in the S. wilsonii mitogenome, including 33 protein-coding genes (PCGs), 22 tRNAs, and 3 rRNAs. In the S. wilsonii mitogenome, four genes (mttB, nad3, nad4, and sdh4) were found to play important roles in its evolution through selection pressure analysis. Collinearity analysis of six Salix mitogenomes revealed high structural variability. To determine the evolutionary position of S. wilsonii, we conducted a phylogenetic analysis of the mitogenomes of S. wilsonii and 12 other species in the order Malpighiales. Results strongly supported the segregation of S. wilsonii and other five Salix species with 100% bootstrap support. The comparative analysis of the S. wilsonii mitogenome not only sheds light on the functional and structural features of S. wilsonii but also provides essential information for genetic studies of the genus Salix.

Complete Chloroplast Genome Sequence of the Endemic and Endangered Plant Dendropanax oligodontus: Genome Structure, Comparative and Phylogenetic Analysis

Dendropanax oligodontus, which belongs to the family Araliaceae, is an endemic and endangered species of Hainan Island, China. It has potential economic and medicinal value owing to the presence of phenylpropanoids, flavonoids, triterpenoids, etc. The analysis of the structure and characteristics of the D. oligodontus chloroplast genome (cpDNA) is crucial for understanding the genetic and phylogenetic evolution of this species. In this study, the cpDNA of D. oligodontus was sequenced for the first time using next-generation sequencing methods, assembled, and annotated. We observed a circular quadripartite structure comprising a large single-copy region (86,440 bp), a small single-copy region (18,075 bp), and a pair of inverted repeat regions (25,944 bp). The total length of the cpDNA was 156,403 bp, and the GC% was 37.99%. We found that the D. oligodontus chloroplast genome comprised 131 genes, with 86 protein-coding genes, 8 rRNA genes, and 37 tRNAs. Furthermore, we identified 26,514 codons, 13 repetitive sequences, and 43 simple sequence repeat sites in the D. oligodontus cpDNA. The most common amino acid encoded was leucine, with a strong A/T preference at the third position of the codon. The prediction of RNA editing sites in the protein-coding genes indicated that RNA editing was observed in 19 genes with a total of 54 editing sites, all of which involved C-to-T transitions. Finally, the cpDNA of 11 species of the family Araliaceae were selected for comparative analysis. The sequences of the untranslated regions and coding regions among 11 species were highly conserved, and minor differences were observed in the length of the inverted repeat regions; therefore, the cpDNAs were relatively stable and consistent among these 11 species. The variable hotspots in the genome included clpP, ycf1, rnK-rps16, rps16-trnQ, atpH-atpI, trnE-trnT, psbM-trnD, ycf3-trnS, and rpl32-trnL, providing valuable molecular markers for species authentication and regions for inferring phylogenetic relationships among them, as well as for evolutionary studies. Evolutionary selection pressure analysis indicated that the atpF gene was strongly subjected to positive environmental selection. Phylogenetic analysis indicated that D. oligodontus and Dendropanax dentiger were the most closely related species within the genus, and D. oligodontus was closely related to the genera Kalopanax and Metapanax in the Araliaceae family. Overall, the cp genomes reported in this study will provide resources for studying the genetic diversity and conservation of the endangered plant D. oligodontus, as well as resolving phylogenetic relationships within the family.

Evaluation of genome size and phylogenetic relationships of the Saccharum complex species

"Saccharum complex" is a hypothetical group of species, which is supposed to be involved in the origin of modern sugarcane, and displays large genomes and complex chromosomal alterations. The utilization of restricted parents in breeding programs of modern cultivated sugarcane has resulted in a genetic blockage, which controlled its improvement because of the limited genetic diversity. The use of wild relatives is an effective way to broaden the genetic composition of cultivated sugarcane. Due to the infrequent characterization of genomes, the potential of wild relatives is diffused in improving the cultivated sugarcane. To characterize the genomes of the wild relatives, the genome size and phylogenetic relationships among eight species, including Saccharum spontaneum, Erianthus arundinaceus, E. fulvus, E. rockii, Narenga porphyrocoma, Miscanthus floridulus, Eulalia quadrinervis, and M. sinensis were evaluated based on flow cytometry, genome surveys, K-mer analysis, chloroplast genome sequencing, and whole-genome SNPs analysis. We observed highly heterozygous genomes of S. spontaneum, E. rockii, and E. arundinaceus and the highly repetitive genome of E. fulvus. The genomes of Eulalia quadrinervis, N. porphyrocoma, M. sinensis, and M. floridulus were highly complex. Phylogenetic results of the two approaches were dissimilar, however, both indicate E. fulvus displayed closer relationships to Miscanthus and Saccharum than other species of Saccharum complex. Eulalia quadrinervis was more closely related to M. floridulus than M. sinensis; E. arundinaceus differ significantly from Miscanthus, Narenga, and Saccharum, but was relatively close to Erianthus. We proved the point of E. rockii and E. fulvus should not be classified as one genus, and E. fulvus should be classified as the Saccharum genus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03338-5.

Characterization of the Chloroplast Genome of Argyranthemum frutescens and a Comparison with Other Species in Anthemideae

Argyranthemum frutescens, which belongs to the Anthemideae (Asteraceae), is widely cultivated as an ornamental plant. In this study, the complete chloroplast genome of A. frutescens was obtained based on the sequences generated by Illumina HiSeq. The chloroplast genome of A. frutescens was 149,626 base pairs (bp) in length, containing a pair of inverted repeats (IR, 24,510 bp) regions separated by a small single-copy (SSC, 18,352 bp) sequence and a large single-copy (LSC, 82,254 bp) sequence. The genome contained 132 genes, consisting of 85 coding DNA sequences, 37 tRNA genes, and 8 rRNA genes, with nineteen genes duplicated in the IR region. A comparison chloroplast genome analysis among ten species from the tribe of Anthemideae revealed that the chloroplast genome size varied, but the genome structure, gene content, and oligonucleotide repeats were highly conserved. Highly divergent regions, e.g., ycf1, trnK-psbK, petN-psbM intronic, were detected. Phylogenetic analysis supported Argyranthemum as a separate genus. The findings of this study will be helpful in the exploration of the phylogenetic relationships of the tribe of Anthemideae and contribute to the breeding improvement of A. frutescens.

Analysis of complete chloroplast genome sequences and insight into the phylogenetic relationships of Ferula L

Background

Ferula L. is one of the largest and most taxonomically complicated genera as well as being an important medicinal plant resource in the family Apiaceae. To investigate the plastome features and phylogenetic relationships of Ferula and its neighboring genera Soranthus Ledeb., Schumannia Kuntze., and Talassia Korovin, we sequenced 14 complete plastomes of 12 species. 

Results

The size of the 14 complete chloroplast genomes ranged from 165,607 to 167,013 base pairs (bp) encoding 132 distinct genes (87 protein-coding, 37 tRNA, and 8 rRNA genes), and showed a typical quadripartite structure with a pair of inverted repeats (IR) regions. Based on comparative analysis, we found that the 14 plastomes were similar in codon usage, repeat sequence, simple sequence repeats (SSRs), and IR borders, and had significant collinearity. Based on our phylogenetic analyses, Soranthus, Schumannia, and Talassia should be considered synonymous with Ferula. Six highly divergent regions (rps16/trnQ-UUG, trnS-UGA/psbZ, psbH/petB, ycf1/ndhF, rpl32, and ycf1) were also detected, which may represent potential molecular markers, and combined with selective pressure analysis, the weak positive selection gene ccsA may be a discriminating DNA barcode for Ferula species.

Conclusion

Plastids contain abundant informative sites for resolving phylogenetic relationships. Combined with previous studies, we suggest that there is still much room for improvement in the classification of Ferula. Overall, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of this genus.

Supplementary Information

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

Comparative Analysis of Complete Chloroplast Genomes of Nine Species of Litsea (Lauraceae): Hypervariable Regions, Positive Selection, and Phylogenetic Relationships

Litsea is a group of evergreen trees or shrubs in the laurel family, Lauraceae. Species of the genus are widely used for a wide range of medicinal and industrial aspects. At present, most studies related to the gene resources of Litsea are restricted to morphological analyses or features of individual genomes, and currently available studies of select molecular markers are insufficient. In this study, we assembled and annotated the complete chloroplast genomes of nine species in Litsea, carried out a series of comparative analyses, and reconstructed phylogenetic relationships within the genus. The genome length ranged from 152,051 to 152,747 bp and a total of 128 genes were identified. High consistency patterns of codon bias, repeats, divergent analysis, single nucleotide polymorphisms (SNP) and insertions and deletions (InDels) were discovered across the genus. Variations in gene length and the presence of the pseudogene ycf1^Ψ, resulting from IR contraction and expansion, are reported. The hyper-variable gene rpl16 was identified for its exceptionally high Ka/Ks and Pi values, implying that those frequent mutations occurred as a result of positive selection. Phylogenetic relationships were recovered for the genus based on analyses of full chloroplast genomes and protein-coding genes. Overall, both genome sequences and potential molecular markers provided in this study enrich the available genomic resources for species of Litsea. Valuable genomic resources and divergent analysis are also provided for further research of the evolutionary patterns, molecular markers, and deeper phylogenetic relationships of Litsea.

Application of chloroplast genome in the identification of Traditional Chinese Medicine Viola philippica

BACKGROUND

Viola philippica Cav. is the only source plant of "Zi Hua Di Ding", which is a Traditional Chinese Medicine (TCM) that is utilized as an antifebrile and detoxicant agent for the treatment of acute pyogenic infections. Historically, many Viola species with violet flowers have been misused in "Zi Hua Di Ding". Viola have been recognized as a taxonomically difficult genera due to their highly similar morphological characteristics. Here, all common V. philippica adulterants were sampled. A total of 24 complete chloroplast (cp) genomes were analyzed, among these 5 cp genome sequences were downloaded from GenBank and 19 cp genomes, including 2 "Zi Hua Di Ding" purchased from a local TCM pharmacy, were newly sequenced.

RESULTS

The Viola cp genomes ranged from 156,483 bp to 158,940 bp in length. A total of 110 unique genes were annotated, including 76 protein-coding genes, 30 tRNAs, and four rRNAs. Sequence divergence analysis screening identified 16 highly diverged sequences; these could be used as markers for the identification of Viola species. The morphological, maximum likelihood and Bayesian inference trees of whole cp genome sequences and highly diverged sequences were divided into five monophyletic clades. The species in each of the five clades were identical in their positions within the morphological and cp genome tree. The shared morphological characters belonging to each clade was summarized. Interestingly, unique variable sites were found in ndhF, rpl22, and ycf1 of V. philippica, and these sites can be selected to distinguish V. philippica from samples all other Viola species, including its most closely related species. In addition, important morphological characteristics were proposed to assist the identification of V. philippica. We applied these methods to examine 2 "Zi Hua Di Ding" randomly purchased from the local TCM pharmacy, and this analysis revealed that the morphological and molecular characteristics were valid for the identification of V. philippica.

CONCLUSIONS

This study provides invaluable data for the improvement of species identification and germplasm of V. philippica that may facilitate the application of a super-barcode in TCM identification and enable future studies on phylogenetic evolution and safe medical applications.

Plastid Phylogenomics and Plastomic Diversity of the Extant Lycophytes

Although extant lycophytes represent the most ancient surviving lineage of early vascular plants, their plastomic diversity has long been neglected. The ancient evolutionary history and distinct genetic diversity patterns of the three lycophyte families, each with its own characteristics, provide an ideal opportunity to investigate the interfamilial relationships of lycophytes and their associated patterns of evolution. To compensate for the lack of data on Lycopodiaceae, we sequenced and assembled 14 new plastid genomes (plastomes). Combined with other lycophyte plastomes available online, we reconstructed the phylogenetic relationships of the extant lycophytes based on 93 plastomes. We analyzed, traced, and compared the plastomic diversity and divergence of the three lycophyte families (Isoëtaceae, Lycopodiaceae, and Selaginellaceae) in terms of plastomic diversity by comparing their plastome sizes, GC contents, substitution rates, structural rearrangements, divergence times, ancestral states, RNA editings, and gene losses. Comparative analysis of plastid phylogenomics and plastomic diversity of three lycophyte families will set a foundation for further studies in biology and evolution in lycophytes and therefore in vascular plants.

Comparative Analyses of Six Complete Chloroplast Genomes from the Genus Cupressus and Juniperus (Cupressaceae)

Cupressaceae is a conifer family distributed around the world. Cupressus and Juniperus are the main genera of the Cupressaceae family and have important medicinal value. This leads to confusion between Cupressus and Juniperus due to similar morphologies. Here, the complete cp genomes of two Cupressus (C. duclouxiana and C. funebri) and four Juniperus (J. chinensis, J. gaussenii J. pingii and J. procumbens) were sequenced. The results revealed that the length of the cp genomes ranged from 126,996 bp to 129,959 bp, with 119 genes comprising 82 protein-coding genes, 33 transfer RNAs and 4 ribosomal RNAs. All chloroplast genomes of Cupressus and Juniperus lost whole IR regions, which is consistent with gymnosperm cp genome studies. In addition, the number of SSRs per species ranged from 54 to 73 and was dominated by mononucleotide repeats. In the six cp genomes of Cupressus and Juniperus, five highly divergent regions, including accD, accD-rpl2, ycf1, ycf2 and rrn23-rrn4.5, can be used as DNA barcodes of interspecific relationships and potential genetic markers. We compared the gene selection pressures (C. chengiana as reference species), and 6 genes underwent positive selection, the majority of which were related to photosynthesis. Phylogenetic results showed that the monophyly of Cupressus and Juniperus supported most bootstrap support. Cupressus funebris and J. chinensis were resolved to be early diverging species within Cupressus and Juniperus, and the two genera were sister groups to each other. This research revealed a new understanding of the structural pluralism and phylogenetic relationships of Cupressaceae cp genomes. These results will facilitate comprehension of the complexity and diversity of conifer cp genomes. SIGNIFICANCE:: Phylogenetic relationships among Cupressus, Juniperus, and their closest relatives are controversial, and generic delimitations have been in flux for the past decade. To address relationships and attempt to produce a more robust classification, we sequenced 6 new plastid genomes (plastomes) from the two variously described genera in this complex (Cupressus and Juniperus) and compared them with additional plastomes from diverse members of Cupressaceae. Our study corroborated the accD of Cupressophytes have a tendency to expand in size and strongly supported a sister relationship between Cupressus and Juniperus. The disparity in these results could be traced to the facts that the chloroplast genome is uniparentally inherited, also the usage of the whole chloroplast genome for this research is of a better advantage compared to usage of selected genes or portion of the plastome. The complete CP genomic data will provide useful information for studying genetic diversity and species identification, which is important for the overarching goal of biodiversity conservation.

Characterization of the complete chloroplast genome of Wisteriopsis reticulata (Fabaceae): an IRLC legumes

The inverted repeat‐lacking clade (IRLC) species are characterized by the loss of an IR region in their plastomes, which has long been of great interest. Wisteriopsis reticulata is one of the members of the tribe Wisterieae, which belongs to Fabaceae and is well-known as IRLC. Here, we reported and characterized the complete chloroplast genome of W. reticulata using the genome skimming approach. The chloroplast genome is 132,477 bp in length and lacks one copy of IR region. The genome encoded 112 unique genes including 89 protein-coding genes, 29 transfer RNA genes and four ribosomal RNA genes. Phylogenetic results supported the monophyly of the tribe Wisterieae (IRLC) and confirmed that W. reticulata belongs to the genus Wisteriopsis.

Infrageneric Plastid Genomes of Cotoneaster (Rosaceae): Implications for the Plastome Evolution and Origin of C. wilsonii on Ulleung Island

Cotoneaster is a taxonomically and ornamentally important genus in the family Rosaceae; however, phylogenetic relationships among its species are complicated owing to insufficient morphological diagnostic characteristics and hybridization associated with polyploidy and apomixis. In this study, we sequenced the complete plastomes of seven Cotoneaster species (C. dielsianus, C. hebephyllus, C. integerrimus, C. mongolicus, C. multiflorus, C. submultiflorus, and C. tenuipes) and included the available complete plastomes in a phylogenetic analysis to determine the origin of C. wilsonii, which is endemic to Ulleung Island, Korea. Furthermore, based on 15 representative lineages within the genus, we carried out the first comparative analysis of Cotoneaster plastid genomes to gain an insight into their molecular evolution. The plastomes were highly conserved, with sizes ranging from 159,595 bp (C. tenuipes) to 160,016 bp (C. hebephyllus), and had a GC content of 36.6%. The frequency of codon usage showed similar patterns among the 15 Cotoneaster species, and 24 of the 35 protein-coding genes were predicted to undergo RNA editing. Eight of the 76 common protein-coding genes, including ccsA, matK, ndhD, ndhF, ndhK, petA, rbcL, and rpl16, were positively selected, implying their potential roles in adaptation and speciation. Of the 35 protein-coding genes, 24 genes (15 photosynthesis-related, seven self-replications, and three others) were found to harbor RNA editing sites. Furthermore, several mutation hotspots were identified, including trnG-UCC/trnR-UCU/atpA and trnT-UGU/trnL-UAA. Maximum likelihood analysis based on 57 representative plastomes of Cotoneaster and two Heteromeles plastomes as outgroups revealed two major lineages within the genus, which roughly correspond to two subgenera, Chaenopetalum and Cotoneaster. The Ulleung Island endemic, C. wilsonii, shared its most recent common ancestor with two species, C. schantungensis and C. zabelii, suggesting its potential origin from geographically close members of the subgenus Cotoneaster, section Integerrimi.

Comparative analysis of chloroplast genomes reveals phylogenetic relationships and intraspecific variation in the medicinal plant Isodon rubescens

Isodon rubescens (Hemsley) H. Hara (Lamiaceae) is a traditional Chinese medicine plant that has been used to treat various human diseases and conditions such as inflammation, respiratory and gastrointestinal bacterial infections, and malignant tumors. However, the contents of the main active components of I. rubescens from different origins differ significantly, which greatly affected its quality. Therefore, a molecular method to identify and classify I. rubescens is needed. Here, we report the DNA sequence of the chloroplast genome of I. rubescens collected from Lushan, Henan province. The genome is 152,642 bp in length and has a conserved structure that includes a pair of IR regions (25,726 bp), a LSC region (83,527 bp) and a SSC region (17,663 bp). The chloroplast genome contains 113 unique genes, four rRNA genes, 30 tRNA genes, and 79 protein-coding genes, 23 of which contain introns. The protein-coding genes account for a total of 24,412 codons, and most of them are A/T biased usage. We identified 32 simple sequence repeats (SSRs) and 48 long repeats. Furthermore, we developed valuable chloroplast molecular resources by comparing chloroplast genomes from three Isodon species, and both mVISTA and DnaSP analyses showed that rps16-trnQ, trnS-trnG, and ndhC-trnM are candidate regions that will allow the identification of intraspecific differences within I. rubescens. Also 14 candidate fragments can be used to identify interspecific differences between species in Isodon. A phylogenetic analysis of the complete chloroplast genomes of 24 species in subfamily Nepetoideae was performed using the maximum likelihood method, and shows that I. rubescens clustered closer to I. serra than I. lophanthoides. Interestingly, our analysis showed that I. rubescens (MW018469.1) from Xianyang, Shaanxi Province (IR-X), is closer to I. serra than to the other two I. rubescens accessions. These results strongly indicate that intraspecific diversity is present in I. rubescens. Therefore, our results provide further insight into the phylogenetic relationships and interspecific diversity of species in the genus Isodon.

Nucleotide Evolution, Domestication Selection, and Genetic Relationships of Chloroplast Genomes in the Economically Important Crop Genus Gossypium

Gossypium hirsutum (upland cotton) is one of the most economically important crops worldwide, which has experienced the long terms of evolution and domestication process from wild species to cultivated accessions. However, nucleotide evolution, domestication selection, and the genetic relationship of cotton species remain largely to be studied. In this study, we used chloroplast genome sequences to determine the evolutionary rate, domestication selection, and genetic relationships of 72 cotton genotypes (36 cultivated cotton accessions, seven semi-wild races of G. hirsutum, and 29 wild species). Evolutionary analysis showed that the cultivated tetraploid cotton genotypes clustered into a single clade, which also formed a larger lineage with the semi-wild races. Substitution rate analysis demonstrated that the rates of nucleotide substitution and indel variation were higher for the wild species than the semi-wild and cultivated tetraploid lineages. Selection pressure analysis showed that the wild species might have experienced greater selection pressure, whereas the cultivated cotton genotypes underwent artificial and domestication selection. Population clustering analysis indicated that the cultivated cotton accessions and semi-wild races have existed the obviously genetic differentiation. The nucleotide diversity was higher in the semi-wild races compared with the cultivated genotypes. In addition, genetic introgression and gene flow occurred between the cultivated tetraploid cotton and semi-wild genotypes, but mainly via historical rather than contemporary gene flow. These results provide novel molecular mechanisms insights into the evolution and domestication of economically important crop cotton species.

Comparative and phylogenetic analyses of six Kenya Polystachya (Orchidaceae) species based on the complete chloroplast genome sequences

BACKGROUND

Polystachya Hook. is a large pantropical orchid genus (c. 240 species) distributed in Africa, southern Asia and the Americas, with the center of diversity in Africa. Previous studies on species of this genus have not obtained the complete chloroplast genomes, structures and variations. Additionally, the phylogenetic position of the genus in the Orchidaceae is still controversial and uncertain. Therefore, in this study, we sequenced the complete plastomes of six Kenya Polystachya species based on genome skimming, subjected them to comparative genomic analysis, and reconstructed the phylogenetic relationships with other Orchidaceae species.

RESULTS

The results exhibited that the chloroplast genomes had a typical quadripartite structure with conserved genome arrangement and moderate divergence. The plastomes of the six Polystachya species ranged from 145,484 bp to 149,274 bp in length and had an almost similar GC content of 36.9-37.0%. Gene annotation revealed 106-109 single-copy genes. In addition, 19 genes are duplicated in the inverted regions, and 16 genes each possessd one or more introns. Although no large structural variations were observed among the Polystachya plastomes, about 1 kb inversion was found in Polystachya modesta and all 11 ndh genes in the Polystachya plastomes were lost or pseudogenized. Comparative analysis of the overall sequence identity among six complete chloroplast genomes confirmed that for both coding and non-coding regions in Polystachya, SC regions exhibit higher sequence variation than IRs. Furthermore, there were various amplifications in the IR regions among the six Polystachya species. Most of the protein-coding genes of these species had a high degree of codon preference. We screened out SSRs and found seven relatively highly variable loci. Moreover, 13 genes were discovered with significant positive selection. Phylogenetic analysis showed that the six Polystachya species formed a monophyletic clade and were more closely related to the tribe Vandeae. Phylogenetic relationships of the family Orchidaceae inferred from the 85 chloroplast genome sequences were generally consistent with previous studies and robust.

CONCLUSIONS

Our study is the initial report of the complete chloroplast genomes of the six Polystachya species, elucidates the structural characteristics of the chloroplast genome of Polystachya, and filters out highly variable sequences that can contribute to the development of DNA markers for use in the study of genetic variability and evolutionary studies in Polystachya. In addition, the phylogenetic results strongly support that the genus of Polystachya is a part of the tribe Vandeae.

The complete nucleotide sequence of Viburnum odoratissimum chloroplast genome

The complete chloroplast genome (cpDNA) Viburnum odoratissimumis sequenced and assembled from the whole genome data. The cpDNA of V. odoratissimum is 158,744bp in lengthwith the overall GC content of 38.1%. It consists of a pair of inverted repeats (IRa and IRb, 26,494bp) which separate a large single copy (LCS, 87,348 bp) and small single copy (SSC,18,408 bp). The complete chloroplast genome contains 129 genes, including 84 protein-coding genes, 8 rRNA genes and 37 tRNA genes. A phylogenetic tree of DNA sequences of barcoding regions, including rbcL, matK, psbA-trnHfrom 8 species of the Genus Viburnum shows that V. odoratissimum is closely related to V. furcatum and V. burejaeticum.

The chloroplast genome sequence and phylogenetic analysis of Apocynum venetum L

Apocynum venetum L. (Apocynaceae) is valuable for its medicinal compounds and fiber content. Native A. venetum populations are threatened and require protection. Wild A. venetum resources are limited relative to market demand and a poor understanding of the composition of A. venetum at the molecular level. The chloroplast genome contains genetic markers for phylogenetic analysis, genetic diversity evaluation, and molecular identification. In this study, the entire genome of the A. venetum chloroplast was sequenced and analyzed. The A. venetum cp genome is 150,878 bp, with a pair of inverted repeat regions (IRA and IRB). Each inverted repeat region is 25,810 bp, which consist of large (LSC, 81,951 bp) and small (SSC, 17,307 bp) single copy areas. The genome-wide GC content was 38.35%, LSC made up 36.49%, SSC made up 32.41%, and IR made up 43.3%. The A. venetum chloroplast genome encodes 131 genes, including 86 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. This study identified the unique characteristics of the A. venetum chloroplast genome, which will help formulate effective conservation and management strategies as well as molecular identification approaches for this important medicinal plant.

The complete chloroplast genome of Desmodium styracifolium

The complete chloroplast genome (cpDNA) of Desmodium styracifolium, an important medicinal herb for urolithiasis treatment, was sequenced and assembled from the whole genome data. The cpDNA of D. styracifolium is 149,155 bp in length with GC content of 35.2%. The genome has a quadripartite structure that is composed of a large single-copy (LCS, 82,476 bp) and small single-copy (SSC, 18,439 bp) separated by a pair of inverted repeats (IRa and IRb, 24,120 bp each). There are 128 genes in the chloroplast genome, including 83 protein-coding genes, 8 rRNA genes and 37 tRNA genes.

The complete chloroplast genome of Onobrychis gaubae (Fabaceae-Papilionoideae): comparative analysis with related IR-lacking clade species

BACKGROUND

Plastome (Plastid genome) sequences provide valuable markers for surveying evolutionary relationships and population genetics of plant species. Papilionoideae (papilionoids) has different nucleotide and structural variations in plastomes, which makes it an ideal model for genome evolution studies. Therefore, by sequencing the complete chloroplast genome of Onobrychis gaubae in this study, the characteristics and evolutionary patterns of plastome variations in IR-loss clade were compared.

RESULTS

In the present study, the complete plastid genome of O. gaubae, endemic to Iran, was sequenced using Illumina paired-end sequencing and was compared with previously known genomes of the IRLC species of legumes. The O. gaubae plastid genome was 122,688 bp in length and included a large single-copy (LSC) region of 81,486 bp, a small single-copy (SSC) region of 13,805 bp and one copy of the inverted repeat (IR_b) of 29,100 bp. The genome encoded 110 genes, including 76 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes and possessed 83 simple sequence repeats (SSRs) and 50 repeated structures with the highest proportion in the LSC. Comparative analysis of the chloroplast genomes across IRLC revealed three hotspot genes (ycf1, ycf2, clpP) which could be used as DNA barcode regions. Moreover, seven hypervariable regions [trnL(UAA)-trnT(UGU), trnT(GGU)-trnE(UUC), ycf1, ycf2, ycf4, accD and clpP] were identified within Onobrychis, which could be used to distinguish the Onobrychis species. Phylogenetic analyses revealed that O. gaubae is closely related to Hedysarum. The complete O. gaubae genome is a valuable resource for investigating evolution of Onobrychis species and can be used to identify related species.

CONCLUSIONS

Our results reveal that the plastomes of the IRLC are dynamic molecules and show multiple gene losses and inversions. The identified hypervariable regions could be used as molecular markers for resolving phylogenetic relationships and species identification and also provide new insights into plastome evolution across IRLC.

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.