Molecular evolution of rbcL in the mycoheterotrophic coralroot orchids (Corallorhiza Gagnebin, Orchidaceae)

The RuBisCO large subunit gene (rbcL) has been the focus of numerous plant phylogenetic studies and studies on molecular evolution in parasitic plants. However, there has been a lack of investigation of photosynthesis gene molecular evolution in fully mycoheterotrophic plants. These plants invade pre-existing mutualistic associations between ectomycorrhizal trees and fungi, from which they obtain fixed carbon and nutrients. The mycoheterotrophic orchid Corallorhiza contains both green (photosynthetic) and non-green (putatively nonphotosynthetic) species. We sequenced rbcL from 31 accessions of eight species of Corallorhiza and hypothesized that some lineages would have pseudogenes resulting from relaxation of purifying selection on RuBisCO's carboxylase function. Phylogenetic analysis of rbcL+ITS gave high jackknife support for relationships among species. We found evidence of pseudogene formation in all lineages of the Corallorhiza striata complex and in some lineages of the C. maculata complex. Evidence includes: stop codons, frameshifts, decreased d(S)/d(N) ratios, replacements not observed in photosynthetic species, rate heterogeneity, and high likelihood of neutral evolution. The evolution of rbcL in Corallorhiza may serve as an exemplary system in which to study the effects of relaxed evolutionary constraints on photosynthesis genes for >400 documented fully mycoheterotrophic plant species.     

Characterization of Mitochondrial Small-Subunit Ribosomal RNAs from Holoparasitic Plants

Mitochondrial small-subunit (19S) rDNA sequences were obtained from 10 angiosperms to further characterize sequence divergence levels and structural variation in this molecule. These sequences were derived from seven holoparasitic (nonphotosynthetic) angiosperms as well as three photosynthetic plants. 19S rRNA is composed of a conservative core region (ca. 1450 nucleotides) as well as two variable regions (V1 and V7). In pairwise comparisons of photosynthetic angiosperms to Glycine, the core 19S rDNA sequences differed by less than 1.4%, thus supporting the observation that variation in mitochondrial rDNA is 3–4 times lower than seen in protein coding and rDNA genes of other subcellular organelles. Sequences representing four distinct lineages of nonasterid holoparasites showed significantly increased numbers of substitutions in their core 19S rDNA sequences (2.3–7.6%), thus paralleling previous findings that showed accelerated rates in nuclear (18S) and plastid (16S) rDNA from the same plants. Relative rate tests confirmed the accelerated nucleotide substitution rates in the holoparasites whereas rates in nonparasitic plants were not significantly increased. Among comparisons of both parasitic and nonparasitic plants, transversions outnumbered transitions, in many cases more than two to one. The core 19S rRNA is conserved in sequence and structure among all nonparasitic angiosperms whereas 19S rRNA from members of holoparasitic Balanophoraceae have unique extensions to the V5 and V6 variable domains. Substitution and insertion/deletion mutations characterized the V1 and V7 regions of the nonasterid holoparasites. The V7 sequence of one holoparasite (Scybalium) contained repeat motifs. The cause of substitution rate increases in the holoparasites does not appear to be a result of RNA editing, hence the underlying molecular mechanism remains to be fully documented. Received: 18 May 1997 / Accepted: 11 July 1997     

Phylogenetic affinity of arbuscular mycorrhizal symbionts in <Emphasis Type="Italic">Psilotum nudum</Emphasis>

Many lineages of land plants (from lycopsids to angiosperms) have non-photosynthetic life cycle phases that involve obligate mycoheterotrophic arbuscular mycorrhizal (AM) associations where the plant host gains organic carbon through glomalean symbionts. Our goal was to isolate and phylogenetically identify the AM fungi associated with both the autotrophic and underground mycoheterotrophic life cycle phases of Psilotum nudum. Phylogenetic analyses recovered 11 fungal phylotypes in four diverse clades of Glomus A that form AM associations with P. nudum mycoheterotrophic gametophytes and autotrophic sporophytes, and angiosperm roots found in the same greenhouse pots. The correspondence of identities of AM symbionts in P. nudum sporophytes, gametophytes and neighboring angiosperms provides compelling evidence that photosynthetic heterospecific and conspecific plants can serve as the ultimate sources of fixed carbon for mycoheterotrophic gametophytes of P. nudum, and that the transfer of carbon occurs via shared fungal networks. Moreover, broader phylogenetic analyses suggest greenhouse Psilotum populations, like field-surveyed populations of mycoheterotrophic plants, form AM associations with restricted clades of Glomus A. The phylogenetic affinities and distribution of Glomus A symbionts indicate that P. nudum greenhouse populations have the potential to be exploited as an experimental system to further study the physiology, ecology and evolution of mycoheterotrophic AM associations. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

