Plastome Evolution in Hemiparasitic Mistletoes

Phylogenetic placement of the monotypic Baolia (Amaranthaceae s.l.) based on morphological and molecular evidence

BACKGROUND

Baolia H.W.Kung & G.L.Chu is a monotypic genus only known in Diebu County, Gansu Province, China. Its systematic position is contradictory, and its morphoanatomical characters deviate from all other Chenopodiaceae. Recent study has regarded Baolia as a sister group to Corispermoideae. We therefore sequenced and compared the chloroplast genomes of this species, and resolved its phylogenetic position based on both chloroplast genomes and marker sequences.

RESULTS

We sequenced 18 chloroplast genomes of 16 samples from two populations of Baolia bracteata and two Corispermum species. These genomes of Baolia ranged in size from 152,499 to 152,508 bp. Simple sequence repeats (SSRs) were primarily located in the LSC region of Baolia chloroplast genomes, and most of them consisted of single nucleotide A/T repeat sequences. Notably, there were differences in the types and numbers of SSRs between the two populations of B. bracteata. Our phylogenetic analysis, based on both complete chloroplast genomes from 33 species and a combination of three markers (ITS, rbcL, and matK) from 91 species, revealed that Baolia and Corispermoideae (Agriophyllum, Anthochlamys, and Corispermum) form a well-supported clade and sister to Acroglochin. According to our molecular dating results, a major divergence event between Acroglochin, Baolia, and Corispermeae occurred during the Middle Eocene, approximately 44.49 mya. Ancestral state reconstruction analysis showed that Baolia exhibited symplesiomorphies with those found in core Corispermoideae characteristics including pericarp and seed coat.

CONCLUSIONS

Comparing the chloroplast genomes of B. bracteata with those of eleven typical Chenopodioideae and Corispermoideae species, we observed a high overall similarity and a one notable noteworthy case of inversion of approximately 3,100 bp. of DNA segments only in two Atriplex and four Chenopodium species. We suggest that Corispermoideae should be considered in a broader sense, it includes Corispermeae (core Corispermoideae: Agriophyllum, Anthochlamys, and Corispermum), as well as two new monotypic tribes, Acroglochineae (Acroglochin) and Baolieae (Baolia).

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.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

The complete chloroplast genome sequence of Psittacanthus schiedeanus (Cham. & Schltdl.) G.Don. (Santalales: Loranthaceae), the first plastome of a mistletoe species in the Psittacantheae tribe

Psittacanthus schiedeanus (Cham. & Schltdl.) G.Don., 1834, is a mistletoe species in the Loranthaceae, characteristic of the canopy in cloud forest edges and widely distributed in northern Mesoamerica. Here, we report the complete chloroplast genome sequence of P. schiedeanus, the first for a species in the Psittacantheae tribe. The circularized quadripartite structure of the P. schiedeanus chloroplast genome was 122,586 bp in length and included a large single-copy region of 72,507 bp and two inverted repeats of 21,283 bp separated by a small single-copy region of 7,513 bp. The genome contained 112 genes, of which 96 are unique, including 65 protein-coding genes, 27 transfer RNA, and four ribosomal RNA. The overall GC content in the plastome of P. schiedeanus is 36.9%. Based on 43 published complete chloroplast genome sequences for species in the families Loranthaceae and Santalaceae (Santalales), the maximum-likelihood phylogenetic tree with high-support bootstrap values indicated that P. schiedeanus in the Psittacantheae tribe is sister to the tribe Lorantheae. The chloroplast genome provided in this study represents a valuable resource for genetic, phylogenetic and conservation studies of Psittacanthus species, and an important advance for unraveling the evolutionary history of these hemiparasitic plants.

The Viscum album Gene Space database

The hemiparasitic flowering plant Viscum album (European mistletoe) is known for its very special life cycle, extraordinary biochemical properties, and extremely large genome. The size of its genome is estimated to be 30 times larger than the human genome and 600 times larger than the genome of the model plant Arabidopsis thaliana. To achieve insights into the Gene Space of the genome, which is defined as the space including and surrounding protein-coding regions, a transcriptome project based on PacBio sequencing has recently been conducted. A database resulting from this project contains sequences of 39,092 different open reading frames encoding 32,064 distinct proteins. Based on 'Benchmarking Universal Single-Copy Orthologs' (BUSCO) analysis, the completeness of the database was estimated to be in the range of 78%. To further develop this database, we performed a transcriptome project of V. album organs harvested in summer and winter based on Illumina sequencing. Data from both sequencing strategies were combined. The new V. album Gene Space database II (VaGs II) contains 90,039 sequences and has a completeness of 93% as revealed by BUSCO analysis. Sequences from other organisms, particularly fungi, which are known to colonize mistletoe leaves, have been removed. To evaluate the quality of the new database, proteome data of a mitochondrial fraction of V. album were re-analyzed. Compared to the original evaluation published five years ago, nearly 1000 additional proteins could be identified in the mitochondrial fraction, providing new insights into the Oxidative Phosphorylation System of V. album. The VaGs II database is available at https://viscumalbum.pflanzenproteomik.de/. Furthermore, all V. album sequences have been uploaded at the European Nucleotide Archive (ENA).

Habitat-related plastome evolution in the mycoheterotrophic Neottia listeroides complex (Orchidaceae, Neottieae)

BACKGROUND

Mycoheterotrophs, acquiring organic carbon and other nutrients from mycorrhizal fungi, have evolved repeatedly with substantial plastid genome (plastome) variations. To date, the fine-scale evolution of mycoheterotrophic plastomes at the intraspecific level is not well-characterized. A few studies have revealed unexpected plastome divergence among species complex members, possibly driven by various biotic/abiotic factors. To illustrate evolutionary mechanisms underlying such divergence, we analyzed plastome features and molecular evolution of 15 plastomes of Neottia listeroides complex from different forest habitats.

RESULTS

These 15 samples of Neottia listeroides complex split into three clades according to their habitats approximately 6 million years ago: Pine Clade, including ten samples from pine-broadleaf mixed forests, Fir Clade, including four samples from alpine fir forests and Fir-willow Clade with one sample. Compared with those of Pine Clade members, plastomes of Fir Clade members show smaller size and higher substitution rates. Plastome size, substitution rates, loss and retention of plastid-encoded genes are clade-specific. We propose to recognized six species in N. listeroides complex and slightly modify the path of plastome degradation.

CONCLUSIONS

Our results provide insight into the evolutionary dynamics and discrepancy of closely related mycoheterotrophic orchid lineages at a high phylogenetic resolution.

Understanding parasitism in Loranthaceae: Insights from plastome and mitogenome of Helicanthes elastica

