Complete Chloroplast Genome Sequence of Holoparasite Cistanche deserticola (Orobanchaceae) Reveals Gene Loss and Horizontal Gene Transfer from Its Host Haloxylon ammodendron (Chenopodiaceae)

Background

The central function of chloroplasts is to carry out photosynthesis, and its gene content and structure are highly conserved across land plants. Parasitic plants, which have reduced photosynthetic ability, suffer gene losses from the chloroplast (cp) genome accompanied by the relaxation of selective constraints. Compared with the rapid rise in the number of cp genome sequences of photosynthetic organisms, there are limited data sets from parasitic plants.

Principal Findings/Significance

Here we report the complete sequence of the cp genome of Cistanche deserticola, a holoparasitic desert species belonging to the family Orobanchaceae. The cp genome of C. deserticola is greatly reduced both in size (102,657 bp) and in gene content, indicating that all genes required for photosynthesis suffer from gene loss and pseudogenization, except for psbM. The striking difference from other holoparasitic plants is that it retains almost a full set of tRNA genes, and it has lower dN/dS for most genes than another close holoparasitic plant, E. virginiana, suggesting that Cistanche deserticola has undergone fewer losses, either due to a reduced level of holoparasitism, or to a recent switch to this life history. We also found that the rpoC2 gene was present in two copies within C. deserticola. Its own copy has much shortened and turned out to be a pseudogene. Another copy, which was not located in its cp genome, was a homolog of the host plant, Haloxylon ammodendron (Chenopodiaceae), suggesting that it was acquired from its host via a horizontal gene transfer.     

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.      

Rampant gene loss in the underground orchid Rhizanthella gardneri highlights evolutionary constraints on plastid genomes

Since the endosymbiotic origin of chloroplasts from cyanobacteria 2 billion years ago, the evolution of plastids has been characterized by massive loss of genes. Most plants and algae depend on photosynthesis for energy and have retained ～110 genes in their chloroplast genome that encode components of the gene expression machinery and subunits of the photosystems. However, nonphotosynthetic parasitic plants have retained a reduced plastid genome, showing that plastids have other essential functions besides photosynthesis. We sequenced the complete plastid genome of the underground orchid, Rhizanthella gardneri. This remarkable parasitic subterranean orchid possesses the smallest organelle genome yet described in land plants. With only 20 proteins, 4 rRNAs, and 9 tRNAs encoded in 59,190 bp, it is the least gene-rich plastid genome known to date apart from the fragmented plastid genome of some dinoflagellates. Despite numerous differences, striking similarities with plastid genomes from unrelated parasitic plants identify a minimal set of protein-encoding and tRNA genes required to reside in plant plastids. This prime example of convergent evolution implies shared selective constraints on gene loss or transfer.     

Phytochrome evolution in green and nongreen plants

Photoreceptors are critical molecules that function at the interface between organism and environment. Plants use specific light signals to determine their place in time and space, allowing them to synchronize their growth, metabolism, and development to the environments in which they occur. Thus, innovation in light sensing mechanisms is expected to coincide with adaptation and diversification. Three studies involving the well-characterized phytochrome photoreceptor system in plants indicate that much work is yet needed to test this expectation. In early diverging flowering plants, episodic positive selection influenced the evolution of phytochrome A, but little of the functional data needed to link molecular adaptation with a change in gene function are available. In the model plant Arabidopsis thaliana, known functional differences between a recently duplicated gene pair remain difficult to characterize at the sequence level. In parasitic plants, patterns of development that in autotrophs are under the control of light signals are highly modified, suggesting that phytochromes and other photoreceptors function differently in nonphotosynthetic plants. Analyses of phytochrome A coding sequences indicate that they are evolving under relaxed constraints in nonphotosynthetic Orobanchaceae, consistent with the expectation of functional change. Further work is needed to determine which of the processes mediated by phyA may have been altered, a line of investigation that may improve our understanding of divergence points in downstream signaling pathways.     

Rapid evolution of the plastid translational apparatus in a nonphotosynthetic plant: Loss or accelerated sequence evolution of tRNA and ribosomal protein genes

