Comparative Analysis of Chloroplast Genomes within Saxifraga (Saxifragaceae) Takes Insights into Their Genomic Evolution and Adaption to the High-Elevation Environment

Comprehensive comparative analysis and development of molecular markers for Lasianthus species based on complete chloroplast genome sequences

BACKGROUND

Lasianthus species are widely used in traditional Chinese folk medicine with high medicinal value. However, source materials and herbarium specimens are often misidentified due to morphological characteristics and commonly used DNA barcode fragments are not sufficient for accurately identifying Lasianthus species. To improve the molecular methods for distinguishing among Lasianthus species, we report the complete chloroplast (CP) genomes of Lasianthus attenuatus, Lasianthus henryi, Lasianthus hookeri, Lasianthus sikkimensis, obtained via high-throughput Illumina sequencing.

RESULTS

These showed CP genomes size of 160164-160246 bp and a typical quadripartite structure, including a large single-copy region (86675-86848 bp), a small single-copy region (17177-17326 bp), and a pair of inverted repeats (28089-28135 bp). As a whole, the gene order, GC content and IR/SC boundary structure were remarkably similar among of the four Lasianthus CP genomes, the partial gene length and IR, LSC and SSC regions length are still different. The average GC content of the CP genomes was 36.71-36.75%, and a total of 129 genes were detected, including 83 different protein-coding genes, 8 different rRNA genes and 38 different tRNA genes. Furthermore, we compared our 4 complete CP genomes data with publicly available CP genome data from six other Lasianthus species, and we initially screened eleven highly variable region fragments were initially screened. We then evaluated the identification efficiency of eleven highly variable region fragments and 5 regular barcode fragments. Ultimately, we found that the optimal combination fragment' ITS2 + psaI-ycf4' could authenticated the Lasianthus species well. Additionally, the results of genome comparison of Rubiaceae species showed that the coding region is more conservative than the non-coding region, and the ycf1 gene shows the most significant variation. Finally, 49 species of CP genome sequences belonging to 16 genera of the Rubiaceae family were used to construct phylogenetic trees.

CONCLUSIONS

Our research is the first to analyze the chloroplast genomes of four species of Lasianthus in detail and we ultimately determined that the combination fragment' ITS2 + psaI-ycf4' is the optimal barcode combination for identifying the genus of Lasianthus. Meanwhile, we gathered the available CP genome sequences from the Rubiaceae and used them to construct the most comprehensive phylogenetic tree for the Rubiaceae family. These investigations provide an important reference point for further studies in the species identification, genetic diversity, and phylogenetic analyses of Rubiaceae species.

Factors contributing to organelle genomes size variation and the intracellular DNA transfer in Polygonaceae

The use of complete organelle genomes, including chloroplast and mitochondrial genomes, is a powerful molecular method for studying biological evolution and gene transfer. However, in the case of Polygonaceae, an important family with numerous edible, medicinal, and ornamental species, the mitochondrial genomes of only three species have been sequenced and analyzed. In this study, we present the mitochondrial and chloroplast genomes of two important Tibetan medicinal plants, Bistorta viviparum and B. macrophyllum. All the organelle genomes are assembled into a single circular structure and contain a common set of 32 protein-coding genes (PCGs). Some genes such as rps2 and ndhF were found to have high nucleotide polymorphism (Pi) in the chloroplast genomes, while cox1, mttB and rps12 showed pronounced Pi values in the mitochondrial genomes. Furthermore, our analysis revealed that most chloroplast genes and mitochondrial PCGs in Polygonaceae plants are under purifying selection. However, a few genes, including the chloroplast gene psaJ and the mitochondrial genes ccmFc, atp8 and nad4, showed positive selection in certain Polygonaceae plants, as indicated by a Ka/Ks ratio greater than one. Structural variation analysis revealed a wealth of differences between the mitochondrial genomes of five Polygonaceae species, with a particularly notable large-scale inversion observed between Reynoutria japonica and Fallopia aubertii. Furthermore, an analysis of the homologous sequences in the chloroplast and mitochondrial genomes revealed that the rps7 has been transferred from the chloroplast to the mitochondrial genome in all five Polygonaceae species. Finally, ecological niche models were constructed for B. viviparum and B. macrophyllum, indicating that mean annual temperature and altitude are the main climatic factors influencing the distribution of both species. Although the current distribution of B. viviparum is significantly wider than that of B. macrophyllum, projections suggest that the optimal growth ranges of both species will expand in the future, with B. macrophyllum potentially exceeding B. viviparum. This study not only contributes to the plastid genome database for Polygonaceae plants, but also provides theoretical insights into the adaptive evolution of these plants.