Loranthaceae is the largest family of the order Santalales and includes root and stem hemiparasites. The parasites are known to exhibit reductions in the genomic features as well as relaxed or intensified selection shifts. In this study, we report plastome and mitogenome sequence of Helicanthes elastica (subtribe Amyeminae, tribe Lorantheae), an endemic, monotypic genus of Western Ghats, India growing on remarkably diverse host range. The length of plastome sequence was 1,28,805 bp while that of mitogenome was 1,65,273 bp. This is the smallest mitogenome from Loranthaceae reported till date. The plastome of Helicanthes exhibited loss of ndh genes (ψndhB), ψinfA, rps15, rps16, rpl32, trnK-UUU, trnG-UCC, trnV-UAC and trnA-UGC while mitogenome exhibited pseudogenized cox2, nad1 and nad4 genes. The comparative study of Loranthaceae plastomes revealed that the pseudogenization or loss of genes was not specific to any genus or tribe and variation was noted in the number of introns of clpP gene in the family. Several photosynthetic genes have undergone relaxed selection supporting lower photosynthetic rates in parasitic plants while some respiratory genes exhibited intensified selection supporting the idea of host-parasite arm race in Loranthaceae. The plastome gene content was found conserved in root hemiparasites compared to stem hemiparasites. The atp1 gene of mitogenome was chimeric and part of it exhibited similarities with Lamiales members. The phylogenetic analysis based on plastid genes placed Helicanthes sister to the members of subtribe Dendrophthoinae.

Gene Losses and Homology of the Chloroplast Genomes of Taxillus and Phacellaria Species

Research on the chloroplast genome of parasitic plants is limited. In particular, the homology between the chloroplast genomes of parasitic and hyperparasitic plants has not been reported yet. In this study, three chloroplast genomes of Taxillus (Taxillus chinensis, Taxillus delavayi, and Taxillus thibetensis) and one chloroplast genome of Phacellaria (Phacellaria rigidula) were sequenced and analyzed, among which T. chinensis is the host of P. rigidula. The chloroplast genomes of the four species were 119,941-138,492 bp in length. Compared with the chloroplast genome of the autotrophic plant Nicotiana tabacum, all of the ndh genes, three ribosomal protein genes, three tRNA genes and the infA gene were lost in the three Taxillus species. Meanwhile, in P. rigidula, the trnV-UAC gene and the ycf15 gene were lost, and only one ndh gene (ndhB) existed. The results of homology analysis showed that the homology between P. rigidula and its host T. chinensis was low, indicating that P. rigidula grows on its host T. chinensis but they do not share the chloroplast genome. In addition, horizontal gene transfer was not found between P. rigidula and its host T. chinensis. Several candidate highly variable regions in the chloroplast genomes of Taxillus and Phacellaria species were selected for species identification study. Phylogenetic analysis revealed that the species of Taxillus and Scurrula were closely related and supported that Scurrula and Taxillus should be treated as congeneric, while species in Phacellaria had a close relationship with that in Viscum.

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

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

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.

Comparing complete organelle genomes of holoparasitic Christisonia kwangtungensis (Orabanchaceae) with its close relatives: how different are they?

BACKGROUND

Orobanchaceae is the only flowering plant family with species from free-living nonparasite, hemi-parasite to holoparasite, making it an ideal system for studying the evolution of parasitism. However, both plastid and mitochondrial genome have been sequenced in only few parasitic species in Orobanchaceae. Therefore, further comparative study is wanted to investigate the impact of holoparasitism on organelle genomes evolution between close relatives. Here, we sequenced organelle genomes and transcriptome of holoparasitic Christisonia kwangtungensis and compared it with its closely related groups to analyze similarities and differences in adaption strategies to the holoparasitic lifestyle.

RESULTS

The plastid genome of C. kwangtungensis has undergone extensive pseudogenization and gene loss, but its reduction pattern is different from that of Aeginetia indica, the close relative of C. kwangtungensis. Similarly, the gene expression detected in the photosynthetic pathway of these two genera is different. In Orobanchaceae, holoparasites in Buchnereae have more plastid gene loss than Rhinantheae, which reflects their longer history of holoparasitism. Distinct from severe degradation of the plastome, protein-coding genes in the mitochondrial genome of C. kwangtungensis are relatively conserved. Interestingly, besides intracellularly transferred genes which are still retained in its plastid genome, we also found several horizontally transferred genes of plastid origin from diverse donors other than their current hosts in the mitochondrial genome, which probably indicate historical hosts.

CONCLUSION

Even though C. kwangtungensis and A. indica are closely related and share severe degradation of plastome, they adapt organelle genomes to the parasitic lifestyle in different ways. The difference between their gene loss and gene expression shows they ultimately lost photosynthetic genes but through different pathways. Our study exemplifies how parasites part company after achieving holoparasitism.