The plastid genome of the mycoheterotrophic Corallorhiza striata (Orchidaceae) is in the relatively early stages of degradation

? Premise of the study: Plastid genomes of nonphotosynthetic, mycoheterotrophic plants represent apt systems in which to study effects of relaxed evolutionary constraints. The few mycoheterotrophic angiosperm plastomes sequenced to date display drastic patterns of degradation/reduction relative to those of photosynthetic relatives. The goal of this study was to focus on a mycoheterotrophic orchid hypothesized to be in the "early" stages of plastome degradation, to provide perspective on this process. ? Methods: Short-read sequencing was used to generate a complete plastome sequence for Corallorhiza striata var. vreelandii, a mycoheterotrophic orchid, to investigate the extent of plastome degradation. Patterns of nonsynonymous/synonymous mutations were also assessed, and comparisons were made between Corallorhiza and other heterotrophic plant lineages. ? Key results: Corallorhiza yielded a plastome of 137505 bp, with several photosynthesis-related genes either lost or pseudogenized. Members of all major photosynthesis complexes, except ATP-synthase genes, were affected. "Housekeeping" genes were intact, despite the loss of a single tRNA. Intact photosynthesis genes (excluding atp genes) together displayed elevated nonsynonymous changes, while housekeeping genes did not. ? Conclusions: The Corallorhiza plastome is not drastically reduced in overall size (～6% reduction relative to that of photosynthetic Oncidium), but displays a pattern congruent with a loss of photosynthetic function. Comparing Corallorhiza with other heterotrophs allows some emergent evolutionary patterns to be inferred, but these remain as hypotheses to be tested, especially at lower taxonomic levels, and in lineages illustrating transitions from autotrophy to heterotrophy. The independent, unique processes of plastome modification among mycoheterotrophic lineages illustrate the urgency of their conservation.     

Plastid genome evolution in mycoheterotrophic Ericaceae

Unlike parasitic plants, which are linked to their hosts directly through haustoria, mycoheterotrophic (MHT) plants derive all or part of their water and nutrients from autothrophs via fungal mycorrhizal intermediaries. Ericaceae, the heather family, are a large and diverse group of plants known to form elaborate symbiotic relationships with mycorrhizal fungi. Using PHYA sequence data, we first investigated relationships among mycoheterotrophic Ericaceae and their close autotrophic relatives. Phylogenetic results suggest a minimum of two independent origins of MHT within this family. Additionally, a comparative investigation of plastid genomes (plastomes) grounded within this phylogenetic framework was conducted using a slot-blot Southern hybridization approach. This survey encompassed numerous lineages of Ericaceae with different life histories and trophic levels, including multiple representatives from mixotrophic Pyroleae and fully heterotrophic Monotropeae and Pterosporeae. Fifty-four probes derived from all categories of protein coding genes typically found within the plastomes of flowering plants were used. Our results indicate that the holo-mycoheterotrophic Ericaceae exhibit extensive loss of genes relating to photosynthetic function and expression of the plastome but retain genes with possible functions outside photosynthesis. Mixotrophic taxa tend to retain most genes relating to photosynthetic functions but are varied regarding the plastid ndh gene content. This investigation extends previous inferences that the loss of the NDH complex occurs prior to becoming holo-heterotrophic and it shows that the pattern of gene losses among mycoheterotrophic Ericaceae is similar to that of haustorial parasites. Additionally, we identify the most desirable candidate species for entire plastome sequencing.     

Mycoheterotrophic interactions are not limited to a narrow phylogenetic range of arbuscular mycorrhizal fungi

