The First Complete Chloroplast Genome Sequences in Actinidiaceae: Genome Structure and Comparative Analysis

Enhancing Festuca chloroplast genome resources

Festuca L. is a perennial plant of the Poaceae family, known for its good palatability and strong tolerance. The characteristics of the chloroplast genomes of six native Festuca species on the Qinghai-Tibet Plateau were analyzed, including four F. sinensis germplasms from Gansu (1), Qinghai (57), Sichuan (111), and Tibet (10-85); F. rubra subsp. arctica; F. kryloviana; F. nitidula; F. valesiaca subsp. sulcata; and F. elata. The chloroplast genome consists of 127 genes, including 4 ribosomal RNA (rRNA) genes, 38 transfer RNA (tRNA) genes, and 83 protein-coding genes. The number of identified SSRs based on 16 shared optimal ciphers of the four F. sinensis germplasms were 48, 50, 49, and 38, respectively, with the majority being single nucleotide repeat sequences. The phylogenetic tree shows that F. nitidula, F. elata, and one germplasm of F. sinensis from Tibet are grouped into the same clade, while three germplasms of F. sinensis from Gansu, Qinghai, and Sichuan are grouped into another clade. F. valesiaca subsp. sulcata, F. ovina, and F. brevipila are grouped into the same clade. F. rubra subsp. arctica and F. rubra are clustered into one branch. This study clarifies the taxonomy of F. sinensis within the genus, providing a reference for subsequent research on the genetic diversity, genetic structure, and phylogeny of Festuca plants.

The complete chloroplast genome of Stachyurus himalaicus

The complete chloroplast genome of Stachyurus himalaicus Hook. f. et Thoms. ex Benth. was sequenced and analyzed. The circular cp genome of S. himalaicus is 162,642 bp in length, consisting of inverted repeats (IRs; 26,437 bp), small single-copy (SSC; 18,762 bp), and large single-copy regions (LSC; 91,006 bp). The genome contains 131 genes (86 protein-coding, 37 tRNA, 8 rRNA) with 37.0% GC content. Phylogenetic analysis revealed the monophyly of Stachyuraceae, and S. himalaicus was sister to the other 3 Stachyurus species examined. Our results provide useful genetic resources for further studies on the evolution and phylogeny of Stachyuraceae.

The complete chloroplast genome sequence and phylogenetic analysis of Actinidia suberifolia C.Y. Wu (actinidiaceae)

Actinidia suberifolia C.Y. Wu (Actinidiaceae), an endemic species of Yunnan province in China, exhibits substantial therapeutic importance in daily life. Given its narrow distribution and small population size, complete genome sequence is needed to reveal its phylogenetic position in Actinidia. In this study, we assembled and annotated the complete chloroplast genome of A. suberifolia and conducted the phylogenetic analysis among the genus Actinidia. The results showed that A. suberifolia had a typical quadripartite structure, exhibiting a total length of 156,716 bp and consisting of two inverted repeats (IRs) of 23,805 bp separated by a large single-copy (LSC) and a small single-copy (SSC) of 88,437 bp and 20,669 bp. A maximum-likelihood (ML) phylogenetic tree including A. suberifolia and 21 related species indicated that it was close to Actinidia latifolia. The study will offer valuable genetic resources and improve the phylogenetic resolution of Actinidia.

Comparative Chloroplast Genomics of Actinidia deliciosa Cultivars: Insights into Positive Selection and Population Evolution

The chloroplast genome, as an important evolutionary marker, can provide a new breakthrough direction for the population evolution of plant species. Actinidia deliciosa represents one of the most economically significant and widely cultivated fruit species in the genus Actinidia. In this study, we sequenced and analyzed the complete chloroplast genomes of seven cultivars of Actinidia. deliciosa to detect the structural variation and population evolutionary characteristics. The total genome size ranged from 156,404 bp (A. deliciosa cv. Hayward) to 156,495 bp (A. deliciosa cv. Yate). A total of 321 simple sequence repeats (SSRs) and 1335 repetitive sequences were identified. Large-scale repeat sequences may facilitate indels and substitutions, molecular variations in A. deliciosa varieties' chloroplast genomes. Additionally, four polymorphic chloroplast DNA loci (atpF-atpH, atpH-atpI, atpB, and accD) were detected, which could potentially provide useful molecular genetic markers for further population genetics studies within A. deliciosa varieties. Site-specific selection analysis revealed that six genes (atpA, rps3, rps7, rpl22, rbcL, and ycf2) underwent protein sequence evolution. These genes may have played key roles in the adaptation of A. deliciosa to various environments. The population evolutionary analysis suggested that A. deliciosa cultivars were clustered into an individual evolutionary branch with moderate-to-high support values. These results provided a foundational genomic resource that will be a major contribution to future studies of population genetics, adaptive evolution, and genetic improvement in Actinidia.

The complete chloroplast genome sequence of Enkianthus perulatus (Miq.) C.K.Schneid. (Ericaceae)

Enkianthus perulatus (Miq.) C.K.Schneid. is a popular ornamental plant in Japan, for which the chloroplast genome has not been previously reported. Its complete chloroplast genome was therefore determined, using a high-throughput sequencing technique. The total length of the E. perulatus chloroplast genome was 157,256 bp, comprising a large single copy region of 89,680 bp and a small single copy region of 20,472 bp, separated by a pair of 23,552 bp inverted repeat regions. A phylogenetic analysis, based on the chloroplast genome, indicated that E. perulatus is the earliest diverged species within Ericaceae. The chloroplast genome structure of E. perulatus is similar to species from families within Ericales, but not other species of Ericaceae. The presented chloroplast genome will be useful for further systematic analyses and planning conservation strategies.

Comparative genomics and phylogenetic analysis of six Malvaceae species based on chloroplast genomes

The Malvaceae family, comprising 9 subfamilies and 4,225 species, includes economically significant taxa, such as Ceiba pentandra, Gossypium ekmanianum, Gossypium stephensii, Kokia drynarioides, Talipariti hamabo, and Durio zibethinus. Chloroplast (cp) genome research is crucial for elucidating the evolutionary divergence and species identification within this family. In this study, we assembled and annotated cp genomes of six Malvaceae species, conducting comprehensive comparative genomic and phylogenomic analyses. The assembled genomes range from 160,495 to 163,970 bp in size, with 125-129 genes annotated. Notable differences were observed in the IR (inverted repeat) regions, and SSR analysis revealed that Durio zibethinus has the highest number of specific variation sites. Among the six species, Talipariti hamabo uniquely exhibits more palindromic repeats than forward repeats. Seven highly mutated regions were identified, offering potential markers for species identification. Phylogenetic reconstruction using the maximum likelihood method revealed two primary clades within Malvaceae: Byttneriina and Malvadendrina. Within Malvadendrina, the subfamily Helicteroideae represents the earliest divergence, followed by Sterculioideae. This study provides a robust phylogenetic framework and valuable insights into the classification and evolutionary history of Malvaceae species.

Comparative chloroplast genome analysis of five widespread species (Zanthoxylum L.) and development of molecular markers for their discrimination

Background

Zanthoxylum L., an important genus in the Rutaceae family, has great edible and medical values. However, the high degree of morphological similarity among Zanthoxylum species and the lack of sufficient chloroplast (cp) genomic resources have greatly impeded germplasm identification and phylogenetic analyses of Zanthoxylum.

Methods

Here we assembled cp genomes of five widespread species (Zanthoxylum bungeanum, Z. armatum, Z. nitidum, Z. ailanthoides and Z. piasezkii) in China as a case study, comparative analysis of these assembled cp genomes.

Results

Each of them, ranging from 157,231 to 158,728 bp, has a quadripartite structure. Except for one extra gene in Z. piasezkii, 132 genes were identified in each species, including 87 encode protein genes, 37 transfer ribose nucleic acid (tRNA) genes, and eight ribosomal RNA (rRNA) genes. Substantial variation was observed among these five cp genome sequences in the IR/SC boundary regions. Variation in insertions and deletions were observed in the cp genomes of the five species over three different intervals, and a large number of single-nucleotide polymorphism variants were detected in the rps3-rpl22-rps19 region. Phylogenetic analysis of complete cp genome sequences revealed the evolutionary relationships among 23 Zanthoxylum species (29 samples).

Conclusion

