NDH expression marks major transitions in plant evolution and reveals coordinate intracellular gene loss

Chloroplast NADH dehydrogenase-like complex-mediated cyclic electron flow is the main electron transport route in C4 bundle sheath cells

The superior productivity of C4 plants is achieved via a metabolic C4 cycle which acts as a CO2 pump across mesophyll and bundle sheath (BS) cells and requires an additional input of energy in the form of ATP. The importance of chloroplast NADH dehydrogenase-like complex (NDH) operating cyclic electron flow (CEF) around Photosystem I (PSI) for C4 photosynthesis has been shown in reverse genetics studies but the contribution of CEF and NDH to cell-level electron fluxes remained unknown. We have created gene-edited Setaria viridis with null ndhO alleles lacking functional NDH and developed methods for quantification of electron flow through NDH in BS and mesophyll cells. We show that CEF accounts for 84% of electrons reducing PSI in BS cells and most of those electrons are delivered through NDH while the contribution of the complex to electron transport in mesophyll cells is minimal. A decreased leaf CO2 assimilation rate and growth of plants lacking NDH cannot be rescued by supplying additional CO2. Our results indicate that NDH-mediated CEF is the primary electron transport route in BS chloroplasts highlighting the essential role of NDH in generating ATP required for CO2 fixation by the C3 cycle in BS cells.

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.

Distinctive plastome evolution in carnivorous angiosperms

BACKGROUND

Independent origins of carnivory in multiple angiosperm families are fabulous examples of convergent evolution using a diverse array of life forms and habitats. Previous studies have indicated that carnivorous plants have distinct evolutionary trajectories of plastid genome (plastome) compared to their non-carnivorous relatives, yet the extent and general characteristics remain elusive.

RESULTS

We compared plastomes from 9 out of 13 carnivorous families and their non-carnivorous relatives to assess carnivory-associated evolutionary patterns. We identified inversions in all sampled Droseraceae species and four species of Utricularia, Pinguicula, Darlingtonia and Triphyophyllum. A few carnivores showed distinct shifts in inverted repeat boundaries and the overall repeat contents. Many ndh genes, along with some other genes, were independently lost in several carnivorous lineages. We detected significant substitution rate variations in most sampled carnivorous lineages. A significant overall substitution rate acceleration characterizes the two largest carnivorous lineages of Droseraceae and Lentibulariaceae. We also observe moderate substitution rates acceleration in many genes of Cephalotus follicularis, Roridula gorgonias, and Drosophyllum lusitanicum. However, only a few genes exhibit significant relaxed selection.

CONCLUSION

Our results indicate that the carnivory of plants have different effects on plastome evolution across carnivorous lineages. The complex mechanism under carnivorous habitats may have resulted in distinctive plastome evolution with conserved plastome in the Brocchinia hechtioides to strongly reconfigured plastomes structures in Droseraceae. Organic carbon obtained from prey and the efficiency of utilizing prey-derived nutrients might constitute possible explanation.

More than a spiny morphology: plastome variation in the prickly pear cacti (Opuntieae)

BACKGROUND

Plastid genomes (plastomes) have long been recognized as highly conserved in their overall structure, size, gene arrangement and content among land plants. However, recent studies have shown that some lineages present unusual variations in some of these features. Members of the cactus family are one of these lineages, with distinct plastome structures reported across disparate lineages, including gene losses, inversions, boundary movements or loss of the canonical inverted repeat (IR) region. However, only a small fraction of cactus diversity has been analysed so far.

METHODS

Here, we investigated plastome features of the tribe Opuntieae, the remarkable prickly pear cacti, which represent one of the most diverse and important lineages of Cactaceae. We assembled de novo the plastome of 43 species, representing a comprehensive sampling of the tribe, including all seven genera, and analysed their evolution in a phylogenetic comparative framework. Phylogenomic analyses with different datasets (full plastome sequences and genes only) were performed, followed by congruence analyses to assess signals underlying contentious nodes.

KEY RESULTS

Plastomes varied considerably in length, from 121 to 162 kbp, with striking differences in the content and size of the IR region (contraction and expansion events), including a lack of the canonical IR in some lineages and the pseudogenization or loss of some genes. Overall, nine different types of plastomes were reported, deviating in the presence of the IR region or the genes contained in the IR. Overall, plastome sequences resolved phylogenetic relationships within major clades of Opuntieae with high bootstrap values but presented some contentious nodes depending on the dataset analysed (e.g. whole plastome vs. genes only). Congruence analyses revealed that most plastidial regions lack phylogenetic resolution, while few markers are supporting the most likely topology. Likewise, alternative topologies are driven by a handful of plastome markers, suggesting recalcitrant nodes in the phylogeny.

CONCLUSIONS

Our study reveals a dynamic nature of plastome evolution across closely related lineages, shedding light on peculiar features of plastomes. Variation of plastome types across Opuntieae is remarkable in size, structure and content and can be important for the recognition of species in some major clades. Unravelling connections between the causes of plastome variation and the consequences for species biology, physiology, ecology, diversification and adaptation is a promising and ambitious endeavour in cactus research. Although plastome data resolved major phylogenetic relationships, the generation of nuclear genomic data is necessary to confront these hypotheses and assess the recalcitrant nodes further.

Characterization and comparative analysis of the complete plastid genomes of four Astragalus species

Astragalus is the largest flowering plant genus. We assembled the plastid genomes of four Astragalus species (Astragalus iranicus, A. macropelmatus, A. mesoleios, A. odoratus) using next-generation sequencing and analyzed their plastomes including genome organization, codon usage, nucleotide diversity, prediction of RNA editing and etc. The total length of the newly sequenced Astragalus plastomes ranged from 121,050 bp to 123,622 bp, with 110 genes comprising 76 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes. Comparative analysis of the chloroplast genomes of Astragalus revealed several hypervariable regions comprising three non-coding sites (trnQ(UUG)-accD, rps7 -trnV(GAC) and trnR(ACG)-trnN(GUU)) and four protein-coding genes (ycf1, ycf2, accD and clpP), which have potential as molecular markers. Positive selection signatures were found in five genes in Astragalus species including rps11, rps15, accD, clpP and ycf1. The newly sequenced species, A. macropelmatus, has an approximately 13-kb inversion in IR region. Phylogenetic analysis based on 75 protein-coding gene sequences confirmed that Astragalus form a monophyletic clade within the tribe Galegeae and Oxytropis is sister group to the Coluteoid clade. The results of this study may helpful in elucidating the chloroplast genome structure, understanding the evolutionary dynamics at genus Astragalus and IRLC levels and investigating the phylogenetic relationships. Moreover, the newly plastid genomes sequenced have been increased the plastome data resources on Astragalus that can be useful in further phylogenomic studies.

Comparative plastome genomics, taxonomic delimitation and evolutionary divergences of Tetraena hamiensis var. qatarensis and Tetraena simplex (Zygophyllaceae)

The Zygophyllum and Tetraena genera are intriguingly important ecologically and medicinally. Based on morphological characteristics, T. hamiensis var. qatarensis, and T. simplex were transferred from Zygophyllum to Tetraena with the least genomic datasets available. Hence, we sequenced the T. hamiensis and T. simplex and performed in-depth comparative genomics, phylogenetic analysis, and estimated time divergences. The complete plastomes ranged between 106,720 and 106,446 bp-typically smaller than angiosperms plastomes. The plastome circular genomes are divided into large single-copy regions (~ 80,964 bp), small single-copy regions (~ 17,416 bp), and two inverted repeats regions (~ 4170 bp) in both Tetraena species. An unusual shrinkage of IR regions 16-24 kb was identified. This resulted in the loss of 16 genes, including 11 ndh genes which encode the NADH dehydrogenase subunits, and a significant size reduction of Tetraena plastomes compared to other angiosperms. The inter-species variations and similarities were identified using genome-wide comparisons. Phylogenetic trees generated by analyzing the whole plastomes, protein-coding genes, matK, rbcL, and cssA genes exhibited identical topologies, indicating that both species are sisters to the genus Tetraena and may not belong to Zygophyllum. Similarly, based on the entire plastome and proteins coding genes datasets, the time divergence of Zygophyllum and Tetraena was 36.6 Ma and 34.4 Ma, respectively. Tetraena stem ages were 31.7 and 18.2 Ma based on full plastome and protein-coding genes. The current study presents the plastome as a distinguishing and identification feature among the closely related Tetraena and Zygophyllum species. It can be potentially used as a universal super-barcode for identifying plants.

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 comparative study across the parasitic plants of Cuscuta subgenus Grammica (Convolvulaceae) reveals a possible loss of the plastid genome in its section Subulatae