Complete chloroplast genome of the desert date (Balanites aegyptiaca (L.) Del. comparative analysis, and phylogenetic relationships among the members of Zygophyllaceae

Background

Balanites aegyptiaca (L.) Delile, commonly known as desert date, is a thorny evergreen tree belonging to the family Zygophyllaceae and subfamily Tribuloideae that is widespread in arid and semiarid regions. This plant is an important source of food and medicines and plays an important role in conservation strategies for restoring degraded desert ecosystems.

Results

In the present study, we sequenced the complete plastome of B. aegyptiaca. The chloroplast genome was 155,800 bp, with a typical four-region structure: a large single copy (LSC) region of 86,562 bp, a small single copy (SSC) region of 18,102 bp, and inverted repeat regions (IRa and IRb) of 25,568 bp each. The GC content was 35.5%. The chloroplast genome of B. aegyptiaca contains 107 genes, 75 of which coding proteins, 28 coding tRNA, and 4 coding rRNA.

We did not observe a large loss in plastid genes or a reduction in the genome size in B. aegyptiaca, as found previously in some species belonging to the family Zygophyllaceae. However, we noticed a divergence in the location of certain genes at the IR-LSC and IR-SSC boundaries and loss of ndh genes relative to other species. Furthermore, the phylogenetic tree constructed from the complete chloroplast genome data broadly supported the taxonomic classification of B. aegyptiaca as belonging to the Zygophyllaceae family. The plastome of B. aegyptiaca was found to be rich in single sequence repeats (SSRs), with a total of 240 SSRs.

Conclusions

The genomic data available from this study could be useful for developing molecular markers to evaluate population structure, investigate genetic variation, and improve production programs for B. aegyptiaca. Furthermore, the current data will support future investigation of the evolution of the family Zygophyllaceae.

Gene Losses and Plastome Degradation in the Hemiparasitic Species Plicosepalus acaciae and Plicosepalus curviflorus: Comparative Analyses and Phylogenetic Relationships among Santalales Members

The Plicosepalus genus includes hemiparasitic mistletoe and belongs to the Loranthaceae family, and it has several medicinal uses. In the present study, we sequenced the complete plastomes of two species, Plicosepalus acaciae and Plicosepalus curviflorus, and compared them with the plastomes of photosynthetic species (hemiparasites) and nonphotosynthetic species (holoparasites) in the order Santalales. The complete chloroplast genomes of P. acaciae and P. curviflorus are circular molecules with lengths of 120,181 bp and 121,086 bp, respectively, containing 106 and 108 genes and 63 protein-coding genes, including 25 tRNA and 4 rRNA genes for each species. We observed a reduction in the genome size of P. acaciae and P. curviflorus and the loss of certain genes, although this reduction was less than that in the hemiparasite and holoparasitic cp genomes of the Santalales order. Phylogenetic analysis supported the taxonomic state of P. acaciae and P. curviflorus as members of the family Loranthaceae and tribe Lorantheae; however, the taxonomic status of certain tribes of Loranthaceae must be reconsidered and the species that belong to it must be verified. Furthermore, available chloroplast genome data of parasitic plants could help to strengthen efforts in weed management and encourage biotechnology research to improve host resistance.

Transcriptional Basis for Haustorium Formation and Host Establishment in Hemiparasitic Psittacanthus schiedeanus Mistletoes

The mistletoe Psittacanthus schiedeanus, a keystone species in interaction networks between plants, pollinators, and seed dispersers, infects a wide range of native and non-native tree species of commercial interest. Here, using RNA-seq methodology we assembled the whole circularized quadripartite structure of P. schiedeanus chloroplast genome and described changes in the gene expression of the nuclear genomes across time of experimentally inoculated seeds. Of the 140,467 assembled and annotated uniGenes, 2,000 were identified as differentially expressed (DEGs) and were classified in six distinct clusters according to their expression profiles. DEGs were also classified in enriched functional categories related to synthesis, signaling, homoeostasis, and response to auxin and jasmonic acid. Since many orthologs are involved in lateral or adventitious root formation in other plant species, we propose that in P. schiedeanus (and perhaps in other rootless mistletoe species), these genes participate in haustorium formation by complex regulatory networks here described. Lastly, and according to the structural similarities of P. schiedeanus enzymes with those that are involved in host cell wall degradation in fungi, we suggest that a similar enzymatic arsenal is secreted extracellularly and used by mistletoes species to easily parasitize and break through tissues of the host.

Characterization of the complete chloroplast genome of Zephyranthes phycelloides (Amaryllidaceae, tribe Hippeastreae) from Atacama region of Chile

Sporadic rains in the Atacama Desert reveal a high biodiversity of plant species that only occur there. One of these rare species is the “Red añañuca” (Zephyranthes phycelloides), formerly known as Rhodophiala phycelloides. Many species of Zephyranthes in the Atacama Desert are dangerously threatened, due to massive extraction of bulbs and cutting of flowers. Therefore, studies of the biodiversity of these endemic species, which are essential for their conservation, should be conducted sooner rather than later. There are some chloroplast genomes available for Amaryllidaceae species, however there is no complete chloroplast genome available for any of the species of Zephyranthes subgenus Myostemma. The aim of the present work was to characterize and analyze the chloroplast of Z. phycelloides by NGS sequencing. The chloroplast genome of the Z. phycelloides consists of 158,107 bp, with typical quadripartite structures: a large single copy (LSC, 86,129 bp), a small single copy (SSC, 18,352 bp), and two inverted repeats (IR, 26,813 bp). One hundred thirty-seven genes were identified: 87 coding genes, 8 rRNA, 38 tRNA and 4 pseudogenes. The number of SSRs was 64 in Z. phycelloides and a total of 43 repeats were detected. The phylogenetic analysis of Z. phycelloides shows a distinct subclade with respect to Z. mesochloa. The average nucleotide variability (Pi) between Z. phycelloides and Z. mesochloa was of 0.02000, and seven loci with high variability were identified: psbA, trnS^GCU-trnG^UCC, trnD^GUC-trnY^GUA, trnL^UAA-trnF^GAA, rbcL, psbE-petL and ndhG-ndhI. The differences between the species are furthermore confirmed by the high amount of SNPs between these two species. Here, we report for the first time the complete cp genome of one species of the Zephyranthes subgenus Myostemma, which can be used for phylogenetic and population genomic studies.

The gene space of European mistletoe (Viscum album)

European mistletoe (Viscum album) is a hemiparasitic flowering plant that is known for its very special life cycle and extraordinary biochemical properties. Particularly, V. album has an unusual mode of cellular respiration that takes place in the absence of mitochondrial complex I. However, insights into the molecular biology of V. album so far are very limited. Since the genome of V. album is extremely large (estimated 600 times larger than the genome of the model plant Arabidopsis thaliana) it has not been sequenced up to now. We here report sequencing of the V. album gene space (defined as the space including and surrounding genic regions, encompassing coding as well as 5' and 3' non-coding regions). mRNA fractions were isolated from different V. album organs harvested in summer or winter and were analyzed via single-molecule real-time sequencing. We determined sequences of 39 092 distinct open reading frames encoding 32 064 V. album proteins (designated V. album protein space). Our data give new insights into the metabolism and molecular biology of V. album, including the biosynthesis of lectins and viscotoxins. The benefits of the V. album gene space information are demonstrated by re-evaluating mass spectrometry-based data of the V. album mitochondrial proteome, which previously had been evaluated using the A. thaliana genome sequence. Our re-examination allowed the additional identification of nearly 200 mitochondrial proteins, including four proteins related to complex I, which all have a secondary function not related to respiratory electron transport. The V. album gene space sequences are available at the NCBI.

Genes from oxidative phosphorylation complexes II-V and two dual-function subunits of complex I are transcribed in Viscum album despite absence of the entire mitochondrial holo-complex I

Mistletoes (Viscum) and close relatives are unique among flowering plants in having a drastically altered electron transport chain. Lack of complex I genes has previously been reported for the mitochondrial genome, and here we report an almost complete absence of nuclear-encoded complex I genes in the transcriptome of Viscum album. Compared to Arabidopsis with approximately 40 nuclear complex I genes, we recover only transcripts of two dual-function genes: gamma carbonic anhydrase and L-galactono-1,4-lactone dehydrogenase. The complement of genes belonging to complexes II-V of the oxidative phosphorylation pathway appears to be in accordance with other vascular plants. Additionally, transcripts encoding alternative NAD(P)H dehydrogenases and alternative oxidase were found. Despite sequence divergence, structural modeling suggests that the encoded proteins are structurally conserved. Complex I loss is a special feature in Viscum species and relatives, as all other parasitic flowering plants investigated to date seem to have a complete OXPHOS system. Hence, Viscum offers a unique system for specifically investigating molecular consequences of complex I absence, such as the role of complex I subunits involved in secondary functions.

Plastid NDH Pseudogenization and Gene Loss in a Recently Derived Lineage from the Largest Hemiparasitic Plant Genus Pedicularis (Orobanchaceae)

The plastid genome (plastome) is highly conserved in both gene order and content and has a lower mutation rate than the nuclear genome. However, the plastome is more variable in heterotrophic plants. To date, most such studies have investigated just a few species or only holoheterotrophic groups, and few have examined plastome evolution in recently derived lineages at an early stage of transition from autotrophy to heterotrophy. In this study, we investigated the evolutionary dynamics of plastomes in the monophyletic and recently derived Pedicularis sect. Cyathophora (Orobanchaceae). We obtained 22 new plastomes, 13 from the six recognized species of section Cyathophora, six from hemiparasitic relatives and three from autotrophic relatives. Comparative analyses of gene content, plastome structure and selection pressure showed dramatic differences among species in section Cyathophora and in Pedicularis as a whole. In comparison with autotrophic relatives and other Pedicularis spp., we found that the inverted repeat (IR) region in section Cyathophora had expansions to the small single-copy region, with a large expansion event and two independent contraction events. Moreover, NA(D)H dehydrogenase, accD and ccsA have lost function multiple times, with the function of accD being replaced by nuclear copies of an accD-like gene in Pedicularis spp. The ccsA and ndhG genes may have evolved under selection in association with IR expansion/contraction events. This study is the first to report high plastome variation in a recently derived lineage of hemiparasitic plants and therefore provides evidence for plastome evolution in the transition from autotrophy to heterotrophy.

Plastome variation and phylogeny of Taxillus (Loranthaceae)

Several molecular phylogenetic studies of the mistletoe family Loranthaceae have been published such that now the general pattern of relationships among the genera and their biogeographic histories are understood. Less is known about species relationships in the larger (> 10 species) genera. This study examines the taxonomically difficult genus Taxillus composed of 35–40 Asian species. The goal was to explore the genetic diversity present in Taxillus plastomes, locate genetically variable hotspots, and test these for their utility as potential DNA barcodes. Using genome skimming, complete plastomes, as well as nuclear and mitochondrial rDNA sequences, were newly generated for eight species. The plastome sequences were used in conjunction with seven publicly available Taxillus sequences and three sequences of Scurrula, a close generic relative. The Taxillus plastomes ranged from 121 to 123 kbp and encoded 90–93 plastid genes. In addition to all of the NADH dehydrogenase complex genes, four ribosomal genes, infA and four intron-containing tRNA genes were lost or pseudogenized in all of the Taxillus and Scurrula plastomes. The topologies of the plastome, mitochondrial rDNA and nuclear rDNA trees were generally congruent, though with discordance at the position of T. chinensis. Several variable regions in the plastomes were identified that have sufficient numbers of parsimony informative sites as to recover the major clades seen in the complete plastome tree. Instead of generating complete plastome sequences, our study showed that accD alone or the concatenation of accD and rbcL can be used in future studies to facilitate identification of Taxillus samples and to generate a molecular phylogeny with robust sampling within the genus.

Genomic reconfiguration in parasitic plants involves considerable gene losses alongside global genome size inflation and gene births

Parasitic plant genomes and transcriptomes reveal numerous genetic innovations, the functional-evolutionary relevance and roles of which open unprecedented research avenues.

Highly degenerate plastomes in two hemiparasitic dwarf mistletoes: Arceuthobium chinense and A. pini (Viscaceae)

The leafless and endophytic habitat may significantly relax the selection pressure on photosynthesis, and plastid transcription and translation, causing the loss/pseudogenization of several essential plastid-encoding genes in dwarf mistletoes. Dwarf mistletoes (Arceuthobium spp., Viscaceae) are the most destructive plant parasites to numerous conifer species worldwide. In this study, the plastid genomes (plastomes) of Arceuthobium chinense Lecomte and A. pini Hawksworth and Wiens were sequenced and characterized. Although dwarf mistletoes are hemiparasites capable of photosynthesis, their plastomes were highly degenerated, as indicated by the smallest plastome size, the lowest GC content, and relatively very few intact genes among the Santalales hemiparasites. Unexpectedly, several essential housekeeping genes (rpoA, rpoB, rpoC1, and rpoC2) and some core photosynthetic genes (psbZ and petL), as well as the rpl33 gene, that is indispensable for plants under stress conditions, were deleted or pseudogenized in the Arceuthobium plastomes. Our data suggest that the leafless and endophytic habit, which heavily relies on the coniferous hosts for nutrients and carbon requirement, may largely relax the selection pressure on photosynthesis, as well as plastid transcription and translation, thus resulting in the loss/pseudogenization of such essential plastid-encoding genes in dwarf mistletoes. Therefore, the higher level of plastome degradation in Arceuthobium species than other Santalales hemiparasites is likely correlated with the evolution of leafless and endophytic habit. A higher degree of plastome degradation in Arceuthobium. These findings provide new insights into the plastome degeneration associated with parasitism in Santalales and deepen our understanding of the biology of dwarf mistletoes.

Chloroplast Genomes of Two Species of Cypripedium: Expanded Genome Size and Proliferation of AT-Biased Repeat Sequences

The size of the chloroplast genome (plastome) of autotrophic angiosperms is generally conserved. However, the chloroplast genomes of some lineages are greatly expanded, which may render assembling these genomes from short read sequencing data more challenging. Here, we present the sequencing, assembly, and annotation of the chloroplast genomes of Cypripedium tibeticum and Cypripedium subtropicum. We de novo assembled the chloroplast genomes of the two species with a combination of short-read Illumina data and long-read PacBio data. The plastomes of the two species are characterized by expanded genome size, proliferated AT-rich repeat sequences, low GC content and gene density, as well as low substitution rates of the coding genes. The plastomes of C. tibeticum (197,815 bp) and C. subtropicum (212,668 bp) are substantially larger than those of the three species sequenced in previous studies. The plastome of C. subtropicum is the longest one of Orchidaceae to date. Despite the increase in genome size, the gene order and gene number of the plastomes are conserved, with the exception of an ∼75 kb large inversion in the large single copy (LSC) region shared by the two species. The most striking is the record-setting low GC content in C. subtropicum (28.2%). Moreover, the plastome expansion of the two species is strongly correlated with the proliferation of AT-biased non-coding regions: the non-coding content of C. subtropicum is in excess of 57%. The genus provides a typical example of plastome expansion induced by the expansion of non-coding regions. Considering the pros and cons of different sequencing technologies, we recommend hybrid assembly based on long and short reads applied to the sequencing of plastomes with AT-biased base composition.

The Organelle Genomes in the Photosynthetic Red Algal Parasite Pterocladiophila hemisphaerica (Florideophyceae, Rhodophyta) Have Elevated Substitution Rates and Extreme Gene Loss in the Plastid Genome

Comparative organelle genome studies of parasites can highlight genetic changes that occur during the transition from a free-living to a parasitic state. Our study focuses on a poorly studied group of red algal parasites, which are often closely related to their red algal hosts and from which they presumably evolved. Most of these parasites are pigmented and some show photosynthetic capacity. Here, we assembled and annotated the complete organelle genomes of the photosynthetic red algal parasite, Pterocladiophila hemisphaerica. The plastid genome is the smallest known red algal plastid genome at 68,701 bp. The plastid genome has many genes missing, including all photosynthesis-related genes. In contrast, the mitochondrial genome is similar in architecture to that of other free-living red algae. Both organelle genomes show elevated mutation rates and significant changes in patterns of selection, measured as dN/dS ratios. This caused phylogenetic analyses, even of multiple aligned proteins, to be unresolved or give contradictory relationships. Full plastid datasets interfered by selected best gene evolution models showed the supported relationship of P. hemisphaerica within the Ceramiales, but the parasite was grouped with support as sister to the Gracilariales when interfered under the GHOST model. Nuclear rDNA showed a supported grouping of the parasite within a clade containing several red algal orders including the Gelidiales. This photosynthetic parasite, which is unable to photosynthesize with its own plastid due to the total loss of all photosynthesis genes, raises intriguing questions on parasite-host organelle genome capabilities and interactions.

Comparative Plastome Analysis of Root- and Stem-Feeding Parasites of Santalales Untangle the Footprints of Feeding Mode and Lifestyle Transitions

In plants, parasitism triggers the reductive evolution of plastid genomes (plastomes). To disentangle the molecular evolutionary associations between feeding on other plants below- or aboveground and general transitions from facultative to obligate parasitism, we analyzed 34 complete plastomes of autotrophic, root- and stem-feeding hemiparasitic, and holoparasitic Santalales. We observed inexplicable losses of housekeeping genes and tRNAs in hemiparasites and dramatic genomic reconfiguration in holoparasitic Balanophoraceae, whose plastomes have exceptionally low GC contents. Genomic changes are related primarily to the evolution of hemi- or holoparasitism, whereas the transition from a root- to a stem-feeding mode plays no major role. In contrast, the rate of molecular evolution accelerates in a stepwise manner from autotrophs to root- and then stem-feeding parasites. Already the ancestral transition to root-parasitism coincides with a relaxation of selection in plastomes. Another significant selectional shift in plastid genes occurs as stem-feeders evolve, suggesting that this derived form coincides with trophic specialization despite the retention of photosynthetic capacity. Parasitic Santalales fill a gap in our understanding of parasitism-associated plastome degeneration. We reveal that lifestyle-genome associations unfold interdependently over trophic specialization and feeding mode transitions, where holoparasitic Balanophoraceae provide a system for exploring the functional realms of plastomes.

The Complete Plastomes of Five Hemiparasitic Plants (Osyris wightiana, Pyrularia edulis, Santalum album, Viscum liquidambaricolum, and V. ovalifolium): Comparative and Evolutionary Analyses Within Santalales

Most species of Santalales (the sandalwood order) are hemiparasites, including both facultative and obligate hemiparasites. Despite its rich diversity, only a small fraction of the species in the sandalwood order have sequenced plastomes. The evolution of parasitism-associated plastome reduction in Santalales remains under-studied. Here, we report the complete plastomes of three facultative hemiparasites (Pyrularia edulis, Cervantesiaceae; Osyris wightiana, and Santalum album, Santalaceae), and two obligate hemiparasites (Viscum liquidambaricolum and Viscum ovalifolium, Viscaceae). Coupled with publicly available data, we investigated the dynamics of plastome degradation in Santalales hemiparasites. Our results indicate that these hemiparasites can be characterized by various degrees of plastome downsizing, structural rearrangement, and gene loss. The loss or pseudogenization of ndh genes was commonly observed in Santalales hemiparasites, which may be correlated to the lifestyle shift from photoautotroph to hemiparasitism. However, the obligate hemiparasites did not exhibit a consistently higher level of gene loss or pseudogenization compared to facultative hemiparasites, which suggests that the degree of plastome reduction is not correlated with the trophic level facultative or obligate hemiparasitism. Instead, closely related taxa tend to possess highly similar plastome size, structure, and gene content. This implies the parasitism-associated plastome degradation in Santalales may evolve in a lineage-specific manner.

Photosynthetic cyclic electron transport provides ATP for homeostasis during trap closure in Dionaea muscipula

BACKGROUND AND AIMS

The processes connected with prey capture and the early consumption of prey by carnivorous Dionaea muscipula require high amounts of energy. The aim of the present study was to identify processes involved in flytrap energy provision and ATP homeostasis under these conditions.

METHODS

We determined photosynthetic CO2 uptake and chlorophyll fluorescence as well as the dynamics of ATP contents in the snap traps upon closure with and without prey.

KEY RESULTS

The results indicate that upon prey capture, a transient switch from linear to cyclic electron transport mediates a support of ATP homeostasis. Beyond 4 h after prey capture, prey resources contribute to the traps' ATP pool and, 24 h after prey capture, export of prey-derived resources to other plant organs may become preferential and causes a decline in ATP contents.

CONCLUSIONS

Apparently, the energy demand of the flytrap for prey digestion and nutrient mining builds on both internal and prey-derived resources.

The First Plastid Genome of the Holoparasitic Genus Prosopanche (Hydnoraceae)

Plastomes of parasitic and mycoheterotrophic plants show different degrees of reduction depending on the plants’ level of heterotrophy and host dependence in comparison to photoautotrophic sister species, and the amount of time since heterotrophic dependence was established. In all but the most recent heterotrophic lineages, this reduction involves substantial decrease in genome size and gene content and sometimes alterations of genome structure. Here, we present the first plastid genome of the holoparasitic genus Prosopanche, which shows clear signs of functionality. The plastome of Prosopanche americana has a length of 28,191 bp and contains only 24 unique genes, i.e., 14 ribosomal protein genes, four ribosomal RNA genes, five genes coding for tRNAs and three genes with other or unknown function (accD, ycf1, ycf2). The inverted repeat has been lost. Despite the split of Prosopanche and Hydnora about 54 MYA ago, the level of genome reduction is strikingly congruent between the two holoparasites although highly dissimilar nucleotide sequences are observed. Our results lead to two possible evolutionary scenarios that will be tested in the future with a larger sampling: 1) a Hydnoraceae plastome, similar to those of Hydnora and Prosopanche today, existed already in the most recent common ancestor and has not changed much with respect to gene content and structure, or 2) the genome similarities we observe today are the result of two independent evolutionary trajectories leading to almost the same endpoint. The first hypothesis would be most parsimonious whereas the second would point to taxon dependent essential gene sets for plants released from photosynthetic constraints.

