The evolution of the plastid chromosome in land plants: gene content, gene order, gene function

Plastid genome evolution in leafless members of the orchid subfamily Orchidoideae, with a focus on Degranvillea dermaptera

PREMISE

Leafless, heterotrophic plants are prime examples of organismal modification, the genomic consequences of which have received considerable interest. In particular, plastid genomes (plastomes) are being sequenced at a high rate, allowing continual refinement of conceptual models of reductive evolution in heterotrophs. However, numerous sampling gaps exist, hindering the ability to conduct comprehensive phylogenomic analyses in these plants.

METHODS

Using floral tissue from an herbarium specimen, we sequenced and analyzed the plastome of Degranvillea dermaptera, a rarely collected, leafless orchid species from South America about which little is known, including its phylogenetic affinities.

RESULTS

The plastome is the most reduced of those sequenced among the orchid subfamily Orchidoideae. In Degranvillea, it has lost the majority of genes found in leafy autotrophic species, is structurally rearranged, and has similar gene content to the most reduced plastomes among the orchids. We found strong evidence for the placement of Degranvillea within the subtribe Spiranthinae using models that explicitly account for heterotachy, or lineage-specific evolutionary rate variation over time. We further found evidence of relaxed selection on several genes and of correlations among substitution rates and several other "traits" of the plastome among leafless members of orchid subfamily Orchidoideae.

CONCLUSIONS

Our findings advance knowledge on the phylogenetic relationships and paths of plastid genome evolution among the orchids, which have experienced more independent transitions to heterotrophy than any other plant family. This study demonstrates the importance of herbarium collections in comparative genomics of poorly known species of conservation concern.

Characterization of the complete chloroplast genome of Rhodiola sachalinensis and comparative analysis with its congeneric plants

Rhodiola, belonging to the Crassulaceae family, is a perennial herbaceous plant genus. There are about 90 Rhodiola species worldwide, some of which have been reported to have medicinal properties. Rhodiola sachalinensis is a perennial medicinal herb within this genus and, in the present study, its chloroplast genome was sequenced, assembled, annotated and compared with 24 other Rhodiola species. The results obtained show that the chloroplast genome of R. sachalinensis is 151 595 bp long and has a CG content of 37.7%. The inverted repeats (IR) region of the Rhodiola chloroplast genome is the most conserved region, with the main differences being observed in the ycf1 and ndhF genes at the IRb-small single copy boundary, and rps19 and trnH genes at the IRa-large single copy boundary. Phylogenetic analysis showed that Rhodiola species form two major clades, and species with recorded medicinal properties, clustered together in one branch except for R. dumulosa. Within the genus, R. sachalinensis is most closely related to Rhodiola rosea, although comparative analyses showed that only R. sachalinensis and Rhodiola subopposita contained the psbZ gene, which encodes a highly conserved protein subunit of the Photosystem II core complex. Overall, the present study contributes to the understanding of the chloroplast genome of Rhodiola species, and provides a theoretical basis for the study of their genetic diversity and possible use as medicinal plants.

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.

Plant organellar genomes: much done, much more to do

Plastids and mitochondria are the only organelles that possess genomes of endosymbiotic origin. In recent decades, advances in sequencing technologies have contributed to a meteoric rise in the number of published organellar genomes, and have revealed greatly divergent evolutionary trajectories. In this review, we quantify the abundance and distribution of sequenced plant organellar genomes across the plant tree of life. We compare numerous genomic features between the two organellar genomes, with an emphasis on evolutionary trajectories, transfers, the current state of organellar genome editing by transcriptional activator-like effector nucleases (TALENs), transcription activator-like effector (TALE)-mediated deaminase, and clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas), as well as genetic transformation. Finally, we propose future research to understand these different evolutionary trajectories, and genome-editing strategies to promote functional studies and eventually improve organellar genomes.

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.

Insights into U-to-C RNA editing from the lycophyte Phylloglossum drummondii

The lycophyte Phylloglossum drummondii is the sole inhabitant of its genus in the Huperzioideae group and one of a small minority of plants which perform uridine to cytidine RNA editing. We assembled the P. drummondii chloroplast and mitochondrial genomes and used RNA sequence data to build a comprehensive profile of organellar RNA editing events. In addition to many C-to-U editing events in both organelles, we found just four U-to-C editing events in the mitochondrial transcripts cob, nad1, nad5 and rpl2. These events are conserved in related lycophytes in the genera Huperzia and Phlegmariurus. De novo transcriptomes for three of these lycophytes were assembled to search for putative U-to-C RNA editing enzymes. Four putative U-to-C editing factors could be matched to the four mitochondrial U-to-C editing sites. Due to the unusually few numbers of U-to-C RNA editing sites, P. drummondii and related lycophytes are useful models for studying this poorly understood mechanism.

Phylogeography of Morella nana: The Wumeng Mountains as a natural geographical isolation boundary on the Yunnan-Guizhou Plateau

The Yunnan-Guizhou Plateau (YGP) is characterized by the distinctive isolated habitat of the limestone Karst Islands and features the Wumeng Mountains, which divide the YGP into the two Plateaus of Yunnan and Guizhou. This study aimed to assess the effects of geographic isolation and past climate fluctuation on the distribution of flora in the YGP. To achieve this, we carried out the phylogeographical pattern and genetic structure based on chloroplast and nuclear ribosomal DNA sequence in relation to past (Last Glacial Maximum) and present distributions based on ecological niche modeling for Morella nana, an important wild plant resource and endemic to the YGP once considered a vulnerable species. The results suggest that the genetic and chlorotype network structures of M. nana are divided into at least two groups: cpDNA chlorotype H2 (or dominant nrDNA haplotypes h1 and h2), distributed primarily to the east of the Wumeng Mountains, and cpDNA chlorotypes H1 and H3-H10 (or dominant nrDNA haplotype h2 and h3), distributed to the west of the Wumeng Mountains. A deep genetic split was noted within the two groups to reach 25 steps, especially for the cpDNA fragment variation. This east-west divergence reveals the existence of a natural geographical isolation boundary in the form of the Wumeng Mountains, and supports the existence of at least two glacial refuges during the Quaternary glacial period, along with two genetic diversity center, and at least two large geographic protection units for the important species of M. nana. This study indicates that the phylogeographical pattern of M. nana can be attributed to geographic/environmental isolation caused by the Wumeng Mountains and climate fluctuation during the last glacial maximum, and proposes an effective strategy to protecting this important plant resource.

The complete chloroplast genome of the Syzygium buxifolium Hook. Et Arn.1833 and its phylogenetic analysis

Syzygium buxifolium. Hook. Et Arn.1833 is a member of the Myrtaceae family. This species is used in traditional Chinese medicines. It possesses numerous synonyms, reflecting the ambiguity in its taxonomy. The chloroplast genome has been widely used for species identification and phylogenetic analysis. Regrettably, there is a lack of information regarding the chloroplast genome of S. buxifolium. Here, we intend to obtain the chloroplast genome of S. buxifolium to resolve its classification problems. In particular, we utilized Illumina sequencing technology to sequence, GetOrganelle to assemble, and CPGAVAS2 to characterize the chloroplast genome of S. buxifolium. The chloroplast genome of S. buxifolium had a length of 158,581 bp and consisted of 111 unique genes, comprising 78 protein-coding genes, 29 transfer RNA (tRNA) genes, and four ribosomal RNA (rRNA) genes. In addition, we identified 86 Simple Sequence Repeats, 345 tandem repetitive sequences, and 34 dispersed repetitive sequences using modules implemented in CPGAVAS2. Lastly, we carried out phylogenetic analysis using Phylosuite. The results indicated a close relationship between S. buxifolium and S. grijsii. This study offers novel genetic data for the molecular identification and subsequent phylogenetic analysis of the Syzygium genus.

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.

Complete chloroplast genomes and phylogenetic relationships of Pedicularis chinensis and Pedicularis kansuensis