Most species in Cuscuta subgenus Grammica retain many photosynthesis-related plastid genes, generally under purifying selection. A group of holoparasitic species in section Subulatae may have lost their plastid genomes entirely. The c. 153 species of plants belonging to Cuscuta subgenus Grammica are all obligate stem parasites. However, some have completely lost the ability to conduct photosynthesis while others retain photosynthetic machinery and genes. The plastid genome that primarily encodes key photosynthesis genes functions as a bellwether for how reliant plants are on primary production. This research assembles and analyses 17 plastomes across Cuscuta subgenus Grammica with the aim of characterizing the state of the plastome in each of its sections. By comparing the structure and content of plastid genomes across the subgenus, as well as by quantifying the selection acting upon each gene, we reconstructed the patterns of plastome change within the phylogenetic context for this group. We found that species in 13 of the 15 sections that comprise Grammica retain the bulk of plastid photosynthesis genes and are thus hemiparasitic. The complete loss of photosynthesis can be traced to two clades: the entire section Subulatae and a complex of three species within section Ceratophorae. We were unable to recover any significant plastome sequences from section Subulatae, suggesting that plastomes in these species are either drastically reduced or lost entirely.

Comparative Analyses of Chloroplast Genomes for Parasitic Species of Santalales in the Light of Two Newly Sequenced Species, Taxillus nigrans and Scurrula parasitica

When a flowering plant species changes its life history from self-supply to parasite, its chloroplast genomes may have experienced functional physical reduction, and gene loss. Most species of Santalales are hemiparasitic and few studies focus on comparing the chloroplast genomes of the species from this order. In this study, we collected and compared chloroplast genomes of 12 species of Santalales and sequenced the chloroplast genomes of Taxillus nigrans and Scurrula parasitica for the first time. The chloroplast genomes for these species showed typical quadripartite structural organization. Phylogenetic analysis suggested that these 12 species of Santalales clustered into three clades: Viscum (4 spp.) and Osyris (1 sp.) in the Santalaceae and Champereia (1 sp.) in the Opiliaceae formed one clade, while Taxillus (3 spp.) and Scurrula (1 sp.) in the Loranthaceae and Schoepfia (1 sp.) in the Schoepfiaceae formed another clade. Erythropalum (1 sp.), in the Erythropalaceae, appeared as a third, most distant, clade within the Santalales. In addition, both Viscum and Taxillus are monophyletic, and Scurrula is sister to Taxillus. A comparative analysis of the chloroplast genome showed differences in genome size and the loss of genes, such as the ndh genes, infA genes, partial ribosomal genes, and tRNA genes. The 12 species were classified into six categories by the loss, order, and structure of genes in the chloroplast genome. Each of the five genera (Viscum, Osyris, Champereia, Schoepfia, and Erythropalum) represented an independent category, while the three Taxillus species and Scurrula were classified into a sixth category. Although we found that different genes were lost in various categories, most genes related to photosynthesis were retained in the 12 species. Hence, the genetic information accorded with observations that they are hemiparasitic species. Our comparative genomic analyses can provide a new case for the chloroplast genome evolution of parasitic species.

Ascorbate peroxidase postcold regulation of chloroplast NADPH dehydrogenase activity controls cold memory

Exposure of Arabidopsis (Arabidopsis thaliana) to 4°C imprints a cold memory that modulates gene expression in response to a second (triggering) stress stimulus applied several days later. Comparison of plastid transcriptomes of cold-primed and control plants directly before they were exposed to the triggering stimulus showed downregulation of several subunits of chloroplast NADPH dehydrogenase (NDH) and regulatory subunits of ATP synthase. NDH is, like proton gradient 5 (PGR5)-PGR5-like1 (PGRL1), a thylakoid-embedded, ferredoxin-dependent plastoquinone reductase that protects photosystem I and stabilizes ATP synthesis by cyclic electron transport (CET). Like PGRL1A and PGRL1B transcript levels, ndhA and ndhD transcript levels decreased during the 24-h long priming cold treatment. PGRL1 transcript levels were quickly reset in the postcold phase, but expression of ndhA remained low. The transcript abundances of other ndh genes decreased within the next days. Comparison of thylakoid-bound ascorbate peroxidase (tAPX)-free and transiently tAPX-overexpressing or tAPX-downregulating Arabidopsis lines demonstrated that ndh expression is suppressed by postcold induction of tAPX. Four days after cold priming, when tAPX protein accumulation was maximal, NDH activity was almost fully lost. Lack of the NdhH-folding chaperonin Crr27 (Cpn60β4), but not lack of the NDH activity modulating subunits NdhM, NdhO, or photosynthetic NDH subcomplex B2 (PnsB2), strengthened priming regulation of zinc finger of A. thaliana 10, which is a nuclear-localized target gene of the tAPX-dependent cold-priming pathway. We conclude that cold-priming modifies chloroplast-to-nucleus stress signaling by tAPX-mediated suppression of NDH-dependent CET and that plastid-encoded NdhH, which controls subcomplex A assembly, is of special importance for memory stabilization.

Analysis of the complete plastomes and nuclear ribosomal DNAs from Euonymus hamiltonianus and its relatives sheds light on their diversity and evolution

Euonymus hamiltonianus and its relatives (Celastraceae family) are used for ornamental and medicinal purposes. However, species identification in Euonymus is difficult due to their morphological diversity. Using plastid genome (plastome) data, we attempt to reveal phylogenetic relationship among Euonymus species and develop useful markers for molecular identification. We assembled the plastome and nuclear ribosomal DNA (nrDNA) sequences from five Euonymus lines collected from South Korea: three Euonymus hamiltonianus accessions, E. europaeus, and E. japonicus. We conducted an in-depth comparative analysis using ten plastomes, including other publicly available plastome data for this genus. The genome structures, gene contents, and gene orders were similar in all Euonymus plastomes in this study. Analysis of nucleotide diversity revealed six divergence hotspots in their plastomes. We identified 339 single nucleotide polymorphisms and 293 insertion or deletions among the four E. hamiltonianus plastomes, pointing to abundant diversity even within the same species. Among 77 commonly shared genes, 9 and 33 were identified as conserved genes in the genus Euonymus and E. hamiltonianus, respectively. Phylogenetic analysis based on plastome and nrDNA sequences revealed the overall consensus and relationships between plastomes and nrDNAs. Finally, we developed six barcoding markers and successfully applied them to 31 E. hamiltonianus lines collected from South Korea. Our findings provide the molecular basis for the classification and molecular taxonomic criteria for the genus Euonymus (at least in Korea), which should aid in more objective classification within this genus. Moreover, the newly developed markers will be useful for understanding the species delimitation of E. hamiltonianus and closely related species.

New plastome structural rearrangements discovered in core Tillandsioideae (Bromeliaceae) support recently adopted taxonomy

Full plastome sequences for land plants have become readily accessible thanks to the development of Next Generation Sequencing (NGS) techniques and powerful bioinformatic tools. Despite this vast amount of genomic data, some lineages remain understudied. Full plastome sequences from the highly diverse (>1,500 spp.) subfamily Tillandsioideae (Bromeliaceae, Poales) have been published for only three (i.e., Guzmania, Tillandsia, and Vriesea) out of 22 currently recognized genera. Here, we focus on core Tillandsioideae, a clade within subfamily Tillandsioideae, and explore the contribution of individual plastid markers and data categories to inform deep divergences of a plastome phylogeny. We generated 37 high quality plastome assemblies and performed a comparative analysis in terms of plastome structure, size, gene content and order, GC content, as well as number and type of repeat motifs. Using the obtained phylogenetic context, we reconstructed the evolution of these plastome attributes and assessed if significant shifts on the evolutionary traits' rates have occurred in the evolution of the core Tillandsioideae. Our results agree with previously published phylogenetic hypotheses based on plastid data, providing stronger statistical support for some recalcitrant nodes. However, phylogenetic discordance with previously published nuclear marker-based hypotheses was found. Several plastid markers that have been consistently used to address phylogenetic relationships within Tillandsioideae were highly informative for the retrieved plastome phylogeny and further loci are here identified as promising additional markers for future studies. New lineage-specific plastome rearrangements were found to support recently adopted taxonomic groups, including large inversions, as well as expansions and contractions of the inverted repeats. Evolutionary trait rate shifts associated with changes in size and GC content of the plastome regions were found across the phylogeny of core Tillandsioideae.

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.