The complete chloroplast genome of Magnolia polytepala: Comparative analyses offer implication for genetics and phylogeny of Yulania

Magnoliaceae is a primitive taxon in the angiosperms, comprising approximately 240 species in 2-17 genera. Many of them have been widely cultivated due to their horticultural and medicinal value. However, there are uncertainties and controversies about the delimitation of the genera except Liriodendron L. in this family. The Yulania taxa is also the focus of dispute at the genus and section levels. In this study, we compared ten Yulania plastomes, including the newly sequenced M. polytepala. The plastome-wide comparative analysis demonstrated that 1) Yulania cp genomes were highly conserved, and the majority differences existed in IR regions with the loss/retention of trnV-GAC or ycf15 gene, 2) mutational hotspots with high levels of nucleotide diversity (Pi > 0.02) existed in both coding (rpoA, and ycf1) and no-coding (ccsA-ndhD, ndhE-ndhG, ndhF-rpl32, petA-psbJ, rpl32-trnL, rps3-rps19, and trnH-psbA) regions among the genus Yulania. Combined with other data from Magnoliaceae plastomes, our reconstructed molecular phylogenetic tree revealed that Yulania is monophyletic, separated from the genus Magnolia L. (=Magnolia subg. Magnolia L.), but seems a sister of Michelia L. Moreover, M. polytepala which belongs to the genus Yulania is most closely related to M. liliiflora. All these results indicated that plastome data may contribute to investigating taxonomy, population genetics and phylogeny of Yulania.

