The evolutionary fate of rpl32 and rps16 losses in the Euphorbia schimperi (Euphorbiaceae) plastome

Taken to extremes: Loss of plastid rpl32 in Streptophyta and Cuscuta's unconventional solution for its replacement

The evolution of plant genomes is riddled with exchanges of genetic material within one plant (endosymbiotic gene transfer/EGT) and between unrelated plants (horizontal gene transfer/HGT). These exchanges have left their marks on plant genomes. Parasitic plants with their special evolutionary niche provide ample examples for these processes because they are under a reduced evolutionary pressure to maintain autotrophy and thus to conserve their plastid genomes. On the other hand, the close physical connections with different hosts enabled them to acquire genetic material from other plants. Based on an analysis of an extensive dataset including the parasite Cuscuta campestris and other parasitic plant species, we identified a unique evolutionary history of rpl32 genes coding for an essential plastid ribosomal subunit in Cuscuta. Our analysis suggests that the gene was most likely sequestered by HGT from a member of the Oxalidales order serving as host to an ancestor of the Cuscuta subgenus Grammica. Oxalidales had suffered an ancestral EGT of rpl32 predating the evolution of the genus Cuscuta. The HGT subsequently relieved the plastid rpl32 from its evolutionary constraint and led to its loss from the plastid genome. The HGT-based acquisition in Cuscuta is supported by a high sequence similarity of the mature L32 protein between species of the subgenus Grammica and representatives of the Oxalidales, and by a surprisingly conserved transit peptide, whose functionality in Cuscuta was experimentally verified. The findings are discussed in view of an overall pattern of EGT events for plastid ribosomal subunits in Streptophyta.

Elucidating the evolutionary dynamics of parasitism in Cuscuta: in-depth phylogenetic reconstruction and extensive plastomes reduction

BACKGROUND

The genus Cuscuta L. (Convolvulaceae), commonly known as dodder, is a holoparasite plant that relies on host plants for nutrition, leading to significant genomic changes, particularly in plastomes. This dependency has led to significant reductions and modifications in their plastomes compared to autotrophic plants. In contrast to the well-conserved plastomes of photosynthetic plants, Cuscuta exhibits substantial genomic reductions reflecting the loss of photosynthetic functions and associated genes.

RESULT

This study examines eight plastomes within Cuscuta and reconstructs the phylogenetic relationships among 40 Cuscuta taxa using five other genera as an outgroup. The size of plastid genome varies significantly, with the smallest being 60 kb and the largest 121 kb, highlighting extensive genomic reduction. In special cases, the subgenera Cuscuta exhibit the loss of inverted repeats, distinguishing from them other subge within the Cuscuta genus. This reduction is most pronounced in genes related to photosynthesis, such as atp, pet, psa, psb, and ycf genes, particularly in the subg. Grammica (Lour.) Peter. The study also notes the frequent and independent loss of the plastid genes infA, rpl23, rpl32, rps15, and rps16 across various angiosperm lineages, often involving transfer to the nuclear genome. In parasitic plants like Cuscuta, the ndh genes, crucial for photosynthesis, are often lost. The study also highlights that in the subg. Grammica, the matK and rpo genes, along with trnR-ACG genes, are lost in parallel, indicating that these parasitic plants do not need matK and rpo genes after the loss of ndh genes for survival. Analysis of selective relaxation pressure on plastid genes shows a reductive trend, with genes such as atp, pet, psa, psb, rpo, and ycf progressively becoming pseudogenes over time, with housekeeping genes like rpl and rps expected to follow. However, the pseudogenization process is specific to the subg. Grammica, Pachystigma (Engelm.) Baker & C.H.Wright, and Cuscuta, rather than in the subg. Monogynella (Des Moul.) Peter, Engl. & Prantl (ancient clade species).

CONCLUSION