The complete cp genomes of Pedicularis chinensis (GenBank accession number: OQ587614) and Pedicularis kansuensis (GenBank accession number: OQ587613) were sequenced, assembled, and annotated. Their chloroplast (cp) genome lengths were 146,452 bp, and 146,852 bp, respectively; 120 and 116 genes were identified, comprising 75 and 72 protein-coding genes (PCGs), 37 and 36 transfer RNA (tRNA) genes, and 8 and 8 ribosomal RNA (rRNA) genes, for P. chinensis and P. kansuensis, respectively. A simple sequence repeat (SSR) analysis revealed that the repetitive sequences were mainly composed of mononucleotide repeats (A/T motif) and dinucleotide repeats (AT/TA motif). Comparative genomics identified several variant genes (rpl22, rps19, rpl12, ycf1, trnH, psbA, and ndhH) and variant regions (trnS-GGA, trnV-UAC, ndhJ-trnV, ycf4-cemA, ndhE-nhdG, and rpl32-trnL) with a high Pi, indicating the potential to serve as deoxyribo nucleic acid (DNA) barcodes for Pedicularis species identification. The results show that the cp genomes of P. chinensis and P. kansuensis were the same as those of other plants in Pedicularis, with different degrees of AT preference for codons. Large differences in the number of SSRs and the expansion of the inverted repeat (IR) region showed strong variability and interspecific differentiation between these two species and other species represented in the genus Pedicularis. A phylogenetic analysis showed that P. kansuensis had the closest relationship with P. oliveriana, and P. chinensis had the closest relationship with P. aschistorhyncha. These results will facilitate the study of the phylogenetic classification and interspecific evolution of Pedicularis plants.

Characterization of the complete chloroplast genome of the rare medicinal plant: Mandragora caulescens (Solanaceae)

In this study, we assembled high-quality chloroplast genomes of Mandragora caulescens through a reference-guided approach using high-throughput Illumina sequencing reads. The resulting chloroplast genome assembly displayed a typical quadripartite structural organization, comprising a large single-copy (LSC) region of 85,233 bp, two inverted repeat (IR) regions of 25,685 bp each, and a small single-copy (SSC) region of 18,207 bp. The chloroplast genome harbored 141 complete genes, and its overall GC content was 38.0%. In maximum-likelihood (ML) and Bayesian inference (BI) trees, the 19 Solanaceae species formed a monophyletic group, dividing into two main clades. M. caulescens and Nicandra physalodes formed a monophyletic group, suggesting a close relationship between the two species. The M. caulescens cp genome presented in this study lays a good foundation for further genetic and genomic studies of the Solanaceae.

The chloroplast genome sequences of Ipomoea alba and I. obscura (Convolvulaceae): genome comparison and phylogenetic analysis

Ipomoea species have diverse uses as ornamentals, food, and medicine. However, their genomic information is limited; I. alba and I. obscura were sequenced and assembled. Their chloroplast genomes were 161,353 bp and 159,691 bp, respectively. Both genomes exhibited a quadripartite structure, consisting of a pair of inverted repeat (IR) regions, which are separated by the large single-copy (LSC) and small single-copy (SSC) regions. The overall GC content was 37.5% for both genomes. A total of 104 and 93 simple sequence repeats, 50 large repeats, and 30 and 22 short tandem repeats were identified in the two chloroplast genomes, respectively. G and T were more preferred than C and A at the third base position based on the Parity Rule 2 plot analysis, and the neutrality plot revealed correlation coefficients of 0.126 and 0.105, indicating the influence of natural selection in shaping the codon usage bias in most protein-coding genes (CDS). Genome comparative analyses using 31 selected Ipomoea taxa from Thailand showed that their chloroplast genomes are rather conserved, but the presence of expansion or contraction of the IR region was identified in some of these Ipomoea taxa. A total of five highly divergent regions were identified, including the CDS genes accD, ndhA, and ndhF, as well as the intergenic spacer regions psbI-atpA and rpl32-ccsA. Phylogenetic analysis based on both the complete chloroplast genome sequence and CDS datasets of 31 Ipomoea taxa showed that I. alba is resolved as a group member for series (ser.) Quamoclit, which contains seven other taxa, including I. hederacea, I. imperati, I. indica, I. nil, I. purpurea, I. quamoclit, and I. × sloteri, while I. obscura is grouped with I. tiliifolia, both of which are under ser. Obscura, and is closely related to I. biflora of ser. Pes-tigridis. Divergence time estimation using the complete chloroplast genome sequence dataset indicated that the mean age of the divergence for Ipomoeeae, Argyreiinae, and Astripomoeinae, was approximately 29.99 Mya, 19.81 Mya, and 13.40 Mya, respectively. The node indicating the divergence of I. alba from the other members of Ipomoea was around 10.06 Mya, and the split between I. obscura and I. tiliifolia is thought to have happened around 17.13 Mya. The split between the I. obscura accessions from Thailand and Taiwan is thought to have taken place around 0.86 Mya.

The chloroplast genome of Chrozophora sabulosa Kar. & Kir. and its exploration in the evolutionary position uncertainty of genus Chrozophora

Chrozophora sabulosa Kar. & Kir. is a biennial herbaceous plant that belongs to the Euphorbiaceae family and has medicinal properties. This research aimed to identify the genetic characteristics and phylogenetic position of the Chrozophora genus within the Euphorbiaceae family. The evolutionary position of the Chrozophora genus was previously unknown due to insufficient research. Therefore, to determine the evolutionary link between C. sabulosa and other related species, we conducted a study using the NGS Illumina platform to sequence the C. sabulosa chloroplast (cp.) genome. The study results showed that the genome was 156,488 bp in length. It had a quadripartite structure consisting of two inverted repeats (IRb and IRa) of 24,649-bp, separated by an 87,696-bp LSC region and a 19,494-bp SSC region. The CP genome contained 113 unique genes, including four rRNA genes, 30 tRNA genes, and 79 CDS genes. In the second copy of the inverted repeat, there were 18 duplicated genes. The C. sabulosa lacks the petD, petB, rpl2, and rps16 intron. The analysis of simple sequence repeats (SSRs) revealed 93 SSR loci of 22 types and 78 oligonucleotide repeats of four kinds. The phylogenetic investigation showed that the Chrozophora genus evolved paraphyletically from other members of the Euphorbiaceae family. To support the phylogenetic findings, we selected species from the Euphorbiaceae and Phyllanthaceae families to compare with C. sabulosa for Ks and Ka substitution rates, InDels investigation, IR contraction and expansion, and SNPs analysis. The results of these comparative studies align with the phylogenetic findings. We identified six highly polymorphic regions shared by both families, which could be used as molecular identifiers for the Chrozophora genus (rpl33-rps18, rps18-rpl20, rps15-ycf1, ndhG-ndhI, psaI-ycf4, petA-psbJ). The cp. genome sequence of C. sabulosa reveals the evolution of plastid sequences in Chrozophora species. This is the first time the cp. genome of a Chrozophora genus has been sequenced, serving as a foundation for future sequencing of other species within the Chrozophoreae tribe and facilitating in-depth taxonomic research. The results of this research will also aid in identifying new Chrozophora species.

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.

Phylogenetic analysis and divergence time estimation of Lycium species in China based on the chloroplast genomes

BACKGROUND

Lycium is an economically and ecologically important genus of shrubs, consisting of approximately 70 species distributed worldwide, 15 of which are located in China. Despite the economic and ecological importance of Lycium, its phylogeny, interspecific relationships, and evolutionary history remain relatively unknown. In this study, we constructed a phylogeny and estimated divergence time based on the chloroplast genomes (CPGs) of 15 species, including subspecies, of the genus Lycium from China.

RESULTS

We sequenced and annotated 15 CPGs in this study. Comparative analysis of these genomes from these Lycium species revealed a typical quadripartite structure, with a total sequence length ranging from 154,890 to 155,677 base pairs (bp). The CPGs was highly conserved and moderately differentiated. Through annotation, we identified a total of 128-132 genes. Analysis of the boundaries of inverted repeat (IR) regions showed consistent positioning: the junctions of the IRb/LSC region were located in rps19 in all Lycium species, IRb/SSC between the ycf1 and ndhF genes, and SSC/IRa within the ycf1 gene. Sequence variation in the SSC region exceeded that in the IR region. We did not detect major expansions or contractions in the IR region or rearrangements or insertions in the CPGs of the 15 Lycium species. Comparative analyses revealed five hotspot regions in the CPG: trnR(UCU), atpF-atpH, ycf3-trnS(GGA), trnS(GGA), and trnL-UAG, which could potentially serve as molecular markers. In addition, phylogenetic tree construction based on the CPG indicated that the 15 Lycium species formed a monophyletic group and were divided into two typical subbranches and three minor branches. Molecular dating suggested that Lycium diverged from its sister genus approximately 17.7 million years ago (Mya) and species diversification within the Lycium species of China primarily occurred during the recent Pliocene epoch.