Organellomic data sets confirm a cryptic consensus on (unrooted) land-plant relationships and provide new insights into bryophyte molecular evolution

PREMISE

Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution.

METHODS

We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements.

RESULTS

Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses.

CONCLUSIONS

A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure.

Gene Losses and Variations in Chloroplast Genome of Parasitic Plant Macrosolen and Phylogenetic Relationships within Santalales

Macrosolen plants are parasitic shrubs, several of which are important medicinal plants, that are used as folk medicine in some provinces of China. However, reports on Macrosolen are limited. In this study, the complete chloroplast genome sequences of Macrosolen cochinchinensis, Macrosolen tricolor and Macrosolen bibracteolatus are reported. The chloroplast genomes were sequenced by Illumina HiSeq X. The length of the chloroplast genomes ranged from 129,570 bp (M. cochinchinensis) to 126,621 bp (M. tricolor), with a total of 113 genes, including 35 tRNA, eight rRNA, 68 protein-coding genes, and two pseudogenes (ycf1 and rpl2). The simple sequence repeats are mainly comprised of A/T mononucleotide repeats. Comparative genome analyses of the three species detected the most divergent regions in the non-coding spacers. Phylogenetic analyses using maximum parsimony and maximum likelihood strongly supported the idea that Loranthaceae and Viscaceae are monophyletic clades. The data obtained in this study are beneficial for further investigations of Macrosolen in respect to evolution and molecular identification.