The study of Cuscuta plastomes reveals the profound impact of parasitism on genome evolution, highlighting the complex interplay of gene retention and loss through phylogenomic approaches. This research enriches our understanding of plant genome evolution and the intricate host-parasite relationships. It also sheds light on the evolutionary history and genomic adaptations of Cuscuta, illustrating the diverse strategies enabling subg. Grammica, Pachystigma, Cuscuta, and Monogynella thrive as parasitic species. These findings provide valuable insights into the molecular mechanisms underlying parasitism and its impact on plastid genome organization.

Chloroplast genomics unravels taxonomic status and evolutionary relationships in Euphorbia fauriei and E. pekinensis

Euphorbia fauriei, a perennial plant endemic to South Korea, exhibits both morphological characteristics and intricate genetic identities akin to E. pekinensis. This study aimed to provide clarity on the taxonomic status of E. faurirei and E. pekinensis through a comprehensive chloroplast genome analysis. Additionally, we sequenced the Acalypha australis chloroplast genome as an outgroup for the construction of a phylogenetic tree with other Euphorbia species. The three chloroplast genomes, ranging from 162,834 bp to 168,832 bp, displayed typical quadripartite structures. The Euphorbia genomes contained 111 unique genes, whereas the A. australis genome contained two additional protein-coding genes (PCGs), rpl32 and rps16. Comparative analysis unveiled the loss of rpl32 and rps16 as synapomorphic characteristics in Euphorbia, whereas the loss of infA occurred across Euphorbiaceae. High collinearity and sequence similarity were observed among Euphorbia species, accompanied by significant inversion regions in the E. pekinensis chloroplast genomes from China and Japan, indicating regional genetic variability. Nucleotide substitution analysis revealed that the ndh group exhibited the highest KA/KS values (0.224), whereas the atp, psa, psb and rps groups had the lowest. Phylogenomic analysis utilizing whole genomes, PCGs, and intron and intergenic regions consistently demonstrated that E. pekinensis from South Korea clusters closely with E. fauriei. These findings challenge the current taxonomic distinction between E. pekinensis and E. fauriei in Korea, suggesting that while they exhibit distinct characteristics, E. fauriei should be considered a closely related subspecies rather than the same species as E. pekinensis. This study emphasizes the need for population studies to clarify the taxonomic relationships between E. pekinensis and E. fauriei.

The plastomes of Lepismium cruciforme (Vell.) Miq and Schlumbergera truncata (Haw.) Moran reveal tribe-specific rearrangements and the first loss of the trnT-GGU gene in Cactaceae

BACKGROUND

Recent studies have revealed atypical features in the plastomes of the family Cactaceae, the largest lineage of succulent species adapted to arid and semi-arid regions. Most plastomes sequenced to date are from short-globose and cylindrical cacti, while little is known about plastomes of epiphytic cacti. Published cactus plastomes reveal reduction and complete loss of IRs, loss of genes, pseudogenization, and even degeneration of tRNA structures. Aiming to contribute with new insights into the plastid evolution of Cactaceae, particularly within the tribe Rhipsalideae, we de novo assembled and analyzed the plastomes of Lepismium cruciforme and Schlumbergera truncata, two South American epiphytic cacti.

METHODS AND RESULTS

Our data reveal many gene losses in both plastomes and the first loss of functionality of the trnT-GGU gene in Cactaceae. The trnT-GGU is a pseudogene in L. cruciforme plastome and appears to be degenerating in the tribe Rhipsalideae. Although the plastome structure is conserved among the species of the tribe Rhipsalideae, with tribe-specific rearrangements, we mapped around 200 simple sequence repeats and identified nine nucleotide polymorphism hotspots, useful to improve the phylogenetic resolutions of the Rhipsalideae. Furthermore, our analysis indicated high gene divergence and rapid evolution of RNA editing sites in plastid protein-coding genes in Cactaceae.

CONCLUSIONS

Our findings show that some characteristics of the Rhipsalideae tribe are conserved, such as plastome structure with IRs containing only the ycf2 and two tRNA genes, structural degeneration of the trnT-GGU gene and ndh complex, and lastly, pseudogenization of rpl33 and rpl23 genes, both plastid translation-related genes.

