CPGAVAS2, an integrated plastome sequence annotator and analyzer

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.

Complete chloroplast genome sequence of Karelinia caspia (Pall.) Less. (Compositae)

Karelinia caspia (Compositae) is a perennial herbaceous plant owning high economic, feeding and medicinal values. It is widely distributed in desertification and saline alkali areas. The complete chloroplast genome was firstly reported in this study. The chloroplast genome of K. caspia with a total size of 151,239 bp consists of two inverted repeats separated by a large single-copy region and a small single-copy region. Its chloroplast genome contains 129 genes, including 85 protein-coding genes, 36 tRNA genes, and 8 rRNA genes. Also, a total of 62 simple sequence repeats were identified. These results will be useful for study on the evolution and genetic diversity of K. caspia in the future.

Saito M, Iwasaki T. The complete chloroplast genome of Viola vaginata (Violaceae), an endemic species of the snowy region in Japan

Viola vaginata, a perennial herb in subsection Stolonosae, is endemic to the snowy mountainous regions on the Sea of Japan side of Japan. Its complete chloroplast genome was 156,056 bp in length, comprising one large single-copy region (86,407 bp), one small single-copy region (17,301 bp), and two inverted repeat regions (27,174 bp each). It contained 111 unique genes, including 77 protein-coding genes, 30 transfer RNA genes, and 4 ribosomal RNA genes. Phylogenetic analysis placed V. vaginata in a clade with subsection Biobatae species and some Patellares species, while other Patellares species formed a distinct clade, contrasting with previous nuclear ITS results. These findings highlight the phylogenetic complexity within Viola.

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.

Characterization of the complete mitochondrial genome of the rice bean (Vigna umbellata)

BACKGROUND

Rice bean (Vigna umbellata), an underrated legume crop, demonstrates strong adaptability to poor soil fertility and has significant potential to enhance global food security. It is valuable both as a vegetable and fodder crop due to its high protein content, essential fatty acids, and micronutrients. Despite the sequencing of a high-quality genome of rice bean, its mitochondrial genome (mitogenome) sequence has not yet been reported.

RESULTS

For the first time, the rice bean mitogenome was assembled and annotated using PacBio HiFi sequencing and Geseq software. The mitogenome is a circular molecule with a length of 404,493 bp, containing 32 protein-coding genes, 17 tRNAs, and 3 rRNAs. Codon usage and sequence repeats were also determined. Six gene migration events from the chloroplast to the mitogenome were detected in rice bean. A phylogenetic analysis, including the rice bean mitogenome and 25 other taxa (23 of which are Fabales species), clarified the evolutionary and taxonomic status of rice bean. Additionally, a collinearity analysis of seven Fabales mitogenomes revealed high structural variability. In total, 473 RNA editing sites in protein-coding genes were identified.

CONCLUSIONS

This study presents the first sequencing, assembly, annotation, and analysis of the rice bean mitogenome, providing valuable background information for understanding the evolution of this species. These findings lay the groundwork for future genetic studies and molecular breeding efforts aimed at improving rice bean.

Exploring the phylogenetic framework and trait evolution of Impatiens through chloroplast genome analysis

BACKGROUND

The genus Impatiens, which includes both annual and perennial herbs, holds considerable ornamental, economic, and medicinal value. However, it posed significant challenges for taxonomic and systematic reconstruction. This was largely attributed to its high intraspecific diversity and low interspecific variation in morphological characteristics. In this study, we sequenced samples from 12 Impatiens species native to China and assessed their phylogenetic resolution using the complete chloroplast genome, in conjunction with published samples of Impatiens. In addition, a comparative analysis of chloroplast genomes were conducted to explore the evolution of the chloroplast genome in Impatiens.

RESULTS

The chloroplast genomes of 12 Impatiens species exhibited high similarity to previously published samples in terms of genome size, gene content, and sequence. The chloroplast genome of Impatiens exhibited a typical four-part structure, with lengths ranging from 146,987 bp(I. morsei)- 152,872 bp(I. jinpingensis). Our results identified 10 mutant hotspot regions (rps16, rps16-trnG, trnS-trnR, and rpoB-trnC) that could serve as effective molecular markers for phylogenetic analyses and species identification within the Impatiens. Phylogenetic analyses supported the classification of Impatiens as a monophyletic taxon. The identified affinities supported the taxonomic classification of the subgenus Clavicarpa within the Impatiens, with subgenus Clavicarpa being the first taxon to diverge. In phylogenetic tree,the Impatiens was divided into eight distinct clades. The results of ancestral trait reconstruction suggested that the ancestral traits of Impatiens included a perennial life cycle, four sepals and three pollen grooves. However, the ancestral morphology regarding fruit shape, flower colour, and spacing length remained ambiguous.

CONCLUDE

Our study largely supported the family-level taxonomic treatment of Impatiens species in China and demonstrated the utility of whole chloroplast genome sequences for phylogenetic resolution. Comparative analysis of the chloroplast genomes of Impatiens facilitated the development of molecular markers.The results of ancestral trait reconstruction showed that the ancestor type of habit was perennial, the number of sepals was 4, and morphology and number of aperture was 3 colpus. The traits of capsule shape, flower colour, and spur length underwent a complex evolutionary process. Our results provided data support for further studies and some important new insights into the evolution of the Impatiens.

Intraspecific Chloroplast Genome Genetic Polymorphism of Pinellia ternata (Xi Junecry) and Its Revelation of a Single Origin in Phylogeny

Background: Xi Junecry (Pinellia ternata), a perennial herb of the Araceae family, is indigenous to Xinxian County, Henan Province, China, and is regarded as a premium variety among similar medicinal materials. However, the lack of comprehensive genetic information on Xi Junecry germplasm resources has constrained the cultivation and identification of high-quality varieties. Methods: In this study, six chloroplast genomes of Xi Junecry were assembled and annotated using high-throughput sequencing. Subsequently, comparative analyses were conducted, and a phylogenetic tree was constructed. Results: The six Xi Junecry chloroplast genome lengths ranged from 157,456 to 158,406 bp, and the GC content was between 36.0% and 36.2%. A total of 265 single nucleotide polymorphism sites were identified across the six genomes, with a whole-genome nucleotide diversity (Pi) value of 0.00084. Among the four genomic regions, the small single-copy region exhibited the highest Pi, followed by the large single-copy region, while the inverted repeat region showed the lowest. Nucleotide polymorphism in coding regions was significantly lower than in non-coding regions. Nine hypervariable regions were identified, as follows: ndhE-ndhG, trnN-GUU-ndhF, trnS-GCU-trnG-UCC, atpB-rbcL, psaI, accD-ycf4, psbE-petL, psaC-ndhE, and psbI-trnG-UCC. Positive selection sites were detected in the accD and rbcL genes. Phylogenetic analysis clustered the six Xi Junecry samples into a distinct clade, separating them from other regional Pinellia samples. Conclusions: These findings elucidate the genetic variation levels in Xi Junecry and provide high-variability loci for population history inference, genetic diversity assessment, species domestication studies, and new cultivar development.

Analysis of the complete mitochondrial genome of Panax quinquefolius reveals shifts from cis-splicing to trans-splicing of intron cox2i373

