PGA: a software package for rapid, accurate, and flexible batch annotation of plastomes

Comparative plastome analysis of Arundinelleae (Poaceae, Panicoideae), with implications for phylogenetic relationships and plastome evolution

BACKGROUND

Arundinelleae is a small tribe within the Poaceae (grass family) possessing a widespread distribution that includes Asia, the Americas, and Africa. Several species of Arundinelleae are used as natural forage, feed, and raw materials for paper. The tribe is taxonomically cumbersome due to a paucity of clear diagnostic morphological characters. There has been scant genetic and genomic research conducted for this group, and as a result the phylogenetic relationships and species boundaries within Arundinelleae are poorly understood.

RESULTS

We compared and analyzed 11 plastomes of Arundinelleae, of which seven plastomes were newly sequenced. The plastomes range from 139,629 base pairs (bp) (Garnotia tenella) to 140,943 bp (Arundinella barbinodis), with a standard four-part structure. The average GC content was 38.39%, but varied in different regions of the plastome. In all, 110 genes were annotated, comprising 76 protein-coding genes, 30 tRNA genes, and four rRNA genes. Furthermore, 539 simple sequence repeats, 519 long repeats, and 10 hyper-variable regions were identified from the 11 plastomes of Arundinelleae. A phylogenetic reconstruction of Panicoideae based on 98 plastomes demonstrated the monophyly of Arundinella and Garnotia, but the circumscription of Arundinelleae remains unresolved.

CONCLUSION

Complete chloroplast genome sequences can improve phylogenetic resolution relative to single marker approaches, particularly within taxonomically challenging groups. All phylogenetic analyses strongly support the monophyly of Arundinella and Garnotia, respectively, but the monophylly of Arundinelleae was not well supported. The intergeneric phylogenetic relationships within Arundinelleae require clarification, indicating that more data is necessary to resolve generic boundaries and evaluate the monophyly of Arundinelleae. A comprehensive taxonomic revision for the tribe is necessary. In addition, the identified hyper-variable regions could function as molecular markers for clarifying phylogenetic relationships and potentially as barcoding markers for species identification in the future.

Comparative Phylogenomic Study of Malaxidinae (Orchidaceae) Sheds Light on Plastome Evolution and Gene Divergence

Malaxidinae is one of the most confusing groups in the Orchidaceae classification. Previous phylogenetic analyses have revealed that the relationships between the taxa in Malaxidinae have not yet been reliably established, using only a few plastome regions and nuclear ribosomal internal transcribed spacer (nrITS). In the present study, the complete plastomes of Oberonia integerrima and Crepidium purpureum were assembled using high-throughput sequencing. Combined with publicly available complete plastome data, this resulted in a dataset of 19 plastomes, including 17 species of Malaxidinae. The plastome features and phylogenetic relationships were compared and analyzed. The results showed the following: (1) Malaxidinae species plastomes possess the quadripartite structure of typical angiosperms, with sizes ranging from 142,996 to 158,787 bp and encoding from 125 to 133 genes. The ndh genes were lost or pseudogenized to varying degrees in six species. An unusual inversion was detected in the large single-copy region (LSC) of Oberonioides microtatantha. (2) Eight regions, including ycf1, matK, rps16, rpl32, ccsA-ndhD, clpP-psbB, trnFGAA-ndhJ, and trnSGCU-trnGUCC, were identified as mutational hotspots. (3) Based on complete plastomes, 68 protein-coding genes, and 51 intergenic regions, respectively, our phylogenetic analyses revealed the genus-level relationships in this subtribe with strong support. The Liparis was supported as non-monophyletic.

Assessing Genetic Diversity in Endangered Plant Orchidantha chinensis: Chloroplast Genome Assembly and Simple Sequence Repeat Marker-Based Evaluation

Orchidantha chinensis T. L. Wu, an endemic species in China, is listed as a key protected wild plant in Guangdong Province. However, the lack of reports on the chloroplast genome and simple sequence repeat (SSR) markers has hindered the assessment of its genetic diversity and conservation strategies. The limited number of molecular markers to assess the genetic diversity of this species, and thus develop proper conservation strategies, highlighted the urgent need to develop new ones. This study developed new SSR markers and investigated genetic variation using 96 samples of O. chinensis from seven populations. Through high-throughput sequencing, a complete chloroplast genome of 134,407 bp was assembled. A maximum-likelihood phylogenetic tree, based on the chloroplast genome, showed that O. chinensis is closely related to Ravenala madagascariensis. The study identified 52 chloroplast SSRs (cpSSRs) and 5094 expressed sequence tag SSRs (EST-SSRs) loci from the chloroplast genome and leaf transcriptome, respectively. Twenty-one polymorphic SSRs (seven cpSSRs and fourteen EST-SSRs) were selected to evaluate the genetic variation in 96 accessions across seven populations. Among these markers, one cpSSR and 11 EST-SSRs had high polymorphism information content (>0.5). Cluster, principal coordinate, and genetic structure analyses indicated that groups G1 and G6 were distinct from the other five groups. However, an analysis of molecular variance showed greater variation within groups than among groups. The genetic distance among the populations was significantly positively correlated with geographical distance. These findings provide new markers for studying the genetic variability of O. chinensis and offer a theoretical foundation for its conservation strategies.

The complete chloroplast genome sequence of Primula medogensis (Primulaceae) and its phylogeny

Primula medogensis W.B Ju, B. Xu & X.F. Gao 2023, a new species categorized under P. sect. Cordifoliae, was officially described in 2023. Given its recent classification, the genetic resources for this species are currently very limited. Here, we sequenced and assembled the first complete chloroplast genome of P. medogensis using Illumina sequencing technology. The complete chloroplast genome of P. medogensis is 151,486 bp in length, exhibiting a typical quadripartite structure. It consists of a large single-copy region (LSC; 83,407 bp) and a small single-copy region (SSC;17675 bp), separated by a pair of inverted repeat regions (IRs; 25202 bp). A total of 131 genes were annotated, including 86 protein-coding, 37 tRNA, and eight rRNA genes. The overall GC content was 37.1%. Phylogenetic analysis of 59 Primula species revealed a close relationship between P. medogensis and P. calliantha subsp. bryophila.

Characterization of the complete chloroplast genome of Liparis gigantea (Orchidaceae)

Liparis gigantea is a Chinese traditional medicinal herb in the Orchidaceae family. It is a rare and special Liparis species that exhibits relatively large flowers. To illuminate its phylogenetic status and augment genomic resources, the complete chloroplast (cp) genome of L. gigantea was first sequenced and assembled using whole genome next-generation sequencing in this study. The cp genome size is 158,462 bp with a total GC content of 36.9%. Characterized by a quadripartite structure, the genome consists of a large single-copy (LSC) region of 86,032 bp, a small single-copy (SSC) region of 18,322 bp, which is separated by a pair of 27,054 bp inverted repeat regions (IRs). A total of 133 genes were annotated, including 87 protein-coding genes, 38 tRNA genes, and 8 rRNA genes. Phylogenetic analysis strongly supported L. gigantea as the sister to two closely related terrestrial species, Liparis nervosa and L. vivipara. The results of this study provide genomic information for future research and application of this medicinal herb.

From light into shadow: comparative plastomes in Petrocosmea and implications for low light adaptation

BACKGROUND

Plastids originated from an ancient endosymbiotic event and evolved into the photosynthetic organelles in plant cells. They absorb light energy and carbon dioxide, converting them into chemical energy and oxygen, which are crucial for plant development and adaptation. However, little is known about the plastid genome to light adaptation. Petrocosmea, a member of the Gesneriaceae family, comprises approximately 70 species with diverse light environment, serve as an ideal subject for studying plastomes adapt to light.

RESULTS

