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

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.

The complete chloroplast genome of Rhododendronambiguum and comparative genomics of related species

Rhododendron Linnaeus, 1753, the largest genus of woody plants in the Northern Hemisphere, includes some of the most significant species in horticulture. Rhododendronambiguum Hemsl, 1911, a member of subsection Triflora Sleumer 1947, exemplifies typical alpine Rhododendron species. The analysis of the complete chloroplast genome of R.ambiguum offers new insights into the evolution of Rhododendron species and enhances the resolution of phylogenetic relationships. This genome is composed of 207,478 base pairs, including a pair of inverted repeats (IRs) of 47,249 bp each, separated by a large single-copy (LSC) region of 110,367 bp and a small single-copy (SSC) region of 2,613 bp. It contains 110 genes: 77 protein-coding genes, 29 tRNAs, four unique rRNAs (4.5S, 5S, 16S, and 23S), with 16 genes duplicated in the IRs. Comparative analyses reveal substantial diversity in the Rhododendron chloroplast genome structures, identifying a fourth variant pattern. Specifically, four highly divergent regions (trnI-rpoB, ndhE-psaC, rpl32-ndhF, rrn16S-trnI) were noted in the intergenic spacers. Additionally, 76 simple sequence repeats were identified. Positive selection signals were detected in four genes (cemA, rps4, rpl16, and rpl14), evidenced by high Ka/Ks ratios. Phylogenetic reconstruction based on two datasets (shared protein-coding genes and complete chloroplast genomes) suggests that R.ambiguum is closely related to R.concinnum Hemsley, 1889. However, the phylogenetic positions of subsection Triflora Pojarkova, 1952 species remain unresolved, indicating that the use of complete chloroplast genomes for phylogenetic research in Rhododendron requires careful consideration. Overall, our findings provide valuable genetic information that will enhance understanding of the evolution, molecular biology, and genetic improvement of Rhododendron spieces.

The complete chloroplast genome sequences of four Liparis species (Orchidaceae) and phylogenetic implications

Liparis Richard (Malaxideae, Epidendroideae) is a large and diverse genus of the family Orchidaceae, the taxonomy of which is complicated and controversial. In this study, we sequenced, assembled and analyzed four complete chloroplast genomes of Liparis species including L. kumokiri, L. makinoana, L. pauliana, and L. viridiflora, and evaluated their phylogenetic relationships with related species for the first time. These four chloroplast genomes (size range 153,095 to 158,239 bp) possess typical quadripartite structures that consist of a large single copy (LSC, 83,533-86,752 bp), a small single copy (SSC, 17,938-18,156 bp) and a pair of inverted repeats (IRs, 26,421-26,933 bp). The genomes contain 133 genes, including 87 protein coding genes, 38 tRNAs and 8 rRNA genes. The genome arrangements, gene contents, gene order, long repeats and simple sequence repeats were similar with small differences observed among these four chloroplast genomes. Five highly variable regions including ycf1, ndhA, ndhF, trnQ and trnK were identified from the comparative analysis with other nine related Liparis species, which had the potential to be used as DNA markers for species identification and phylogenetic studies of Liparis species. Phylogenetic analysis based on the complete chloroplast genome sequences strongly supported the polyphyly of Liparis and its further division into three branches. These results provided valuable information to illustrate the complicated taxonomy, phylogeny and evolution process of the Liparis genus.

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.

Plastome evolution and phylogenomics of Impatiens (Balsaminaceae)

This study reported seven new plastomes from Impatiens and observed three highly variable regions for phylogeny and DNA barcoding, which resolved the relationships among sections of subgenus Impatiens. Impatiens L. (Balsaminaceae, Ericales) is one of the largest and most diverse genera of angiosperms, widely known for its taxonomic difficulty. In this study, we reevaluated the infrageneric relationships within the genus Impatiens, using complete plastome sequence data. Seven complete plastomes of Impatiens (representing 6 species) were newly sequenced and characterized along with 20 previously published plastomes of other Impatiens species, plus 2 plastomes of outgroups (Hydrocera triflora, Balsaminaceae; Marcgravia coriacea, Marcgraviaceae). The total size of these 29 plastomes ranged from 151,538 bp to 152,917 bp, except 2 samples of Impatiens morsei, which exhibited a shorter length and lost some genes encoding NADH dehydrogenase subunits. Moreover, the number of simple sequence repeats (SSRs) ranged from 51 to 113, and the number of long repeats from 17 to 26. In addition, three highly variable regions were identified (trnG-GCC (The previous one), ndhF-rpl32-trnL-UGA-ccsA, and ycf1). Our phylogenomic analysis based on 80 plastome-derived protein-coding genes strongly supported the monophyly of Impatiens and its two subgenera (Clavicarpa and Impatiens), and fully resolved relationships among the six (out of seven) sampled sections of subgenus Impatiens. Overall, the plastome DNA markers and phylogenetic results reported in this study will facilitate future identification, taxonomic and DNA barcoding studies in Impatiens as well as evolutionary studies in Balsaminaceae.