Mycoheterotrophic Epirixanthes (Polygalaceae) has a typical angiosperm mitogenome but unorthodox plastid genomes

BACKGROUND AND AIMS

Fully mycoheterotrophic plants derive carbon and other nutrients from root-associated fungi and have lost the ability to photosynthesize. While mycoheterotroph plastomes are often degraded compared with green plants, the effect of this unusual symbiosis on mitochondrial genome evolution is unknown. By providing the first complete organelle genome data from Polygalaceae, one of only three eudicot families that developed mycoheterotrophy, we explore how both organellar genomes evolved after loss of photosynthesis.

METHODS

We sequenced and assembled four complete plastid genomes and a mitochondrial genome from species of Polygalaceae, focusing on non-photosynthetic Epirixanthes. We compared these genomes with those of other mycoheterotroph and parasitic plant lineages, and assessed whether organelle genes in Epirixanthes experienced relaxed or intensified selection compared with autotrophic relatives.

KEY RESULTS

Plastomes of two species of Epirixanthes have become substantially degraded compared with that of autotrophic Polygala. Although the lack of photosynthesis is presumably homologous in the genus, the surveyed Epirixanthes species have marked differences in terms of plastome size, structural rearrangements, gene content and substitution rates. Remarkably, both apparently replaced a canonical plastid inverted repeat with large directly repeated sequences. The mitogenome of E. elongata incorporated a considerable number of fossilized plastid genes, by intracellular transfer from an ancestor with a less degraded plastome. Both plastid and mitochondrial genes in E. elongata have increased substitution rates, but the plastid genes of E. pallida do not. Despite this, both species have similar selection patterns operating on plastid housekeeping genes.

CONCLUSIONS

Plastome evolution largely fits with patterns of gene degradation seen in other heterotrophic plants, but includes highly unusual directly duplicated regions. The causes of rate elevation in the sequenced Epirixanthes mitogenome and of rate differences in plastomes of related mycoheterotrophic species are not currently understood.

Extensive plastome reduction and loss of photosynthesis genes in Diphelypaea coccinea, a holoparasitic plant of the family Orobanchaceae

BACKGROUND

Parasitic plants have the ability to obtain nutrients from their hosts and are less dependent on their own photosynthesis or completely lose this capacity. The reduction in plastid genome size and gene content in parasitic plants predominantly results from loss of photosynthetic genes. Plants from the family Orobanchaceae are used as models for studying plastid genome evolution in the transition from an autotrophic to parasitic lifestyle. Diphelypaea is a poorly studied genus of the Orobanchaceae, comprising two species of non-photosynthetic root holoparasites. In this study, we sequenced the plastid genome of Diphelypaea coccinea and compared it with other Orobanchaceae, to elucidate patterns of plastid genome evolution. In addition, we used plastid genome data to define the phylogenetic position of Diphelypaea spp.

METHODS

The complete nucleotide sequence of the plastid genome of D. coccinea was obtained from total plant DNA, using pyrosequencing technology.

RESULTS

The D. coccinea plastome is only 66,616 bp in length, and is highly rearranged; however, it retains a quadripartite structure. It contains only four rRNA genes, 25 tRNA genes and 25 protein-coding genes, being one of the most highly reduced plastomes among the parasitic Orobanchaceae. All genes related to photosynthesis, including the ATP synthase genes, had been lost, whereas most housekeeping genes remain intact. The plastome contains two divergent, but probably intact clpP genes. Intron loss had occurred in some protein-coding and tRNA genes. Phylogenetic analysis yielded a fully resolved tree for the Orobanchaceae, with Diphelypaea being a sister group to Orobanche sect. Orobanche.

Plastome Reduction in the Only Parasitic Gymnosperm Parasitaxus Is Due to Losses of Photosynthesis but Not Housekeeping Genes and Apparently Involves the Secondary Gain of a Large Inverted Repeat

Plastid genomes (plastomes) of parasitic plants undergo dramatic reductions as the need for photosynthesis relaxes. Here, we report the plastome of the only known heterotrophic gymnosperm Parasitaxus usta (Podocarpaceae). With 68 unique genes, of which 33 encode proteins, 31 tRNAs, and four rRNAs in a plastome of 85.3-kb length, Parasitaxus has both the smallest and the functionally least capable plastid genome of gymnosperms. Although the heterotroph retains chlorophyll, all genes for photosynthesis are physically or functionally lost, making photosynthetic energy gain impossible. The pseudogenization of the three plastome-encoded light-independent chlorophyll biosynthesis genes chlB, chlL, and chlN implies that Parasitaxus relies on either only the light-dependent chlorophyll biosynthesis pathway or another regulation system. Nesting within a group of gymnosperms known for the absence of the large inverted repeat regions (IRs), another unusual feature of the Parasitaxus plastome is the existence of a 9,256-bp long IR. Its short length and a gene composition that completely differs from those of IR-containing gymnosperms together suggest a regain of this critical, plastome structure-stabilizing feature. In sum, our findings highlight the particular path of lifestyle-associated reductive plastome evolution, where structural features might provide additional cues of a continued selection for plastome maintenance.