Comparative analysis revealed that rps3-rpl22-rps19 is a highly variable divergent region in Zanthoxylum that could be developed as candidate markers for phylogenetic studies and species identification. This study identified a pair of molecular markers from hypervariable regions that can be used to distinguish between the five Zanthoxylum species and validated their utility. Overall, the results of this study provide new insights into the genetic breeding, germplasm exploration, and phylogeny of Zanthoxylum species.

Comparative analysis of the complete chloroplast genome of Pueraria provides insights for species identification, phylogenetic relationships, and taxonomy

BACKGROUND

Pueraria is an edible and medicinal raw material, which is of great value to the pharmaceutical and food industries. Nonetheless, due to morphological diversity and complex domestication history, the classification of Pueraria plants is ambiguous. As the varieties on the market are mixed, the species are difficult to distinguish, and their morphological characteristics are similar to the physical and chemical properties. It is difficult to accurately identify them by traditional identification methods. Chloroplast (cp) genomes are widely used in species identification and phylogenetic studies to achieve accurate identification of medicinal plants, and can also provide more reference information for phylogenetic studies. Based on interspecific and intraspecific sampling, the cp genomes of eight species or varieties of Pueraria plants were examined in this study.

RESULTS

The study unveiled that the cp genome size varied from 151,555 to 153,668 base pairs (bp), with the total GC content ranging from 35.4 to 37.0%. Moreover, it was discerned that the cp genome contained between 128 and 135 genes. Comparative analysis indicated that the highest number of Simple Sequence Repeats (SSRs) was identified in P. montana and P. alopecuroides, with a preponderance of these SSRs being rich in Adenine (A) and Thymine (T) nucleotides. Complete comparison and sliding window analysis of the cp genome established that the non-coding region exhibited greater sequence differences than the coding region, and that the large single copy (LSC) region demonstrated higher nucleotide polymorphism levels. Fourteen highly variable loci such as rpoB,ycf1,rbcL,trnF-GAA-trnL,psbC-psbD, and ycf4-cemA were detected as potential molecular markers for Pueraria species identification. Moreover, the phylogenetic tree demonstrated that other Pueraria species had the most distant relationship with Haymondia wallichii and Toxicopueraria peduncularis, thereby offering fresh perspectives into the species classification of Pueraria. The molecular clock analysis results indicate that the divergence time of Pueraria may occur at ∼6.46 Ma. It is speculated that the cold climate may be the cause of Pueraria species diversity and promote the radiation of the genus.

CONCLUSION

This research provides theoretical backing and serves as a reference point for the identification and taxonomical classification of Pueraria species. The findings will prove beneficial in future studies on the preservation of medicinal resources, phylogenetic relationships, and genetic engineering of Pueraria plants.

New insights on the phylogeny, evolutionary history, and ecological adaptation mechanism in cycle-cup oaks based on chloroplast genomes

Cycle-cup oaks (Quercus section Cyclobalanopsis) are one of the principal components of forests in the tropical and subtropical climates of East and Southeast Asia. They have experienced relatively recent increases in the diversification rate, driven by changing climates and the Himalayan orogeny. However, the evolutionary history and adaptive mechanisms at the chloroplast genome level in cycle-cup oaks remain largely unknown. Therefore, we studied this problem by conducting chloroplast genomics on 50 of the ca. 90 species. Comparative genomics and other analyses showed that Quercus section Cyclobalanopsis had a highly conserved chloroplast genome structure. Highly divergent regions, such as the ndhF and ycf1 gene regions and the petN-psbM and rpoB-trnC-GCA intergenic spacer regions, provided potential molecular markers for subsequent analysis. The chloroplast phylogenomic tree indicated that Quercus section Cyclobalanopsis was not monophyletic, which mixed with the other two sections of subgenus Cerris. The reconstruction of ancestral aera inferred that Palaeotropics was the most likely ancestral range of Quercus section Cyclobalanopsis, and then dispersed to Sino-Japan and Sino-Himalaya. Positive selection analysis showed that the photosystem genes had the lowest ω values among the seven functional gene groups. And nine protein-coding genes containing sites for positive selection: ndhA, ndhD, ndhF, ndhH, rbcL, rpl32, accD, ycf1, and ycf2. This series of analyses together revealed the phylogeny, evolutionary history, and ecological adaptation mechanism of the chloroplast genome of Quercus section Cyclobalanopsis in the long river of earth history. These chloroplast genome data provide valuable information for deep insights into phylogenetic relationships and intraspecific diversity in Quercus.

Phylogenetic relationships, selective pressure and molecular markers development of six species in subfamily Polygonoideae based on complete chloroplast genomes

The subfamily Polygonoideae encompasses a diverse array of medicinal and horticultural plants that hold significant economic value. However, due to the lack of a robust taxonomy based on phylogenetic relationships, the classification within this family is perplexing, and there is also a scarcity of reports on the chloroplast genomes of many plants falling under this classification. In this study, we conducted a comprehensive analysis by sequencing and characterizing the complete chloroplast genomes of six Polygonoideae plants, namely Pteroxygonum denticulatum, Pleuropterus multiflorus, Pleuropterus ciliinervis, Fallopia aubertii, Fallopia dentatoalata, and Fallopia convolvulus. Our findings revealed that these six plants possess chloroplast genomes with a typical quadripartite structure, averaging 162,931 bp in length. Comparative chloroplast analysis, codon usage analysis, and repetitive sequence analysis demonstrated a high level of conservation within the chloroplast genomes of these plants. Furthermore, phylogenetic analysis unveiled a distinct clade occupied by P. denticulatum, while P. ciliinrvis displayed a closer relationship to the three plants belonging to the Fallopia genus. Selective pressure analysis based on maximum likelihood trees showed that a total of 14 protein-coding genes exhibited positive selection, with psbB and ycf1 having the highest number of positive amino acid sites. Additionally, we identified four molecular markers, namely petN-psbM, psal-ycf4, ycf3-trnS-GGA, and trnL-UAG-ccsA, which exhibit high variability and can be utilized for the identification of these six plants.

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.

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome

BACKGROUND

Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear.

RESULTS

We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species.

CONCLUSIONS

These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.

Chloroplast genome of Calamus tetradactylus revealed rattan phylogeny

BACKGROUND

Calamus tetradactylus, a species primarily distributed in Vietnam, Laos, and southern China, is highly valued for its utilization as a small-diameter rattan material. While its physical and mechanical properties have been extensively studied, the genomic characteristics of C. tetradactylus remain largely unexplored.

RESULTS

To gain a better understanding of its chloroplast genomic features and evolutionary relationships, we conducted sequencing and assembly of the chloroplast genome of C. tetradactylus. The complete chloroplast genome exhibited the typical highly conserved quartile structure, with specific variable regions identified in the single-copy region (like psbF-psbE, π = 0.10327, ndhF-rpl32, π = 0.10195), as well as genes such as trnT-GGU (π = 0.05764) and ycf1 (π = 0.03345) and others. We propose that these regions and genes hold potential as markers for species identification. Furthermore, phylogenetic analysis revealed that C. tetradactylus formed a distinct clade within the phylogenetic tree, alongside other Calamus species, and C. tetradactylus was most closely related to C. walkeri, providing support for the monophyly of the genus.

CONCLUSION

The analysis of the chloroplast genome conducted in this study provides valuable insights that can contribute to the improvement of rattan breeding programs and facilitate sustainable development in the future.

The complete plastome sequences of invasive weed Parthenium hysterophorus: genome organization, evolutionary significance, structural features, and comparative analysis

Parthenium hysterophorus, a globally widespread weed, poses a significant threat to agricultural ecosystems due to its invasive nature. We investigated the chloroplast genome of P. hysterophorus in this study. Our analysis revealed that the chloroplast genome of P. hysterophorus spans a length of 151,881 base pairs (bp). It exhibits typical quadripartite structure commonly found in chloroplast genomes, including inverted repeat regions (IR) of 25,085 bp, a small single copy (SSC) region of 18,052 bp, and a large single copy (LSC) region of 83,588 bp. A total of 129 unique genes were identified in P. hysterophorus chloroplast genomes, including 85 protein-coding genes, 36 tRNAs, and eight rRNAs genes. Comparative analysis of the P. hysterophorus plastome with those of related species from the tribe Heliantheae revealed both conserved structures and intriguing variations. While many structural elements were shared among the species, we identified a rearrangement in the large single-copy region of P. hysterophorus. Moreover, our study highlighted notable gene divergence in several specific genes, namely matK, ndhF, clpP, rps16, ndhA, rps3, and ndhD. Phylogenetic analysis based on the 72 shared genes placed P. hysterophorus in a distinct clade alongside another species, P. argentatum. Additionally, the estimated divergence time between the Parthenium genus and Helianthus (sunflowers) was approximately 15.1 million years ago (Mya). These findings provide valuable insights into the evolutionary history and genetic relationships of P. hysterophorus, shedding light on its divergence and adaptation over time.