CONCLUSION

The divergence time estimation presented in this study will facilitate future research on Lycium, aid in species differentiation, and facilitate diverse investigations into this economically and ecologically important genus.

Piece and parcel of gymnosperm organellar genomes

MAIN CONCLUSION

Significant past, present, and potential future research into the organellar (plastid and mitochondrial) genomes of gymnosperms that can provide insight into the unknown origin and evolution of plants is highlighted. Gymnosperms are vascular seed plants that predominated the ancient world before their sister clade, angiosperms, took over during the Late Cretaceous. The divergence of gymnosperms and angiosperms took place around 300 Mya, with the latter evolving into the diverse group of flowering plants that dominate the plant kingdom today. Although gymnosperms have reportedly made some evolutionary innovations, the literature on their genome advances, particularly their organellar (plastid and mitochondrial) genomes, is relatively scattered and fragmented. While organellar genomes can shed light on plant origin and evolution, they are frequently overlooked, due in part to their limited contribution to gene expression and lack of evolutionary dynamics when compared to nuclear genomes. A better understanding of gymnosperm organellar genomes is critical because they reveal genetic changes that have contributed to their unique adaptations and ecological success, potentially aiding in plant survival, enhancement, and biodiversity conservation in the face of climate change. This review reveals significant information and gaps in the existing knowledge base of organellar genomes in gymnosperms, as well as the challenges and research needed to unravel their complexity.

A comparative plastome approach enhances the assessment of genetic variation in the Melilotus genus

BACKGROUND

Melilotus, a member of the Fabaceae family, is a pivotal forage crop that is extensively cultivated in livestock regions globally due to its notable productivity and ability to withstand abiotic stress. However, the genetic attributes of the chloroplast genome and the evolutionary connections among different Melilotus species remain unresolved.

RESULTS

In this study, we compiled the chloroplast genomes of 18 Melilotus species and performed a comprehensive comparative analysis. Through the examination of protein-coding genes, we successfully established a robust phylogenetic tree for these species. This conclusion is further supported by the phylogeny derived from single-nucleotide polymorphisms (SNPs) across the entire chloroplast genome. Notably, our findings revealed that M. infestus, M. siculus, M. sulcatus, and M. speciosus formed a distinct subgroup within the phylogenetic tree. Additionally, the chloroplast genomes of these four species exhibit two shared inversions. Moreover, inverted repeats were observed to have reemerged in six species within the IRLC. The distribution patterns of single-nucleotide polymorphisms (SNPs) and insertions/deletions (InDels) within protein-coding genes indicated that ycf1 and ycf2 accumulated nonconservative alterations during evolutionary development. Furthermore, an examination of the evolutionary rate of protein-coding genes revealed that rps18, rps7, and rpl16 underwent positive selection specifically in Melilotus.

CONCLUSIONS

We present a comparative analysis of the complete chloroplast genomes of Melilotus species. This study represents the most thorough and detailed exploration of the evolution and variability within the genus Melilotus to date. Our study provides valuable chloroplast genomic information for improving phylogenetic reconstructions and making biogeographic inferences about Melilotus and other Papilionoideae species.

Guiding seed movement: environmental heterogeneity drives genetic differentiation in Plathymenia reticulata, providing insights for restoration

Forest and landscape restoration is one of the main strategies for overcoming the environmental crisis. This activity is particularly relevant for biodiversity-rich areas threatened by deforestation, such as tropical forests. Efficient long-term restoration requires understanding the composition and genetic structure of native populations, as well as the factors that influence these genetic components. This is because these populations serve as the seed sources and, therefore, the gene reservoirs for areas under restoration. In the present study, we investigated the influence of environmental, climatic and spatial distance factors on the genetic patterns of Plathymenia reticulata, aiming to support seed translocation strategies for restoration areas. We collected plant samples from nine populations of P. reticulata in the state of Bahia, Brazil, located in areas of Atlantic Forest and Savanna, across four climatic types, and genotyped them using nine nuclear and three chloroplast microsatellite markers. The populations of P. reticulata evaluated generally showed low to moderate genotypic variability and low haplotypic diversity. The populations within the Savanna phytophysiognomy showed values above average for six of the eight evaluated genetic diversity parameters. Using this classification based on phytophysiognomy demonstrated a high predictive power for genetic differentiation in P. reticulata. Furthermore, the interplay of climate, soil and geographic distance influenced the spread of alleles across the landscape. Based on our findings, we propose seed translocation, taking into account the biome, with restricted use of seed sources acquired or collected from the same environment as the areas to be restored (Savanna or Atlantic Forest).

Comparative Analysis of Complete Chloroplast Genomes of Rubus in China: Hypervariable Regions and Phylogenetic Relationships

With more than 200 species of native Rubus, China is considered a center of diversity for this genus. Due to a paucity of molecular markers, the phylogenetic relationships for this genus are poorly understood. In this study, we sequenced and assembled the plastomes of 22 out of 204 Chinese Rubus species (including varieties) from three of the eight sections reported in China, i.e., the sections Chamaebatus, Idaeobatus, and Malachobatus. Plastomes were annotated and comparatively analyzed with the inclusion of two published plastomes. The plastomes of all 24 Rubus species were composed of a large single-copy region (LSC), a small single-copy region (SSC), and a pair of inverted repeat regions (IRs), and ranged in length from 155,464 to 156,506 bp. We identified 112 unique genes, including 79 protein-coding genes, 29 transfer RNAs, and four ribosomal RNAs. With highly consistent gene order, these Rubus plastomes showed strong collinearity, and no significant changes in IR boundaries were noted. Nine divergent hotspots were identified based on nucleotide polymorphism analysis: trnH-psbA, trnK-rps16, rps16-trnQ-psbK, petN-psbM, trnT-trnL, petA-psbJ, rpl16 intron, ndhF-trnL, and ycf1. Based on whole plastome sequences, we obtained a clearer phylogenetic understanding of these Rubus species. All sampled Rubus species formed a monophyletic group; however, sections Idaeobatus and Malachobatus were polyphyletic. These data and analyses demonstrate the phylogenetic utility of plastomes for systematic research within Rubus.

Comparative chloroplast genomics, phylogenetic relationships and molecular markers development of Aglaonema commutatum and seven green cultivars of Aglaonema

Aglaonema commutatum is a famous species in the Aglaonema genus, which has important ornamental and economic value. However, its chloroplast genome information and phylogenetic relationships among popular green cultivars of Aglaonema in southern China have not been reported. Herein, chloroplast genomes of one variety of A. commutatum and seven green cultivars of Aglaonema, namely, A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Sapphire', 'Silver Queen', 'Snow White', 'White Gem', and 'White Horse Prince', were sequenced and assembled for comparative analysis and phylogeny. These eight genomes possessed a typical quadripartite structure that consisted of a LSC region (90,799-91,486 bp), an SSC region (20,508-21,137 bp) and a pair of IR regions (26,661-26,750 bp). Each genome contained 112 different genes, comprising 79 protein-coding genes, 29 tRNA genes and 4 rRNA genes. The gene orders, GC contents, codon usage frequency, and IR/SC boundaries were highly conserved among these eight genomes. Long repeats, SSRs, SNPs and indels were analyzed among these eight genomes. Comparative analysis of 15 Aglaonema chloroplast genomes identified 7 highly variable regions, including trnH-GUG-exon1-psbA, trnS-GCU-trnG-UCC-exon1, trnY-GUA-trnE-UUC, psbC-trnS-UGA, trnF-GAA-ndhJ, ccsA-ndhD, and rps15-ycf1-D2. Reconstruction of the phylogenetic trees based on chloroplast genomes, strongly supported that Aglaonema was a sister to Anchomanes, and that the Aglaonema genus was classified into two sister clades including clade I and clade II, which corresponded to two sections, Aglaonema and Chamaecaulon, respectively. One variety and five cultivars, including A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Silver Queen', 'Snow White', and 'White Horse Prince', were classified into clade I; and the rest of the two cultivars, including 'Sapphire' and 'White Gem', were classified into clade II. Positive selection was observed in 34 protein-coding genes at the level of the amino acid sites among 77 chloroplast genomes of the Araceae family. Based on the highly variable regions and SSRs, 4 DNA markers were developed to differentiate the clade I and clade II in Aglaonema. In conclusion, this study provided chloroplast genomic resources for Aglaonema, which were useful for its classification and phylogeny.

Testing plastomes and nuclear ribosomal DNA sequences as the next-generation DNA barcodes for species identification and phylogenetic analysis in Acer

