Transcriptome map of plant mitochondria reveals islands of unexpected transcribed regions

Evolution of cox2 introns in angiosperm mitochondria and efficient splicing of an elongated cox2i691 intron

In angiosperms, the mitochondrial cox2 gene harbors up to two introns, commonly referred to as cox2i373 and cox2i691. We studied the cox2 from 222 fully-sequenced mitogenomes from 30 angiosperm orders and analyzed the evolution of their introns. Unlike cox2i373, cox2i691 shows a distribution among plants that is shaped by frequent intron loss events driven by localized retroprocessing. In addition, cox2i691 exhibits sporadic elongations, frequently in domain IV of introns. Such elongations are poorly related to repeat content and two of them showed the presence of LINE transposons, suggesting that increasing intron size is very likely due to nuclear intracelular DNA transfer followed by incorporation into the mitochondrial DNA. Surprisingly, we found that cox2i691 is erroneously annotated as absent in 30 mitogenomes deposited in public databases. Although each of the cox2 introns is ∼1.5 kb in length, a cox2i691 of 4.2 kb has been reported in Acacia ligulata (Fabaceae). It is still unclear whether its unusual length is due to a trans-splicing arrangement or the loss of functionality of the interrupted cox2. Through analyzing short-read RNA sequencing of Acacia with a multi-step computational strategy, we found that the Acacia cox2 is functional and its long intron is spliced in cis in a very efficient manner despite its length.

Transcriptional Landscape and Splicing Efficiency in Arabidopsis Mitochondria

Plant mitochondrial transcription is initiated from multiple promoters without an apparent motif, which precludes their identification in other species based on sequence comparisons. Even though coding regions take up only a small fraction of plant mitochondrial genomes, deep RNAseq studies uncovered that these genomes are fully or nearly fully transcribed with significantly different RNA read depth across the genome. Transcriptomic analysis can be a powerful tool to understand the transcription process in diverse angiosperms, including the identification of potential promoters and co-transcribed genes or to study the efficiency of intron splicing. In this work, we analyzed the transcriptional landscape of the Arabidopsis mitochondrial genome (mtDNA) based on large-scale RNA sequencing data to evaluate the use of RNAseq to study those aspects of the transcription process. We found that about 98% of the Arabidopsis mtDNA is transcribed with highly different RNA read depth, which was elevated in known genes. The location of a sharp increase in RNA read depth upstream of genes matched the experimentally identified promoters. The continuously high RNA read depth across two adjacent genes agreed with the known co-transcribed units in Arabidopsis mitochondria. Most intron-containing genes showed a high splicing efficiency with no differences between cis and trans-spliced introns or between genes with distinct splicing mechanisms. Deep RNAseq analyses of diverse plant species will be valuable to recognize general and lineage-specific characteristics related to the mitochondrial transcription process.

Horizontal gene transfers dominate the functional mitochondrial gene space of a holoparasitic plant

Although horizontal gene transfer (HGT) is common in angiosperm mitochondrial DNAs (mtDNAs), few cases of functional foreign genes have been identified. The one outstanding candidate for large-scale functional HGT is the holoparasite Lophophytum mirabile, whose mtDNA has lost most native genes but contains intact foreign homologs acquired from legume host plants. To investigate the extent to which this situation results from functional replacement of native by foreign genes, functional mitochondrial gene transfer to the nucleus, and/or loss of mitochondrial biochemical function in the context of extreme parasitism, we examined the Lophophytum mitochondrial and nuclear transcriptomes by deep paired-end RNA sequencing. Most foreign mitochondrial genes in Lophophytum are highly transcribed, accurately spliced, and efficiently RNA edited. By contrast, we found no evidence for functional gene transfer to the nucleus or loss of mitochondrial functions in Lophophytum. Many functional replacements occurred via the physical replacement of native genes by foreign genes. Some of these events probably occurred as the final act of HGT itself. Lophophytum mtDNA has experienced an unprecedented level of functional replacement of native genes by foreign copies. This raises important questions concerning population-genetic and molecular regimes that underlie such a high level of foreign gene takeover.

Non-coding RNAs and plant male sterility: current knowledge and future prospects