Panax quinquefolius is a perennial plant with medicinal values. In this study, we assembled the complete mitochondrial genome (mitogenome) of P. quinquefolius using PMAT assembler. The total length of P. quinquefolius mitogenome is 573,154 bp. We annotated a total of 34 protein-coding genes (PCGs), 35 tRNA genes, and 6 rRNA genes in this mitogenome. The analysis of repetitive elements shows that there are 153 SSRs, 24 tandem repeats and 242 pairs of dispersed repeats this mitogenome. Also, we found 24 homologous sequences with a total length of 64,070 bp among its mitogenome and plastome, accounting for 41.05 % of the plastome, and 11.18 % of the mitogenome, showing a remarkable frequent sequence dialogue between plastome and mitogenomes. Besides, a total of 583 C to U RNA editing sites on 34 PCGs of high confidence were predicted by using Deepred-mt. We also inferred the phylogenetic relationships of P. quinquefolius and other angiosperms based on mitochondrial PCGs. Finally, we observed a shift from cis- to trans-splicing in P. quinquefolius for two mitochondrial introns, namely cox2i373 and nad1i728, and a pair of 48 bp short repetitive sequences may be associated with the breaking and rearrangement of the cox2i373 intron. The fragmentation of the cox2i373 intron was further confirmed by our PCR amplification experiments. In summary, our report on the P. quinquefolius mitogenome provides a new perspective on the intron evolution of the mitogenome.

Revealing the Dynamic History of Parasitic Plant Plastomes via Structural Characterization, Comparative Analysis, and Phylogenomics

Background: The shift to a parasitic lifestyle in plants often leaves distinct marks on their plastid genomes, given the central role plastids play in photosynthesis. Studying these unique adaptations in parasitic plants is essential for understanding the mechanisms and evolutionary patterns driving plastome reduction in angiosperms. By exploring these changes, we can gain deeper insights into how parasitism reshapes the genomic architecture of plants. Method: This study analyzed and compared the plastomes of 113 parasitic plants from different families. Results: The Orobanchaceae family (hemiparasitic plants) displayed the largest plastome size, while Apodanthaceae exhibited the shortest. Additionally, Orobanchaceae showcased little to no gene loss in their plastomes. However, holoparasitic species typically had reduced plastome sizes. Convolvulaceae exhibited significantly reduced plastome sizes due to high gene loss, and Apodanthaceae retained only a few genes. Gene divergence among different families was also investigated, and rps15, rps18, and rpl33 in Orobanchaceae; accD and ycf1 in Convolvulaceae; atpF and ccsA in Loranthaceae; and rpl32 in Santalaceae showed greater divergence. Additionally, Orobanchaceae had the highest numbers of all repeat types, whereas Loranthaceae and Convolvulaceae exhibited the lowest repeat numbers. Similarly, more simple sequence repeats were reported in Loranthaceae and Santalaceae. Our phylogenetic analysis also uncovered a distinct clade comprising Loranthaceae, with a single Schoepfiaceae species clustering nearby. Contrary to expectations, parasitic and hemiparasitic plants formed mixed groupings instead of segregating into separate clades. Conclusions: These findings offer insights into parasitic plants' evolutionary relationships, revealing shared and divergent genomic features across diverse lineages.

Complete chloroplast genome sequence of the medicinal plant Oxyceros horridus (Rubiaceae) and phylogenetic analysis

Oxyceros horridus Lour. 1790, valued for ornamental and medicinal properties, has been extensively utilized in traditional medicinal in Vietnam. This study assembled and characterized the first chloroplast of Oxyceros horridus, consisted of 155,333 bp with a GC content of 37.5%. The genome included 130 annotated genes (including 85 protein-coding genes, 37 tRNA genes, and 8 rRNA genes), exhibited a typical quadripartite structure consistent with those found in other members of the Rubiaceae family. Phylogenetic analysis placed O. horridus in Ixoroideae subfamily, closely related to Fosergia shweliensis and Gardenia tenophylla. This study enriches cp genome data and phylogenetic insights within Rubiaceae.

The complete chloroplast genome sequence of Melothria scabra (Cucurbitaceae)

Melothria scabra has gradually become an economically important plant worldwide. The complete chloroplast genome of M. scabra has a length of 156,744 bp, contains a large single-copy (LSC) region (86,387 bp), a small single-copy (SSC) region (18,055 bp), and two inverted repeats (IRs) with the same length of 26,151 bp. In total, 126 genes were detected, including 83 protein-encoding genes, 35 transfer RNA (tRNA) genes, and eight ribosomal RNA (rRNA) genes. For phylogenetic analysis, M. scabra has a closer genetic relationship with Cucumis sativus and Citrullus lanatus. The complete chloroplast genome sequence of M. scabra would promote the germplasm exploration, phylogenetic relationships, and molecular biology researches in Melothria.

Comprehensive analysis of the first complete mitogenome and plastome of a traditional Chinese medicine Viola diffusa

BACKGROUND

Viola diffusa is used in the formulation of various Traditional Chinese Medicines (TCMs), including antiviral, antimicrobial, antitussive, and anti-inflammatory drugs, due to its richness in flavonoids and triterpenoids. The biosynthesis of these compounds is largely mediated by cytochrome P450 enzymes, which are primarily located in the membranes of mitochondria and the endoplasmic reticulum.

RESULTS

This study presents the complete assembly of the mitogenome and plastome of Viola diffusa. The circular mitogenome spans 474,721 bp with a GC content of 44.17% and encodes 36 unique protein-coding genes, 21 tRNA, and 3 rRNA. Except for the RSCU values of 1 observed for the start codon (AUG) and tryptophan (UGG), the mitochondrial protein-coding genes exhibited a codon usage bias, with most estimates deviating from 1, similar to patterns observed in closely related species. Analysis of repetitive sequences in the mitogenome demonstrated potential homologous recombination mediated by these repeats. Sequence transfer analysis revealed 24 homologous sequences shared between the mitogenome and plastome, including nine full-length genes. Collinearity was observed among Viola diffusa species within the other members of Malpighiales order, indicated by the presence of homologous fragments. The length and arrangement of collinear blocks varied, and the mitogenome exhibited a high frequency of gene rearrangement.

CONCLUSIONS

We present the first complete assembly of the mitogenome and plastome of Viola diffusa, highlighting its implications for pharmacological, evolutionary, and taxonomic studies. Our research underscores the multifaceted importance of comprehensive mitogenome analysis.

Cladopusyangjiangensis (Podostemaceae), a new species from Guangdong, South China, redefining the phylogenetic relationships within Cladopus

This paper introduces Cladopusyangjiangensis, a newly identified species that enriches our understanding of the diversity of the Podostemaceae in East Asia. Distinctive in its morphological traits, this species is characterized by the region's longest flowering shoots and exhibits a high number of elongated leaves per cluster, along with relatively slender roots. Phylogenetic analyses using Maximum Likelihood and Bayesian Inference methods on plastome and matK sequences confirm C.yangjiangensis as a distinct species. It forms a clade with C.fukienensis, its closest relative, together branching off from C.austrosinensis. The plastome of C.yangjiangensis is 132,818 bp in length, comprising two inverted repeat regions of 20,881 bp, which are separated by large and small single-copy regions of 78,713 and 12,343 bp, respectively. Genetic analysis reveals the extensive loss of the ycf1 and ycf2 genes in the chloroplast genome, a trait common to the Podostemaceae, suggesting adaptations to environmental conditions or gene transfers to nuclear or mitochondrial genomes. This study improves the clarity of phylogenetic relationships in previous studies and underscores the importance of continued taxonomic and phylogenetic research.

Phylogenomic analyses of the pantropical Platycerium Desv. (Platycerioideae) reveal their complex evolution and historical biogeography