BACKGROUND

Acer is a taxonomically intractable and speciose genus that contains over 150 species. It is challenging to distinguish Acer species only by morphological method due to their abundant variations. Plastome and nuclear ribosomal DNA (nrDNA) sequences are recommended as powerful next-generation DNA barcodes for species discrimination. However, their efficacies were still poorly studied. The current study will evaluate the application of plastome and nrDNA in species identification and perform phylogenetic analyses for Acer.

RESULT

Based on a collection of 83 individuals representing 55 species (c. 55% of Chinese species) from 13 sections, our barcoding analyses demonstrated that plastomes exhibited the highest (90.47%) species discriminatory power among all plastid DNA markers, such as the standard plastid barcodes matK + rbcL + trnH-psbA (61.90%) and ycf1 (76.19%). And the nrDNA (80.95%) revealed higher species resolution than ITS (71.43%). Acer plastomes show abundant interspecific variations, however, species identification failure may be due to the incomplete lineage sorting (ILS) and chloroplast capture resulting from hybridization. We found that the usage of nrDNA contributed to identifying those species that were unidentified by plastomes, implying its capability to some extent to mitigate the impact of hybridization and ILS on species discrimination. However, combining plastome and nrDNA is not recommended given the cytonuclear conflict caused by potential hybridization. Our phylogenetic analysis covering 19 sections (95% sections of Acer) and 128 species (over 80% species of this genus) revealed pervasive inter- and intra-section cytonuclear discordances, hinting that hybridization has played an important role in the evolution of Acer.

CONCLUSION

Plastomes and nrDNA can significantly improve the species resolution in Acer. Our phylogenetic analysis uncovered the scope and depth of cytonuclear conflict in Acer, providing important insights into its evolution.

Comparative and phylogenetic analyses of Loranthaceae plastomes provide insights into the evolutionary trajectories of plastome degradation in hemiparasitic plants

BACKGROUND

The lifestyle transition from autotrophy to heterotrophy often leads to extensive degradation of plastomes in parasitic plants, while the evolutionary trajectories of plastome degradation associated with parasitism in hemiparasitic plants remain poorly understood. In this study, phylogeny-oriented comparative analyses were conducted to investigate whether obligate Loranthaceae stem-parasites experienced higher degrees of plastome degradation than closely related facultative root-parasites and to explore the potential evolutionary events that triggered the 'domino effect' in plastome degradation of hemiparasitic plants.

RESULTS

Through phylogeny-oriented comparative analyses, the results indicate that Loranthaceae hemiparasites have undergone varying degrees of plastome degradation as they evolved towards a heterotrophic lifestyle. Compared to closely related facultative root-parasites, all obligate stem-parasites exhibited an elevated degree plastome degradation, characterized by increased downsizing, gene loss, and pseudogenization, thereby providing empirical evidence supporting the theoretical expectation that evolution from facultative parasitism to obligate parasitism may result in a higher degree of plastome degradation in hemiparasites. Along with infra-familial divergence in Loranthaceae, several lineage-specific gene loss/pseudogenization events occurred at deep nodes, whereas further independent gene loss/pseudogenization events were observed in shallow branches.

CONCLUSIONS

The findings suggest that in addition to the increasing levels of nutritional reliance on host plants, cladogenesis can be considered as another pivotal evolutionary event triggering the 'domino effect' in plastome degradation of hemiparasitic plants. These findings provide new insights into the evolutionary trajectory of plastome degradation in hemiparasitic plants.

Characterization of the Plastid Genomes of Four Caroxylon Thunb. Species from Kazakhstan

The family Chenopodiaceae Vent. (Amaranthaceae s.l.) is known for its taxonomic complexity, comprising species of significant economic and ecological importance. Despite its significance, the availability of plastid genome data for this family remains limited. This study involved assembling and characterizing the complete plastid genomes of four Caroxylon Thunb. species within the tribe Salsoleae s.l., utilizing next-generation sequencing technology. We compared genome features, nucleotide diversity, and repeat sequences and conducted a phylogenetic analysis of ten Salsoleae s.l. species. The size of the plastid genome varied among four Caroxylon species, ranging from 150,777 bp (C. nitrarium) to 151,307 bp (C. orientale). Each studied plastid genome encoded 133 genes, including 114 unique genes. This set of genes includes 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Eight divergent regions (accD, atpF, matK, ndhF-ndhG, petB, rpl20-rpl22, rpoC2, and ycf3) were identified in ten Salsoleae s.l. plastid genomes, which could be potential DNA-barcoding markers. Additionally, 1106 repeat elements were detected, consisting of 814 simple sequence repeats, 92 tandem repeats, 88 forward repeats, 111 palindromic repeats, and one reverse repeat. The phylogenetic analysis provided robust support for the relationships within Caroxylon species. These data represent a valuable resource for future phylogenetic studies within the genus.

Comparative chloroplast genomes of Argentina species: genome evolution and phylogenomic implications

The genus Argentina Hill belongs to the tribe Potentilleae Sweet and contains approximately 75 species predominantly distributed in the Sino-Himalayan region and the Malesian archipelago. So far we have less knowledge on the phylogenetic relationships within Argentina owing to limited sampling of Argentina taxa or gene fragments in previous studies. Moreover, to date there is no phylogenetic study on Argentina from the perspective of comparative chloroplast (cp) genomics. Here we performed comparative genomic analyses on the cp genomes of 39 accessions representing 18 taxa of Argentina. The Argentina cp genomes presented the typical quadripartite structure, with the sizes ranging from 155 096 bp to 157 166 bp. The 39 Argentina cp genomes contained a set of 112 unique genes, comprising four ribosomal RNA (rRNA) genes, 30 transfer RNA (tRNA) genes, as well as 78 protein-coding genes (PCGs). The cp genome organization, gene content and order in Argentina were highly conserved, but some visible divergences were present in IR/SC boundary regions. Ten regions (trnH-GUG-psbA, trnG-GCC-trnfM-CAU, trnD-GUC-trnY-GUA, rpl32-trnL-UAG, atpH-atpI, rps16-trnQ-UUG, trnS-GCU-trnG-UCC, ndhF-rpl32, trnR-UCU-atpA, and accD-psaI) were identified as excellent candidate DNA markers for future studies on species identification, population genetics and phylogeny of Argentina. Our results indicated that Argentina is monophyletic. In the current sampling, the A. smithiana - A. anserina clade was sister to the remainder of Argentina. Our results corroborated the previous taxonomic treatments to transfer A. phanerophlebia and A. micropetala from the genus Sibbaldia L. to Argentina. Our results showed close relationships among A. stenophylla, A. microphylla, A. taliensis, and A. tatsienluensis, congruent with previous studies based on the morphology of these species. Twenty-six genes (rps3, rps15, rps16, rps19, rpl16, rpl20, rpl22, rpoA, rpoB, rpoC1, rpoC2, atpA, atpF, psbB, psbF, ndhA, ndhB, ndhC, ndhD, ndhF, rbcL, accD, ccsA, matK, ycf1, ycf2) were with sites under positive selection, and adaptive evolution of these genes might have played crucial roles in Argentina species adaptation to the harsh mountain environment. This study will facilitate future work on taxonomy, phylogenetics, and adaptive evolution of Argentina.

Unprecedented variation pattern of plastid genomes and the potential role in adaptive evolution in Poales

BACKGROUND

The plastid is the photosynthetic organelle in plant cell, and the plastid genomes (plastomes) are generally conserved in evolution. As one of the most economically and ecologically important order of angiosperms, Poales was previously documented to exhibit great plastomic variation as an order of photoautotrophic plants.

RESULTS

We acquired 93 plastomes, representing all the 16 families and 5 major clades of Poales to reveal the extent of their variation and evolutionary pattern. Extensive variation including the largest one in monocots with 225,293 bp in size, heterogeneous GC content, and a wide variety of gene duplication and loss were revealed. Moreover, rare occurrences of three inverted repeat (IR) copies in angiosperms and one IR loss were observed, accompanied by short IR (sIR) and small direct repeat (DR). Widespread structural heteroplasmy, diversified inversions, and unusual genomic rearrangements all appeared in Poales, occasionally within a single species. Extensive repeats in the plastomes were found to be positively correlated with the observed inversions and rearrangements. The variation all showed a "small-large-moderate" trend along the evolution of Poales, as well as for the sequence substitution rate. Finally, we found some positively selected genes, mainly in C4 lineages, while the closely related lineages of those experiencing gene loss tended to have undergone more relaxed purifying selection.