In this study, we selected ten representative species of Petrocosmea from diverse light environments, assembled their plastid genomes, and conducted a comparative genomic analysis. We found that the plastid genome of Petrocosmea is highly conserved in both structure and gene content. The phylogenetic relationships reconstructed based on the plastid genes were divided into five clades, which is consistent with the results of previous studies. The vast majority of plastid protein-coding genes were under purifying selection, with only the rps8 and rps16 genes identified under positive selection in different light environments. Notably, significant differences of evolutionary rate were observed in NADH dehydrogenase, ATPase ribosome, and RNA polymerase between Clade A and the other clades. Additionally, we identified ycf1 and several intergenic regions (trnH-psbA, trnK-rps16, rpoB-trnC, petA-psbJ, ccsA-trnL, rps16-trnQ, and trnS-trnG) as candidate barcodes for this emerging ornamental horticulture.

CONCLUSION

We newly assembled ten plastid genomes of Petrocosmea and identified several hypervariable regions, providing genetic resources and candidate markers for this promising emerging ornamental horticulture. Furthermore, our study suggested that rps8 and rps16 were under positive selection and that the evolutionary patterns of NADH dehydrogenase, ATPase ribosome, and RNA polymerase were related to the diversity light environment in Petrocosmea. This revealed an evolutionary scenario for light adaptation of the plastid genome in plants.

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.

The complete chloroplast genome of Camellia huulungensis Rosmann et Ninh, a golden Camellia species endemic to Vietnam

Camellia huulungensis Rosmann & Ninh 1997, belonging to the sect. Chrysantha, holds important ornamental value and medicinal value. In this study, the complete chloroplast genome sequence of C. huulungensis was assembled using high-throughput sequencing technology. The entire length of chloroplast genome is 156,546 bp and contains a small single-copy region (18,257 bp), a large single-copy region (86,219 bp), and a pair of inverted repeat regions (26,035 bp). A total of 133 genes were annotated, including 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. The overall GC content is 37.33%. The phylogenetic analysis showed that C. huulungensis is sister to C. aurea. The results can provide genetic data for further phylogenetic studies of Camellia.

Resurrection of Sorbustapashana (Rosaceae) based on molecular and morphological evidence

Sorbustapashana (Rosaceae) from Taibai Shan, Qinling, China, has been treated as a synonym of S.tianschanica. Both species belong to a distinctive group characterized by white tomentose buds, relatively large flowers, and red fruits. However, these two species do not cluster together in the plastome-based phylogenetic analysis. Morphologically, S.tapashana differs from S.tianschanica by its persistent white tomentose on the peduncle, pedicels, rachis, both sides of the midrib on abaxial surface, its leaflets with 31-51 teeth on each side and much smaller corymbs and fruits. Therefore, S.tapashana is reinstated as a distinct species here.

The complete chloroplast genome of Tabebuia rosea (Bignoniaceae)

Tabebuia rosea is a world-renowned woody plant with colorful flowers in full bloom. In addition to its high ornamental value, it also has ecological and medicinal value. In this study, the complete circular chloroplast genome of T. rosea was reconstructed and annotated using Illumina sequencing. The chloroplast genome was 158,919 bp in size with GC content of 38.21%, including a large single-copy region of 85,823 bp, a small single-copy region of 12,816 bp, and a pair of inverted repeats of 30,140 bp. It encoded 132 genes, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Based on current available chloroplast genome sequences, the phylogenetic analysis indicated that T. rosea was clustered with T. nodosa and H. chrysanthus. This study provided insights into the evolutionary relationships among different species of Bignoniaceae.

The complete chloroplast genome sequence of Lilium saccatum S. Yun Liang (Liliaceae, Lilieae) and its phylogenetic analysis

Lilium saccatum is a species of ornamental plant found in southeastern Xizang, China. In the present study, the complete chloroplast (cp) genome of L. saccatum was sequenced using next-generation sequencing (NGS). The de novo assembled cp genome was 151,839 bp in length, including a pair of inverted repeat regions (IRs; 26,421 bp), a small single-copy region (SSC; 17,528 bp), and a large single-copy region (LSC; 81,469 bp). The cp genome encodes 113 unique genes, including 79 protein-coding genes (PCGs), 30 tRNA genes, and four rRNA genes. The total GC content of the cp genome was 37.0%. Phylogenetic analysis of 24 cp genomes revealed that L. saccatum was closely related to L. souliei. This study could provide fundamental information for the phylogenomics and utilization of Lilium.

Deciphering the evolution and biogeography of Ant-ferns Lecanopteris s.s

Southeast Asia is a biodiversity hotspot characterized by a complex paleogeography, and its Polypodiopsida flora is particularly diverse. While hybridization is recognized as common in ferns, further research is needed to investigate the relationship between hybridization events and fern diversity. Lecanopteris s.s., an ant-associated fern, has been subject to debate regarding species delimitations primarily due to limited DNA markers and species sampling. Our study integrates 22 newly generated plastomes, 22 transcriptomes, and flow cytometry of all native species along with two cultivated hybrids. Our objective is to elucidate the reticulate evolutionary history within Lecanopteris s.s. through the integration of phylobiogeographic reconstruction, gene flow inference, and genome size estimation. Key findings of our study include: (1) An enlarged plastome size (178-187 Kb) in Lecanopteris, attributed to extreme expansion of the Inverted Repeat (IR) regions; (2) The traditional 'pumila' and 'crustacea' groups are paraphyletic; (3) Significant cytonuclear discordance attributed to gene flow; (4) Natural hybridization and introgression in the 'pumila' and 'darnaedii' groups; (5) L. luzonensis is the maternal parent of L. 'Yellow Tip', with L. pumila suggested as a possible paternal parent; (6) L. 'Tatsuta' is a hybrid between L. luzonensis and L. crustacea; (7) Lecanopteris first diverged during the Neogene and then during the middle Miocene climatic optimum in the Indochina and Sundaic regions. In conclusion, the biogeographic history and speciation of Lecanopteris have been profoundly shaped by past climate changes and geodynamics of Southeast Asia. Dispersals, hybridization and introgression between species act as pivotal factors in the evolutionary trajectory of Lecanopteris. This research provides a robust framework for further exploration and understanding of the complex dynamics driving the diversification and distribution patterns within Polypodiaceae subfamily Microsoroideae.

Plant species diversification in the Himalaya-Hengduan Mountains region: an example from an endemic lineage of Pedicularis (Orobanchaceae) in the role of floral specializations and rapid range expansions

The c. 270 endemic species of Pedicularis in the Himalaya-Hengduan Mountains (HHM) region exhibit high diversity in geographic distribution, elevational range and floral morphology. Many of these, including the species with the longest corolla tubes and beaked galeas, are monophyletic and represent a putative in situ radiation. In this study, we focus on the representative Clade 3 within the HHM region. We integrate the plastid phylogeny of this clade with environmental data and species distributions to infer environmental correlates of species diversity. We estimate macroevolutionary rates and reconstructed ancestral states for geographic ranges and corolla traits, and analyse patterns of range overlap and niche evolution to assess drivers of diversification in the HHM region. Our results show that the region from northwest Yunnan to southwest Sichuan is the centre of diversity for this clade of Pedicularis. Rates of diversification are associated with precipitation and multiple environmental factors. Multiple range expansions from the Sanjiang (Three Parallel Rivers) region, followed by allopatric speciation across the HHM region, contributed to early rapid diversification. Corolla traits are not significantly associated with species diversification. This study highlights the importance of integrated evidence for understanding species diversification dynamics and contributes to our understanding of the origins of the remarkable richness of plant species in the HHM region.

The complete chloroplast genome of Dryas octopetala var. asiatica (Dryadoideae, Rosaceae) and phylogenetic analysis