Comparative analysis of the complete chloroplast genomes of thirteen Bougainvillea cultivars from South China with implications for their genome structures and phylogenetic relationships

Bougainvillea spp., belonging to the Nyctaginaceae family, have high economic and horticultural value in South China. Despite the high similarity in terms of leaf appearance and hybridization among Bougainvillea species, especially Bougainvillea × buttiana, their phylogenetic relationships are very complicated and controversial. In this study, we sequenced, assembled and analyzed thirteen complete chloroplast genomes of Bougainvillea cultivars from South China, including ten B. × buttiana cultivars and three other Bougainvillea cultivars, and identified their phylogenetic relationships within the Bougainvillea genus and other species of the Nyctaginaceae family for the first time. These 13 chloroplast genomes had typical quadripartite structures, comprising a large single-copy (LSC) region (85,169-85,695 bp), a small single-copy (SSC) region (18,050-21,789 bp), and a pair of inverted-repeat (IR) regions (25,377-25,426 bp). These genomes each contained 112 different genes, including 79 protein-coding genes, 29 tRNAs and 4 rRNAs. The gene content, codon usage, simple sequence repeats (SSRs), and long repeats were essentially conserved among these 13 genomes. Single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were detected among these 13 genomes. Four divergent regions, namely, trnH-GUG_psbA, trnS-GCU_trnG-UCC-exon1, trnS-GGA_rps4, and ccsA_ndhD, were identified from the comparative analysis of 16 Bougainvillea cultivar genomes. Among the 46 chloroplast genomes of the Nyctaginaceae family, nine genes, namely, rps12, rbcL, ndhF, rpoB, rpoC2, ndhI, psbT, ycf2, and ycf3, were found to be under positive selection at the amino acid site level. Phylogenetic relationships within the Bougainvillea genus and other species of the Nyctaginaceae family based on complete chloroplast genomes and protein-coding genes revealed that the Bougainvillea genus was a sister to the Belemia genus with strong support and that 35 Bougainvillea individuals were divided into 4 strongly supported clades, namely, Clades Ⅰ, Ⅱ, Ⅲ and Ⅳ. Clade Ⅰ included 6 individuals, which contained 2 cultivars, namely, B. × buttiana 'Gautama's Red' and B. spectabilis 'Flame'. Clades Ⅱ only contained Bougainvillea spinosa. Clade Ⅲ comprised 7 individuals of wild species. Clade Ⅳ included 21 individuals and contained 11 cultivars, namely, B. × buttiana 'Mahara', B. × buttiana 'California Gold', B. × buttiana 'Double Salmon', B. × buttiana 'Double Yellow', B. × buttiana 'Los Banos Beauty', B. × buttiana 'Big Chitra', B. × buttiana 'San Diego Red', B. × buttiana 'Barbara Karst', B. glabra 'White Stripe', B. spectabilis 'Splendens' and B. × buttiana 'Miss Manila' sp. 1. In conclusion, this study not only provided valuable genome resources but also helped to identify Bougainvillea cultivars and understand the chloroplast genome evolution of the Nyctaginaceae family.

Phylogenomic analysis and dynamic evolution of chloroplast genomes of Clematis nannophylla

Clematis nannophylla is a perennial shrub of Clematis with ecological, ornamental, and medicinal value, distributed in the arid and semi-arid areas of northwest China. This study successfully determined the chloroplast (cp) genome of C. nannophylla, reconstructing a phylogenetic tree of Clematis. This cp genome is 159,801 bp in length and has a typical tetrad structure, including a large single-copy, a small single-copy, and a pair of reverse repeats (IRa and IRb). It contains 133 unique genes, including 89 protein-coding, 36 tRNA, and 8 rRNA genes. Additionally, 66 simple repeat sequences, 50 dispersed repeats, and 24 tandem repeats were found; many of the dispersed and tandem repeats were between 20-30 bp and 10-20 bp, respectively, and the abundant repeats were located in the large single copy region. The cp genome was relatively conserved, especially in the IR region, where no inversion or rearrangement was observed, further revealing that the coding regions were more conserved than the noncoding regions. Phylogenetic analysis showed that C. nannophylla is more closely related to C. fruticosa and C. songorica. Our analysis provides reference data for molecular marker development, phylogenetic analysis, population studies, and cp genome processes to better utilise C. nannophylla.