Platycerium is a genus of pantropical epiphytic ferns consisting of ca. 18 species and are highly sought after by horticultural enthusiasts. Although the monophyly of this genus has been well supported in previous molecular studies, as an intercontinentally disjunct genus, the origin and distribution pattern of Platycerium were elusive and controversial. This is mainly due to limited taxon sampling, a plastid representing only a single coalescent history, the lack of fossil evidence, and so on. Here, by utilizing genome-skimming sequencing, transcriptome sequencing, and flow cytometry, we integrated chloroplast genomes, data of single-copy nuclear genes, ploidy levels, morphology, and geographic distribution to understand the species phylogeny and the evolutionary and biogeographic history of Platycerium. Our major results include: (1) based on both plastid and nuclear datasets, Platycerium is consistently resolved into three fully supported clades: the Afro-American (AA) clade, the Javan-Australian (JA) clade, and the Malayan-Asian (MA) clade. The AA clade and MA clade are further divided into three and two subclades, respectively; (2) a large amount of gene tree conflict, as well as cytonuclear discordance, was found and can be explained by hybridization and incomplete lineage sorting, and most of the hybridization hypotheses represented ancient hybridization events; (3) through molecular dating, the crown age of Platycerium is determined to be at approximately 32.79 Ma based on the plastid dataset or 29.08 Ma based on the nuclear dataset in the Middle Oligocene; (4) ancestral area reconstruction analysis from different datasets showed that Platycerium most likely originated from Indochina; (5) current distribution patterns are resultant from long-distance dispersals, ancient orogeny, and an ancient climate event; and (6) species diversification was driven by polyploidization, dispersal, and hybridization. This study presented here will help understand the evolution of tropical plant flora and provide a reference for the cultivation and breeding of staghorn ferns.

Mitochondrial Genome Assembly and Comparative Analysis of Chionanthus Retusus (Oleaceae)

Background/Objectives: Chionanthus retusus Lindl. & Paxton is an ornamental tree species native to North China. Research on the mitochondrial genome can elucidate the evolution and biological characteristics of C. retusus and better protect this important species. Methods and Results: This work aimed to clarify the evolutionary and phylogenetic links by sequencing, assembling, annotating, and analyzing the entire mitochondrial genome of C. retusus. The single-loop structure that made up the mitochondrial genome had a total length of 657,640 bp and a GC content of 44.52%. In total, 37 unique protein-coding genes, 20 tRNA genes, and 3 rRNA genes were identified. Numerous repeat sequences and migrating fragments of chloroplast sequences were found. Using the mitochondrial protein-coding genes to construct evolutionary trees, it was found that the closest relative of C. retusus is C. rupicola (Lingelsh.) Kiew. Conclusions: This research represents the first comprehensive set of data on the mitochondrial genome of an ancient (>500 yr) C. retusus specimen. In addition to elucidating the biological characteristics of C. retusus. The findings contribute to the Oleaceae mitochondrial genome database and offer valuable insights for future studies in molecular breeding, evolutionary biology, and genetic diversity conservation.

Structural and evolutionary analyses of the mitochondrial genome of Spuriopimpinella brachycarpa

Introduction

Spuriopimpinella brachycarpa (Kom.) Kitag., a member of the Apiaceae family, is a perennial aromatic herb native to Northeast Asia with applications in culinary and traditional medicine. Despite its significance, most studies on S. brachycarpa have primarily focused on its phytochemical properties, with limited insights into its molecular and genomic characteristics.

Methods

This study presents the sequencing and assembly of the mitochondrial genome (mitogenome) of S. brachycarpa using second- and third-generation high-throughput sequencing technologies. Comprehensive analyses were performed on its structural organization, RNA editing sites, relative synonymous codon usage (RSCU), and repeat sequences. Comparative analyses with closely related species were also conducted.

Results

The mitogenome exhibited a multi-branched structure, with a total length of 523,512 bp and a GC content of 43.37%. Annotation revealed 30 unique protein-coding genes, 21 tRNA genes, and three rRNA genes. Comparative analysis indicated that the S. brachycarpa mitogenome contains structural variations but shares collinear features with other Apiaceae species. We identified 618 potential RNA editing sites involving C-to-U conversions and discovered 59 homologous fragments between the mitogenome and plastome, comprising 8.13% of the mitogenome.

Discussion

These results enrich the genomic database of Apiaceae, providing valuable insights into the evolutionary relationships and genetic diversity within the family.

Chloroplast genome of four Amorphophallus species: genomic features,comparative analysis, and phylogenetic relationships among Amorphophallus species

BACKGROUND

The genus Amorphophallus (Araceae) contains approximately 250 species, most of which have high ecological and economic significance. The chloroplast genome data and the comprehensive analysis of the chloroplast genome structure of Amorphophallus is limited. In this study, four chloroplast genomes of Amorphophallus were sequenced and assembled. For the first time, comparative analyses of chloroplast genomes were conducted on the 13 Amorphophallus species in conjunction with nine published sequences.

RESULTS

The Amorphophallus chloroplast genomes exhibited typical quadripartite structures with lengths ranging from 164,417 to 177,076 bp. These structures consisted of a large single copy (LSC, 90,705 - 98,561 bp), a small single copy (SSC, 14,172 - 21,575 bp), and a pair of inverted repeats (IRs, 26,225 - 35,204 bp). The genomes contain 108 - 113 unique genes, including 76 - 79 protein-coding genes, 28 - 29 tRNA genes, and 4 rRNA genes. The molecular structure, gene order, content, codon usage, long repeats, and simple sequence repeats (SSRs) within Amorphophallus were generally conserved. However, several variations in intron loss and gene expansion on the IR-SSC boundary regions were found among these 13 genomes. Four mutational hotspot regions, including trnM-atpE, atpB, atpB-rbcL and ycf1 were identified. They could identify and phylogeny future species in the genus Amorphophallus. Positive selection was found for rpl36, ccsA, rpl16, rps4, rps8, rps11, rps12, rps14, clpP, rps3, ycf1, rpl20, rps2, rps18, rps19, atpA, atpF, rpl14, rpoA, rpoC1, rpoC2 and rps15 based on the analyses of Ka/Ks ratios. Phylogenetic inferences based on the complete chloroplast genomes revealed a sister relationship between Amorphophallus and Caladieae. All Amorphophallus species formed a monophyletic evolutionary clade and were divided into three groups, including CA-II, SEA, and CA-I. Amorphophallus albus, A. krausei, A. kachinensis and A. konjac were clustered into the CA-II clade, A. paeoniifolius and A. titanum were clustered into the SEA clade, A. muelleri 'zhuyajin1', Amorphophallus sp, A. coaetaneus, A. tonkinensis and A. yunnanensis were clustered into CA- I clade.

CONCLUSIONS

The genome structure and gene content of Amorphophallus chloroplast genomes are consistent across various species. In this study, the structural variation and comparative genome of chloroplast genomes of Amorphophallus were comprehensively analyzed for the first time. The results provide important genetic information for species classification, identification, molecular breeding, and evolutionary exploration of the genus Amorphophallus.

Complete Chloroplast Genomes and Phylogenetic Analysis of Woody Climbing Genus Phanera (Leguminosae)

BACKGROUND

Phanera Lour., a genus in the subfamily Cercidoideae of the family Leguminosae, is characterized by woody liana habit, tendrils, and distinctive bilobate or bifoliolate leaves. The genus holds important medicinal value and constitutes a complex group characterized by morphological diversity and unstable taxonomic boundaries. However, limited information on the chloroplast genomes of this genus currently available constrains our understanding of its species diversity. Hence, it is necessary to obtain more chloroplast genome information to uncover the genetic characteristics of this genus.

METHODS