Summary The vestigial plastid genome of Epifagus virginiana (beechdrops), a nonphotosynthetic parasitic flowering plant, is functional but lacks six ribosomal protein and 13 tRNA genes found in the chloroplast DNAs of photosynthetic flowering plants. Import of nuclear gene products is hypothesized to compensate for many of these losses. Codon usage and amino acid usage patterns in Epifagus plastic genes have not been affected by the tRNA gene losses, though a small shift in the base composition of the whole genome (toward A + T -richness) is apparent. The ribosomal protein and tRNA genes that remain have had a high rate of molecular evolution, perhaps due to relaxation of constraints on the translational apparatus. Despite the compactness and extensive gene loss, one translational gene (infA, encoding initiation factor 1) that is a pseudogene in tobacco has been maintained intact in Epifagus.Offprint requests to: J.D. Palmer     

Power analysis of tests for loss of selective constraint in cave crayfish and nonphotosynthetic plant lineages

Loss of selective constraint on a gene may be expected following changes in the environment or life history that render its function unnecessary. The long-term persistence of protein-coding genes after the loss of known functional necessity can occur by chance or because of selective maintenance of an unknown gene function. The selective maintenance of an alternative gene function is not demonstrated by the failure of statistical tests to reject the hypothesis that there has been no change in the degree of constraint on the evolution of coding genes. Maintenance may be inferred, however, when power analyses of such tests demonstrate that there has been a sufficient number of nucleotide substitutions to detect the loss of selective constraint. Here, we describe a power analysis for tests of loss of constraint on protein-coding genes. The power analysis was applied to loss-of-constraint tests for opsin gene evolution in cave-dwelling crayfish and rbcL evolution in nonphotosynthetic parasitic plants. The power of previously applied tests for loss of constraint on cave crayfish opsin genes was insufficient to distinguish between chance retention and selective maintenance of opsin genes. However, the power of codon-based likelihood ratio tests for change in d(N)/d(S) (=omega) (nonsynonymous to synonymous change) did have sufficient power to detect a loss of constraint on rbcL associated with a loss of photosynthesis in most examples but failed to detect such a change in three independent lineages. We conclude that rbcL has been selectively maintained in these holoparasitic plant lineages. This conclusion suggests that either these taxa are photosynthetic for at least a part of their life or rbcL may have an unknown function in these plants unrelated to photosynthesis.     

Purifying selection detected in the plastid gene matK and flanking ribozyme regions within a group II intron of nonphotosynthetic plants

In a striking contrast, matK is one of the most rapidly evolving plastid genes and also one of the few plastid genes to be retained in all nonphotosynthetic plants examined to date. DNA sequences of this region were obtained from photosynthetic and nonphotosynthetic plants of Orobanchaceae and their relatives. The resulting plastid DNA phylogeny was congruent with that recently obtained from analyses of rps2 and provided much better resolution. This phylogeny was then used to examine the relative degrees of evolutionary constraint of both the matK gene and the non-protein-coding regions that flank it inside the trnK intron. The method of subtree contrasts was introduced to compare levels of constraint. matK has evolved with a low but significant level of constraint on its amino acid sequence in both photosynthetic and nonphotosynthetic plants. Constraint is greater in photosynthetic than in nonphotosynthetic plants of this group. Domain X, thought to contain the active site of the protein, is not significantly more constrained than the rest of the protein. The portions of the flanking regions that are thought to form paired stem structures also show constraint, but in this case, there is no significant difference in degree of constraint between photosynthetic and nonphotosynthetic plants.     

Discovery of photosynthesis genes through whole-genome sequencing of acetate-requiring mutants of Chlamydomonas reinhardtii