Molecular Phylogeny and Historical Biogeography of Goodyera R. Br. (Orchidaceae): A Case of the Vicariance Between East Asia and North America

Understanding of intercontinental distribution in the Northern Hemisphere has attracted a lot of attention from botanists. However, although Orchidaceae is the largest group of angiosperms, biogeographical studies on the disjunctive pattern have not been sufficient for this family. Goodyera R. Br. (tribe Cranichideae, subfamily Orchidoideae, family Orchidaceae) is widely distributed in temperate and tropical regions. Although the phylogenetic relationship of Goodyera inferred from both morphological and molecular data has been conducted, the sampled taxa were mainly distributed in Asia regions that resulted in non-monophyly of this genus. In this study, the complete plastid genomes of Goodyera, generated by next-generation sequencing (NGS) technique and sampled in East Asia and North America, were used to reconstruct phylogeny and explore the historical biogeography. A total of 18 Goodyera species including seven newly sequenced species were analyzed. Based on 79 protein-coding genes, the phylogenetic analysis revealed that Goodyera could be subdivided into four subclades with high support values. The polyphyletic relationships among Goodyera taxa were confirmed, and the unclear position of G. foliosa was also resolved. The datasets that are composed of the 14 coding sequences (CDS) (matK, atpF, ndhK, accD, cemA, clpP, rpoA, rpl22, ndhF, ccsA, ndhD, ndhI, ndhA, and ycf 1) showed the same topology derived from 79 protein-coding genes. Molecular dating analyses revealed the origin of Goodyera in the mid-Miocene (15.75 Mya). Nearctic clade of Goodyera was diverged at 10.88 Mya from their most recent common ancestor (MRCA). The biogeographical reconstruction suggests that subtropical or tropical Asia is the origin of Goodyera and it has subsequently spread to temperate Asia during the Miocene. In addition, Nearctic clade is derived from East Asian species through Bering Land Bridge (BLB) during the Miocene. The speciation of Goodyera is most likely to have occurred during Miocene, and climatic and geological changes are thought to have had a part in this diversification. Our findings propose both origin and vicariance events of Goodyera for the first time and add an example for the biogeographical history of the Northern Hemisphere.

Plastid RNA editing reduction accompanied with genetic variations in Cymbidium, a genus with diverse lifestyle modes

Recent sequencing efforts have broadly uncovered the evolutionary trajectory of plastid genomes (plastomes) of flowering plants in diverse habitats, yet our knowledge of the evolution of plastid posttranscriptional modifications is limited. In this study, we generated 11 complete plastomes and performed ultra-deep transcriptome sequencing to investigate the co-evolution of plastid RNA editing and genetic variation in Cymbidium, a genus with diverse trophic lifestyles. Genome size and gene content is reduced in terrestrial and green mycoheterotrophic orchids relative to their epiphytic relatives. This could be partly due to extensive losses and pseudogenization of ndh genes for the plastid NADH dehydrogenase-like complex, but independent pseudogenization of ndh genes has also occurred in the epiphyte C. mannii, which was reported to use strong crassulacean acid metabolism photosynthesis. RNA editing sites are abundant but variable in number among Cymbidium plastomes. The nearly twofold variation in editing abundance is mainly due to extensive reduction of ancestral editing sites in ndh transcripts of terrestrial, mycoheterotrophic, and C. mannii plastomes. The co-occurrence of editing reduction and pseudogenization in ndh genes suggests functional constraints on editing machinery may be relaxed, leading to nonrandom loss of ancestral edited sites via reduced editing efficiency. This study represents the first systematic examination of RNA editing evolution linked to plastid genome variation in a single genus. We also propose an explanation for how genomic and posttranscriptional variations might be affected by lifestyle-associated ecological adaptation strategies in Cymbidium.

The complete chloroplast genome sequencing analysis revealed an unusual IRs reduction in three species of subfamily Zygophylloideae

Tetraena mongolica, Zygophyllum xanthoxylon, and Z. fabago are three typical dryland plants with important ecological values in subfamily Zygophylloideae of Zygophyllaceae. Studies on the chloroplast genomes of them are favorable for understanding the diversity and phylogeny of Zygophyllaceae. Here, we sequenced and assembled the whole chloroplast genomes of T. mongolica, Z. xanthoxylon, and Z. fabago, and performed comparative genomic and phylogenetic analysis. The total size, structure, gene content and orders of these three chloroplast genomes were similar, and the three chloroplast genomes exhibited a typical quadripartite structure with a large single-copy region (LSC; 79,696–80,291 bp), a small single-copy region (SSC; 16,462–17,162 bp), and two inverted repeats (IRs; 4,288–4,413 bp). A total of 107 unique genes were identified from the three chloroplast genomes, including 70 protein-coding genes, 33 tRNAs, and 4 rRNAs. Compared with other angiosperms, the three chloroplast genomes were significantly reduced in overall length due to an unusual 16–24 kb shrinkage of IR regions and loss of the 11 genes which encoded subunits of NADH dehydrogenase. Genome-wide comparisons revealed similarities and variations between the three species and others. Phylogenetic analysis based on the three chloroplast genomes supported the opinion that Zygophyllaceae belonged to Zygophyllales in Fabids, and Z. xanthoxylon and Z. fabago belonged to Zygophyllum. The genome-wide comparisons revealed the similarity and variations between the chloroplast genomes of the three Zygophylloideae species and other plant species. This study provides a valuable molecular biology evidence for further studies of phylogenetic status of Zygophyllaceae.

The complete chloroplast genome of Onobrychis gaubae (Fabaceae-Papilionoideae): comparative analysis with related IR-lacking clade species

BACKGROUND

Plastome (Plastid genome) sequences provide valuable markers for surveying evolutionary relationships and population genetics of plant species. Papilionoideae (papilionoids) has different nucleotide and structural variations in plastomes, which makes it an ideal model for genome evolution studies. Therefore, by sequencing the complete chloroplast genome of Onobrychis gaubae in this study, the characteristics and evolutionary patterns of plastome variations in IR-loss clade were compared.

RESULTS

In the present study, the complete plastid genome of O. gaubae, endemic to Iran, was sequenced using Illumina paired-end sequencing and was compared with previously known genomes of the IRLC species of legumes. The O. gaubae plastid genome was 122,688 bp in length and included a large single-copy (LSC) region of 81,486 bp, a small single-copy (SSC) region of 13,805 bp and one copy of the inverted repeat (IR_b) of 29,100 bp. The genome encoded 110 genes, including 76 protein-coding genes, 30 transfer RNA (tRNA) genes and four ribosome RNA (rRNA) genes and possessed 83 simple sequence repeats (SSRs) and 50 repeated structures with the highest proportion in the LSC. Comparative analysis of the chloroplast genomes across IRLC revealed three hotspot genes (ycf1, ycf2, clpP) which could be used as DNA barcode regions. Moreover, seven hypervariable regions [trnL(UAA)-trnT(UGU), trnT(GGU)-trnE(UUC), ycf1, ycf2, ycf4, accD and clpP] were identified within Onobrychis, which could be used to distinguish the Onobrychis species. Phylogenetic analyses revealed that O. gaubae is closely related to Hedysarum. The complete O. gaubae genome is a valuable resource for investigating evolution of Onobrychis species and can be used to identify related species.

CONCLUSIONS

Our results reveal that the plastomes of the IRLC are dynamic molecules and show multiple gene losses and inversions. The identified hypervariable regions could be used as molecular markers for resolving phylogenetic relationships and species identification and also provide new insights into plastome evolution across IRLC.

Comparative Analyses of 3,654 Plastid Genomes Unravel Insights Into Evolutionary Dynamics and Phylogenetic Discordance of Green Plants

