Plastid ndh genes in plant evolution

Cyclic electron flow and Photosystem II-less photosynthesis

Oxygenic photosynthesis is characterised by the cooperation of two photo-driven complexes, Photosystem II (PSII) and Photosystem I (PSI), sequentially linked through a series of redox-coupled intermediates. Divergent evolution has resulted in photosystems exhibiting complementary redox potentials, spanning the range necessary to oxidise water and reduce CO2 within a single system. Catalysing nature's most oxidising reaction to extract electrons from water is a highly specialised task that limits PSII's metabolic function. In contrast, potential electron donors in PSI span a range of redox potentials, enabling it to accept electrons from various metabolic processes. This metabolic flexibility of PSI underpins the capacity of photosynthetic organisms to balance energy supply with metabolic demands, which is key for adaptation to environmental changes. Here, we review the phenomenon of 'PSII-less photosynthesis' where PSI functions independently of PSII by operating cyclic electron flow using electrons derived from non-photochemical reactions. PSII-less photosynthesis enables supercharged ATP production and is employed, for example, by cyanobacteria's heterocysts to host nitrogen fixation and by bundle sheath cells of C4 plants to boost CO2 assimilation. We discuss the energetic benefits of this arrangement and the prospects of utilising it to improve the productivity and stress resilience of photosynthetic organisms.

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

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

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

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

MORF9-dependent specific plastid RNA editing inhibits root growth under sugar starvation in Arabidopsis

Multiple organellar RNA editing factor (MORF) complex was shown to be highly associated with C-to-U RNA editing of vascular plant editosome. However, mechanisms by which MORF9-dependent plastid RNA editing controls plant development and responses to environmental alteration remain obscure. In this study, we found that loss of MORF9 function impaired PSII efficiency, NDH activity, and carbohydrate production, rapidly promoted nuclear gene expression including sucrose transporter and sugar/energy responsive genes, and attenuated root growth under sugar starvation conditions. Sugar repletion increased MORF9 and MORF2 expression in wild-type seedlings and reduced RNA editing of matK-706, accD-794, ndhD-383 and ndhF-290 in the morf9 mutant. RNA editing efficiency of ndhD-383 and ndhF-290 sites was diminished in the gin2/morf9 double mutants, and that of matK-706, accD-794, ndhD-383 and ndhF-290 sites were significantly diminished in the snrk1/morf9 double mutants. In contrast, overexpressing HXK1 or SnRK1 promoted RNA editing rate of matK-706, accD-794, ndhD-383 and ndhF-290 in leaves of morf9 mutants, suggesting that HXK1 partially impacts MORF9 mediated ndhD-383 and ndhF-290 editing, while SnRK1 may only affect MORF9-mediated ndhF-290 site editing. Collectively, these findings suggest that sugar and/or its intermediary metabolites impair MORF9-dependent plastid RNA editing resulting in derangements of plant root development.

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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.

Genomic insights into the clonal reproductive Opuntia cochenillifera: mitochondrial and chloroplast genomes of the cochineal cactus for enhanced understanding of structural dynamics and evolutionary implications

Background

The cochineal cactus (Opuntia cochenillifera), notable for its substantial agricultural and industrial applications, predominantly undergoes clonal reproduction, which presents significant challenges in breeding and germplasm innovation. Recent developments in mitochondrial genome engineering offer promising avenues for introducing heritable mutations, potentially facilitating selective sexual reproduction through the creation of cytoplasmic male sterile genotypes. However, the lack of comprehensive mitochondrial genome information for Opuntia species hinders these efforts. Here, we intended to sequence and characterize its mitochondrial genome to maximize the potential of its genomes for evolutionary studies, molecular breeding, and molecular marker developments.

Results

We sequenced the total DNA of the O. cochenillifera using DNBSEQ and Nanopore platforms. The mitochondrial genome was then assembled using a hybrid assembly strategy using Unicycler software. We found that the mitochondrial genome of O. cochenillifera has a length of 1,156,235 bp, a GC content of 43.06%, and contains 54 unique protein-coding genes and 346 simple repeats. Comparative genomic analysis revealed 48 homologous fragments shared between mitochondrial and chloroplast genomes, with a total length of 47,935 bp. Additionally, the comparison of mitochondrial genomes from four Cactaceae species highlighted their dynamic nature and frequent mitogenomic reorganizations.

Conclusion

Our study provides a new perspective on the evolution of the organelle genome and its potential application in genetic breeding. These findings offer valuable insights into the mitochondrial genetics of Cactaceae, potentially facilitating future research and breeding programs aimed at enhancing the genetic diversity and adaptability of O. cochenillifera by leveraging its unique mitochondrial genome characteristics.

The Complete Chloroplast Genomes of Bulbophyllum (Orchidaceae) Species: Insight into Genome Structure Divergence and Phylogenetic Analysis

