Extreme plastomes in holoparasitic Balanophoraceae are not the norm.

BACKGROUND: Balanophoraceae plastomes are known for their highly condensed and re-arranged nature alongside the most extreme nucleotide compositional bias known to date, culminating in two independent reconfigurations of their genetic code. Currently, a large portion of the Balanophoraceae diversity remains unexplored, hindering, among others, evolutionary pattern recognition. Here, we explored newly sequenced plastomes of Sarcophyte sanguinea and Thonningia sanguinea. The reconstructed plastomes were analyzed using various methods of comparative genomics based on a representative taxon sampling. RESULTS: Sarcophyte, recovered sister to the other sampled Balanophoraceae s. str., has plastomes up to 50% larger than those currently published. Its gene set contains five genes lost in any other species, including matK. Five cis-spliced introns are maintained. In contrast, the Thonningia plastome is similarly reduced to published Balanophoraceae and retains only a single cis-spliced intron. Its protein-coding genes show a more biased codon usage compared to Sarcophyte, with an accumulation of in-frame TAG stop codons. Structural plastome comparison revealed multiple, previously unknown, structural rearrangements within Balanophoraceae. CONCLUSIONS: For the "minimal plastomes" of Thonningia, we propose a genetic code change identical to sister genus Balanophora. Sarcophyte however differs drastically from our current understanding on Balanophoraceae plastomes. With a less-extreme nucleotide composition, there is no evidence for an altered genetic code. Using comparative genomics, we identified a hotspot for plastome reconfiguration in Balanophoraceae. Based on previously published and newly identified structural reconfigurations, we propose an updated model of evolutionary plastome trajectories for Balanophoraceae, illustrating a much greater plastome diversity than previously known.

Cytonuclear coevolution in a holoparasitic plant with highly disparate organellar genomes.

KEY MESSAGE: Contrasting substitution rates in the organellar genomes of Lophophytum agree with the DNA repair, replication, and recombination gene content. Plastid and nuclear genes whose products form multisubunit complexes co-evolve. The organellar genomes of the holoparasitic plant Lophophytum (Balanophoraceae) show disparate evolution. In the plastid, the genome has been severely reduced and presents a > 85% AT content, while in the mitochondria most protein-coding genes have been replaced by homologs acquired by horizontal gene transfer (HGT) from their hosts (Fabaceae). Both genomes carry genes whose products form multisubunit complexes with those of nuclear genes, creating a possible hotspot of cytonuclear coevolution. In this study, we assessed the evolutionary rates of plastid, mitochondrial and nuclear genes, and their impact on cytonuclear evolution of genes involved in multisubunit complexes related to lipid biosynthesis and proteolysis in the plastid and those in charge of the oxidative phosphorylation in the mitochondria. Genes from the plastid and the mitochondria (both native and foreign) of Lophophytum showed extremely high and ordinary substitution rates, respectively. These results agree with the biased loss of plastid-targeted proteins involved in angiosperm organellar repair, replication, and recombination machinery. Consistent with the high rate of evolution of plastid genes, nuclear-encoded subunits of plastid complexes showed disproportionate increases in non-synonymous substitution rates, while those of the mitochondrial complexes did not show different rates than the control (i.e. non-organellar nuclear genes). Moreover, the increases in the nuclear-encoded subunits of plastid complexes were positively correlated with the level of physical interaction they possess with the plastid-encoded ones. Overall, these results suggest that a structurally-mediated compensatory factor may be driving plastid-nuclear coevolution in Lophophytum, and that mito-nuclear coevolution was not altered by HGT.

Native and foreign mitochondrial and nuclear encoded proteins conform the OXPHOS complexes of a holoparasitic plant.

The intimate contact between the holoparasitic plant Lophophytum mirabile (Balanophoraceae) and its host plant (Fabaceae) facilitates the exchange of genetic information, increasing the frequency of horizontal gene transfer (HGT). Lophophytum stands out because it acquired a large number of mitochondrial genes (greater than 20) from its legume host that replaced the majority of the native homologs. These foreign genes code for proteins that form multisubunit enzyme complexes, such as those in the oxidative phosphorylation system (OXPHOS) and cytochrome c maturation (ccm) system, together with dozens of nuclear-encoded subunits. However, the existence and the origin of the nuclear subunits that form the major part of the OXPHOS and ccm system in Lophophytum remain unknown. It was proposed that nuclear-encoding genes whose products interact with foreign mitochondrial proteins are also foreign, minimizing the incompatibilities that could arise in the assembly and functioning of these multiprotein complexes. We identified a nearly complete set of OXPHOS and ccm system subunits evolving under selective constraints in the transcriptome of Lophophytum, indicating that OXPHOS is functional and resembles that of free-living angiosperms. Maximum Likelihood phylogenetic analyses revealed a single case of HGT in the nuclear genes, which results in mosaic OXPHOS and ccm system in Lophophytum. These observations raise new questions about the evolution and physiology of this parasitic plant. A putative case of cooperation between two foreign (one mitochondrial and one nuclear) genes is presented.