Comprehensive analysis of complete chloroplast genome sequence of Plantago asiatica L. (Plantaginaceae)

Plantago asiatica L. is a representative individual species of Plantaginaceae, whose high reputation is owed to its edible and medicinal values. However, the phylogeny and genes of the P. asiatica chloroplast have not yet been well described. Here we report the findings of a comprehensive analysis of the P. asiatica chloroplast genome. The P. asiatica chloroplast genome is 164,992 bp, circular, and has a GC content of 37.98%. The circular genome contains 141 genes, including 8 rRNAs, 38 tRNAs, and 95 protein-coding genes. Seventy-two simple sequence repeats are detected. Comparative chloroplast genome analysis of six related species suggests that a higher similarity exists in the coding region than the non-coding region, and differences in the degree of preservation is smaller between P. asiatica and Plantago depressa than among others. Our phylogenetic analysis illustrates P. asiatica has a relatively close relationship with P. depressa, which was also divided into different clades with Plantago ovata and Plantago lagopus in the genus Plantago. This analysis of the P. asiatica chloroplast genome contributes to an improved deeply understanding of the evolutionary relationships among Plantaginaceae.

Chloroplast Genome Comparison and Phylogenetic Analysis of the Commercial Variety Actinidia chinensis 'Hongyang'

Actinidia chinensis 'Hongyang', also known as red yangtao (red heart kiwifruit), is a vine fruit tree native to China possessing significant nutritional and economic value. However, information on its genetic diversity and phylogeny is still very limited. The first chloroplast (cp) genome of A. chinensis 'Hongyang' cultivated in China was sequenced using de novo technology in this study. A. chinensis 'Hongyang' possesses a cp genome that spans 156,267 base pairs (bp), exhibiting an overall GC content of 37.20%. There were 132 genes that were annotated, with 85 of them being protein-coding genes, 39 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. A total of 49 microsatellite sequences (SSRs) were detected, mainly single nucleotide repeats, mostly consisting of A or T base repeats. Compared with 14 other species, the cp genomes of A. chinensis 'Hongyang' were biased towards the use of codons containing A/U, and the non-protein coding regions in the A. chinensis 'Hongyang' cpDNA showed greater variation than the coding regions. The nucleotide polymorphism analysis (Pi) yielded nine highly variable region hotspots, most in the large single copy (LSC) region. The cp genome boundary analysis revealed a conservative order of gene arrangement in the inverted repeats (IRs) region of the cp genomes of 15 Actinidia plants, with small expansions and contractions of the boundaries. Furthermore, phylogenetic tree indicated that A. chinensis 'Hongyang' was the closest relative to A. indochinensis. This research provides a useful basis for future genetic and evolutionary studies of A. chinensis 'Hongyang', and enriches the biological information of Actinidia species.

Comparative Analysis of Plastomes in Elsholtzieae: Phylogenetic Relationships and Potential Molecular Markers

The Elsholtzieae, comprising ca. 7 genera and 70 species, is a small tribe of Lamiaceae (mint family). Members of Elsholtzieae are of high medicinal, aromatic, culinary, and ornamentals value. Despite the rich diversity and value of Elsholtzieae, few molecular markers or plastomes are available for phylogenetics. In the present study, we employed high-throughput sequencing to assemble two Mosla plastomes, M. dianthera and M. scabra, for the first time, and compared with other plastomes of Elsholtzieae. The plastomes of Elsholtzieae exhibited a quadripartite structure, ranging in size from 148,288 bp to 152,602 bp. Excepting the absence of the pseudogene rps19 in Elsholtzia densa, the exhaustive tally revealed the presence of 132 genes (113 unique genes). Among these, 85 protein-coding genes (CDS), 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (rps19 and ycf1) were annotated. Comparative analyses showed that the plastomes of these species have minor variations at the gene level. Notably, the E. eriostchya plastid genome exhibited increased GC content regions in the LSC and SSC, resulting in an increased overall GC content of the entire plastid genome. The E. densa plastid genome displayed modified boundaries due to inverted repeat (IR) contraction. The sequences of CDS and intergenic regions (IGS) with elevated variability were identified as potential molecular markers for taxonomic inquiries within Elsholtzieae. Phylogenetic analysis indicated that four genera formed monophyletic entities, with Mosla and Perilla forming a sister clade. This clade was, in turn, sister to Collinsonia, collectively forming a sister group to Elsholtzia. Both CDS, and CDS + IGS could construct a phylogenetic tree with stronger support. These findings facilitate species identification and DNA barcoding investigations in Elsholtzieae and provide a foundation for further exploration and resource utilization within this tribe.

Complete chloroplast genome structural characterization of two Phalaenopsis (Orchidaceae) species and comparative analysis with their alliance

BACKGROUND

The taxonomy and infrageneric delimitation of Phalaenopsis Blume has been significantly disputed due to some overlapping morphological features between species related, which needed further evidence for clarification. The structural characterization of complete chloroplast genomes of P. storbatiana and P. wilsonii were analyzed and compared with those of related taxa to provide a better understanding of their genomic information on taxonomy and phylogeny.

RESULTS

It was shown that chloroplast genomes of Phalaenopsis storbatiana and P. wilsonii had a typical quadripartite structure with conserved genome arrangements and moderate divergence. The chloroplast genomes of P. storbatiana and P. wilsonii were 145,885 bp and 145,445 bp in length, respectively, and shared a similar GC content of 36.8%. Gene annotations of two species revealed 109 single-copy genes consistently. In addition, 20 genes duplicated in the inverted regions, 16 genes each possessed one or more introns, and five ndh (NA (D)H dehydrogenase) genes were observed in both. Comparative analysis of the total cp genomes of P. storbatiana and P. wilsonii with those of other six related Phalaenopsis species confirmed the stable sequence identity for coding and non-coding regions and higher sequence variation in SC regions than IR regions. Most of their protein-coding genes had a high degree of codon preference. Moreover, 45 genes were discovered with significantly positive selection. However, different amplifications in IR regions were observed in these eight species. Phylogenetic analysis based on CDS from 60 species representing main clades in Orchidaceae indicated that Phalaenopsis species including P. stobartiana and P. wilsonii formed a monophyletic clade with high bootstrap nested in tribe Vandeae of Epidendroideae, which was consistent with those from previous studies.

CONCLUSIONS

The results could provide insight into understanding the plastome evolution and phylogenetic relationships of Phalaenopsis.

The complete chloroplast genome sequence of Eugenia klotzschiana O. Berg unveils the evolutionary dynamics in plastomes of Myrteae DC. tribe (Myrtaceae)

Myrteae is the most diversified tribe in the Myrtaceae family and has great ecological and economic importance. Here, we performed the assembly and annotation of the chloroplast genome of Eugenia klotzschiana O. Berg and used this in a comparative analysis with other 13 species from the Myrteae tribe. The E. klotzschiana plastome exhibited a length of 158,977 bp and a very conserved structure and gene composition when compared with other Myrteae genomes. We identified 34 large repetitive sequences and 94 SSR repeats in E. klotzschiana plastome. The trnT-trnL, rpl32-trnL, ndhF-rpl32, psbE-petL, and ycf1 regions were identified as mutational hotspots. A negative selection signal was detected in 74 protein-coding genes while neutral selection was detected in two genes (rps12 and psaI). Furthermore, 222 RNA editing sites were identified in the E. klotzschiana plastome. We also obtained a plastome-based Myrtales phylogenetic tree, including E. klotzschiana for the first time in a molecular phylogeny, recovering its sister relationship for all other Eugenia species. Our results illuminate how evolution shaped the chloroplast genome structure and composition in the Myrteae tribe, especially in the E. klotzschiana plastome.

The First Complete Chloroplast Genome of Cordia monoica: Structure and Comparative Analysis