We collected and assembled the complete chloroplast genomes of nine representative Phanera plants, including Phanera erythropoda, Phanera vahlii, Phanera aureifolia, Phanera bidentata, Phanera japonica, Phanera saigonensis, Phanera championii, Phanera yunnanensis, and Phanera apertilobata. We then conducted a comparative analysis of these genomes and constructed phylogenetic trees.

RESULTS

These species are each characterized by a typical quadripartite structure. A total of 130-135 genes were annotated, and the GC content ranged from 39.25-42.58%. Codon usage analysis indicated that codons encoding alanine were dominant. We found 82-126 simple sequence repeats, along with 5448 dispersed repeats, mostly in the form of forward repeats. Phylogenetic analysis revealed that 16 Phanera species form a well-supported monophyletic group, suggesting a possible monophyletic genus. Furthermore, 10 hypervariable regions were detected for identification and evolutionary studies.

CONCLUSIONS

We focused on comparing chloroplast genome characteristics among nine Phanera species and conducted phylogenetic analyses, laying the foundation for further phylogenetic research and species identification of Phanera.

The complete chloroplast genome of Erodium stephanianum (Geraniaceae)

Erodium stephanianum (Erodium stephanianum Willd. Sp. Pl. 1800) has not had its complete chloroplast genome reported, which limits our understanding of its genetics and evolution. In this study, we assembled and annotated its chloroplast genome, revealing a quadripartite structure with 76 protein-coding genes. Repeat analysis indicated the presence of simple sequence repeats. Phylogenetic analysis confirmed E. stephanianum's placement within the genus Erodium of the Geraniaceae family. These findings offer valuable genomic resources for comparative studies in Erodium and Geraniaceae, aiding genetic diversity and phylogenetic analyses.

Comparative analysis of chloroplast genomes and phylogenetic relationships in the endemic Chinese bamboo Gelidocalamus (Bambusoideae)

Introduction

Gelidocalamus Wen is a small yet taxonomically challenging genus within the Arundinarieae tribe. Recent molecular studies have suggested it may not be monophyletic. However, limited species sampling and insufficient molecular marker information have resulted in poorly resolved phylogenetic relationships within this genus.

Methods

The complete chloroplast genomes covering all 16 species and one variant of Gelidocalamus were sequenced, and comparative analyses were conducted. Phylogenetic analyses were performed using different molecular markers, including chloroplast data, the nuclear ribosomal DNA (nrDNA) repeats region, and 29 mitochondrial protein-coding genes. Additionally, the divergence times of Gelidocalamus were estimated to reveal their evolutionary history.

Results

The plastomes of Gelidocalamus ranged in size from 139,500 bp to 139,801 bp, with a total of 137 identified genes, including 90 protein-coding genes, 39 tRNA genes, and 8 rRNA genes. The size of the nrDNA repeats ranged from 5,802 bp to 5,804 bp. Phylogenetic analysis based on chloroplast data revealed that Gelidocalamus is polyphyletic, with different subclades distributed within the IV and V clades. However, phylogenetic analysis based on nrDNA and mitochondrial genes did not effectively resolve the relationships within the genus.

Discussion

Comparative analysis of chloroplast genomes indicated that Gelidocalamus shares a high degree of similarity with closely related genera in terms of chloroplast genome collinearity, codon usage bias, and repetitive sequences. Divergence time estimation suggests that it is a relatively young group, with all members appearing successively over the past four million years. The complex phylogenetic patterns may arise from the rapid radiation of Arundinarieae. This study provides a preliminary foundation for further in-depth research on the phylogeny, genomic structural features, and divergence times of this genus.

Comparative analysis of chloroplast genome of Lonicera japonica cv. Damaohua

Lonicera japonica is a well-known medicinal plant, and the Damaohua cultivar is one of the oldest known honeysuckle cultivars in China. The 155,151 bp chloroplast genome of this cultivar was obtained through Illumina sequencing. The genome includes a pair of inverted repeats (IRa and IRb; 23,789 bp each), a large single-copy region (88,924 bp), and a small single-copy (SSC) region (18,649 bp). In total, 127 unique genes were identified: 80 protein-coding, 39 tRNA, and 8 rRNA genes. Only ycf3 contained two introns. Eighty-nine large repetitive sequences and 54 simple sequence repeats were detected. Fifty potential RNA editing sites were predicted. Adaptive evolution analysis revealed that infA, matK, petB, petD, rbcL, rpl16, rpl2, rps3, ycf1, and ycf2 were positively selected, possibly reflecting the specific environmental adaptations of this cultivar. Sequence alignment and analysis revealed several candidate fragments for Lonicera species identification, such as the intergenic regions rpoB-petN, rbcL-accD, and psaA-ycf3. The IR region boundary and phylogenetic analysis revealed that the L. japonica cv. Damaohua chloroplast genome was most closely related to the L. japonica genome, but there were five distinct differences between the two. There are four sites with high variability between L. japonica and L. japonica cv. Damaohua with nucleotide variability (Pi) greater than 0.002, including rps2-rpoC2, atpB-rbcL, ycf1, and ycf1-trnN GUU. The differences between L. japonica and L. japonica cv. Damaohua were further confirmed by the single nucleotide polymorphism sites between these two species. Therefore, this study revealed that the chloroplast genome can serve as a universal super barcode for plant identification, which can identify differences and help distinguish Lonicera japonica from related species. An understanding of Lonicera japonica cv. Damaohua chloroplast genomics and a comparative analysis of Lonicera species will provide a scientific basis for breeding, species identification, systematic evolution analysis, and chloroplast genetic engineering research on medicinal honeysuckle plants.

The complete genome sequence of the chloroplast of Bidens aurea

Bidens aurea (Asteraceae), a native of tropical America is now widespread in Asia and the Americas. We explored the B. aurea chloroplast genome and conducted a phylogenetic analysis. The chloroplast genome was circular, consisting of a large single copy (LSC) of 83,909 base pairs (bp), a small single copy (SSC) of 18,407 bp, and two inverted repeat regions (IR) of 24,729 bp each. Phylogenetic analysis showed that the 19 Bidens taxa were divided into five major clades, and B. aurea was most closely related to two species. Our findings offer a high-quality B. aurea chloroplast genome, aiding DNA barcode development and evolutionary history reconstruction.

The complete plastid genome of Polygonatum gracile P. Y. Li (Asparagaceae): characterization and phylogeny

Polygonatum gracile P. Y. Li is a perennial rhizomatous medicinal herb endemic to China with recognized medicinal value. Here, the first complete plastid genome of P. gracile is reported and characterized, and the phylogenetic relationships of P. gracile with other Polygonatum species are clarified by phylogenetic inferences based on complete plastid genome sequences. The complete plastid genome of P. gracile is 155,734 bp in length and exhibits the typical quadripartite circular structure consisting of a large single-copy (LSC, 84,385 bp) region, a small single-copy (SSC, 18,519 bp) region, and a pair of inverted repeat regions (IRa/b, 26,415 bp). A total of 112 unique genes, including 78 protein-coding genes, 30 tRNA genes, and four rRNA genes were identified. Seven protein-coding genes, eight tRNA genes, and four rRNA genes were duplicated in the IR regions. Phylogenetic reconstructions indicated that P. gracile belongs to sect. Verticillata and is most closely related to P. zanlanscianense Pamp. This study provides a basis for further phylogenetic reconstructions, conservation and utilization of Polygonatum species.

Pan-Chloroplast Genomes Reveal the Accession-Specific Marker for Gastrodia elata f. glauca