Molecular and morphological evidence supports the resurrection of Chrysospleniumguangxiense H.G.Ye & Gui C.Zhang (Saxifragaceae)

Chrysospleniumguangxiense H.G.Ye & Gui C.Zhang was first described as a new species in 1994 but later synonymized in the Flora of China treatment with C.glossophyllum H.Hara. Plastid genomes and nrDNA sequences were used to infer the phylogenetic relationships of selected taxa in Chrysosplenium. Our phylogenetic analyses revealed that C.guangxiense belongs to sect. Alternifolia, is closely related to Chrysospleniumhydrocotylifolium H.Lév. & Vaniot but distant from C.glossophyllum. Morphologically, C.guangxiense could be easily distinguished from C.glossophyllum by having robust rhizomes, basal leaves with a long cuneate base and fewer teeth in the margin, curled sepal margins, and red, larger seeds. It could also be easily distinguished from C.hydrocotylifolium by possessing long elliptic leaves and a long cuneate leaf base. Along with the phylogenetic studies, the complete plastid genome of C.guangxiense was also reported. The plastid genome was 154,004 bp in length and comprised two inverted repeats (IRs) of 28,120 bp, separated by a large single-copy of 80,646 bp and a small single-copy of 17,118 bp. A total of 111 functional genes were discovered, comprising 78 protein-coding genes, 29 tRNA genes, and four rRNA genes. Based on assessment of morphological and molecular data Chrysospleniumguangxiense H.G.Ye & Gui C.Zhang is resurrected from C.glossophyllum H.Hara at species level. A global conservation assessment classifies C.guangxiense as Vulnerable (VU).

The complete chloroplast genome of Meconopsis torquata (Papaveraceae), a traditional Tibetan medicine

Meconopsis torquata Prain 1906, a national second-class rare and endangered plant, is reported here for the first time for its complete chloroplast genome. The genome is 153,290 bp in length, comprising a large single-copy region (LSC, 83,918 bp), a small single-copy region (SSC, 17,740 bp), and two inverted repeat sequences (IRa and IRb, each 25,816 bp). The overall GC content is 38.7%, with the IR region having the highest content (43.1%). The genome is annotated with 112 unique genes, including 4 rRNA genes, 29 tRNA genes, and 79 protein-coding genes. Analysis of codon usage bias reveals that codons ending in A/T account for 96.7% of those with a Relative Synonymous Codon Usage (RSCU) value above 1. This predominance of A/T-ending codons might be indicative of M. torquata adaptation to high-altitude environments. Phylogenetic analysis reveals a close kinship between M. torquata and M. pinnatifolia and M. paniculata, indicating that the ancestral groups of these species might have a complex evolutionary history. This study uncovers the genetic characteristics and adaptive evolution of M. torquata, offering a new perspective in understanding the phylogenetic relationships within the genus. The findings not only provide a solid theoretical foundation for the conservation and sustainable use of this rare and endangered species but also offer significant scientific support for the conservation of biodiversity.

Chloroplast genome data of five Amygdalus species: Clarifying genome structure and phylogenetic relationships

Amygdalus species have considerable ecological and economic value, however, the phylogenetic relationships among Amygdalus remain controversy. In this study, we sequenced and assembled the chloroplast (cp) genomes of five Amygdalus species: Prunus communis, P. mongolica, P. pedunculata, P. triloba, and P. mira. We then conducted comparative genomic analyses and constructed their phylogenetic relationships. The genome length ranged from 157,870 to 158,451 bp, and 131 genes were annotated (86 protein-coding genes, 37 tRNAs, and 8 rRNAs). Additionally, 49-57 simple sequence repeats were detected, with most in the large single-copy region and with AT base preferences. Comparative genomic analyses revealed high similarities in structure, order, and gene content. However, we identified four highly divergent sequences: trnR-UCU-atpA, nbdhC-trnV-UAC, ycf4-cemA, and rpl32-trnL-UAG. The phylogenomic relationship analysis suggested that the Amygdalus species were grouped together, in which P. pedunculata, P. triloba, and Prunus tangutica were categorized into a branch, P. mongolica and Prunus davidiana were clustered a branch. This study provides an improved understanding of the genetic relationships among the Amygdalus and provides a basis for the development and utilization of Amygdalus resources.