Latest outcomes assign functional role to non-coding (nc) RNA molecules in regulatory networks that confer male sterility to plants. Male sterility in plants offers great opportunity for improving crop performance through application of hybrid technology. In this respect, cytoplasmic male sterility (CMS) and sterility induced by photoperiod (PGMS)/temperature (TGMS) have greatly facilitated development of high-yielding hybrids in crops. Participation of non-coding (nc) RNA molecules in plant reproductive development is increasingly becoming evident. Recent breakthroughs in rice definitively associate ncRNAs with PGMS and TGMS. In case of CMS, the exact mechanism through which the mitochondrial ORFs exert influence on the development of male gametophyte remains obscure in several crops. High-throughput sequencing has enabled genome-wide discovery and validation of these regulatory molecules and their target genes, describing their potential roles performed in relation to CMS. Discovery of ncRNA localized in plant mtDNA with its possible implication in CMS induction is intriguing in this respect. Still, conclusive evidences linking ncRNA with CMS phenotypes are currently unavailable, demanding complementing genetic approaches like transgenics to substantiate the preliminary findings. Here, we review the recent literature on the contribution of ncRNAs in conferring male sterility to plants, with an emphasis on microRNAs. Also, we present a perspective on improved understanding about ncRNA-mediated regulatory pathways that control male sterility in plants. A refined understanding of plant male sterility would strengthen crop hybrid industry to deliver hybrids with improved performance.

The mitochondrial genome of the terrestrial carnivorous plant Utricularia reniformis (Lentibulariaceae): Structure, comparative analysis and evolutionary landmarks

The carnivorous plants of the family Lentibulariaceae have attained recent attention not only because of their interesting lifestyle, but also because of their dynamic nuclear genome size. Lentibulariaceae genomes span an order of magnitude and include species with the smallest genomes in angiosperms, making them a powerful system to study the mechanisms of genome expansion and contraction. However, little is known about mitochondrial DNA (mtDNA) sequences of this family, and the evolutionary forces that shape this organellar genome. Here we report the sequencing and assembly of the complete mtDNA from the endemic terrestrial Brazilian species Utricularia reniformis. The 857,234bp master circle mitochondrial genome encodes 70 transcriptionaly active genes (42 protein-coding, 25 tRNAs and 3 rRNAs), covering up to 7% of the mtDNA. A ltrA-like protein related to splicing and mobility and a LAGLIDADG homing endonuclease have been identified in intronic regions, suggesting particular mechanisms of genome maintenance. RNA-seq analysis identified properties with putative diverse and important roles in genome regulation and evolution: 1) 672kbp (78%) of the mtDNA is covered by full-length reads; 2) most of the 243kbp intergenic regions exhibit transcripts; and 3) at least 69 novel RNA editing sites in the protein-coding genes. Additional genomic features are hypothetical ORFs (48%), chloroplast insertions, including truncated plastid genes that have been lost from the chloroplast DNA (5%), repeats (5%), relics of transposable elements mostly related to LTR retrotransposons (5%), and truncated mitovirus sequences (0.4%). Phylogenetic analysis based on 32 different Lamiales mitochondrial genomes corroborate that Lentibulariaceae is a monophyletic group. In summary, the U. reniformis mtDNA represents the eighth largest plant mtDNA described to date, shedding light on the genomic trends and evolutionary characteristics and phylogenetic history of the family Lentibulariaceae.

Non-coding RNA may be associated with cytoplasmic male sterility in Silene vulgaris

Cytoplasmic male sterility (CMS) is a widespread phenomenon in flowering plants caused by mitochondrial (mt) genes. CMS genes typically encode novel proteins that interfere with mt functions and can be silenced by nuclear fertility-restorer genes. Although the molecular basis of CMS is well established in a number of crop systems, our understanding of it in natural populations is far more limited. To identify CMS genes in a gynodioecious plant, Silene vulgaris, we constructed mt transcriptomes and compared transcript levels and RNA editing patterns in floral bud tissue from female and hermaphrodite full siblings. The transcriptomes from female and hermaphrodite individuals were very similar overall with respect to variation in levels of transcript abundance across the genome, the extent of RNA editing, and the order in which RNA editing and intron splicing events occurred. We found only a single genomic region that was highly overexpressed and differentially edited in females relative to hermaphrodites. This region is not located near any other transcribed elements and lacks an open-reading frame (ORF) of even moderate size. To our knowledge, this transcript would represent the first non-coding mt RNA associated with CMS in plants and is, therefore, an important target for future functional validation studies.

Cytoplasmic male sterility (CMS) in hybrid breeding in field crops