Cordia monoica is a member of the Boraginaceae family. This plant is widely distributed in tropical regions and has a great deal of medical value as well as economic importance. In the current study, the complete chloroplast (cp) genome of C. monoica was sequenced, assembled, annotated, and reported. This circular chloroplast genome had a size of 148,711 bp, with a quadripartite structure alternating between a pair of repeated inverted regions (26,897-26,901 bp) and a single copy region (77,893 bp). Among the 134 genes encoded by the cp genome, there were 89 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. A total of 1387 tandem repeats were detected, with the hexanucleotides class making up 28 percent of the repeats. Cordia monoica has 26,303 codons in its protein-coding regions, and leucine amino acid was the most frequently encoded amino acid in contrast to cysteine. In addition, 12 of the 89 protein-coding genes were found to be under positive selection. The phyloplastomic taxonomical clustering of the Boraginaceae species provides further evidence that chloroplast genome data are reliable not only at family level but also in deciphering the phylogeny at genus level (e.g., Cordia).

Complete chloroplast genome and phylogenetic relationship of Nymphaea nouchali (Nymphaeaceae), a rare species of water lily in China

Nymphaea nouchali is a native species of Chinese water lily with important ornamental, economical, and medicinal purposes. However, due to the serious disturbance by alien biological invasion and human factors, N. nouchali is in an endangered state in China and urgently needs to be protected. Here, we reported the complete chloroplast genome of N. nouchali for the first time, and we found that its plastome is 159 978 bp long, comprising large and small single copies and two inverted repeats (90 001, 19 603, and 50 374 bp, respectively), indicating a typical tetrad structure. In total, 130 genes were identified, including 85 protein-coding genes, 37 transfer RNAs, and 8 ribosomal RNAs. Additionally, 136 simple sequence repeat sites were identified, composed mainly of single nucleotide (46.32%) and trinucleotide (47.05%) sequences. Five highly variable sites (psaI, rps19, ndhF, rps15, and ycf1) with a high Pi value were identified as potential molecular markers. Phylogenetic analysis showed that N. nouchali and N. ampla are closely related, and further validated previous water lily classification results based on morphological characteristics, which divided water lilies into five subgenera: Nymphaea, Brachyceras, Anecphya, Hydrocallis, and Lotos. These results are valuable for the identification and the formulation of protection strategies of N. nouchali, as well as contributing to understanding the evolutionary relationships among Nymphaeaceae species.

Comparative chloroplast genome analyses of diverse Phoebe (Lauraceae) species endemic to China provide insight into their phylogeographical origin

The genus Phoebe (Lauraceae) includes about 90 evergreen tree species that are an ideal source of timber. Habitat destruction and deforestation have resulted in most of them being endemic to China. The accurate identification of endangered Phoebe species in China is necessary for their conservation. Chloroplast genome sequences can play an important role in species identification. In this study, comparative chloroplast genome analyses were conducted on diverse Phoebe species that are primarily distributed in China. Despite the conserved nature of chloroplast genomes, we detected some highly divergent intergenic regions (petA-psbE, ndhF-rpl32, and psbM-trnD-GUC) as well as three highly divergent genes (rbcL, ycf1, and ycf2) that have potential applications in phylogenetics and evolutionary analysis. The phylogenetic analysis indicated that various Phoebe species in China were divided into three clades. The complete chloroplast genome was better suited for phylogenetic analysis of Phoebe species. In addition, based on the phylogeographical analysis of Phoebe species in China, we inferred that the Phoebe species in China first originated in Yunnan and then spread to other southern areas of the Yangtze River. The results of this research will add to existing case studies on the phylogenetic analysis of Phoebe species and have the potential to contribute to the conservation of Phoebe species that are in danger of extinction.

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.

Evolutionary Comparison of the Complete Chloroplast Genomes in Convallaria Species and Phylogenetic Study of Asparagaceae

The genus Convallaria (Asparagaceae) comprises three herbaceous perennial species that are widely distributed in the understory of temperate deciduous forests in the Northern Hemisphere. Although Convallaria species have high medicinal and horticultural values, studies related to the phylogenetic analysis of this genus are few. In the present study, we assembled and reported five complete chloroplast (cp) sequences of three Convallaria species (two of C. keiskei Miq., two of C. majalis L., and one of C. montana Raf.) using Illumina paired-end sequencing data. The cp genomes were highly similar in overall size (161,365–162,972 bp), and all consisted of a pair of inverted repeats (IR) regions (29,140–29,486 bp) separated by a large single-copy (LSC) (85,183–85,521 bp) and a small single-copy (SSC) region (17,877–18,502 bp). Each cp genome contained the same 113 unique genes, including 78 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. Gene content, gene order, AT content and IR/SC boundary structure were nearly identical among all of the Convallaria cp genomes. However, their lengths varied due to contraction/expansion at the IR/LSC borders. Simple sequence repeat (SSR) analyses indicated that the richest SSRs are A/T mononucleotides. Three highly variable regions (petA-psbJ, psbI-trnS and ccsA-ndhD) were identified as valuable molecular markers. Phylogenetic analysis of the family Asparagaceae using 48 cp genome sequences supported the monophyly of Convallaria, which formed a sister clade to the genus Rohdea. Our study provides a robust phylogeny of the Asparagaceae family. The complete cp genome sequences will contribute to further studies in the molecular identification, genetic diversity, and phylogeny of Convallaria.

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.

Complete Chloroplast Genomes of 14 Subspecies of D. glomerata: Phylogenetic and Comparative Genomic Analyses

Orchardgrass (Dactylis glomerata L.) is a species in the Gramineae family that is highly important economically and valued for its role in ecology. However, the phylogeny and taxonomy of D. glomerata are still controversial based on current morphological and molecular evidence. The study of chloroplast (cp) genomes has developed into a powerful tool to develop molecular markers for related species and reveal the relationships between plant evolution and phylogenetics. In this study, we conducted comparative genomic analyses and phylogenetic inferences on 14 cp genomes of D. glomerata originating from the Mediterranean and Eurasia. The genome size ranged from 134,375 bp to 134,993 bp and exhibited synteny of gene organization and order. A total of 129–131 genes were identified, including 85–87 protein coding genes, 38 tRNA genes and 8 rRNA genes. The cp sequences were highly conserved, and key sequence variations were detected at the junctions of inverted repeats (IRs)/small single–copy (SSC) regions. Moreover, nine highly variable regions were identified among the subspecies based on a sequence divergence analysis. A total of 285 RNA editing sites were detected that were relevant to 52 genes, where rpoB exhibited the most abundant RNA editing sites. The phylogenetic analysis revealed that all Dactylis subspecies clustered into a monophyletic group and most branches provided a high support bootstrap. The main divergence time of D. glomerata was dated to the Miocene era, and this could have been due to changes in the climate. These findings will provide useful insights for further studies on phylogeny, the identification of subspecies and the development of hypotheses for the evolutionary history of the genus Dactylis and of the Gramineae family.

Taxonomy, comparative genomics of Mullein (Verbascum, Scrophulariaceae), with implications for the evolution of Verbascum and Lamiales

BACKGROUND

The genus Verbascum L. (Scrophulariaceae) is distributed in Africa, Europe, and parts of Asia, with the Mediterranean having the most species variety. Several researchers have already worked on the phylogenetic and taxonomic analysis of Verbascum by using ITS data and chloroplast genome fragments and have produced different conclusions. The taxonomy and phylogenetic relationships of this genus are unclear.

RESULTS

The complete plastomes (cp) lengths for V. chaixii, V. songaricum, V. phoeniceum, V. blattaria, V. sinaiticum, V. thapsus, and V. brevipedicellatum ranged from 153,014 to 153,481 bp. The cp coded 114 unique genes comprising of 80 protein-coding genes, four ribosomal RNA (rRNA), and 30 tRNA genes. We detected variations in the repeat structures, gene expansion on the inverted repeat, and single copy (IR/SC) boundary regions. The substitution rate analysis indicated that some genes were under purifying selection pressure. Phylogenetic analysis supported the sister relationship of (Lentibulariaceae + Acanthaceae + Bignoniaceae + Verbenaceae + Pedaliaceae) and (Lamiaceae + Phyrymaceae + Orobanchaceae + Paulowniaceae + Mazaceae) in Lamiales. Within Scrophulariaceae, Verbascum was sister to Scrophularia, while Buddleja formed a monophyletic clade from (Scrophularia + Verbascum) with high bootstrap support values. The relationship of the nine species within Verbascum was highly supported.