Gastrodia elata rhizomes have been applied as traditional medicinal materials for thousands of years. In China, G. elata f. elata (red flower and stem, Ge), G. elata f. viridis (green, Gv), and G. elata f. glauca (black, Gg) represent the primary cultivars in artificial cultivation. Although the annual output of G. elata amounts to 150,000 tons, only 20% is Gg. The long production period, low yield, and high quality of Gg led to its extremely high market prices. However, an effective method to identify this crude drug based solely on its morphological or chemical characteristics is lacking. In this study, the complete chloroplast genomes of three G. elata variants were sequenced using the Illumina HiSeq 2500 platform. Another 21 chloroplast genomes from Gastrodia species, which have published in previous reports, were combined and analyzed together. Our results showed that larger genomic sizes, fewer long tandem repeats, and more simple sequence repeats were the major features of the Gg chloroplast genomes. Phylogenetic analysis showed that the Gg samples were separately clustered in a subclade. Moreover, an accession-specific marker was successfully developed and validated for distinguishing additional rhizome samples. Our study provides investigations of the taxonomic relationships of Gastrodia species. The molecular marker will be useful for differentiating Gastrodia products on the market.

Characterization and phylogenetic analysis of the chloroplast genome of Diospyros aff. oleifera

The Diospyros genus (Ebenaceae) has significant economic value. During field surveys, we discovered a Diospyros specimen showing morphological overlap with both D. oleifera and D. kaki var. silvestris, provisionally named Diospyros aff. oleifera. To resolve its taxonomy, we sequenced and analyzed its chloroplast genome. The complete chloroplast genome is 157,732 bp with a quadripartite structure. mVISTA analysis revealed unique sequence variations compared to related species. Phylogenetic analysis using 75 protein-coding genes grouped it with D. oleifera, indicating their close relationship. Our findings suggest this specimen likely represents a novel, undescribed species. This study provides insights into Diospyros diversity and a foundation for future research.

Comparative analysis of complete chloroplast genomes of 14 Asteraceae species

BACKGROUND

The Asteraceae family, the largest and one of the most diverse families of angiosperms, presents significant challenges in taxonomic classification and systematic research due to its vast species diversity and morphological complexity. A comprehensive understanding of the chloroplast genomes within this family is essential for refining taxonomic classifications and advancing phylogenetic studies.

METHODS AND RESULTS

In this study, we sequenced the complete chloroplast genomes of 14 Asteraceae species and conducted a thorough bioinformatic analysis of their characteristics. The chloroplast genomes, ranging from 150,907 bp to 152,858 bp, exhibit a typical quadripartite structure: a large single-copy (LSC) region (83,044 bp to 84,625 bp), a small single-copy (SSC) region (18,223 bp to 18,673 bp), and a pair of inverted repeats (IRs) (24,806 bp to 25,201 bp). These genomes encode 87 to 89 protein-coding genes (PCGs), 36 to 37 tRNA genes, and 8 rRNA genes, with high conservation in size, structure, gene content, and order. Comparative analysis with other Asteraceae species' chloroplast genomes revealed notable similarities and structural variations, particularly in the IR regions. Nucleotide polymorphism analysis identified four genes-trnY-GUA, trnE-UUC, ycf1, and rrn23-with higher Pi values, suggesting potential hotspots for evolutionary studies. Phylogenetic analysis using maximum likelihood (ML) and Bayesian inference (BI) approaches provided new insights, proposing the reclassification of Himalaiella auriculata and Jacobaea raphanifolia as independent genera, distinct from Saussurea and Senecio.

CONCLUSIONS

This study presents a comprehensive analysis of the chloroplast genome structures and phylogenetic relationships of 14 Asteraceae species, offering critical data for future molecular identification, evolutionary biology, and population genetics research. The findings hold significant implications for the ongoing refinement of Asteraceae taxonomic classifications and enhance our understanding of the evolutionary dynamics within this complex family.

Mitochondrial Genome Assembly and Structural Characteristics Analysis of Gentiana rigescens

Gentiana rigescens, an alpine plant with significant medicinal value, possesses a complex genetic background. However, comprehensive genomic research on G. rigescens is still lacking, particularly concerning its organelle genome. In this study, G. rigescens was studied to sequence the mitochondrial genome (mitogenome) and ascertain the assembly, informational content, and developmental expression of the mitogenome. The mitogenome of G. rigescens was 393,595 bp in length and comprised four circular chromosomes ranging in size from 6646 bp to 362,358 bp. The GC content was 43.73%. The mitogenome featured 30 distinct protein-coding genes, 26 tRNA genes, and 3 rRNA genes. The mitogenome of G. rigescens also revealed 70 SSRs, which were mostly tetra-nucleotides. In addition, 48 homologous fragments were found between the mitogenome and the chloroplast genome, with the longest measuring 23,330 bp. The documentation of the mitochondrial genome of G. rigescens is instrumental in advancing the understanding of its physiological development. Decoding the G. rigescens mitogenome will offer valuable genetic material for phylogenetic research on Gentianaceae and enhance the use of species germplasm resources.

The Vitextrifolia complex (Lamiaceae) in the Philippines

The Vitextrifolia complex in the Philippines comprises economically important medicinal species, but its taxonomic status has become ambiguous due to numerous historical nomenclatural revisions. We assembled the complete chloroplast genomes of five species belonging to this species complex to provide additional evidence for their species delimitation. Based on a detailed analysis of specimens which combined molecular and morphological data, we propose reinstating V.elmeri Moldenke and delineating V.arvensis Gentallan, Sengun & M.B.Bartolome as a new endemic species belonging to this complex. The new species is a putative hybrid between V.bicolor Willd. and V.elmeri Moldenke. The specific epithet arvensis reflects its predominantly cultivated nature, both on a commercial scale and in home gardens, as a valued medicinal plant. We also provided a key to identify the five species belonging to the V.trifolia complex in the Philippines.

Repeat-mediated recombination results in Complex DNA structure of the mitochondrial genome of Trachelospermum jasminoides

BACKGROUND

Trachelospermum jasminoides has medicinal and ornamental value and is widely distributed in China. Although the chloroplast genome has been documented, the mitochondrial genome has not yet been studied.

RESULTS

The mitochondrial genome of T. jasminoides was assembled and functionally annotated using Illumina and nanopore reads. The mitochondrial genome comprises a master circular molecular structure of 605,764 bp and encodes 65 genes: 39 protein-coding genes, 23 transfer RNA (tRNA) genes and 3 ribosomal RNA genes. In addition to the single circular conformation, we found many alternative conformations of the T. jasminoides mitochondrial genome mediated by 42 repetitive sequences. Six repetitive sequences (DRS01-DRS06) were supported by nanopore long reads, polymerase chain reaction (PCR) amplifications, and Sanger sequencing of the PCR products. Eleven homologous fragments were identified by comparing the mitochondrial and chloroplast genome sequences, including three complete tRNA genes. Moreover, 531 edited RNA sites were identified in the protein-coding sequences based on RNA sequencing data, with nad4 having the highest number of sites (54).

CONCLUSION

To our knowledge, this is the first description of the mitochondrial genome of T. jasminoides. Our results demonstrate the existence of multiple conformations. These findings lay a foundation for understanding the genetics and evolutionary dynamics of Apocynaceae.

Characterization of the complete chloroplast genome sequence of Artemisia sylvatica Maximowicz 1859 (Asteraceae)