The majority of achlorophyllous mycoheterotrophic plant species associate with arbuscular mycorrhizal fungi (AMF). Previous studies have shown that some species are highly specialized towards narrow lineages of AMF and have suggested that only particular lineages of these fungi are targeted by mycoheterotrophic plants. To test this hypothesis, we analyzed all available partial SSU sequences of AMF associated with mycoheterotrophic plants including data from 13 additional specimens from French Guiana, Gabon and Australia. Sequences were assigned to 'virtual taxa' (VT) according to the MaarjAM database. We found that 20% of all known Glomeromycota VT are involved in mycoheterotrophic interactions and the majority of associations involve Glomeraceae (Glomus Group A) fungi. While some mycoheterotrophic plant species have been found growing with only a single VT, many species are able to associate with a wide range of AMF. We calculated significant phylogenetic clustering of Glomeromycota VT involved in mycoheterotrophic interactions, suggesting that associations between mycoheterotrophic plants and AMF are influenced by the phylogenetic relationships of the fungi. Our results demonstrate that many lineages of AMF are prone to exploitation by mycoheterotrophic plants. However, mycoheterotrophs from different plant lineages and different geographical regions tend to be dependent on lineages of AMF that are phylogenetically related.     

Photosynthetic evolution in parasitic plants: insight from the chloroplast genome

Despite the enormous diversity in plant form, structure and growth environment across the seed-bearing plants (angiosperms and gymnosperms), the chloroplast genome has, with few exceptions, remained remarkably conserved. This conservation suggests the existence of universal evolutionary selection pressures associated with photosynthesis-the primary function of chloroplasts. The stark exceptions to this conservation occur in parasitic angiosperms, which have escaped the dominant model by evolving the capacity to obtain some or all of their carbon (and nutrients) from their plant hosts. The consequence of this evolution to parasitism is a relaxation of the evolutionary constraints associated with the need to maintain photosynthetic function, the very function that drove early stages of the ancient symbiotic relationship that produced the contemporary chloroplast. Extreme examples of reductionism among parasitic angiosperms reveals major alterations in chloroplast function with the loss of photosynthetic capacity and, with that, massive alterations in chloroplast genome content. This review highlights emerging patterns in reported gene loss and gene retention in the chloroplast genomes of parasitic plants. Some gene losses appear to occur in the early stages of parasitic evolution, even before the loss of photosynthetic capacity, like the chlororespiratory (ndh) genes. This contrasts with unexpected gene retentions, like that of the rbcL gene responsible for photosynthetic carbon dioxide fixation, and belies current understanding of gene function. The review relates gene retention to current knowledge of protein function and gene processing that has implications to broader aspects of genome conservation in organelles.     

A chromosome-scale Gastrodia elata genome and large-scale comparative genomic analysis indicate convergent evolution by gene loss in mycoheterotrophic and parasitic plants

Mycoheterotrophic and parasitic plants are heterotrophic and parasitize on fungi and plants, respectively, to obtain nutrients. Large-scale comparative genomics analysis has not been conducted in mycoheterotrophic or parasitic plants or between these two groups of parasites. We assembled a chromosome-level genome of the fully mycoheterotrophic plant Gastrodia elata (Orchidaceae) and performed comparative genomic analyses on the genomes of G. elata and four orchids (initial mycoheterotrophs), three parasitic plants (Cuscuta australis, Striga asiatica, and Sapria himalayana), and 36 autotrophs from various angiosperm lineages. It was found that while in the hemiparasite S. asiatica and initial mycoheterotrophic orchids, approximately 4–5% of the conserved orthogroups were lost, the fully heterotrophic G. elata and C. australis both lost approximately 10% of the conserved orthogroups, indicating that increased heterotrophy is positively associated with gene loss. Importantly, many genes that are essential for autotrophs, including those involved in photosynthesis, the circadian clock, flowering time regulation, immunity, nutrient uptake, and root and leaf development, were convergently lost in both G. elata and C. australis. The high-quality genome of G. elata will facilitate future studies on the physiology, ecology, and evolution of mycoheterotrophic plants, and our findings highlight the critical role of gene loss in the evolution of plants with heterotrophic lifestyles.     

A bi-organellar phylogenomic study of Pandanales: inference of higher-order relationships and unusual rate-variation patterns