The minicircular and extremely heteroplasmic mitogenome of the holoparasitic plant Rhopalocnemis phalloides.

The plastid and nuclear genomes of parasitic plants exhibit deeply altered architectures,1-13 whereas the few examined mitogenomes range from deeply altered to conventional.14-20 To provide further insight on mitogenome evolution in parasitic plants, we report the highly modified mitogenome of Rhopalocnemis phalloides, a holoparasite in Balanophoraceae. Its mitogenome is uniquely arranged in 21 minicircular chromosomes that vary in size from 4,949 to 7,861 bp, with a total length of only 130,713 bp. All chromosomes share an identical 896 bp conserved region, with a large stem-loop that acts as the origin of replication, flanked on each side by hypervariable and semi-conserved regions. Similar minicircular structures with shared and unique regions have been observed in parasitic animals and free-living protists,21-24 suggesting convergent structural evolution. Southern blots confirm both the minicircular structure and the replication origin of the mitochondrial chromosomes. PacBio reads provide evidence for chromosome recombination and rolling-circle replication for the R. phalloides mitogenome. Despite its small size, the mitogenome harbors a typical set of genes and introns within the unique regions of each chromosome, yet introns are the smallest among seed plants and ferns. The mitogenome also exhibits extreme heteroplasmy, predominantly involving short indels and more complex variants, many of which cause potential loss-of-function mutations for some gene copies. All heteroplasmic variants are transcribed, and functional and nonfunctional protein-coding variants are spliced and RNA edited. Our findings offer a unique perspective into how mitogenomes of parasitic plants can be deeply altered and shed light on plant mitogenome replication.

Plastomes in the holoparasitic family Balanophoraceae: Extremely high AT content, severe gene content reduction, and two independent genetic code changes.

The transition to a heterotrophic lifestyle in angiosperms is characterized by convergent evolutionary changes. Plastid genome remodeling includes dramatic functional and physical reductions with the highest degrees observed in fully heterotrophic plants. Genes related to photosynthesis are generally absent or pseudogenized, while a few genes related to other metabolic processes that take place within the plastid are almost invariably maintained. The family Balanophoraceae consists of root holoparasites that present reduced plastid genomes with an extraordinarily elevated AT content and the single genetic code change ever documented in land plant plastomes (the stop codon TAG now codes for tryptophan). Here, we studied the plastomes of Lophophytum leandri and Ombrophytum subterraneum (Balanophoraceae) that showed the remarkable absence of the gene trnE, a highly biased nucleotide composition, and an independent genetic code change (the standard stop codon TGA codes for tryptophan). This is the second genetic code change identified in land plant plastomes. Analysis of the transcriptome of Lophophytum indicated that the entire C5 pathway typical of plants is conserved despite the lack of trnE in its plastome. A hypothetical model of plastome evolution in the Balanophoraceae is presented.

Multichromosomal structure and foreign tracts in the Ombrophytum subterraneum (Balanophoraceae) mitochondrial genome.

Horizontal gene transfer (HGT) is frequent in parasitic plant mitochondria as a result of vascular connections established in host-parasite relationships. Recent studies of the holoparasitic plant Lophophytum mirabile (Balanophoraceae) revealed the unprecedented acquisition of a large amount of mitochondrial sequences from its legume host. We focused on a close relative, the generalist holoparasite Ombrophytum subterraneum, to examine the incidence of HGT events in the mitochondrial genome (mtDNA). The mtDNA of O. subterraneum assembles into 54 circular chromosomes, only 34 of which contain the 51 full-length coding regions. Numerous foreign tracts (totaling almost 100 kb, ~ 14% of the mtDNA), including 12 intact genes, were acquired by HGT from the Asteraceae hosts. Nine chromosomes concentrate most of those regions and eight are almost entirely foreign. Native homologs of each foreign gene coexist in the mtDNA and are potentially functional. A large proportion of shorter regions were related to the Fabaceae (a total of ~ 110 kb, 15.4%), some of which were shared with L. mirabile. We also found evidence of foreign sequences donated by angiosperm lineages not reported as hosts (Apocynaceae, Euphorbiaceae, Lamiaceae, and Malvales). We propose an evolutionary hypothesis that involves ancient transfers from legume hosts in the common ancestor of Ombrophytum and Lophophytum followed by more recent transfer events in L. mirabile. Besides, the O. subterraneum mtDNA was also subjected to additional HGT events from diverse angiosperm lineages, including large and recent transfers from the Asteraceae, and also from Lamiaceae.