Dryas octopetala L. var. asiatica (Nakai) Nakai 1918 is a dwarf shrub that mainly grow in alpine and arctic zones of the Northern Hemisphere, representing an endemic variety in Asia. In the present study, the complete chloroplast (cp) genome of D. octopetala var. asiatica was first characterized and used for its phylogenetic analysis. The cp genome span 158,271 bp with an overall GC content of 36.5%. A total of 129 genes were identified, including 84 protein-coding genes (PCGs), 37 tRNA genes, and 8 rRNA genes. In addition, repetitive sequences and microsatellites were detected within this species. Phylogenetic analysis involving 39 cp genomes from Rosaceae family indicated that D. octopetala var. asiatica was sister to the clade of Amygdaloideae. This study contributes fundamental insights into the cp genome of Dryas octopetala var. asiatica, which will have expanded its use in photosynthesis and evolutionary study.

Complete mitochondrial genome of Agropyron cristatum reveals gene transfer and RNA editing events

BACKGROUND

As an important forage in arid and semi-arid regions, Agropyron cristatum provides livestock with exceptionally high nutritional value. Additionally, A. cristatum exhibits outstanding genetic characteristics to endure drought and disease. Therefore, rich genetic diversity serves as a cornerstone for the improvement of major food crops. The purposes of this study were to systematically describe mitogenome of A.cristatum and preliminarily analyze its internal variations.

RESULT

The A. cristatum mitogenome was a single-ring molecular structure of 381,065 bp that comprised 52 genes, including 35 protein-coding, 3 rRNA and 14 tRNA genes. Among these, two pseudoprotein-coding genes and multiple copies of tRNA genes were observed. A total of 320 repetitive sequences was found to cover more than 10% of the mitogenome (105 simple sequences, 185 dispersed and 30 tandem repeats), which led to a large number of fragment rearrangements in the mitogenome of A. cristatum. Leucine was the most frequent amino acid (n = 1087,10.8%) in the protein-coding genes of A. cristatum mitogenome, and the highest usage codon was ATG (initiation codon). The number of A/T changes at the third base of the codon was much higher than that of G/C. Among 23 PCGs, the range of Pi values is from 0.0021 to 0.0539, with an average of 0.013. Additionally, 81 RNA editing sites were predicted, which were considerably fewer than those reported in other plant mitogenomes. Most of the RNA editing site base positions were concentrated at the first and second codon bases, which were C to T transitions. Moreover, we identified 95 sequence fragments (total length of 34, 343 bp) that were transferred from the chloroplast to mitochondria genes, introns, and intergenic regions. The stability of the tRNA genes was maintained during this process. Selection pressure analysis of 23 protein-coding genes shared by 15 Poaceae plants, showed that most genes were subjected to purifying selection during evolution, whereas rps4, cob, mttB, and ccmB underwent positive selection in different plants. Finally, a phylogenetic tree was constructed based on 22 plant mitogenomes, which showed that Agropyron plants have a high degree of independent heritability in Triticeae.

CONCLUSION

The findings of this study provide new data for a better understanding of A. cristatum genes, and demonstrate that mitogenomes are suitable for the study of plant classifications, such as those of Agropyron. Moreover, it provides a reference for further exploration of the phylogenetic relationships within Agropyron species, and establishes a theoretical basis for the subsequent development and utilization of A. cristatum plant germplasm resources.

Phylogenomics analysis of Scutellaria (Lamiaceae) of the world

BACKGROUND

Scutellaria, a sub-cosmopolitan genus, stands as one of the Lamiaceae family's largest genera, encompassing approximately 500 species found in both temperate and tropical montane regions. Recognized for its significant medicinal properties, this genus has garnered attention as a research focus, showcasing anti-cancer, anti-inflammatory, antioxidant, and hepatoprotective qualities. Additionally, it finds application in agriculture and horticulture. Comprehending Scutellaria's taxonomy is pivotal for its effective utilization and conservation. However, the current taxonomic frameworks, primarily based on morphological characteristics, are inadequate. Despite several phylogenetic studies, the species relationships and delimitations remain ambiguous, leaving the genus without a stable and reliable classification system.

RESULTS

This study analyzed 234 complete chloroplast genomes, comprising 220 new and 14 previously published sequences across 206 species, subspecies, and varieties worldwide. Phylogenetic analysis was conducted using six data matrices through Maximum Likelihood and Bayesian Inference, resulting in a robustly supported phylogenetic framework for Scutellaria. We propose three subgenera, recommending the elevation of Section Anaspis to subgeneric rank and the merging of Sections Lupulinaria and Apeltanthus. The circumscription of Subgenus Apeltanthus and Section Perilomia needs to be reconsidered. Comparative analysis of chloroplast genomes highlighted the IR/SC boundary feature as a significant taxonomic indicator. We identified a total of 758 SSRs, 558 longer repetitive sequences, and ten highly variable regions, including trnK-rps16, trnC-petN, petN-psbM, accD-psaI, petA-psbJ, rpl32-trnL, ccsA-ndhD, rps15-ycf1, ndhF, and ycf1. These findings serve as valuable references for future research on species identification, phylogeny, and population genetics.

CONCLUSIONS

The phylogeny of Scutellaria, based on the most comprehensive sample collection to date and complete chloroplast genome analysis, has significantly enhanced our understanding of its infrageneric relationships. The extensive examination of chloroplast genome characteristics establishes a solid foundation for the future development and utilization of Scutellaria, an important medicinal plant globally.

Plastome Evolution and Comparative Analyses of a Recently Radiated Genus Vanda (Aeridinae, Orchidaceae)

Vanda R.Br. is an epiphytic orchid genus with significant horticultural and ornamental value. Previous molecular studies expanded Vanda including some members from five other genera. However, the interspecific relationships of this recently radiated genus have remained unclear based on several DNA markers until now. In this study, the complete plastome has been used to infer the phylogenetic relationships of Vanda s.l. The five newly obtained plastomes ranged from 146,340 bp to 149,273 bp in length, with a GC content ranging from 36.5% to 36.7%. The five plastomes contained 74 protein-coding genes (CDSs), 38 tRNAs, and 8 rRNAs, and their ndh genes underwent loss or pseudogenization. Comparative plastome analyses of 13 Vanda species revealed high conservation in terms of genome size, structure, and gene order, except for a large inversion from trnGGCC to ycf3 in V. coerulea. Moreover, six CDSs and five non-CDSs were selected as candidate DNA barcodes. Our phylogenetic analyses demonstrated that Vanda s.l. is a monophyletic group with high supporting values based on five different datasets (complete plastome with one IR, 68 CDSs, LSC, five hypervariable non-CDSs, and six hypervariable CDSs), while the phylogenetic relationships among species were fully resolved based on the complete plastome with one IR dataset. Our results confirmed that the complete plastome has a great power in resolving the phylogenetic relationships of recently radiated lineages.

Rampant intraspecific variation of plastid genomes in Gentiana section Chondrophyllae

Exploring the level of intraspecific diversity in taxa experienced radiation is helpful to understanding speciation and biodiversity assembly. Gentiana section Chondrophyllae sensu lato encompasses more than 180 species and occupies more a half of species in the genus. In this study, we collected samples across the range of three species (Gentiana aristata, G. crassuloides and G. haynaldii) in section Chondrophyllae s.l., and recovered the intra-species variation by comparing with closely related taxon. Using 25 newly sequenced plastid genomes together with previously published data, we compared structural differences, quantified the variations in plastome size, and measured nucleotide diversity in various regions. Our results showed that the plastome size variation in the three Chondrophyllae species ranged from 285 to 628 bp, and the size variation in LSC, IR and SSC ranged from 236 to 898 bp, 52 to 393 bp and 135 to 356 bp, respectively. Nucleotide diversity of plastome or any of the four regions was much higher than the control species. The average nucleotide diversity in plastomes of the three species ranged from 0.0010 to 0.0023 in protein coding genes, and from 0.0023 to 0.0061 in intergenic regions. More repeat sequence variations were detected within the three Chondrophyllae species than the control species. Various plastid sequence matrixes resulted in different backbone topology in two target species, showed uncertainty in phylogenetic relationship based inference. In conclusion, our results recovered that species of G. section Chondrophyllae s.l. has high intraspecific plastome variation, and provided insights into the radiation in this speciose lineage.