Bulbophyllum is one of the largest genera and presents some of the most intricate taxonomic problems in the family Orchidaceae, including species of ornamental and medical importance. The lack of knowledge regarding the characterization of Bulbophyllum chloroplast (cp) genomes has imposed current limitations on our study. Here, we report the complete cp genomes of seven Bulbophyllum species, including B. ambrosia, B. crassipes, B. farreri, B. hamatum, B. shanicum, B. triste, and B. violaceolabellum, and compared with related taxa to provide a better understanding of their genomic information on taxonomy and phylogeny. A total of 28 Bulbophyllum cp genomes exhibit typical quadripartite structures with lengths ranging from 145,092 bp to 165,812 bp and a GC content of 36.60% to 38.04%. Each genome contained 125-132 genes, encompassing 74-86 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The genome arrangements, gene contents, and length were similar, with differences observed in ndh gene composition. It is worth noting that there were exogenous fragment insertions in the IR regions of B. crassipes. A total of 18-49 long repeats and 38-80 simple sequence repeats (SSRs) were detected and the single nucleotide (A/T) was dominant in Bulbophyllum cp genomes, with an obvious A/T preference. An analysis of relative synonymous codon usage (RSCU) revealed that leucine (Leu) was the most frequently used codon, while cysteine (Cys) was the least used. Six highly variable regions (rpl32-trnLUAG > trnTUGU-trnLUAA > trnFGAA-ndhJ > rps15-ycf1 > rbcL-accD > psbI-trnSGCU) and five coding sequences (ycf1 > rps12 > matK > psbK > rps15) were identified as potential DNA markers based on nucleotide diversity. Additionally, 31,641 molecular diagnostic characters (MDCs) were identified in complete cp genomes. A phylogenetic analysis based on the complete cp genome sequences and 68 protein-coding genes strongly supported that 28 Bulbophyllum species can be divided into four branches, sects. Brachyantha, Cirrhopetalum, and Leopardinae, defined by morphology, were non-monophyly. Our results enriched the genetic resources of Bulbophyllum, providing valuable information to illustrate the complicated taxonomy, phylogeny, and evolution process of the genus.

Complete Chloroplast Genomes of Four Oaks from the Section Cyclobalanopsis Improve the Phylogenetic Analysis and Understanding of Evolutionary Processes in the Genus Quercus

Quercus is a valuable genus ecologically, economically, and culturally. They are keystone species in many ecosystems. Species delimitation and phylogenetic studies of this genus are difficult owing to frequent hybridization. With an increasing number of genetic resources, we will gain a deeper understanding of this genus. In the present study, we collected four Quercus section Cyclobalanopsis species (Q. poilanei, Q. helferiana, Q. camusiae, and Q. semiserrata) distributed in Southeast Asia and sequenced their complete genomes. Following analysis, we compared the results with those of other species in the genus Quercus. These four chloroplast genomes ranged from 160,784 bp (Q. poilanei) to 161,632 bp (Q. camusiae) in length, with an overall guanine and cytosine (GC) content of 36.9%. Their chloroplast genomic organization and order, as well as their GC content, were similar to those of other Quercus species. We identified seven regions with relatively high variability (rps16, ndhk, accD, ycf1, psbZ-trnG-GCC, rbcL-accD, and rpl32-trnL-UAG) which could potentially serve as plastid markers for further taxonomic and phylogenetic studies within Quercus. Our phylogenetic tree supported the idea that the genus Quercus forms two well-differentiated lineages (corresponding to the subgenera Quercus and Cerris). Of the three sections in the subgenus Cerris, the section Ilex was split into two clusters, each nested in the other two sections. Moreover, Q. camusiae and Q. semiserrata detected in this study diverged first in the section Cyclobalanopsis and mixed with Q. engleriana in the section Ilex. In particular, 11 protein coding genes (atpF, ndhA, ndhD, ndhF, ndhK, petB, petD, rbcL, rpl22, ycf1, and ycf3) were subjected to positive selection pressure. Overall, this study enriches the chloroplast genome resources of Quercus, which will facilitate further analyses of phylogenetic relationships in this ecologically important tree genus.

CLAW: An automated Snakemake workflow for the assembly of chloroplast genomes from long-read data

Chloroplasts are photosynthetic organelles in algal and plant cells that contain their own genome. Chloroplast genomes are commonly used in evolutionary studies and taxonomic identification and are increasingly becoming a target for crop improvement studies. As DNA sequencing becomes more affordable, researchers are collecting vast swathes of high-quality whole-genome sequence data from laboratory and field settings alike. Whole tissue read libraries sequenced with the primary goal of understanding the nuclear genome will inadvertently contain many reads derived from the chloroplast genome. These whole-genome, whole-tissue read libraries can additionally be used to assemble chloroplast genomes with little to no extra cost. While several tools exist that make use of short-read second generation and third-generation long-read sequencing data for chloroplast genome assembly, these tools may have complex installation steps, inadequate error reporting, poor expandability, and/or lack scalability. Here, we present CLAW (Chloroplast Long-read Assembly Workflow), an easy to install, customise, and use Snakemake tool to assemble chloroplast genomes from chloroplast long-reads found in whole-genome read libraries (https://github.com/aaronphillips7493/CLAW). Using 19 publicly available reference chloroplast genome assemblies and long-read libraries from algal, monocot and eudicot species, we show that CLAW can rapidly produce chloroplast genome assemblies with high similarity to the reference assemblies. CLAW was designed such that users have complete control over parameterisation, allowing individuals to optimise CLAW to their specific use cases. We expect that CLAW will provide researchers (with varying levels of bioinformatics expertise) with an additional resource useful for contributing to the growing number of publicly available chloroplast genome assemblies.

Comparative analyses of complete chloroplast genomes reveal interspecific difference and intraspecific variation of Tripterygium genus

The genus Tripterygium was of great medicinal value and attracted much attention on the taxonomic study using morphological and molecular methods. In this study, we assembled 12 chloroplast genomes of Tripterygium to reveal interspecific difference and intraspecific variation. The sequence length (156,692-157,061 bp) and structure of Tripterygium were conserved. Comparative analyses presented abundant variable regions for further study. Meanwhile, we determined the ndhB gene under positive selection through adaptive evolution analysis. And the phylogenetic analyses based on 15 chloroplast genomes supported the monophyly of Tripterygium hypoglaucum and the potential sister relationship between Tripterygium wilfordii and Tripterygium regelii. Molecular dating analysis indicated that the divergence time within Tripterygium was approximately 5.99 Ma (95% HPD = 3.11-8.68 Ma). The results in our study provided new insights into the taxonomy, evolution process, and phylogenetic construction of Tripterygium using complete plastid genomes.

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.

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.

A Machine Learning Framework Identifies Plastid-Encoded Proteins Harboring C3 and C4 Distinguishing Sequence Information

C4 photosynthesis is known to have at least 61 independent origins across plant lineages making it one of the most notable examples of convergent evolution. Of the >60 independent origins, a predicted 22-24 origins, encompassing greater than 50% of all known C4 species, exist within the Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) clade of the Poaceae family. This clade is therefore primed with species ideal for the study of genomic changes associated with the acquisition of the C4 photosynthetic trait. In this study, we take advantage of the growing availability of sequenced plastid genomes and employ a machine learning (ML) approach to screen for plastid genes harboring C3 and C4 distinguishing information in PACMAD species. We demonstrate that certain plastid-encoded protein sequences possess distinguishing and informative sequence information that allows them to train accurate ML C3/C4 classification models. Our RbcL-trained model, for example, informs a C3/C4 classifier with greater than 99% accuracy. Accurate prediction of photosynthetic type from individual sequences suggests biologically relevant, and potentially differing roles of these sequence products in C3 versus C4 metabolism. With this ML framework, we have identified several key sequences and sites that are most predictive of C3/C4 status, including RbcL, subunits of the NAD(P)H dehydrogenase complex, and specific residues within, further highlighting their potential significance in the evolution and/or maintenance of C4 photosynthetic machinery. This general approach can be applied to uncover intricate associations between other similar genotype-phenotype relationships.