The plastid organelle is essential for many vital cellular processes and the growth and development of plants. The availability of a large number of complete plastid genomes could be effectively utilized to understand the evolution of the plastid genomes and phylogenetic relationships among plants. We comprehensively analyzed the plastid genomes of Viridiplantae comprising 3,654 taxa from 298 families and 111 orders and compared the genomic organizations in their plastid genomic DNA among major clades, which include gene gain/loss, gene copy number, GC content, and gene blocks. We discovered that some important genes that exhibit similar functions likely formed gene blocks, such as the psb family presumably showing co-occurrence and forming gene blocks in Viridiplantae. The inverted repeats (IRs) in plastid genomes have doubled in size across land plants, and their GC content is substantially higher than non-IR genes. By employing three different data sets [all nucleotide positions (nt123), only the first and second codon positions (nt12), and amino acids (AA)], our phylogenomic analyses revealed Chlorokybales + Mesostigmatales as the earliest-branching lineage of streptophytes. Hornworts, mosses, and liverworts forming a monophylum were identified as the sister lineage of tracheophytes. Based on nt12 and AA data sets, monocots, Chloranthales and magnoliids are successive sister lineages to the eudicots + Ceratophyllales clade. The comprehensive taxon sampling and analysis of different data sets from plastid genomes recovered well-supported relationships of green plants, thereby contributing to resolving some long-standing uncertainties in the plant phylogeny.

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.

Unprecedented Intraindividual Structural Heteroplasmy in Eleocharis (Cyperaceae, Poales) Plastomes

Plastid genomes (plastomes) of land plants have a conserved quadripartite structure in a gene-dense unit genome consisting of a large inverted repeat that separates two single copy regions. Recently, alternative plastome structures were suggested in Geraniaceae and in some conifers and Medicago the coexistence of inversion isomers has been noted. In this study, plastome sequences of two Cyperaceae, Eleocharis dulcis (water chestnut) and Eleocharis cellulosa (gulf coast spikerush), were completed. Unlike the conserved plastomes in basal groups of Poales, these Eleocharis plastomes have remarkably divergent features, including large plastome sizes, high rates of sequence rearrangements, low GC content and gene density, gene duplications and losses, and increased repetitive DNA sequences. A novel finding among these features was the unprecedented level of heteroplasmy with the presence of multiple plastome structural types within a single individual. Illumina paired-end assemblies combined with PacBio single-molecule real-time sequencing, long-range polymerase chain reaction, and Sanger sequencing data identified at least four different plastome structural types in both Eleocharis species. PacBio long read data suggested that one of the four E. dulcis plastome types predominates.

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.

Highly variable chloroplast genome from two endangered Papaveraceae lithophytes Corydalis tomentella and Corydalis saxicola

The increasingly wide application of chloroplast (cp) genome super-barcode in taxonomy and the recent breakthrough in cp genetic engineering make the development of new cp gene resources urgent and significant. Corydalis is recognized as the most genotypes complicated and taxonomically challenging plant taxa in Papaveraceae. However, there currently are few reports about cp genomes of the genus Corydalis. In this study, we sequenced four complete cp genomes of two endangered lithophytes Corydalis saxicola and Corydalis tomentella in Corydalis, conducted a comparison of these cp genomes among each other as well as with others of Papaveraceae. The cp genomes have a large genome size of 189,029-190,247 bp, possessing a quadripartite structure and with two highly expanded inverted repeat (IR) regions (length: 41,955-42,350 bp). Comparison between the cp genomes of C. tomentella, C. saxicola, and Papaveraceae species, five NADH dehydrogenase-like genes (ndhF, ndhD, ndhL, ndhG, and ndhE) with psaC, rpl32, ccsA, and trnL-UAG normally located in the SSC region have migrated to IRs, resulting in IR expansion and gene duplication. An up to 9 kb inversion involving five genes (rpl23, ycf2, ycf15, trnI-CAU, and trnL-CAA) was found within IR regions. The accD gene was found to be absent and the ycf1 gene has shifted from the IR/SSC border to the SSC region as a single copy. Phylogenetic analysis based on the sequences of common CDS showed that the genus Corydalis is quite distantly related to the other genera of Papaveraceae, it provided a new clue for recent advocacy to establish a separate Fumariaceae family. Our results revealed one special cp genome structure in Papaveraceae, provided a useful resources for classification of the genus Corydalis, and will be valuable for understanding Papaveraceae evolutionary relationships.

Lineage-specific plastid degradation in subtribe Gentianinae (Gentianaceae)

The structure and sequence of plastid genomes is highly conserved across most land plants, except for a minority of lineages that show gene loss and genome degradation. Understanding the early stages of plastome degradation may provide crucial insights into the repeatability and predictability of genomic evolutionary trends. We investigated these trends in subtribe Gentianinae of the Gentianaceae, which encompasses ca. 450 species distributed around the world, particularly in alpine and subalpine environments. We sequenced, assembled, and annotated the plastomes of 41 species, representing all six genera in subtribe Gentianinae and all main sections of the species-rich genus Gentiana L. We reconstructed the phylogeny, estimated divergence times, investigated the phylogenetic distribution of putative gene losses, and related these to substitution rate shifts and species' habitats. We obtained a strongly supported topology consistent with earlier studies, with all six genera in Gentianinae recovered as monophyletic and all main sections of Gentiana having full support. While closely related species have very similar plastomes in terms of size and structure, independent gene losses, particularly of the ndh complex, have occurred in multiple clades across the phylogeny. Gene loss was usually associated with a shift in the boundaries of the small single-copy and inverted repeat regions. Substitution rates were variable between clades, with evidence for both elevated and decelerated rate shifts. Independent lineage-specific loss of ndh genes occurred at a wide range of times, from Eocene to Pliocene. Our study illustrates that diverse degradation patterns shape the evolution of the plastid in this species-rich plant group.

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.

Plastid Genomes of Flowering Plants: Essential Principles

The plastid genome (plastome ) has proved a valuable source of data for evaluating evolutionary relationships among angiosperms. Through basic and applied approaches, plastid transformation technology offers the potential to understand and improve plant productivity, providing food, fiber, energy, and medicines to meet the needs of a burgeoning global population. The growing genomic resources available to both phylogenetic and biotechnological investigations is allowing novel insights and expanding the scope of plastome research to encompass new species. In this chapter, we present an overview of some of the seminal and contemporary research that has contributed to our current understanding of plastome evolution and attempt to highlight the relationship between evolutionary mechanisms and the tools of plastid genetic engineering.

Regulation of electron transport is essential for photosystem I stability and plant growth

Photosynthetic electron transport is regulated by cyclic and pseudocyclic electron flow (CEF and PCEF) to maintain the balance between light availability and metabolic demands. CEF transfers electrons from photosystem I to the plastoquinone pool with two mechanisms, dependent either on PGR5/PGRL1 or on the type I NADH dehydrogenase-like (NDH) complex. PCEF uses electrons from photosystem I to reduce oxygen and in many groups of photosynthetic organisms, but remarkably not in angiosperms, it is catalyzed by flavodiiron proteins (FLVs). In this study, Physcomitrella patens plants depleted in PGRL1, NDH and FLVs in different combinations were generated and characterized, showing that all these mechanisms are active in this moss. Surprisingly, in contrast to flowering plants, Physcomitrella patens can cope with the simultaneous inactivation of PGR5- and NDH-dependent CEF but, when FLVs are also depleted, plants show strong growth reduction and photosynthetic activity is drastically reduced. The results demonstrate that mechanisms for modulation of photosynthetic electron transport have large functional overlap but are together indispensable to protect photosystem I from damage and they are an essential component for photosynthesis in any light regime.

Extensive variation in nucleotide substitution rate and gene/intron loss in mitochondrial genomes of Pelargonium

Geraniaceae organelle genomes have been shown to exhibit several highly unusual features compared to most other photosynthetic angiosperms. This includes massively rearranged plastomes with considerable size variation, extensive gene and intron loss, accelerated rates of nucleotide substitutions in both mitogenomes and plastomes, and biparental inheritance and cytonuclear incompatibility of the plastome. Most previous studies have focused on plastome evolution with mitogenome comparisons limited to only a few taxa or genes. In this study, mitogenomes and transcriptomes were examined for 27 species of Geraniales, including 13 species of Pelargonium. Extensive gene and intron losses were detected across the Geraniales with Pelargonium representing the most gene depauperate lineage in the family. Plotting these events on the Geraniaceae phylogenetic tree showed that gene losses occurred multiple times, whereas intron losses more closely reflected the relationships among taxa. In addition, P. australe acquired an intron by horizontal transfer. Comparisons of nucleotide substitution rates in Pelargonium showed that synonymous changes in nuclear genes were much lower than in mitochondrial genes. This is in contrast to the previously published studies that indicated that nuclear genes have 16 fold higher rates than mitochondrial genes across angiosperms. Elevated synonymous substitutions occurred for each mitochondrial gene in Pelargonium with the highest values 783 and 324 times higher than outgroups and other Geraniaceae, respectively. Pelargonium is one of four unrelated genera of angiosperms (Ajuga, Plantago and Silene) that have experienced highly accelerated nucleotide substitutions in mitogenomes. It is distinct from most angiosperms in also having elevated substitution rates in plastid genes but the cause of rate accelerations in Pelargonium plastomes and mitogenomes may be different.