Comparative analysis the chloroplast genomes of Celastrus (Celastraceae) species: Provide insights into molecular evolution, species identification and phylogenetic relationships

BACKGROUND

The genus Celastrus is an important medicinal plant resource. The similarity of morphology and the lack of complete chloroplast genome analysis have significantly impeded the exploration of species identification, molecular evolution and phylogeny of Celastrus.

PURPOSE

In order to resolve the phylogenic controversy of Celastrus species, the chloroplast genome comparative analysis was performed to provide genetic evidence.

METHODS

In this study, we collected and sequenced ten chloroplast genomes of Celastrus species from China and downloaded three chloroplast genomes from the databases. The chloroplast genomes were compared and analyzed to explore their characteristics and evolution. Furthermore, the phylogenetic relationships of Celastrus species were inferred based on the whole chloroplast genomes and protein-coding genes.

RESULTS

All the 13 Celastrus species chloroplast genomes showed a typical quadripartite structure with genome sizes ranging from 155,113 to 157,366 bp. The intron loss of the rps16 gene occurred in all the 13 Celastrus species. The GC content, gene sequence, repeat types and codon bias pattern were highly conserved. Ten highly variation regions were identified, which can be used as potential DNA markers in molecular identification of Celastrus species. Eight genes, including accD, atp4, ndhB, rpoC1, rbcL, rpl2, rpl20 and ycf1, were detected to experience positive selection. Phylogenetic analysis showed that Celastrus was a monophyletic group and Tripterygium was the closest sister-group. Noteworthy, C. gemmatus Loes. and C. orbiculatus Thunb. can be discriminated using the chloroplast genome as a super barcode. The comparative and phylogenetic analysis results proposed that C. tonkinensis Pitard. was the synonym of C. hindsii Benth.

CONCLUSION

The comparative analysis of the Celastrus chloroplast genomes can provide comprehensive genetic evidence for molecular evolution, species identification and phylogenetic relationships.

Decoding the Chloroplast Genome of Tetrastigma (Vitaceae): Variations and Phylogenetic Selection Insights

Tetrastigma (Vitaceae) is known for its ornamental, medicinal, and ecological significance. However, the structural and variational characteristics of the Tetrastigma chloroplast genome and their impact on phylogenetic relationships remain underexplored. This study utilized bioinformatics methods to assemble and annotate the chloroplast genomes of 10 Tetrastigma species and compare them with five previously sequenced species. This study analyzed gene composition, simple sequence repeats, and codon usage patterns, revealing a high A/T content, uniquely identified pentanucleotide repeats in five species and several preferred codons. In addition, comparative analyses were conducted of the chloroplast genomes of 15 Tetrastigma species, examining their structural differences and identifying polymorphic hotspots (rps16, rps16-trnQ, trnS, trnD, psbC-trnS-psbZ, accD-psaI, psbE-petL-petG, etc.) suitable for DNA marker development. Furthermore, phylogenetic and selective pressure analyses were performed based on the chloroplast genomes of these 15 Tetrastigma species, validating and elucidating intra-genus relationships within Tetrastigma. Futhermore, several genes under positive selection, such as atpF and accD, were identified, shedding light on the adaptive evolution of Tetrastigma. Utilizing 40 Vitaceae species, the divergence time of Tetrastigma was estimated, clarifying the evolutionary relationships within Tetrastigma relative to other genera. The analysis revealed diverse divergences of Tetrastigma in the Miocene and Pliocene, with possible ancient divergence events before the Eocene. Furthermore, family-level selective pressure analysis identified key features distinguishing Tetrastigma from other genera, showing a higher degree of purifying selection. This research enriches the chloroplast genome data for Tetrastigma and offers new insights into species identification, phylogenetic analysis, and adaptive evolution, enhancing our understanding of the genetic diversity and evolutionary history of these species.

Characterization, comparison, and phylogenetic analyses of chloroplast genomes of Euphorbia species