Novel genetic code and record-setting AT-richness in the highly reduced plastid genome of the holoparasitic plant Balanophora.

Plastid genomes (plastomes) vary enormously in size and gene content among the many lineages of nonphotosynthetic plants, but key lineages remain unexplored. We therefore investigated plastome sequence and expression in the holoparasitic and morphologically bizarre Balanophoraceae. The two Balanophora plastomes examined are remarkable, exhibiting features rarely if ever seen before in plastomes or in any other genomes. At 15.5 kb in size and with only 19 genes, they are among the most reduced plastomes known. They have no tRNA genes for protein synthesis, a trait found in only three other plastid lineages, and thus Balanophora plastids must import all tRNAs needed for translation. Balanophora plastomes are exceptionally compact, with numerous overlapping genes, highly reduced spacers, loss of all cis-spliced introns, and shrunken protein genes. With A+T contents of 87.8% and 88.4%, the Balanophora genomes are the most AT-rich genomes known save for a single mitochondrial genome that is merely bloated with AT-rich spacer DNA. Most plastid protein genes in Balanophora consist of ≥90% AT, with several between 95% and 98% AT, resulting in the most biased codon usage in any genome described to date. A potential consequence of its radical compositional evolution is the novel genetic code used by Balanophora plastids, in which TAG has been reassigned from stop to tryptophan. Despite its many exceptional properties, the Balanophora plastome must be functional because all examined genes are transcribed, its only intron is correctly trans-spliced, and its protein genes, although highly divergent, are evolving under various degrees of selective constraint.

Genome-scale transfer of mitochondrial DNA from legume hosts to the holoparasite Lophophytum mirabile (Balanophoraceae).

Angiosperm mitochondrial horizontal gene transfer (HGT) has been widely reported during the past decades. With a few exceptions, foreign sequences are mitochondrial genes or intronic regions from other plants, indicating that HGT has played a major role in shaping mitochondrial genome evolution. Host-parasite relationships are a valuable system to study this phenomenon due to the high frequency of HGT. In particular, the interaction between mimosoid legumes and holoparasites of the genus Lophophytum represents an outstanding opportunity to discern HGT events. The mitochondrial genome of the holoparasite L. mirabile has remarkable properties, the most extraordinary of which is the presence of 34 out of 43 mitochondrial protein genes acquired from its legume host, with the stunning replacement of up to 26 native homologs. However, the origin of the intergenic sequences that represent the majority (>90%) of the L. mirabile mtDNA remains largely unknown. The lack of mitochondrial sequences available from the donor angiosperm lineage (mimosoid legumes) precluded a large-scale evolutionary study. We sequenced and assembled the mitochondrial genome of the mimosoid Acacia ligulata and performed genome wide comparisons with L. mirabile. The A. ligulata mitochondrial genome is almost 700 kb in size, encoding 60 genes. About 60% of the L. mirabile mtDNA had greatest affinity to members of the family Fabaceae (∼49% to mimosoids in particular) with an average sequence identity of ∼96%, including genes but mostly intergenic regions. These findings strengthen the mitochondrial fusion compatibility model for angiosperm mitochondrion-to-mitochondrion HGT.

Rate heterogeneity in six protein-coding genes from the holoparasite Balanophora (Balanophoraceae) and other taxa of Santalales.