Complete Chloroplast Genome Determination of Ranunculus sceleratus from Republic of Korea (Ranunculaceae) and Comparative Chloroplast Genomes of the Members of the Ranunculus Genus

Ranunculus sceleratus (family: Ranunculaceae) is a medicinally and economically important plant; however, gaps in taxonomic and species identification limit its practical applicability. This study aimed to sequence the chloroplast genome of R. sceleratus from Republic of Korea. Chloroplast sequences were compared and analyzed among Ranunculus species. The chloroplast genome was assembled from Illumina HiSeq 2500 sequencing raw data. The genome was 156,329 bp and had a typical quadripartite structure comprising a small single-copy region, a large single-copy region, and two inverted repeats. Fifty-three simple sequence repeats were identified in the four quadrant structural regions. The region between the ndhC and trnV-UAC genes could be useful as a genetic marker to distinguish between R. sceleratus populations from Republic of Korea and China. The Ranunculus species formed a single lineage. To differentiate between Ranunculus species, we identified 16 hotspot regions and confirmed their potential using specific barcodes based on phylogenetic tree and BLAST-based analyses. The ndhE, ndhF, rpl23, atpF, rps4, and rpoA genes had a high posterior probability of codon sites in positive selection, while the amino acid site varied between Ranunculus species and other genera. Comparison of the Ranunculus genomes provides useful information regarding species identification and evolution that could guide future phylogenetic analyses.

Backbone phylogeny and adaptive evolution of Pleurospermum s. l.: New insights from phylogenomic analyses of complete plastome data

Pleurospermum is a taxonomically challenging taxon of Apiaceae, as its circumscription and composition remain controversial for morphological similarities with several related genera, leading to a dispute between Pleurospermum in the broad sense and strict sense. While evidence from previous molecular studies recognized plural branching lineages within the Pleurospermum s. l., it did not support the latest delimitation of Pleurospermum s. str. by only two closely related northern species. So far, no proper delimitation for Pleurospermum has come up, and many of the plural taxa in Pleurospermum s. l. remain unresolved, which may be due to poor phylogenetic resolution yielded barely from ITS sequences. Herein, we newly assembled 40 complete plastomes from 36 species of Pleurospermum s. l. and related genera, 34 of which were first reported and generated a well-resolved backbone phylogeny in a framework of the subfamily Apioideae. From the phylogeny with greatly improved resolution, a total of six well-supported monophyletic lineages within Pleurospermum s. l. were recognized falling in different major clades of Apioideae. Combining morphological characteristics with phylogenetic inference, we suggested to re-delimit the Pleurospermum s. str. by introducing nine species mainly from the Himalayan regions and proposed its boundary features; the remaining species were suggested to be excluded from Pleurospermum to incorporate into their more related taxa being revealed. On this basis, the plastome comparison revealed not only the high conservatism but also the mild differences among lineages in plastome structure and gene evolution. Overall, our study provided a backbone phylogeny essential for further studies of the taxonomically difficult taxa within Pleurospermum s. l.

Evolution of 101 Apocynaceae plastomes and phylogenetic implications

Apocynaceae are one of the ten species-richest angiosperm families. However, the backbone phylogeny of the family is yet less well supported, and the evolution of plastome structure has not been thoroughly studied for the whole family. Herein, a total of 101 complete plastomes including 35 newly sequenced, 24 reassembled from public raw data and the rest from the NCBI GenBank database, representing 26 of 27 tribes of Apocynaceae, were used for comparative plastome analysis. Phylogenetic analyses were conducted using a combined plastid data matrix of 77 protein-coding genes from 162 taxa, encompassing all tribes and 41 of 49 subtribes of Apocynaceae. Plastome lengths ranged from 150,897 bp in Apocynum venetum to 178,616 bp in Hoya exilis. Six types of boundaries between the inverted repeat (IR) regions and single copy (SC) regions were identified. Different sizes of IR expansion were found in three lineages, including Alyxieae, Ceropegieae and Marsdenieae, suggesting multiple expansion events of the IRs over the SC regions in Apocynaceae. The IR regions of Marsdenieae evolved in two ways: expansion towards the large single copy (LSC) region in Lygisma + Stephanotis + Ruehssia + Gymnema (Cosmopolitan clade), and expansion towards both LSC and small single copy (SSC) region in Dischidia-Hoya alliance and Marsdenia (Asia-Pacific clade). Six coding genes and five non-coding regions were identified as highly variable, including accD, ccsA-ndhD, clpP, matK, ndhF, ndhG-ndhI, trnG(GCC)-trnfM(CAU), trnH(GUG)-psbA, trnY(GUA)-trnE(UUC), ycf1, and ycf2. Maximum likelihood and Bayesian phylogenetic analyses resulted in nearly identical tree topologies and produced a well-resolved backbone comprising 15 consecutive dichotomies that subdivided Apocynaceae into 15 clades. The subfamily Periplocoideae were embedded in the Apocynoid grade and were sister to the Echiteae-Odontadenieae-Mesechiteae clade with high support values. Three tribes (Melodineae, Vinceae, and Willughbeieae), the subtribe Amphineuriinae, and four genera (Beaumontia, Ceropegia, Hoya, and Stephanotis) were not resolved as monophyletic. Our work sheds light on the backbone phylogenetic relationships in the family Apocynaceae and offers insights into the evolution of Apocynaceae plastomes using the most densely sampled plastome dataset to date.