CONCLUSIONS

The variation of plastomes in Poales may be related to its successful diversification into diverse habitats and multiple photosynthetic pathway transitions. Our order-scale analyses revealed unusual evolutionary scenarios for plastomes in the photoautotrophic order of Poales and provided new insights into the plastome evolution in angiosperms as a whole.

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

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

Decoding Evolution of Rubioideae: Plastomes Reveal Sweet Secrets of Codon Usage, Diagnostides, and Superbarcoding

Galium genus belongs to the Rubiaceae family, which consists of approximately 14,000 species. In comparison to its well-known relatives, the plastomes of the Galium genus have not been explored so far. The plastomes of this genus have a typical, quadripartite structure, but differ in gene content, since the infA gene is missing in Galium palustre and Galium trfidum. An evaluation of the effectiveness of using entire chloroplast genome sequences as superbarcodes for accurate plant species identification revealed the high potential of this method for molecular delimitation within the genus and tribe. The trnE-UUC-psbD region showed the biggest number of diagnostides (diagnostic nucleotides) which might be new potential barcodes, not only in Galium, but also in other closely related genera. Relative synonymous codon usage (RSCU) appeared to be connected with the phylogeny of the Rubiaceae family, showing that during evolution, plants started preferring specific codons over others.

Insights into comparative genomics, structural features, and phylogenetic relationship of species from Eurasian Aster and its related genera (Asteraceae: Astereae) based on complete chloroplast genome

Aster L. is an economically and phylogenetically important genus in the tribe Astereae. Here, the complete plastomes of the eight Aster species were assembled and characterized using next-generation sequencing datasets. The results indicated the complete plastomes of Aster had a quadripartite structure. These genomes were 152,045-152,729 bp in length and contained 132-133 genes, including 87 protein-coding genes, 37-38 tRNA genes, and eight rRNA genes. Expansion or contraction of inverted repeat regions and forward, palindromic, complement, and reverse repeats were detected in the eight Aster species. Additionally, our analyses showed the richest type of simple sequence repeats was A/T mononucleotides, and 14 highly variable regions were discovered by analyzing the border regions, sequence divergence, and hotspots. Phylogenetic analyses indicated that 27 species in Astereae were clustered into six clades, i.e., A to D, North American, and outgroup clades, and supported that the genera Heteropappus, Kalimeris, and Heteroplexis are nested within Aster. The results indicated the clades B to D might be considered as genera. Divergence time estimate showed the clades A, B, C, and D diverged at 23.15 Mya, 15.13 Mya, 24.29 Mya, and 21.66 Mya, respectively. These results shed light on the phylogenetic relationships of Aster and provided new information on species identification of Aster and its related genera.

Complete Chloroplast Genome Sequence Structure and Phylogenetic Analysis of Kohlrabi (Brassica oleracea var. gongylodes L.)

Kohlrabi is an important swollen-stem cabbage variety belonging to the Brassicaceae family. However, few complete chloroplast genome sequences of this genus have been reported. Here, a complete chloroplast genome with a quadripartite cycle of 153,364 bp was obtained. A total of 132 genes were identified, including 87 protein-coding genes, 37 transfer RNA genes and eight ribosomal RNA genes. The base composition analysis showed that the overall GC content was 36.36% of the complete chloroplast genome sequence. Relative synonymous codon usage frequency (RSCU) analysis showed that most codons with values greater than 1 ended with A or U, while most codons with values less than 1 ended with C or G. Thirty-five scattered repeats were identified and most of them were distributed in the large single-copy (LSC) region. A total of 290 simple sequence repeats (SSRs) were found and 188 of them were distributed in the LSC region. Phylogenetic relationship analysis showed that five Brassica oleracea subspecies were clustered into one group and the kohlrabi chloroplast genome was closely related to that of B. oleracea var. botrytis. Our results provide a basis for understanding chloroplast-dependent metabolic studies and provide new insight for understanding the polyploidization of Brassicaceae species.

The complete chloroplast genome and phylogentic results support the species position of Swertia banzragczii and Swertia marginata (Gentianaceae) in Mongolia

BACKGROUND

Swertia banzragczii and S. marginata are important medicinal species in Mongolia. However, their taxonomic positions and genetic backgrounds remain unknown. In this study, we explored the complete chloroplast genomes and DNA barcoding of these species and compared them with those of closely related species within the subgenus to determine their taxonomic positions and phylogenetic relationships.

RESULT

The chloroplast genomes of S. banzragczii and S. marginata encoded 114 genes, including 80 protein-coding genes, 30 tRNA genes, and 4 rRNA genes. Among them, 16 genes contained a single intron, and 2 genes had two introns. Closely related species had a conserved genome structure and gene content. Only differences in genome length were noticed, which were caused by the expansion and contraction of the inverted repeat (IR) region and loss of exons in some genes. The trnH-GUG-psbA and trnD-GUC-trnY-GUA intergenic regions had high genetic diversity within Swertia plastomes. Overall, S. banzragczii and S. marginata are true species and belong to the subgenus Swertia.

CONCLUSIONS

These results provide valuable genetic and morphological information on rare and subendemic Swertia species in Mongolia, which can be used for further advanced studies on the Swertia genus.

Phylogenomics and plastomics offer new evolutionary perspectives on Kalanchoideae (Crassulaceae)

BACKGROUND AND AIMS

Kalanchoideae is one of three subfamilies within Crassulaceae and contains four genera. Despite previous efforts, the phylogeny of Kalanchoideae remains inadequately resolved with persistent issues including low support, unstructured topologies and polytomies. This study aimed to address two central objectives: (1) resolving the pending phylogenetic questions within Kalanchoideae by using organelle-scale 'barcodes' (plastomes) and nuclear data; and (2) investigating interspecific diversity patterns among Kalanchoideae plastomes.

METHODS

To explore the plastome evolution in Kalanchoideae, we newly sequenced 38 plastomes representing all four constituent genera (Adromischus, Cotyledon, Kalanchoe and Tylecodon). We performed comparative analyses of plastomic features, including GC and gene contents, gene distributions at the IR (inverted repeat) boundaries, nucleotide divergence, plastomic tRNA (pttRNA) structures and codon aversions. Additionally, phylogenetic inferences were inferred using both the plastomic dataset (79 genes) and nuclear dataset (1054 genes).

KEY RESULTS

Significant heterogeneities were observed in plastome lengths among Kalanchoideae, strongly correlated with LSC (large single copy) lengths. Informative diversities existed in the gene content at SSC/IRa (small single copy/inverted repeat a), with unique patterns individually identified in Adromischus leucophyllus and one major Kalanchoe clade. The ycf1 gene was assessed as a shared hypervariable region among all four genera, containing nine lineage-specific indels. Three pttRNAs exhibited unique structures specific to Kalanchoideae and the genera Adromischus and Kalanchoe. Moreover, 24 coding sequences revealed a total of 41 lineage-specific unused codons across all four constituent genera. The phyloplastomic inferences clearly depicted internal branching patterns in Kalanchoideae. Most notably, by both plastid- and nuclear-based phylogenies, our research offers the first evidence that Kalanchoe section Eukalanchoe is not monophyletic.

CONCLUSIONS

This study conducted comprehensive analyses on 38 newly reported Kalanchoideae plastomes. Importantly, our results not only reconstructed well-resolved phylogenies within Kalanchoideae, but also identified highly informative unique markers at the subfamily, genus and species levels. These findings significantly enhance our understanding of the evolutionary history of Kalanchoideae.

The size diversity of the Pteridaceae family chloroplast genome is caused by overlong intergenic spacers

BACKGROUND

While the size of chloroplast genomes (cpDNAs) is often influenced by the expansion and contraction of inverted repeat regions and the enrichment of repeats, it is the intergenic spacers (IGSs) that appear to play a pivotal role in determining the size of Pteridaceae cpDNAs. This provides an opportunity to delve into the evolution of chloroplast genomic structures of the Pteridaceae family. This study added five Pteridaceae species, comparing them with 36 published counterparts.

RESULTS

Poor alignment in the non-coding regions of the Pteridaceae family was observed, and this was attributed to the widespread presence of overlong IGSs in Pteridaceae cpDNAs. These overlong IGSs were identified as a major factor influencing variations in cpDNA size. In comparison to non-expanded IGSs, overlong IGSs exhibited significantly higher GC content and were rich in repetitive sequences. Species divergence time estimations suggest that these overlong IGSs may have already existed during the early radiation of the Pteridaceae family.