Omics Approaches in Uncovering Molecular Evolution and Physiology of Botanical Carnivory

Systems biology has been increasingly applied with multiple omics for a holistic comprehension of complex biological systems beyond the reductionist approach that focuses on individual molecules. Different high-throughput omics approaches, including genomics, transcriptomics, metagenomics, proteomics, and metabolomics have been implemented to study the molecular mechanisms of botanical carnivory. This covers almost all orders of carnivorous plants, namely Caryophyllales, Ericales, Lamiales, and Oxalidales, except Poales. Studies using single-omics or integrated multi-omics elucidate the compositional changes in nucleic acids, proteins, and metabolites. The omics studies on carnivorous plants have led to insights into the carnivory origin and evolution, such as prey capture and digestion as well as the physiological adaptations of trap organ formation. Our understandings of botanical carnivory are further enhanced by the discoveries of digestive enzymes and transporter proteins that aid in efficient nutrient sequestration alongside dynamic molecular responses to prey. Metagenomics studies revealed the mutualistic relationships between microbes and carnivorous plants. Lastly, in silico analysis accelerated the functional characterization of new molecules from carnivorous plants. These studies have provided invaluable molecular data for systems understanding of carnivorous plants. More studies are needed to cover the diverse species with convergent evolution of botanical carnivory.

The photosynthesis apparatus of European mistletoe (Viscum album)

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

The complete mitochondrial genome of carnivorous Genlisea tuberosa (Lentibulariaceae): Structure and evolutionary aspects

Sequenced genomic data for carnivorous plants are scarce, especially regarding the mitogenomes (MTs) and further studies are crucial to obtain a better understanding of the topic. In this study, we sequenced and characterized the mitochondrial genome of the tuberous carnivorous plant Genlisea tuberosa, being the first of its genus to be sequenced. The genome comprises 729,765 bp, encoding 80 identified genes of which 36 are protein-coding, 40 tRNA, four rRNA genes, and three pseudogenes. An intronic region from the cox1 gene was identified that encodes an endonuclease enzyme that is present in the other sequenced species of Lentibulariaceae. Chloroplast genes (pseudogene and complete) inserted in the MT genome were identified, showing possible horizontal transfer between organelles. In addition, 50 pairs of long repeats from 94 to 274 bp are present, possibly playing an important role in the maintenance of the MT genome. Phylogenetic analysis carried out with 34 coding mitochondrial genes corroborated the positioning of the species listed here within the family. The molecular dynamism in the mitogenome (e.g. the loss or pseudogenization of genes, insertion of foreign genes, the long repeats as well as accumulated mutations) may be reflections of the carnivorous lifestyle where a significant part of cellular energy was shifted for the adaptation of leaves into traps molding the mitochondrial DNA. The sequence and annotation of G. tuberosa's MT will be useful for further studies and serve as a model for evolutionary and taxonomic clarifications of the group as well as improving our comprehension of MT evolution.

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.

Calophyllaceae plastomes, their structure and insights in relationships within the clusioids

A complete chloroplast genome is not yet available for numerous species of plants. Among the groups that lack plastome information is the clusioid clade (Malpighiales), which includes five families: Bonnetiaceae, Calophyllaceae, Clusiaceae, Hypericaceae, and Podostemaceae. With around 2200 species, it has few published plastomes and most of them are from Podostemaceae. Here we assembled and compared six plastomes from members of the clusioids: five from Calophyllaceae (newly sequenced) and one from Clusiaceae. Putative regions for evolutionary studies were identified and the newly assembled chloroplasts were analyzed with other available chloroplasts for the group, focusing on Calophyllaceae. Our results mostly agree with recent studies which found a general conserved structure, except for the two Podostemaceae species that have a large inversion (trnK-UUU–rbcL) and lack one intron from ycf3. Within Calophyllaceae we observed a longer LSC and reduced IRs in Mahurea exstipulata, resulting in some genic rearrangement, and a short inversion (psbJ–psbE) in Kielmeyera coriacea. Phylogenetic analyses recovered the clusioids and the five families as monophyletic and revealed that conflicts in relationships reported in the literature for the group agree with nodes concentrating uninformative or conflicting gene trees. Our study brings new insights about clusioid plastome architecture and its evolution.