Artemisia sylvatica Maximowicz 1859 is one of the medicinal herbs in Artemisia. This study presents the complete chloroplast genome of A. sylvatica, sequenced using the Illumina NovaSeq platform. The genome is 151,161 bp in length, featuring a GC content of 38%. It consists of a large single-copy (LSC) region of 82,892 bp, a small single-copy (SSC) region of 18,353 bp, and two inverted repeat (IR) regions of 24,958 bp each. In total, the genome contains 132 genes, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Phylogenetic analysis positions A. sylvatica within the subgenus Artemisia, highlighting its evolutionary relationships within this diverse genus. The first chloroplast genome of A. sylvatica was reported in this work contributes to the enrichment of genomic data for the genus Artemisia.

Deciphering the multi- partite mitochondrial genome of Crataegus pinnatifida: insights into the evolution and genetics of cultivated Hawthorn

Flowering plant (angiosperm) mitochondrial genomes are remarkably dynamic in their structures. We present the complete mitochondrial genome of hawthorn (Crataegus pinnatifida Bunge), a shrub that bears fruit and is celebrated for its extensive medicinal history. We successfully assembled the hawthorn mitogenome utilizing the PacBio long-read sequencing technique, which yielded 799,862 reads, and the Illumina novaseq6000 sequencing platform, which producing 6.6 million raw paired reads. The C. pinnatifida mitochondria sequences encompassed a total length of 440,295 bp with a GC content of 45.42%. The genome annotates 54 genes, including 34 that encode proteins, 17 that encode tRNA, and three genes for rRNA. A fascinating interplay was observed between the chloroplast and mitochondrial genomes, which share 17 homologous sequences sequences that rotal 1,933 bp. A total of 134 SSRs, 22 tandem repeats and 42 dispersed repeats were identified in the mitogenome. Four conformations of C. pinnatifida mitochondria sequences recombination were verified through PCR experiments and Sanger sequencing, and C. pinnatifida mitogenome is more likely to be assembled into three circular-mapping chromosomes. All the RNA editing sites that were identified C-U edits, which predominantly occurred at the first and second positions of the codons. Phylogenetic and collinearity analyses identified the evolutionary trajectory of C. pinnatifida, which reinforced the genetic identity of the hawthorn section. This unveiling of the unique multi-partite structure of the hawthorn mitogenome offers a foundational reference for future study into the evolution and genetics of C. pinnatifida.

The complete chloroplast genome of Calyptothecium philippinense Broth. (Pterobryaceae, Hypnales)

The genus Calyptothecium, currently comprising ca. 30 species worldwide, is the largest genus within the family Pterobryaceae. However, a comprehensive taxonomic revision of this genus is lacking. Calyptothecium philippinense Broth. 1899, a moss species widely found in the tropical regions of Asia, is characterized by the unique rugose leaves and large auriculate leaf bases. In this study, we sequenced the complete chloroplast genome (CPG) of C. philippinense using the Illumina NovaSeq 6000 platform. The length of the CPG of C. philippinense was determined to be 124,513 bp, with an AT content of 74%. The CPG of C. philippinense exhibited a standard quadripartite structure, consisting of one small single-copy (SSC) region (18,541 bp), one large single-copy region (LSC) (87,222 bp), and two inverted repeat (IR) regions (9375 bp each). A total of 126 genes from the CPG of C. philippinense were annotated, including 82 protein-coding genes, eight ribosomal RNA genes, and 36 transfer RNA genes. Phylogenetic analysis based on the CPGs of 25 bryophyte taxa revealed that the three Pterobryaceae species C. philippinense, Calyptothecium hookeri (Mitt.) Broth. and Pterobryopsis orientalis (Müll. Hal.) M. Fleisch. formed a robust clade. The findings could facilitate more accurate classification and help elucidate evolutionary relationships within Calyptothecium.

Assembly and comparative analysis of the first complete mitochondrial genome of Salix psammophila, a good windbreak and sand fixation shrub

Salix psammophila, commonly known as the sandlive willow, is a vital shrub species within the Salicaceae family, particularly significant for its ecological role in regions susceptible to desertification and sandy soils. In this study, we assembled the complete S. psammophila mitochondrial genome using Pacbio HiFi third-generation sequencing data. The genome was found to be a typical single circular structure, with a total length of 715,555 bp and a GC content of 44.89%. We annotated 33 unique protein-coding genes (PCGs), which included 24 core mitochondrial genes and 9 variable genes, as well as 18 tRNA genes (5 of which were multicopy genes) and 3 rRNA genes. Comparative analysis of the PCGs from the mitochondrial genomes of S. psammophila, Populus deltoides, Populus simonii, Salix wilsonii, and Salix suchowensis revealed that these genes are relatively conserved within the Salicaceae family, with variability primarily occurring in the ribosomal protein genes. The absence of the rps14, which encodes a ribosomal protein, may have played a role in the evolution of stress tolerance in Salicaceae plants. Additionally, we identified 232 SSRs, 19 tandem repeat sequences, and 236 dispersed repeat sequences in the S. psammophila mitochondrial genome, with palindromic and forward repeats being the most abundant. The longest palindromic repeat measured 260 bp, while the longest forward repeat was 86,068 bp. Furthermore, 324 potential RNA editing sites were discovered, all involving C-to-U edits, with the nad4 having the highest number of edits. These findings provide valuable insights into the phylogenetic and genetic research of Salicaceae plants.

The complete chloroplast genome of Buxus sinica var. parvifolia (Buxaceae) and its phylogenetic analysis

Buxus sinica var. parvifolia is a shrub or small arbor of the Buxaceae family, rich in various medicinal alkaloids and of great horticultural value. In this study, we sequenced, assembled, and annotated the complete chloroplast genome of B. sinica var. parvifolia for the first time. The length of the chloroplast genome is 158,995 bp with 38.1% overall GC content. It includes a large single-copy (LSC) region of 88,140 bp, a small single-copy (SSC) region of 17,761 bp, and two inverted repeat regions of 26,547 bp. Additionally, 132 functional genes in the genome are identified, including 87 protein-coding genes, eight ribosomal RNA genes, and 37 transfer RNA genes. Phylogenetic analysis showed that B. sinica var. parvifolia is closely related to Buxus microphylla. The complete chloroplast genome sequence of B. sinica var. parvifolia and its phylogenetic analysis provides useful genomic information for the further study of B. sinica var. parvifolia and other Buxus species.

The complete chloroplast genome of Clematis serratifolia (Ranunculaceae) from Jilin province, China

Clematis serratifolia has high medicinal and ornamental value. In this study, we characterize and report, for the first time, the complete chloroplast genome sequence of C. serratifolia based on high-throughput sequence dates. The whole chloroplast genome of C. serratifolia is a circular molecule of 159,648 bp in length, consisting of a large single-copy (LSC) region of 79,394 bp, a small single-copy (SSC) region of 18,112 bp, and two inverted repeat (IR) regions of 31,071 bp. The overall GC content of the chloroplast genome is 38%, while that in the LSC, SSC, and IR regions is 36.3%, 31.3%, and 42.1%, respectively. The chloroplast genome of C. serratifolia contains 133 genes, including 89 coding genes, 8 ribosomal RNAs, and 36 transfer RNAs. Among them, 14 protein-coding genes have a single intron, and 2 genes have two introns. The phylogenetic analysis showed a close relationship between C. serratifolia and C. heracleifolia.

Characterization and phylogenetic analysis of the chloroplast genome of Solanum pseudocapsicum (Solanaceae)