CONCLUSIONS

This study reveals new insights into the genetic variation, evolutionary history, and dynamic changes in the cpDNA structure of the Pteridaceae family, providing a fundamental resource for further exploring its evolutionary research.

Variations and reduction of plastome are associated with the evolution of parasitism in Convolvulaceae

Parasitic lifestyle can often relax the constraint on the plastome, leading to gene pseudogenization and loss, and resulting in diverse genomic structures and rampant genome degradation. Although several plastomes of parasitic Cuscuta have  been reported, the evolution of parasitism in the family Convolvulaceae which is linked to structural variations and reduction of plastome has not been well investigated. In this study, we assembled and collected 40 plastid genomes belonging to 23 species representing four subgenera of Cuscuta and ten species of autotrophic Convolvulaceae. Our findings revealed nine types of structural variations and six types of inverted repeat (IR) boundary variations in the plastome of Convolvulaceae spp. These structural variations were associated with the shift of parasitic lifestyle, and IR boundary shift, as well as the abundance of long repeats. Overall, the degradation of Cuscuta plastome proceeded gradually, with one clade exhibiting an accelerated degradation rate. We observed five stages of gene loss in Cuscuta, including NAD(P)H complex → PEP complex → Photosynthesis-related → Ribosomal protein subunits → ATP synthase complex. Based on our results, we speculated that the shift of parasitic lifestyle in early divergent time promoted relaxed selection on plastomes, leading to the accumulation of microvariations, which ultimately resulted in the plastome reduction. This study provides new evidence towards a better understanding of plastomic evolution, variation, and reduction in the genus Cuscuta.

Chloroplast genome structure analysis of Equisetum unveils phylogenetic relationships to ferns and mutational hotspot region

Equisetum is one of the oldest extant group vascular plants and is considered to be the key to understanding vascular plant evolution. Equisetum is distributed almost all over the world and has a high degree of adaptability to different environments. Despite the fossil record of horsetails (Equisetum, Equisetaceae) dating back to the Carboniferous, the phylogenetic relationship of this genus is not well, and the chloroplast evolution in Equisetum remains poorly understood. In order to fill this gap, we sequenced, assembled, and annotated the chloroplast genomes of 12 species of Equisetum, and compared them to 13 previously published vascular plants chloroplast genomes to deeply examine the plastome evolutionary dynamics of Equisetum. The chloroplast genomes have a highly conserved quadripartite structure across the genus, but these chloroplast genomes have a lower GC content than other ferns. The size of Equisetum plastomes ranges from 130,773 bp to 133,684 bp and they encode 130 genes. Contraction/expansion of IR regions and the number of simple sequences repeat regions underlie large genomic variations in size among them. Comparative analysis revealed we also identified 13 divergence hotspot regions. Additionally, the genes accD and ycf1 can be used as potential DNA barcodes for the identification and phylogeny of the genus Equisetum. Twelve photosynthesis-related genes were specifically selected in Equisetum. Comparative genomic analyses implied divergent evolutionary patterns between Equisetum and other ferns. Phylogenomic analyses and molecular dating revealed a relatively distant phylogenetic relationship between Equisetum and other ferns, supporting the division of pteridophyte into Lycophytes, Equisetaceae and ferns. The results show that the chloroplast genome can be used to solve phylogenetic problems within or between Equisetum species, and also provide genomic resources for the study of Equisetum systematics and evolution.

Chloroplast genomes of Caragana tibetica and Caragana turkestanica: structures and comparative analysis

BACKGROUND

The genus Caragana encompasses multiple plant species that possess medicinal and ecological value. However, some species of Caragana are quite similar in morphology, so identifying species in this genus based on their morphological characteristics is considerably complex. In our research, illumina paired-end sequencing was employed to investigate the genetic organization and structure of Caragana tibetica and Caragana turkestanica, including the previously published chloroplast genome sequence of 7 Caragana plants.

RESULTS

The lengths of C. tibetica and C. turkestanica chloroplast genomes were 128,433 bp and 129,453 bp, respectively. The absence of inverted repeat sequences in these two species categorizes them under the inverted repeat loss clade (IRLC). They encode 110 and 111 genes (4 /4 rRNA genes, 30 /31tRNA genes, and 76 /76 protein-coding genes), respectively. Comparison of the chloroplast genomes of C. tibetica and C. turkestanica with 7 other Caragana species revealed a high overall sequence similarity. However, some divergence was observed between certain intergenic regions (matK-rbcL, psbD-psbM, atpA-psbI, and etc.). Nucleotide diversity (π) analysis revealed the detection of five highly likely variable regions, namely rps2-atpI, accD-psaI-ycf4, cemA-petA, psbN-psbH and rpoA-rps11. Phylogenetic analysis revealed that C. tibetica's sister species is Caragana jubata, whereas C. turkestanica's closest relative is Caragana arborescens.

CONCLUSIONS

The present study provides worthwhile information about the chloroplast genomes of C. tibetica and C. turkestanica, which aids in the identification and classification of Caragana species.

Complete plastid genome of Iris orchioides and comparative analysis with 19 Iris plastomes

Iris is a cosmopolitan genus comprising approximately 280 species distributed throughout the Northern Hemisphere. Although Iris is the most diverse group in the Iridaceae, the number of taxa is debatable owing to various taxonomic issues. Plastid genomes have been widely used for phylogenetic research in plants; however, only limited number of plastid DNA markers are available for phylogenetic study of the Iris. To understand the genomic features of plastids within the genus, including its structural and genetic variation, we newly sequenced and analyzed the complete plastid genome of I. orchioides and compared it with those of 19 other Iris taxa. Potential plastid markers for phylogenetic research were identified by computing the sequence divergence and phylogenetic informativeness. We then tested the utility of the markers with the phylogenies inferred from the markers and whole-plastome data. The average size of the plastid genome was 152,926 bp, and the overall genomic content and organization were nearly identical among the 20 Iris taxa, except for minor variations in the inverted repeats. We identified 10 highly informative regions (matK, ndhF, rpoC2, ycf1, ycf2, rps15-ycf, rpoB-trnC, petA-psbJ, ndhG-ndhI and psbK-trnQ) and inferred a phylogeny from each region individually, as well as from their concatenated data. Remarkably, the phylogeny reconstructed from the concatenated data comprising three selected regions (rpoC2, ycf1 and ycf2) exhibited the highest congruence with the phylogeny derived from the entire plastome dataset. The result suggests that this subset of data could serve as a viable alternative to the complete plastome data, especially for molecular diagnoses among closely related Iris taxa, and at a lower cost.

Comparative Analysis of Chloroplast Genomes for the Genus Manglietia Blume (Magnoliaceae): Molecular Structure and Phylogenetic Evolution

Manglietia Blume, belonging to the Magnoliaceae family and mainly distributed in tropical and subtropical regions of Asia, has great scientific and economic value. In this study, we employed next-generation sequencing followed by de novo assembly to investigate the adaptive evolution of Manglietia using plastid genetic information. We newly sequenced the complete or nearly complete plastomes of four Manglietia species (Manglietia aromatica, Manglietia calcarea, Manglietia kwangtungensis, and Manglietia glauca) and conducted comparative analysis with seventeen published plastomes to examine the evolutionary pattern within this genus. The plastomes of these five newly sequenced Manglietia species range from 157,093 bp (M. calcarea2) to 160,493 bp (M. kwangtungensis), all exhibiting circular structures when mapped. Nucleotide diversity was observed across the plastomes, leading us to identify 13 mutational hotspot regions, comprising eight intergenic spacer regions and five gene regions. Our phylogenetic analyses based on 77 protein-coding genes generated phylogenetic relationships with high support and resolution for Manglietia. This genus can be divided into three clades, and the previously proposed infrageneric classifications are not supported by our studies. Furthermore, the close affinity between M. aromatica and M. calcarea is supported by the present work, and further studies are necessary to conclude the taxonomic treatment for the latter. These results provide resources for the comparative plastome, breeding, and plastid genetic engineering of Magnoliaceae and flowering plants.

Full-length transcriptome analysis of Ophioglossum vulgatum: effects of experimentally identified chloroplast gene clusters on expression and evolutionary patterns