Phylogenomic analyses revealed widely occurring hybridization events across Elsholtzieae (Lamiaceae)

Obtaining a robust phylogeny proves challenging due to the intricate evolutionary history of species, where processes such as hybridization and incomplete lineage sorting can introduce conflicting signals, thereby complicating phylogenetic inference. In this study, we conducted comprehensive sampling of Elsholtzieae, with a particular focus on its largest genus, Elsholtzia. We utilized 503 nuclear loci and complete plastome sequences obtained from 99 whole-genome sequencing datasets to elucidate the interspecific relationships within the Elsholtzieae. Additionally, we explored various sources of conflicts between gene trees and species trees. Fully supported backbone phylogenies were recovered, and the monophyly of Elsholtzia and Keiskea was not supported. Significant gene tree heterogeneity was observed at numerous nodes, particularly regarding the placement of Vuhuangia and the E. densa clade. Further investigations into potential causes of this discordance revealed that incomplete lineage sorting (ILS), coupled with hybridization events, has given rise to substantial gene tree discordance. Several species, represented by multiple samples, exhibited a closer association with geographical distribution rather than following a strictly monophyletic pattern in plastid trees, suggesting chloroplast capture within Elsholtzieae and providing evidence of hybridization. In conclusion, this study provides phylogenomic insights to untangle taxonomic problems in the tribe Elsholtzieae, especially the genus Elsholtzia.

Complete chloroplast genome sequence and phylogenetic analysis of Populus maximowiczii Henry 1913 (Salicaceae Mirb.)

The classification and identification of species in Populus has remained a formidable challenge due to widespread interspecies hybridization. The complete chloroplast genome of Populus maximowiczii was obtained by Illumina high-throughput sequencing technology, with a typical quadripartite structure and 37.0% GC content. The chloroplast genome of P. maximowiczii was 156,892 in length, including a large single-copy region (LSC: 84,988 bp), a small single-copy region (SSC: 16,630 bp), and a pair of inverted repeats (IRs: each 27,637 bp in length). A total of 131 genes were annotated, including 86 protein-coding genes, 37 tRNAs, and 8 rRNAs. The phylogenetic analysis indicated that 43 species belonging to Populus were classified into monophyly, with P. cathayana being the closest relatives to P. maximowiczii. In conclusion, this study provides valuable insights into understanding the phylogeny of Populus.

The complete chloroplast genome of Camellia flava (Pitard) Sealy, a golden camellia of Vietnam

Camellia flava (Pit.) Sealy 1949 is a rare and precious species with golden flowers, which hold important ornamental and breeding values. In this study, the complete chloroplast genome of C. flava is reported for the first time. The chloroplast genome exhibits a typical quadripartite structure with a total length of 156,670 bp and a GC content of 37.32%, including a large single-copy region (86,250 bp), a small single-copy region (18,292 bp), and a pair of inverted repeat regions (26,064 bp). A total of 133 genes, including 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes were annotated. The phylogenetic analysis revealed a close relationship between C. flava and C. tamdaoensis. The chloroplast genome sequence of C. flava serves as a valuable resource for further breeding research and genetic phylogenetic studies.

The complete chloroplast genome of Semiaquilegia danxiashanensis (Ranunculaceae), a rare species endemic to Danxia landform in Guangdong, China

Semiaquilegia danxiashanensis is currently known only from the type locality, Danxia Mountain, characterized by its spectacular red sandstone cliffscape. In this study, we assembled the complete chloroplast genome sequence of S. danxiashanensis and inferred its phylogenetic relationships. Total length of the chloroplast genome was 160,548 bp, with an overall GC content of 39%. The chloroplast genome had typical quadripartite structure and contained one LSC region (89,882 bp) and one SSC region (17,386 bp), which were separated by two IRs regions (26,640 bp, respectively). It comprised 133 genes, including 84 protein coding genes, 41 tRNA genes and eight rRNA genes. The maximum likelihood phylogenetic analysis indicated that S. danxiashanensis was sister to S. adoxoides; meanwhile, Semiaquilegia was closely related to both Urophysa and Aquilegia in Ranunculaceae. This study sheds light on the evolutionary history of Semiaquilegia and provides preliminary data for future comparative analysis of chloroplast genomes.

Analysis of chloroplast genome structure and phylogeny of the traditional medicinal of Ardisia crispa (Myrsinaceae)

Ardisia crispa(Myrsinaceae) is an ethnomedicine with horticultural and important medicinal values. Its morphology is complex, and its identification is difficult. We analyse the chloroplast genome characteristics and phylogenetic position of A. crispa to provide basic research data for the identification of A. crispa species and resource conservation. This study assemble and annotate the chloroplast genome of A. crispa and to compare it with the chloroplast genome within Ardisia. The A. crispa chloroplast genome is 156,785 bp in length, with a typical quadripartite structure containing 131 genes, including 86 protein-coding genes, 37 tRNA genes, and 8 rRNA genes; a total of 59 SSRs sites were identified, and the codon preference of this chloroplast genome is greater in A/U than in G/C, and leucine is the amino acid with the highest frequency of use. The chloroplast genomes of the nine Ardisia species are conserved in gene content and number, with more stable boundaries and less variation. In the phylogenetic tree, A. crispa is clustered on a branch with A. crispa var dielsii, and is closely related to A. mamillata and A. pedalis. In this study, we constructed and analyzed the chloroplast genome structure of A. crispa, and conducted phylogenetic analysis using the whole chloroplast genome sequence data of Ardisia plants, which is of great significance in understanding the genetic basis of A. crispa and adaptive evolution in Ardisia plants, and this will lay the foundation for the future research on A. crispa resource conservation and species identification.

Phylogeny, adaptive evolution, and taxonomy of Acronema (Apiaceae): evidence from plastid phylogenomics and morphological data

Introduction

The genus Acronema, belonging to Apiaceae, includes approximately 25 species distributed in the high-altitude Sino-Himalayan region from E Nepal to SW China. This genus is a taxonomically complex genus with often indistinct species boundaries and problematic generic delimitation with Sinocarum and other close genera, largely due to the varied morphological characteristics.

Methods

To explore the phylogenetic relationships and clarify the limits of the genus Acronema and its related genera, we reconstructed a reliable phylogenetic framework with high support and resolution based on two molecular datasets (plastome data and ITS sequences) and performed morphological analyses.

Results

Both phylogenetic analyses robustly supported that Acronema was a non-monophyletic group that fell into two clades: Acronema Clade and East-Asia Clade. We also newly sequenced and assembled sixteen Acronema complete plastomes and performed comprehensively comparative analyses for this genus. The comparative results showed that the plastome structure, gene number, GC content, codon bias patterns were high similarity, but varied in borders of SC/IR and we identified six different types of SC/IR border. The SC/IR boundaries of Acronema chienii were significantly different from the other Acronema members which was consistent with the type VI pattern in the genus Tongoloa. We also identified twelve potential DNA barcode regions (ccsA, matK, ndhF, ndhG, psaI, psbI, rpl32, rps15, ycf1, ycf3, psaI-ycf4 and psbM-trnD) for species identification in Acronema. The molecular evolution of Acronema was relatively conservative that only one gene (petG) was found to be under positive selection (ω = 1.02489).

Discussion

The gene petG is one of the genes involved in the transmission of photosynthetic electron chains during photosynthesis, which plays a crucial role in the process of photosynthesis in plants. This is also a manifestation of the adaptive evolution of plants in high-altitude areas to the environment. In conclusion, our study provides novel insights into the plastome adaptive evolution, phylogeny, and taxonomy of genus Acronema.

Comparative analysis of complete chloroplast genomes of Synotis species (Asteraceae, Senecioneae) for identification and phylogenetic analysis