Large-scale mutant libraries have been indispensable for genetic studies, and the development of next-generation genome sequencing technologies has greatly advanced efforts to analyze mutants. In this work, we sequenced the genomes of 660 Chlamydomonas reinhardtii acetate-requiring mutants, part of a larger photosynthesis mutant collection previously generated by insertional mutagenesis with a linearized plasmid. We identified 554 insertion events from 509 mutants by mapping the plasmid insertion sites through paired-end sequences, in which one end aligned to the plasmid and the other to a chromosomal location. Nearly all (96%) of the events were associated with deletions, duplications, or more complex rearrangements of genomic DNA at the sites of plasmid insertion, and together with deletions that were unassociated with a plasmid insertion, 1470 genes were identified to be affected. Functional annotations of these genes were enriched in those related to photosynthesis, signaling, and tetrapyrrole synthesis as would be expected from a library enriched for photosynthesis mutants. Systematic manual analysis of the disrupted genes for each mutant generated a list of 253 higher-confidence candidate photosynthesis genes, and we experimentally validated two genes that are essential for photoautotrophic growth, CrLPA3 and CrPSBP4. The inventory of candidate genes includes 53 genes from a phylogenomically defined set of conserved genes in green algae and plants. Altogether, 70 candidate genes encode proteins with previously characterized functions in photosynthesis in Chlamydomonas, land plants, and/or cyanobacteria; 14 genes encode proteins previously shown to have functions unrelated to photosynthesis. Among the remaining 169 uncharacterized genes, 38 genes encode proteins without any functional annotation, signifying that our results connect a function related to photosynthesis to these previously unknown proteins. This mutant library, with genome sequences that reveal the molecular extent of the chromosomal lesions and resulting higher-confidence candidate genes, will aid in advancing gene discovery and protein functional analysis in photosynthesis.     

Nonessential plastid-encoded ribosomal proteins in tobacco: a developmental role for plastid translation and implications for reductive genome evolution

Plastid genomes of higher plants contain a conserved set of ribosomal protein genes. Although plastid translational activity is essential for cell survival in tobacco (Nicotiana tabacum), individual plastid ribosomal proteins can be nonessential. Candidates for nonessential plastid ribosomal proteins are ribosomal proteins identified as nonessential in bacteria and those whose genes were lost from the highly reduced plastid genomes of nonphotosynthetic plastid-bearing lineages (parasitic plants, apicomplexan protozoa). Here we report the reverse genetic analysis of seven plastid-encoded ribosomal proteins that meet these criteria. We have introduced knockout alleles for the corresponding genes into the tobacco plastid genome. Five of the targeted genes (ribosomal protein of the large subunit22 [rpl22], rpl23, rpl32, ribosomal protein of the small subunit3 [rps3], and rps16) were shown to be essential even under heterotrophic conditions, despite their loss in at least some parasitic plastid-bearing lineages. This suggests that nonphotosynthetic plastids show elevated rates of gene transfer to the nuclear genome. Knockout of two ribosomal protein genes, rps15 and rpl36, yielded homoplasmic transplastomic mutants, thus indicating nonessentiality. Whereas Δrps15 plants showed only a mild phenotype, Δrpl36 plants were severely impaired in photosynthesis and growth and, moreover, displayed greatly altered leaf morphology. This finding provides strong genetic evidence that chloroplast translational activity influences leaf development, presumably via a retrograde signaling pathway.     

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.     

The complete plastomes of two flowering epiparasites (Phacellaria glomerata and P. compressa): Gene content,organization, and plastome degradation

A plant parasite obligately parasitizing another plant parasite is referred to as epiparasite, which is extremely rare in angiosperms, and their complete plastome sequences have not been characterized to date. In this study, the complete plastomes of two flowering epiparasites: Phacellaria compressa and P. glomerata (Amphorogynaceae, Santalales) were sequenced. The plastomes of both species are of similar size, structure, gene content, and arrangement of genes to other hemiparasites in Santalales. Their plastomes were characterized by the functional loss of plastid-encoded NAD(P)H-dehydrogenase and infA genes, which strongly coincides with the general pattern of plastome degradation observed in Santalales hemiparasites. Our study demonstrates that the relatively higher level of nutritional reliance on the host plants and the reduced vegetative bodies of P. compressa and P. glomerata do not appear to cause any unique plastome degradation compared with their closely related hemiparasites.     

Evolution of glutamine synthetase in heterokonts: evidence for endosymbiotic gene transfer and the early evolution of photosynthesis