Characteristics of the complete chloroplast genome of Saxifragaceae species Bergenia purpurascens (Hook. f. et Thoms.) Engl

Bergenia purpurascens (Hook. f. et Thoms.) Engl. is one species of traditional Chinese medicinal plant. This is the first publication of its complete chloroplast (cp) genome. The whole cp genome has 157,246 base pairs in length with 132 annotated genes, of which were 87 protein-coding genes, 37 tRNAs, and 8 rRNAs. According to the phylogenetic study, B. purpurascens and Bergenia scopulosa T. P. Wang. 1974 had a sister relationship. This genomic data and conclusions from B. purpurascens phylogenetic research will provide useful information and throw light on more in-depth investigations of the systematics and evolutionary patterns of Saxifragaceae.

Complete chloroplast genome assembly and phylogenetic analysis of blackcurrant (Ribes nigrum), red and white currant (Ribes rubrum), and gooseberry (Ribes uva-crispa) provide new insights into the phylogeny of Grossulariaceae

Background

Blackcurrant (Ribes nigrum), red currant (R. rubrum), white currant (R. rubrum), and gooseberry (R. uva-crispa) belong to Grossulariaceae and are popular small-berry crops worldwide. The lack of genomic data has severely limited their systematic classification and molecular breeding.

Methods

The complete chloroplast (cp) genomes of these four taxa were assembled for the first time using MGI-DNBSEQ reads, and their genome structures, repeat elements and protein-coding genes were annotated. By genomic comparison of the present four and previous released five Ribes cp genomes, the genomic variations were identified. By phylogenetic analysis based on maximum-likelihood and Bayesian methods, the phylogeny of Grossulariaceae and the infrageneric relationships of the Ribes were revealed.

Results

The four cp genomes have lengths ranging from 157,450 to 157,802 bp and 131 shared genes. A total of 3,322 SNPs and 485 Indels were identified from the nine released Ribes cp genomes. Red currant and white currant have 100% identical cp genomes partially supporting the hypothesis that white currant (R. rubrum) is a fruit color variant of red currant (R. rubrum). The most polymorphic genic and intergenic region is ycf1 and trnT-psbD, respectively. The phylogenetic analysis demonstrated the monophyly of Grossulariaceae in Saxifragales and the paraphyletic relationship between Saxifragaceae and Grossulariaceae. Notably, the Grossularia subgenus is well nested within the Ribes subgenus and shows a paraphyletic relationship with the co-ancestor of Calobotrya and Coreosma sections, which challenges the dichotomous subclassification of the Ribes genus based on morphology (subgenus Ribes and subgenus Grossularia). These data, results, and insights lay a foundation for the phylogenetic research and breeding of Ribes species.

Codon Usage Analyses Reveal the Evolutionary Patterns among Plastid Genes of Saxifragales at a Larger-Sampling Scale

Saxifragales is a 15-family order of early-divergent Eudicots with a rich morphological diversity and an ancient rapid radiation. Codon usage bias (CUB) analyses have emerged as an essential tool for understanding the evolutionary dynamics in genes. Thus far, the codon utilization patterns had only been reported in four separate genera within Saxifragales. This study provides a comprehensive assessment of the codon manipulation based on 50 plastid genes, covering 11 constituent families at a larger sampling scale. Our results first showed a high preference for AT bases and AT-ending codons. We then used effective number of codons (ENC) to assess a range of codon bias levels in the plastid genes. We also detected high-informative intrafamilial differences of ENC in three families. Subsequently, parity rule 2 (PR2) plot analyses revealed both family-unique and order-shared bias patterns. Most importantly, the ENC plots and neutrality analyses collectively supported the dominant roles of selection in the CUB of Saxifragales plastid genes. Notably, the phylogenetic affinities inferred by both ML and BI methods were consistent with each other, and they all comprised two primary clades and four subclades. These findings significantly enhance our understanding of the evolutionary processes of the Saxifrage order, and could potentially inspire more CUB analyses at higher taxonomic levels.