The genus Euphorbia (Euphorbiaceae) has near-cosmopolitan distribution and serves as a significant resource for both ornamental and medicinal purposes. Despite its economic importance, Euphorbia's taxonomy has long been challenged by the intricate nature of morphological traits exhibiting high levels of convergence. While molecular markers are essential for phylogenetic studies, their availability for Euphorbia has been limited. To address this gap, we conducted comparative analyses focusing on the chloroplast (CP) genomes of nine Euphorbia species, incorporating three newly sequenced and annotated accessions. In addition, phylogenetic informativeness and nucleotide diversity were computed to identify candidate markers for phylogenetic analyses among closely related taxa in the genus. Our investigation revealed relatively conserved sizes and structures of CP genomes across the studied species, with notable interspecific variations observed primarily in non-coding regions and IR/SC borders. By leveraging phylogenetic informativeness and nucleotide diversity, we identified rpoB gene as the optimal candidate for species delimitation and shallow-level phylogenetic inference within the genus. Through this comprehensive analysis of CP genomes across multiple taxa, our study sheds light on the evolutionary dynamics and taxonomic intricacies of Euphorbia, offering valuable insights into its CP genome evolution and taxonomy.

Intraspecific differentiation of Lindera obtusiloba as revealed by comparative plastomic and evolutionary analyses

Lindera obtusiloba Blume is the northernmost tree species in the family Lauraceae, and it is a key species in understanding the evolutionary history of this family. The species of L. obtusiloba in East Asia has diverged into the Northern and Southern populations, which are geographically separated by an arid belt. Though the morphological differences between populations have been observed and well documented, intraspecific variations at the plastomic level have not been systematically investigated to date. Here, ten chloroplast genomes of L. obtusiloba individuals were sequenced and analyzed along with three publicly available plastomes. Comparative plastomic analysis suggests that both the Northern and the Southern populations share similar overall structure, gene order, and GC content in their plastomes although the size of the plasome and the level of intraspecific variability do vary between the two populations. The Northern have relatively larger plastomes while the Southern population possesses higher intraspecific variability, which could be attributed to the complexity of the geological environments in the South. Phylogenomic analyses also support the split of the Northern and Southern clades among L. obtusiloba individuals. However, there is no obvious species boundary between var. obtusiloba and var. heterophylla in the Southern population, indicating that gene flow could still occur between these two varieties, and this could be used as a good example of reticulate evolution. It is also found that a few photosynthesis-related genes are under positive selection, which is mainly related to the geological and environmental differences between the Northern and the Southern regions. Our results provide a reference for phylogenetic analysis within species and suggest that phylogenomic analyses with a sufficient number of nuclear and chloroplast genomic target loci from widely distributed individuals could provide a deeper understanding of the population evolution of the widespread species.

Insight into chloroplast genome structural variation of the Mongolian endemic species Adonis mongolica (Ranunculaceae) in the Adonideae tribe

Adonis mongolica is a threatened species that is endemic to Mongolia. It is a medicinal plant from the Adonis genus and has been used to treat heart diseases. However, the genomics and evolution of this species have not been thoroughly studied. We sequenced the first complete plastome of A. mongolica and compared it with ten Adonideae species to describe the plastome structure and infer phylogenetic relationships. The complete plastome of A. mongolica was 157,521 bp long and had a typical quadripartite structure with numerous divergent regions. The plastomes of Adonideae had relatively constant genome structures and sizes, except for those of Adonis. The plastome structure was consistent across Adonis. We identified a 44.8 kb large-scale inversion within the large single-copy region and rpl32 gene loss in the Adonis plastomes compared to other members of the Adonideae tribe. Additionally, Adonis had a smaller plastome size (156,917-157,603 bp) than the other genera within the tribe (159,666-160,940 bp), which was attributed to deletions of intergenic regions and partial and complete gene losses. These results suggested that an intramolecular mutation occurred in the ancestor of the Adonis genus. Based on the phylogenetic results, Adonis separated earlier than the other genera within the Adonideae tribe. The genome structures and divergences of specific regions in the Adonis genus were unique to the Adonideae tribe. This study provides fundamental knowledge for further genomic research in Mongolia and a better understanding of the evolutionary history of endemic plants.