Although the endosymbiotic evolution of chloroplasts through primary and secondary associations is well established, the evolutionary timing and stability of the secondary endosymbiotic events is less well resolved. Heterokonts include both photosynthetic and nonphotosynthetic members and the nonphotosynthetic lineages branch basally in phylogenetic reconstructions. Molecular and morphological data indicate that heterokont chloroplasts evolved via a secondary endosymbiosis, involving a heterotrophic host cell and a photosynthetic ancestor of the red algae and this endosymbiotic event may have preceded the divergence of heterokonts and alveolates. If photosynthesis evolved early in this lineage, nuclear genomes of the nonphotosynthetic groups may contain genes that are not essential to photosynthesis but were derived from the endosymbiont genome through gene transfer. These genes offer the potential to trace the evolutionary history of chloroplast gains and losses within these lineages. Glutamine synthetase (GS) is essential for ammonium assimilation and glutamine biosynthesis in all organisms. Three paralogous gene families (GSI, GSII, and GSIII) have been identified and are broadly distributed among prokaryotic and eukaryotic lineages. In diatoms (Heterokonta), the nuclear-encoded chloroplast and cytosolic-localized GS isoforms are encoded by members of the GSII and GSIII family, respectively. Here, we explore the evolutionary history of GSII in both photosynthetic and nonphotosynthetic heterokonts, red algae, and other eukaryotes. GSII cDNA sequences were obtained from two species of oomycetes by polymerase chain reaction amplification. Additional GSII sequences from eukaryotes and bacteria were obtained from publicly available databases and genome projects. Bayesian inference and maximum likelihood phylogenetic analyses of GSII provided strong support for the monophyly of heterokonts, rhodophytes, chlorophytes, and plants and strong to moderate support for the Opisthokonts. Although the phylogeny is reflective of the unikont/bikont division of eukaryotes, we propose based on the robustness of the phylogenetic analyses that the heterokont GSII gene evolved via endosymbiotic gene transfer from the nucleus of the red-algal endosymbiont to the nucleus of the host. The lack of GSIII sequences in the oomycetes examined here further suggests that the GSIII gene that functions in the cytosol of photosynthetic heterokonts was replaced by the endosymbiont-derived GSII gene.     

Lineage‐specific plastid degradation in subtribe Gentianinae (Gentianaceae)

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

Piecing together the puzzle of parasitic plant plastome evolution

The importance of photosynthesis as a mode of energy production has put plastid genomes of plants under a constant purifying selection. This has shaped the characteristic features of plastid genomes across the entire spectrum of photosynthetic plants and has led to a highly uniform and conserved plastid genome with respect to structure, size, gene order, intron and editing site positions and coding capacity. Parasitic species that have dropped photosynthesis as the main energy provider share striking deviations from the plastid genome norm: multiple rearrangements within the circular chromosome, pseudogenization and gene deletions, promoter losses, intron losses as well as the extensive loss of mRNA editing competence have been reported. The collective loss of larger sets of functionally related genes like those for the plastid NADH–dehydrogenase complex and concomitant losses of RNA polymerase genes together with their target promoters point to “domino effects” where an initial loss might have triggered others. An example, which will be discussed in more detail, is the concomitant loss of the intron maturase gene matK and all introns that are supposedly subject to MatK-dependent splicing in two Cuscuta species.     

Targeted disruption of chloroplast genes in Chlamydomonas reinhardtii.

Summary We have developed an efficient procedure for the disruption of Chlamydomonas chloroplast genes. Wild-type C. reinhardtii cells were bombarded with microprojectiles coated with a mixture of two plasmids, one encoding selectable, antibiotic-resistance mutations in the 16S ribosomal RNA gene and the other containing either the atpB or rbcL photosynthetic gene inactivated by an insertion of 0.48 kb of yeast DNA in the coding sequence. Antibiotic-resistant transformants were selected under conditions permissive for growth of nonphotosynthetic mutants. Approximately half of these transformants were initially heteroplasmic for copies of the disrupted atpB or rbcL genes integrated into the recipient chloroplast genome but still retained photosynthetic competence. A small fraction of the transformants (1.1% for atpB; 4.3% for rbcL) were nonphotosynthetic and homoplasmic for the disrupted gene at the time they were isolated. Single cell cloning of the initially heteroplasmic transformants also yielded nonphotosynthetic segregants that were homoplasmic for the disrupted gene. Polypeptide products of the disrupted atpB and rbcL genes could not be detected using immunoblotting techniques. We believe that any nonessential Chlamydomonas chloroplast gene, such as those involved in photosynthesis, should be amenable to gene disruption by cotransformation. The method should prove useful for the introduction of site-specific mutations into chloroplast genes and flanking regulatory sequences with a view to elucidating their function.     