Structural Features of Carnivorous Plant (Genlisea, Utricularia) Tubers as Abiotic Stress Resistance Organs

Carnivorous plants from the Lentibulariaceae form a variety of standard and novel vegetative organs and survive unfavorable environmental conditions. Within Genlisea, only G. tuberosa, from the Brazilian Cerrado, formed tubers, while Utricularia menziesii is the only member of the genus to form seasonally dormant tubers. We aimed to examine and compare the tuber structure of two taxonomically and phylogenetically divergent terrestrial carnivorous plants: Genlisea tuberosa and Utriculariamenziesii. Additionally, we analyzed tubers of U. mannii. We constructed phylogenetic trees using chloroplast genes matK/trnK and rbcL and used studied characters for ancestral state reconstruction. All examined species contained mainly starch as histologically observable reserves. The ancestral state reconstruction showed that specialized organs such as turions evolved once and tubers at least 12 times from stolons in Lentibulariaceae. Different from other clades, tubers probably evolved from thick stolons for sect. Orchidioides and both structures are primarily water storage structures. In contrast to species from section Orchidioides, G. tuberosa, U.menziesii and U. mannii form starchy tubers. In G. tuberosa and U. menziesii, underground tubers provide a perennating bud bank that protects the species in their fire-prone and seasonally desiccating environments.

Plastome comparative genomics in maples resolves the infrageneric backbone relationships

Maples (Acer) are among the most diverse and ecologically important tree genera of the north-temperate forests. They include species highly valued as ornamentals and as a source of timber and sugar products. Previous phylogenetic studies employing plastid markers have not provided sufficient resolution, particularly at deeper nodes, leaving the backbone of the maple plastid tree essentially unresolved. We provide the plastid genome sequences of 16 species of maples spanning the sectional diversity of the genus and explore the utility of these sequences as a source of information for genetic and phylogenetic studies in this group. We analyzed the distribution of different types of repeated sequences and the pattern of codon usage, and identified variable regions across the plastome. Maximum likelihood and Bayesian analyses using two partitioning strategies were performed with these and previously published sequences. The plastomes ranged in size from 155,212 to 157,023 bp and had structure and gene content except for Acer palmatum (sect. Palmata), which had longer inverted repeats and an additional copy of the rps19 gene. Two genes, rps2 and rpl22, were found to be truncated at different positions and might be non-functional in several species. Most dispersed repeats, SSRs, and overall variation were detected in the non-coding sequences of the LSC and SSC regions. Fifteen loci, most of which have not been used before in the genus, were identified as the most variable and potentially useful as molecular markers for barcoding and genetic studies. Both ML and Bayesian analyses produced similar results irrespective of the partitioning strategy used. The plastome-based tree largely supported the topology inferred in previous studies using cp markers while providing resolution to the backbone relationships but was highly incongruous with a recently published nuclear tree presenting an opportunity for further research to investigate the causes of discordance, and particularly the role of hybridization in the diversification of the genus. Plastome sequences are valuable tools to resolve deep-level relationships within Acer. The variable loci and SSRs identified in this study will facilitate the development of markers for ecological and evolutionary studies in the genus. This study underscores the potential of plastid genome sequences to improve our understanding of the evolution of maples.

Genlisea hawkingii (Lentibulariaceae), a new species from Serra da Canastra, Minas Gerais, Brazil

Genlisea hawkingii, which is a new species of Genlisea subgen. Tayloria (Lentibulariaceae) from cerrado in southwest Brazil, is described and illustrated. This species has been found in only one locality thus far, in the Serra da Canastra, which is located in the Delfinópolis municipality in Minas Gerais, Brazil. The new species is morphologically similar to Genlisea violacea and G. flexuosa, but differs from them in having a corolla with a conical and curved spur along with sepals with an acute apex and reproductive organs that only have glandular hairs. Moreover, it is similar to G. uncinata's curved spur. G. hawkingii is nested within the subgen. Tayloria clade as a sister group to all the other species of this subgenus. Therefore, both morphological and phylogenetic results strongly support G. hawkingii as a new species in the subgen. Tayloria.

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.