A comprehensive understanding of CMS/Rf system enabled by modern omics tools and technologies considerably improves our ability to harness hybrid technology for enhancing the productivity of field crops. Harnessing hybrid vigor or heterosis is a promising approach to tackle the current challenge of sustaining enhanced yield gains of field crops. In the context, cytoplasmic male sterility (CMS) owing to its heritable nature to manifest non-functional male gametophyte remains a cost-effective system to promote efficient hybrid seed production. The phenomenon of CMS stems from a complex interplay between maternally-inherited (mitochondrion) and bi-parental (nucleus) genomic elements. In recent years, attempts aimed to comprehend the sterility-inducing factors (orfs) and corresponding fertility determinants (Rf) in plants have greatly increased our access to candidate genomic segments and the cloned genes. To this end, novel insights obtained by applying state-of-the-art omics platforms have substantially enriched our understanding of cytoplasmic-nuclear communication. Concomitantly, molecular tools including DNA markers have been implicated in crop hybrid breeding in order to greatly expedite the progress. Here, we review the status of diverse sterility-inducing cytoplasms and associated Rf factors reported across different field crops along with exploring opportunities for integrating modern omics tools with CMS-based hybrid breeding.

Mitochondrial ORF79 levels determine pollen abortion in cytoplasmic male sterile rice

Cytoplasmic male sterility (CMS) is an important agricultural trait characterized by lack of functional pollen, and caused by ectopic and defective mitochondrial gene expression. The pollen function in CMS plants is restored by the presence of nuclear-encoded restorer of fertility (Rf) genes. Previously, we cloned Rf2, which restores the fertility of Lead Rice (LD)-type CMS rice. However, neither the function of Rf2 nor the identity of the mitochondrial gene causing CMS has been determined in LD-CMS rice. Here, we show that the mitochondrial gene orf79 acts as a CMS-associated gene in LD-CMS rice, similar to its role in BT-CMS rice originating from Chinsurah Boro II, and Rf2 weakly restores fertility in BT-CMS rice. We also show that RF2 promotes degradation of atp6-orf79 RNA in a different manner from that of RF1, which is the Rf gene product in BT-CMS rice. The amount of ORF79 protein in LD-CMS rice was one-twentieth of the amount in BT-CMS rice. The difference in ORF79 protein levels probably accounts for the mild and severe pollen defects in LD-CMS and BT-CMS rice, respectively. In the presence of Rf2, accumulation of ORF79 was reduced to almost zero and 25% in LD-CMS and BT-CMS rice, respectively, which probably accounts for the complete and weak fertility restoration abilities of Rf2 in LD-CMS and BT-CMS rice, respectively. These observations indicate that the amount of ORF79 influences the pollen fertility in two strains of rice in which CMS is induced by orf79.

High transcript abundance, RNA editing, and small RNAs in intergenic regions within the massive mitochondrial genome of the angiosperm Silene noctiflora

BACKGROUND

Species within the angiosperm genus Silene contain the largest mitochondrial genomes ever identified. The enormity of these genomes (up to 11 Mb in size) appears to be the result of increased non-coding DNA, which represents >99 % of the genome content. These genomes are also fragmented into dozens of circular-mapping chromosomes, some of which contain no identifiable genes, raising questions about if and how these 'empty' chromosomes are maintained by selection. To assess the possibility that they contain novel and unannotated functional elements, we have performed RNA-seq to analyze the mitochondrial transcriptome of Silene noctiflora.

RESULTS

We identified regions of high transcript abundance in almost every chromosome in the mitochondrial genome including those that lack any annotated genes. In some cases, these transcribed regions exhibited higher expression levels than some core mitochondrial protein-coding genes. We also identified RNA editing sites throughout the genome, including 97 sites that were outside of protein-coding gene sequences and found in pseudogenes, introns, UTRs, and transcribed intergenic regions. Unlike in protein-coding sequences, however, most of these RNA editing sites were only edited at intermediate frequencies. Finally, analysis of mitochondrial small RNAs indicated that most were likely degradation products from longer transcripts, but we did identify candidates for functional small RNAs that mapped to intergenic regions and were not associated with longer RNA transcripts.

CONCLUSIONS

Our findings demonstrate transcriptional activity in many localized regions within the extensive intergenic sequence content in the S. noctiflora mitochondrial genome, supporting the possibility that the genome contains previously unidentified functional elements. However, transcription by itself is not proof of functional importance, and we discuss evidence that some of the observed transcription and post-transcriptional modifications are non-adaptive. Therefore, further investigations are required to determine whether any of the identified transcribed regions have played a functional role in the proliferation and maintenance of the enormous non-coding regions in Silene mitochondrial genomes.