BACKGROUND

The Synotis (C. B. Clarke) C. Jeffrey & Y. L. Chen is an ecologically important genus of the tribe Senecioneae, family Asteraceae. Because most species of the genus bear similar morphology, traditional morphological identification methods are very difficult to discriminate them. Therefore, it is essential to develop a reliable and effective identification method for Synotis species. In this study, the complete chloroplast (cp.) genomes of four Synotis species, S. cavaleriei (H.Lév.) C. Jeffrey & Y.L. Chen, S. duclouxii (Dunn) C. Jeffrey & Y.L. Chen, S. nagensium (C.B. Clarke) C. Jeffrey & Y.L. Chen and S. erythropappa (Bureau & Franch.) C. Jeffrey & Y. L. Chen had been sequenced using next-generation sequencing technology and reported here.

RESULTS

These four cp. genomes exhibited a typical quadripartite structure and contained the large single-copy regions (LSC, 83,288 to 83,399 bp), the small single-copy regions (SSC, 18,262 to 18,287 bp), and the inverted repeat regions (IR, 24,837 to 24,842 bp). Each of the four cp. genomes encoded 134 genes, including 87 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (ycf1 and rps19). The highly variable regions (trnC-GCA-petN, ccsA-psaC, trnE-UUC-rpoB, ycf1, ccsA and petN) may be used as potential molecular barcodes. The complete cp. genomes sequence of Synotis could be used as the potentially effective super-barcode to accurately identify Synotis species. Phylogenetic analysis demonstrated that the four Synotis species were clustered into a monophyletic group, and they were closed to the Senecio, Crassocephalum and Dendrosenecio in tribe Senecioneae.

CONCLUSIONS

This study will be useful for further species identification, evolution, genetic diversity and phylogenetic studies within this genus Synotis and the tribe Senecioneae.

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

Comparative analysis of 12 water lily plastid genomes reveals genomic divergence and evolutionary relationships in early flowering plants

The aquatic plant Nymphaea, a model genus of the early flowering plant lineage Nymphaeales and family Nymphaeaceae, has been extensively studied. However, the availability of chloroplast genome data for this genus is incomplete, and phylogenetic relationships within the order Nymphaeales remain controversial. In this study, 12 chloroplast genomes of Nymphaea were assembled and analyzed for the first time. These genomes were 158,290-160,042 bp in size and contained 113 non-repeat genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. We also report on codon usage, RNA editing sites, microsatellite structures, and new repetitive sequences in this genus. Comparative genomics revealed that expansion and contraction of IR regions can lead to changes in the gene numbers. Additionally, it was observed that the highly variable regions of the chloroplast genome were mainly located in intergenic regions. Furthermore, the phylogenetic tree showed the order Nymphaeales was divided into three families, and the genus Nymphaea can be divided into five (or three) subgenera, with the subgenus Nymphaea being the oldest. The divergence times of nymphaealean taxa were analyzed, with origins of the order Nymphaeales and family Nymphaeaceae being about 194 and 131 million years, respectively. The results of the phylogenetic analysis and estimated divergence times will be useful for future evolutionary studies of basal angiosperm lineages.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-024-00242-0.

Polyphyly of Boehmeria (Urticaceae) congruent with plastome structural variation

Boehmeria is a taxonomically challenging group within the nettle family (Urticaceae). The polyphyly of the genus has been proposed by previous studies with respect to five genera (Debregeasia, Cypholophus, Sarcochlamys, Archiboehmeria, and Astrothalamus). Extensive homoplasy of morphological characters has made generic delimitation problematic. Previous studies in other plant groups suggest that plastome structural variations have the potential to provide characters useful in reconstructing evolutionary relationships. We aimed to test this across Boehmeria and its allied genera by mapping plastome structural variations onto a resolved strongly supported phylogeny. In doing so, we expanded the sampling of the plastome to include Cypholophus, Sarcochlamys, Archiboehmeria, and Astrothalamus for the first time. The results of our phylogenomic analyses provide strong support for Sarcochlamys as being more closely related to Leucosyke puya than to Boehmeria and for the clustering of Boehmeria s.l. into four subclades. The sizes of the plastomes in Boehmeria s.l. ranged from 142,627 bp to 170,958 bp. The plastomes recovered a typical quadripartite structure comprising 127~146 genes. We observe several obvious structural variations across the taxa such as gene loss and multiple gene duplication, inverted repeat (IR) contraction and wide expansions, and inversions. Moreover, we recover a trend for these variations that the early clades were relatively conserved in evolution, whereas the later diverging clades were variable. We propose that the structural variations documented may be linked to the adaptation of Boehmeria s.l. to a wide range of habitats, from moist broadleaf forests in Asia to xeric shrublands and deserts in Africa. This study confirms that variation in plastome gene loss/duplication, IR contraction/expansion, and inversions can provide evidence useful for the reconstruction of evolutionary relationships.

Characterization and phylogenetic analysis of the complete chloroplast genome of Carpesium lipskyi (Asteraceae, Inuleae)

The species of Carpesium lipskyi C.Winkl. 1998 is an important traditional Chinese medicine in China. In this study, the complete chloroplast (cp) genome of C. lipskyi was determined and analyzed. The result showed that the complete cp genome of C. lipskyi was 151,244 bp in length, consisting of a large single-copy (LSC) region of 82,908 bp, a small single-copy (SSC) region of 18,430 bp, and a pair of inverted repeats (IRs) of 24,953 bp. The overall GC content of the C. lipskyi is 37.68%. The species of C. lipskyi possessed 127 genes, including 83 protein-coding genes, 36 transfer RNA genes, and eight ribosomal RNA genes. The present study found that Inula is sister groups with the closest genetic relationship. The obtained knowledge could provide useful information for future phylogenetic, taxonomic, and evolutionary studies on Inuleae.

The complete chloroplast genome of Leonurus sibiricus Linnaeus (Labiatae, Leonurus Miller)

Leonurus sibiricus Linnaeus 1753, an annual or biennial herb found in northern China, Mongolia, and Russia, typically grows in stony, sandy grasslands, and pine forests. This study sequenced and reported the complete chloroplast genome of L. sibiricus for the first time. The entire circular genome measures 151,689 bp in length, with a GC content of 38.4%. A total of 133 genes were annotated, including 88 protein-coding genes, 37 tRNAs, and eight rRNAs. The genome exhibits a typical quadripartite structure, comprising a large single-copy (LSC 82,820 bp) region, a small single-copy (SSC 17,619 bp) region, and a pair of inverted repeat (IR 25,625 bp each) regions. Phylogenetic analysis using the maximum-likelihood method indicates that L. sibiricus is most closely related to L. japonicus Houttuyn. This study provides valuable genomic resources for further research on the phylogenetics and biodiversity of the genus Leonurus.

Characteristics of the chloroplast genome and genetic divergence of Tamarix hispida Willd. 1816 (Tamaricaceae)

Tamarix hispida Willd. 1816, a crucial native plant species in the arid desert region of northwestern China, plays a significant role in maintaining ecological stability. It is instrumental in addressing soil salinity-alkalinity and heavy metal pollution. This research aims to analyze the phylogenetic divergence pattern and evolutionary history of T. hispida by comparing chloroplast genome structures across different populations. Despite the minimal differences in chloroplast genome structure due to conserved genes and junction regions, sequencing was conducted using the Illumina NovaSeq platform to verify the historical evolutionary processes between different populations, followed by assembly and annotation. The results revealed that the T. hispida chloroplast genome is approximately 156,164-156,186 bp in length, with a quadripartite structure and 131 annotated genes. Phylogenetic analysis indicated two lineages within T. hispida, with a divergence time of 3.15 Ma. These findings emphasize the low genetic diversity in T. hispida and offer valuable insights into its evolutionary past. To effectively protect and manage this species, increased scientific research and monitoring of its genetic diversity are necessary. This study underscores the importance of comprehending the genetic mechanisms behind species divergence to develop informed conservation strategies.