New Insights Into The Evolution of Chloroplast Genomes in Ochna Species (Ochnaceae, Malpighiales)

Ochnaceae DC. includes more than 600 species that exhibit potential values for environmental ecology, ornamental, pharmaceutical, and timber industries. Although studies on phylogeny and phytochemicals have been intensively conducted, chloroplast genome data of Ochnaceae species have not been fully explored. In this study, the next-generation sequencing method was used to sequence the chloroplast genomes of Ochna integerrima and Ochna serrulata which were 157 329 and 157 835 bp in length, respectively. These chloroplast genomes had a quadripartite structure and contained 78 protein-coding genes, 30 tRNAs, and 4 rRNAs. Comparative analysis revealed 8 hypervariable regions, including trnK_UUU-trnQ_UUG, rpoB-psbM, trnS_GGA-rps4, accD-psaI, rpl33-rps18, rpl14-rpl16, ndhF-trnL_UAG, and rps15-ycf1 among 6 Ochnaceae taxa. Additionally, there were shared and unique repeats among 6 examined chloroplast genomes. The notable changes were the loss of rpl32 in Ochna species and the deletion of rps16 exon 2 in O. integerrima compared to other taxa. This study is the first comprehensive comparative genomic analysis of complete chloroplast genomes of Ochna species and related taxa in Ochnaceae. Consequently, the current study provides initial results for further research on genomic evolution, population genetics, and developing molecular markers in Ochnaceae and related taxa.

Complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens: genome structures, comparative and phylogenetic analysis

The genus Cerastium includes about 200 species that are mostly found in the temperate climates of the Northern Hemisphere. Here we report the complete chloroplast genomes of Cerastium alpinum, C. arcticum and C. nigrescens. The length of cp genomes ranged from 147,940 to 148,722 bp. Their quadripartite circular structure had the same gene organization and content, containing 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Repeat sequences varied from 16 to 23 per species, with palindromic repeats being the most frequent. The number of identified SSRs ranged from 20 to 23 per species and they were mainly composed of mononucleotide repeats containing A/T units. Based on Ka/Ks ratio values, most genes were subjected to purifying selection. The newly sequenced chloroplast genomes were characterized by a high frequency of RNA editing, including both C to U and U to C conversion. The phylogenetic relationships within the genus Cerastium and family Caryophyllaceae were reconstructed based on the sequences of 71 protein-coding genes. The topology of the phylogenetic tree was consistent with the systematic position of the studied species. All representatives of the genus Cerastium were gathered in a single clade with C. glomeratum sharing the least similarity with the others.

Plastomes of Garcinia mangostana L. and Comparative Analysis with Other Garcinia Species

The two varieties of mangosteen (Garcinia mangostana L.) cultivated in Malaysia are known as Manggis and Mesta. The latter is preferred for its flavor, texture, and seedlessness. Here, we report a complete plastome (156,580 bp) of the Mesta variety that was obtained through a hybrid assembly approach using PacBio and Illumina sequencing reads. It encompasses a large single-copy (LSC) region (85,383 bp) and a small single-copy (SSC) region (17,137 bp) that are separated by 27,230 bp of inverted repeat (IR) regions at both ends. The plastome comprises 128 genes, namely, 83 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The plastome of the Manggis variety (156,582 bp) obtained from reference-guided assembly of Illumina reads was found to be nearly identical to Mesta except for two indels and the presence of a single-nucleotide polymorphism (SNP). Comparative analyses with other publicly available Garcinia plastomes, including G. anomala, G. gummi-gutta, G. mangostana var. Thailand, G. oblongifolia, G. paucinervis, and G. pedunculata, found that the gene content, gene order, and gene orientation were highly conserved among the Garcinia species. Phylogenomic analysis divided the six Garcinia plastomes into three groups, with the Mesta and Manggis varieties clustered closer to G. anomala, G. gummi-gutta, and G. oblongifolia, while the Thailand variety clustered with G. pedunculata in another group. These findings serve as future references for the identification of species or varieties and facilitate phylogenomic analysis of lineages from the Garcinia genus to better understand their evolutionary history.