Solanum pseudocapsicum Linnaeus 1753, a popular indoor potted plant known for its ornamental fruits, had its chloroplast genome sequenced in this study to determine its phylogenetic relationship with other related species and to construct a phylogenetic analysis tree. The research findings are as follows: 1. The chloroplast genome of S. pseudocapsicum comprises a large single-copy (LSC) region of 86,260 base pairs, a small single-copy (SSC) region of 18,325 base pairs, and two inverted repeat (IR) regions, each measuring 25,390 base pairs in length. 2. The G + C content of the entire chloroplast genome is 37.59%, with the highest G + C content found in the IR regions, reaching 43.03%; followed by the LSC region, which has a G + C content of 35.68%; and the lowest in the SSC region, with a G + C content of 31.53%. 3. The genome contains 127 genes, including 82 protein-coding genes, 37 tRNA genes, and 8 rRNA genes, with 18 genes duplicated in the IR regions. 4. Phylogenetic analysis revealed that S. pseudocapsicum, Solanum betaceum, Solanum laciniatum, and Solanum nitidum are genetically closely related and are located on the same branch of the phylogenetic tree, indicating a close relationship among them. This study provides a foundation for the identification, classification, and exploration of genetic diversity within the Solanum genus.

Characterization and phylogenetic analysis of the chloroplast genome in Elaeocarpus duclouxii

Elaeocarpus duclouxii, an evergreen tree species, is renowned for its fruits rich in flavonoids exhibiting potent antioxidant properties. Despite its significance, the chloroplast genome of this plant has remained unexplored until now. Our study presents the first comprehensive sequencing and analysis of the E. duclouxii chloroplast genome, revealing a circular DNA molecule of 158,148 base pairs. This genome comprises a large single-copy region of 85,700 base pairs, a small single-copy region of 17,672 base pairs, and a pair of inverted repeat regions totaling 27,388 base pairs. The genome encodes 132 genes, including 87 protein-coding genes, 37 transfer RNA genes, and 8 ribosomal RNA genes. Phylogenetic analyses indicate a close evolutionary relationship between E. duclouxii and E. sylvestris. This study not only represents the first phylogenetic investigation of E. duclouxii but also establishes a crucial genomic foundation for future research area such as conservation genetics, evolutionary biology, and potential biotechnological applications.

The complete chloroplast genome sequence of Begonia pedatifida

Begonia pedatifida has persistently been utilized as a traditional folk herbal medicine. This study has sequenced the chloroplast genome of B. pedatifida to establish its genomic characteristics and to discern its phylogenetic relationships with other closely related species. The chloroplast genome structure of B. pedatifida reveals a circular molecule with a length of 169,606 bp, including a large single copy (LSC) region of 76,086 bp, a small single copy (SSC) region of 18,314 bp, and a pair of inverted repeats (IRS) region of 37,603 bp. The entire genome contains 138 genes, which consist of 88 protein-coding genes, 42 tRNA genes, and 8 rRNA genes. Phylogenetic analysis suggests that B. pedatifida is closely related to Begonia emeiensis, Begonia jinyunensis, and Begonia pulchrifolia, sharing a common ancestor and forming sister lineages. This research provides genetic information for further study on B. pedatifida.

Complete mitochondrial genome assembly of Juglans regia unveiled its molecular characteristics, genome evolution, and phylogenetic implications

BACKGROUND

The Persian walnut (Juglans regia), an economically vital species within the Juglandaceae family, has seen its mitochondrial genome sequenced and assembled in the current study using advanced Illumina and Nanopore sequencing technology.

RESULTS

The 1,007,576 bp mitogenome of J. regia consisted of three circular chromosomes with a 44.52% GC content encoding 39 PCGs, 47 tRNA, and five rRNA genes. Extensive repetitive sequences, including 320 SSRs, 512 interspersed, and 83 tandem repeats, were identified, contributing to genomic complexity. The protein-coding sequences (PCGs) favored A/T-ending codons, and the codon usage bias was primarily shaped by selective pressure. Intracellular gene transfer occurred among the mitogenome, chloroplast, and nuclear genomes. Comparative genomic analysis unveiled abundant structure and sequence variation among J. regia and related species. The results of selective pressure analysis indicated that most PCGs underwent purifying selection, whereas the atp4 and ccmB genes had experienced positive selection between many species pairs. In addition, the phylogenetic examination, grounded in mitochondrial genome data, precisely delineated the evolutionary and taxonomic relationships of J. regia and its relatives. We identified a total of 539 RNA editing sites, among which 288 were corroborated by transcriptome sequencing data. Furthermore, expression profiling under temperature stress highlighted the complex regulation pattern of 28 differently expressed PCGs, wherein NADH dehydrogenase and ATP synthase genes might be critical in the mitochondria response to cold stress.

CONCLUSIONS

Our results provided valuable molecular resources for understanding the genetic characteristics of J. regia and offered novel perspectives for population genetics and evolutionary studies in Juglans and related woody species.

Comparative Plastomics of Plantains (Plantago, Plantaginaceae) as a Tool for the Development of Species-Specific DNA Barcodes

Plantago (plantains, Plantaginaceae) is a cosmopolitan genus including over 250 species used as functional foods, forage, and traditional medicine. Among them, Plantago lanceolata is commonly used as an ingredient of herbal products, but the close similarity to other Plantago species can cause misidentifications with potentially serious consequences for product safety/quality. To test the possibility of developing species-specific barcoding markers, we de novo assembled plastome sequences of individuals of Plantago argentea, Plantago atrata, P. lanceolata, and Plantago maritima. These genomes were characterized in comparison with both previously sequenced conspecific accessions and other publicly available plastomes, thus providing an assessment of both intraspecific and interspecific genetic variation in Plantago plastomes. Additionally, molecular evolutionary analyses indicated that eleven protein-coding genes involved in different plastid functions in Plantago plastomes underwent positive selection, suggesting they might have contributed to enhancing species' adaptation during the evolutionary history of Plantago. While the most variable mutational hotspots in Plantago plastomes were not suitable for the development of species-specific molecular markers, species-specific polymorphisms could discriminate P. lanceolata from its closest relatives. Taken together, these results highlight the potential of plastome sequencing for the development of molecular markers to improve the identification of species with relevance in herbal products.

The complete plastome of Rorippa palustris Besser 1821 and its phylogenetic analysis

Rorippa palustris Besser 1821, a species of Brassicaceae, is widely distributed around the world and used for both food and traditional Chinese medicinal purposes. Despite the plant's significance, its genetic diversity must be better understood. In this study, we have successfully assembled and characterized a complete plastome of R. palustris, marking a significant advancement toward comprehending its genetic composition. The plastome is 154,674 bp long and harbors 128 genes, including 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Our phylogenomic analysis indicated that R. palustris is closely related to R. curvipes. These findings are crucial for conserving and utilizing this important plant species. They also highlight the potential for future research into the evolution and preservation of R. palustris, which could be advantageous in pharmaceutical applications.

Characteristics and phylogenetic analysis of the complete chloroplast genome of Mesembryanthemum cordifolium L. F. (Aizoaceae)

Mesembryanthemum cordifolium, a perennial plant with crassulacean acid metabolism (CAM) in the Aizoaceae family, has significant ornamental and medicinal values. In this study, we reported the first complete chloroplast genome sequence of this species. The total genome size was 153,734 bp in length, including a large single-copy (LSC) region of 85,692 bp, a small single-copy (SSC) region of 18,212 bp, and a pair of inverted repeat (IR) regions of 24,915 bp by each. The overall GC content of the M. cordifolium chloroplast genome was 37.08%. The genome encodes 131 genes, comprising 87 protein-coding genes (PCGs), 36 transfer RNA genes (tRNAs), and eight ribosomal RNA genes (rRNAs). Phylogenetic analysis shows this species was relatively close to M. crystallinum. This chloroplast genome sequence will be valuable for species discrimination and for understanding phylogenetic relationships within the genus Mesembryanthemum.

