Loss of all ndh genes as determined by sequencing the entire chloroplast genome of the black pine Pinus thunbergii

Plastome phylogenomics unravels the evolutionary relationships and biogeographic history of Chloranthaceae

BACKGROUND

Chloranthaceae is a pantropical family of flowering plants distributed mainly across Central and South America, East Asia, and the Pacific, encompassing approximately 73 species belonging to four extant genera-Ascarina, Chloranthus, Hedyosmum, and Sarcandra. As one of the most ancient lineages of extant angiosperms, Chloranthaceae holds substantial value in traditional medicine globally and offers major insights into the evolutionary history of flowering plants. However, phylogenetic relationships within this family remain partially resolved, and its origin continues to be debated. We here sequenced, assembled, annotated, and compared the chloroplast genomes (plastomes) of 22 Chloranthaceae species to investigate their plastome evolution and reconstructed the family's phylogeny by using both plastome and nuclear ribosomal DNA (nrDNA). Additionally, we employed the divergence-extinction-cladogenesis model to infer the ancestral area and historical migration patterns of Chloranthaceae.

RESULTS

Results The plastomes of Chloranthaceae exhibit a typical quadripartite structure, with genome sizes ranging from 157,449 bp to 159,218 bp and encoding 113 unique genes, including 17-18 genes duplicated in inverted repeat (IR) regions. Interestingly, the genus Hedyosmum demonstrated notable IR expansions, resulting in trnHGUG duplication in five species. A comparative analysis unveiled substantial variations in simple sequence repeats (SSRs) and tandem repeats across different species, whereas interspersed repeats remained relatively stable. Correlation analyses revealed positive correlations between plastome size and both IR length and SSR count. Furthermore, our phylogenetic reconstruction based on plastome and nrDNA (18 S-ITS1-5.8 S-ITS2-26 S) data provided robust support for the intergeneric relationships, confirming the monophyly of each genus and resolving long-standing uncertainties in the intrageneric relationships of Chloranthus and Hedyosmum. Notably, H. orientale, the sole Hedyosmum species found in Asia, was phylogenetically embedded within the genus, rather than occupying a basal position contrary to previous results. Our biogeographical reconstruction, incorporating both extant and fossil evidence, supports a Laurasian origin of Chloranthaceae and its subsequent dispersal to the tropics.

CONCLUSION

Conclusions This study provides the most comprehensive plastome-based phylogenomic analysis of Chloranthaceae to date, offering novel insights into its plastome evolution, resolving long-standing ambiguities in interspecific and phylogenetic relationships and reconstructing the probable biogeographic origin and historical migrations of this family. The proposed resources establish a robust framework for future evolutionary investigations within the realm of Chloranthaceae.

Chloroplast Genome Sequencing and Comparative Analysis of Six Medicinal Plants of Polygonatum

The genus Polygonatum boasts abundant germplasm resources and comprises numerous species. Among these, medicinal plants of this genus, which have a long history, have garnered attention of scholars. This study sequenced and analyzed the chloroplast genomes of six species of Polygonatum medicinal plants (P. zanlanscianense, P. kingianum, P. sibiricum, P. cyrtonema, P. filipes, and P. odoratum, respectively) to explore their interspecific relationships. The sequence length (154, 578-155, 807 bp) and genome structure were conserved among the six Polygonatum species, with a typical tetrad structure. Among the 127-131 genes contained in the genomes, 84-85 are protein-coding genes, 37-38 are transfer RNA genes, and 6-8 are ribosomal RNA genes. The genomes contained 64-76 simple sequence repeats (SSRs) and 36-62 long repetitive sequences. Codon bias patterns tended to use codons ending in A/T. In 30 types of codons with RSCU > 1, 93.3% ended in A/T of the six species. Twenty-one highly variable plastid regions were identified in the chloroplast genomes of the six medicinal plants. Furthermore, a phylogenetic analysis encompassing these and 53 other chloroplast genomes of Polygonatum species revealed that P. cyrtonema, P. odoratum, and P. filipes clustered together on one clade, whereas P. kingianum and P. zanlanscianense formed separate clades. Notably, P. sibiricum emerged as a standalone clade, and our phylogenetic tree reinforces the classification of P. sibiricum as forming a monophyly. This study provides a novel basis for intragenus taxonomy and DNA barcoding molecular identification within the genus Polygonatum medicinal plants.

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.

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.

Plastid phylogenomics and cytonuclear discordance in Rubioideae, Rubiaceae

In this study of evolutionary relationships in the subfamily Rubioideae (Rubiaceae), we take advantage of the off-target proportion of reads generated via previous target capture sequencing projects based on nuclear genomic data to build a plastome phylogeny and investigate cytonuclear discordance. The assembly of off-target reads resulted in a comprehensive plastome dataset and robust inference of phylogenetic relationships, where most intratribal and intertribal relationships are resolved with strong support. While the phylogenetic results were mostly in agreement with previous studies based on plastome data, novel relationships in the plastid perspective were also detected. For example, our analyses of plastome data provide strong support for the SCOUT clade and its sister relationship to the remaining members of the subfamily, which differs from previous results based on plastid data but agrees with recent results based on nuclear genomic data. However, several instances of highly supported cytonuclear discordance were identified across the Rubioideae phylogeny. Coalescent simulation analysis indicates that while ILS could, by itself, explain the majority of the discordant relationships, plastome introgression may be the better explanation in some cases. Our study further indicates that plastomes across the Rubioideae are, with few exceptions, highly conserved and mainly conform to the structure, gene content, and gene order present in the majority of the flowering plants.

Complete Chloroplast Genome of Alternanthera sessilis and Comparative Analysis with Its Congeneric Invasive Weed Alternanthera philoxeroides

Alternanthera sessilis is considered the closest relative to the invasive weed Alternanthera philoxeroides in China, making it an important native species for studying the invasive mechanisms and adaptations of A. philoxeroides. Chloroplasts play a crucial role in a plant's environmental adaptation, with their genomes being pivotal in the evolution and adaptation of both invasive and related species. However, the chloroplast genome of A. sessilis has remained unknown until now. In this study, we sequenced and assembled the complete chloroplast genome of A. sessilis using high-throughput sequencing. The A. sessilis chloroplast genome is 151,935 base pairs long, comprising two inverted repeat regions, a large single copy region, and a small single copy region. This chloroplast genome contains 128 genes, including 8 rRNA-coding genes, 37 tRNA-coding genes, 4 pseudogenes, and 83 protein-coding genes. When compared to the chloroplast genome of the invasive weed A. philoxeroides and other Amaranthaceae species, we observed significant variations in the ccsA, ycf1, and ycf2 regions in the A. sessilis chloroplast genome. Moreover, two genes, ccsA and accD, were found to be undergoing rapid evolution due to positive selection pressure. The phylogenetic trees were constructed for the Amaranthaceae family, estimating the time of independent species formation between A. philoxeroides and A. sessilis to be approximately 3.5186-8.8242 million years ago. These findings provide a foundation for understanding the population variation within invasive species among the Alternanthera 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.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae)

BACKGROUND

Corydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes.

RESULTS

We obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit β (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera.

CONCLUSIONS

The C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

Plastid genomes provide insights into the phylogeny and chloroplast evolution of the paper daisy tribe Gnaphalieae (Asteraceae)

Chloroplast genomes, as an essential source of phylogenetic information, are increasingly utilized in the evolutionary study of angiosperms. Gnaphalieae is a medium-sized tribe of the sunflower family of Asteraceae, with about 2,100 species in 178 genera distributed in temperate habitats worldwide. There has been considerable progress in our understanding of their phylogenetic evolution using both nuclear and chloroplast sequences, but no focus on chloroplast genomic data. In this study, we performed sequencing, assembly, and annotation of 16 representative chloroplast genomes from all the major lineages of Gnaphalieae. Our results showed that the plastomes exhibited a typical circular tetrad structure with similar genomic structure gene content. But there were differences in genome size, SSRs, and codon usage within the tribe. Phylogenetic analysis revealed Relhania clade is the earliest diverged lineages with the Lasiopogon clade and the Gnaphalium s.s. clade diverged subsequently. The core group includes FLAG clade sister to the HAP and Australasian group. Compared with the outgroup species, chloroplast genome size of the FLAG clade is much reduced whereas those of Australasian, HAP, Gnaphalium s.s., Lasiopogon and Relhania clades are relatively expanded. Insertions and deletions in the intergenic regions associated with repetitive sequence variations are supposed to be the main factor leading to length variations in the chloroplast genomes of Gnaphalieae. The comparative analyses of chloroplast genomes would provide useful implications into understanding the taxonomic and evolutionary history of Gnaphalieae.

Complete chloroplast genome sequence of Pinus tabuliformis var. henryi (Mast.) C.T.Kuan 1983 (Pinaceae)

Pinus tabuliformis var. henryi (Mast.) C.T.Kuan 1983 is an endemic and rare subtropical pine, mainly distributed in central China. In this study, we sequenced the complete chloroplast (cp) genome of P. tabuliformis var. henryi and reported it for the first time. The cp genome was 119,634 base pairs (bp) in total length, including two inverted repeats (IRs, 495 bp), separated by a large single-copy region (LSC, 65,600 bp) and a small single-copy region (SSC, 53,044 bp). There are 114 different genes in the cp genome of P. tabuliformis var. henryi, including 74 protein-coding genes, 36 transfer RNA genes, and four ribosomal RNA genes. The overall GC content of the cp genome was 38.5%. Our phylogenetic analysis of P. tabuliformis var. henryi demonstrated that it was closely related to P. tabuliformis and could be used to identify and analyze its genetic diversity, which was expected to provide new data for taxonomic and phylogenetic studies of Pinus.

Complete Chloroplast Genomes and the Phylogenetic Analysis of Three Native Species of Paeoniaceae from the Sino-Himalayan Flora Subkingdom

Paeonia delavayi var. lutea, Paeonia delavayi var. angustiloba, and Paeonia ludlowii are Chinese endemics that belong to the Paeoniaceae family and have vital medicinal and ornamental value. It is often difficult to classify Paeoniaceae plants based on their morphological characteristics, and the limited genomic information has strongly hindered molecular evolution and phylogenetic studies of Paeoniaceae. In this study, we sequenced, assembled, and annotated the chloroplast genomes of P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii. The chloroplast genomes of these strains were comparatively analyzed, and their phylogenetic relationships and divergence times were inferred. These three chloroplast genomes exhibited a typical quadripartite structure and were 152,687-152,759 bp in length. Each genome contains 126-132 genes, including 81-87 protein-coding genes, 37 transfer RNAs, and 8 ribosomal RNAs. In addition, the genomes had 61-64 SSRs, with mononucleotide repeats being the most abundant. The codon bias patterns of the three species tend to use codons ending in A/U. Six regions of high variability were identified (psbK-psbL, trnG-UCC, petN-psbM, psbC, rps8-rpl14, and ycf1) that can be used as DNA molecular markers for phylogenetic and taxonomic analysis. The Ka/Ks ratio indicates positive selection for the rps18 gene associated with self-replication. The phylogenetic analysis of 99 chloroplast genomes from Saxifragales clarified the phylogenetic relationships of Paeoniaceae and revealed that P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii are monophyletic groups and sisters to P. delavayi. Divergence time estimation revealed two evolutionary divergences of Paeoniaceae species in the early Oligocene and Miocene. Afterward, they underwent rapid adaptive radiation from the Pliocene to the early Pleistocene when P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii formed. The results of this study enrich the chloroplast genomic information of Paeoniaceae and reveal new insights into the phylogeny of Paeoniaceae.

The complete chloroplast genome sequence of Zygophyllum brachypterum (Zygophyllaceae) reveals its distinctive characteristics and evolutionary implication

Zygophyllum brachypterum Karelin & Kirilov belongs to Zygophyllaceae and is mainly distributed in the desert regions of Central Asia, Mongolia, and Northwest China. The species is valuable in exploring the adaptations of Zygophyllaceae plants to salt stress in ecological environments. In this study, we report the complete chloroplast (cp) genome of Z. brachypterum. The entire cp genome was 104590 bp in length, with a large single-copy region (LSC, 79170 bp), a small single-copy region (SSC, 16778 bp), and two inverted repeats (IRa/IRb) of 4321 bp each. A total of 106 genes were detected, among which seven were located in the IRs, and 65, 30, and 4 were protein-coding, tRNA, and rRNA genes, respectively. Notably, eleven genes encoding the subunits of NAD(P)H dehydrogenase complex (NDH) were absent. Phylogenetic analysis indicated that Z. brachypterum belonged to Zygophylloideae (Zygophyllaceae). Furthermore, it was closely related to Z. fabago and Z. kansuense.

Loss of plastid ndh genes in an autotrophic desert plant

Plant plastid genomes are highly conserved with most flowering plants having the same complement of essential plastid genes. Here, we report the loss of five of the eleven NADH dehydrogenase subunit genes (ndh) in the plastid of a desert plant jojoba (Simmondsia chinensis). The plastid genome of jojoba was 156,496 bp with one large single copy region (LSC), a very small single copy region (SSC) and two expanded inverted repeats (IRA + IRB). The NADH dehydrogenase (NDH) complex is comprised of several protein subunits, encoded by the ndh genes of the plastome and the nucleus. The ndh genes are critical to the proper functioning of the photosynthetic electron transport chain and protection of plants from oxidative stress. Most plants are known to contain all eleven ndh genes. Plants with missing or defective ndh genes are often heterotrophs either due to their complete or holo- or myco- parasitic nature. Plants with a defective NDH complex, caused by the deletion/pseudogenisation of some or all the ndh genes, survive in milder climates suggesting the likely extinction of plant lineages lacking these genes under harsh climates. Interestingly, some autotrophic plants do exist without ndh gene/s and can cope with high or low light. This implies that these plants are protected from oxidative stress by mechanisms excluding ndh genes. Jojoba has evolved mechanisms to cope with a non-functioning NDH complex and survives in extreme desert conditions with abundant sunlight and limited water.

Analysis of the RNA Editing Sites and Orthologous Gene Function of Transcriptome and Chloroplast Genomes in the Evolution of Five Deutzia Species

In this study, the chloroplast genomes and transcriptomes of five Deutzia genus species were sequenced, characterized, combined, and analyzed. A phylogenetic tree was constructed, including 32 other chloroplast genome sequences of Hydrangeoideae species. The results showed that the five Deutzia chloroplast genomes were typical circular genomes 156,860-157,025 bp in length, with 37.58-37.6% GC content. Repeat analysis showed that the Deutzia species had 41-45 scattered repeats and 199-201 simple sequence repeats. Comparative genomic and pi analyses indicated that the genomes are conservative and that the gene structures are stable. According to the phylogenetic tree, Deutzia species appear to be closely related to Kirengeshoma palmata and Philadelphus. By combining chloroplast genomic and transcriptomic analyses, 29-31 RNA editing events and 163-194 orthologous genes were identified. The ndh, rpo, rps, and atp genes had the most editing sites, and all RNA editing events were of the C-to-U type. Most of the orthologous genes were annotated to the chloroplast, mitochondria, and nucleus, with functions including energy production and conversion, translation, and protein transport. Genes related to the biosynthesis of monoterpenoids and flavonoids were also identified from the transcriptome of Deutzia spp. Our results will contribute to further studies of the genomic information and potential uses of the Deutzia spp.

The complete chloroplast genome of Ulmus mianzhuensis with insights into structural variations, adaptive evolution, and phylogenetic relationships of Ulmus (Ulmaceae)

BACKGROUND

Ulmus mianzhuensis is an endemic tree species in China with high ornamental and economic value. Currently, little is known regarding its genomic architecture, phylogenetic position, or adaptive evolution. Here, we sequenced the complete chloroplast genome (cp genome) of U. mianzhuensis and further compared the variations in gene organization and structure within Ulmus species to define their genomic evolution, then reconstructed the phylogenomic relationship of 31 related Ulmus species to explore the systematic position of U. mianzhuensis and the utility of cp genome for resolving phylogenetics among Ulmus species.

RESULTS

Our results revealed that all the Ulmus species exhibited a typical quadripartite structure, with a large single copy (LSC) region of 87,170 - 88,408 bp, a small single copy (SSC) region of 18,650 - 19,038 bp and an inverted repeat (IR) region of 26,288 - 26,546 bp. Within Ulmus species, gene structure and content of cp genomes were highly conserved, although slight variations were found in the boundary of SC/IR regions. Moreover, genome-wide sliding window analysis uncovered the variability of ndhC-trnV-UAC, ndhF-rpl32, and psbI-trnS-GCU were higher among 31 Ulmus that may be useful for the population genetics and potential DNA barcodes. Two genes (rps15 and atpF) were further detected under a positive selection of Ulmus species. Comparative phylogenetic analysis based on the cp genome and protein-coding genes revealed consistent topology that U. mianzhuensis is a sister group to U. parvifolia (sect. Microptelea) with a relatively low-level nucleotide variation of the cp genome. Additionally, our analyses also found that the traditional taxonomic system of five sections in Ulmus is not supported by the current phylogenomic topology with a nested evolutionary relationship between sections.

CONCLUSIONS

Features of the cp genome length, GC content, organization, and gene order were highly conserved within Ulmus. Furthermore, molecular evidence from the low variation of the cp genome suggested that U. mianzhuensis should be merged into U. parvifolia and regarded as a subspecies of U. parvifolia. Overall, we demonstrated that the cp genome provides valuable information for understanding the genetic variation and phylogenetic relationship in Ulmus.

Insights into phylogenetic relationships in Pinus inferred from a comparative analysis of complete chloroplast genomes

BACKGROUND

Pinus is the largest genus of Pinaceae and the most primitive group of modern genera. Pines have become the focus of many molecular evolution studies because of their wide use and ecological significance. However, due to the lack of complete chloroplast genome data, the evolutionary relationship and classification of pines are still controversial. With the development of new generation sequencing technology, sequence data of pines are becoming abundant. Here, we systematically analyzed and summarized the chloroplast genomes of 33 published pine species.

RESULTS

Generally, pines chloroplast genome structure showed strong conservation and high similarity. The chloroplast genome length ranged from 114,082 to 121,530 bp with similar positions and arrangements of all genes, while the GC content ranged from 38.45 to 39.00%. Reverse repeats showed a shrinking evolutionary trend, with IRa/IRb length ranging from 267 to 495 bp. A total of 3,205 microsatellite sequences and 5,436 repeats were detected in the studied species chloroplasts. Additionally, two hypervariable regions were assessed, providing potential molecular markers for future phylogenetic studies and population genetics. Through the phylogenetic analysis of complete chloroplast genomes, we offered novel opinions on the genus traditional evolutionary theory and classification.

CONCLUSION

We compared and analyzed the chloroplast genomes of 33 pine species, verified the traditional evolutionary theory and classification, and reclassified some controversial species classification. This study is helpful in analyzing the evolution, genetic structure, and the development of chloroplast DNA markers in Pinus.

Complete Chloroplast Genome of Corethrodendron fruticosum (Papilionoideae: Fabaceae): Comparative and Phylogenetic Analysis

Corethrodendron fruticosum is an endemic forage grasses in China with high ecological value. In this study, the complete chloroplast genome of C. fruticosum was sequenced using Illumina paired-end sequencing. The C. fruticosum chloroplast genome was 123,100 bp and comprised 105 genes, including 74 protein-coding genes, 4 rRNA-coding genes, and 27 tRNA-coding genes. The genome had a GC content of 34.53%, with 50 repetitive sequences and 63 simple repeat repetitive sequences that did not contain reverse repeats. The simple repeats included 45 single-nucleotide repeats, which accounted for the highest proportion and primarily comprised A/T repeats. A comparative analysis of C. fruticosum, C. multijugum, and four Hedysarum species revealed that the six genomes were highly conserved, with differentials primarily located in the conserved non-coding regions. Moreover, the accD and clpP genes in the coding regions exhibited high nucleotide variability. Accordingly, these genes may serve as molecular markers for the classification and phylogenetic analysis of Corethrodendron species. Phylogenetic analysis further revealed that C. fruticosum and C. multijugum appeared in different clades than the four Hedysarum species. The newly sequenced chloroplast genome provides further insights into the phylogenetic position of C. fruticosum, which is useful for the classification and identification of Corethrodendron.

Flavodiiron-mediated O2 photoreduction at photosystem I acceptor-side provides photoprotection to conifer thylakoids in early spring

Green organisms evolve oxygen (O₂) via photosynthesis and consume it by respiration. Generally, net O₂ consumption only becomes dominant when photosynthesis is suppressed at night. Here, we show that green thylakoid membranes of Scots pine (Pinus sylvestris L) and Norway spruce (Picea abies) needles display strong O₂ consumption even in the presence of light when extremely low temperatures coincide with high solar irradiation during early spring (ES). By employing different electron transport chain inhibitors, we show that this unusual light-induced O₂ consumption occurs around photosystem (PS) I and correlates with higher abundance of flavodiiron (Flv) A protein in ES thylakoids. With P700 absorption changes, we demonstrate that electron scavenging from the acceptor-side of PSI via O₂ photoreduction is a major alternative pathway in ES. This photoprotection mechanism in vascular plants indicates that conifers have developed an adaptative evolution trajectory for growing in harsh environments.

The First Complete Chloroplast Genome of Cordia monoica: Structure and Comparative Analysis

Cordia monoica is a member of the Boraginaceae family. This plant is widely distributed in tropical regions and has a great deal of medical value as well as economic importance. In the current study, the complete chloroplast (cp) genome of C. monoica was sequenced, assembled, annotated, and reported. This circular chloroplast genome had a size of 148,711 bp, with a quadripartite structure alternating between a pair of repeated inverted regions (26,897-26,901 bp) and a single copy region (77,893 bp). Among the 134 genes encoded by the cp genome, there were 89 protein-coding genes, 37 transfer RNA (tRNA) genes, and 8 ribosomal RNA (rRNA) genes. A total of 1387 tandem repeats were detected, with the hexanucleotides class making up 28 percent of the repeats. Cordia monoica has 26,303 codons in its protein-coding regions, and leucine amino acid was the most frequently encoded amino acid in contrast to cysteine. In addition, 12 of the 89 protein-coding genes were found to be under positive selection. The phyloplastomic taxonomical clustering of the Boraginaceae species provides further evidence that chloroplast genome data are reliable not only at family level but also in deciphering the phylogeny at genus level (e.g., Cordia).

The phylogenomics and evolutionary dynamics of the organellar genomes in carnivorous Utricularia and Genlisea species (Lentibulariaceae)

Utricularia and Genlisea are highly specialized carnivorous plants whose phylogenetic history has been poorly explored using phylogenomic methods. Additional sampling and genomic data are needed to advance our phylogenetic and taxonomic knowledge of this group of plants. Within a comparative framework, we present a characterization of plastome (PT) and mitochondrial (MT) genes of 26 Utricularia and six Genlisea species, with representatives of all subgenera and growth habits. All PT genomes maintain similar gene content, showing minor variation across the genes located between the PT junctions. One exception is a major variation related to different patterns in the presence and absence of ndh genes in the small single copy region, which appears to follow the phylogenetic history of the species rather than their lifestyle. All MT genomes exhibit similar gene content, with most differences related to a lineage-specific pseudogenes. We find evidence for episodic positive diversifying selection in PT and for most of the Utricularia MT genes that may be related to the current hypothesis that bladderworts' nuclear DNA is under constant ROS oxidative DNA damage and unusual DNA repair mechanisms, or even low fidelity polymerase that bypass lesions which could also be affecting the organellar genomes. Finally, both PT and MT phylogenetic trees were well resolved and highly supported, providing a congruent phylogenomic hypothesis for Utricularia and Genlisea clade given the study sampling.

A phylogenomic study of Iridaceae Juss. based on complete plastid genome sequences

The plastid genome has proven to be an effective tool for examining deep correlations in plant phylogenetics, owing to its highly conserved structure, uniparental inheritance, and limited variation in evolutionary rates. Iridaceae, comprising more than 2,000 species, includes numerous economically significant taxa that are frequently utilized in food industries and medicines and for ornamental and horticulture purposes. Molecular studies on chloroplast DNA have confirmed the position of this family in the order Asparagales with non-asparagoids. The current subfamilial classification of Iridaceae recognizes seven subfamilies-Isophysioideae, Nivenioideae, Iridoideae, Crocoideae, Geosiridaceae, Aristeoideae, and Patersonioideae-which are supported by limited plastid DNA regions. To date, no comparative phylogenomic studies have been conducted on the family Iridaceae. We assembled and annotated (de novo) the plastid genomes of 24 taxa together with seven published species representing all the seven subfamilies of Iridaceae and performed comparative genomics using the Illumina MiSeq platform. The plastomes of the autotrophic Iridaceae represent 79 protein-coding, 30 tRNA, and four rRNA genes, with lengths ranging from 150,062 to 164,622 bp. The phylogenetic analysis of the plastome sequences based on maximum parsimony, maximum likelihood, and Bayesian inference analyses suggested that Watsonia and Gladiolus were closely related, supported by strong support values, which differed considerably from recent phylogenetic studies. In addition, we identified genomic events, such as sequence inversions, deletions, mutations, and pseudogenization, in some species. Furthermore, the largest nucleotide variability was found in the seven plastome regions, which can be used in future phylogenetic studies. Notably, three subfamilies-Crocoideae, Nivenioideae, and Aristeoideae-shared a common ycf2 gene locus deletion. Our study is a preliminary report of a comparative study of the complete plastid genomes of 7/7 subfamilies and 9/10 tribes, elucidating the structural characteristics and shedding light on plastome evolution and phylogenetic relationships within Iridaceae. Additionally, further research is required to update the relative position of Watsonia within the tribal classification of the subfamily Crocoideae.

The chloroplast genome of black pepper (Piper nigrum L.) and its comparative analysis with related Piper species

Piper nigrum, also known as black pepper, is an economically and ecologically important crop of the genus Piper. It has been titled as the king of spices due to its wide consumption throughout the world. In the present investigation, the chloroplast genome of P. nigrum has been assembled from a whole genome sequence by integrating the short and long reads generated through Illumina and PacBio platforms, respectively. The chloroplast genome was observed to be 161,522 bp in size, having a quadripartite structure with a large single copy (LSC) region of 89,153 bp and a small single copy (SSC) region of 18,255 bp separated by a copy of inverted repeats (IRs), each 27,057 bp in length. Taking into consideration all the duplicated genes, a total of 131 genes were observed, which included 81 protein-coding genes, 37 tRNAs, 4 rRNAs, and 1 pseudogene. Individually, the LSC region consisted of 83 genes, the SSC region had 13 genes, and 18 genes were present in each IR region. Additionally, 216 SSRs were detected and 11 of these were validated through amplification in 12 species of Piper. The features of the chloroplast genome have been compared with those of the genus Piper. Our results provide useful insights into evolutionary and molecular studies of black pepper which will contribute to its further genetic improvement and breeding.

Phylogenetic relationship and comparative analysis of the main Bupleuri Radix species in China

Background

Bupleuri Radix (Chaihu) is a famous traditional Chinese medicine derived from Bupleurum, Apiaceae. The origin of cultivated Chaihu germplasm in China is unclear, which has led to unstable Chaihu quality. In this study, we reconstructed the phylogeny of the main Chaihu germplasm species in China and identified potential molecular markers to authenticate its origin.

Methods

Three Bupleurum species (eight individuals), B. bicaule, B. chinense, and B. scorzonerifolium, were selected for genome skimming. Published genomes from B. falcatum and B. marginatum var. stenophyllum were used for comparative analysis.

Results

Sequences of the complete plastid genomes were conserved with 113 identical genes ranging from 155,540 to 155,866 bp in length. Phylogenetic reconstruction based on complete plastid genomes resolved intrageneric relationships of the five Bupleurum species with high support. Conflicts between the plastid and nuclear phylogenies were observed, which were mainly ascribed to introgressive hybridization. Comparative analysis showed that noncoding regions of the plastomes had most of the variable sequences. Eight regions (atpF-atpH, petN-psbM, rps16-psbK, petA-psbJ, ndhC-trnV/UAC and ycf1) had high divergence values in Bupleurum species and could be promising DNA barcodes for Chaihu authentication. A total of seven polymorphic cpSSRs and 438 polymorphic nSSRs were detected across the five Chaihu germplasms. Three photosynthesis-related genes were under positive selection, of which accD reflected the adaptation fingerprint of B. chinense to different ecological habitats. Our study provides valuable genetic information for phylogenetic investigation, germplasm authentication, and molecular breeding of Chaihu species.

Comparison of the chloroplast genomes and phylogenomic analysis of Elaeocarpaceae

Background

Elaeocarpaceae is a vital family in tropical and subtropical forests. Compared with the important position of Elaeocarpaceae species in forest ecosystem and the concern of medicinal value, the most research on Elaeocarpaceae are classification and taxonomy. Molecular systematics has corrected the morphological misjudgment, and it belongs to Oxalidales. Phylogenetic and divergence time estimates of Elaeocarpaceae is mostly constructed by using chloroplast gene fragments. At present, although there are reports on the chloroplast structure of Elaeocarpaceae, a comprehensive analysis of the chloroplast structure of Elaeocarpaceae is lacking.

Methods

To understand the variation in chloroplast sequence size and structure in Elaeocarpaceae, the chloroplast genomes of nine species were sequenced using the Illumina HiSeq 2500 platform and further assembled and annotated with Elaeocarpus japonicus and Sloanea sinensis (family Elaeocarpaceae) as references. A phylogenomic tree was constructed based on the complete chloroplast genomes of the 11 species representing five genera of Elaeocarpaceae. Chloroplast genome characteristics were examined by using Circoletto and IRscope software.

Results

The results revealed the following: (a) The 11 sequenced chloroplast genomes ranged in size from 157,546 to 159,400 bp. (b) The chloroplast genomes of Elaeocarpus, Sloanea, Crinodendron and Vallea lacked the rpl32 gene in the small single-copy (SSC) region. The large single-copy (LSC) region of the chloroplast genomes lacked the ndhK gene in Elaeocarpus, Vallea stipularis, and Aristotelia fruticosa. The LSC region of the chloroplast genomes lacked the infA gene in genus Elaeocarpus and Crinodendron patagua. (c) Through inverted repeat (IR) expansion and contraction analysis, a significant difference was found between the LSC/IRB and IRA/LSC boundaries among these species. Rps3 was detected in the neighboring regions of the LSC and IRb regions in Elaeocarpus. (d) Phylogenomic analysis revealed that the genus Elaeocarpus is closely related to Crinodendron patagua on an independent branch and Aristotelia fruticosa is closely related to Vallea stipularis, forming a clade with the genus Sloanea. Structural comparisons showed that Elaeocarpaceae diverged at 60 Mya, the genus Elaeocarpus diverged 53 Mya and that the genus Sloanea diverged 0.44 Mya. These results provide new insight into the evolution of the Elaeocarpaceae.

Dynamic evolution of the plastome in the Elm family (Ulmaceae)

This study compared the plastomes of Ulmaceae allowing analyses of the dynamic evolution, including genome structure, codon usage bias, repeat sequences, molecular mutation rates, and phylogenetic inferences. Ulmaceae is a small family in the order Rosales. This family consists of seven genera, including Ulmus, Zelkova, Planera, Hemiptelea, Phyllostylon, Ampelocera, and Holoptelea. Ulmaceae is an interesting lineage from plant biogeographic, systematic, evolutionary, and paleobotanic perspectives. It is also a good model to investigate the evolution of the plastomes in woody plants. In this study, we sequenced and assembled the complete plastomes of the six Ulmaceae genera to compare genomic structures and reveal the molecular evolutionary patterns. The size of the quadripartite plastomes ranged from 158,290 bp to 161,886 bp. The genomes contained 131 genes, including 87 coding genes, 36 tRNA, and 8 rRNA. The gene number, gene content, and genomic structure were highly consistent among the Ulmaceae genera. Nine variable regions including ndhA intron, ndhF-rpl32, ycf1, psbK-trnS, rps16-trnQ, trnT-trnL, trnT-psbD, trnS-trnG, and rpl32-trnL, were identified in Ulmaceae plastomes according to the nucleotide diversity values. Condon usage was biased among the genes and showed consistent trends in the seven genera. Molecular evolution analyses revealed that most of the genes and all gene groups were under widespread purifying selection. Twelve genes (ccsA, matK, psbH, psbK, rbcL, rpl22, rpl32, rpoA, rps12, rps15, rps16, and ycf2) were under positive selection. Phylogenetic analyses supported that Ulmaceae should be divided into two main clades, such as the temperate clade, including Ulmus, Zelkova, Planera, and Hemiptelea and the tropical clade, including Phyllostylon, Ampelocera and Holoptelea. This study reports the structure and evolutionary characteristics of the Elm family. These new genomic data will benefit assessments of genomic evolution and provide information to elucidate the phylogenetic relationships among Ulmaceae species.

Assembly and Annotation of Red Spruce (Picea rubens) Chloroplast Genome, Identification of Simple Sequence Repeats, and Phylogenetic Analysis in Picea

We have sequenced the chloroplast genome of red spruce (Picea rubens) for the first time using the single-end, short-reads (44 bp) Illumina sequences, assembled and functionally annotated it, and identified simple sequence repeats (SSRs). The contigs were assembled using SOAPdenovo2 following the retrieval of chloroplast genome sequences using the black spruce (Picea mariana) chloroplast genome as the reference. The assembled genome length was 122,115 bp (gaps included). Comparatively, the P. rubens chloroplast genome reported here may be considered a near-complete draft. Global genome alignment and phylogenetic analysis based on the whole chloroplast genome sequences of Picea rubens and 10 other Picea species revealed high sequence synteny and conservation among 11 Picea species and phylogenetic relationships consistent with their known classical interrelationships and published molecular phylogeny. The P. rubens chloroplast genome sequence showed the highest similarity with that of P. mariana and the lowest with that of P. sitchensis. We have annotated 107 genes including 69 protein-coding genes, 28 tRNAs, 4 rRNAs, few pseudogenes, identified 42 SSRs, and successfully designed primers for 26 SSRs. Mononucleotide A/T repeats were the most common followed by dinucleotide AT repeats. A similar pattern of microsatellite repeats occurrence was found in the chloroplast genomes of 11 Picea species.

Phylogenomic analysis and development of molecular markers for the determination of twelve plum cultivars (Prunus, Rosaceae)

Background

Plums are one of the most important economic crops of the Rosaceae family and are produced all over the world. China has many local varieties, but the genomic information is limited for genetic studies. Here, we first sequenced, assembled, and analyzed the plastomes of twelve plum cultivars and developed molecular markers to distinguish them.

Results

The twelve plastomes of plum cultivars have a circular structure of 157,863–157,952 bp containing a large single-copy region (LSC) of 86,109–86,287 bp, a small copy region (SSC) of 18,927–19,031 bp, and two inverted repeats (IR) of 26,353–26,387 bp each. The plastomes of plum cultivars encode 131 genes, including 86 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. We detected 50, 54, 54, 53, 53, 50, 54, 54, 54, 49, 50, 54 SSRs in the twelve analyzed varieties, respectively. For repeat sequences, we identified 553 tandem repeats, 204 direct repeats, and 270 palindromic repeats. We also analyzed the expansion/contraction of IR regions. The genes rpl22, rps19, rpl2, ycf1, ndhF, and the trnH span on or near the boundary of IR and single-copy regions. Phylogenetic analysis showed that the twelve cultivars were clustered with the P. salicina and P. domestica. We developed eight markers LZ01 to LZ08 based on whole plastomes and nuclear genes and validated them successfully with six repetitions.

Conclusions

The results obtained here could fill in the blanks of the plastomes of these twelve plum cultivars and provide a wider perspective based on the basis of the plastomes of Prunus to the molecular identification and phylogenetic construction accurately. The analysis from this study provides an important and valuable resource for studying the genetic basis for agronomic and adaptive differentiation of the Prunus species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08965-z.

The photosynthesis apparatus of European mistletoe (Viscum album)

European mistletoe (Viscum album) is known for its special mode of cellular respiration. It lacks the mitochondrial NADH dehydrogenase complex (Complex I of the respiratory chain) and has restricted capacities to generate mitochondrial adenosine triphosphate (ATP). Here, we present an investigation of the V. album energy metabolism taking place in chloroplasts. Thylakoids were purified from young V. album leaves, and membrane-bound protein complexes were characterized by Blue native polyacrylamide gel electrophoresis as well as by the complexome profiling approach. Proteins were systematically identified by label-free quantitative shotgun proteomics. We identified >1,800 distinct proteins (accessible at https://complexomemap.de/va_leaves), including nearly 100 proteins forming part of the protein complexes involved in the light-dependent part of photosynthesis. The photosynthesis apparatus of V. album has distinct features: (1) comparatively low amounts of Photosystem I; (2) absence of the NDH complex (the chloroplast pendant of mitochondrial Complex I involved in cyclic electron transport (CET) around Photosystem I); (3) reduced levels of the proton gradient regulation 5 (PGR5) and proton gradient regulation 5-like 1 (PGRL1) proteins, which offer an alternative route for CET around Photosystem I; (4) comparable amounts of Photosystem II and the chloroplast ATP synthase complex to other seed plants. Our data suggest a restricted capacity for chloroplast ATP biosynthesis by the photophosphorylation process. This is in addition to the limited ATP supply by the mitochondria. We propose a view on mistletoe's mode of life, according to which its metabolism relies to a greater extent on energy-rich compounds provided by the host trees.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Comparative chloroplast genome analysis of Ficus (Moraceae): Insight into adaptive evolution and mutational hotspot regions

As the largest genus in Moraceae, Ficus is widely distributed across tropical and subtropical regions and exhibits a high degree of adaptability to different environments. At present, however, the phylogenetic relationships of this genus are not well resolved, and chloroplast evolution in Ficus remains poorly understood. Here, we sequenced, assembled, and annotated the chloroplast genomes of 10 species of Ficus, downloaded and assembled 13 additional species based on next-generation sequencing data, and compared them to 46 previously published chloroplast genomes. We found a highly conserved genomic structure across the genus, with plastid genome sizes ranging from 159,929 bp (Ficus langkokensis) to 160,657 bp (Ficus religiosa). Most chloroplasts encoded 113 unique genes, including a set of 78 protein-coding genes, 30 transfer RNA (tRNA) genes, four ribosomal RNA (rRNA) genes, and one pseudogene (infA). The number of simple sequence repeats (SSRs) ranged from 67 (Ficus sagittata) to 89 (Ficus microdictya) and generally increased linearly with plastid size. Among the plastomes, comparative analysis revealed eight intergenic spacers that were hotspot regions for divergence. Additionally, the clpP, rbcL, and ccsA genes showed evidence of positive selection. Phylogenetic analysis indicated that none of the six traditionally recognized subgenera of Ficus were monophyletic. Divergence time analysis based on the complete chloroplast genome sequences showed that Ficus species diverged rapidly during the early to middle Miocene. This research provides basic resources for further evolutionary studies of Ficus.

Complete Chloroplast Genomes of 14 Subspecies of D. glomerata: Phylogenetic and Comparative Genomic Analyses

Orchardgrass (Dactylis glomerata L.) is a species in the Gramineae family that is highly important economically and valued for its role in ecology. However, the phylogeny and taxonomy of D. glomerata are still controversial based on current morphological and molecular evidence. The study of chloroplast (cp) genomes has developed into a powerful tool to develop molecular markers for related species and reveal the relationships between plant evolution and phylogenetics. In this study, we conducted comparative genomic analyses and phylogenetic inferences on 14 cp genomes of D. glomerata originating from the Mediterranean and Eurasia. The genome size ranged from 134,375 bp to 134,993 bp and exhibited synteny of gene organization and order. A total of 129–131 genes were identified, including 85–87 protein coding genes, 38 tRNA genes and 8 rRNA genes. The cp sequences were highly conserved, and key sequence variations were detected at the junctions of inverted repeats (IRs)/small single–copy (SSC) regions. Moreover, nine highly variable regions were identified among the subspecies based on a sequence divergence analysis. A total of 285 RNA editing sites were detected that were relevant to 52 genes, where rpoB exhibited the most abundant RNA editing sites. The phylogenetic analysis revealed that all Dactylis subspecies clustered into a monophyletic group and most branches provided a high support bootstrap. The main divergence time of D. glomerata was dated to the Miocene era, and this could have been due to changes in the climate. These findings will provide useful insights for further studies on phylogeny, the identification of subspecies and the development of hypotheses for the evolutionary history of the genus Dactylis and of the Gramineae family.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

Background

Calanthe(Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of theCalanthealliance (Calanthe,Cephalantheropsis, andPhaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed ade novoassembly, and did a comparative analysis of 8Calanthegroup species' plastomes: 6Calantheand 2Phaiusspecies. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

Results

The complete plastomes of theCalanthegroup species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost inC.delavayi(P. delavayi), includingndhC,ndhF, andndhKgenes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in theCalanthegroup, and confirmed the position ofPhaius delavayiin the genusCalantheas opposed to its previous placement inPhaius.

Conclusion

This study provides a report on the complete plastomes of 6Calantheand 2Phaiusspecies and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

Plastome structure of 8 Calanthe s.l. species (Orchidaceae): comparative genomics, phylogenetic analysis

BACKGROUND

Calanthe (Epidendroideae, Orchidaceae) is a pantropical genus distributed in Asia and Africa. Its species are of great importance in terms of economic, ornamental and medicinal values. However, due to limited and confusing delimitation characters, the taxonomy of the Calanthe alliance (Calanthe, Cephalantheropsis, and Phaius) has not been sufficiently resolved. Additionally, the limited genomic information has shown incongruences in its systematics and phylogeny. In this study, we used illumina platform sequencing, performed a de novo assembly, and did a comparative analysis of 8 Calanthe group species' plastomes: 6 Calanthe and 2 Phaius species. Phylogenetic analyses were used to reconstruct the relationships of the species as well as with other species of the family Orchidaceae.

RESULTS

The complete plastomes of the Calanthe group species have a quadripartite structure with varied sizes ranging between 150,105bp-158,714bp, including a large single-copy region (LSC; 83,364bp- 87,450bp), a small single-copy region (SSC; 16,297bp -18,586bp), and a pair of inverted repeat regions (IRs; 25,222bp - 26,430bp). The overall GC content of these plastomes ranged between 36.6-36.9%. These plastomes encoded 131-134 differential genes, which included 85-88 protein-coding genes, 37-38 tRNA genes, and 8 rRNA genes. Comparative analysis showed no significant variations in terms of their sequences, gene content, gene order, sequence repeats and the GC content hence highly conserved. However, some genes were lost in C. delavayi (P. delavayi), including ndhC, ndhF, and ndhK genes. Compared to the coding regions, the non-coding regions had more sequence repeats hence important for species DNA barcoding. Phylogenetic analysis revealed a paraphyletic relationship in the Calanthe group, and confirmed the position of Phaius delavayi in the genus Calanthe as opposed to its previous placement in Phaius.

CONCLUSION

This study provides a report on the complete plastomes of 6 Calanthe and 2 Phaius species and elucidates the structural characteristics of the plastomes. It also highlights the power of plastome data to resolve phylogenetic relationships and clarifies taxonomic disputes among closely related species to improve our understanding of their systematics and evolution. Furthermore, it also provides valuable genetic resources and a basis for studying evolutionary relationships and population genetics among orchid species.

Comparative Analysis of the Complete Chloroplast Genomes in Allium Section Bromatorrhiza Species (Amaryllidaceae): Phylogenetic Relationship and Adaptive Evolution

With the development of molecular sequencing approaches, many taxonomic and phylogenetic problems of the genus Allium L. have been solved; however, the phylogenetic relationships of some subgenera or sections, such as section Bromatorrhiza, remain unresolved, which has greatly impeded our full understanding of the species relationships among the major clades of Allium. In this study, the complete chloroplast (cp) genomes of nine species in the Allium sect. Bromatorrhiza were determined using the Illumina paired-end sequencing, the NOVOPlasty de novo assembly strategy, and the PGA annotation method. The results showed that the cp genome exhibited high conservation and revealed a typical circular tetrad structure. Among the sect. Bromatorrhiza species, the gene content, SSRs, codon usage, and RNA editing site were similar. The genome structure and IR regions’ fluctuation were investigated while genes, CDSs, and non-coding regions were extracted for phylogeny reconstruction. Evolutionary rates (Ka/Ks values) were calculated, and positive selection analysis was further performed using the branch-site model. Five hypervariable regions were identified as candidate molecular markers for species authentication. A clear relationship among the sect. Bromatorrhiza species were detected based on concatenated genes and CDSs, respectively, which suggested that sect. Bromatorrhiza is monophyly. In addition, there were three genes with higher Ka/Ks values (rps2, ycf1, and ycf2), and four genes (rpoC2, atpF, atpI, and rpl14) were further revealed to own positive selected sites. These results provide new insights into the plastome component, phylogeny, and evolution of Allium species.

Genome-Wide Insights Into the Organelle Translocation of Photosynthetic NDH-1 Genes During Evolution

Translocation of chloroplast-located genes to mitochondria or nucleus is considered to be a safety strategy that impedes mutation of photosynthetic genes and maintains their household function during evolution. The organelle translocation strategy is also developed in photosynthetic NDH-1 (pNDH-1) genes but its understanding is still far from complete. Here, we found that the mutation rate of the conserved pNDH-1 genes was gradually reduced but their selection pressure was maintained at a high level during evolution from cyanobacteria to angiosperm. By contrast, oxygenic photosynthesis-specific (OPS) pNDH-1 genes had an opposite trend, explaining the reason why they were transferred from the reactive oxygen species (ROS)-enriched chloroplast to the ROS-barren nucleus. Further, genome-wide sequence analysis supported the possibility that all conserved pNDH-1 genes lost in chloroplast genomes of Chlorophyceae and Pinaceae were transferred to the ROS-less mitochondrial genome as deduced from their truncated pNDH-1 gene fragments. Collectively, we propose that the organelle translocation strategy of pNDH-1 genes during evolution is necessary to maintain the function of the pNDH-1 complex as an important antioxidant mechanism for efficient photosynthesis.

Genome skimming approach reveals the gene arrangements in the chloroplast genomes of the highly endangered Crocus L. species: Crocus istanbulensis (B.Mathew) Rukšāns

Crocus istanbulensis (B.Mathew) Rukšāns is one of the most endangered Crocus species in the world and has an extremely limited distribution range in Istanbul. Our recent field work indicates that no more than one hundred individuals remain in the wild. In the present study, we used genome skimming to determine the complete chloroplast (cp) genome sequences of six C. istanbulensis individuals collected from the locus classicus. The cp genome of C. istanbulensis has 151,199 base pairs (bp), with a large single-copy (LSC) (81,197 bp), small single copy (SSC) (17,524 bp) and two inverted repeat (IR) regions of 26,236 bp each. The cp genome contains 132 genes, of which 86 are protein-coding (PCGs), 8 are rRNA and 38 are tRNA genes. Most of the repeats are found in intergenic spacers of Crocus species. Mononucleotide repeats were most abundant, accounting for over 80% of total repeats. The cp genome contained four palindrome repeats and one forward repeat. Comparative analyses among other Iridaceae species identified one inversion in the terminal positions of LSC region and three different gene (psbA, rps3 and rpl22) arrangements in C. istanbulensis that were not reported previously. To measure selective pressure in the exons of chloroplast coding sequences, we performed a sequence analysis of plastome-encoded genes. A total of seven genes (accD, rpoC2, psbK, rps12, ccsA, clpP and ycf2) were detected under positive selection in the cp genome. Alignment-free sequence comparison showed an extremely low sequence diversity across naturally occurring C. istanbulensis specimens. All six sequenced individuals shared the same cp haplotype. In summary, this study will aid further research on the molecular evolution and development of ex situ conservation strategies of C. istanbulensis.

Plastid and mitochondrial phylogenomics reveal correlated substitution rate variation in Koenigia (Polygonoideae, Polygonaceae) and a reduced plastome for Koenigia delicatula including loss of all ndh genes

Koenigia, a genus proposed by Linnaeus, has a contentious taxonomic history. In particular, relationships among species and the circumscription of the genus relative to Aconogonon remain uncertain. To explore phylogenetic relationships of Koenigia with other members of tribe Persicarieae and to establish the timing of major evolutionary diversification events, genome skimming of organellar sequences was used to assemble plastomes and mitochondrial genes from 15 individuals representing 13 species. Most Persicarieae plastomes exhibit a conserved structure and content relative to other flowering plants. However, Koenigia delicatula has lost functional copies of all ndh genes and the intron from atpF. In addition, the rpl32 gene was relocated in the K. delicatula plastome, which likely occurred via overlapping inversions or differential expansion and contraction of the inverted repeat. The highly supported but conflicting relationships between plastome and mitochondrial trees and among gene trees complicates the circumscription of Koenigia, which could be caused by rapid diversification within a short period. Moreover, the plastome and mitochondrial trees revealed correlated variation in substitution rates among Persicarieae species, suggesting a shared underlying mechanism promoting evolutionary rate variation in both organellar genomes. The divergence of dwarf K. delicatula from other Koenigia species may be associated with the well-known Eocene Thermal Maximum 2 or Early Eocene Climatic Optimum event, while diversification of the core-Koenigia clade associates with the Mid-Miocene Climatic Optimum and the uplift of Qinghai-Tibetan Plateau and adjacent areas.

Phylogeography and Population Genetics of Rosa chinensis var. spontanea and R. lucidissima Complex, the Important Ancestor of Modern Roses

Rosa chinensis var. spontanea and R. lucidissima complex are the morphologically very similar key ancestors of modern roses with high importance in rose research and breeding. Although widely distributed in subtropical central and southwestern China, these two taxa are highly endangered. We sampled a total of 221 specimens and 330 DNA samples from 25 populations across the two taxa's whole range. Leaf morphological traits were compared. Two chloroplast DNA intergenic spacers (trnG-trnS, petL-psbE) and ITS were used for population genetics and phylogenetic study to delimit the boundary between the two taxa, assess the genetic variation, uncover the possible evolutionary mechanism responsible for the differentiation within the complex, and make the conservation recommendations. The complex exhibited high levels of genetic variation (h _TcpDNA = 0.768, h _TITS = 0.726) and high population differentiation even over small geographic distance. We suggest R. chinensis var. spontanea and R. lucidissma be treated as independent taxa, and the northern populations around and within the Sichuan Basin being R. chinensis var. spontanea, having broader leaflets and paler full-blooming flowers, while those in the middle and southern Yunnan-Guizhou Plateau and the adjacent regions being R. lucidissma, having narrower leaflets and darker full-blooming flowers. Transitional areas between the southeastern Sichuan Basin and northeastern Guizhou are the contact or the hybridization zone of the two taxa. Ancestral haplotypes of the complex (R. lucidissma) evolved at about 1.21-0.86 Mya in southeastern Yunnan-Guizhou Plateau and its adjacent regions and survived there during the Quaternary Oscillation. Ancestral haplotypes of R. chinensis var. spontanea deviated from R. lucidissma at about 0.022-0.031 Mya at the transitional areas (Daloushan and Wulingshan Mountains) between the northeastern edge of Yunnan-Guizhou Plaeteau and the southeastern border of Sichuan Basin, where they survived the LGM. The evolution of the complex included spatial isolation and inter-species hybridization. The complex's endangered status might be the result of over-exploitation for its ornamental and medical value, or due to reforestation of some originally open habitats. We provide specific recommendations for the two taxa's in situ and ex situ conservation.

Comparative and phylogenetic analyses of six Kenya Polystachya (Orchidaceae) species based on the complete chloroplast genome sequences

BACKGROUND

Polystachya Hook. is a large pantropical orchid genus (c. 240 species) distributed in Africa, southern Asia and the Americas, with the center of diversity in Africa. Previous studies on species of this genus have not obtained the complete chloroplast genomes, structures and variations. Additionally, the phylogenetic position of the genus in the Orchidaceae is still controversial and uncertain. Therefore, in this study, we sequenced the complete plastomes of six Kenya Polystachya species based on genome skimming, subjected them to comparative genomic analysis, and reconstructed the phylogenetic relationships with other Orchidaceae species.

RESULTS

The results exhibited that the chloroplast genomes had a typical quadripartite structure with conserved genome arrangement and moderate divergence. The plastomes of the six Polystachya species ranged from 145,484 bp to 149,274 bp in length and had an almost similar GC content of 36.9-37.0%. Gene annotation revealed 106-109 single-copy genes. In addition, 19 genes are duplicated in the inverted regions, and 16 genes each possessd one or more introns. Although no large structural variations were observed among the Polystachya plastomes, about 1 kb inversion was found in Polystachya modesta and all 11 ndh genes in the Polystachya plastomes were lost or pseudogenized. Comparative analysis of the overall sequence identity among six complete chloroplast genomes confirmed that for both coding and non-coding regions in Polystachya, SC regions exhibit higher sequence variation than IRs. Furthermore, there were various amplifications in the IR regions among the six Polystachya species. Most of the protein-coding genes of these species had a high degree of codon preference. We screened out SSRs and found seven relatively highly variable loci. Moreover, 13 genes were discovered with significant positive selection. Phylogenetic analysis showed that the six Polystachya species formed a monophyletic clade and were more closely related to the tribe Vandeae. Phylogenetic relationships of the family Orchidaceae inferred from the 85 chloroplast genome sequences were generally consistent with previous studies and robust.

CONCLUSIONS

Our study is the initial report of the complete chloroplast genomes of the six Polystachya species, elucidates the structural characteristics of the chloroplast genome of Polystachya, and filters out highly variable sequences that can contribute to the development of DNA markers for use in the study of genetic variability and evolutionary studies in Polystachya. In addition, the phylogenetic results strongly support that the genus of Polystachya is a part of the tribe Vandeae.

Plastome structure, evolution, and phylogeny of Selaginella

As one of the earliest land plant lineages, Selaginella is important for studying land plant evolution. It is the largest genus of lycophytes containing 700-800 species. Some unique characters of Selaginella plastomes have been reported, but based only on 20 species. There have been no plastome phylogenies of Selaginella based on a relatively large sampling, and no efforts have been made to resolve the phylogeny of the enigmatic Sinensis group whose relationships have been unclear based on small datasets. Here we investigated the structures of 59 plastomes representing 51 species covering all six subgenera and 18 sections of Selaginella except two sections and including the intriguing Sinensis group for the first time. Our major results include: (1) the plastome size of Selaginella ranges tremendously from 78,492 bp to 187,632 bp; (2) there are numerous gene losses in Selaginella comparing with other lycophytes, Isoëtaceae and Lycopodiaceae; (3) the gene contents and plastome structures in Selaginella vary lineage-specifically and all infrageneric taxa are well supported in the plastome phylogeny; (4) the ndh gene family tends to lose or pseudogenize in those species with DR structure and without other short or medium repeats; (5) the short and medium repeat regions in SC mediate many conformations causing diverse and complex plastome structures, and six new conformations are discovered; (6) forty-eight species sampled have high GC content (>50%) but three species in the Sinensis group have ∼30% GC content in plastomes, similar to most vascular plants; (7) the Sinensis group is monophyletic, includes at least two subgroups, and has the smallest plastomes in land plants except some parasitic plants, and their plastomes do not contain any tRNAs; (8) the younger lineages in Selaginella tend to have higher GC content, whereas the older lineages tend to have lower GC content; and (9) because of incomplete genomic data and abnormal structures or some unknown reasons, even the concatenated plastomes could not well resolve the phylogenetic relationships in Selaginella with confidence, highlighting the difficulty in resolving the phylogeny and evolution of this particularly important land plant lineage.

On the Edge of Dispensability, the Chloroplast ndh Genes

The polypeptides encoded by the chloroplast ndh genes and some nuclear genes form the thylakoid NADH dehydrogenase (Ndh) complex, homologous to the mitochondrial complex I. Except for Charophyceae (algae related to higher plants) and a few Prasinophyceae, all eukaryotic algae lack ndh genes. Among vascular plants, the ndh genes are absent in epiphytic and in some species scattered among different genera, families, and orders. The recent identification of many plants lacking plastid ndh genes allows comparison on phylogenetic trees and functional investigations of the ndh genes. The ndh genes protect Angiosperms under various terrestrial stresses, maintaining efficient photosynthesis. On the edge of dispensability, ndh genes provide a test for the natural selection of photosynthesis-related genes in evolution. Variable evolutionary environments place Angiosperms without ndh genes at risk of extinction and, probably, most extant ones may have lost ndh genes recently. Therefore, they are evolutionary endpoints in phylogenetic trees. The low number of sequenced plastid DNA and the long lifespan of some Gymnosperms lacking ndh genes challenge models about the role of ndh genes protecting against stress and promoting leaf senescence. Additional DNA sequencing in Gymnosperms and investigations into the molecular mechanisms of their response to stress will provide a unified model of the evolutionary and functional consequences of the lack of ndh genes.

Comparative plastomics of Amaryllidaceae: inverted repeat expansion and the degradation of the ndh genes in Strumaria truncata Jacq

Amaryllidaceae is a widespread and distinctive plant family contributing both food and ornamental plants. Here we present an initial survey of plastomes across the family and report on both structural rearrangements and gene losses. Most plastomes in the family are of similar gene arrangement and content however some taxa have shown gains in plastome length while in several taxa there is evidence of gene loss. Strumaria truncata shows a substantial loss of ndh family genes while three other taxa show loss of cemA, which has been reported only rarely. Our sparse sampling of the family has detected sufficient variation to suggest further sampling across the family could be a rich source of new information on plastome variation and evolution.

Heat-induced down-regulation of photosystem II protects photosystem I in honeysuckle (Lonicera japonica)

Honeysuckle (Lonicera japonica Thunb.) is a traditional medicinal plant in China which is often threatened by high temperature at midday during summer. Heat-induced effects on the photosynthetic apparatus in honeysuckle are associated with a depression of the photosystem II (PSII) photochemical efficiency. However, very limited information is available on regulation of photosynthetic electron flow in PSI photoprotection in heat-stressed honeysuckle. Simultaneous analyses of chlorophyll fluorescence and the change in absorbance of P700 showed that energy transformation and electron transfer activity in PSII decreased under heat stress, but the fraction of photo-oxidizable PSI (Pm) remained stable. With treatments at 38 and 42 °C, the photochemical electron transport in PSII was suppressed, whereas the cyclic electron flow (CEF) around PSI was induced. In addition, the levels of high energy state quenching (qE) and P700 oxidation increased significantly with increasing temperature. However, a decline of qE in antimycin A (AA)- or 3-(3,4-dichlorophenyl)-1,1-dimethylurea (DCMU)-treated leaves after heat treatment was observed, while P700 oxidation decreased only in the presence of AA. The results indicate that heat-induced inhibition of PSII and induction of CEF cooperatively protect PSI from ROS damages through moderate down-regulation of photosynthetic electron flow from PSII to PSI.

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

Comparative Chloroplast Genomes of Zosteraceae Species Provide Adaptive Evolution Insights Into Seagrass

Seagrasses are marine flowering plants found in tropical and sub-tropical areas that live in coastal regions between the sea and land. All seagrass species evolved from terrestrial monocotyledons, providing the opportunity to study plant adaptation to sea environments. Here, we sequenced the chloroplast genomes (cpGenomes) of three Zostera species, then analyzed and compared their cpGenome structures and sequence variations. We also performed a phylogenetic analysis using published seagrass chloroplasts and calculated the selection pressure of 17 species within seagrasses and nine terrestrial monocotyledons, as well as estimated the number of shared genes of eight seagrasses. The cpGenomes of Zosteraceae species ranged in size from 143,877 bp (Zostera marina) to 152,726 bp (Phyllospadix iwatensis), which were conserved and displayed similar structures and gene orders. Additionally, we found 17 variable hotspot regions as candidate DNA barcodes for Zosteraceae species, which will be helpful for studying the phylogenetic relationships and interspecies differences between seagrass species. Interestingly, nine genes had positive selection sites, including two ATP subunit genes (atpA and atpF), two ribosome subunit genes (rps4 and rpl20), two DNA-dependent RNA polymerase genes (rpoC1 and rpoC2), as well as accD, clpP, and ycf2. These gene regions may have played key roles in the seagrass adaptation to diverse environments. The Branch model analysis showed that seagrasses had a higher rate of evolution than terrestrial monocotyledons, suggesting that seagrasses experienced greater environmental pressure. Moreover, a branch-site model identified positively selected sites (PSSs) in ccsA, suggesting their involvement in the adaptation to sea environments. These findings are valuable for further investigations on Zosteraceae cpGenomes and will serve as an excellent resource for future studies on seagrass adaptation to sea environments.

Gene Loss, Pseudogenization in Plastomes of Genus Allium (Amaryllidaceae), and Putative Selection for Adaptation to Environmental Conditions

Amaryllidaceae is a large family with more than 1,600 species, belonging to 75 genera. The largest genus—Allium—is vast, comprising about a thousand species. Allium species (as well as other members of the Amaryllidaceae) are widespread and diversified, they are adapted to a wide range of habitats from shady forests to open habitats like meadows, steppes, and deserts. The genes present in chloroplast genomes (plastomes) play fundamental roles for the photosynthetic plants. Plastome traits could thus be associated with geophysical abiotic characteristics of habitats. Most chloroplast genes are highly conserved and are used as phylogenetic markers for many families of vascular plants. Nevertheless, some studies revealed signatures of positive selection in chloroplast genes of many plant families including Amaryllidaceae. We have sequenced plastomes of the following nine Allium (tribe Allieae of Allioideae) species: A. zebdanense, A. moly, A. victorialis, A. macleanii, A. nutans, A. obliquum, A. schoenoprasum, A. pskemense, A. platyspathum, A. fistulosum, A. semenovii, and Nothoscordum bivalve (tribe Leucocoryneae of Allioideae). We compared our data with previously published plastomes and provided our interpretation of Allium plastome genes’ annotations because we found some noteworthy inconsistencies with annotations previously reported. For Allium species we estimated the integral evolutionary rate, counted SNPs and indels per nucleotide position as well as compared pseudogenization events in species of three main phylogenetic lines of genus Allium to estimate whether they are potentially important for plant physiology or just follow the phylogenetic pattern. During examination of the 38 species of Allium and the 11 of other Amaryllidaceae species we found that rps16, rps2, infA, ccsA genes have lost their functionality multiple times in different species (regularly evolutionary events), while the pseudogenization of other genes was stochastic events. We found that the “normal” or “pseudo” state of rps16, rps2, infA, ccsA genes correlates well with the evolutionary line of genus the species belongs to. The positive selection in various NADH dehydrogenase (ndh) genes as well as in matK, accD, and some others were found. Taking into account known mechanisms of coping with excessive light by cyclic electron transport, we can hypothesize that adaptive evolution in genes, coding subunits of NADH-plastoquinone oxidoreductase could be driven by abiotic factors of alpine habitats, especially by intensive light and UV radiation.

Evolution of an assembly factor-based subunit contributed to a novel NDH-PSI supercomplex formation in chloroplasts

Chloroplast NADH dehydrogenase-like (NDH) complex is structurally related to mitochondrial Complex I and forms a supercomplex with two copies of Photosystem I (the NDH-PSI supercomplex) via linker proteins Lhca5 and Lhca6. The latter was acquired relatively recently in a common ancestor of angiosperms. Here we show that NDH-dependent Cyclic Electron Flow 5 (NDF5) is an NDH assembly factor in Arabidopsis. NDF5 initiates the assembly of NDH subunits (PnsB2 and PnsB3) and Lhca6, suggesting that they form a contact site with Lhca6. Our analysis of the NDF5 ortholog in Physcomitrella and angiosperm genomes reveals the subunit PnsB2 to be newly acquired via tandem gene duplication of NDF5 at some point in the evolution of angiosperms. Another Lhca6 contact subunit, PnsB3, has evolved from a protein unrelated to NDH. The structure of the largest photosynthetic electron transport chain complex has become more complicated by acquiring novel subunits and supercomplex formation with PSI.

The chloroplast NDH complex interacts with Photosystem I to form the NDH-PSI supercomplex. Here the authors show that Arabidopsis NDF5 shares a common ancestor with the NDH subunit PnsB2 and acts as an NDH assembly factor initiating the assembly of PnsB2 and the evolutionarily distinct PnsB3.

Plastomes from tribe Plantagineae (Plantaginaceae) reveal infrageneric structural synapormorphies and localized hypermutation for Plantago and functional loss of ndh genes from Littorella

Tribe Plantagineae (Plantaginaceae) comprises ~ 270 species in three currently recognized genera (Aragoa, Littorella, Plantago), of which Plantago is most speciose. Plantago plastomes exhibit several atypical features including large inversions, expansions of the inverted repeat, increased repetitiveness, intron losses, and gene-specific increases in substitution rate, but the prevalence of these plastid features among species and subgenera is unknown. To assess phylogenetic relationships and plastomic evolutionary dynamics among Plantagineae genera and Plantago subgenera, we generated 25 complete plastome sequences and compared them with existing plastome sequences from Plantaginaceae. Using whole plastome and partitioned alignments, our phylogenomic analyses provided strong support for relationships among major Plantagineae lineages. General plastid features-including size, GC content, intron content, and indels-provided additional support that reinforced major Plantagineae subdivisions. Plastomes from Plantago subgenera Plantago and Coronopus have synapomorphic expansions and inversions affecting the size and gene order of the inverted repeats, and particular genes near the inversion breakpoints exhibit accelerated nucleotide substitution rates, suggesting localized hypermutation associated with rearrangements. The Littorella plastome lacks functional copies of ndh genes, which may be related to an amphibious lifestyle and partial reliance on CAM photosynthesis.

Plastid Genomes of the Early Vascular Plant Genus Selaginella Have Unusual Direct Repeat Structures and Drastically Reduced Gene Numbers

The early vascular plants in the genus Selaginella, which is the sole genus of the Selaginellaceae family, have an important place in evolutionary history, along with ferns, as such plants are valuable resources for deciphering plant evolution. In this study, we sequenced and assembled the plastid genome (plastome) sequences of two Selaginella tamariscina individuals, as well as Selaginella stauntoniana and Selaginella involvens. Unlike the inverted repeat (IR) structures typically found in plant plastomes, Selaginella species had direct repeat (DR) structures, which were confirmed by Oxford Nanopore long-read sequence assembly. Comparative analyses of 19 lycophytes, including two Huperzia and one Isoetes species, revealed unique phylogenetic relationships between Selaginella species and related lycophytes, reflected by structural rearrangements involving two rounds of large inversions that resulted in dynamic changes between IR and DR blocks in the plastome sequence. Furthermore, we present other uncommon characteristics, including a small genome size, drastic reductions in gene and intron numbers, a high GC content, and extensive RNA editing. Although the 16 Selaginella species examined may not fully represent the genus, our findings suggest that Selaginella plastomes have undergone unique evolutionary events yielding genomic features unparalleled in other lycophytes, ferns, or seed plants.