The complete chloroplast genome of Viola grypoceras (Violaceae)

We constructed and characterized the chloroplast genome of Viola grypoceras via de novo assembly of Illumina data. The complete circular chloroplast genome is 158,357 bp long and contains four parts: a large single-copy (LSC) region of 86,764 bp, a small single-copy (SSC) region of 17,345 bp, and two inverted-repeat regions (IRa and IRb) of 27,124 bp each. Genome annotation predicted that this genome harbors 111 genes, comprising 77 protein-coding genes, 30 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes. Phylogenetic analysis demonstrated that V. grypoceras shares a close systematic relationship with V. mirabilis and V. websteri by forming a basal clade in the genus Viola.

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

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

Gene loss, genome rearrangement, and accelerated substitution rates in plastid genome of Hypericum ascyron (Hypericaceae)

BACKGROUND

Comparative genomic analysis exhibits dynamic evolution of plastid genome (plastome) in the clusioid clade of Malpighiales, which comprise five families, including multiple inversions and gene losses. Little is known about the plastome evolution in Hypericaceae, a large family in the clade. Only the plastome of one species, Cratoxylum cochinchinense, has been published.

RESULTS

We generated a complete plastome sequence for Hypericum ascyron, providing the first complete plastome from the tribe Hypericeae (Hypericaceae). The H. ascyron plastome exhibits dynamic changes in gene and intron content, structure, and sequence divergence compared to the C. cochinchinense plastome from the tribe Cratoxyleae (Hypericaceae). Transcriptome data determined the evolutionary fate of the missing plastid genes infA, rps7, rps16, rpl23, and rpl32 in H. ascyron. Putative functional transfers of infA, rps7, and rpl32 were detected to the nucleus, whereas rps16 and rpl23 were substituted by nuclear-encoded homologs. The plastid rpl32 was integrated into the nuclear-encoded SODcp gene. Our findings suggested that the transferred rpl32 had undergone subfunctionalization by duplication rather than alternative splicing. The H. ascyron plastome rearrangements involved seven inversions, at least three inverted repeat (IR) boundary shifts, which generated gene relocations and duplications. Accelerated substitution rates of plastid genes were observed in the H. ascyron plastome compared with that of C. cochinchinense plastid genes. The higher substitution rates in the accD and clpP were correlated with structural change, including a large insertion of amino acids and losses of two introns, respectively. In addition, we found evidence of positive selection of the clpP, matK, and rps3 genes in the three branches related to H. ascyron. In particular, the matK gene was repeatedly under selection within the family Hypericaceae. Selective pressure in the H. ascyron matK gene was associated with the loss of trnK-UUU and relocation into the IR region.

CONCLUSIONS

The Hypericum ascyron plastome sequence provides valuable information for improving the understanding of plastome evolution among the clusioid of the Malpighiales. Evidence for intracellular gene transfer from the plastid to the nucleus was detected in the nuclear transcriptome, providing insight into the evolutionary fate of plastid genes in Hypericaceae.

An integrated study of Violae Herba (Viola philippica) and five adulterants by morphology, chemical compositions and chloroplast genomes: insights into its certified plant origin

BACKGROUND

Viola philippica Cav. is the only original plant for Violae Herba, as described in the Chinese Pharmacopoeia. The quality of this crude drug is affected by several adulterants from congeneric Viola species, and the authentic plant origin of Violae Herba is still controversial. Genome-based identification offers abundant genetic information and potential molecular markers that can be used for the authentication of closely related species. This study aims to investigate the certified origin of Violae Herba and to develop more effective markers for these easily confused species at the genetic level.