Analysis of codon usage patterns in complete plastomes of four medicinal Polygonatum species (Asparagaceae)

Polygonati Rhizoma and Polygonati odorati Rhizoma, known as "Huangjing" and "Yuzhu" in China, are medicinal Polygonatum species resources with top-grade medical and edible properties. The chloroplast (cp) genome has been used to study species diversity, evolution, and breeding of species for applications in genetic engineering. Codon usage bias (CUB), a common and complex natural phenomenon, is essential for studies of codon optimization of exogenous genes, genetic engineering, and molecular evolution. However, the CUB of medicinal Polygonatum species chloroplast genomes has not been systematically studied. In our study, a detailed analysis of CUB was performed in the medicinal Polygonatum species chloroplast genomes. We investigated the codon bias of 204 plastid protein-coding genes (PCGs) in 4 medicinal Polygonatum species using CodonW and CUSP online software. Through the analysis of the codon bias index, we found that the medicinal Polygonatum species chloroplast genomes had weak codon usage bias. In addition, our results also showed a high preference for AT bases in medicinal Polygonatum species chloroplast genomes, and the preference to use AT-ending codons was observed in these species chloroplast genomes. The neutrality plot, ENC plot, PR2-Bias plot, and correspondence analysis showed that compared with mutation pressure, natural selection was the most important factor of CUB. Based on the comparative analysis of high-frequency codons and high expression codons, we also determined the 10-11 optimal codons of investigative medicinal Polygonatum species. Furthermore, the result of RSCU-based cluster analysis showed that the genetic relationship between different medicinal Polygonatum species could be well reflected. This study provided an essential understanding of CUB and evolution in the medicinal Polygonatum species chloroplast genomes.

Insights into Aquilaria phylogenetics through comparative plastomic resources

The plastid is an essential organelle for its role in photosynthesis and energy production and its genomic information is always employed as important evolutionary markers to explore the relationship among species. Agarwood (Aquilaria), prized for its aromatic blend, finds extensive use in various cultures as incense and perfume. Despite its high economic importance, the phylogenetic status among Aquilaria based on plastomes remains ambiguous due to the lack of available plastomic resources. To bridge this knowledge gap, 22 Aquilaria plastomes were newly sequenced, similar variation patterns in this genus were determined, including a shared 16 bp extension of the rps19 gene and seven highly variable regions. The analysis highlighted the highest prevalence of the A/T motif among simple sequence repeats in these plastomes. Further phylogenetic analysis revealed Aquilaria's phylogenetic implications with an expanded dataset. This comprehensive plastomic resource not only enhances our understanding of Aquilaria evolution but also presents potential molecular markers for DNA barcoding.

De novo assembly and comprehensive analysis of the mitochondrial genome of Taxus wallichiana reveals different repeats mediate recombination to generate multiple conformations

Taxus plants are the exclusive source of paclitaxel, an anticancer drug with significant medicinal and economic value. Interspecies hybridization and gene introgression during evolution have obscured distinctions among Taxus species, complicating their phylogenetic classification. While the chloroplast genome of Taxus wallichiana, a widely distributed species in China, has been sequenced, its mitochondrial genome (mitogenome) remains uncharacterized.We sequenced and assembled the T. wallichiana mitogenome using BGI short reads and Nanopore long reads, facilitating comparisons with other gymnosperm mitogenomes. The T. wallichiana mitogenome spanning 469,949 bp, predominantly forms a circular configuration with a GC content of 50.51%, supplemented by 3 minor configurations mediated by one pair of LRs and two pairs of IntRs. It includes 32 protein-coding genes, 7 tRNA genes, and 3 rRNA genes, several of which exist in multiple copies.We detailed the mitogenome's structure, codon usage, RNA editing, and sequence migration between organelles, constructing a phylogenetic tree to elucidate evolutionary relationships. Unlike typical gymnosperm mitochondria, T. wallichiana shows no evidence of mitochondrial-plastid DNA transfer (MTPT), highlighting its unique genomic architecture. Synteny analysis indicated extensive genomic rearrangements in T. wallichiana, likely driven by recombination among abundant repetitive sequences. This study offers a high-quality T. wallichiana mitogenome, enhancing our understanding of gymnosperm mitochondrial evolution and supporting further cultivation and utilization of Taxus species.

Structural analysis of the mitochondrial genome of Santalum album reveals a complex branched configuration

BACKGROUND

Santalum album L. is an evergreen tree which is mainly distributes throughout tropical and temperate regions. And it has a great medicinal and economic value.

RESULTS

In this study, the complete mitochondrial genome of S. album were assembled and annotated, which could be descried by a complex branched structure consisting of three contigs. The lengths of these three contigs are 165,122 bp, 93,430 bp and 92,491 bp. We annotated 34 genes coding for proteins (PCGs), 26 tRNA genes, and 4 rRNA genes. The analysis of repeated elements shows that there are 89 SSRs and 242 pairs of dispersed repeats in S. album mitochondrial genome. Also we found 20 MTPTs among the chloroplast and mitochondria. The 20 MTPTs sequences span a combined length of 22,353 bp, making up 15.52 % of the plastome, 6.37 % of the mitochondrial genome. Additionally, by using the Deepred-mt tool, we found 628 RNA editing sites in 34 PCGs. Moreover, significant genomic rearrangement is observed between S. album and its associated mitochondrial genomes. Finally, based on mitochondrial genome PCGs, we deduced the phylogenetic ties between S. album and other angiosperms.

CONCLUSIONS

We reported the mitochondrial genome from Santalales for the first time, which provides a crucial genetic resource for our study of the evolution of mitochondrial genome.

Assembly and comparative analysis of the complete mitochondrial genome of Vaccinium carlesii Dunn

Vaccinium L. is an important fruit tree with nutritional, medicinal, and ornamental values. However, the mitochondrial (mt) genome of Vaccinium L. remains largely unexplored. Vaccinium carlesii Dunn is an endemic wild resource in China, which is crucial for blueberry breeding. The V. carlesii mt genomes were sequenced using Illumina and Nanopore, which total length was 636,904 bp with 37 protein coding genes, 20 tRNA genes, and three rRNA genes. We found four pairs of long repeat fragments homologous recombination mediated the generation of substructures in the V. carlesii mt genome. We predicted 383 RNA editing sites, all converting cytosine (C) to uracil (U). According to the phylogenetic analysis, V. carlesii and V. macrocarpon of the Ericaceae exhibited the closest genetic relationship. This study provides a theoretical basis for understanding the evolution of higher plants, species classification and identification, and will also be useful for further utilization of Vaccinium germplasm resources.

The complete plastome of Saxifraga giraldiana Engler and its phylogenetic analysis

Saxifraga giraldiana Engler is a common subalpine and alpine plant belonging to Saxifragaceae. However, the genetic diversity of this species has remained to be explored. In this study, we have assembled and characterized the complete chloroplast genome of S. giraldiana, filling this knowledge gap and uncovering its genetic composition. The chloroplast genome is 147,267 bp long and contains 131 genes, including 85 protein-coding genes, 38 tRNA genes, and eight rRNA genes. Furthermore, we have performed a phylogenetic analysis of 19 representative species within Saxifraga. As a result, we have found that S. giraldiana, together with S. implicans and S. stellariifolia, forms a monophyletic group. These findings have implications for the conservation and utilization of S. giraldiana.