Do nonasterid holoparasitic flowering plants have plastid genomes?

Past work involving the plastid genome (plastome) of holoparasitic plants has been confined to Scrophulariaceae (or Orobanchaceae) which have truncated plastomes owing to loss of photosynthetic and other genes. Nonasterid holoparasites from Balanophoraceae (Corynaea), Hydnoraceae (Hydnora) and Cytinaceae (Cytinus) were tested for the presence of plastid genes and a plastome. Using PCR, plastid 16S rDNA was successfully amplified and sequenced from the above three holoparasites. The sequence of Cytinus showed 121 single base substitutions relative to Nicotiana (8% of the molecule) whereas higher sequence divergence was observed in Hydnora and Corynaea (287 and 513 changes, respectively). Secondary structural models for these 16S rRNAs show that most changes are compensatory, thus suggesting they are functional. Probes constructed for 16S rDNA and for four plastid-encoded ribosomal protein genes (rps2, rps4, rps7 and rpl16) were used in Southern blots of digested genomic DNA from the three holoparasites. Positive hybridizations were obtained using each of the five probes only for Cytinus. For SmaI digests, all plastid gene probes hybridized to a common fragment ca. 20 kb in length in this species. Taken together, these data provide preliminary evidence suggestive of the retention of highly diverged and truncated plastid genome in Cytinus. The greater sequence divergence for 16S rDNA and the negative hybridization results for Hydnora and Corynaea suggests two possibilities: the loss of typically conserved elements of their plastomes or the complete absence of a plastome.     

Herbarium phylogenomics: Resolving the generic status of the enigmatic Pseudobartsia (Orobanchaceae)

The millions of herbarium specimens in collections around the world provide historical resources for phylogenomics and evolutionary studies. Many rare and endangered species exist only as historical specimens. Here, we report a case study of the monotypic Pseudobartsia yunnanensis D. Y. Hong (=Pseudobartsia glandulosa[Bentham] W. B. Yu & D. Z. Li: Orobanchaceae) known from a single Chinese collection taken in 1940. We obtained genomic data of Pseudobartsia glandulosa using high-throughput short-read sequencing, and then assembled a complete chloroplast genome and nuclear ribosome DNA region in this study. We found that the newly assembled three plastid DNA regions (atpB-rbcL, rpl16, and trnS-G) and nuclear ribosomal internal transcribed spacer (nrITS) of Pseudobartsia glandulosa were more than 99.98% similar to published sequences obtained by target sequencing. Phylogenies of Orobanchaceae using 30 plastomes (including 10 new plastomes), using both supermatrix and multispecies coalescent approaches following a novel plastid phylogenomic workflow, recovered seven recognized tribes and two unranked groups, both of which were proposed as new tribes, that is, Brandisieae and Pterygielleae. Within Pterygielleae, all analyses strongly supported Xizangia D. Y. Hong as the first diverging genus, with Pseudobartsia D. Y. Hong as sister to Pterygiella Oliver + Phtheirospermum Bunge (excluding Phtheirospermum japonicum [Thunberg] Kanitz); this supports reinstatement of Pseudobartsia and Xizangia. Although elements of Buchnereae-Cymbarieae-Orobancheae and Brandisieae-Pterygielleae-Rhinantheae showed incongruence among gene trees, the topology of the supermatrix tree was congruent with the majority of gene trees and functional-group trees. Therefore, most plastid genes are evolving as a linkage group, allowing the supermatrix tree approach to yield internally consistent phylogenies for Orobanchaceae.