METHODS

We compared the morphology and chemical composition of 18 batches of commercial samples and six widespread medicinal Viola plants used as Violae Herba or its substitutes by TLC and HPLC-Triple-TOF-MS/MS analyses. The complete chloroplast genomes of these species were sequenced and analyzed, including the general features, repeat sequences, mutational hotspots and phylogeny. The complete chloroplast genomes used as superbarcodes and some specific barcodes screened from mutational hotspots were tested for their ability to distinguish Viola species.

RESULTS

A comparative study showed that Violae Herba is a multi-origin traditional Chinese medicine. Commercial decoction pieces and the standard reference drug were mainly derived from V. prionantha, clashing with the record in the Chinese Pharmacopoeia. Chloroplast genome analyses of V. philippica and five adulterants indicated that sequence divergence was relatively low within Viola species. By tree-based approaches, the complete chloroplast genomes showed a better discrimination ability and phylogenetic resolution for each Viola species. These results indicate that the whole chloroplast genomes can be used as superbarcodes to differentiate Viola medicinal plants. More specific DNA barcodes could be further developed from the Viola chloroplast genomes for more efficient and rapid identification of commercial Violae Herba and its adulterants.

CONCLUSIONS

This study has implications for chloroplast genome-based phylogenetic analysis and the authentication of multiple Viola species used as Violae Herba. The legal origin recorded in the Chinese Pharmacopoeia should be further revised to V. prionantha, in line with the commercial Violae Herba in the TCM markets.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

Chloroplast Phylogenomic Analyses Reveal a Maternal Hybridization Event Leading to the Formation of Cultivated Peanuts

Peanuts (Arachis hypogaea L.) offer numerous healthy benefits, and the production of peanuts has a prominent role in global food security. As a result, it is in the interest of society to improve the productivity and quality of peanuts with transgenic means. However, the lack of a robust phylogeny of cultivated and wild peanut species has limited the utilization of genetic resources in peanut molecular breeding. In this study, a total of 33 complete peanut plastomes were sequenced, analyzed and used for phylogenetic analyses. Our results suggest that sect. Arachis can be subdivided into two lineages. All the cultivated species are contained in Lineage I with AABB and AA are the two predominant genome types present, while species in Lineage II possess diverse genome types, including BB, KK, GG, etc. Phylogenetic studies also indicate that all allotetraploid cultivated peanut species have been derived from a possible maternal hybridization event with one of the diploid Arachis duranensis accessions being a potential AA sub-genome ancestor. In addition, Arachis monticola, a tetraploid wild species, is placed in the same group with all the cultivated peanuts, and it may represent a transitional species, which has been through the recent hybridization event. This research could facilitate a better understanding of the taxonomic status of various Arachis species/accessions and the evolutionary relationship among them, and assists in the correct and efficient use of germplasm resources in breeding efforts to improve peanuts for the benefit of human beings.

Complete chloroplast genome sequences of five Bruguiera species (Rhizophoraceae): comparative analysis and phylogenetic relationships

Bruguiera is a genus of true mangroves that are mostly distributed in the Indo-West Pacific region. However, the number of published whole chloroplast genome sequences of Bruguiera species are limited. Here, the complete chloroplast sequences of five Bruguiera species were sequenced and assembled using Illumina data. The chloroplast genomes of B. gymnorhiza, B. hainesii, B. cylindrica, B. parviflora and B. sexangula were assembled into 161,195, 164,295, 164,297, 163,228 and 164,170 bp, respectively. All chloroplast genomes contain 37 tRNA and eight rRNA genes, with either 84 or 85 protein-coding genes. A comparative analysis of these genomes revealed high similarity in gene structure, gene order and boundary position of the LSC, SSC and two IR regions. Interestingly, B. gymnorhiza lost a rpl32 gene in the SSC region. In addition, a ndhF gene in B. parviflora straddles both the SSC and IRB boundary regions. These genes reveal differences in chloroplast evolution among Bruguiera species. Repeats and SSRs in the chloroplast genome sequences were found to be highly conserved between B. cylindrica and B. hainesii as well as B. gymnorhiza and B. sexangula indicating close genetic relationships based on maternal inheritance. Notably, B. hainesii, which is considered a hybrid between B. gymnorhiza and B. cylindrica, appears to have inherited the chloroplast from B. cylindrica. Investigating the effects of selection events on shared protein-coding genes showed a positive selection in rps7 and rpl36 genes in all species compared to land-plant species. A phylogenetic analysis, based on 59 conserved chloroplast protein-coding genes, showed strong support that all Bruguiera species are in the clade Rhizophoraceae. This study provides valuable genetic information for the study of evolutionary relationships and population genetics in Bruguiera and other mangrove species.