Genes with similar or related functions in chloroplasts are often arranged in close proximity, forming clusters on chromosomes. These clusters are transcribed coordinated to facilitate the expression of genes with specific function. Our previous study revealed a significant negative correlation between the chloroplast gene expression level of the rare medicinal fern Ophioglossum vulgatum and its evolutionary rates as well as selection pressure. Therefore, in this study, we employed a combination of SMRT and Illumina sequencing technology to analyze the full-length transcriptome sequencing of O. vulgatum for the first time. In particular, we experimentally identified gene clusters based on transcriptome data and investigated the effects of chloroplast gene clustering on expression and evolutionary patterns. The results revealed that the total sequenced data volume of the full-length transcriptome of O. vulgatum amounted to 71,950,652,163 bp, and 110 chloroplast genes received transcript coverage. Nine different types of gene clusters were experimentally identified in their transcripts. The chloroplast cluster genes may cause a decrease in non-synonymous substitution rate and selection pressure, as well as a reduction in transversion rate, transition rate, and their ratio. While expression levels of chloroplast cluster genes in leaf, sporangium, and stem would be relatively elevated. The Mann-Whitney U test indicated statistically significant in the selection pressure, sporangia and leaves groups (P < 0.05). We have contributed novel full-length transcriptome data resources for ferns, presenting new evidence on the effects of chloroplast gene clustering on expression land evolutionary patterns, and offering new theoretical support for transgenic research through gene clustering.

Comparative plastome analysis of the sister genera Ceratocephala and Myosurus (Ranunculaceae) reveals signals of adaptive evolution to arid and aquatic environments

BACKGROUND

Expansion and contraction of inverted repeats can cause considerable variation of plastid genomes (plastomes) in angiosperms. However, little is known about whether structural variations of plastomes are associated with adaptation to or occupancy of new environments. Moreover, adaptive evolution of angiosperm plastid genes remains poorly understood. Here, we sequenced the complete plastomes for four species of xerophytic Ceratocephala and hydrophytic Myosurus, as well as Ficaria verna. By an integration of phylogenomic, comparative genomic, and selection pressure analyses, we investigated evolutionary patterns of plastomes in Ranunculeae and their relationships with adaptation to dry and aquatic habitats.

RESULTS

Owing to the significant contraction of the boundary of IRA/LSC towards the IRA, plastome sizes and IR lengths of Myosurus and Ceratocephala are smaller within Ranunculeae. Compared to other Ranunculeae, the Myosurus plastome lost clpP and rps16, one copy of rpl2 and rpl23, and one intron of rpoC1 and rpl16, and the Ceratocephala plastome added an infA gene and lost one copy of rpl2 and two introns of clpP. A total of 11 plastid genes (14%) showed positive selection, two genes common to Myosurus and Ceratocephala, seven in Ceratocephala only, and two in Myosurus only. Four genes showed strong signals of episodic positive selection. The rps7 gene of Ceratocephala and the rpl32 and ycf4 genes of Myosurus showed an increase in the rate of variation close to 3.3 Ma.

CONCLUSIONS

The plastomic structure variations as well as the positive selection of two plastid genes might be related to the colonization of new environments by the common ancestor of Ceratocephala and Myosurus. The seven and two genes under positive selection might be related to the adaptation to dry and aquatic habitats in Ceratocephala and Myosurus, respectively. Moreover, intensified aridity and frequent sea-level fluctuations, as well as global cooling, might have favored an increased rate of change in some genes at about 3.3 Ma, associated with adaptation to dry and aquatic environments, respectively. These findings suggest that changing environments might have influenced structural variations of plastomes and fixed new mutations arising on some plastid genes owing to adaptation to specific habitats.

The newly assembled chloroplast genome of Aeluropus littoralis: molecular feature characterization and phylogenetic analysis with related species

Aeluropus littoralis, a halophyte grass, is widely distributed from the Mediterranean to the Indian subcontinent through the Mongolian Gobi. This model halophyte has garnered increasing attention owing to its use as forage and its high tolerance to environmental stressors. The chloroplast genomes of many plants have been extensively examined for molecular, phylogenetic and transplastomic applications. However, no published research on the A. littoralis chloroplast (cp) genome was discovered. Here, the entire chloroplast genome of A. littoralis was assembled implementing accurate long-read sequences. The entire chloroplast genome, with an estimated length of 135,532 bp (GC content: 38.2%), has a quadripartite architecture and includes a pair of inverted repeat (IR) regions, IRa and IRb (21,012 bp each), separated by a large and a small single-copy regions (80,823 and 12,685 bp, respectively). The features of A. littoralis consist of 133 genes that synthesize 87 peptides, 38 transfer RNAs, and 8 ribosomal RNAs. Of these genes, 86 were unique, whereas 19 were duplicated in IR regions. Additionally, a total of forty-six simple sequence repeats, categorized into 32-mono, four-di, two-tri, and eight-tetranucleotides, were discovered. Furthermore, ten sets of repeats greater than 20 bp were located primarily in the LSC region. Evolutionary analysis based on chloroplast sequence data revealed that A. littoralis with A. lagopoides and A. sinensis belong to the Aeluropodinae subtribe, which is a sister to the Eleusininae in the tribe Cynodonteae and the subfamily Chloridoideae. This subfamily belongs to the PACMAD clade, which contains the majority of the C4 photosynthetic plants in the Poaceae. The newly constructed A. littoralis cp genome offers valuable knowledge for DNA barcoding, phylogenetic, transplastomic research, and other biological studies.

Three complete chloroplast genomes from two north American Rhus species and phylogenomics of Anacardiaceae

BACKGROUND

The suamc genus Rhus (sensu stricto) includes two subgenera, Lobadium (ca. 25 spp.) and Rhus (ca. 10 spp.). Their members, R. glabra and R. typhina (Rosanae: Sapindales: Anacardiaceae), are two economic important species. Chloroplast genome information is of great significance for the study of plant phylogeny and taxonomy.

RESULTS

The three complete chloroplast genomes from two Rhus glabra and one R. typhina accessions were obtained with a total of each about 159k bp in length including a large single-copy region (LSC, about 88k bp), a small single-copy regions (SSC, about 19k bp) and a pair of inverted repeats regions (IRa/IRb, about 26k bp), to form a canonical quadripartite structure. Each genome contained 88 protein-coding genes, 37 transfer RNA genes, eight ribosomal RNA genes and two pseudogenes. The overall GC content of the three genomes all were same (37.8%), and RSCU values showed that they all had the same codon prefers, i.e., to use codon ended with A/U (93%) except termination codon. Three variable hotspots, i.e., ycf4-cemA, ndhF-rpl32-trnL and ccsA-ndhD, and a total of 152-156 simple sequence repeats (SSR) were identified. The nonsynonymous (Ka)/synonymous (Ks) ratio was calculated, and cemA and ycf2 genes are important indicators of gene evolution. The phylogenetic analyses of the family Anacardiaceae showed that the eight genera were grouped into three clusters, and supported the monophyly of the subfamilies and all the genera. The accessions of five Rhus species formed four clusters, while, one individual of R. typhina grouped with the R. glabra accessions instead of clustering into the two other individuals of R. typhina in the subgenus Rhus, which showed a paraphyletic relationship.

CONCLUSIONS

Comparing the complete chloroplast genomes of the Rhus species, it was found that most SSRs were A/T rich and located in the intergenic spacer, and the nucleotide divergence exhibited higher levels in the non-coding region than in the coding region. The Ka/Ks ratio of cemA gene was > 1 for species collected in America, while it was < 1 for other species in China, which dedicated that the Rhus species from North America and East Asia have different evolutionary pressure. The phylogenetic analysis of the complete chloroplast genome clarified the Rhus placement and relationship. The results obtained in this study are expected to provide valuable genetic resources to perform species identification, molecular breeding, and intraspecific diversity of the Rhus species.

Complete chloroplast genome and comparison of herbicides toxicity on Aeschynomene indica (Leguminosae) in upland direct-seeding paddy field

BACKGROUND

Indian jointvetch (Aeschynomene indica) is a common and pernicious weed found in the upland direct-seeding rice fields in the lower reaches of the Yangtze River in China. However, there are few reports on the degree of harm, genetic characteristics, and management methods of this weed. The purpose of this study is to clarify the harm of Indian jointvetch to upland direct-seeding rice, analyze the genetic characteristics of this weed based on chloroplast genomics and identify its related species, and screen herbicides that are effective in managing this weed in upland direct-seeding rice fields.

RESULTS