Comparative analysis of 11 Brassicales mitochondrial genomes and the mitochondrial transcriptome of Brassica oleracea

To elucidate the evolution of mitochondrial genomic diversity within a single order of angiosperms, we sequenced seven Brassicales genomes and the transcriptome of Brassica oleracea. In the common ancestor of Brassicaceae, several genes of known function were lost and the ccmFN gene was split into two independent genes, which also coincides with a trend of genome reduction towards the smallest sequenced angiosperm genomes of Brassica. For most ORFs of unknown function, the lack of conservation throughout Brassicales and the generally low expression and absence of RNA editing in B. oleracea argue against functionality. However, two chimeric ORFs were expressed and edited in B. oleracea, suggesting a potential role in cytoplasmic male sterility in certain nuclear backgrounds. These results demonstrate how frequent shifts in size, structure, and content of plant mitochondrial genomes can occur over short evolutionary time scales.

RNA metabolism in plant mitochondria

Mitochondria are essential for the eukaryotic cell and are derived from the endosymbiosis of an α-proteobacterial ancestor. Compared to other eukaryotes, RNA metabolism in plant mitochondria is complex and combines bacterial-like traits with novel features that evolved in the host cell. These complex RNA processes are regulated by families of nucleus-encoded RNA-binding proteins. Transcription is particularly relaxed and is initiated from multiple promoters covering the entire genome. The variety of RNA precursors accumulating in mitochondria highlights the importance of post-transcriptional processes to determine the size and abundance of transcripts. Here we review RNA metabolism in plant mitochondria, from RNA transcription to translation, with a special focus on their unique features that are controlled by trans-factors.

The genome and transcriptome of perennial ryegrass mitochondria

Background

Perennial ryegrass (Lolium perenne L.) is one of the most important forage and turf grass species of temperate regions worldwide. Its mitochondrial genome is inherited maternally and contains genes that can influence traits of agricultural importance. Moreover, the DNA sequence of mitochondrial genomes has been established and compared for a large number of species in order to characterize evolutionary relationships. Therefore, it is crucial to understand the organization of the mitochondrial genome and how it varies between and within species. Here, we report the first de novo assembly and annotation of the complete mitochondrial genome from perennial ryegrass.

Results

Intact mitochondria from perennial ryegrass leaves were isolated and used for mtDNA extraction. The mitochondrial genome was sequenced to a 167-fold coverage using the Roche 454 GS-FLX Titanium platform, and assembled into a circular master molecule of 678,580 bp. A total of 34 proteins, 14 tRNAs and 3 rRNAs are encoded by the mitochondrial genome, giving a total gene space of 48,723 bp (7.2%). Moreover, we identified 149 open reading frames larger than 300 bp and covering 67,410 bp (9.93%), 250 SSRs, 29 tandem repeats, 5 pairs of large repeats, and 96 pairs of short inverted repeats. The genes encoding subunits of the respiratory complexes – nad1 to nad9, cob, cox1 to cox3 and atp1 to atp9 – all showed high expression levels both in absolute numbers and after normalization.

Conclusions

The circular master molecule of the mitochondrial genome from perennial ryegrass presented here constitutes an important tool for future attempts to compare mitochondrial genomes within and between grass species. Our results also demonstrate that mitochondria of perennial ryegrass contain genes crucial for energy production that are well conserved in the mitochondrial genome of monocotyledonous species. The expression analysis gave us first insights into the transcriptome of these mitochondrial genes in perennial ryegrass.

Inhibition of a basal transcription factor 3-like gene Osj10gBTF3 in rice results in significant plant miniaturization and typical pollen abortion

BTF3, which was originally recognized as a basal transcription factor, has been known to be involved in transcription initiation, translational regulation and protein localization in many eukaryotic organisms. However, its function remains largely unknown in plant species. In the present study, we analyzed a BTF3-related sequence in Oryza sativa L. subsp. japonica, which shares the conserved domain of a nascent polypeptide-associated complex with human BTF3, and was referred to as Osj10gBTF3. The expression of Osj10gBTF3 was primarily constitutive and generally modulated by salt, high temperature and exogenous phytohormone stress. The Osj10gBTF3::EGFP (enhanced green fluorescence protein) fusion protein was localized in both the nucleus and cytoplasmic membrane system. Inhibition of Osj10gBTF3 led to significant morphological changes in all detected tissues and organs, with a reduced size of between 25% and 52%. Furthermore, the pollen that developed was completely sterile, which was correlated with the altered expression of two Rf (fertility restorer)-like genes that encode pentatricopeptide repeat-containing proteins OsPPR676 and OsPPR920, translational initiation factors OseIF3e and OseIF3h, and the heat shock protein OsHSP82. These findings were verified through a yeast two-hybrid assay using a Nipponbare callus cDNA library as bait followed by the reverse transcription-PCR analysis of total leaf or anther RNAs. Our demonstration of the important role of Osj10gBTF3 in rice growth and development provides new insights showing that more complex regulatory functions are associated with BTF3 in plants.