The chloroplast NADH dehydrogenase-like complex influences the photosynthetic activity of the moss Physcomitrella patens

Alternative electron pathways contribute to regulation of photosynthetic light reactions to adjust to metabolic demands in dynamic environments. The chloroplast NADH dehydrogenase-like (NDH) complex mediates the cyclic electron transport pathway around PSI in different cyanobacteria, algae, and plant species, but it is not fully conserved in all photosynthetic organisms. In order to assess how the physiological role of this complex changed during plant evolution, we isolated Physcomitrella patens lines knocked out for the NDHM gene that encodes a subunit fundamental for the activity of the complex. ndhm knockout mosses indicated high PSI acceptor side limitation upon abrupt changes in illumination. In P. patens, pseudo-cyclic electron transport mediated by flavodiiron proteins (FLVs) was also shown to prevent PSI over-reduction in plants exposed to light fluctuations. flva ndhm double knockout mosses had altered photosynthetic performance and growth defects under fluctuating light compared with the wild type and single knockout mutants. The results showed that while the contribution of NDH to electron transport is minor compared with FLV, NDH still participates in modulating photosynthetic activity, and it is critical to avoid PSI photoinhibition, especially when FLVs are inactive. The functional overlap between NDH- and FLV-dependent electron transport supports PSI activity and prevents its photoinhibition under light variations.

The plastid NAD(P)H dehydrogenase-like complex: structure, function and evolutionary dynamics

The thylakoid NAD(P)H dehydrogenase-like (NDH) complex is a large protein complex that reduces plastoquinone and pumps protons into the lumen generating protonmotive force. In plants, the complex consists of both nuclear and chloroplast-encoded subunits. Despite its perceived importance for stress tolerance and ATP generation, chloroplast-encoded NDH subunits have been lost numerous times during evolution in species occupying seemingly unrelated environmental niches. We have generated a phylogenetic tree that reveals independent losses in multiple phylogenetic lineages, and we use this tree as a reference to discuss possible evolutionary contexts that may have relaxed selective pressure for retention of ndh genes. While we are still yet unable to pinpoint a singular specific lifestyle that negates the need for NDH, we are able to rule out several long-standing explanations. In light of this, we discuss the biochemical changes that would be required for the chloroplast to dispense with NDH functionality with regards to known and proposed NDH-related reactions.

Degradation of key photosynthetic genes in the critically endangered semi-aquatic flowering plant Saniculiphyllum guangxiense (Saxifragaceae)

BACKGROUND

Plastid gene loss and pseudogenization has been widely documented in parasitic and mycoheterotrophic plants, which have relaxed selective constraints on photosynthetic function. More enigmatic are sporadic reports of pseudogenization and loss of important photosynthesis genes in lineages thought to be fully photosynthetic. Here we report the complete plastid genome of Saniculiphyllum guangxiense, a critically endangered and phylogenetically isolated plant lineage, along with genomic evidence of reduced chloroplast function. We also report 22 additional plastid genomes representing the diversity of its containing clade Saxifragales, characterizing gene content and placing variation in a broader phylogenetic context.

RESULTS

We find that the plastid genome of Saniculiphyllum has experienced pseudogenization of five genes of the ndh complex (ndhA, ndhB, ndhD, ndhF, and ndhK), previously reported in flowering plants with an aquatic habit, as well as the surprising pseudogenization of two genes more central to photosynthesis (ccsA and cemA), contrasting with strong phylogenetic conservatism of plastid gene content in all other sampled Saxifragales. These genes participate in photooxidative protection, cytochrome synthesis, and carbon uptake. Nuclear paralogs exist for all seven plastid pseudogenes, yet these are also unlikely to be functional.

CONCLUSIONS

Saniculiphyllum appears to represent the greatest degree of plastid gene loss observed to date in any fully photosynthetic lineage, perhaps related to its extreme habitat specialization, yet plastid genome length, structure, and substitution rate are within the variation previously reported for photosynthetic plants. These results highlight the increasingly appreciated dynamism of plastid genomes, otherwise highly conserved across a billion years of green plant evolution, in plants with highly specialized life history traits.

Insights Into Chloroplast Genome Evolution Across Opuntioideae (Cactaceae) Reveals Robust Yet Sometimes Conflicting Phylogenetic Topologies

Chloroplast genomes (plastomes) are frequently treated as highly conserved among land plants. However, many lineages of vascular plants have experienced extensive structural rearrangements, including inversions and modifications to the size and content of genes. Cacti are one of these lineages, containing the smallest plastome known for an obligately photosynthetic angiosperm, including the loss of one copy of the inverted repeat (∼25 kb) and the ndh gene suite, but only a few cacti from the subfamily Cactoideae have been sufficiently characterized. Here, we investigated the variation of plastome sequences across the second-major lineage of the Cactaceae, the subfamily Opuntioideae, to address (1) how variable is the content and arrangement of chloroplast genome sequences across the subfamily, and (2) how phylogenetically informative are the plastome sequences for resolving major relationships among the clades of Opuntioideae. Our de novo assembly of the Opuntia quimilo plastome recovered an organelle of 150,347 bp in length with both copies of the inverted repeat and the presence of all the ndh gene suite. An expansion of the large single copy unit and a reduction of the small single copy unit was observed, including translocations and inversion of genes, as well as the putative pseudogenization of some loci. Comparative analyses among all clades within Opuntioideae suggested that plastome structure and content vary across taxa of this subfamily, with putative independent losses of the ndh gene suite and pseudogenization of genes across disparate lineages, further demonstrating the dynamic nature of plastomes in Cactaceae. Our plastome dataset was robust in resolving three tribes with high support within Opuntioideae: Cylindropuntieae, Tephrocacteae and Opuntieae. However, conflicting topologies were recovered among major clades when exploring different assemblies of markers. A plastome-wide survey for highly informative phylogenetic markers revealed previously unused regions for future use in Sanger-based studies, presenting a valuable dataset with primers designed for continued evolutionary studies across Cactaceae. These results bring new insights into the evolution of plastomes in cacti, suggesting that further analyses should be carried out to address how ecological drivers, physiological constraints and morphological traits of cacti may be related with the common rearrangements in plastomes that have been reported across the family.

The highly efficient NDH-dependent photosystem I cyclic electron flow pathway in the marine angiosperm Zostera marina

Zostera marina, a fully submerged marine angiosperm with a unique evolutionary history associated with its terrestrial origin, has distinct photochemical characteristics caused by its oxygen-evolving complex (OEC) being prone to deactivation in visible light. Based on the present phylogenetic analysis, the chloroplast NADPH dehydrogenase-like (NDH) complex was found to be completed in of Z. marina, unlike other marine plants, suggesting its crucial role. Thus, the responses of electron transport to irradiation were investigated through multiple chlorophyll fluorescence techniques and Western blot analysis. Moreover, the respective contribution of the two photosystem I cyclic electron flow (PSI-CEF) pathways to the generation of trans-thylakoid proton gradient (∆pH) was also examined using inhibitors. The contributions of the two PSI-CEF pathways to ∆pH were similar; furthermore, there was a trade-off between the two pathways under excess irradiation: the PGR5/L1-dependent PSI-CEF decreased gradually following its activation during the initial illumination, while NDH-dependent PSI-CEF was activated gradually with exposure duration. OEC inactivation was continuously under excess irradiation, which exhibits a positive linear correlation with the activation of NDH-dependent PSI-CEF. We suggest that PGR5/L1-dependent PSI-CEF was preferentially activated to handle the excess electron caused by the operation of OEC during the initial illumination. Subsequently, the increasing OEC inactivation with exposure duration resulted in a deficit of electrons. Limited electrons from PSI might preferentially synthesize NADPH, which could support the function of NDH-dependent PSI-CEF to generate ∆pH and ATP via reducing ferredoxin, thereby maintaining OEC stability.

Characterization of sequence variability hotspots in Cranichideae plastomes (Orchidaceae, Orchidoideae)