A comparison of 25 complete chloroplast genomes between sister mangrove species Kandelia obovata and Kandelia candel geographically separated by the South China Sea

In 2003, Kandelia obovata was identified as a new mangrove species differentiated from Kandelia candel. However, little is known about their chloroplast (cp) genome differences and their possible ecological significance. In this study, 25 whole cp genomes, with seven samples of K. candel from Malaysia, Thailand, and Bangladesh and 18 samples of K. obovata from China, were sequenced for comparison. The cp genomes of both species encoded 128 genes, namely 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes, but the cp genome size of K. obovata was ~2 kb larger than that of K. candle due to the presence of more and longer repeat sequences. Of these, tandem repeats and simple sequence repeats exhibited great differences. Principal component analysis based on indels, and phylogenetic tree analyses constructed with homologous protein genes from the single-copy genes, as well as 38 homologous pair genes among 13 mangrove species, gave strong support to the separation of the two species within the Kandelia genus. Homologous genes ndhD and atpA showed intraspecific consistency and interspecific differences. Molecular dynamics simulations of their corresponding proteins, NAD(P)H dehydrogenase chain 4 (NDH-D) and ATP synthase subunit alpha (ATP-A), predicted them to be significantly different in the functions of photosynthetic electron transport and ATP generation in the two species. These results suggest that the energy requirement was a pivotal factor in their adaptation to differential environments geographically separated by the South China Sea. Our results also provide clues for future research on their physiological and molecular adaptation mechanisms to light and temperature.

Chloroplast phylogenomics and the taxonomy of Saxifraga section Ciliatae (Saxifragaceae)

Comprising ca. 200 species, Saxifraga sect. Ciliatae is the most species-rich section of Saxifraga s.str., whose center of diversity is in the Tibeto-Himalayan region. The infra-sectional classification of sect. Ciliatae is still in debate due to the high level of species richness, as well as remarkable variations of habitat, morphology, physiology and life cycles. Subdivisions of this section proposed in various taxonomic systems have not been adequately tested in previous phylogenetic studies, partly due to low taxonomic sampling density, but also to the use of few DNA markers. In order to achieve a more robust infra-sectional classification of sect. Ciliatae, complete chloroplast genomes of 94 taxa from this section were analyzed, of which 93 were newly sequenced, assembled and annotated. The length of the 94 plastomes of sect. Ciliatae taxa range from 143,479 to 159,938 bp, encoding 75 to 79 unique protein-coding genes (PCGs). Analyses of the 94 plastomes revealed high conservation in structural organization, gene arrangement, and gene content. Gene loss and changes of IR boundaries were detected but in extremely low frequency. The molecular phylogenetic tree from concatenated PCGs and complete chloroplast genome sequences exhibited high resolution and support values and confirms that sect. Ciliatae is monophyletic. Three well-supported clades were revealed within the section that agree relatively well with the subsectional taxonomy of Gornall (1987), but some minor modifications should be made. Firstly, the monotypic subsection Cinerascentes should be abandoned and its constituent species, S. cinerascens, assigned to subsect. Gemmiparae. Secondly, subsections Rosulares and Serpyllifoliae should be merged and become subsect. Rosulares. Section Ciliatae thus comprises: subsect. Hirculoideae Engl. & Irmsch.; subsect. Rosulares Gornall; subsect. Gemmiparae Engl. & Irmsch.; subsect. Flagellares (C. B. Clarke) Engl. & Irmsch. and subsect. Hemisphaericae (Engl. & Irmsch.) Gornall.

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.

The ecological adaptation of the unparalleled plastome character evolution in slipper orchids

Plastomes may have undergone adaptive evolution in the process of plant adaptation to diverse environments, whereby species may differ in plastome characters. Cypripedioideae successfully colonized distinct environments and could be an ideal group for studying the interspecific variation and adaptive evolution of plastomes. Comparative study of plastomes, ancestral state reconstruction, phylogenetic-based analysis, ecological niche modelling, and selective pressure analysis were conducted to reveal the evolutionary patterns of plastomes in Cypripedioideae and their relationship with environmental factors. The plastomes of the three evolved genera had reduced plastome size, increased GC content, and compacted gene content compared to the basal group. Variations in plastome size and GC content are proved to have clear relationships with climate regions. Furthermore, ecological niche modelling revealed that temperature and water factors are important climatic factors contributing to the distributional difference which is directly correlated with the climate regions. The temperature-sensitive genes ndh genes, infA, and rpl20 were found to be either lost/pseudogenized or under positive selection in the evolved groups. Unparalleled plastome character variations were discovered in slipper orchids. Our study indicates that variations in plastome characters have adaptive consequences and that temperature and water factors are important climatic factors that affect plastome evolution. This research highlights the expectation that plants can facilitate adaptation to different environmental conditions with the changes in plastome and has added critical insight for understanding the process of plastome evolution in plants.