We used a bi-organellar phylogenomic approach to address higher-order relationships in Pandanales, including the first molecular phylogenetic study of the panama-hat family, Cyclanthaceae. Our genus-level study of plastid and mitochondrial gene sets includes a comprehensive sampling of photosynthetic lineages across the order, and provides a framework for investigating clade ages, biogeographic hypotheses and organellar molecular evolution. Using multiple inference methods and both organellar genomes, we recovered mostly congruent and strongly supported relationships within and between families, including the placement of fully mycoheterotrophic Triuridaceae. Cyclanthaceae and Pandanaceae plastomes have slow substitution rates, contributing to weakly supported plastid-based relationships in Cyclanthaceae. While generally slowly evolving, mitochondrial genomes exhibit sporadic rate elevation across the order. However, we infer well-supported relationships even for slower evolving mitochondrial lineages in Cyclanthaceae. Clade age estimates across photosynthetic lineages are largely consistent with previous studies, are well correlated between the two organellar genomes (with slightly younger inferences from mitochondrial data), and support several biogeographic hypotheses. We show that rapidly evolving non-photosynthetic lineages may bias age estimates upwards at neighbouring photosynthetic nodes, even using a relaxed clock model. Finally, we uncovered new genome structural variants in photosynthetic taxa at plastid inverted repeat boundaries that show promise as interfamilial phylogenetic markers.     

Arbuscular mycorrhizal symbionts in Botrychium (Ophioglossaceae)

Many plant species are characterized by a life cycle with a long-lived, subterranean phase that is completely dependent on mycorrhizal fungal symbionts for fixed carbon. This type of life cycle is both phylogenetically and ecologically widespread and is found in diverse vascular plant lineages from the tropics to subalpine meadows. Here we report on the molecular identities of the arbuscular mycorrhizal fungi associated with the autotrophic and underground mycoheterotrophic life cycle phases of the ferns Botrychium crenulatum and B. lanceolatum. We show that the Glomus taxa found in the mycoheterotrophic life cycle phases of B. crenulatum and B. lanceolatum are also found in conspecific and heterospecific photosynthetic neighboring plants. From our DNA sequence data, we infer carbon flow from photosynthetic plants to mycoheterotrophic plants through shared glomalean fungal networks. Finally, our phylogenetic analyses identify a major Glomus clade that forms associations with mycoheterotrophic life cycle stages of B. crenulatum and B. lanceolatum.     

Conservation and Structural Divergence of Organellar DNA and Gene Expression in Non-Photosynthetic Plastids During Ontogenetic Differentiation and Phylogenetic Adaptation

Plastids of non-photosynthetic cells or tissues, such as chromoplasts or leukoplasts, which develop during the course of ontogenetic differentiation contain DNA which is identical to chloroplast DNA with respect to size, organization and gene content. Also in ribosome-deficient bleached plastids, produced in leaves by experimental treatments or mutation, chloroplast DNA remains unaltered. The chloroplast DNA of various bleached mutant strains of Euglena has suffered major deletions or rearrangements, but is, however, never totally lost. Also leukoplasts of parasitic higher plants contain DNA. In the organellar DNA of several parasitic plants photosynthetic genes are conserved. In the heterotrophic flagellate Astasia and in the holoparasite Epifagus virginiana (Orobanchaceae) the size of the plastid DNA is greatly reduced by major deletions and most or all photosynthetic genes or genes related to the chloroplastic respiratory chain are lost. The residual plastid genomes have, however, retained genes for RNAs, tRNAs and ribosomal polypeptides and these are transcribed, although plastidic RNA-polymerase genes are lost in Epifagus. These findings demand the existence of a nuclear-encoded RNA-polymerase. The relevance of the conservation of plastid DNA and of plastidic gene expression in non-photosynthetic cells is discussed, remains, however, at present elusive. Open reading frames of unknown function might be of particular significance for non-photosynthetic plastids.     

Changing partners in the dark: isotopic and molecular evidence of ectomycorrhizal liaisons between forest orchids and trees