Thirteen New Plastid Genomes from Mixotrophic and Autotrophic Species Provide Insights into Heterotrophy Evolution in Neottieae Orchids

Mixotrophic species use both organic and mineral carbon sources. Some mixotrophic plants combine photosynthesis and a nutrition called mycoheterotrophy, where carbon is obtained from fungi forming mycorrhizal symbiosis with their roots. These species can lose photosynthetic abilities and evolve full mycoheterotrophy. Besides morphological changes, the latter transition is associated with a deep alteration of the plastid genome. Photosynthesis-related genes are lost first, followed by housekeeping genes, eventually resulting in a highly reduced genome. Whether relaxation of selective constraints already occurs for the plastid genome of mixotrophic species, which remain photosynthetic, is unclear. This is partly due to the difficulty of comparing plastid genomes of autotrophic, mixotrophic, and mycoheterotrophic species in a narrow phylogenetic framework. We address this question in the orchid tribe Neottieae, where this large assortment of nutrition types occurs. We sequenced 13 new plastid genomes, including 9 mixotrophic species and covering all 6 Neottieae genera. We investigated selective pressure on plastid genes in each nutrition type and conducted a phylogenetic inference of the group. Surprisingly, photosynthesis-related genes did not experience selection relaxation in mixotrophic species compared with autotrophic relatives. Conversely, we observed evidence for selection intensification for some plastid genes. Photosynthesis is thus still under purifying selection, maybe because of its role in fruit formation and thus reproductive success. Phylogenetic analysis resolved most relationships, but short branches at the base of the tree suggest an evolutionary radiation at the beginning of Neottieae history, which, we hypothesize, may be linked to mixotrophy emergence.

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

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

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

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

Mitochondrial genome evolution in parasitic plants

Background

Parasitic plants rely on their host to cover their nutritional requirements either for their entire life or a smaller part of it. Depending on the level of parasitism, a proportional reduction on the plastid genome has been found. However, knowledge on gene loss and evolution of the mitogenome of parasitic plants is only available for four hemiparasitic Viscum species (Viscaceae), which lack many of the mitochondrial genes, while the remaining genes exhibit very fast molecular evolution rates. In this study, we include another genus, Phoradendron, from the Viscaceae, as well as 10 other hemiparasitic or holoparasitic taxa from across the phylogeny of the angiosperms to investigate how fast molecular evolution works on their mitogenomes, and the extent of gene loss.

Results

Our observations from Viscum were replicated in Phoradendron liga, whereas the remaining parasitic plants in the study have a complete set of the core mitochondrial genes and exhibit moderate or only slightly raised substitution rates compared to most autotrophic taxa, without any statistically significant difference between the different groups (autotrophs, hemiparasites and holoparasites). Additionally, further evidence is provided for the placement of Balanophoraceae within the order Santalales, while the exact placement of Cynomoriaceae still remains elusive.

Conclusions

We examine the mitochondrial gene content of 11 hemiparasitic and holoparasitic plants and confirm previous observations in Viscaceae. We show that the remaining parasitic plants do not have significantly higher substitution rates than autotrophic plants in their mitochondrial genes. We provide further evidence for the placement of Balanophoraceae in the Santalales.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1401-8) contains supplementary material, which is available to authorized users.

Genome Reports: Contracted Genes and Dwarfed Plastome in Mycoheterotrophic Sciaphila thaidanica (Triuridaceae, Pandanales)

With a reduced need for photosynthesis, the plastome of parasitic and mycoheterotrophic plants degrades. In the tiny, fully mycoheterotrophic plant Sciaphila thaidanica, we find one of the smallest plastomes yet encountered. Its size is just 12,780 bp and it contains only 20 potentially functional housekeeping genes. Thus S. thaidanica fits the proposed model of gene loss in achlorophyllous plants. The most astonishing feature of the plastome is its extremely compact nature, with more than half of the genes having overlapping reading frames. Additionally, intergenic sequences have been reduced to a bare minimum, and the retained genes have been reduced in length both compared with the orthologous genes in another mycoheterotrophic species of Sciaphila and in the autotrophic relative Carludovica.

Convergent Plastome Evolution and Gene Loss in Holoparasitic Lennoaceae

The Lennoaceae, a small monophyletic plant family of root parasites endemic to the Americas, are one of the last remaining independently evolved lineages of parasitic angiosperms lacking a published plastome. In this study, we present the assembled and annotated plastomes of two species spanning the crown node of Lennoaceae, Lennoa madreporoides and Pholisma arenarium, as well as their close autotrophic relative from the sister family Ehretiaceae, Tiquilia plicata. We find that the plastomes of L. madreporoides and P. arenarium are similar in size and gene content, and substantially reduced compared to T. plicata, consistent with trends seen in other holoparasitic lineages. In particular, most plastid genes involved in photosynthesis function have been lost, whereas housekeeping genes (ribosomal protein-coding genes, rRNAs, and tRNAs) are retained. One notable exception is the persistence of a rbcL open reading frame in P. arenarium but not L. madreporoides suggesting a nonphotosynthetic function for this gene. Of the retained coding genes, d_N/d_S ratios indicate that some remain under purifying selection, whereas others show relaxed selection. Overall, this study supports the mounting evidence for convergent plastome evolution in flowering plants following the shift to heterotrophy.

Punctuated plastome reduction and host-parasite horizontal gene transfer in the holoparasitic plant genus Aphyllon

Foundational studies of chloroplast genome (plastome) evolution in parasitic plants have focused on broad trends across large clades, particularly among the Orobanchaceae, a species-rich and ecologically diverse family of root parasites. However, the extent to which such patterns and processes of plastome evolution, such as stepwise gene loss following the complete loss of photosynthesis (shift to holoparasitism), are detectable at shallow evolutionary time scale is largely unknown. We used genome skimming to assemble eight chloroplast genomes representing complete taxonomic sampling of Aphyllon sect. Aphyllon, a small clade within the Orobanchaceae that evolved approximately 6 Ma, long after the origin of holoparasitism. We show substantial plastome reduction occurred in the stem lineage, but subsequent change in plastome size, gene content, and structure has been relatively minimal, albeit detectable. This lends additional fine-grained support to existing models of stepwise plastome reduction in holoparasitic plants. Additionally, we report phylogenetic evidence based on an rbcL gene tree and assembled 60+ kb fragments of the Aphyllon epigalium mitochondrial genome indicating host-to-parasite horizontal gene transfers (hpHGT) of several genes originating from the plastome of an ancient Galium host into the mitochondrial genome of a recent common ancestor of A. epigalium Ecologically, this evidence of hpHGT suggests that the host-parasite associations between Galium and A. epigalium have been stable at least since its subspecies diverged hundreds of thousands of years ago.

Understanding plastome evolution in Hemiparasitic Santalales: Complete chloroplast genomes of three species, Dendrotrophe varians, Helixanthera parasitica, and Macrosolen cochinchinensis