Comparative chloroplast genomics provides insights into the genealogical relationships of endangered Tetraena mongolica and the chloroplast genome evolution of related Zygophyllaceae species

A comprehensive understanding of genetic background for rare species will provide an important theoretical basis for the future species management, monitoring and conservation. Tetraena mongolica is restrictedly distributed in the western Ordos plateau of China and has been listed as a national protected plant. We generated 13 chloroplast (cp) genomes of T. mongolica (size range of 106,062-106,230 bp) and conducted a series of comparative analyses of six Zygophyllaceae cp genomes. T. mongolica cp genome exhibited a quadripartite structure with drastically reduced inverted repeats (IRs, 4,315 bp) and undergone the loss of a suit of ndh genes and a copy of rRNAs. Furthermore, all the T. mongolica populations were divided into two genetic groups based on complete cp phylogenomics. In addition, notably variable genome size, gene order and structural changes had been observed among the six Zygophyllaceae cp genomes. Overall, our findings provide insights into the cp genome evolution mode and intraspecific relationships of T. mongolica, and provide a molecular basis for scientific conservation of this endangered plant.

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.

An ancient function of PGR5 in iron delivery?

In all phototrophic organisms, the photosynthetic apparatus must be protected from light-induced damage. One important mechanism that mitigates photodamage in plants is antimycin A (AA)-sensitive cyclic electron flow (CEF), the evolution of which remains largely obscure. Here we show that proton gradient regulation 5 (PGR5), a key protein involved in AA-sensitive CEF, displays intriguing commonalities - including sequence and structural features - with a group of ferritin-like proteins. We therefore propose that PGR5 may originally have been involved in prokaryotic iron mobilization and delivery, which facilitated a primordial type of CEF as a side effect. The abandonment of the bacterioferritin system during the transformation of cyanobacterial endosymbionts into chloroplasts might have allowed PGR5 to functionally specialize in CEF.

Decoding RNA Editing Sites Through Transcriptome Analysis in Rice Under Alkaline Stress

Ribonucleic acid editing (RE) is a post-transcriptional process that altered the genetics of RNA which provide the extra level of gene expression through insertion, deletions, and substitutions. In animals, it converts nucleotide residues C-U. Similarly in plants, the role of RNA editing sites (RES) in rice under alkaline stress is not fully studied. Rice is a staple food for most of the world population. Alkaline stress cause reduction in yield. Here, we explored the effect of alkaline stress on RES in the whole mRNA from rice chloroplast and mitochondria. Ribonucleic acid editing sites in both genomes (3336 RESs) including chloroplast (345 RESs) and mitochondria (2991 RESs) with average RES efficiency ∼55% were predicted. Our findings showed that majority of editing events found in non-synonymous codon changes and change trend in amino acids was hydrophobic. Four types of RNA editing A-G (A-I), C-T (C-U), G-A, and T-C were identified in treated and untreated samples. Overall, RNA editing efficiency was increased in the treated samples. Analysis of Gene Ontology revealed that mapped genes were engaged in many biological functions and molecular processes. We also checked the expression of pentatricopeptide repeat (PPR), organelle zinc-finger (OZI), and multiple organellar RNA editing factors/RNA editing factor interacting proteins genes in control and treatment, results revealed upregulation of PPR and OZ1 genes in treated samples. This induction showed the role of these genes in RNA editing. The current findings report that RNA editing increased under alkaline stress which may contribute in adaptation for rice by changing amino acids in edited genes (88 genes). These findings will provide basis for identification of RES in other crops and also will be useful in alkaline tolerance development in rice.

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.

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.

Molecular evolution and phylogenetic relationships of Ligusticum (Apiaceae) inferred from the whole plastome sequences

Background

The genus Ligusticum belongs to Apiaceae, and its taxonomy has long been a major difficulty. A robust phylogenetic tree is the basis of accurate taxonomic classification of Ligusticum. We herein used 26 (including 14 newly sequenced) plastome-scale data to generate reliable phylogenetic trees to explore the phylogenetic relationships of Chinese Ligusticum.

Results

We found that these plastid genomes exhibited diverse plastome characteristics across all four currently identified clades in China, while the plastid protein-coding genes were conserved. The phylogenetic analyses by the concatenation and coalescent methods obtained a more robust molecular phylogeny than prior studies and showed the non-monophyly of Chinese Ligusticum. In the concatenation-based phylogeny analyses, the two datasets yielded slightly different topologies that may be primarily due to the discrepancy in the number of variable sites.

Conclusions

Our plastid phylogenomics analyses emphasized that the current circumscription of the Chinese Ligusticum should be reduced, and the taxonomy of Ligusticum urgently needs revision. Wider taxon sampling including the related species of Ligusticum will be necessary to explore the phylogenetic relationships of this genus. Overall, our study provided new insights into the taxonomic classification of Ligusticum and would serve as a framework for future studies on taxonomy and delimitation of Ligusticum from the perspective of the plastid genome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02010-z.

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.

Tightly Constrained Genome Reduction and Relaxation of Purifying Selection during Secondary Plastid Endosymbiosis

