Trans splicing in Oenothera mitochondria: nad1 mRNAs are edited in exon and trans-splicing group II intron sequences

Mitochondrial genome annotation with MFannot: a critical analysis of gene identification and gene model prediction

Compared to nuclear genomes, mitochondrial genomes (mitogenomes) are small and usually code for only a few dozen genes. Still, identifying genes and their structure can be challenging and time-consuming. Even automated tools for mitochondrial genome annotation often require manual analysis and curation by skilled experts. The most difficult steps are (i) the structural modelling of intron-containing genes; (ii) the identification and delineation of Group I and II introns; and (iii) the identification of moderately conserved, non-coding RNA (ncRNA) genes specifying 5S rRNAs, tmRNAs and RNase P RNAs. Additional challenges arise through genetic code evolution which can redefine the translational identity of both start and stop codons, thus obscuring protein-coding genes. Further, RNA editing can render gene identification difficult, if not impossible, without additional RNA sequence data. Current automated mito- and plastid-genome annotators are limited as they are typically tailored to specific eukaryotic groups. The MFannot annotator we developed is unique in its applicability to a broad taxonomic scope, its accuracy in gene model inference, and its capabilities in intron identification and classification. The pipeline leverages curated profile Hidden Markov Models (HMMs), covariance (CMs) and ERPIN models to better capture evolutionarily conserved signatures in the primary sequence (HMMs and CMs) as well as secondary structure (CMs and ERPIN). Here we formally describe MFannot, which has been available as a web-accessible service (https://megasun.bch.umontreal.ca/apps/mfannot/) to the research community for nearly 16 years. Further, we report its performance on particularly intron-rich mitogenomes and describe ongoing and future developments.

Graph-based models of the Oenothera mitochondrial genome capture the enormous complexity of higher plant mitochondrial DNA organization

Plant mitochondrial genomes display an enormous structural complexity, as recombining repeat-pairs lead to the generation of various sub-genomic molecules, rendering these genomes extremely challenging to assemble. We present a novel bioinformatic data-processing pipeline called SAGBAC (Semi-Automated Graph-Based Assembly Curator) that identifies recombinogenic repeat-pairs and reconstructs plant mitochondrial genomes. SAGBAC processes assembly outputs and applies our novel ISEIS (Iterative Sequence Ends Identity Search) algorithm to obtain a graph-based visualization. We applied this approach to three mitochondrial genomes of evening primrose (Oenothera), a plant genus used for cytoplasmic genetics studies. All identified repeat pairs were found to be flanked by two alternative and unique sequence-contigs defining so-called 'double forks', resulting in four possible contig-repeat-contig combinations for each repeat pair. Based on the inferred structural models, the stoichiometry of the different contig-repeat-contig combinations was analyzed using Illumina mate-pair and PacBio RSII data. This uncovered a remarkable structural diversity of the three closely related mitochondrial genomes, as well as substantial phylogenetic variation of the underlying repeats. Our model allows predicting all recombination events and, thus, all possible sub-genomes. In future work, the proposed methodology may prove useful for the investigation of the sub-genome organization and dynamics in different tissues and at various developmental stages.

Plant mitochondria - past, present and future

The study of plant mitochondria started in earnest around 1950 with the first isolations of mitochondria from animal and plant tissues. The first 35 years were spent establishing the basic properties of plant mitochondria and plant respiration using biochemical and physiological approaches. A number of unique properties (compared to mammalian mitochondria) were observed: (i) the ability to oxidize malate, glycine and cytosolic NAD(P)H at high rates; (ii) the partial insensitivity to rotenone, which turned out to be due to the presence of a second NADH dehydrogenase on the inner surface of the inner mitochondrial membrane in addition to the classical Complex I NADH dehydrogenase; and (iii) the partial insensitivity to cyanide, which turned out to be due to an alternative oxidase, which is also located on the inner surface of the inner mitochondrial membrane, in addition to the classical Complex IV, cytochrome oxidase. With the appearance of molecular biology methods around 1985, followed by genomics, further unique properties were discovered: (iv) plant mitochondrial DNA (mtDNA) is 10-600 times larger than the mammalian mtDNA, yet it only contains approximately 50% more genes; (v) plant mtDNA has kept the standard genetic code, and it has a low divergence rate with respect to point mutations, but a high recombinatorial activity; (vi) mitochondrial mRNA maturation includes a uniquely complex set of activities for processing, splicing and editing (at hundreds of sites); (vii) recombination in mtDNA creates novel reading frames that can produce male sterility; and (viii) plant mitochondria have a large proteome with 2000-3000 different proteins containing many unique proteins such as 200-300 pentatricopeptide repeat proteins. We describe the present and fairly detailed picture of the structure and function of plant mitochondria and how the unique properties make their metabolism more flexible allowing them to be involved in many diverse processes in the plant cell, such as photosynthesis, photorespiration, CAM and C4 metabolism, heat production, temperature control, stress resistance mechanisms, programmed cell death and genomic evolution. However, it is still a challenge to understand how the regulation of metabolism and mtDNA expression works at the cellular level and how retrograde signaling from the mitochondria coordinates all those processes.

Organelle Genomes and Transcriptomes of Nymphaea Reveal the Interplay between Intron Splicing and RNA Editing

Posttranscriptional modifications, including intron splicing and RNA editing, are common processes during regulation of gene expression in plant organelle genomes. However, the intermediate products of intron-splicing, and the interplay between intron-splicing and RNA-editing were not well studied. Most organelle transcriptome analyses were based on the Illumina short reads which were unable to capture the full spectrum of transcript intermediates within an organelle. To fully investigate the intermediates during intron splicing and the underlying relationships with RNA editing, we used PacBio DNA-seq and Iso-seq, together with Illumina short reads genome and transcriptome sequencing data to assemble the chloroplast and mitochondrial genomes of Nymphaea 'Joey Tomocik' and analyze their posttranscriptional features. With the direct evidence from Iso-seq, multiple intermediates partially or fully intron-spliced were observed, and we also found that both cis- and trans-splicing introns were spliced randomly. Moreover, by using rRNA-depleted and non-Oligo(dT)-enrichment strand-specific RNA-seq data and combining direct SNP-calling and transcript-mapping methods, we identified 98 and 865 RNA-editing sites in the plastome and mitogenome of N. 'Joey Tomocik', respectively. The target codon preference, the tendency of increasing protein hydrophobicity, and the bias distribution of editing sites are similar in both organelles, suggesting their common evolutionary origin and shared editing machinery. The distribution of RNA editing sites also implies that the RNA editing sites in the intron and exon regions may splice synchronously, except those exonic sites adjacent to intron which could only be edited after being intron-spliced. Our study provides solid evidence for the multiple intermediates co-existing during intron-splicing and their interplay with RNA editing in organelle genomes of a basal angiosperm.

Rapid shifts in mitochondrial tRNA import in a plant lineage with extensive mitochondrial tRNA gene loss

In most eukaryotes, transfer RNAs (tRNAs) are one of the very few classes of genes remaining in the mitochondrial genome, but some mitochondria have lost these vestiges of their prokaryotic ancestry. Sequencing of mitogenomes from the flowering plant genus Silene previously revealed a large range in tRNA gene content, suggesting rapid and ongoing gene loss/replacement. Here, we use this system to test longstanding hypotheses about how mitochondrial tRNA genes are replaced by importing nuclear-encoded tRNAs. We traced the evolutionary history of these gene loss events by sequencing mitochondrial genomes from key outgroups (Agrostemma githago and Silene [=Lychnis] chalcedonica). We then performed the first global sequencing of purified plant mitochondrial tRNA populations to characterize the expression of mitochondrial-encoded tRNAs and the identity of imported nuclear-encoded tRNAs. We also confirmed the utility of high-throughput sequencing methods for the detection of tRNA import by sequencing mitochondrial tRNA populations in a species (Solanum tuberosum) with known tRNA trafficking patterns. Mitochondrial tRNA sequencing in Silene revealed substantial shifts in the abundance of some nuclear-encoded tRNAs in conjunction with their recent history of mt-tRNA gene loss and surprising cases where tRNAs with anticodons still encoded in the mitochondrial genome also appeared to be imported. These data suggest that nuclear-encoded counterparts are likely replacing mitochondrial tRNAs even in systems with recent mitochondrial tRNA gene loss, and the redundant import of a nuclear-encoded tRNA may provide a mechanism for functional replacement between translation systems separated by billions of years of evolutionary divergence.

Molecular and functional diversity of organelle RNA editing mediated by RNA recognition motif-containing protein ORRM4 in tomato

Plant organellar RNA editing is a distinct type of post-transcriptional RNA modification that is critical for plant development. We showed previously that the RNA editing factor SlORRM4 is required for mitochondrial function and fruit ripening in tomato (Solanum lycopersicum). However, a comprehensive atlas of the RNA editing mediated by SlORRM4 is lacking. We observed that SlORRM4 is targeted to both chloroplasts and mitochondria, and its knockout results in pale-green leaves and delayed fruit ripening. Using high-throughput sequencing, we identified 12 chloroplast editing sites and 336 mitochondrial editing sites controlled by SlORRM4, accounting for 23% of chloroplast sites in leaves and 61% of mitochondrial sites in fruits, respectively. Analysis of native RNA immunoprecipitation sequencing revealed that SlORRM4 binds to 31 RNA targets; 19 of these targets contain SlORRM4-dependent editing sites. Large-scale analysis of putative SlORRM4-interacting proteins identified SlRIP1b, a RIP/MORF protein. Moreover, functional characterization demonstrated that SlRIP1b is involved in tomato fruit ripening. Our results indicate that SlORRM4 binds to RNA targets and interacts with SlRIP1b to broadly affect RNA editing in tomato organelles. These results provide insights into the molecular and functional diversity of RNA editing factors in higher plants.

Plant organellar RNA editing: what 30 years of research has revealed

The central dogma in biology defines the flow of genetic information from DNA to RNA to protein. Accordingly, RNA molecules generally accurately follow the sequences of the genes from which they are transcribed. This rule is transgressed by RNA editing, which creates RNA products that differ from their DNA templates. Analyses of the RNA landscapes of terrestrial plants have indicated that RNA editing (in the form of C-U base transitions) is highly prevalent within organelles (that is, mitochondria and chloroplasts). Numerous C→U conversions (and in some plants also U→C) alter the coding sequences of many of the organellar transcripts and can also produce translatable mRNAs by creating AUG start sites or eliminating premature stop codons, or affect the RNA structure, influence splicing and alter the stability of RNAs. RNA-binding proteins are at the heart of post-transcriptional RNA expression. The C-to-U RNA editing process in plant mitochondria involves numerous nuclear-encoded factors, many of which have been identified as pentatricopeptide repeat (PPR) proteins that target editing sites in a sequence-specific manner. In this review we report on major discoveries on RNA editing in plant organelles, since it was first documented 30 years ago.

The Whole Genome Assembly and Comparative Genomic Research of Thellungiella parvula (Extremophile Crucifer) Mitochondrion

The complete nucleotide sequences of the mitochondrial (mt) genome of an extremophile species Thellungiella parvula (T. parvula) have been determined with the lengths of 255,773 bp. T. parvula mt genome is a circular sequence and contains 32 protein-coding genes, 19 tRNA genes, and three ribosomal RNA genes with a 11.5% coding sequence. The base composition of 27.5% A, 27.5% T, 22.7% C, and 22.3% G in descending order shows a slight bias of 55% AT. Fifty-three repeats were identified in the mitochondrial genome of T. parvula, including 24 direct repeats, 28 tandem repeats (TRs), and one palindromic repeat. Furthermore, a total of 199 perfect microsatellites have been mined with a high A/T content (83.1%) through simple sequence repeat (SSR) analysis and they were distributed unevenly within this mitochondrial genome. We also analyzed other plant mitochondrial genomes' evolution in general, providing clues for the understanding of the evolution of organelles genomes in plants. Comparing with other Brassicaceae species, T. parvula is related to Arabidopsis thaliana whose characters of low temperature resistance have been well documented. This study will provide important genetic tools for other Brassicaceae species research and improve yields of economically important plants.

RNA editing in plant mitochondria—connecting RNA target sequences and acting proteins

RNA editing changes several hundred cytidines to uridines in the mRNAs of mitochondria in flowering plants. The target cytidines are identified by a subtype of PPR proteins characterized by tandem modules which each binds with a specific upstream nucleotide. Recent progress in correlating repeat structures with nucleotide identities allows to predict and identify target sites in mitochondrial RNAs. Additional proteins have been found to play a role in RNA editing; their precise function still needs to be elucidated. The enzymatic activity performing the C to U reaction may reside in the C-terminal DYW extensions of the PPR proteins; however, this still needs to be proven. Here we update recent progress in understanding RNA editing in flowering plant mitochondria.

The life of plant mitochondrial complex I

The mitochondrial NADH dehydrogenase complex (complex I) of the respiratory chain has several remarkable features in plants: (i) particularly many of its subunits are encoded by the mitochondrial genome, (ii) its mitochondrial transcripts undergo extensive maturation processes (e.g. RNA editing, trans-splicing), (iii) its assembly follows unique routes, (iv) it includes an additional functional domain which contains carbonic anhydrases and (v) it is, indirectly, involved in photosynthesis. Comprising about 50 distinct protein subunits, complex I of plants is very large. However, an even larger number of proteins are required to synthesize these subunits and assemble the enzyme complex. This review aims to follow the complete "life cycle" of plant complex I from various molecular perspectives. We provide arguments that complex I represents an ideal model system for studying the interplay of respiration and photosynthesis, the cooperation of mitochondria and the nucleus during organelle biogenesis and the evolution of the mitochondrial oxidative phosphorylation system.

RNA editing in mitochondrial trans-introns is required for splicing

In plant mitochondria, gene expression of translatable mRNAs is a complex process with two critical steps, RNA editing and splicing. We studied the role of RNA editing on non-coding regions of the mat-r-nad1e-nad5c transcript from wheat mitochondria. This RNA contains two trans-introns, 3'-nad1-I4 and 3'-nad5-I2, involved in different trans-splicing events, ensuring the association of nad1d-nad1e and nad5b-nad5c exons from nad1 and nad5 mRNAs respectively. The C-to-U editing changes studied here affect homologous positions on 3'-nad1-I4 and 3'-nad5-I2. It is proposed that these base changes are necessary to place an Adenosine residue in a bulging conformation characteristic of domain VI (D6) from group II introns. In this work, we investigated the role of RNA editing events on 3'-nad1-I4 and 3'-nad5-I2 in the trans-splicing process using in vivo and in organello approaches. When the branched intermediates formed during the splicing process were analyzed, the C residues from D6 intron domains from 3'-nad1-I4 and 3'-nad5-I2 were found changed to U, suggesting that RNA editing of these residues could be mandatory for splicing. This assumption was tested by expressing recombinant mat-r-nad1e transgenes introduced into mitochondria by electroporation. Mutation of the editing target residue dramatically affected trans-splicing. Interestingly, the exon joining efficiency was not recovered by compensatory mutations, suggesting that the role of RNA editing is not confined to the restoration of the secondary structure of domain D6 of the intron. Our results strongly support the hypothesis that RNA editing in trans-introns precedes maturation, and is required for the splicing reaction. In addition, this is the first report using an in organello approach to study the trans-splicing process, opening the way to future studies of this peculiar mechanism.

Intron RNA editing is essential for splicing in plant mitochondria

Most plant mitochondria messenger RNAs (mRNAs) undergo editing through C-to-U conversions located mainly in exon sequences. However, some RNA editing events are found in non-coding regions at critical positions in the predicted secondary and tertiary structures of introns, suggesting that RNA editing could be important for splicing. Here, we studied the relationships between editing and splicing of the mRNA encoding the ribosomal protein S10 (rps10), which has a group II intron and five editing sites. Two of them, C2 and C3, predicted to stabilize the folded structure of the intron necessary for splicing, were studied by using rps10 mutants introduced into isolated potato mitochondria by electroporation. While mutations of C2 involved in EBS2/IBS2 interactions did not affect splicing, probably by the presence of an alternative EBS2′ region in domain I of the intron, the edition of site C3 turned out to be critical for rps10 mRNA splicing; only the edited (U) form of the transcript was processed. Interestingly, RNA editing was strongly reduced in transcripts from two different intronless genes, rps10 from potato and cox2 from wheat, suggesting that efficient RNA processing may require a close interaction of factors engaged in different maturation processes. This is the first report linking editing and splicing in conditions close to the in vivo situation.

The long and winding road to RNA editing in plant mitochondria: the Tübingen-Berlin chapter

It took several independent observations of C-to-T differences between genomic mtDNA sequences and corresponding complementary DNA (cDNA) sequences before RNA editing in plant mitochondria was accepted as a fact by the group at Tübingen and later Berlin (Hiesel et al., Science246 (1989) 1632-1634). The first such deviating sequence runs were critically viewed in the lab as being errors of some kind, most likely cloning artifacts, which occur only too frequently. Several such cDNA-mtDNA differences identified in independent cDNA clones in different libraries and finally CGG to TGG codon changes dispelled the skeptical view, and this phenomenon was finally recognized as plant mitochondrial RNA editing of a type similar to the apolipoprotein B RNA editing in mammals.

Trans-splicing of organelle introns--a detour to continuous RNAs

In eukaryotes, RNA trans-splicing is an important RNA-processing form for the end-to-end ligation of primary transcripts that are derived from separately transcribed exons. So far, three different categories of RNA trans-splicing have been found in organisms as diverse as algae to man. Here, we review one of these categories: the trans-splicing of discontinuous group II introns, which occurs in chloroplasts and mitochondria of lower eukaryotes and plants. Trans-spliced exons can be predicted from DNA sequences derived from a large number of sequenced organelle genomes. Further molecular genetic analysis of mutants has unravelled proteins, some of which being part of high-molecular-weight complexes that promote the splicing process. Based on data derived from the alga Chlamydomonas reinhardtii, a model is provided which defines the composition of an organelle spliceosome. This will have a general relevance for understanding the function of RNA-processing machineries in eukaryotic organelles.

A trans-splicing group I intron and tRNA-hyperediting in the mitochondrial genome of the lycophyte Isoetes engelmannii

Plant mitochondrial genomes show much more evolutionary plasticity than those of animals. We analysed the first mitochondrial DNA (mtDNA) of a lycophyte, the quillwort Isoetes engelmannii, which is separated from seed plants by more than 350 million years of evolution. The Isoetes mtDNA is particularly rich in recombination events, and chloroplast as well as nuclear DNA inserts document the incorporation of foreign sequences already in this most ancestral vascular plant lineage. On the other hand, particularly small group II introns and short intergenic regions reveal a tendency of evolution towards a compact mitochondrial genome. RNA editing reaches extreme levels exceeding 100 pyrimidine exchanges in individual mRNAs and, hitherto unobserved in such frequency, also in tRNAs with 18 C-to-U conversions in the tRNA for proline. In total, some 1500 sites of RNA editing can be expected for the Isoetes mitochondrial transcriptome. As a unique molecular novelty, the Isoetes cox1 gene requires trans-splicing via a discontinuous group I intron demonstrating disrupted, but functional, RNAs for yet another class of natural ribozymes.

Angiosperm mitochondrial genomes and mutations

Flowering plants harbor the largest mitochondrial genomes reported so far. At present, the nucleotide sequences of 15 mitochondrial genomes from seven angiosperm species are available, making detailed comparative analysis feasible. The gene content is variable among the species, but the most striking feature is the fluidity of intergenic regions, where species-specific sequences predominate. Additionally, angiosperm mitochondrial genomes, even within a species, show a remarkable amount of rearrangement. We also review mitochondrial mutants in angiosperms from a genomic viewpoint, and discuss how they have arisen. The involvement of nuclear genes in mitochondrial genome stability and organization is currently being revealed through the analysis of mutants.

Cis- and trans-splicing of group II introns in plant mitochondria

Group II-type introns in the mitochondrial genes of flowering plants belong to the ribozymic, mobile retroelement family, but not all exhibit conventional structural features and some follow unusual splicing pathways. Moreover, several introns have been disrupted by DNA rearrangements, so that separately-transcribed precursors undergo splicing in trans. RNA processing in plant mitochondria has the added complexity of C-to-U RNA editing which also sometimes occurs within core intron structures or at exon sites very close to introns. It appears that mitochondrial introns in flowering plants have followed quite different evolutionary pathways than other group II introns.

Evolutionary origin of a plant mitochondrial group II intron from a reverse transcriptase/maturase-encoding ancestor

Group II introns are widespread in plant cell organelles. In vivo, most if not all plant group II introns do not self-splice but require the assistance of proteinaceous splicing factors. In some cases, a splicing factor (also referred to as maturase) is encoded within the intronic sequence and produced by translation of the (excised) intron RNA. However, most present-day group II introns in plant organellar genomes do not contain open reading frames (ORFs) for splicing factors, and their excision may depend on proteins encoded by other organellar introns or splicing factors encoded in the nuclear genome. Whether or not the ancestors of all of these noncoding organellar introns originally contained ORFs for maturases is currently unknown. Here we show that a noncoding intron in the mitochondrial cox2 gene of seed plants is likely to be derived from an ancestral reverse transcriptase/maturase-encoding form. We detected remnants of maturase and reverse transcriptase sequences in the 2.7 kb cox2 intron of Ginkgo biloba, the only living species of an ancient gymnosperm lineage, suggesting that the intron originally harbored a splicing factor. This finding supports the earlier proposed hypothesis that the ancient group II introns that invaded organellar genomes were autonomous genetic entities in that they encoded the factor(s) required for their own excision and mobility.

Multiple physical forms of excised group II intron RNAs in wheat mitochondria

Plant mitochondrial group II introns do not all possess hallmark ribozymic features such as the bulged adenosine involved in lariat formation. To gain insight into their splicing pathways, we have examined the physical form of excised introns in germinating wheat embryos. Using RT–PCR and cRT–PCR, we observed conventional lariats consistent with a two-step transesterification pathway for introns such as nad2 intron 4, but this was not the case for the cox2 intron or nad1 intron 2. For cox2, we detected full-length linear introns, which possess non-encoded 3′terminaladenosines, as well as heterogeneous circular introns, which lack 3′ nucleotide stretches. These observations are consistent with hydrolytic splicing followed by polyadenylation as well as an in vivo circularization pathway, respectively. The presence of both linear and circular species in vivo is supported by RNase H analysis. Furthermore, the nad1 intron 2, which lacks a bulged nucleotide at the branchpoint position, comprised a mixed population of precisely full-length molecules and circular ones which also include a short, discrete block of non-encoded nucleotides. The presence of these various linear and circular forms of excised intron molecules in plant mitochondria points to multiple novel group II splicing mechanisms in vivo.

Gene expression studies in isolated mitochondria: Solanum tuberosum rps10 is recognized by cognate potato but not by the transcription, splicing and editing machinery of wheat mitochondria

The complex gene expression mechanisms that occur in plant mitochondria, such as RNA editing and splicing, are not yet well understood. RNA editing in higher plant mitochondria is a highly specific process which modifies mRNA sequences by C-to-U conversions. It has been suggested that in some cases this process is required for splicing. Here, we use an experimental model based on the introduction of DNA into isolated mitochondria by electroporation to study organellar gene expression events. Our aim was to compare processing and editing of potato small ribosomal protein 10 gene (rps10) transcripts in heterologous (wheat mitochondria) and homologous (potato mitochondria) contexts. rps10 is a suitable model because it contains a group II intron, is absent in wheat mitochondria but is actively expressed in potato mitochondria, where transcripts are spliced and undergo five C-to-U editing events. For this purpose, conditions for electroporating isolated potato mitochondria were established. rps10 was placed under the control of either potato or wheat cox2 promoters. We found that rps10 was only transcribed under the control of a cognate promoter. In wheat mitochondria, rps10 transcripts were neither spliced nor edited while they are correctly processed in potato mitochondria. Interestingly, a wheat editing site grafted into rps10 was not recognized by wheat mitochondria but was correctly edited in potato mitochondria. Taken together, these results suggest that editing might occur only when the transcripts are engaged in processing and that they would not be available to editing factors outside of a putative RNA maturation machinery complex.

Group II intron retroelements: function and diversity

Group II introns are a class of retroelements capable of carrying out both self-splicing and retromobility reactions. In recent years, the number of known group II introns has increased dramatically, particularly in bacteria, and the new information is altering our understanding of these intriguing elements. Here we review the basic properties of group II introns, and summarize the differences between the organellar and bacterial introns with regard to structures, insertion patterns and inferred behaviors. We also discuss the evolution of group II introns, as they are the putative ancestors of spliceosomal introns and possibly non-LTR retroelements, and may have played an important role in the development of eukaryote genomes.

Regulation of TRPV1 by a novel renally expressed rat TRPV1 splice variant

The capsaicin receptor and transient receptor potential channel TRPV1 senses heat, protons, and vanilloid agonists in peripheral sensory ganglia. Abundant data have suggested the presence of potentially novel splice variants in the kidney. We report a novel rat TRPV1 splice variant, TRPV1(VAR), cloned from kidney papilla. TRPV1(VAR) cDNA was identified in multiple kidney tissues. Its sequence was fully compatible with potential splice donor and acceptor sites in the rat TRPV1 gene. TRPV1(VAR) is predicted to encode a truncated form of TRPV1 consisting of the NH2-terminal 248 residues of TRPV1 (all within the NH2-terminal intracellular domain) followed by five nonconsensus amino acids (Arg-Glu-Ala-Met-Trp) and a stop codon. The variant utilizes the same consensus Kozak sequence as canonical TRPV1. A band of the appropriate molecular mass was identified in rat kidney papillary (but not medullary) lysates immunoblotted with an antibody directed against the NH2 terminus of TRPV1, whereas an antibody recognizing the TRPV1 COOH terminus failed to detect it. Upon heterologous expression in HEK 293 cells, TRPV1(VAR) potentiated the ability of cotransfected TRPV1 to confer calcium influx in response to resiniferatoxin. TRPV1(VAR) did not influence expression or cell surface localization of cotransfected TRPV1. TRPV1(VAR) protein product associated with the NH2 terminus of canonical TRPV1. Interestingly, when expressed in the COS-7 epithelial cell line, TRPV1(VAR) functioned in a dominant-negative acting capacity, partially blocking TRPV1-dependent resiniferatoxin responsiveness. We conclude that TRPV1(VAR) is one of perhaps several TRPV1 splice variants expressed in rat kidney and that it may serve to modulate TRPV1 responsiveness in some tissues.

Many independent origins of trans splicing of a plant mitochondrial group II intron

We examined the cis- vs. trans-splicing status of the mitochondrial group II intron nad1i728 in 439 species (427 genera) of land plants, using both Southern hybridization results (for 416 species) and intron sequence data from the literature. A total of 164 species (157 genera), all angiosperms, was found to have a trans-spliced form of the intron. Using a multigene land plant phylogeny, we infer that the intron underwent a transition from cis to trans splicing 15 times among the sampled angiosperms. In 10 cases, the intron was fractured between its 5' end and the intron-encoded matR gene, while in the other 5 cases the fracture occurred between matR and the 3' end of the intron. The 15 intron fractures took place at different time depths during the evolution of angiosperms, with those in Nymphaeales, Austrobaileyales, Chloranthaceae, and eumonocots occurring early in angiosperm evolution and those in Syringodium filiforme, Hydrocharis morsus- ranae, Najas, and Erodium relatively recently. The trans-splicing events uncovered in Austrobaileyales, eumonocots, Polygonales, Caryophyllales, Sapindales, and core Rosales reinforce the naturalness of these major clades of angiosperms, some of which have been identified solely on the basis of recent DNA sequence analyses.

Nanoarchaeum equitans creates functional tRNAs from separate genes for their 5'- and 3'-halves

Analysis of the genome sequence of the small hyperthermophilic archaeal parasite Nanoarchaeum equitans has not revealed genes encoding the glutamate, histidine, tryptophan and initiator methionine transfer RNA species. Here we develop a computational approach to genome analysis that searches for widely separated genes encoding tRNA halves that, on the basis of structural prediction, could form intact tRNA molecules. A search of the N. equitans genome reveals nine genes that encode tRNA halves; together they account for the missing tRNA genes. The tRNA sequences are split after the anticodon-adjacent position 37, the normal location of tRNA introns. The terminal sequences can be accommodated in an intervening sequence that includes a 12-14-nucleotide GC-rich RNA duplex between the end of the 5' tRNA half and the beginning of the 3' tRNA half. Reverse transcriptase polymerase chain reaction and aminoacylation experiments of N. equitans tRNA demonstrated maturation to full-size tRNA and acceptor activity of the tRNA(His) and tRNA(Glu) species predicted in silico. As the joining mechanism possibly involves tRNA trans-splicing, the presence of an intron might have been required for early tRNA synthesis.

Distribution of introns in the mitochondrial gene nad1 in land plants: phylogenetic and molecular evolutionary implications

Forty-six species of diverse land plants were investigated by sequencing for their intron content in the mitochondrial gene nad1. A total of seven introns, all belonging to group II, were found, and two were newly discovered in this study. All 13 liverworts examined contain no intron, the same condition as in green algae. Mosses and hornworts, however, share one intron by themselves and another one with vascular plants. These intron distribution patterns are consistent with the hypothesis that liverworts represent the basal-most land plants and that the two introns were gained in the common ancestor of mosses-hornworts-vascular plants after liverworts had diverged. Hornworts also possess a unique intron of their own. A fourth intron was found only in Equisetum L., Marattiaceae, Ophioglossum L., Osmunda L., Asplenium L., and Adiantum L., and was likely acquired in their common ancestor, which supports the monophyly of moniliformopses. Three introns that were previously characterized in angiosperms and a few pteridophytes are now all extended to lycopods, and were likely gained in the common ancestor of vascular plants. Phylogenetic analyses of the intron sequences recovered topologies mirroring those of the plants, suggesting that the introns have all been vertically inherited. All seven nad1 group II introns show broad phylogenetic distribution patterns, with the narrowest being in moniliformopses and hornworts, lineages that date back to at least the Devonian (345 million years ago) and Silurian (435 million years ago), respectively. Hence, these introns must have invaded the genes via ancient transpositional events during the early stage of land plant evolution. Potentially heavy RNA editing was observed in nad1 of Haplomitrium Dedecek, Takakia Hatt. & Inoue, hornworts, Isoetes L., Ophioglossum, and Asplenium. A new nomenclature is proposed for group II introns.

The mitochondrial DNA of land plants: peculiarities in phylogenetic perspective

Land plants exhibit a significant evolutionary plasticity in their mitochondrial DNA (mtDNA), which contrasts with the more conservative evolution of their chloroplast genomes. Frequent genomic rearrangements, the incorporation of foreign DNA from the nuclear and chloroplast genomes, an ongoing transfer of genes to the nucleus in recent evolutionary times and the disruption of gene continuity in introns or exons are the hallmarks of plant mtDNA, at least in flowering plants. Peculiarities of gene expression, most notably RNA editing and trans-splicing, are significantly more pronounced in land plant mitochondria than in chloroplasts. At the same time, mtDNA is generally the most slowly evolving of the three plant cell genomes on the sequence level, with unique exceptions in only some plant lineages. The slow sequence evolution and a variable occurrence of introns in plant mtDNA provide an attractive reservoir of phylogenetic information to trace the phylogeny of older land plant clades, which is as yet not fully resolved. This review attempts to summarize the unique aspects of land plant mitochondrial evolution from a phylogenetic perspective.

Horizontal gene transfer from flowering plants to Gnetum

Although horizontal gene transfer is well documented in microbial genomes, no case has been reported in higher plants. We discovered horizontal transfer of the mitochondrial nad1 intron 2 and adjacent exons b and c from an asterid to Gnetum (Gnetales, gymnosperms). Gnetum has two copies of intron 2, a group II intron, that differ in their exons, nucleotide composition, domain lengths, and structural characteristics. One of the copies, limited to an Asian clade of Gnetum, is almost identical to the homologous locus in angiosperms, and partial sequences of its exons b and c show characteristic substitutions unique to angiosperms. Analyses of 70 seed plant nad1 exons b and c and intron 2 sequences, including representatives of all angiosperm clades, support that this copy originated from a euasterid and was horizontally transferred to Gnetum. Molecular clock dating, using calibrations provided by gnetalean macrofossils, suggests an age of 5 to 2 million years for the Asian clade that received the horizontal transfer.

Inferences from mitochondrial DNA patterns on the domestication history of alfalfa (Medicago sativa)

The evolutionary history of a wild-cultivated complex is the outcome of different factors, among which are the domestication of the cultivated form, its geographical and demographic expansion, as well as gene flow between natural and cultivated populations. To clarify this history for alfalfa (Medicago sativa), we adopted a phylogeographical approach based on a sample of natural and cultivated populations of the M. sativa species complex. This sample was characterized for mtDNA variation through restriction fragment length polymorphisms. Twenty-two mitotypes were identified in the whole data set (155 individuals). In the wild pool from the presumed area of origin (Near East to Central Asia), the diversity was high but the absence of geographical differentiation hinders a more precise location of the centre of domestication. Within the cultivated alfalfa, the geographical structure suggests strongly the existence of at least two independent routes of dissemination of alfalfa from its centre of origin. In particular, original mitotypes detected in some regions suggest that wild populations not sampled in this study, and located outside the presumed area of domestication, may have contributed to the diversity of the cultivated pool. Finally, analysis of the mtDNA variation shows that natural populations collected in Spain constitute an endemic wild pool now introgressed partially by cultivated alfalfa.

Lariat formation and a hydrolytic pathway in plant chloroplast group II intron splicing

Lariat formation has been studied intensively only with a few self-splicing group II introns, and little is known about how the numerous diverse introns in plant organelles are excised. Several of these introns have branch-points that are not a single bulge but are adjoined by A:A, A:C, A:G and G:G pairs. Using a highly sensitive in vivo approach, we demonstrate that all but one of the barley chloroplast introns splice via the common pathway that produces a branched product. RNA editing does not improve domain 5 and 6 structures of these introns. The conserved branch-point in tobacco rpl16 is chosen even if an adjacent unpaired adenosine is available, suggesting that spatial arrangements in domain 6 determine correct branch-point selection. Lariats were not detected for the chloroplast trnV intron, which lacks an unpaired adenosine in domain 6. Instead, this intron is released as linear molecules that undergo further polyadenylation. trnV, which is conserved throughout plant evolution, constitutes the first example of naturally occurring hydrolytic group II intron splicing in vivo.

The RNA world of plant mitochondria

Mitochondria are well known as the cellular power factory. Much less is known about these organelles as a genetic system. This is particularly true for mitochondria of plants, which subsist with respect to attention by the scientific community in the shadow of the chloroplasts. Nevertheless the mitochondrial genetic system is essential for the function of mitochondria and thus for the survival of the plant. In plant mitochondria the pathway from the genetic information encoded in the DNA to the functional protein leads through a very diverse RNA world. How the RNA is generated and what kinds of regulation and control mechanisms are operative in transcription are current topics in research. Furthermore, the modes of posttranscriptional alterations and their consequences for RNA stability and thus for gene expression in plant mitochondria are currently objects of intensive investigations. In this article current results obtained in the examination of plant mitochondrial transcription, RNA processing, and RNA stability are illustrated. Recent developments in the characterization of promoter structure and the respective transcription apparatus as well as new aspects of RNA processing steps including mRNA 3' processing and stability, mRNA polyadenylation, RNA editing, and tRNA maturation are presented. We also consider new suggestions concerning the endosymbiont hypothesis and evolution of mitochondria. These novel considerations may yield important clues for the further analysis of the plant mitochondrial genetic system. Conversely, an increasing knowledge about the mechanisms and components of the organellar genetic system might reveal new aspects of the evolutionary history of mitochondria.

RNA editing of the ribosomal protein S13 transcripts in magnolia and sunflower mitochondria

To obtain information on ribosomal protein gene content and RNA editing in plant mitochondria, we investigated the conservation of the gene for ribosomal protein S13 (rps13) in the mitochondrial genomes of Magnolia spp. and Helianthus annuus. A complete rps13 open reading frame was identified in the mitochondrial genome of Magnolia and Helianthus. Comparison of genomic and cDNA sequences revealed that transcripts from the mitochondrial rps13 locus in both Magnolia and Helianthus undergo RNA editing suggesting the functionality of the related gene. Six and four editing events in the rps13 transcript population of Magnolia and Helianthus mitochondria, respectively, cause critical changes of codons and, consequently, induce 5 and 3.4% amino acid modifications in the respective genomically-encoded S13 polypeptides. Both editing pattern and efficiency of RNA editing differ in the reverse transcription (RT)-PCR-derived cDNA populations from the rps13 locus from the analyzed plant species with only the Magnolia rps13 being fully edited.

Functions and mechanisms of RNA editing

RNA editing can be broadly defined as any site-specific alteration in an RNA sequence that could have been copied from the template, excluding changes due to processes such as RNA splicing and polyadenylation. Changes in gene expression attributed to editing have been described in organisms from unicellular protozoa to man, and can affect the mRNAs, tRNAs, and rRNAs present in all cellular compartments. These sequence revisions, which include both the insertion and deletion of nucleotides, and the conversion of one base to another, involve a wide range of largely unrelated mechanisms. Recent advances in the development of in vitro editing and transgenic systems for these varied modifications have provided a better understanding of similarities and differences between the biochemical strategies, regulatory sequences, and cellular factors responsible for such RNA processing events.

From gene to protein in higher plant mitochondria

Higher plant mitochondria contain a genetic system with a genome, transcription and translation processes, which have to be logistically integrated with the two other genomes in the nucleus and the plastid. In plant mitochondria, after transcripts have been synthesised, at least in some cases by a phage-type RNA polymerase, they have to go through a complex processing apparatus, which depends on protein factors imported from the cytosol. Processing involves cis- and trans-splicing, internal RNA editing and maturation at the transcript termini, these steps often occurring in parallel. Transcript life is terminated by RNA degradation mechanisms, one of which involves polyadenylation. RNA metabolism seems to be a key element of the regulation of gene expression in higher plant mitochondria.

Heterogeneous Sp1 mRNAs in human HepG2 cells include a product of homotypic trans-splicing

Sp1 is one of the well documented transcription factors, but the whole structure of human Sp1 has not been determined yet. In the present study, we isolated several cDNAs representing two forms of human Sp1 mRNA with different 5'-terminal structures in HepG2 cells. Isolation of a genomic clone established that one of the cDNAs represents the mRNA having consecutive alignment of exons, which allowed deducing the complete amino acid sequence for human Sp1. Another cDNA clone had a surprising structure that possessed an alignment of exons 3-2-3. Both reverse transcriptase-polymerase chain reaction and RNase protection assays confirmed accumulation of the two forms of Sp1 mRNA in HepG2 cells. Because Southern blot analysis suggested that exon 3 is of a single copy in the genome, the cDNA clone having the duplicated sequences for exon 3 appeared to reflect the trans-splicing between pre-mRNAs of human Sp1.

Expression of mitochondrial protein-coding genes in Tetrahymena pyriformis

In the ciliate protozoon, Tetrahymena pyriformis, mitochondrial protein-coding genes are highly divergent in sequence, and in a number of cases they lack AUG initiation codons. We asked whether RNA editing might be acting to generate protein sequences that are more conventional than those inferred from the corresponding gene sequences, and/or to create standard AUG initiation codons where these are absent. However, comparison of genomic and cDNA sequences (the latter generated by reverse transcriptase sequencing of T. pyriformis mitochondrial mRNAs) yielded no evidence of mitochondrial RNA editing in this organism. To delineate the 5' ends of mitochondrial protein-coding transcripts, primer extension experiments were conducted. In all cases, 5' termini were found to map within a few nucleotides of potential initiation codons, indicating that T. pyriformis mitochondrial mRNAs have little or no 5' untranslated leader sequence. The pattern of strong primer extension stops suggested that both standard (AUG) and non-standard (AUU, AUA, GUG, UUG) initiation codons are utilized by the Tetrahymena mitochondrial translation system. We also investigated expression of the nad1 gene, which in both T. pyriformis and Paramecium aurelia is split into two portions that are encoded by and transcribed from different DNA strands. Northern hybridization analysis showed that the corresponding transcripts are not trans-spliced, implying that separate N-terminal and C-terminal portions of Nad1 are made in this system. Finally, in a search for primary transcripts, we isolated from a T. pyriformis mitochondrial fraction several small RNAs that were reproducibly labeled by incubation in the presence of [alpha-(32)P]GTP and guanylyltransferase. Partial sequence information revealed that none of these cappable RNAs is encoded in the T. pyriformis mitochondrial genome.

RNA maturation of the rice SPK gene may involve trans-splicing

A gene encoding a calcium-dependent seed-specific protein kinase (SPK) is abundantly expressed in developing rice seeds (Kawasaki, T et al. Gene (1993) 129, 183-189). Rice genomic clones encoding SPK were isolated using the entire cDNA fragment as a probe. Physical mapping of these genomic clones indicated that the genomic region corresponding to the entire cDNA was divided into two different regions, SPK-A and SPK-B, located on different rice chromosomes. The results of RACE-PCR analyses showed that the respective transcripts from SPK-A and SPK-B contained additional sequences which were not found in the SPK cDNA, and that these sequences were removed like introns during maturation of the SPK mRNA. These results suggest that two different RNAs were independently transcribed from SPK-A and SPK-B and joined, possibly by trans-splicing.

In the Nicotiana sylvestris CMSII mutant, a recombination-mediated change 5' to the first exon of the mitochondrial nad1 gene is associated with lack of the NADH:ubiquinone oxidoreductase (complex I) NAD1 subunit

We previously reported that the Nicotiana sylvestris CMSII mutant mitochondrial DNA carried a large deletion. Several expressed sequences, most of which are duplicated, and the unique copy of the nad7 gene encoding the NAD7 subunit of the NADH:ubiquinone oxidoreductase complex (complex I) are found in the deletion. Here, we show that the orf87-nad3-nad1/A cotranscription unit transcribed from a unique promoter element in the wild-type, is disrupted in CMSII. Nad3, orf87 and the promoter element are part of the deleted sequence, whilst the nad1/A sequence is present and transcribed from a new promoter brought by the recombination event, as indicated by Northern and primer extension experiments. However, Western analyses of mitochondrial protein fractions and of complex I purified using anti-NAD9 affinity columns, revealed that NAD1 is lacking in CMSII mitochondria. Our results suggest that translation of nad1 transcripts rather than transcription itself could be altered in the mutant. Consequences of lack of this submit belonging the membrane arm of complex I and thought to contain the ubiquinone-binding site, are discussed.

Physiological, biochemical and molecular aspects of mitochondrial complex I in plants

Respiratory complex I of plant mitochondria has to date been investigated with respect to physiological function, biochemical properties and molecular structure. In the respiratory chain complex I is the major entry gate for low potential electrons from matrix NADH, reducing ubiquinone and utilizing the released energy to pump protons across the inner membrane. Plant complex I is active against a background of several other NAD(P)H dehydrogenases, which do not contribute in proton pumping, but permit and establish several different routes of shuttling electrons from NAD(P)H to ubiquinone. Identification of the corresponding molecular structures, that is the proteins and genes of the different NADH dehydrogenases, will allow more detailed studies of this interactive regulatory network in plant mitochondria. Present knowledge of the structure of complex I and the respective mitochondrial and nuclear genes encoding various subunits of this complex in plants is summarized here. Copyright 1998 Elsevier Science B.V.

Assays for investigating RNA editing in plant mitochondria

In plant organelles transcripts are modified posttranscriptionally by RNA editing. This modification process changes almost every protein-coding RNA at specific cytidine and uridine positions. Therefore, mitochondrially encoded protein sequences differ from the genomically fixed information and show, after editing, a higher conservation. To investigate this unusual processing step in plant mitochondria, several assays have been developed. However, compared with the progress made in other RNA editing fields, knowledge about the factors involved in plant mitochondrial editing is limited. One reason for this is the lack of a reliable in vitro system for mitochondria. To reveal the biochemical nature of the RNA editing reaction in plant mitochondria, we developed an in vitro system by which we were able to show that cytidine is specifically modified to uridine by a deamination or transamination process. Here we describe the development of a pea in vitro system and discuss assays to follow the editing process.

Splicing and intron-internal RNA editing of trnK-matK transcripts in barley plastids: support for MatK as an essential splice factor

Group II introns frequently require assistance by specific factors, maturases, for folding and effective splicing in vivo. The only putative maturase of higher plant chloroplasts is encoded by matK, located in the intron of trnK. We show that in barley matK transcripts are modified at a first codon base by C-to-U RNA editing. The resulting H --> Y substitution restores a sequence motif that is present in maturases of yeast and plant mitochondria and of Lactococcus ltrA and that is positioned within the X domain. Processing of trnK-matK transcripts was further investigated in plastids lacking functional ribosomes due to a mutation. Absence of the intron-encoded matK gene product in these plastids is correlated with the accumulation of precursor transcripts for tRNALys(UUU)-matK, processed to different degrees, and by the lack of mature and spliced tRNA molecules. These results suggest an essential role of MatK for splicing of its own transcript in vivo. Processing of the 5' end of trnK exon 1 was found to proceed efficiently also in the mutant plastids although the two tRNA exons were separated by the 2481 nt intron. Consequently, presence of the intron does not interfere with the formation of mature 5' termini.

Evolution of trans-splicing plant mitochondrial introns in pre-Permian times

Trans-splicing in angiosperm plant mitochondria connects exons from independent RNA molecules by means of group II intron fragments. Homologues of trans-splicing introns in the angiosperm mitochondrial nad2 and nad5 genes are now identified as uninterrupted group II introns in the ferns Asplenium nidus and Marsilea drummondii. These fern introns are correctly spliced from the pre-mRNA at the sites predicted from their well-conserved secondary structures. The flanking exon sequences of the nad2 and nad5 genes in the ferns require RNA editing, including the removal of in-frame stop codons by U-to-C changes for correct expression of the genetic information. We conclude that cis-splicing introns like the ones now identified in ferns are the ancestors of trans-splicing introns in angiosperm mitochondria. Intron disruption is apparently due to a size increase of the structurally variable group II intron domain IV followed by DNA recombination in the plant mitochondrial genome.

RNA editing status of nad7 intron domains in wheat mitochondria

The most highly conserved structures of group II introns are the helical domains V and VI near the 3'splice site. Within this region of each of the four introns in the wheat mitochondrial nad7 gene encoding NADH dehydrogenase subunit 7, there are A-C mispairs. To determine whether C-to-U type RNA editing restores conventional A-U pairing, we sequenced RT-PCR products from partially-spliced nad7 template RNA and gel-fractionated, excised intron RNA. We examined transcripts from germinating wheat embryos and seedlings because these two stages of development show pronounced differences in steady state levels of nad7 intronic RNAs. We observed editing at only two of the six predicted sites, and they were located at homologous positions within domain V of the third and fourth introns. A third site was found to be edited within the unmodelled domain VI loop of the fourth intron. Similar patterns of RNA editing were seen in wheat embryos and seedlings. These observations, and the presence of other non-conventional base pairs particularly within domain V of plant mitochondrial introns, indicate weaker helical core structure than in ribozymic group II introns. Moreover, the incompleteness or absence of editing in wheat nad7 excised intron RNA suggests that, although editing may contribute to splicing efficiency, it is not essential for splicing.

The mitochondrial genome of Arabidopsis thaliana contains 57 genes in 366,924 nucleotides

We have determined the complete sequence of the mitochondrial DNA in the model plant species Arabidopsis thaliana, affording access to the first of its three genomes. The 366,924 nucleotides code for 57 identified genes, which cover only 10% of the genome. Introns in these genes add about 8%, open reading frames larger than 100 amino acids represent 10% of the genome, duplications account for 7%, remnants of retrotransposons of nuclear origin contribute 4% and integrated plastid sequences amount to 1%-leaving 60% of the genome unaccounted for. With the significant contribution of duplications, imported foreign DNA and the extensive background of apparently functionless sequences, the mosaic structure of the Arabidopsis thaliana mitochondrial genome features many aspects of size-relaxed nuclear genomes.