CONCLUSION

Based on the phylogenetic results, we proposed to reinstate the species status of V. brevipedicellatum (Engl.) Hub.-Mor. Additionally, three genera (Mazus, Lancea, and Dodartia) placed in the Phyrymaceae family formed a separate clade within Lamiaceae. The classification of the three genera was supported by previous studies. Thus, the current study also suggests the circumscription of these genera as documented previously to be reinstated. The divergence time of Lamiales was approximated to be 86.28 million years ago (Ma) (95% highest posterior density (HPD), 85.12-89.91 Ma). The complete plastomes sequence data of the Verbascum species will be important for understanding the Verbascum phylogenetic relationships and evolution in order Lamiales.

A Comparative Analysis of the Chloroplast Genomes of Four Polygonum Medicinal Plants

Polygonum is a generalized genus of the Polygonaceae family that includes various herbaceous plants. In order to provide aid in understanding the evolutionary and phylogenetic relationship in Polygonum at the chloroplast (cp) genome-scale level, we sequenced and annotated the complete chloroplast genomes of four Polygonum species using next-generation sequencing technology and CpGAVAS. Then, repeat sequences, IR contractions, and expansion and transformation sites of chloroplast genomes of four Polygonum species were studied, and a phylogenetic tree was built using the chloroplast genomes of Polygonum. The results indicated that the chloroplast genome construction of Polygonum also displayed characteristic four types of results, comparable to the published chloroplast genome of recorded angiosperms. The chloroplast genomes of the four Polygonum plants are highly consistent in genome size (159,015 bp-163,461 bp), number of genes (112 genes, including 78 protein-coding genes, 30 tRNA genes, and 4 rRNA genes), gene types, gene order, codon usage, and repeat sequence distribution, which identifies the high preservation among the Polygonum chloroplast genomes. The Polygonum phylogenetic tree was recreated by a full sequence of the chloroplast genome, which illustrates that the P. bistorta, P. orientale, and P. perfoliatum are divided into the same branch, and P. aviculare belongs to Fallopia. The precise system site of lots base parts requires further verification, but the study would provide a basis for developing the available genetic resources and evolutionary relationships of Polygonum.

Comprehensive analysis of plastid gene expression during fruit development and ripening of kiwifruit

Global survey of plastid gene expression during fruit ripening in kiwifruit provides cis-elements for the future engineering of the plastid genome of kiwifruit. A limitation in the application of plastid biotechnology for molecular farming is the low-level expression of transgenes in non-green plastids compared with photosynthetically active chloroplasts. Unlike other fruits, not all chloroplasts are transformed into chromoplasts during ripening of red-fleshed kiwifruit (Actinidia chinensis cv. Hongyang) fruits, which may make kiwifruit an ideal horticultural plant for recombinant protein production by plastid engineering. To identify cis-elements potentially triggering high-level transgene expression in edible tissues of the 'Hongyang' kiwifruit, here we report a comprehensive analysis of kiwifruit plastid gene transcription in green leaves and fruits at three different developmental stages. While transcripts of a few photosynthesis-related genes and most genetic system genes were substantially upregulated in green fruits compared with leaves, nearly all plastid genes were significantly downregulated at the RNA level during fruit development. Expression of a few genes remained unchanged, including psbA, the gene encoding the D1 polypeptide of photosystem II. However, PsbA protein accumulation decreased continuously during chloroplast-to-chromoplast differentiation. Analysis of post-transcriptional steps in mRNA maturation, including intron splicing and RNA editing, revealed that splicing and editing may contribute to regulation of plastid gene expression. Altogether, 40 RNA editing sites were verified, and 5 of them were newly discovered. Taken together, this study has generated a valuable resource for the analysis of plastid gene expression and provides cis-elements for future efforts to engineer the plastid genome of kiwifruit.

Comparative Analysis of Codon Bias in the Chloroplast Genomes of Theaceae Species

Theaceae species are dicotyledonous angiosperms with extremely high ornamental and economic value. The chloroplast genome is traditionally used to study species evolution, expression of chloroplast genes and chloroplast transformation. Codon usage bias (CUB) analysis is beneficial for investigations of evolutionary relationships and can be used to improve gene expression efficiency in genetic transformation research. However, there are relatively few systematic studies of the CUB in the chloroplast genomes of Theaceae species. In this study, CUB and nucleotide compositions parameters were determined by the scripts written in the Perl language, CodonW 1.4.2, CU.Win2000, RStudio and SPSS 23.0. The chloroplast genome data of 40 Theaceae species were obtained to analyse the codon usage (CU) characteristics of the coding regions and the influence of the source of variation on CUB. To explore the relationship between the CUB and gene expression levels in these 40 Theaceae plastomes, the synonymous codon usage order (SCUO) and measure independent of length and composition (MILC) values were determined. Finally, phylogenetic analysis revealed the genetic evolutionary relationships among these Theaceae species. Our results showed that based on the chloroplast genomes of these 40 Theaceae species, the CUB was for codons containing A/T bases and those that ended with A/T bases. Moreover, there was great commonality in the CUB of the Theaceae species according to comparative analysis of relative synonymous codon usage (RSCU) and relative frequency of synonymous codon (RFSC): these species had 29 identical codons with bias (RSCU > 1), and there were 19 identical high-frequency codons. The CUB of Theaceae species is mainly affected by natural selection. The SCUO value of the 40 Theaceae species was 0.23 or 0.24, and the chloroplast gene expression level was moderate, according to MILC values. Additionally, we observed a positive correlation between the SCUO and MILC values, which indicated that CUB might affect gene expression. Furthermore, the phylogenetic analysis showed that the evolutionary relationships in these 40 Theaceae species were relatively conserved. A systematic study on the CUB and expression of Theaceae species provides further evidence for their evolution and phylogeny.

The comparative studies of complete chloroplast genomes in Actinidia (Actinidiaceae): novel insights into heterogenous variation, clpP gene annotation and phylogenetic relationships

The genus Actinidia, also called kiwifruit, is characterized with abundant balanced nutritional metabolites, including exceptionally high vitamin C content. However, the traditional classification could not fully reflect the actual Actinidia species' relationships, which need further revision through more accurate approaches. Compared to the nuclear genome, the chloroplast genome has simple heredity characteristics, conserved genome structure and small size, suitable for deciphering complicated species' phylogenetic relationships. Here, the genome-wide comprehensive comparative analyses were performed over 29 independent chloroplast genomes' sequences derived from 25 Actinidia taxa. The average genome size is 156,673.38 bp, with an average 37.20% GC content. The long repeat sequences rather than SSRs (simple sequence repeats) in Actinidia were revealed to be the causal agent leading to the chloroplast genome size expansion. The clpP gene sequences with exon merge and intron deletion were annotated in all the 29 chloroplast genomes tested, which has been previously reported to be lost in Actinidia species. Comprehensive sequence analyses indicated the distinct variation at the clpP gene locus was Actinidiaceae-specific, emerging after the Actinidiaceae-other Ericales species divergence. Four highly divergent sequences (i.e., rps16 ~ trnQ-UUG, rps4 ~ trnT-UGU, petA ~ psbJ, and rps12 ~ psbB) evolved in the LSC (large single-copy) and SSC (small single-copy) regions embodying rps12 ~ psbB (including clpP gene and its up/downstream noncoding sequence) were identified as variation hot spots in Actinidia species. Based on either LSC region alone, combined sequences of LSC and SSC or the whole chloroplast genome sequences, three identical phylogenetic trees of the 25 Actinidia taxa with relatively improved resolution were reconstructed, consistently supporting the reticulate evolutionary lineage in Actinidia. Our findings could help to better understand the evolution characteristics of chloroplast genomes and phylogenetic relationships among Actinidia species.

Rapid sequence evolution is associated with genetic incompatibilities in the plastid Clp complex

KEY MESSAGE