In the mycorrhizal symbiosis, plants exchange photosynthates for mineral nutrients acquired by fungi from the soil. This mutualistic arrangement has been subverted by hundreds of mycorrhizal plant species that lack the ability to photosynthesize. The most numerous examples of this behaviour are found in the largest plant family, the Orchidaceae. Although these non-photosynthetic orchid species are known to be highly specialized exploiters of the ectomycorrhizal symbiosis, photosynthetic orchids are thought to use free-living saprophytic, or pathogenic, fungal lineages. However, we present evidence that putatively photosynthetic orchids from five species which grow in the understorey of forests: (i) form mycorrhizas with ectomycorrhizal fungi of forest trees; and (ii) have stable isotope signatures indicating distinctive pathways for nitrogen and carbon acquisition approaching those of non-photosynthetic orchids that associate with ectomycorrhizal fungi of forest trees. These findings represent a major shift in our understanding of both orchid ecology and evolution because they explain how orchids can thrive in low-irradiance niches and they show that a shift to exploiting ectomycorrhizal fungi precedes viable losses of photosynthetic ability in orchid lineages.     

Relative rates of synonymous substitutions in the mitochondrial, chloroplast and nuclear genomes of seed plants

Previous studies have estimated that, in angiosperms, the synonymous substitution rate of chloroplast genes is three times higher than that of mitochondrial genes and that of nuclear genes is twelve times higher than that of mitochondrial genes. Here we used 12 genes in 27 seed plant species to investigate whether these relative rates of substitutions are common to diverse seed plant groups. We find that the overall relative rate of synonymous substitutions of mitochondrial, chloroplast and nuclear genes of all seed plants is 1:3:10, that these ratios are 1:2:4 in gymnosperms but 1:3:16 in angiosperms and that they go up to 1:3:20 in basal angiosperms. Our results show that the mitochondrial, chloroplast and nuclear genomes of seed plant groups have different synonymous substitutions rates, that these rates are different in different seed plant groups and that gymnosperms have smaller ratios than angiosperms.     

Ins and outs of plastid genome evolution

Recent findings have established cracks in the straight-laced image of the plastid genome as a molecule whose sole function is photosynthesis and whose gene content is highly conserved. Genes for numerous non-photosynthetic functions have been identified. Algal plastid genomes contain many genes with no homologs in angiosperms, and the recent transfer of genes from the plastid to the nuclear genome has been described. Wholesale abandonment of genes encoding photosynthetic and gene-expression functions has occurred in the plastid genomes of a non-green plant and alga. The origins of plastid DNA, its use in phylogenetic studies, and the origins of plastid introns are also reviewed.     

The effect of relaxed functional constraints on the photosynthetic gene rbcL in photosynthetic and nonphotosynthetic parasitic plants

The photosynthetic gene rbcL has been lost or dramatically altered in some lineages of nonphotosynthetic parasitic plants, but the dynamics of these events following loss of photosynthesis and whether rbcL has sustained functionally significant changes in photosynthetic parasitic plants are unknown. To assess the changes to rbcL associated with the loss of functional constraints for photosynthesis, nucleotide sequences from nonparasitic and parasitic plants of Scrophulariales were used for phylogeny reconstruction and character analysis. Plants in this group display a broad range of parasitic abilities, from photosynthetic ("hemiparasites") to nonphotosynthetic ("holoparasites"). With the exception of Conopholis (Orobanchaceae), the rbcL locus is present in all parasitic plants of Scrophulariales examined. Several holoparasitic genera included in this study, including Boschniakia, Epifagus, Orobanche, and Hyobanche, have rbcL pseudogenes. However, the holoparasites Alectra orobanchoides, Harveya capensis, Harveya purpurea, Lathraea clandestina, Orobanche corymbosa, O. fasciculata, and Striga gesnerioides have intact open reading frames (ORFs) for the rbcL gene. Phylogenetic hypotheses based on rbcL are largely in agreement with those based on sequences of the nonphotosynthetic genes rps2 and matK and show a single origin of parasitism, and loss of photosynthesis and pseudogene formation have been independently derived several times in Scrophulariales. The mutations in rbcL in nonparasitic and hemiparasitic plants would result in largely conservative amino acid substitutions, supporting the hypothesis that functional proteins can experience only a limited range of changes, even in minimally photosynthetic plants. In contrast, ORFs in some holoparasites had many previously unobserved missense substitutions at functionally important amino acid residues, suggesting that rbcL genes in these plants have evolved under relaxed or altered functional constraints.      