Comparative Analysis of Plastome Sequences of Seven Tulipa L. (Liliaceae Juss.) Species from Section Kolpakowskianae Raamsd. Ex Zonn and Veldk

Tulipa L. is a genus of significant economic, environmental, and cultural importance in several parts of the world. The exact number of species in the genus remains uncertain due to inherent taxonomic challenges. We utilized next-generation sequencing technology to sequence and assemble the plastid genomes of seven Tulipa species collected in Kazakhstan and conducted a comparative analysis. The total number of annotated genes was 136 in all seven studied Tulipa species, 114 of which were unique, including 80 protein-coding, 30 tRNA, and 4 rRNA genes. Nine regions (petD, ndhH, ycf2-ycf3, ndhA, rpl16, clpP, ndhD-ndhF, rpoC2, and ycf1) demonstrated significant nucleotide variability, suggesting their potential as molecular markers. A total of 1388 SSRs were identified in the seven Tulipa plastomes, with mononucleotide repeats being the most abundant (60.09%), followed by dinucleotide (34.44%), tetranucleotide (3.90%), trinucleotide (1.08%), pentanucleotide (0.22%), and hexanucleotide (0.29%). The Ka/Ks values of the protein-coding genes ranged from 0 to 3.9286, with the majority showing values <1. Phylogenetic analysis based on a complete plastid genome and protein-coding gene sequences divided the species into three major clades corresponding to their subgenera. The results obtained in this study may contribute to understanding the phylogenetic relationships and molecular taxonomy of Tulipa species.

Chloroplast genome of plantago major, a medicinal plant in China

Plantago major (Plantaginaceae) is a medicinal plant in Chinese folk culture, known for its famous medicinal components such as plantagomain. In this study, we conducted genome sequencing of P. major using Illumina sequencing technology. The assembled complete chloroplast genome had a length of 165,044 bp, comprising a large single-copy regions (82,963 bp), a small single-copy regions (4,633 bp), and a pair of inverted repeat regions (38,724 bp). A total of 140 genes were detected, including 94 CDS, 38 tRNA, and 8 rRNA. Phylogenetic analysis revealed a close genetic relationship between P. major and P. rigida. These findings provide valuable data for a comprehensive understanding of the biological characteristics of P. major.

Complete chloroplast genome and phylogenetic analysis of Amorphophallus paeoniifolius (Araceae)

Amorphophallus paeoniifolius (Dennst.) Nicolson, 1885, often known as elephant foot yam, is a tropical tuber crop that originates from south-east Asia and belongs to the Araceae family. It is known for its high production potential and popularity as a medicinal plant. However, the phylogeny and genes for this species are still unavailable. In this study, the first complete chloroplast genome of A. paeoniifolius was reported and phylogenetic analysis was conducted with Araceae species. The chloroplast genome was 176,258 bp in length with 34.80% overall GC content and includes a large single-copy (LSC) region (93,951 bp), a small single-copy (SSC) region (15,013 bp), and a pair of inverted repeat (IRs) regions (33,647 bp). The chloroplast genome has 130 genes, which include 85 protein-coding genes, 37 tRNA genes, and eight rRNA genes. A maximum-likelihood (ML) phylogenetic analysis indicated that all Amorphophallus species formed a single monophyletic clade with a high bootstrap value and A. paeoniifolius was closely related to A. konjac, A. albus, A. krausei, and A. titanum. The chloroplast genome reported in this study will be useful for further taxonomic and evolutionary studies of Amorphophallus.

Phylogenomics identifies parents of naturally occurring tetraploid bananas

BACKGROUND

Triploid bananas are almost sterile. However, we succeeded in harvesting seeds from two edible triploid banana individuals (Genotype: ABB) in our conservation repository where various wild diploid bananas were also grown. The resulting rare offspring survived to seedling stages. DNA content analyses reveal that they are tetraploid. Since bananas contain maternally inherited plastids and paternally inherited mitochondria, we sequenced and assembled plastomes and mitogenomes of these seedlings to trace their hybridization history.

RESULTS

The coding sequences of both organellar genomic scaffolds were extracted, aligned, and concatenated for constructing phylogenetic trees. Our results suggest that these tetraploid seedlings be derived from hybridization between edible triploid bananas and wild diploid Musa balbisiana (BB) individuals. We propose that generating female triploid gametes via apomeiosis may allow the triploid maternal bananas to produce viable seeds.

CONCLUSIONS

Our study suggests a practical avenue towards expanding genetic recombination and increasing genetic diversity of banana breeding programs. Further cellular studies are needed to understand the fusion and developmental processes that lead to formation of hybrid embryos in banana reproduction, polyploidization, and evolution.

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

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.

Assembly and comparative analysis of the complete mitochondrial and chloroplast genome of Cyperus stoloniferus (Cyperaceae), a coastal plant possessing saline-alkali tolerance

BACKGROUND

Cyperus stoloniferus is an important species in coastal ecosystems and possesses economic and ecological value. To elucidate the structural characteristics, variation, and evolution of the organelle genome of C. stoloniferus, we sequenced, assembled, and compared its mitochondrial and chloroplast genomes.

RESULTS

We assembled the mitochondrial and chloroplast genomes of C. stoloniferus. The total length of the mitochondrial genome (mtDNA) was 927,413 bp, with a GC content of 40.59%. It consists of two circular DNAs, including 37 protein-coding genes (PCGs), 22 tRNAs, and five rRNAs. The length of the chloroplast genome (cpDNA) was 186,204 bp, containing 93 PCGs, 40 tRNAs, and 8 rRNAs. The mtDNA and cpDNA contained 81 and 129 tandem repeats, respectively, and 346 and 1,170 dispersed repeats, respectively, both of which have 270 simple sequence repeats. The third high-frequency codon (RSCU > 1) in the organellar genome tended to end at A or U, whereas the low-frequency codon (RSCU < 1) tended to end at G or C. The RNA editing sites of the PCGs were relatively few, with only 9 and 23 sites in the mtDNA and cpDNA, respectively. A total of 28 mitochondrial plastid DNAs (MTPTs) in the mtDNA were derived from cpDNA, including three complete trnT-GGU, trnH-GUG, and trnS-GCU. Phylogeny and collinearity indicated that the relationship between C. stoloniferus and C. rotundus are closest. The mitochondrial rns gene exhibited the greatest nucleotide variability, whereas the chloroplast gene with the greatest nucleotide variability was infA. Most PCGs in the organellar genome are negatively selected and highly evolutionarily conserved. Only six mitochondrial genes and two chloroplast genes exhibited Ka/Ks > 1; in particular, atp9, atp6, and rps7 may have undergone potential positive selection.

CONCLUSION

We assembled and validated the mtDNA of C. stoloniferus, which contains a 15,034 bp reverse complementary sequence. The organelle genome sequence of C. stoloniferus provides valuable genomic resources for species identification, evolution, and comparative genomic research in Cyperaceae.

Plastome structure, phylogeny and evolution of plastid genes in Reevesia (Helicteroideae, Malvaceae)

Reevesia is an eastern Asian-eastern North American disjunction genus in the family Malvaceae s.l. and comprises approximately 25 species. The relationships within the genus are not well understood. Here, 15 plastomes representing 12 Reevesia species were compared, with the aim of better understanding the species circumscription and phylogenetic relationships within the genus and among genera in the family Malvaceae s.l. The 11 newly sequenced plastomes range between 161,532 and 161, 945 bp in length. The genomes contain 114 unique genes, 18 of which are duplicated in the inverted repeats (IRs). Gene content of these plastomes is nearly identical. All the protein-coding genes are under purifying selection in the Reevesia plastomes compared. The top ten hypervariable regions, SSRs, and the long repeats identified are potential molecular markers for future population genetic and phylogenetic studies. Phylogenetic analysis based on the whole plastomes confirmed the monophyly of Reevesia and a close relationship with Durio (traditional Bombacaceae) in subfamily Helicteroideae, but not with the morphologically similar genera Pterospermum and Sterculia (both of traditional Sterculiaceae). Phylogenetic relationships within Reevesia suggested that two species, R. pubescens and R. thyrsoidea, as newly defined, are not monophyletic. Six taxa, R. membranacea, R. xuefengensis, R. botingensis, R. lofouensis, R. longipetiolata and R. pycnantha, are suggested to be recognized.