Replacing the native clpP1 gene in the Nicotiana plastid genome with homologs from different donor species showed that the extent of genetic incompatibilities depended on the rate of sequence evolution. The plastid caseinolytic protease (Clp) complex plays essential roles in maintaining protein homeostasis and comprises both plastid-encoded and nuclear-encoded subunits. Despite the Clp complex being retained across green plants with highly conserved protein sequences in most species, examples of extremely accelerated amino acid substitution rates have been identified in numerous angiosperms. The causes of these accelerations have been the subject of extensive speculation but still remain unclear. To distinguish among prevailing hypotheses and begin to understand the functional consequences of rapid sequence divergence in Clp subunits, we used plastome transformation to replace the native clpP1 gene in tobacco (Nicotiana tabacum) with counterparts from another angiosperm genus (Silene) that exhibits a wide range in rates of Clp protein sequence evolution. We found that antibiotic-mediated selection could drive a transgenic clpP1 replacement from a slowly evolving donor species (S. latifolia) to homoplasmy but that clpP1 copies from Silene species with accelerated evolutionary rates remained heteroplasmic, meaning that they could not functionally replace the essential tobacco clpP1 gene. These results suggest that observed cases of rapid Clp sequence evolution are a source of epistatic incompatibilities that must be ameliorated by coevolutionary responses between plastid and nuclear subunits.

The complete chloroplast genomes of three Hamamelidaceae species: Comparative and phylogenetic analyses

Hamamelidaceae is an important group that represents the origin and early evolution of angiosperms. Its plants have many uses, such as timber, medical, spice, and ornamental uses. In this study, the complete chloroplast genomes of Loropetalum chinense (R. Br.) Oliver, Corylopsis glandulifera Hemsl., and Corylopsis velutina Hand.-Mazz. were sequenced using the Illumina NovaSeq 6000 platform. The sizes of the three chloroplast genomes were 159,402 bp (C. glandulifera), 159,414 bp (C. velutina), and 159,444 bp (L. chinense), respectively. These chloroplast genomes contained typical quadripartite structures with a pair of inverted repeat (IR) regions (26,283, 26,283, and 26,257 bp), a large single-copy (LSC) region (88,134, 88,146, and 88,160 bp), and a small single-copy (SSC) region (18,702, 18,702, and 18,770 bp). The chloroplast genomes encoded 132-133 genes, including 85-87 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. The coding regions were composed of 26,797, 26,574, and 26,415 codons, respectively, most of which ended in A/U. A total of 37-43 long repeats and 175-178 simple sequence repeats (SSRs) were identified, and the SSRs contained a higher number of A + T than G + C bases. The genome comparison showed that the IR regions were more conserved than the LSC or SSC regions, while the noncoding regions contained higher variability than the gene coding regions. Phylogenetic analyses revealed that species in the same genus tended to cluster together. Chunia Hung T. Chang, Mytilaria Lecomte, and Disanthus Maxim. may have diverged early and Corylopsis Siebold & Zucc. was closely related to Loropetalum R. Br. This study provides valuable information for further species identification, evolution, and phylogenetic studies of Hamamelidaceae plants.

Uncovering the first complete plastome genomics, comparative analyses, and phylogenetic dispositions of endemic medicinal plant Ziziphus hajarensis (Rhamnaceae)

BACKGROUND

Ziziphus hajarensis is an endemic plant species well-distributed in the Western Hajar mountains of Oman. Despite its potential medicinal uses, little is known regarding its genomic architecture, phylogenetic position, or evolution. Here we sequenced and analyzed the entire chloroplast (cp) genome of Z. hajarensis to understand its genetic organization, structure, and phylogenomic disposition among Rhamnaceae species.

RESULTS

The results revealed the genome of Z. hajarensis cp comprised 162,162 bp and exhibited a typical quadripartite structure, with a large single copy (LSC) region of 895,67 bp, a small single copy (SSC) region of 19,597 bp and an inverted repeat (IR) regions of 26,499 bp. In addition, the cp genome of Z. hajarensis comprises 126 genes, including 82 protein-coding genes, eight rRNA genes, and 36 tRNA genes. Furthermore, the analysis revealed 208 microsatellites, 96.6% of which were mononucleotides. Similarly, a total of 140 repeats were identified, including 11 palindromic, 24 forward, 14 reverse, and 104 tandem repeats. The whole cp genome comparison of Z. hajarensis and nine other species from family Rhamnaceae showed an overall high degree of sequence similarity, with divergence among some intergenic spacers. Comparative phylogenetic analysis based on the complete cp genome, 66 shared genes and matK gene revealed that Z. hajarensis shares a clade with Z. jujuba and that the family Rhamnaceae is the closest family to Barbeyaceae and Elaeagnaceae.

CONCLUSION

All the genome features such as genome size, GC content, genome organization and gene order were highly conserved compared to the other related genomes. The whole cp genome of Z. hajarensis gives fascinating insights and valuable data that may be used to identify related species and reconstruct the phylogeny of the species.

Structural characterization of Platanthera ussuriensis chloroplast genome and comparative analyses with other species of Orchidaceae

BACKGROUND

The genus Tulotis has been classified into the genus Platanthera in the present taxonomic studies since the morphological characteristics of this genus is very similar to that of Platanthera. Platanthera ussuriensis, formerly named as Tulotis ussuriensis, is a small terrestrial orchid species and has been listed as wild plant under State protection (category II) in China. An improved understanding of the genomic information will enable future applications of conservation strategy as well as phylogenetic studies for this rare orchid species. The objective of this research was to characterize and compare the chloroplast genome of P. ussuriensis with other closely related species of Orchidaceae.

RESULTS

The chloroplast genome sequence of P. ussuriensis is 155,016 bp in length, which included a pair of inverted repeats (IRs) of 26,548 bp that separated a large single copy (LSC) region of 83,984 bp and a small single copy (SSC) region of 17,936 bp. The annotation contained a total of 132 genes, including 86 protein-coding genes, 38 tRNA genes and 8 rRNA genes. The simple sequence repeat (SSR) analysis showed that there were 104 SSRs in the chloroplast genome of P. ussuriensis. RNA editing sites recognition indicated 72 RNA editing events occurred, and all codon changes were C to T conversions. Comparative genomics showed that the chloroplast sequence of Platanthera related species were relatively conserved, while there were still some high variation regions that could be used as molecular markers. Moreover, Platanthera related species showed similar IR/SSC and IR/LSC borders. The phylogenetic analysis suggested that P. ussuriensis had a closer evolutionary relationship with P. japonica followed by the remaining Platanthera species.

CONCLUSION

Orchidaceae is a key group of biodiversity protection and also a hot spot group in the plant taxonomy and evolution studies due to their characteristics of high specialization and rapid evolution. This research determined the complete chloroplast genome of P. ussuriensis for the first time, and compared the sequence with other closely related orchid species. These results provide a foundation for future genomic and molecular evolution of the Orchidaceae species, and provide insights into the development of conservation strategy for Platanthera species.

Analysis of Chloroplast Genomes Provides Insights Into the Evolution of Agropyron

Plants of the Agropyron genus are important pasture resources, and they also play important roles in the ecological restoration. Chloroplast genomes are inherited from maternal parents, and they are important for studying species taxonomy and evolution. In this study, we sequenced the complete chloroplast genomes of five typical species of the Agropyron genus (eg., A. cristatum × A. desertorum Fisch. Schult, A. desertorum, A. desertorum Fisch. Schult. cv. Nordan, A. michnoi Roshev, and A. mongolicum Keng) using the Illumina NovaSeq platform. We found that these five chloroplast genomes exhibit a typical quadripartite structure with a conserved genome arrangement and structure. Their chloroplast genomes contain the large single-copy regions (LSC, 79,613 bp-79,634 bp), the small single-copy regions (SSC, 12,760 bp-12,768 bp), and the inverted repeat regions (IR, 43,060 bp-43,090 bp). Each of the five chloroplast genomes contains 129 genes, including 38 tRNA genes, eight rRNA genes, and 83 protein-coding genes. Among them, the genes trnG-GCC, matK, petL, ccsA, and rpl32 showed significant nucleotide diversity in these five species, and they may be used as molecular markers in taxonomic studies. Phylogenetic analysis showed that A. mongolicum is closely related to A. michnoi, while others have a closer genetic relationship with the Triticum genus.

Comprehensive analysis of chloroplast genome of Albizia julibrissin Durazz. (Leguminosae sp.)