Endosymbiosis, the establishment of a former free-living prokaryotic or eukaryotic cell as an organelle inside a host cell, can dramatically alter the genomic architecture of the endosymbiont. Plastids or chloroplasts, the light-harvesting organelle of photosynthetic eukaryotes, are excellent models to study this phenomenon because plastid origin has occurred multiple times in evolution. Here, we investigate the genomic signature of molecular processes acting through secondary plastid endosymbiosis-the origination of a new plastid from a free-living eukaryotic alga. We used phylogenetic comparative methods to study gene loss and changes in selective regimes on plastid genomes, focusing on green algae that have given rise to three independent lineages with secondary plastids (euglenophytes, chlorarachniophytes, and Lepidodinium). Our results show an overall increase in gene loss associated with secondary endosymbiosis, but this loss is tightly constrained by the retention of genes essential for plastid function. The data show that secondary plastids have experienced temporary relaxation of purifying selection during secondary endosymbiosis. However, this process is tightly constrained, with selection relaxed only relative to the background in primary plastids. Purifying selection remains strong in absolute terms even during the endosymbiosis events. Selection intensity rebounds to pre-endosymbiosis levels following endosymbiosis events, demonstrating the changes in selection efficiency during different origin phases of secondary plastids. Independent endosymbiosis events in the euglenophytes, chlorarachniophytes, and Lepidodinium differ in their degree of relaxation of selection, highlighting the different evolutionary contexts of these events. This study reveals the selection-drift interplay during secondary endosymbiosis and evolutionary parallels during organellogenesis.

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.

Extensive reorganization of the chloroplast genome of Corydalis platycarpa: A comparative analysis of their organization and evolution with other Corydalis plastomes

Introduction

The chloroplast (cp) is an autonomous plant organelle with an individual genome that encodes essential cellular functions. The genome architecture and gene content of the cp is highly conserved in angiosperms. The plastome of Corydalis belongs to the Papaveraceae family, and the genome is comprised of unusual rearrangements and gene content. Thus far, no extensive comparative studies have been carried out to understand the evolution of Corydalis chloroplast genomes.

Methods

Therefore, the Corydalis platycarpa cp genome was sequenced, and wide-scale comparative studies were conducted using publicly available twenty Corydalis plastomes.

Results

Comparative analyses showed that an extensive genome rearrangement and IR expansion occurred, and these events evolved independently in the Corydalis species. By contrast, the plastomes of its closely related subfamily Papaveroideae and other Ranunculales taxa are highly conserved. On the other hand, the synapomorphy characteristics of both accD and the ndh gene loss events happened in the common ancestor of the Corydalis and sub-clade of the Corydalis lineage, respectively. The Corydalis-sub clade species (ndh lost) are distributed predominantly in the Qinghai-Tibetan plateau (QTP) region. The phylogenetic analysis and divergence time estimation were also employed for the Corydalis species.

Discussion

The divergence time of the ndh gene in the Corydalis sub-clade species (44.31 - 15.71 mya) coincides very well with the uplift of the Qinghai-Tibet Plateau in Oligocene and Miocene periods, and maybe during this period, it has probably triggered the radiation of the Corydalis species.

Conclusion

To the best of the authors' knowledge, this is the first large-scale comparative study of Corydalis plastomes and their evolution. The present study may provide insights into the plastome architecture and the molecular evolution of Corydalis species.

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.

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.

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.

Occurrence, Evolution and Specificities of Iron-Sulfur Proteins and Maturation Factors in Chloroplasts from Algae

Iron-containing proteins, including iron-sulfur (Fe-S) proteins, are essential for numerous electron transfer and metabolic reactions. They are present in most subcellular compartments. In plastids, in addition to sustaining the linear and cyclic photosynthetic electron transfer chains, Fe-S proteins participate in carbon, nitrogen, and sulfur assimilation, tetrapyrrole and isoprenoid metabolism, and lipoic acid and thiamine synthesis. The synthesis of Fe-S clusters, their trafficking, and their insertion into chloroplastic proteins necessitate the so-called sulfur mobilization (SUF) protein machinery. In the first part, we describe the molecular mechanisms that allow Fe-S cluster synthesis and insertion into acceptor proteins by the SUF machinery and analyze the occurrence of the SUF components in microalgae, focusing in particular on the green alga Chlamydomonas reinhardtii. In the second part, we describe chloroplastic Fe-S protein-dependent pathways that are specific to Chlamydomonas or for which Chlamydomonas presents specificities compared to terrestrial plants, putting notable emphasis on the contribution of Fe-S proteins to chlorophyll synthesis in the dark and to the fermentative metabolism. The occurrence and evolutionary conservation of these enzymes and pathways have been analyzed in all supergroups of microalgae performing oxygenic photosynthesis.

Comparative Chloroplast Genomics of Corydalis Species (Papaveraceae): Evolutionary Perspectives on Their Unusual Large Scale Rearrangements

The chloroplast genome (plastome) of angiosperms (particularly photosynthetic members) is generally highly conserved, although structural rearrangements have been reported in a few lineages. In this study, we revealed Corydalis to be another unusual lineage with extensive large-scale plastome rearrangements. In the four newly sequenced Corydalis plastomes that represent all the three subgenera of Corydalis, we detected (1) two independent relocations of the same five genes (trnV-UAC-rbcL) from the typically posterior part of the large single-copy (LSC) region to the front, downstream of either the atpH gene in Corydalis saxicola or the trnK-UUU gene in both Corydalis davidii and Corydalis hsiaowutaishanensis; (2) relocation of the rps16 gene from the LSC region to the inverted repeat (IR) region in Corydalis adunca; (3) uniform inversion of an 11-14 kb segment (ndhB-trnR-ACG) in the IR region of all the four Corydalis species (the same below); (4) expansions (>10 kb) of IR into the small single-copy (SSC) region and corresponding contractions of SSC region; and (5) extensive pseudogenizations or losses of 13 genes (accD, clpP, and 11 ndh genes). In addition, we also found that the four Corydalis plastomes exhibited elevated GC content in both gene and intergenic regions and high number of dispersed repeats. Phylogenomic analyses generated a well-supported topology that was consistent with the result of previous studies based on a few DNA markers but contradicted with the morphological character-based taxonomy to some extent. This study provided insights into the evolution of plastomes throughout the three Corydalis subgenera and will be of value for further study on taxonomy, phylogeny, and evolution of Corydalis.

Plastomes of eight Ligusticum species: characterization, genome evolution, and phylogenetic relationships