BACKGROUND AND AIMS: The holoparasitic flowering plant Balanophora displays extreme floral reduction and was previously found to have enormous rate acceleration in the nuclear 18S rDNA region. So far, it remains unclear whether non-ribosomal, protein-coding genes of Balanophora also evolve in an accelerated fashion and whether the genes with high substitution rates retain their functionality. To tackle these issues, six different genes were sequenced from two Balanophora species and their rate variation and expression patterns were examined. METHODS: Sequences including nuclear PI, euAP3, TM6, LFY and RPB2 and mitochondrial matR were determined from two Balanophora spp. and compared with selected hemiparasitic species of Santalales and autotrophic core eudicots. Gene expression was detected for the six protein-coding genes and the expression patterns of the three B-class genes (PI, AP3 and TM6) were further examined across different organs of B. laxiflora using RT-PCR. KEY RESULTS: Balanophora mitochondrial matR is highly accelerated in both nonsynonymous (d(N)) and synonymous (d(S)) substitution rates, whereas the rate variation of nuclear genes LFY, PI, euAP3, TM6 and RPB2 are less dramatic. Significant d(S) increases were detected in Balanophora PI, TM6, RPB2 and d(N) accelerations in euAP3. All of the protein-coding genes are expressed in inflorescences, indicative of their functionality. PI is restrictively expressed in tepals, synandria and floral bracts, whereas AP3 and TM6 are widely expressed in both male and female inflorescences. CONCLUSIONS: Despite the observation that rates of sequence evolution are generally higher in Balanophora than in hemiparasitic species of Santalales and autotrophic core eudicots, the five nuclear protein-coding genes are functional and are evolving at a much slower rate than 18S rDNA. The mechanism or mechanisms responsible for rapid sequence evolution and concomitant rate acceleration for 18S rDNA and matR are currently not well understood and require further study in Balanophora and other holoparasites.

Microsatellite markers for the endangered root holoparasite Dactylanthus taylorii (Balanophoraceae) from 454 pyrosequencing.

PREMISE OF THE STUDY: Microsatellite loci were isolated and developed as polymorphic markers for the New Zealand endemic root holoparasite Dactylanthus taylorii for use in population and conservation genetics studies. METHODS AND RESULTS: Shotgun 454 pyrosequencing was performed on genomic DNA pooled from three individuals of D. taylorii. From 61709 individual sequence reads, primers for 753 microsatellite loci were developed in silico and 72 of these were tested for consistent amplification and variability. Ten microsatellite loci were found to be polymorphic and consistently scorable when screened in 44 individuals from five geographically distant populations. The number of alleles per locus ranged from four to 16 with an average of 9.7, and average observed heterozygosity per locus was between 0.182 and 0.634. CONCLUSIONS: These polymorphic microsatellite markers establish an important resource for ongoing conservation initiatives and planned population genetic studies of D. taylorii.

Morphology and phylogenetics of two holoparasitic plants, Balanophora japonica and Balanophora yakushimensis (Balanophoraceae), and their hosts in Taiwan and Japan.

Balanophora japonica and B. yakushimensis are two putatively agamospermic taxa previously reported from southern Japan. Their inflorescences superficially represent those of B. laxiflora and B. fungosa. In this study we confirmed their presence in Taiwan by morphological and phylogenetic analysis using nuclear 18S rDNA and nrITS sequences with related taxa. B. japonica, B. yakushimensis, and B. laxiflora formed a well-supported clade that is distinct from other Balanophora. All three taxa also show considerable differences on morphological and nucleotide sequence differences, therefore the name of B. yakushimensis is retained. The results provide new insights on the intrageneric classification of Balanophora and suggest the positioning of female flowers should be down-weighted. We also successfully identify the hosts of B. japonica and B. yakushimensis by amplifying chloroplast matK sequences from the connected root tissues. The results showed that B. japonica parasitizes on Symplocos species, and that B. yakushimensis parasitizes on Distylium racemosum in Japan and Schima superba in Taiwan's population.

Two micro-scale protocols for the isolation of DNA from polysaccharide-rich plant tissue.

The high polysaccharide content of some plant species hinders the successful isolation of their DNA. As an alternative to the macro-extraction methods previously published for polysaccharide-rich plants, we present two techniques (STE/CTAB and HEPES/CTAB), which are performed in microcentrifuge tubes. These protocols are suitable for small amounts of silica gel-preserved plant tissue such as are commonly available from endangered plants. The critical step to remove polysaccharides was performing initial washes in either STE (0.25 M sucrose, 0.03 M Tris, 0.05 M EDTA) or HEPES (2% ß-mercaptoethanol, 0.2% PVP, 0.1 M HEPES, pH 8.0) buffer. Precipitating the DNA at room temperature with isopropanol also aided in decreasing polysaccharide co-precipitation. Of the two protocols we present the STE/CTAB method has the advantages of being more cost-effective and avoiding the use of the hazardous chemical ß-mercaptoethanol.