In a field investigation in upland direct-seeding rice paddies in Shanghai and Jiangsu, we determined that the plant height and maximum lateral distance of Indian jointvetch reached approximately 134.2 cm and 57.9 cm, respectively. With Indian jointvetch present at a density of 4/m2 and 8/m2, the yield of rice decreased by approximately 50% and 70%, respectively. We further obtained the first assembly of the complete chloroplast (cp.) genome sequence of Indian jointvetch (163,613 bp). There were 161 simple sequence repeats, 166 long repeats, and 83 protein-encoding genes. The phylogenetic tree and inverted repeat region expansion and contraction analysis based on cp. genomes demonstrated that species with closer affinity to A. indica included Glycine soja, Glycine max, and Sesbania cannabina. Moreover, a total of 3281, 3840, and 3838 single nucleotide polymorphisms were detected in the coding sequence regions of the cp. genomes of S. cannabina voucher IBSC, G. soja, and G. max compared with the A. indica sequence, respectively. A greenhouse pot experiment indicated that two pre-emergence herbicides, saflufenacil and oxyfluorfen, and two post-emergence herbicides, florpyrauxifen-benzyl and penoxsulam, can more effectively manage Indian jointvetch than other common herbicides in paddy fields. The combination of these two types of herbicides is recommended for managing Indian jointvetch throughout the entire growth period of upland direct-seeding rice.

CONCLUSIONS

This study provides molecular resources for future research focusing on the identification of the infrageneric taxa, phylogenetic resolution, and biodiversity of Leguminosae plants, along with recommendations for reliable management methods to control Indian jointvetch.

A comparative analysis of plastome evolution in autotrophic Piperales

PREMISE

Many plastomes of autotrophic Piperales have been reported to date, describing a variety of differences. Most studies focused only on a few species or a single genus, and extensive, comparative analyses have not been done. Here, we reviewed publicly available plastome reconstructions for autotrophic Piperales, reanalyzed publicly available raw data, and provided new sequence data for all previously missing genera. Comparative plastome genomics of >100 autotrophic Piperales were performed.

METHODS

We performed de novo assemblies to reconstruct the plastomes of newly generated sequence data. We used Sanger sequencing and read mapping to verify the assemblies and to bridge assembly gaps. Furthermore, we reconstructed the phylogenetic relationships as a foundation for comparative plastome genomics.

RESULTS

We identified a plethora of assembly and annotation issues in published plastome data, which, if unattended, will lead to an artificial increase of diversity. We were able to detect patterns of missing and incorrect feature annotation and determined that the inverted repeat (IR) boundaries were the major source for erroneous assembly. Accounting for the aforementioned issues, we discovered relatively stable junctions of the IRs and the small single-copy region (SSC), whereas the majority of plastome variations among Piperales stems from fluctuations of the boundaries of the IR and the large single-copy (LSC) region.

CONCLUSIONS

This study of all available plastomes of autotrophic Piperales, expanded by new data for previously missing genera, highlights the IR-LSC junctions as a potential marker for discrimination of various taxonomic levels. Our data indicates a pseudogene-like status for cemA and ycf15 in various Piperales. Based on a review of published data, we conclude that incorrect IR-SSC boundary identification is the major source for erroneous plastome assembly. We propose a gold standard for assembly and annotation of high-quality plastomes based on de novo assembly methods and appropriate references for gene annotation.

Pan-chloroplast genomes for accession-specific marker development in Hibiscus syriacus

Hibiscus syriacus L. is a renowned ornamental plant. We constructed 95 chloroplast genomes of H. syriacus L. cultivars using a short-read sequencing platform (Illumina) and a long-read sequencing platform (Oxford Nanopore Technology). The following genome assembly, we delineate quadripartite structures encompassing large single-copy, small single-copy, and inverted repeat (IRa and IRb) regions, from 160,231 bp to 161,041 bp. Our comprehensive analyses confirmed the presence of 79 protein-coding genes, 30 tRNA genes, and 4 rRNA genes in the pan-chloroplast genome, consistent with prior research on the H. syriacus chloroplast genome. Subsequent pangenome analysis unveiled widespread genome sequence conservation alongside unique cultivar-specific variant patterns consisting of 193 single-nucleotide polymorphisms and 61 insertions or deletions. The region containing intra-species variant patterns, as identified in this study, has the potential to develop accession-specific molecular markers, enhancing precision in cultivar classification. These findings are anticipated to drive advancements in breeding strategies, augment biodiversity, and unlock the agricultural potential inherent in H. syriacus.

Chloroplast Genome Provides Insights into Molecular Evolution and Species Relationship of Fleabanes (Erigeron: Tribe Astereae, Asteraceae) in the Juan Fernández Islands, Chile

Erigeron represents the third largest genus on the Juan Fernández Islands, with six endemic species, five of which occur exclusively on the younger Alejandro Selkirk Island with one species on both islands. While its continental sister species is unknown, Erigeron on the Juan Fernández Islands appears to be monophyletic and most likely evolved from South American progenitor species. We characterized the complete chloroplast genomes of five Erigeron species, including accessions of E. fernandezia and one each from Alejandro Selkirk and Robinson Crusoe Islands, with the purposes of elucidating molecular evolution and phylogenetic relationships. We found highly conserved chloroplast genomes in size, gene order and contents, and further identified several mutation hotspot regions. In addition, we found two positively selected chloroplast genes (ccsA and ndhF) among species in the islands. The complete plastome sequences confirmed the monophyly of Erigeron in the islands and corroborated previous phylogenetic relationships among species. New findings in the current study include (1) two major lineages, E. turricola-E. luteoviridis and E. fernandezia-E. ingae-E. rupicola, (2) the non-monophyly of E. fernandezia occurring on the two islands, and (3) the non-monophyly of the alpine species E. ingae complex.

Characterization of Firmiana danxiaensis plastomes and comparative analysis of Firmiana: insight into its phylogeny and evolution

BACKGROUND

Firmiana danxiaensis is a critically endangered and ecologically important tree currently only found in four locations in Danxia or Karst habitats in northern Guangdong Province, China. The specialized habitat preference makes it an ideal model species for study of adaptive evolution. Meanwhile, the phylogenetic relationships of F. danxiaensis in four locations under two landforms are unclear. Therefore, we sequenced its complete chloroplast (cp.) genomes and conducted comprehensive interspecific and intrageneric plastome studies.

RESULTS

The F. danxiaensis plastomes in four locations showed a typical quadripartite and circular structure that ranged from 160,832 to 161,206 bp in size, with 112 unique genes encoded. Comparative genomics showed that the plastomes of F. danxiaensis were relatively conserved with high similarity of genome organization, gene number, GC content and SSRs. While the genomes revealed higher biased codon preferences in Karst habitat than those in Danxia habitats. Eighteen and 11 divergent hotpots were identified at interspecific and intrageneric levels for species identification and further phylogenetic studies. Seven genes (clpP, accD, ccsA, ndhH, rpl20, rpoC2, and rps4) were under positive selection and may be related to adaptation. Phylogenetic analysis revealed that F. danxiaensis is sister to F. major and F. simplex. However, the interspecific relationships are not consistent with the habitat types.

CONCLUSIONS

The characteristics and interspecific relationship of F. danxiaensis plastomes provide new insights into further integration of geographical factors, environmental factors, and genetic variations on the genomic study of F. danxiaensis. Together, our study will contribute to the study of species identification, population genetics, and conservation biology of F. danxiaensis.

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

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

The evolution of ephemeral flora in Xinjiang, China: insights from plastid phylogenomic analyses of Brassicaceae

BACKGROUND

The ephemeral flora of northern Xinjiang, China, plays an important role in the desert ecosystems. However, the evolutionary history of this flora remains unclear. To gain new insights into its origin and evolutionary dynamics, we comprehensively sampled ephemeral plants of Brassicaceae, one of the essential plant groups of the ephemeral flora.

RESULTS

We reconstructed a phylogenetic tree using plastid genomes and estimated their divergence times. Our results indicate that ephemeral species began to colonize the arid areas in north Xinjiang during the Early Miocene and there was a greater dispersal of ephemeral species from the surrounding areas into the ephemeral community of north Xinjiang during the Middle and Late Miocene, in contrast to the Early Miocene or Pliocene periods.

CONCLUSIONS

Our findings, together with previous studies, suggest that the ephemeral flora originated in the Early Miocene, and species assembly became rapid from the Middle Miocene onwards, possibly attributable to global climate changes and regional geological events.