BACKGROUND

The genus Ligusticum consists of approximately 60 species distributed in the Northern Hemisphere. It is one of the most taxonomically difficult taxa within Apiaceae, largely due to the varied morphological characteristics. To investigate the plastome evolution and phylogenetic relationships of Ligusticum, we determined the complete plastome sequences of eight Ligusticum species using a de novo assembly approach.

RESULTS

Through a comprehensive comparative analysis, we found that the eight plastomes were similar in terms of repeat sequence, SSR, codon usage, and RNA editing site. However, compared with the other seven species, L. delavayi exhibited striking differences in genome size, gene number, IR/SC borders, and sequence identity. Most of the genes remained under the purifying selection, whereas four genes showed relaxed selection, namely ccsA, rpoA, ycf1, and ycf2. Non-monophyly of Ligusticum species was inferred from the plastomes and internal transcribed spacer (ITS) sequences phylogenetic analyses.

CONCLUSION

The plastome tree and ITS tree produced incongruent tree topologies, which may be attributed to the hybridization and incomplete lineage sorting. Our study highlighted the advantage of plastome with mass informative sites in resolving phylogenetic relationships. Moreover, combined with the previous studies, we considered that the current taxonomy system of Ligusticum needs to be improved and revised. In summary, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of Ligusticum species.

Cultivated Tomato (Solanum lycopersicum L.) Suffered a Severe Cytoplasmic Bottleneck during Domestication: Implications from Chloroplast Genomes

In various crops, genetic bottlenecks occurring through domestication can limit crop resilience to biotic and abiotic stresses. In the present study, we investigated nucleotide diversity in tomato chloroplast genome through sequencing seven plastomes of cultivated accessions from the Campania region (Southern Italy) and two wild species among the closest (Solanum pimpinellifolium) and most distantly related (S. neorickii) species to cultivated tomatoes. Comparative analyses among the chloroplast genomes sequenced in this work and those available in GenBank allowed evaluating the variability of plastomes and defining phylogenetic relationships. A dramatic reduction in genetic diversity was detected in cultivated tomatoes, nonetheless, a few de novo mutations, which still differentiated the cultivated tomatoes from the closest wild relative S. pimpinellifolium, were detected and are potentially utilizable as diagnostic markers. Phylogenetic analyses confirmed that S. pimpinellifolium is the closest ancestor of all cultivated tomatoes. Local accessions all clustered together and were strictly related with other cultivated tomatoes (S. lycopersicum group). Noteworthy, S. lycopersicum var. cerasiforme resulted in a mixture of both cultivated and wild tomato genotypes since one of the two analyzed accessions clustered with cultivated tomato, whereas the other with S. pimpinellifolium. Overall, our results revealed a very reduced cytoplasmic variability in cultivated tomatoes and suggest the occurrence of a cytoplasmic bottleneck during their domestication.

A systematic comparison of chloroplast genome assembly tools

BACKGROUND

Chloroplasts are intracellular organelles that enable plants to conduct photosynthesis. They arose through the symbiotic integration of a prokaryotic cell into an eukaryotic host cell and still contain their own genomes with distinct genomic information. Plastid genomes accommodate essential genes and are regularly utilized in biotechnology or phylogenetics. Different assemblers that are able to assess the plastid genome have been developed. These assemblers often use data of whole genome sequencing experiments, which usually contain reads from the complete chloroplast genome.

RESULTS

The performance of different assembly tools has never been systematically compared. Here, we present a benchmark of seven chloroplast assembly tools, capable of succeeding in more than 60% of known real data sets. Our results show significant differences between the tested assemblers in terms of generating whole chloroplast genome sequences and computational requirements. The examination of 105 data sets from species with unknown plastid genomes leads to the assembly of 20 novel chloroplast genomes.

CONCLUSIONS

We create docker images for each tested tool that are freely available for the scientific community and ensure reproducibility of the analyses. These containers allow the analysis and screening of data sets for chloroplast genomes using standard computational infrastructure. Thus, large scale screening for chloroplasts within genomic sequencing data is feasible.

A systematic comparison of chloroplast genome assembly tools

BACKGROUND

Chloroplasts are intracellular organelles that enable plants to conduct photosynthesis. They arose through the symbiotic integration of a prokaryotic cell into an eukaryotic host cell and still contain their own genomes with distinct genomic information. Plastid genomes accommodate essential genes and are regularly utilized in biotechnology or phylogenetics. Different assemblers that are able to assess the plastid genome have been developed. These assemblers often use data of whole genome sequencing experiments, which usually contain reads from the complete chloroplast genome.

RESULTS

The performance of different assembly tools has never been systematically compared. Here, we present a benchmark of seven chloroplast assembly tools, capable of succeeding in more than 60% of known real data sets. Our results show significant differences between the tested assemblers in terms of generating whole chloroplast genome sequences and computational requirements. The examination of 105 data sets from species with unknown plastid genomes leads to the assembly of 20 novel chloroplast genomes.

CONCLUSIONS

We create docker images for each tested tool that are freely available for the scientific community and ensure reproducibility of the analyses. These containers allow the analysis and screening of data sets for chloroplast genomes using standard computational infrastructure. Thus, large scale screening for chloroplasts within genomic sequencing data is feasible.

Gene Loss and Evolution of the Plastome