The Albizia julibrissin chloroplasts have a classical chloroplast genome structure, containing 93 coding genes and 34 non-coding genes. Our research provides basic data for plant phylogenetic evolutionary studies. There is limited genomic information available for the important Chinese herb Albizia julibrissin Durazz. In this study, we constructed the chloroplast (Cp) genome of A. julibrissin. The length of the assembled Cp genome was 175,922 bp consisting of four conserved regions: a 5145 bp small single-copy (SSC) region, a 91,323 bp large single-copy (LSC) region, and two identical length-inverted repeat (IR) regions (39,725 bp). This Cp genome included 34 non-coding RNAs and 93 unique genes, the former contains 30 transfer and 4 ribosomal RNA genes. Gene annotation indicated some of the coding genes (82) in the A. julibrissin Cp genome classified in the Leguminosae family, with some to other related families (11). The results show that low GC content (36.9%) and codon bias towards A- or T-terminal codons may affect the frequency of gene codon usage. The sequence analysis identified 30 forward, 18 palindrome, and 1 reverse repeat > 30 bp length, and 149 simple sequence repeats (SSR). Fifty-five RNA editing sites in the Cp of A. julibrissin were predicted, most of which are C-to-U conversions. Analysis of the reverse repeat expansion or contraction and divergence area between several species, including A. julibrissin, was performed. The phylogenetic tree revealed that A. julibrissin was most closely related to Albizia odoratissima and Albizia bracteata, followed by Samanea saman, forming an evolutionary branch with Mimosa pudica and Leucaena trichandra. The research results are helpful for breeding and genetic improvement of A. julibrissin, and also provide valuable information for understanding the evolution of this plant.

Characterization of the complete chloroplast genome of Actinidia hemsleyana

The completed chloroplast genome sequence of Actinidia hemsleyana, collected from eastern China, was firstly determined using next-generation sequencing (NGS) and bioinformatic analysis. The length of cp genome was 156,845 bp, including a large single-copy (LSC) region of 88,666 bp and a small single-copy (SSC) region of 20,543 bp, which were separated by a pair of inverted repeats (IRs) of 23,818 bp. The cp genome contained 132 genes, including 83 protein-coding genes, 41 tRNA genes, and 8 rRNA genes. The overall GC content was 37.19%, whereas the corresponding values of the LSC, SSC, and IR regions were 35.46%, 31.09%, and 43.06%, respectively. Phylogenetic analysis based on whole cp genome sequences of 20 Actinidia species showed that A. hemsleyana has the closest relationship with A.latifolia.

Lineage-Specific Variation in IR Boundary Shift Events, Inversions, and Substitution Rates among Caprifoliaceae s.l. (Dipsacales) Plastomes

Caprifoliaceae s.l. plastid genomes (plastomes) show that one inversion and two inverted repeat boundary shifts occurred in the common ancestor of this family, after which the plastomes are generally conserved. This study reports plastome sequences of five additional species, Fedia cornucopiae, Valeriana fauriei, and Valerianella locusta from the subfamily Valerianoideae, as well as Dipsacus japonicus and Scabiosa comosa from the subfamily Dipsacoideae. Combined with the published plastomes, these plastomes provide new insights into the structural evolution of plastomes within the family. Moreover, the three plastomes from the subfamily Valerianoideae exhibited accelerated nucleotide substitution rates, particularly at synonymous sites, across the family. The patterns of accD sequence divergence in the family are dynamic with structural changes, including interruption of the conserved domain and increases in nonsynonymous substitution rates. In particular, the Valeriana accD gene harbors a large insertion of amino acid repeat (AAR) motifs, and intraspecific polymorphism with a variable number of AARs in the Valeriana accD gene was detected. We found a correlation between intron losses and increased ratios of nonsynonymous to synonymous substitution rates in the clpP gene with intensified positive selection. In addition, two Dipsacoideae plastomes revealed the loss of the plastid-encoded rps15, and a potential functional gene transfer to the nucleus was confirmed.

The complete chloroplast genome sequence of Rubus hirsutus Thunb. and a comparative analysis within Rubus species

Rubus hirsutus is a type of tonifying kidney-essence herb that belongs to the Rosaceae family, and has been commonly used to treat multiple diseases, such as polyuria, impotence, and infertility. In this study, we determined the complete chloroplast sequence of R. hirsutus and conduced a comparative analysis within the genus Rubus. The assembled chloroplast (cp.) genome is 156,380 bp in length with a GC content of 37.0% and shares a conserved quadripartite structure within the other cp. genomes in this genus. A total of 132 unique genes were annotated in the cp. genome of R. hirsutus, which contained 87 protein-coding genes, 37 tRNAs, and eight rRNAs. Seventeen duplicated genes were identified in the inverted repeats region. Furthermore, 70 simple sequence repeats and 35 long repeats were detected in total in the R. hirsutus chloroplast genome. Eight mutational hotspots were identified in the cp. genome of this species with higher nucleotide variations in non-coding regions than those of coding regions. Furthermore, the gene order, codon usage, and repeat sequence distribution were highly consistent in Rubus according to the results of a comparative analysis. A phylogenetic analysis indicated that there was a sister relationship between R. hirsutus and R. chingii. Overall, the complete chloroplast genome of R. hirsutus and the comparative analysis will help to further the evolutionary study, conservation, phylogenetic reconstruction, and development of molecular barcodes for the genus Rubus.

Screening and verification of extranuclear genetic markers in green tide algae from the Yellow Sea

Over the past decade, Ulva compressa, a cosmopolitan green algal species, has been identified as a component of green tides in the Yellow Sea, China. In the present study, we sequenced and annotated the complete chloroplast genome of U. compressa (alpha-numeric code: RD9023) and focused on the assessment of genome length, homology, gene order and direction, intron size, selection strength, and substitution rate. We compared the chloroplast genome with the mitogenome. The generated phylogenetic tree was analyzed based on single and aligned genes in the chloroplast genome of Ulva compared to mitogenome genes to detect evolutionary trends. U. compressa and U. mutabilis chloroplast genomes had similar gene queues, with individual genes exhibiting high homology levels. Chloroplast genomes were clustered together in the entire phylogenetic tree and shared several forward/palindromic/tandem repetitions, similar to those in U. prolifera and U. linza. However, U. fasciata and U. ohnoi were more divergent, especially in sharing complementary/palindromic repetitions. In addition, phylogenetic analyses of the aligned genes from their chloroplast genomes and mitogenomes confirmed the evolutionary trends of the extranuclear genomes. From phylogenetic analysis, we identified the petA chloroplast genes as potential genetic markers that are similar to the tufA marker. Complementary/forward/palindromic interval repetitions were more abundant in chloroplast genomes than in mitogenomes. Interestingly, a few tandem repetitions were significant for some Ulva subspecies and relatively more evident in mitochondria than in chloroplasts. Finally, the tandem repetition [GAAATATATAATAATA × 3, abbreviated as TRg)] was identified in the mitogenome of U. compressa and the conspecific strain U. mutabilis but not in other algal species of the Yellow Sea. Owing to the high morphological plasticity of U. compressa, the findings of this study have implications for the rapid non-sequencing detection of this species during the occurrence of green tides in the region.

The chloroplast genome of Prunus zhengheensis: Genome comparative and phylogenetic relationships analysis

Prunus zhengheensis is a novel species originated in Fujian province, China. However, there is no further information available on its classification and molecular biology study. In this study, we first report the complete chloroplast (cp) genome sequence of P. zhengheensis. The cp genome of P. zhengheensis is 158,106 bp and GC content is 36.73%, 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 86,321 bp, 18,999 bp and 26,393 bp, respectively. The cp genome of P. zhengheensis contains 130 genes, and 242 SSRs are identified in the cp genome. The comparative analysis of cp genomes in eight Prunus plants demonstrates the subtle divergences occur in the protein-coding gene rps18, rps12, psbF, rpl33, matK, and rbcL, and that the K_A/K_S nucleotide substitution ratio of the ndhF of P. zhengheensis and P. armeniaca is 1.79636. The phylogenetic results indicate that the P. zhengheensis is closely related to P. mume, compared to other species of Prunus. Our research results provide the important genomic information for molecular phylogeny of P. zhengheensis.

Characterization of the Complete Chloroplast Genome of Buddleja Lindleyana

BACKGROUND

Buddleja lindleyana Fort., which belongs to the Loganiaceae with a distribution throughout the tropics, is widely used as an ornamental plant in China. Buddleja contains several morphologically similar species, which need to be identified by molecular identification. But there is little molecular research on the genus Buddleja.

OBJECTIVE

Using molecular biology techniques to sequence and analyze the complete chloroplast (cp) genome of B. lindleyana.