Santalales is a large order, with over 2200 species, most of which are root or aerial (stem) hemiparasites. In this study, we report the newly assembled chloroplast genome of Dendrotrophe varians (140,666 bp) in the family Amphorogynaceae and the cp genomes of Helixanthera parasitica (124,881 bp) and Macrosolen cochinchinensis (122,986 bp), both in the family Loranthaceae. We compared the cp genomes of 11 Santalales including eight currently available cp genomes. Santalales cp genomes are slightly or not reduced in size (119-147 kb), similar to other hemiparasitic species, when compared with typical angiosperm cp genomes (120-170 kb). In a phylogeny examining gene content, the NADH dehydrogenase gene group is the only one among eight functional gene groups that lost complete functionally in all examined Santalales. This supports the idea that the functional loss of ndh genes is the initial stage in the evolution of the plastome of parasitic plants, but the loss has occurred independently multiple times in angiosperms, while they are not found in some parasites. This suggests that the functional loss of ndh genes is not essential for the transition from autotroph to parasite. We additionally examined the correlation between gene content and type of parasitism (obligate/facultative and stem/root parasites) of all hemiparasitic species in which cp genomes have been reported to date. Correlation was not found in any types of parasitism.

Plastid phylogenomics resolves infrafamilial relationships of the Styracaceae and sheds light on the backbone relationships of the Ericales

Relationships among the genera of the small, woody family Styracaceae and among families of the large, diverse order Ericales have resisted complete resolution with sequences from one or a few genes. We used plastome sequencing to attempt to resolve the backbone relationships of Styracaceae and Ericales and to explore plastome structural evolution. Complete plastomes for 23 species are newly reported here, including 18 taxa of Styracaceae and five of Ericales (including species of Sapotaceae, Clethraceae, Symplocaceae, and Diapensiaceae). Combined with publicly available complete plastome data, this resulted in a data set of 60 plastomes, including 11 of the 12 genera of Styracaceae and 12 of 22 families of Ericales. Styracaceae plastomes were found to possess the quadripartite structure typical of angiosperms, with sizes ranging from 155 to 159 kb. Most of the plastomes were found to possess the full complement of typical angiosperm plastome genes. Unusual structural features were detected in plastomes of Alniphyllum and Bruinsmia, including the presence of a large 20-kb inversion (14 genes) in the Large Single-Copy region, the loss or pseudogenization of the clpP and accD genes in Bruinsmia, and the loss of the first exon of rps16 in B. styracoides. Likewise, the second intron from clpP was found to be lost in Alniphyllum and Huodendron. Phylogenomic analyses including all 79 plastid protein-coding genes provided improved resolution for relationships among the genera of Styracaceae and families of Ericales. Styracaceae was strongly supported as monophyletic, with Styrax, Huodendron, and a clade of Alniphyllum + Bruinsmia successively sister to the remainder of the family, all with strong support. All genera of Styracaceae were recovered as monophyletic, except for Halesia and Pterostyrax, which were each recovered as polyphyletic with strong support. Within Ericales, all families were recovered as monophyletic with strong support, with Balsaminaceae sister to remaining Ericales. Most relationships recovered in plastome analyses are congruent with previous analyses based on smaller data sets. Our results demonstrate the power of plastid phylogenomics to improve phylogenetic hypotheses among genera and families, and provide new insight into plastome evolution across Ericales.

Gene losses and partial deletion of small single-copy regions of the chloroplast genomes of two hemiparasitic Taxillus species

Numerous variations are known to occur in the chloroplast genomes of parasitic plants. We determined the complete chloroplast genome sequences of two hemiparasitic species, Taxillus chinensis and T. sutchuenensis, using Illumina and PacBio sequencing technologies. These species are the first members of the family Loranthaceae to be sequenced. The complete chloroplast genomes of T. chinensis and T. sutchuenensis comprise circular 121,363 and 122,562 bp-long molecules with quadripartite structures, respectively. Compared with the chloroplast genomes of Nicotiana tabacum and Osyris alba, all ndh genes as well as three ribosomal protein genes, seven tRNA genes, four ycf genes, and the infA gene of these two species have been lost. The results of the maximum likelihood and neighbor-joining phylogenetic trees strongly support the theory that Loranthaceae and Viscaceae are monophyletic clades. This research reveals the effect of a parasitic lifestyle on the chloroplast structure and genome content of T. chinensis and T. sutchuenensis, and enhances our understanding of the discrepancies in terms of assembly results between Illumina and PacBio.

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

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

On the brink: the highly reduced plastomes of nonphotosynthetic Ericaceae

Ericaceae (the heather family) is a large and diverse group of plants that forms elaborate symbiotic relationships with mycorrhizal fungi, and includes several nonphotosynthetic lineages. Using an extensive sample of fully mycoheterotrophic (MH) species, we explored inter- and intraspecific variation as well as selective constraints acting on the plastomes of these unusual plants. The plastomes of seven MH genera were analysed in a phylogenetic context with two geographically disparate individuals sequenced for Allotropa, Monotropa, and Pityopus. The plastomes of nonphotosynthetic Ericaceae are highly reduced in size (c. 33-41 kbp) and content, having lost all photosynthesis-related genes, and are reduced to encoding housekeeping genes as well as a protease subunit (clpP)-like and acetyl-CoA carboxylase subunit D (accD)-like open reading frames. Despite an increase in the rate of their nucleotide substitutions, the remaining protein-coding genes are typically under purifying selection in full MHs. We also identified ribosomal proteins under relaxed or neutral selection. These plastomes also exhibit striking structural rearrangements. Intraspecific variation within MH Ericaceae ranges from a few differences (Allotropa) to extensive population divergences (Monotropa, Hypopitys), which indicates that cryptic speciation may be occurring in several lineages. The pattern of gene loss within fully MH Ericaceae plastomes suggests an advanced state of degradation.

Limited mitogenomic degradation in response to a parasitic lifestyle in Orobanchaceae

In parasitic plants, the reduction in plastid genome (plastome) size and content is driven predominantly by the loss of photosynthetic genes. The first completed mitochondrial genomes (mitogenomes) from parasitic mistletoes also exhibit significant degradation, but the generality of this observation for other parasitic plants is unclear. We sequenced the complete mitogenome and plastome of the hemiparasite Castilleja paramensis (Orobanchaceae) and compared them with additional holoparasitic, hemiparasitic and nonparasitic species from Orobanchaceae. Comparative mitogenomic analysis revealed minimal gene loss among the seven Orobanchaceae species, indicating the retention of typical mitochondrial function among Orobanchaceae species. Phylogenetic analysis demonstrated that the mobile cox1 intron was acquired vertically from a nonparasitic ancestor, arguing against a role for Orobanchaceae parasites in the horizontal acquisition or distribution of this intron. The C. paramensis plastome has retained nearly all genes except for the recent pseudogenization of four subunits of the NAD(P)H dehydrogenase complex, indicating a very early stage of plastome degradation. These results lend support to the notion that loss of ndh gene function is the first step of plastome degradation in the transition to a parasitic lifestyle.