Plastome Evolution in the Hyperdiverse Genus Euphorbia (Euphorbiaceae) Using Phylogenomic and Comparative Analyses: Large-Scale Expansion and Contraction of the Inverted Repeat Region

With c. 2,000 species, Euphorbia is one of the largest angiosperm genera, yet a lack of chloroplast genome (plastome) resources impedes a better understanding of its evolution. In this study, we assembled and annotated 28 plastomes from Euphorbiaceae, of which 15 were newly sequenced. Phylogenomic and comparative analyses of 22 plastome sequences from all four recognized subgenera within Euphorbia revealed that plastome length in Euphorbia is labile, presenting a range of variation c. 42 kb. Large-scale expansions of the inverted repeat (IR) region were identified, and at the extreme opposite, the near-complete loss of the IR region (with only 355 bp left) was detected for the first time in Euphorbiaceae. Other structural variations, including gene inversion and duplication, and gene loss/pseudogenization, were also observed. We screened the most promising molecular markers from both intergenic and coding regions for phylogeny-based utilities, and estimated maximum likelihood and Bayesian phylogenies from four datasets including whole plastome sequences. The monophyly of Euphorbia is supported, and its four subgenera are recovered in a successive sister relationship. Our study constitutes the first comprehensive investigation on the plastome structural variation in Euphorbia and it provides resources for phylogenetic research in the genus, facilitating further studies on its taxonomy, evolution, and conservation.

Comparative Chloroplast Genomics of Sophora Species: Evolution and Phylogenetic Relationships in the Early-Diverging Legume Subfamily Papilionoideae (Fabaceae)

The taxonomy and evolutionary history of Sophora L., a genus with high economic and medicinal value, remain uncertain due to the absence of genetic resource (especially in China) and low polymorphism of molecular markers. Our aim was to elucidate the molecular evolution and phylogenetic relationships in chloroplast genomes of Sophora species in the early-diverging legume subfamily Papilionoideae (Fabaceae). We reported nine Sophora chloroplast genome from China using Illumina sequencing. We performed a series of analyses with previously published genomes of Sophora species to investigate their genomic characteristics, identified simple sequence repeats, large repeat sequences, tandem repeats, and highly polymorphic loci. The genomes were 152,953-158,087 bp in length, and contained 111-113 unique genes, including 76-78 protein coding, 31 tRNA, and 4 rRNA. The expansion of inverted repeat boundary of Sophora resulted in rps12 entering into the LSC region and loss of trnT-CGU gene in some species. Also, we found an approximately 23 kb inversion between trnC-GCA and trnF-GAA within the genus. In addition, we identified seven highly polymorphic loci (pi (π) > 0.035) suitable for inferring the phylogeny of Sophora species. Among these, three regions also co-occurred with large repeat sequences and support use of repeats as a proxy for the identification of polymorphic loci. Based on whole chloroplast genome and protein-coding sequences data-set, a well-supported phylogenetic tree of Sophora and related taxa showed that this genus is monophyletic, but sect. Disamaea and sect. Sophora, are incongruent with traditional taxonomic classifications based on fruit morphology. Our finding provides significant genetic resources to support further investigation into the phylogenetic relationship and evolution of the genus Sophora.