METHODS

According to next-generation sequencing to sequence the genome data, a series of bioinformatics software were used to assembly and analysis the molecular structure of cp genome of B. lindleyana.

RESULTS

The complete cp genome of B. lindleyana is a circular 154,487-bp-long molecule with a GC content of 38.1%. It has a familiar quadripartite structure, including a large single-copy region (LSC; 85,489 bp), a small single-copy region (SSC; 17,898bp) and a pair of inverted repeats (IRs; 25,550 bp). A total of 133 genes were identified in the genome, including 86 protein-coding genes, 37 tRNA genes, 8 rRNA genes and 2 pseudogenes.

CONCLUSIONS

These results suggested that B. lindelyana cp genome could be used as a potential genomic resource to resolve the phylogenetic positions and relationships of Loganiaceae, and will offer valuable information for future research in the identification of Buddleja species and will conduce to genomic investigations of these species.

Chloroplast genome sequence of Chongming lima bean (Phaseolus lunatus L.) and comparative analyses with other legume chloroplast genomes

Background

Lima bean (Phaseolus lunatus L.) is a member of subfamily Phaseolinae belonging to the family Leguminosae and an important source of plant proteins for the human diet. As we all know, lima beans have important economic value and great diversity. However, our knowledge of the chloroplast genome level of lima beans is limited.

Results

The chloroplast genome of lima bean was obtained by Illumina sequencing technology for the first time. The Cp genome with a length of 150,902 bp, including a pair of inverted repeats (IRA and IRB 26543 bp each), a large single-copy (LSC 80218 bp) and a small single-copy region (SSC 17598 bp). In total, 124 unique genes including 82 protein-coding genes, 34 tRNA genes, and 8 rRNA genes were identified in the P. lunatus Cp genome. A total of 61 long repeats and 290 SSRs were detected in the lima bean Cp genome. It has a typical 50 kb inversion of the Leguminosae family and an 70 kb inversion to subtribe Phaseolinae. rpl16, accD, petB, rsp16, clpP, ndhA, ndhF and ycf1 genes in coding regions was found significant variation, the intergenic regions of trnk-rbcL, rbcL-atpB, ndhJ-rps4, psbD-rpoB, atpI-atpA, atpA-accD, accD-psbJ, psbE-psbB, rsp11-rsp19, ndhF-ccsA was found in a high degree of divergence. A phylogenetic analysis showed that P. lunatus appears to be more closely related to P. vulgaris, V.unguiculata and V. radiata.

Conclusions

The characteristics of the lima bean Cp genome was identified for the first time, these results will provide useful insights for species identification, evolutionary studies and molecular biology research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07467-8.

Characterization of the complete chloroplast genome sequence of Actinidia kolomikta

The whole complete chloroplast (cp) genome sequence of Actinidia kolomikta (formerly A. maloide) was sequenced and assembled from Illumina paired-end sequencing. The cp genome of A. arguta was 157,425 bp long, containing a large single copy region (LSC) of 88,498 bp and a small single copy region (SSC) of 20,475 bp and a pair of inverted repeat regions (IR) of 242,266 bp. It contained 113 different genes, including 79 protein-coding genes, 30 tRNA genes, 4 ribosomal RNA genes. Phylogenetic analysis using concatenated alignments of whole cp genome sequences revealed that A. kolomikta was the sister group to all other groups of Actinidia

Green giant-a tiny chloroplast genome with mighty power to produce high-value proteins: history and phylogeny

Free-living cyanobacteria were entrapped by eukaryotic cells ~2 billion years ago, ultimately giving rise to chloroplasts. After a century of debate, the presence of chloroplast DNA was demonstrated in the 1960s. The first chloroplast genomes were sequenced in the 1980s, followed by ~100 vegetable, fruit, cereal, beverage, oil and starch/sugar crop chloroplast genomes in the past three decades. Foreign genes were expressed in isolated chloroplasts or intact plant cells in the late 1980s and stably integrated into chloroplast genomes, with typically maternal inheritance shown in the 1990s. Since then, chloroplast genomes conferred the highest reported levels of tolerance or resistance to biotic or abiotic stress. Although launching products with agronomic traits in important crops using this concept has been elusive, commercial products developed include enzymes used in everyday life from processing fruit juice, to enhancing water absorption of cotton fibre or removal of stains as laundry detergents and in dye removal in the textile industry. Plastid genome sequences have revealed the framework of green plant phylogeny as well as the intricate history of plastid genome transfer events to other eukaryotes. Discordant historical signals among plastid genes suggest possible variable constraints across the plastome and further understanding and mitigation of these constraints may yield new opportunities for bioengineering. In this review, we trace the evolutionary history of chloroplasts, status of autonomy and recent advances in products developed for everyday use or those advanced to the clinic, including treatment of COVID-19 patients and SARS-CoV-2 vaccine.

Chloroplast Genome Analysis of Two Medicinal Coelogyne spp. (Orchidaceae) Shed Light on the Genetic Information, Comparative Genomics, and Species Identification

Although the medicinal properties of Coelogyne spp. have been previously studied, there is little genomic information providing a valuable tool for the plant taxonomy, conservation, and utilization of this genus. This study used the next-generation MiSeq sequencing platform to characterize the chloroplast (cp) genomes of Coelogyne fimbriata and Coelogyne ovalis. The Maximum Likelihood (ML) and Bayesian (BI) methods were employed to confirm the phylogenetic position of two Coelogyne species based on the whole chloroplast genome sequences. Additionally, we developed eight new primers based on the two cp genomes' medium variable regions and evaluated the transferability to another 16 Coelogyne species. We constructed phylogenetic trees including 18 Coelogyne species and four outgroup species using the chloroplast fragments with the ML method. Our results showed that the cp genomes of C. fimbriata and C. ovalis contained a small single-copy region (18,839 and 18,851 bp, respectively) and a large single-copy region (87,606 and 87,759 bp, respectively), separated by two same-length inverted-repeat regions (26,675 bp in C. fimbriata and 26,715 bp C. ovalis, respectively). They all contained 86 protein-coding genes, 38 tRNA genes, and eight rRNA genes, revealing strong structure and gene content similarities. The phylogenetic analysis indicated a close relationship between the genera Coelogyne and Pleione. The newly developed primers revealed good transferability among the Coelogyne taxa and provided enough variable sites to distinguish C. fimbriata and C. ovalis. The two complete cp genomes and the eight new primers of Coelogyne provide new genomic data for further studies on phylogenomics, population genetics, and evolutionary history of Coelogyne taxa.

Phylogenetic and Comparative Analyses of Complete Chloroplast Genomes of Chinese Viburnum and Sambucus (Adoxaceae)

Phylogenetic analyses of complete chloroplast genome sequences have yielded significant improvements in our understanding of relationships in the woody flowering genus Viburnum (Adoxaceae, Dipsacales); however, these relationships were evaluated focusing only on Viburnum species within Central and South America and Southeast Asia. By contrast, despite being a hotspot of Viburnum diversity, phylogenetic relationships of Viburnum species in China are less well known. Here, we characterized the complete chloroplast (cp) genomes of 21 Viburnum species endemic to China, as well as three Sambucus species. These 24 plastomes were highly conserved in genomic structure, gene order and content, also when compared with other Adoxaceae. The identified repeat sequences, simple sequence repeats (SSRs) and highly variable plastid regions will provide potentially valuable genetic resources for further population genetics and phylogeographic studies on Viburnum and Sambucus. Consistent with previous combined phylogenetic analyses of 113 Viburnum species, our phylogenomic analyses based on the complete cp genome sequence dataset confirmed the sister relationship between Viburnum and the Sambucus-Adoxa-Tetradoxa-Sinadoxa group, the monophyly of four recognized sections in Flora of China (i.e., Viburnum sect. Tinus, Viburnum sect. Solenotinus, Viburnum sect. Viburnum and Viburnum sect. Pseudotinus) and the nonmonophyly of Viburnum sect. Odontotinus and Viburnum sect. Megalotinus. Additionally, our study confirmed the sister relationships between the clade Valvatotinus and Viburnum sect. Pseudotinus, as well as between Viburnum sect. Opulus and the Odontotinus-Megalotinus group. Overall, our results clearly document the power of the complete cp genomes in improving phylogenetic resolution, and will contribute to a better understanding of plastome evolution in Chinese Adoxaceae.