This study reports complete plastome sequences for six species of Neotropical Cranichideae and focuses on identification of the most variable regions (hotspots) in this group of orchids. These structure of these six plastomes is relatively conserved, exhibiting lengths ranging between 142,599 to 154,562 bp with 36.7% GC on average and exhibiting typical quadripartite arrangement (LSC, SSC and two IRs). Variation detected in the LSC/IR and SSC/IR junctions is explained by the loss of ndhF and ycf1 length variation. For the two genera of epiphytic clade in Spiranthinae, almost whole sets of the ndh-gene family were missing. Eight mutation hotspots were identified based on nucleotide diversity, sequence variability and parsimony-informative sites. Three of them (rps16-trnQ, trnT-trnL, rpl32-trnL) seem to be universal hotspots in the family, and the other five (trnG-trnR, trnR-atpA, trnP-psaJ, rpl32-infA, and rps15-ycf1) are described for the first time as orchid molecular hotspots. These regions have much more variation than all those used previously in phylogenetics of the group and offer useful plastid markers for phylogenetic, barcoding and population genetic studies. The use of whole plastomes or exclusive no-gap matrices also positioned with high support the holomycotrophic Rhizanthella among Orchidoideae plastomes in model-based analyses, showing the utility of plastomes for phylogenetic placement of this unusual genus.

Intraspecific Variation within the Utricularia amethystina Species Morphotypes Based on Chloroplast Genomes

Utricularia amethystina Salzm. ex A.St.-Hil. & Girard (Lentibulariaceae) is a highly polymorphic carnivorous plant taxonomically rearranged many times throughout history. Herein, the complete chloroplast genomes (cpDNA) of three U. amethystina morphotypes: purple-, white-, and yellow-flowered, were sequenced, compared, and putative markers for systematic, populations, and evolutionary studies were uncovered. In addition, RNA-Seq and RNA-editing analysis were employed for functional cpDNA evaluation. The cpDNA of three U. amethystina morphotypes exhibits typical quadripartite structure. Fine-grained sequence comparison revealed a high degree of intraspecific genetic variability in all morphotypes, including an exclusive inversion in the psbM and petN genes in U. amethystina yellow. Phylogenetic analyses indicate that U. amethystina morphotypes are monophyletic. Furthermore, in contrast to the terrestrial Utricularia reniformis cpDNA, the U. amethystina morphotypes retain all the plastid NAD(P)H-dehydrogenase (ndh) complex genes. This observation supports the hypothesis that the ndhs in terrestrial Utricularia were independently lost and regained, also suggesting that different habitats (aquatic and terrestrial) are not related to the absence of Utricularia ndhs gene repertoire as previously assumed. Moreover, RNA-Seq analyses recovered similar patterns, including nonsynonymous RNA-editing sites (e.g., rps14 and petB). Collectively, our results bring new insights into the chloroplast genome architecture and evolution of the photosynthesis machinery in the Lentibulariaceae.

Plastid Genomes of Carnivorous Plants Drosera rotundifolia and Nepenthes × ventrata Reveal Evolutionary Patterns Resembling Those Observed in Parasitic Plants

Carnivorous plants have the ability to capture and digest small animals as a source of additional nutrients, which allows them to grow in nutrient-poor habitats. Here we report the complete sequences of the plastid genomes of two carnivorous plants of the order Caryophyllales, Drosera rotundifolia and Nepenthes × ventrata. The plastome of D. rotundifolia is repeat-rich and highly rearranged. It lacks NAD(P)H dehydrogenase genes, as well as ycf1 and ycf2 genes, and three essential tRNA genes. Intron losses are observed in some protein-coding and tRNA genes along with a pronounced reduction of RNA editing sites. Only six editing sites were identified by RNA-seq in D. rotundifolia plastid genome and at most conserved editing sites the conserved amino acids are already encoded at the DNA level. In contrast, the N. × ventrata plastome has a typical structure and gene content, except for pseudogenization of the ccsA gene. N. × ventrata and D. rotundifolia could represent different stages of evolution of the plastid genomes of carnivorous plants, resembling events observed in parasitic plants in the course of the switch from autotrophy to a heterotrophic lifestyle.

The unique photosynthetic apparatus of Pinaceae: analysis of photosynthetic complexes in Picea abies

Pinaceae are the predominant photosynthetic species in boreal forests, but so far no detailed description of the protein components of the photosynthetic apparatus of these gymnosperms has been available. In this study we report a detailed characterization of the thylakoid photosynthetic machinery of Norway spruce (Picea abies (L.) Karst). We first customized a spruce thylakoid protein database from translated transcript sequences combined with existing protein sequences derived from gene models, which enabled reliable tandem mass spectrometry identification of P. abies thylakoid proteins from two-dimensional large pore blue-native/SDS-PAGE. This allowed a direct comparison of the two-dimensional protein map of thylakoid protein complexes from P. abies with the model angiosperm Arabidopsis thaliana. Although the subunit composition of P. abies core PSI and PSII complexes is largely similar to that of Arabidopsis, there was a high abundance of a smaller PSI subcomplex, closely resembling the assembly intermediate PSI* complex. In addition, the evolutionary distribution of light-harvesting complex (LHC) family members of Pinaceae was compared in silico with other land plants, revealing that P. abies and other Pinaceae (also Gnetaceae and Welwitschiaceae) have lost LHCB4, but retained LHCB8 (formerly called LHCB4.3). The findings reported here show the composition of the photosynthetic apparatus of P. abies and other Pinaceae members to be unique among land plants.

The Unique Evolutionary Trajectory and Dynamic Conformations of DR and IR/DR-Coexisting Plastomes of the Early Vascular Plant Selaginellaceae (Lycophyte)

Both direct repeats (DR) and inverted repeats (IR) are documented in the published plastomes of Selaginella species indicating the unusual and diverse plastome structure in the family Selaginellaceae. In this study, we newly sequenced complete plastomes of seven species from five main lineages of Selaginellaceae and also resequenced three species (Selaginella tamariscina, Selaginella uncinata, and Selaginella moellendorffii) to explore the evolutionary trajectory of Selaginellaceae plastomes. Our results showed that the plastomes of Selaginellaceae vary remarkably in size, gene contents, gene order, and GC contents. Notably, both DR and IR structures existed in the plastomes of Selaginellaceae with DR structure being an ancestral state. The occurrence of DR structure was at ∼257 Ma and remained in most subgenera of Selaginellaceae, whereas IR structure only reoccurred in Selaginella sect. Lepidophyllae (∼143 Ma) and Selaginella subg. Heterostachys (∼19 Ma). The presence of a pair of large repeats psbK-trnQ, together with DR/IR region in Selaginella bisulcata, Selaginella pennata, S. uncinata, and Selaginella hainanensis, could frequently mediate diverse homologous recombination and create approximately equal stoichiometric isomers (IR/DR-coexisting) and subgenomes. High proportion of repeats is presumably responsible for the dynamic IR/DR-coexisting plastomes, which possess a lower synonymous substitution rate (dS) compared with DR-possessing and IR-possessing plastomes. We propose that the occurrence of DR structure, together with few repeats, is possibly selected to keep the stability of plastomes and the IR/DR-coexisting plastomes also reached an equilibrium in plastome organization through highly efficient homologous recombination to maintain stability.

Lycophyte plastid genomics: extreme variation in GC, gene and intron content and multiple inversions between a direct and inverted orientation of the rRNA repeat

Lycophytes are a key group for understanding vascular plant evolution. Lycophyte plastomes are highly distinct, indicating a dynamic evolutionary history, but detailed evaluation is hindered by the limited availability of sequences. Eight diverse plastomes were sequenced to assess variation in structure and functional content across lycophytes. Lycopodiaceae plastomes have remained largely unchanged compared with the common ancestor of land plants, whereas plastome evolution in Isoetes and especially Selaginella is highly dynamic. Selaginella plastomes have the highest GC content and fewest genes and introns of any photosynthetic land plant. Uniquely, the canonical inverted repeat was converted into a direct repeat (DR) via large-scale inversion in some Selaginella species. Ancestral reconstruction identified additional putative transitions between an inverted and DR orientation in Selaginella and Isoetes plastomes. A DR orientation does not disrupt the activity of copy-dependent repair to suppress substitution rates within repeats. Lycophyte plastomes include the most archaic examples among vascular plants and the most reconfigured among land plants. These evolutionary trends correlate with the mitochondrial genome, suggesting shared underlying mechanisms. Copy-dependent repair for DR-localized genes indicates that recombination and gene conversion are not inhibited by the DR orientation. Gene relocation in lycophyte plastomes occurs via overlapping inversions rather than transposase/recombinase-mediated processes.