Evolution of the plastid ribosomal RNA operon in a nongreen parasitic plant: Accelerated sequence evolution,altered promoter structure,and tRNA pseudogenes

The nucleotide sequence of a 7.4 kb region containing the entire plastid ribosomal RNA operon of the nongreen parasitic plant Epifagus virginiana has been determined. Analysis of the sequence indicates that all four rRNA genes are intact and almost certainly functional. In contrast, the split genes for tRNA^Ile and tRNA^Ala present in the 16S-23S rRNA spacer region have become pseudogenes, and deletion upstream of the 16S rRNA gene has removed a tRNA^Val gene and most of the promoter region for the rRNA operon. The rate of nucleotide substitution in 16S and 23S rRNAs is several times higher in Epifagus than in tobacco, a related photosynthetic plant. Possible reasons for this, including relaxed translational constraints, are discussed.     

Angiosperm phylogeny inferred from sequences of four mitochondrial genes

An angiosperm phylogeny was reconstructed in a maximum likelihood analysis of sequences of four mitochondrial genes, atpl, matR, had5, and rps3, from 380 species that represent 376 genera and 296 families of seed plants. It is largely congruent with the phylogeny of angiosperms reconstructed from chloroplast genes atpB, matK, and rbcL, and nuclear 18S rDNA. The basalmost lineage consists of Amborella and Nymphaeales (including Hydatellaceae). Austrobaileyales follow this clade and are sister to the mesangiosperms, which include Chloranthaceae, Ceratophyllum, magnoliids, monocots, and eudicots. With the exception of Chloranthaceae being sister to Ceratophyllum, relationships among these five lineages are not well supported. In eudicots, Ranunculales, Sabiales, Proteales, Trochodendrales, Buxales, Gunnerales, Saxifragales, Vitales, Berberidopsidales, and Dilleniales form a basal grade of lines that diverged before the diversification of rosids and asterids. Within rosids, the COM (Celastrales-Oxalidales-Malpighiales) clade is sister to malvids (or rosid Ⅱ), instead of to the nitrogen-fixing clade as found in all previous large-scale molecular analyses of angiosperms. Santalales and Caryophyllales are members of an expanded asterid clade. This study shows that the mitochondrial genes are informative markers for resolving relationships among genera, families, or higher rank taxa across angiosperms. The low substitution rates and low homoplasy levels of the mitochondrial genes relative to the chloroplast genes, as found in this study, make them particularly useful for reconstructing ancient phylogenetic relationships. A mitochondrial gene-based angiosperm phylogeny provides an independent and essential reference for comparison with hypotheses of angiosperm phylogeny based on chloroplast genes, nuclear genes, and non-molecular data to reconstruct the underlying organismal phylogeny.     

Complex Patterns of Plastid 16S rRNA Gene Evolution in Nonphotosynthetic Green Algae

This study provides a phylogenetic/comparative approach to deciphering the processes underlying the evolution of plastid rRNA genes in genomes under relaxed functional constraints. Nonphotosynthetic green algal taxa that belong to two distinct classes, Chlorophyceae (Polytoma) and Trebouxiophyceae (Prototheca), were investigated. Similar to the situation described previously for plastid 16S rRNA genes in nonphotosynthetic land plants, nucleotide substitution levels, extent of structural variations, and percentage AT values are increased in nonphotosynthetic green algae compared to their closest photosynthetic relatives. However, the mutational processes appear to be different in many respects. First, with the increase in AT content, more transversions are noted in Polytoma and holoparasite angiosperms, while more transitions characterize the evolution of the 16S rDNA sequences in Prototheca. Second, although structural variations do accumulate in both Polytoma and Prototheca (as well as holoparasitic plastid 16S rRNAs), insertions as large as 1.6 kb characterize the plastid 16S rRNA genes in the former, whereas significantly smaller indels (not exceeding 24 bp) seem to be more prevalent in the latter group. The differences in evolutionary rates and patterns within and between lineages might be due to mutations in replication/repair-related genes; slipped-strand mispairing is likely the mechanism responsible for the expansion of insertions in Polytoma plastid 16S rRNA genes. Received: 29 December 2000 / Accepted: 18 May 2001