Chloroplasts are unique organelles within the plant cells and are responsible for sustaining life forms on the earth due to their ability to conduct photosynthesis. Multiple functional genes within the chloroplast are responsible for a variety of metabolic processes that occur in the chloroplast. Considering its fundamental role in sustaining life on the earth, it is important to identify the level of diversity present in the chloroplast genome, what genes and genomic content have been lost, what genes have been transferred to the nuclear genome, duplication events, and the overall origin and evolution of the chloroplast genome. Our analysis of 2511 chloroplast genomes indicated that the genome size and number of coding DNA sequences (CDS) in the chloroplasts genome of algae are higher relative to other lineages. Approximately 10.31% of the examined species have lost the inverted repeats (IR) in the chloroplast genome that span across all the lineages. Genome-wide analyses revealed the loss of the Rbcl gene in parasitic and heterotrophic plants occurred approximately 56 Ma ago. PsaM, Psb30, ChlB, ChlL, ChlN, and Rpl21 were found to be characteristic signature genes of the chloroplast genome of algae, bryophytes, pteridophytes, and gymnosperms; however, none of these genes were found in the angiosperm or magnoliid lineage which appeared to have lost them approximately 203-156 Ma ago. A variety of chloroplast-encoded genes were lost across different species lineages throughout the evolutionary process. The Rpl20 gene, however, was found to be the most stable and intact gene in the chloroplast genome and was not lost in any of the analyzed species, suggesting that it is a signature gene of the plastome. Our evolutionary analysis indicated that chloroplast genomes evolved from multiple common ancestors ~1293 Ma ago and have undergone vivid recombination events across different taxonomic lineages.

Phylogeny of Euclidieae (Brassicaceae) based on plastome and nuclear ribosomal DNA data

Euclidieae, a morphologically diverse tribe in the family Brassicaceae (Cruciferae), consists of 29 genera and more than 150 species distributed mainly in Asia. Prior phylogenetic analyses on Euclidieae are inadequate. In this study, sequence data from the plastid genome and nuclear ribosomal DNA of 72 species in 27 genera of Euclidieae were used to infer the inter- and intra-generic relationships within. The well-resolved and strongly supported plastome phylogenies revealed that Euclidieae could be divided into five clades. Both Cymatocarpus and Neotorularia are polyphyletic in nuclear and plastome phylogenies. Besides, the conflicts of systematic positions of three species of Braya and two species of Solms-laubachia s.l. indicated that hybridization and or introgression might have happened during the evolutionary history of the tribe. Results from divergence-time analyses suggested an early Miocene origin of Euclidieae, and it probably originated from the Central Asia, Pamir Plateau and West Himalaya. In addition, multiple ndh genes loss and pseudogenization were detected in eight species based on comparative genomic study.

Chlororespiration Serves as Photoprotection for the Photo-Inactivated Oxygen-Evolving Complex in Zostera marina, a Marine Angiosperm

As an alternative electron sink, chlororespiration, comprising the NAD(P)H dehydrogenase complex and plastid terminal plastoquinone oxidase, may play a significant role in sustaining the redox equilibrium between stroma and thylakoid membrane. This study identified a distinct role for chlororespiration in the marine angiosperm Zostera marina, whose oxygen-evolving complex (OEC) is prone to photo-inactivation as a result of its inherent susceptibility to excess irradiation. The strong connectivity between OEC peripheral proteins and key chlororespiratory enzymes, as demonstrated in the interaction network of differentially expressed genes, suggested that the recovery of photo-inactivated OEC was connected with chlororespiration. Chlorophyll fluorescence, transcriptome and Western blot data verified a new physiological role for chlororespiration to function as photoprotection and generate a proton gradient across the thylakoid membrane for the recovery of photo-inactivated OEC. Chlororespiration was only activated in darkness following excess irradiation exposure, which might be related to electron deficiency in the electron transport chain because of the continuous impairment of the OEC. The activation of chlororespiration in Z. marina was prone to proactivity, which was also supported by the further activation of the oxidative pentose-phosphate pathway synthesizing NADPH to meet the demand of chlororespiration during darkness. This phenomenon is distinct from the common assumption that chlororespiration is prone to consuming redundant reducing power during the short transition phase from light to dark.

Plastome phylogenomic study of Gentianeae (Gentianaceae): widespread gene tree discordance and its association with evolutionary rate heterogeneity of plastid genes

BACKGROUND

Plastome-scale data have been prevalent in reconstructing the plant Tree of Life. However, phylogenomic studies currently based on plastomes rely primarily on maximum likelihood inference of concatenated alignments of plastid genes, and thus phylogenetic discordance produced by individual plastid genes has generally been ignored. Moreover, structural and functional characteristics of plastomes indicate that plastid genes may not evolve as a single locus and are experiencing different evolutionary forces, yet the genetic characteristics of plastid genes within a lineage remain poorly studied.

RESULTS

We sequenced and annotated 10 plastome sequences of Gentianeae. Phylogenomic analyses yielded robust relationships among genera within Gentianeae. We detected great variation of gene tree topologies and revealed that more than half of the genes, including one (atpB) of the three widely used plastid markers (rbcL, atpB and matK) in phylogenetic inference of Gentianeae, are likely contributing to phylogenetic ambiguity of Gentianeae. Estimation of nucleotide substitution rates showed extensive rate heterogeneity among different plastid genes and among different functional groups of genes. Comparative analysis suggested that the ribosomal protein (RPL and RPS) genes and the RNA polymerase (RPO) genes have higher substitution rates and genetic variations among plastid genes in Gentianeae. Our study revealed that just one (matK) of the three (matK, ndhB and rbcL) widely used markers show high phylogenetic informativeness (PI) value. Due to the high PI and lowest gene-tree discordance, rpoC2 is advocated as a promising plastid DNA barcode for taxonomic studies of Gentianeae. Furthermore, our analyses revealed a positive correlation of evolutionary rates with genetic variation of plastid genes, but a negative correlation with gene-tree discordance under purifying selection.

CONCLUSIONS

Overall, our results demonstrate the heterogeneity of nucleotide substitution rates and genetic characteristics among plastid genes providing new insights into plastome evolution, while highlighting the necessity of considering gene-tree discordance into phylogenomic studies based on plastome-scale data.