Plastome-Wide Rearrangements and Gene Losses in Carnivorous Droseraceae

The plastid genomes of four related carnivorous plants (Drosera regia, Drosera erythrorhiza, Aldrovanda vesiculosa, and Dionaea muscipula) were sequenced to examine changes potentially induced by the transition to carnivory. The plastid genomes of the Droseraceae show multiple rearrangements, gene losses, and large expansions or contractions of the inverted repeat. All the ndh genes are lost or nonfunctional, as well as in some of the species, clpP1, ycf1, ycf2 and some tRNA genes. Uniquely, among land plants, the trnK gene has no intron. Carnivory in the Droseraceae coincides with changes in plastid gene content similar to those induced by parasitism and mycoheterotrophy, suggesting parallel changes in chloroplast function due to the similar switch from autotrophy to (mixo-) heterotrophy. A molecular phylogeny of the taxa based on all shared plastid genes indicates that the "snap-traps" of Aldrovanda and Dionaea have a common origin.

Balancing protection and efficiency in the regulation of photosynthetic electron transport across plant evolution

Contents Summary 105 I. Introduction 105 II. Diversity of molecular mechanisms for regulation of photosynthetic electron transport 106 III. Role of FLVs in the regulation of photosynthesis in eukaryotes 107 IV. Why were FLVs lost in angiosperms? 108 V. Conclusions 108 Acknowledgements 109 References 109 SUMMARY: Photosynthetic electron transport requires continuous modulation to maintain the balance between light availability and metabolic demands. Multiple mechanisms for the regulation of electron transport have been identified and are unevenly distributed among photosynthetic organisms. Flavodiiron proteins (FLVs) influence photosynthetic electron transport by accepting electrons downstream of photosystem I to reduce oxygen to water. FLV activity has been demonstrated in cyanobacteria, green algae and mosses to be important in avoiding photosystem I overreduction upon changes in light intensity. FLV-encoding sequences were nevertheless lost during evolution by angiosperms, suggesting that these plants increased the efficiency of other mechanisms capable of accepting electrons from photosystem I, making the FLV activity for protection from overreduction superfluous or even detrimental for photosynthetic efficiency.

Genomic comparison of two independent seagrass lineages reveals habitat-driven convergent evolution

Seagrasses are marine angiosperms that live fully submerged in the sea. They evolved from land plant ancestors, with multiple species representing at least three independent return-to-the-sea events. This raises the question of whether these marine angiosperms followed the same adaptation pathway to allow them to live and reproduce under the hostile marine conditions. To compare the basis of marine adaptation between seagrass lineages, we generated genomic data for Halophila ovalis and compared this with recently published genomes for two members of Zosteraceae, as well as genomes of five non-marine plant species (Arabidopsis, Oryza sativa, Phoenix dactylifera, Musa acuminata, and Spirodela polyrhiza). Halophila and Zosteraceae represent two independent seagrass lineages separated by around 30 million years. Genes that were lost or conserved in both lineages were identified. All three species lost genes associated with ethylene and terpenoid biosynthesis, and retained genes related to salinity adaptation, such as those for osmoregulation. In contrast, the loss of the NADH dehydrogenase-like complex is unique to H. ovalis. Through comparison of two independent return-to-the-sea events, this study further describes marine adaptation characteristics common to seagrass families, identifies species-specific gene loss, and provides molecular evidence for convergent evolution in seagrass lineages.

Comparative genomic analysis of Genlisea (corkscrew plants-Lentibulariaceae) chloroplast genomes reveals an increasing loss of the ndh genes

In the carnivorous plant family Lentibulariaceae, all three genome compartments (nuclear, chloroplast, and mitochondria) have some of the highest rates of nucleotide substitutions across angiosperms. While the genera Genlisea and Utricularia have the smallest known flowering plant nuclear genomes, the chloroplast genomes (cpDNA) are mostly structurally conserved except for deletion and/or pseudogenization of the NAD(P)H-dehydrogenase complex (ndh) genes known to be involved in stress conditions of low light or CO2 concentrations. In order to determine how the cpDNA are changing, and to better understand the evolutionary history within the Genlisea genus, we sequenced, assembled and analyzed complete cpDNA from six species (G. aurea, G. filiformis, G. pygmaea, G. repens, G. tuberosa and G. violacea) together with the publicly available G. margaretae cpDNA. In general, the cpDNA structure among the analyzed Genlisea species is highly similar. However, we found that the plastidial ndh genes underwent a progressive process of degradation similar to the other terrestrial Lentibulariaceae cpDNA analyzed to date, but in contrast to the aquatic species. Contrary to current thinking that the terrestrial environment is a more stressful environment and thus requiring the ndh genes, we provide evidence that in the Lentibulariaceae the terrestrial forms have progressive loss while the aquatic forms have the eleven plastidial ndh genes intact. Therefore, the Lentibulariaceae system provides an important opportunity to understand the evolutionary forces that govern the transition to an aquatic environment and may provide insight into how plants manage water stress at a genome scale.

The Complete Plastome Sequences of Seven Species in Gentiana sect. Kudoa (Gentianaceae): Insights Into Plastid Gene Loss and Molecular Evolution

The chloroplast (cp) genome is useful in the study of phylogenomics, molecular dating, and molecular evolution. Gentiana sect. Kudoa is a predominantly alpine flowering plant that is valued for its contributions to medicine, ecology, and horticulture. Previous evolutionary studies showed that the plastid gene loss pattern and intra-sectional phylogenetics in sect. Kudoa are still unclear. In this study, we compared 11 Gentiana plastomes, including 7 newly sequenced plastomes from sect. Kudoa, to represent its three serious: ser. Ornatae, ser. Verticillatae, and ser. Monanthae. The cp genome sizes of the seven species ranged from 137,278 to 147,156 bp. The plastome size variation mainly occurred in the small single-copy and long single-copy regions rather than the inverted repeat regions. Compared with sect. Cruciata, the plastomes in ser. Ornatae and ser. Verticillatae had lost approximately 11 kb of sequences containing 11 ndh genes. Conversely, far fewer losses were observed in ser. Monanthae. The phylogenetic tree revealed that sect. Kudoa was not monophyletic and that ser. Monanthae was more closely related to other sections rather than sect. Kudoa. The molecular dating analysis indicated that ser. Monanthae and sect. Kudoa diverged around 8.23 Ma. In ser. Ornatae and ser. Verticillatae, the divergence occurred at around 0.07-1.78 Ma. The nucleotide diversity analysis indicated that the intergenic regions trnH-psbA, trnK-trnQ, ycf3-trnS and rpl32-trnL constituted divergence hotspots in both sect. Kudoa and Gentiana, and would be useful for future phylogenetic and population genetic studies.

Single-Copy Genes as Molecular Markers for Phylogenomic Studies in Seed Plants

Phylogenetic relationships among seed plant taxa, especially within the gymnosperms, remain contested. In contrast to angiosperms, for which several genomic, transcriptomic and phylogenetic resources are available, there are few, if any, molecular markers that allow broad comparisons among gymnosperm species. With few gymnosperm genomes available, recently obtained transcriptomes in gymnosperms are a great addition to identifying single-copy gene families as molecular markers for phylogenomic analysis in seed plants. Taking advantage of an increasing number of available genomes and transcriptomes, we identified single-copy genes in a broad collection of seed plants and used these to infer phylogenetic relationships between major seed plant taxa. This study aims at extending the current phylogenetic toolkit for seed plants, assessing its ability for resolving seed plant phylogeny, and discussing potential factors affecting phylogenetic reconstruction. In total, we identified 3,072 single-copy genes in 31 gymnosperms and 2,156 single-copy genes in 34 angiosperms. All studied seed plants shared 1,469 single-copy genes, which are generally involved in functions like DNA metabolism, cell cycle, and photosynthesis. A selected set of 106 single-copy genes provided good resolution for the seed plant phylogeny except for gnetophytes. Although some of our analyses support a sister relationship between gnetophytes and other gymnosperms, phylogenetic trees from concatenated alignments without 3rd codon positions and amino acid alignments under the CAT + GTR model, support gnetophytes as a sister group to Pinaceae. Our phylogenomic analyses demonstrate that, in general, single-copy genes can uncover both recent and deep divergences of seed plant phylogeny.

Plastome Evolution in the Sole Hemiparasitic Genus Laurel Dodder (Cassytha) and Insights into the Plastid Phylogenomics of Lauraceae