Plastid phylogenomics provides new insights into the systematics, diversification, and biogeography of Cymbidium (Orchidaceae)

Cymbidium (Orchidaceae: Epidendroideae), with around 60 species, is widely-distributed across Southeast Asia, providing a nice system for studying the processes that underlie patterns of biodiversity in the region. However, phylogenetic relationships of Cymbidium have not been well resolved, hampering investigations of species diversification and the biogeographical history of this genus. In this study, we construct a plastome phylogeny of 56 Cymbidium species, with four well-resolved major clades, which provides a framework for biogeographical and diversification rate analyses. Molecular dating and biogeographical analyses show that Cymbidium likely originated in the region spanning northern Indo-Burma to the eastern Himalayas during the early Miocene (∼21.10 Ma). It then rapidly diversified into four major clades in East Asia within approximately a million years during the middle Miocene. Cymbidium spp. migration to the adjacent regions (Borneo, Philippines, and Sulawesi) primarily occurred during the Pliocene-Pleistocene period. Our analyses indicate that the net diversification rate of Cymbidium has decreased since its origin, and is positively associated with changes in temperature and monsoon intensity. Favorable hydrothermal conditions brought by monsoon intensification in the early Miocene possibly contributed to the initial rapid diversification, after which the net diversification rate was reduced with the cooling climate after the middle Miocene. The transition from epiphytic to terrestrial habits may have enabled adaptation to cooler environments and colonization of northern niches, yet without a significant effect on diversification rates. This study provides new insights into how monsoon activity and temperature changes affected the diversification dynamics of plants in Southeast Asia.

Transition of survival strategies under global climate shifts in the grape family

Faced with environmental changes, plants may either move to track their ancestral niches or evolve to adapt to new niches. Vitaceae, the grape family, has evolved diverse adaptive traits facilitating a global expansion in wide-ranging habitats, making it ideal for investigating transition between move and evolve strategies and exploring the underlying mechanisms. Here we inferred the patterns of biogeographic diversification and trait evolution in Vitaceae based on a robust phylogeny with dense sampling including 495 species (~52% of Vitaceae species). Vitaceae probably originated from Asia-the diversity centre of extant genera and the major source of dispersals. Boundaries of the Eocene, Oligocene and Miocene were identified as turning points in shifting strategies. A significant decrease in move strategy was identified during the Oligocene, followed by increases in move and evolve. After the Miocene, evolve began to dominate, during which increased niche opportunities and key trait innovations played important roles.

Complete chloroplast genomes of eight Delphinium taxa (Ranunculaceae) endemic to Xinjiang, China: insights into genome structure, comparative analysis, and phylogenetic relationships

BACKGROUND

Delphinium L. represents a taxonomically intricate genus of significant phylogenetic and economic importance in Ranunculaceae. Despite the existence of few chloroplast genome datasets, a comprehensive understanding of genome structures and selective pressures within the genus remains unknown. Furthermore, several taxa in this genus are exclusively found in Xinjiang, China, a region renowned for its distribution and diversity of Chinese and Central Asian Delphinium species. Therefore, investigating the features of chloroplast genomes in this area will provide valuable insights into the evolutionary processes and phylogenetic relationships of the genus.

RESULTS

In this study, the eight newly completed chloroplast genomes are examined, ranging in length from 153,979 bp to 154,284 bp. Alongside these, analysing six previously reported taxa re-annotated in Delphinium, 111 unique genes are identified across all samples. Genome structure, distributions of simple sequence repeats and short dispersed repeats, as well as gene content are similar among these Delphinium taxa. Nine hypervariable intergenic spacers and protein coding regions, including ndhF-trnL(TAG), rpl16-intron, rpl33, rps15, rps18, trnK(TTT)-trnQ(TTG), trnP(TGG)-psaJ, trnT(GGT)-psbD and ycf1, are identified among 13 perennial Delphinium. Selective pressure and codon usage bias of all the plastid genes are performed within 14 Delphinium taxa. Phylogenetic analysis based on 14 Delphinium plastomes, alongside two Aconitum (Ranunculaceae) species serving as outgroup taxa, reveals the monophyletic nature of Delphinium. Our findings further discern Delphinium into two distinct clades: perennial species (clade I) and annual species (clade II). In addition, compared with the nrDNA ITS topology, cytological data and morphological characters, D. mollifolium and D. maackianum showed potential involvement in hybridization or polyploidization processes. Excluding these two species, the perennial Delphinium (clade I) exhibits a stronger consistency with the morphology-based system that utilized seed morphology.

CONCLUSION

This study represents the first comprehensive analysis of plastomic variations among Delphinium taxa, based on the examination of 14 complete plastomes. The chloroplast genome structure of Delphinium is similar to other angiosperms and possesses the typical quadripartite structure with the conserved genome arrangement and gene features. In addition, the variation of non-coding regions is larger than coding regions of the chloroplast genome. Through DNA sequence divergence across Delphinium plastomes and subsequent phylogenomic analyses ndhF-trnL(TAG) and ycf1 are identified as promising molecular markers. These highly variable loci held significant potential for future phylogenetic and phylogeographic studies on Delphinium. Our phylogenomic analyses based on the whole plastomes, concatenation of 132 unique intergenic spacer regions, concatenation of 77 unique protein-coding genes and nrDNA ITS, all support the monophyly of Delphinium and perennial taxa clusters together into one clade within this genus. These findings provide crucial data for systematic, phylogenomic and evolutionary research in the genus for future studies.

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).

Characterization and Comparative Analysis of Complete Chloroplast Genomes of Four Bromus (Poaceae, Bromeae) Species

Bromus (Poaceae Bromeae) is a forage grass with high adaptability and ecological and economic value. Here, we sequenced Bromus ciliatus, Bromus benekenii, Bromus riparius, and Bromus rubens chloroplast genomes and compared them with four previously described species. The genome sizes of Bromus species ranged from 136,934 bp (Bromus vulgaris) to 137,189 bp (Bromus ciliates, Bromus biebersteinii), with a typical quadripartite structure. The studied species had 129 genes, consisting of 83 protein-coding, 38 tRNA-coding, and 8 rRNA-coding genes. The highest GC content was found in the inverted repeat (IR) region (43.85-44.15%), followed by the large single-copy (LSC) region (36.25-36.65%) and the small single-copy (SSC) region (32.21-32.46%). There were 33 high-frequency codons, with those ending in A/U accounting for 90.91%. A total of 350 simple sequence repeats (SSRs) were identified, with single-nucleotide repeats being the most common (61.43%). A total of 228 forward and 141 palindromic repeats were identified. No reverse or complementary repeats were detected. The sequence identities of all sequences were very similar, especially with respect to the protein-coding and inverted repeat regions. Seven highly variable regions were detected, which could be used for molecular marker development. The constructed phylogenetic tree indicates that Bromus is a monophyletic taxon closely related to Triticum. This comparative analysis of the chloroplast genome of Bromus provides a scientific basis for species identification and phylogenetic studies.

Chloroplast Genomes Evolution and Phylogenetic Relationships of Caragana species