Mapping of wheat mitochondrial mRNA termini and comparison with breakpoints in DNA homology among plants

Mitochondrial DNA rearrangements occur very frequently in flowering plants and when close to genes there must be concomitant acquisition of new regulatory cis-elements. To explore whether there might be limits to such DNA shuffling, we have mapped the termini of mitochondrial mRNAs in wheat, a monocot, and compared them to the known positions for counterpart genes in the eudicot Arabidopsis. Nine genes share homologous 3' UTRs over their full-length and for six of them, the termini map very close to the site of wheat/Arabidopsis DNA rearrangements. Only one such case was seen for comparisons of 5' UTRs, and the 5' ends of mRNAs are typically more heterogeneous than 3' termini. Approximately half of the thirty-one wheat mitochondrial transcriptional units are preceded by CRTA promoter-like motifs, and of the potential stem-loop or tRNA-like structures identified as candidate RNA processing/stability signals near the 5' or 3' ends, several are shared with Arabidopsis. Comparison of the mitochondrial gene flanking sequences from normal fertile wheat (Triticum aestivum) with those of Aegilops kotschyi which is the source of mitochondria present in K-type cytoplasmic male sterile wheat, revealed six cases where mRNAs are precluded from sharing full-length homologous UTRs because of genomic reorganization events, and the presence of short repeats located at the sites of discontinuity points to a reciprocal recombination-mediated mode of rearrangement.

Potential role of tRNAs in wheat and Lolium mitochondrial rps7 transcript processing

The wheat mitochondrial gene for ribosomal protein S7 exhibits multiple transcripts that share the same 3′ terminus but range in overall length from 3.4 to 0.7 kb because of 5′end-maturation events. The longest detectable precursor RNA maps precisely to the 3′ end of a chloroplast-origin tRNA-Phe gene, consistent with it providing signals for endonucleolytic cleavage. Steady-state levels of precursor RNAs were seen to be lower in seedlings than in germinating embryos, although the degree of editing within untranslated regions (UTRs) was higher in seedlings. In another grass, Lolium multiflorum Lam., rps7 displays transcripts of 1.3 and 0.7 kb, and although the distal 5′ UTRs are unrelated in sequence to those of wheat, the 5′ terminus of the longer transcript also maps to a tRNA gene, in this case the native mitochondrial-type tRNA-Ser. Our findings illustrate the plasticity of plant mitochondrial transcriptional units and the recruitment of chloroplast-origin sequences for the expression of mitochondrial genes.

A complete sequence and transcriptomic analyses of date palm (Phoenix dactylifera L.) mitochondrial genome

Based on next-generation sequencing data, we assembled the mitochondrial (mt) genome of date palm (Phoenix dactylifera L.) into a circular molecule of 715,001 bp in length. The mt genome of P. dactylifera encodes 38 proteins, 30 tRNAs, and 3 ribosomal RNAs, which constitute a gene content of 6.5% (46,770 bp) over the full length. The rest, 93.5% of the genome sequence, is comprised of cp (chloroplast)-derived (10.3% with respect to the whole genome length) and non-coding sequences. In the non-coding regions, there are 0.33% tandem and 2.3% long repeats. Our transcriptomic data from eight tissues (root, seed, bud, fruit, green leaf, yellow leaf, female flower, and male flower) showed higher gene expression levels in male flower, root, bud, and female flower, as compared to four other tissues. We identified 120 potential SNPs among three date palm cultivars (Khalas, Fahal, and Sukry), and successfully found seven SNPs in the coding sequences. A phylogenetic analysis, based on 22 conserved genes of 15 representative plant mitochondria, showed that P. dactylifera positions at the root of all sequenced monocot mt genomes. In addition, consistent with previous discoveries, there are three co-transcribed gene clusters–18S-5S rRNA, rps3-rpl16 and nad3-rps12–in P. dactylifera, which are highly conserved among all known mitochondrial genomes of angiosperms.