To date, little is known about the evolution of plastid genomes (plastomes) in Lauraceae. As one of the top five largest families in tropical forests, the Lauraceae contain many species that are important ecologically and economically. Lauraceous species also provide wonderful materials to study the evolutionary trajectory in response to parasitism because they contain both nonparasitic and parasitic species. This study compared the plastomes of nine Lauraceous species, including the sole hemiparasitic and herbaceous genus Cassytha (laurel dodder; here represented by Cassytha filiformis). We found differential contractions of the canonical inverted repeat (IR), resulting in two IR types present in Lauraceae. These two IR types reinforce Cryptocaryeae and Neocinnamomum-Perseeae-Laureae as two separate clades. Our data reveal several traits unique to Cas. filiformis, including loss of IRs, loss or pseudogenization of 11 ndh and rpl23 genes, richness of repeats, and accelerated rates of nucleotide substitutions in protein-coding genes. Although Cas. filiformis is low in chlorophyll content, our analysis based on dN/dS ratios suggests that both its plastid house-keeping and photosynthetic genes are under strong selective constraints. Hence, we propose that short generation time and herbaceous lifestyle rather than reduced photosynthetic ability drive the accelerated rates of nucleotide substitutions in Cas. filiformis.

Complete plastome sequencing of both living species of Circaeasteraceae (Ranunculales) reveals unusual rearrangements and the loss of the ndh gene family

BACKGROUND

Among the 13 families of early-diverging eudicots, only Circaeasteraceae (Ranunculales), which consists of the two monotypic genera Circaeaster and Kingdonia, lacks a published complete plastome sequence. In addition, the phylogenetic position of Circaeasteraceae as sister to Lardizabalaceae has only been weakly or moderately supported in previous studies using smaller data sets. Moreover, previous plastome studies have documented a number of novel structural rearrangements among early-divergent eudicots. Hence it is important to sequence plastomes from Circaeasteraceae to better understand plastome evolution in early-diverging eudicots and to further investigate the phylogenetic position of Circaeasteraceae.

RESULTS

Using an Illumina HiSeq 2000, complete plastomes were sequenced from both living members of Circaeasteraceae: Circaeaster agrestis and Kingdonia uniflora . Plastome structure and gene content were compared between these two plastomes, and with those of other early-diverging eudicot plastomes. Phylogenetic analysis of a 79-gene, 99-taxon data set including exemplars of all families of early-diverging eudicots was conducted to resolve the phylogenetic position of Circaeasteraceae. Both plastomes possess the typical quadripartite structure of land plant plastomes. However, a large ~49 kb inversion and a small ~3.5 kb inversion were found in the large single-copy regions of both plastomes, while Circaeaster possesses a number of other rearrangements, particularly in the Inverted Repeat. In addition, infA was found to be a pseudogene and accD was found to be absent within Circaeaster, whereas all ndh genes, except for ndhE and ndhJ, were found to be either pseudogenized (ΨndhA, ΨndhB, ΨndhD, ΨndhH and ΨndhK) or absent (ndhC, ndhF, ndhI and ndhG) in Kingdonia. Circaeasteraceae was strongly supported as sister to Lardizabalaceae in phylogenetic analyses.

CONCLUSION

The first plastome sequencing of Circaeasteraceae resulted in the discovery of several unusual rearrangements and the loss of ndh genes, and confirms the sister relationship between Circaeasteraceae and Lardizabalaceae. This research provides new insight to characterize plastome structural evolution in early-diverging eudicots and to better understand relationships within Ranunculales .

Concomitant loss of NDH complex-related genes within chloroplast and nuclear genomes in some orchids

The chloroplast NAD(P)H dehydrogenase-like (NDH) complex consists of about 30 subunits from both the nuclear and chloroplast genomes and is ubiquitous across most land plants. In some orchids, such as Phalaenopsis equestris, Dendrobium officinale and Dendrobium catenatum, most of the 11 chloroplast genome-encoded ndh genes (cp-ndh) have been lost. Here we investigated whether functional cp-ndh genes have been completely lost in these orchids or whether they have been transferred and retained in the nuclear genome. Further, we assessed whether both cp-ndh genes and nucleus-encoded NDH-related genes can be lost, resulting in the absence of the NDH complex. Comparative analyses of the genome of Apostasia odorata, an orchid species with a complete complement of cp-ndh genes which represents the sister lineage to all other orchids, and three published orchid genome sequences for P. equestris, D. officinale and D. catenatum, which are all missing cp-ndh genes, indicated that copies of cp-ndh genes are not present in any of these four nuclear genomes. This observation suggests that the NDH complex is not necessary for some plants. Comparative genomic/transcriptomic analyses of currently available plastid genome sequences and nuclear transcriptome data showed that 47 out of 660 photoautotrophic plants and all the heterotrophic plants are missing plastid-encoded cp-ndh genes and exhibit no evidence for maintenance of a functional NDH complex. Our data indicate that the NDH complex can be lost in photoautotrophic plant species. Further, the loss of the NDH complex may increase the probability of transition from a photoautotrophic to a heterotrophic life history.

Contrasting Patterns of Nucleotide Substitution Rates Provide Insight into Dynamic Evolution of Plastid and Mitochondrial Genomes of Geranium

Geraniaceae have emerged as a model system for investigating the causes and consequences of variation in plastid and mitochondrial genomes. Incredible structural variation in plastid genomes (plastomes) and highly accelerated evolutionary rates have been reported in selected lineages and functional groups of genes in both plastomes and mitochondrial genomes (mitogenomes), and these phenomena have been implicated in cytonuclear incompatibility. Previous organelle genome studies have included limited sampling of Geranium, the largest genus in the family with over 400 species. This study reports on rates and patterns of nucleotide substitutions in plastomes and mitogenomes of 17 species of Geranium and representatives of other Geraniaceae. As detected across other angiosperms, substitution rates in the plastome are 3.5 times higher than the mitogenome in most Geranium. However, in the branch leading to Geranium brycei/Geranium incanum mitochondrial genes experienced significantly higher dN and dS than plastid genes, a pattern that has only been detected in one other angiosperm. Furthermore, rate accelerations differ in the two organelle genomes with plastomes having increased dN and mitogenomes with increased dS. In the Geranium phaeum/Geranium reflexum clade, duplicate copies of clpP and rpoA genes that experienced asymmetric rate divergence were detected in the single copy region of the plastome. In the case of rpoA, the branch leading to G. phaeum/G. reflexum experienced positive selection or relaxation of purifying selection. Finally, the evolution of acetyl-CoA carboxylase is unusual in Geraniaceae because it is only the second angiosperm family where both prokaryotic and eukaryotic ACCases functionally coexist in the plastid.

The Complete Plastome Sequences of Four Orchid Species: Insights into the Evolution of the Orchidaceae and the Utility of Plastomic Mutational Hotspots

Orchidaceae (orchids) is the largest family in the monocots, including about 25,000 species in 880 genera and five subfamilies. Many orchids are highly valued for their beautiful and long-lasting flowers. However, the phylogenetic relationships among the five orchid subfamilies remain unresolved. The major dispute centers on whether the three one-stamened subfamilies, Epidendroideae, Orchidoideae, and Vanilloideae, are monophyletic or paraphyletic. Moreover, structural changes in the plastid genome (plastome) and the effective genetic loci at the species-level phylogenetics of orchids have rarely been documented. In this study, we compared 53 orchid plastomes, including four newly sequenced ones, that represent four remote genera: Dendrobium, Goodyera, Paphiopedilum, and Vanilla. These differ from one another not only in their lengths of inverted repeats and small single copy regions but also in their retention of ndh genes. Comparative analyses of the plastomes revealed that the expansion of inverted repeats in Paphiopedilum and Vanilla is associated with a loss of ndh genes. In orchid plastomes, mutational hotspots are genus specific. After having carefully examined the data, we propose that the three loci 5'trnK-rps16, trnS-trnG, and rps16-trnQ might be powerful markers for genera within Epidendroideae, and clpP-psbB and rps16-trnQ might be markers for genera within Cypripedioideae. After analyses of a partitioned dataset, we found that our plastid phylogenomic trees were congruent in a topology where two one-stamened subfamilies (i.e., Epidendroideae and Orchidoideae) were sisters to a multi-stamened subfamily (i.e., Cypripedioideae) rather than to the other one-stamened subfamily (Vanilloideae), suggesting that the living one-stamened orchids are paraphyletic.