Caragana sensu lato (s.l.) includes approximately 100 species that are mainly distributed in arid and semi-arid regions. Caragana species are ecologically valuable for their roles in windbreaking and sand fixation. However, the taxonomy and phylogenetic relationships of the genus Caragana are still unclear. In this study, we sequenced and assembled the chloroplast genomes of representative species of Caragana and reconstructed robust phylogenetic relationships at the section level. The Caragana chloroplast genome has lost the inverted repeat region and wascategorized in the inverted repeat loss clade (IRLC). The chloroplast genomes of the eight species ranged from 128,458 bp to 135,401 bp and contained 110 unique genes. All the Caragana chloroplast genomes have a highly conserved structure and gene order. The number of long repeats and simple sequence repeats (SSRs) showed significant variation among the eight species, indicating heterogeneous evolution in Caragana. Selective pressure analysis of the genes revealed that most of the protein-coding genes evolved under purifying selection. The phylogenetic analyses indicated that each section forms a clade, except the section Spinosae, which was divided into two clades. This study elucidated the evolution of the chloroplast genome within the widely distributed genus Caragana. The detailed information obtained from this study can serve as a valuable resource for understanding the molecular dynamics and phylogenetic relationships within Caragana.

Structure, gene composition, divergence time and phylogeny analysis of the woody desert species Neltuma alba, Neltuma chilensis and Strombocarpa strombulifera

Neltuma alba (Algarrobo blanco), Neltuma chilensis (Algarrobo Chileno) and Strombocarpa strombulifera (Fortuna) are some of the few drought resistant trees and shrubs found in small highly fragmented populations, throughout the Atacama Desert. We reconstructed their plastid genomes using de novo assembly of paired-end reads from total genomic DNA. We found that the complete plastid genomes of N. alba and N. chilensis are larger in size compared to species of the Strombocarpa genus. The Strombocarpa species presented slightly more GC content than the Neltuma species. Therefore, we assume that Strombocarpa species have been exposed to stronger natural selection than Neltuma species. We observed high variation values in the number of cpSSRs (chloroplast simple sequence repeats) and repeated elements among Neltuma and Strombocarpa species. The p-distance results showed a low evolutionary divergence within the genus Neltuma, whereas a high evolutionary divergence was observed between Strombocarpa species. The molecular divergence time found in Neltuma and Strombocarpa show that these genera diverged in the late Oligocene. With this study we provide valuable information about tree species that provide important ecosystem services in hostile environments which can be used to determine these species in the geographically isolated communities, and keep the highly fragmented populations genetically healthy.

Plastid phylogenomics contributes to the taxonomic revision of taxa within the genus Sanicula L. and acceptance of two new members of the genus

Introduction

The genus Sanicula L. is a taxonomically complicated taxa within Apiaceae, as its high variability in morphology. Although taxonomists have performed several taxonomic revisions for this genus, the interspecific relationships and species boundaries have not been satisfactorily resolved, especially for those endemic to China. This study mainly focused on S. giraldii var. ovicalycina, S. tienmuensis var. pauciflora, and S. orthacantha var. stolonifera and also described two new members of the genus.

Methods

We newly sequenced sixteen plastomes from nine Sanicula species. Combined with eleven plastomes previously reported by us and one plastome downloaded, we performed a comprehensively plastid phylogenomics analysis of 21 Sanicula taxa.

Results and Discussion

The comparative results showed that 21 Sanicula plastomes in their structure and features were highly conserved and further justified that two new species were indeed members of Sanicula. Nevertheless, eleven mutation hotspot regions were still identified. Phylogenetic analyses based on plastome data and the ITS sequences strongly supported that these three varieties were clearly distant from three type varieties. The results implied that these three varieties should be considered as three independent species, which were further justified by their multiple morphological characters. Therefore, revising these three varieties into three independent species was reasonable and convincing. Moreover, we also identified and described two new Sanicula species (S. hanyuanensis and S. langaoensis) from Sichuan and Shanxi, China, respectively. Based on their distinct morphological characteristics and molecular phylogenetic analysis, two new species were included in Sanicula. In summary, our study impelled the revisions of Sanicula members and improved the taxonomic system of the genus.

Impatiensyingjingensis (Balsaminaceae), a new species from Sichuan, China

This study describes Impatiensyingjingensis X.Q. Song, B.N. Song & Biao Yang, sp. nov., a new species collected from the Yingjing area of the Giant Panda National Park. This new species is distributed at an altitude of 1400-2100 m, with a plant height of 30-130 cm. The flowers are purple-red or light purple red, with 3-9 flowers on each inflorescence and the dorsal auricle of the lateral united petals is thread-like and about 2 cm long, differing significantly from other species of Impatiens. Furthermore, molecular data, as well as micro-morphological evidence under SEM (of pollens), also support the establishment of the new species.

Comparative and phylogenetic analysis of the complete chloroplast genomes of 10 Artemisia selengensis resources based on high-throughput sequencing

BACKGROUND

Artemisia selengensis, classified within the genus Artemisia of the Asteraceae family, is a perennial herb recognized for its dual utility in culinary and medicinal domains. There are few studies on the chloroplast genome of A. selengensis, and the phylogeographic classification is vague, which makes phylogenetic analysis and evolutionary studies very difficult.

RESULTS

The chloroplast genomes of 10 A. selengensis in this study were highly conserved in terms of gene content, gene order, and gene intron number. The genome lengths ranged from 151,148 to 151,257 bp and were typical of a quadripartite structure with a total GC content of approximately 37.5%. The chloroplast genomes of all species encode 133 genes, including 88 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Due to the contraction and expansion of the inverted repeats (IR), the overlap of ycf1 and ndhF genes occurred at the inverted repeats B (IRB) and short single copy sequence (SSC) boundaries. According to a codon use study, the frequent base in the chloroplast genome of A. selengensis' third codon position was A/T. The number of SSR repeats was 42-44, most of which were single nucleotide A/T repeats. Sequence alignment analysis of the chloroplast genome showed that variable regions were mainly distributed in single copy regions, nucleotide diversity values of 0 to 0.009 were calculated by sliding window analysis, 8 mutation hotspot regions were detected, and coding regions were more conserved than non-coding regions. Analysis of non-synonymous substitution (Ka) and synonymous substitution (Ks) revealed that accD, rps12, petB, and atpF genes were affected by positive selection and no genes were affected by neutral selection. Based on the findings of the phylogenetic analysis, Artemisia selengensis was sister to the genus Artemisia Chrysanthemum and formed a monophyletic group with other Artemisia genera.

CONCLUSIONS

In this research, the present study systematically compared the chloroplast genomic features of A. selengensis and provided important information for the study of the chloroplast genome of A. selengensis and the evolutionary relationships among Asteraceae species.

Phylogeny, Genetic Diversity and Population Structure of Fritillaria cirrhosa and Its Relatives Based on Chloroplast Genome Data

Fritillaria cirrhosa and its relatives have been utilized in traditional Chinese medicine for many years and are under priority protection in China. Despite their medicinal and protective value, research on their phylogeny, genetic diversity, and divergence remains limited. Here, we investigate the chloroplast genome variation architecture of 46 samples of F. cirrhosa and its relatives collected from various regions, encompassing the majority of wild populations across diverse geographical areas. The results indicate abundant variations in 46 accessions including 1659 single-nucleotide polymorphisms and 440 indels. Six variable markers (psbJ, ndhD, ycf1, ndhG, trnT-trnL, and rpl32-trnL) were identified. Phylogenetic and network analysis, population structure analysis, and principal component analysis showed that the 46 accessions formed five clades with significant divergence, which were related to their geographical distribution. The regions spanning from the southern Hengduan Mountains to the Qinghai-Tibet Plateau exhibited the highest levels of genetic diversity. F. cirrhosa and its relatives may have suffered a genetic bottleneck and have a relatively low genetic diversity level. Moreover, geographical barriers and discrete patches may have accelerated population divergence. The study offers novel perspectives on the phylogeny, genetic diversity, and population structure of F. cirrhosa and its relatives, information that can inform conservation and utilization strategies in the future.