Correlation between amino acid residues converted by RNA editing and functional residues in protein three-dimensional structures in plant organelles

Seeing the unseen in characterizing RNA editome during rice endosperm development

Rice (Oryza sativa L.) endosperm is essential to provide nutrients for seed germination and determine grain yield. RNA editing, a post-transcriptional modification essential for plant development, unfortunately, is not fully characterized during rice endosperm development. Here, we perform systematic analyses to characterize RNA editome during rice endosperm development. We find that most editing sites are C-to-U CDS-recoding in mitochondria, leading to increased hydrophobic amino acids and changed structures of mitochondrial proteins. Comparative analysis of RNA editome reveals that CDS-recoding sites present higher editing frequencies with lower variabilities and their resultant recoded amino acids tend to exhibit stronger evolutionary conservation across many land plants. Furthermore, we classify mitochondrial genes into three groups, presenting distinct patterns in terms of CDS-recoding events. Besides, we conduct genome-wide screening to detect pentatricopeptide repeat (PPR) proteins and construct PPR-RNA binding profiles, yielding candidate PPR editing factors related to rice endosperm development. Taken together, our findings provide valuable insights for deciphering fundamental mechanisms of rice endosperm development underlying RNA editing machinery.

Assembly and characterization of the complete mitochondrial genome of Ventilago leiocarpa

KEY MESSAGE

We reported the mitochondrial genome of Ventilago leiocarpa for the first time. Two and one sites lead to the generation of stop and stat codon through editing were verified. Ventilago leiocarpa, a member of the Rhamnaceae family, is frequently utilized in traditional medicine due to the medicinal properties of its roots. In this study, we successfully assembled the mitogenome of V. leiocarpa using both BGI short reads and Nanopore long reads. This mitogenome has a total length of 331,839 bp. The annotated results showed 36 unique protein-coding, 16 tRNA and 3 rRNA genes in this mitogenome. Furthermore, we confirmed the presence of a branched structure through the utilization of long reads mapping, PCR amplification, and Sanger sequencing. Specifically, the ctg1 can form a single circular molecule or combine with ctg4 to form a linear molecule. Likewise, ctg2 can form a single circular molecule or can be connected to ctg4 to form a linear molecule. Subsequently, through a comparative analysis of the mitogenome and cpgenome sequences, we identified ten mitochondrial plastid sequences (MTPTs), including two complete protein-coding genes and five complete tRNA genes. The existence of MTPTs was verified by long reads. Colinear analysis showed that the mitogenomes of Rosales were highly divergent in structure. Finally, we identified 545 RNA editing sites involving 36 protein-coding genes by Deepred-mt. To validate our findings, we conducted PCR amplification and Sanger sequencing, which confirmed the generation of stop codons in atp9-223 and rps10-391, as well as the generation of a start codon in nad4L-2. This project reported the complex structure and RNA editing event of the V. Leiocarpa mitogenome, which will provide valuable information for the study of mitochondrial gene expression.

A Comprehensive Evolutionary Study of Chloroplast RNA Editing in Gymnosperms: A Novel Type of G-to-A RNA Editing Is Common in Gymnosperms

Although more than 9100 plant plastomes have been sequenced, RNA editing sites of the whole plastome have been experimentally verified in only approximately 21 species, which seriously hampers the comprehensive evolutionary study of chloroplast RNA editing. We investigated the evolutionary pattern of chloroplast RNA editing sites in 19 species from all 13 families of gymnosperms based on a combination of genomic and transcriptomic data. We found that the chloroplast C-to-U RNA editing sites of gymnosperms shared many common characteristics with those of other land plants, but also exhibited many unique characteristics. In contrast to that noted in angiosperms, the density of RNA editing sites in ndh genes was not the highest in the sampled gymnosperms, and both loss and gain events at editing sites occurred frequently during the evolution of gymnosperms. In addition, GC content and plastomic size were positively correlated with the number of chloroplast RNA editing sites in gymnosperms, suggesting that the increase in GC content could provide more materials for RNA editing and facilitate the evolution of RNA editing in land plants or vice versa. Interestingly, novel G-to-A RNA editing events were commonly found in all sampled gymnosperm species, and G-to-A RNA editing exhibits many different characteristics from C-to-U RNA editing in gymnosperms. This study revealed a comprehensive evolutionary scenario for chloroplast RNA editing sites in gymnosperms, and reported that a novel type of G-to-A RNA editing is prevalent in gymnosperms.

Comparative chloroplast genome and transcriptome analysis on the ancient genus Isoetes from China

Isoetes is a famous living fossil that plays a significant role in the evolutionary studies of the plant kingdom. To explore the adaptive evolution of the ancient genus Isoetes from China, we focused on Isoetes yunguiensis (Q.F. Wang and W.C. Taylor), I. shangrilaensis (X. Li, Y.Q. Huang, X.K. Dai & X. Liu), I. taiwanensis (DeVol), I. sinensis (T.C. Palmer), I. hypsophila_GHC (Handel-Mazzetti), and I. hypsophila_HZS in this study. We sequenced, assembled, and annotated six individuals' chloroplast genomes and transcriptomes, and performed a series of analyses to investigate their chloroplast genome structures, RNA editing events, and adaptive evolution. The six chloroplast genomes of Isoetes exhibited a typical quadripartite structure with conserved genome sequence and structure. Comparative analyses of Isoetes species demonstrated that the gene organization, genome size, and GC contents of the chloroplast genome are highly conserved across the genus. Besides, our positive selection analyses suggested that one positively selected gene was statistically supported in Isoetes chloroplast genomes using the likelihood ratio test (LRT) based on branch-site models. Moreover, we detected positive selection signals using transcriptome data, suggesting that nuclear-encoded genes involved in the adaption of Isoetes species to the extreme environment of the Qinghai-Tibetan Plateau (QTP). In addition, we identified 291-579 RNA editing sites in the chloroplast genomes of six Isoetes based on transcriptome data, well above the average of angiosperms. RNA editing in protein-coding transcripts results from amino acid changes to increase their hydrophobicity and conservation in Isoetes, which may help proteins form functional three-dimensional structure. Overall, the results of this study provide comprehensive transcriptome and chloroplast genome resources and contribute to a better understanding of adaptive evolutionary and molecular biology in Isoetes.

Homologous Cloning of Potassium Channel Genes From the Superior Apple Rootstock Line 12-2, Which is Tolerant to Apple Replant Disease

Potassium channels are important ion channels that are responsible for the absorption of potassium in the plant nutrient uptake system. In this study, we used homologous molecular cloning to obtain 8 K⁺ channel genes from the superior apple rootstock line 12-2 (self-named): MsAKT1-1, MsKAT3-2, MsKAT1-3, MsK2P3-4, MsK2P3-5, MsK2P5-6, MsK2P3-7, and MsK2P3-8. Their lengths varied from 942 bp (MsK2P5-6) to 2625 bp (MsAKT1-1), and the number of encoded amino acids varied from 314 (MsK2P5-6) to 874 (MsAKT1-1). Subcellular localization predictions showed that MsAKT1-1, MsKAT3-2, and MsKAT1-3 were localized on the plasma membrane, and MsK2P3-4, MsK2P3-5, MsK2P5-6, MsK2P3-7, and MsK2P3-8 were localized on the vacuole and plasma membrane. The 8 K⁺ channel proteins contained α helices, extended strands, β turns, and random coils. MsKAT1-3 had four transmembrane structures, MsKAT3-2 had six, and the other six K⁺ channel genes had five. Protein structure domain analysis showed that MsAKT1-1 contained nine protein domains, followed by MsKAT3-2 with four, MsKAT1-3 with three, and the other five two-pore domain K⁺ channel proteins with two. Semi-quantitative RT-PCR detection of the K⁺ channel genes showed that their expression levels were high in roots. qRT-PCR analysis showed that the relative expression levels of the 8 genes changed after exposure to ARD stress. The above results provide a theoretical basis for further research on the functions of potassium channel genes in 12-2 and a scientific basis for the breeding of ARD-resistant rootstock.

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing

Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016 bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnP^GGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

Identification of Structure-Controlling Rice Biosynthesis Enzymes

The main enzymes controlling the chain-length distributions (CLDs) of starches are starch synthases (SSs), starch branching enzymes (SBEs), and debranching enzymes (DBEs), which have various isoforms, denoted as SSI, SSII-1, etc. Different isozymes dominate the CLD in different ranges of degrees of polymerization (DPs). Models have been developed for the CLDs in terms of the activities of isoforms of these enzymes, in terms of two parameters: β_i, which is the ratio of the activity of SBE to that of SS in set i, and h_i, which is the relative activity of SS in that set. These provide good fits to data but without specifying which isozymes are in set i. Here, CLDs for amylopectin and amylose synthesis in rice endosperm are explored. Molecular weight distributions of the different chains formed in 87 rice varieties were obtained using size-exclusion chromatography following enzymatic debranching (converting a complex branched macromolecule to linear polymers), and fitted by the biosynthesis-based models. The mutants of each isoform among tested rice varieties were identified by amino-acid mutations in coding sequences based on the extraction and analysis of whole gene sequences. The significant differences between mutant groups of different isoforms indicate that SSI, SSII-3, SSIII-1, SSIII-2, and SBEI as well as GBSSI (an isozyme of granule-bound starch synthase) belong to the enzymes sets that control amylose biosynthesis. Further, GBSSI is in the enzyme sets that control amylopectin chains. This enables specification of all isozymes and the DP range, which they dominate, over the entire DP range. As the CLD controls many functional properties of rice, this can help breeders target and develop improved rice species.

Comparative Analysis of Complete Chloroplast Genome Sequences of Wild and Cultivated Bougainvillea (Nyctaginaceae)

Bougainvillea (Nyctaginaceae) is a popular ornamental plant group primarily grown for its striking colorful bracts. However, despite its established horticultural value, limited genomic resources and molecular studies have been reported for this genus. Thus, to address this existing gap, complete chloroplast genomes of four species (Bougainvillea glabra, Bougainvillea peruviana, Bougainvillea pachyphylla, Bougainvillea praecox) and one Bougainvillea cultivar were sequenced and characterized. The Bougainvillea cp genomes range from 153,966 bp to 154,541 bp in length, comprising a large single-copy region (85,159 bp-85,708 bp) and a small single-copy region (18,014 bp-18,078 bp) separated by a pair of inverted repeats (25,377-25,427 bp). All sequenced plastomes have 131 annotated genes, including 86 protein-coding, eight rRNA, and 37 tRNA genes. These five newly sequenced Bougainvillea cp genomes were compared to the Bougainvillea spectabilis cp genome deposited in GeBank. The results showed that all cp genomes have highly similar structures, contents, and organization. They all exhibit quadripartite structures and all have the same numbers of genes and introns. Codon usage, RNA editing sites, and repeat analyses also revealed highly similar results for the six cp genomes. The amino acid leucine has the highest proportion and almost all favored synonymous codons have either an A or U ending. Likewise, out of the 42 predicted RNA sites, most conversions were from serine (S) to leucine (L). The majority of the simple sequence repeats detected were A/T mononucleotides, making the cp genomes A/T-rich. The contractions and expansions of the IR boundaries were very minimal as well, hence contributing very little to the differences in genome size. In addition, sequence variation analyses showed that Bougainvillea cp genomes share nearly identical genomic profiles though several potential barcodes, such as ycf1, ndhF, and rpoA were identified. Higher variation was observed in both B. peruviana and B. pachyphylla cp sequences based on SNPs and indels analysis. Phylogenetic reconstructions further showed that these two species appear to be the basal taxa of Bougainvillea. The rarely cultivated and wild species of Bougainvillea (B. pachyphylla, B. peruviana, B. praecox) diverged earlier than the commonly cultivated species and cultivar (B. spectabilis, B. glabra, B. cv.). Overall, the results of this study provide additional genetic resources that can aid in further phylogenetic and evolutionary studies in Bougainvillea. Moreover, genetic information from this study is potentially useful in identifying Bougainvillea species and cultivars, which is essential for both taxonomic and plant breeding studies.

Detection and Analysis of C-to-U RNA Editing in Rice Mitochondria-Encoded ORFs

Cytidine to uridine (C-to-U) RNA editing is an important type of substitutional RNA modification and is almost omnipresent in plant chloroplasts and mitochondria. In rice mitochondria, 491 C-to-U editing sites have been identified previously, and case studies have elucidated the function of several C-to-U editing sites in rice, but the functional consequence of most C-to-U alterations needs to be investigated further. Here, by means of Sanger sequencing and publicly available RNA-seq data, we identified a total of 569 C-to-U editing sites in rice mitochondria-encoded open reading frames (ORFs), 85.41% of these editing sites were observed on the first or the second base of a codon, resulting in the alteration of encoded amino acid. Moreover, we found some novel editing sites and several inaccurately annotated sites which may be functionally important, based on the highly conserved amino acids encoded by these edited codons. Finally, we annotated all 569 C-to-U RNA editing sites in their biological context. More precise information about C-to-U editing sites in rice mitochondria-encoded ORFs will facilitate our investigation on the function of C-to-U editing events in rice and also provide a valid benchmark from rice for the analysis of mitochondria C-to-U editing in other plant species.

Comparative analysis of chloroplast genomes in Vasconcellea pubescens A.DC. and Carica papaya L

The chloroplast genome is an integral part of plant genomes in a species along with nuclear and mitochondrial genomes, contributing to adaptation, diversification, and evolution of plant lineages. In the family Caricaceae, only the Carica papaya chloroplast genome and its nuclear and mitochondrial genomes were sequenced, and no chloroplast genome-wide comparison across genera was conducted. Here, we sequenced and assembled the chloroplast genome of Vasconcellea pubescens A.DC. using Oxford Nanopore Technology. The size of the genome is 158,712 bp, smaller than 160,100 bp of the C. papaya chloroplast genome. And two structural haplotypes, LSC_IRa_SSCrc_IRb and LSC_IRa_SSC_IRb, were identified in both V. pubescens and C. papaya chloroplast genomes. The insertion-deletion mutations may play an important role in Ycf1 gene evolution in family Caricaceae. Ycf2 is the only one gene positively selected in the V. pubescens chloroplast genome. In the C. papaya chloroplast genome, there are 46 RNA editing loci with an average RNA editing efficiency of 63%. These findings will improve our understanding of the genomes of these two crops in the family Caricaceae and will contribute to crop improvement.

Light/heat effects on RNA editing in chloroplast NADH-plastoquinone oxidoreductase subunit 2 (ndhB) gene of Calotropis (Calotropis procera)

Background

RNA editing is common in terrestrial plants, especially in mitochondria and chloroplast. In the photosynthesis process, NAD dehydrogenase plays a very important role. Subunit 2 of NADH-dehydrogenase is one of the major subunits in NAD dehydrogenase complex. Using desert plant Calotropis (Calotropis procera), this study focuses on the RNA editing activity of ndhB based on light time.

Results

NdhB (NADH-dehydrogenase subunit 2) gene accession no. MK144329 was isolated from Calotropis procera genomic data (PRJNA292713). Additionally, using RNA-seq data, the cDNA of the ndhB gene of C. procera was isolated at three daylight periods, i.e., dawn (accession no. MK165161), at midday (accession no. MK165160), and pre-dusk (accession no. MK165159). Seven RNA editing sites have been found in several different positions (nucleotide no. C467, C586, C611, C737, C746, C830, and C1481) within the ndhB coding region. The rate of these alterations was deferentially edited across the three daylight periods. RNA editing rate of ndhB gene was highest at dawn, (87.5, 79.6, 78.5, 76, 68.6, 39.3, and 96.9%, respectively), less in midday (74.8, 54.1, 62.6, 47.4, 45.5, 47.4, and 93.4%, respectively), and less at pre-dusk (67, 52.6, 56.9, 40.1, 40.7, 33.2, and 90%, respectively), also all these sites were validated by qRT-PCR.

Conclusion

The differential editing of chloroplast ndhB gene across light periods may be led to a somehow relations between the RNA editing and control of photosynthesis.

High Level of Conservation of Mitochondrial RNA Editing Sites Among Four Populus Species

RNA editing occurs in the endosymbiont organelles of higher plants as C-to-U conversions of defined nucleotides. The availability of large quantities of RNA sequencing data makes it possible to identify RNA editing sites and to quantify their editing extent. We have investigated RNA editing in 34 protein-coding mitochondrial transcripts of four Populus species, a genus noteworthy for its remarkably small number of RNA editing sites compared to other angiosperms. 27 of these transcripts were subject to RNA editing in at least one species. In total, 355 RNA editing sites were identified with high confidence, their editing extents ranging from 10 to 100%. The most heavily edited transcripts were ccmB with the highest density of RNA editing sites (53.7 sites / kb) and ccmFn with the highest number of sites (39 sites). Most of the editing events are at position 1 or 2 of the codons, usually altering the encoded amino acid, and are highly conserved among the species, also with regard to their editing extent. However, one SNP was found in the newly sequenced and annotated mitochondrial genome of P. alba resulting in the loss of an RNA editing site compared to P. tremula and P. davidiana This SNP causes a C-to-T transition and an amino acid exchange from Ser to Phe, highlighting the widely discussed role of RNA editing in compensating mutations.

Comparative RNA editing profile of mitochondrial transcripts in cytoplasmic male sterile and fertile pigeonpea reveal significant changes at the protein level

RNA editing is a process which leads to post-transcriptional alteration of the nucleotide sequence of the corresponding mRNA molecule which may or may not lead to changes at the protein level. Apart from its role in providing variability at the transcript and protein levels, sometimes, such changes may lead to abnormal expression of the mitochondrial gene leading to a cytoplasmic male sterile phenotype. Here we report the editing status of 20 major mitochondrial transcripts in both male sterile (AKCMS11) and male fertile (AKPR303) pigeonpea genotypes. The validation of the predicted editing sites was done by mapping RNA-seq reads onto the amplified mitochondrial genes, and 165 and 159 editing sites were observed in bud tissues of the male sterile and fertile plant respectively. Among the resulting amino acid alterations, the most frequent one was the conversion of hydrophilic amino acids to hydrophobic. The alterations thus detected in our study indicates differential editing, but no major change in terms of the abnormal protein structure was detected. However, the above investigation provides an insight into the behaviour of pigeonpea mitochondrial genome in native and alloplasmic state and could hold clues in identification of editing factors and their role in adaptive evolution in pigeonpea.

RNA editing analysis of ATP synthase genes in the cotton cytoplasmic male sterile line H276A

Background

Pollen development is an energy-consuming process that particularly occurs during meiosis. Low levels of adenosine triphosphate (ATP) may cause cell death, resulting in CMS (cytoplasmic male sterility). DNA sequence differences in ATP synthase genes have been revealed between the N- and S-cytoplasms in the cotton CMS system. However, very few data are available at the RNA level. In this study, we compared five ATP synthase genes in the H276A, H276B and fertile F1 (H276A/H268) lines using RNA editing, RNA blotting and quantitative real time-PCR (qRT-PCR) to explore their contribution to CMS. A molecular marker for identifying male sterile cytoplasm (MSC) was also developed.

Results

RNA blotting revealed the absence of any novel orf for the ATP synthase gene sequence in the three lines. Forty-one RNA editing sites were identified in the coding sequences. RNA editing showed that proteins had 32.43% higher hydrophobicity and that 39.02% of RNA editing sites had proline converted to leucine. Two new stop codons were detected in atp6 and atp9 by RNA editing. Real-time qRT-PCR data showed that the atp1, atp6, atp8, and atp9 genes had substantially lower expression levels in H276A compared with those in H276B. By contrast, the expression levels of all five genes were increased in F1 (H276A/H268). Moreover, a molecular marker based on a 6-bp deletion upstream of atp8 in H276A was developed to identify male sterile cytoplasm (MSC) in cotton.

Conclusions

Our data substantially contributes to the understanding of the function of ATP synthase genes in cotton CMS. Therefore, we suggest that ATP synthase genes might be an indirect cause of cotton CMS. Further research is needed to investigate the relationship among ATP synthase genes in cotton CMS.

Electronic supplementary material

The online version of this article (10.1186/s40659-019-0212-0) contains supplementary material, which is available to authorized users.

RNA editing mutants as surrogates for mitochondrial SNP mutants

In terrestrial plants, RNA editing converts specific cytidines to uridines in mitochondrial and plastidic transcripts. Most of these events appear to be important for proper function of organellar encoded genes, since translated proteins from edited mRNAs show higher similarity with evolutionary conserved polypeptide sequences. So far about 100 nuclear encoded proteins have been characterized as RNA editing factors in plant organelles. Respective RNA editing mutants reduce or lose editing activity at different sites and display various macroscopic phenotypes from pale or albino in the case of chloroplasts to growth retardation or even embryonic lethality. Therefore, RNA editing mutants can be a useful resource of surrogate mutants for organellar encoded genes, especially for mitochondrially encoded genes that it is so far unfeasible to manipulate. However, connections between RNA editing defects and observed phenotypes in the mutants are often hard to elucidate, since RNA editing factors often target multiple RNA sites in different genes simultaneously. In this review article, we summarize the physiological aspects of respective RNA editing mutants and discuss them as surrogate mutants for functional analysis of mitochondrially encoded genes.

Towards a comprehensive picture of C-to-U RNA editing sites in angiosperm mitochondria

Our understanding of the dynamic and evolution of RNA editing in angiosperms is in part limited by the few editing sites identified to date. This study identified 10,217 editing sites from 17 diverse angiosperms. Our analyses confirmed the universality of certain features of RNA editing, and offer new evidence behind the loss of editing sites in angiosperms. RNA editing is a post-transcriptional process that substitutes cytidines (C) for uridines (U) in organellar transcripts of angiosperms. These substitutions mostly take place in mitochondrial messenger RNAs at specific positions called editing sites. By means of publicly available RNA-seq data, this study identified 10,217 editing sites in mitochondrial protein-coding genes of 17 diverse angiosperms. Even though other types of mismatches were also identified, we did not find evidence of non-canonical editing processes. The results showed an uneven distribution of editing sites among species, genes, and codon positions. The analyses revealed that editing sites were conserved across angiosperms but there were some species-specific sites. Non-synonymous editing sites were particularly highly conserved (~ 80%) across the plant species and were efficiently edited (80% editing extent). In contrast, editing sites at third codon positions were poorly conserved (~ 30%) and only partially edited (~ 40% editing extent). We found that the loss of editing sites along angiosperm evolution is mainly occurring by replacing editing sites with thymidines, instead of a degradation of the editing recognition motif around editing sites. Consecutive and highly conserved editing sites had been replaced by thymidines as result of retroprocessing, by which edited transcripts are reverse transcribed to cDNA and then integrated into the genome by homologous recombination. This phenomenon was more pronounced in eudicots, and in the gene cox1. These results suggest that retroprocessing is a widespread driving force underlying the loss of editing sites in angiosperm mitochondria.

Phylogenetic Relationships of the Fern Cyrtomium falcatum (Dryopteridaceae) from Dokdo Island Based on Chloroplast Genome Sequencing

Cyrtomium falcatum is a popular ornamental fern cultivated worldwide. Native to the Korean Peninsula, Japan, and Dokdo Island in the Sea of Japan, it is the only fern present on Dokdo Island. We isolated and characterized the chloroplast (cp) genome of C. falcatum, and compared it with those of closely related species. The genes trnV-GAC and trnV-GAU were found to be present within the cp genome of C. falcatum, whereas trnP-GGG and rpl21 were lacking. Moreover, cp genomes of Cyrtomium devexiscapulae and Adiantum capillus-veneris lack trnP-GGG and rpl21, suggesting these are not conserved among angiosperm cp genomes. The deletion of trnR-UCG, trnR-CCG, and trnSeC in the cp genomes of C. falcatum and other eupolypod ferns indicates these genes are restricted to tree ferns, non-core leptosporangiates, and basal ferns. The C. falcatum cp genome also encoded ndhF and rps7, with GUG start codons that were only conserved in polypod ferns, and it shares two significant inversions with other ferns, including a minor inversion of the trnD-GUC region and an approximate 3 kb inversion of the trnG-trnT region. Phylogenetic analyses showed that Equisetum was found to be a sister clade to Psilotales-Ophioglossales with a 100% bootstrap (BS) value. The sister relationship between Pteridaceae and eupolypods was also strongly supported by a 100% BS, but Bayesian molecular clock analyses suggested that C. falcatum diversified in the mid-Paleogene period (45.15 ± 4.93 million years ago) and might have moved from Eurasia to Dokdo Island.

Structural and functional analyses of Barth syndrome-causing mutations and alternative splicing in the tafazzin acyltransferase domain

Tafazzin is a mitochondrial phospholipid transacylase, and its mutations cause Barth syndrome (BTHS). Human tafazzin gene produces four distinct alternatively spliced transcripts. To understand the molecular mechanisms of tafazzin deficiency, we performed an atomic resolution analysis of the influence of the BTHS mutations and of alternative splicing on the structure and function of tafazzin. From the three-dimensional (3D) homology modeling of tafazzin, we identified candidate amino acid residues that contribute to cardiolipin binding and to mitochondrial membrane associations that facilitate acyl-transfer reactions. Primate specific exon 5, which is alternatively spliced, is predicted to correspond to an intrinsically unstructured region in the protein. We proposed that this region should change the substrate-binding affinity and/or contribute to primate-specific molecular interactions. Exon 7, another alternatively spliced exon, encodes a region forming a part of the putative substrate-binding cleft, suggesting that the gene products lacking exon 7 will lose their substrate-binding ability. We demonstrate a clear localization of the BTHS mutations at residues responsible for membrane association, substrate binding, and the conformational stability of tafazzin. These findings provide new insights into the function of defective tafazzin and the pathogenesis of BTHS at the level of protein 3D structure and the evolution of alternatively spliced exons in primates.

Massive gene transfer and extensive RNA editing of a symbiotic dinoflagellate plastid genome

Genome sequencing of Symbiodinium minutum revealed that 95 of 109 plastid-associated genes have been transferred to the nuclear genome and subsequently expanded by gene duplication. Only 14 genes remain in plastids and occur as DNA minicircles. Each minicircle (1.8–3.3 kb) contains one gene and a conserved noncoding region containing putative promoters and RNA-binding sites. Nine types of RNA editing, including a novel G/U type, were discovered in minicircle transcripts but not in genes transferred to the nucleus. In contrast to DNA editing sites in dinoflagellate mitochondria, which tend to be highly conserved across all taxa, editing sites employed in DNA minicircles are highly variable from species to species. Editing is crucial for core photosystem protein function. It restores evolutionarily conserved amino acids and increases peptidyl hydropathy. It also increases protein plasticity necessary to initiate photosystem complex assembly.

Photosynthetic gene expression in higher plants

Within the chloroplasts of higher plants and algae, photosynthesis converts light into biological energy, fueling the assimilation of atmospheric carbon dioxide into biologically useful molecules. Two major steps, photosynthetic electron transport and the Calvin-Benson cycle, require many gene products encoded from chloroplast as well as nuclear genomes. The expression of genes in both cellular compartments is highly dynamic and influenced by a diverse range of factors. Light is the primary environmental determinant of photosynthetic gene expression. Working through photoreceptors such as phytochrome, light regulates photosynthetic genes at transcriptional and posttranscriptional levels. Other processes that affect photosynthetic gene expression include photosynthetic activity, development, and biotic and abiotic stress. Anterograde (from nucleus to chloroplast) and retrograde (from chloroplast to nucleus) signaling insures the highly coordinated expression of the many photosynthetic genes between these different compartments. Anterograde signaling incorporates nuclear-encoded transcriptional and posttranscriptional regulators, such as sigma factors and RNA-binding proteins, respectively. Retrograde signaling utilizes photosynthetic processes such as photosynthetic electron transport and redox signaling to influence the expression of photosynthetic genes in the nucleus. The basic C3 photosynthetic pathway serves as the default form used by most of the plant species on earth. High temperature and water stress associated with arid environments have led to the development of specialized C4 and CAM photosynthesis, which evolved as modifications of the basic default expression program. The goal of this article is to explain and summarize the many gene expression and regulatory processes that work together to support photosynthetic function in plants.

Identification of RNA editing sites in cotton (Gossypium hirsutum) chloroplasts and editing events that affect secondary and three-dimensional protein structures

RNA editing can alter individual nucleotides in primary transcripts, which can cause the amino acids encoded by edited RNA to deviate from the ones predicted from the DNA template. We investigated RNA editing sites of protein-coding genes from the chloroplast genome of cotton. Fifty-four editing sites were identified in 27 transcripts, which is the highest editing frequency found until now in angiosperms. All these editing sites were C-to-U conversion, biased toward ndh genes and U_A context. Examining published editotypes in various angiosperms, we found that RNA editing mostly converts amino acid from hydrophilic to hydrophobic and restores evolutionary conserved amino acids. Using bioinformatics to analyze the effect of editing events on protein secondary and three-dimensional structures, we found that 21 editing sites can affect protein secondary structures and seven editing sites can alter three-dimensional protein structures. These results imply that 24 editing sites in cotton chloroplast transcripts may play an important role in their protein structures and functions.

The nucleocytoplasmic conflict, a driving force for the emergence of plant organellar RNA editing

RNA editing challenges the central dogma of molecular biology by changing the genetic information at the transcript level. In plant organelles, RNAs are modified by deamination of some specific cytosine residues, but the origin of this process remains puzzling. Different from the generally accepted neutral model to explain the emergence of RNA editing in plant organelles, we propose a new hypothesis based on the nucleocytoplasmic conflict theory. We assume that mutations in organellar genomes arose first and spread into the population provided they increased the transmission of their own maternally inherited genome. RNA editing appeared subsequently as a nuclear-encoded correction mechanism to restore the transmission of the nuclear genome. In plants, a well-known consequence of the nucleocytoplasmic conflict is cytoplasmic male sterility (CMS) which is counteracted by the emergence of fertility restorer genes (Rf) belonging to the pentatricopeptide repeat (PPR) protein family. Interestingly, RNA-editing deficiency can lead to CMS, and it now clearly appears that PPR proteins are major players in RNA editing. This striking similarity between the mechanisms of fertility restoration and RNA editing can be explained if both reactions are the consequence of the same driving force, the nucleocytoplasmic conflict. Similarly, the prevalence of RNA editing in eukaryotic organellar genomes could also be a consequence of the genetic antagonism between organellar and nuclear genomes.

RNA editing sites exist in protein-coding genes in the chloroplast genome of Cycas taitungensis

RNA editing is a post-transcriptional process that results in modifications of ribonucleotides at specific locations. In land plants editing can occur in both mitochondria and chloroplasts and most commonly involves C-to-U changes, especially in seed plants. Using prediction and experimental determination, we investigated RNA editing in 40 protein-coding genes from the chloroplast genome of Cycas taitungensis. A total of 85 editing sites were identified in 25 transcripts. Comparison analysis of the published editotypes of these 25 transcripts in eight species showed that RNA editing events gradually disappear during plant evolution. The editing in the first and third codon position disappeared quicker than that in the second codon position. ndh genes have the highest editing frequency while serine and proline codons were more frequently edited than the codons of other amino acids. These results imply that retained RNA editing sites have imbalanced distribution in genes and most of them may function by changing protein structure or interaction. Mitochondrion protein-coding genes have three times the editing sites compared with chloroplast genes of Cycas, most likely due to slower evolution speed.

Are substitution rates and RNA editing correlated?

Background

RNA editing is a post-transcriptional process that, in seed plants, involves a cytosine to uracil change in messenger RNA, causing the translated protein to differ from that predicted by the DNA sequence. RNA editing occurs extensively in plant mitochondria, but large differences in editing frequencies are found in some groups. The underlying processes responsible for the distribution of edited sites are largely unknown, but gene function, substitution rate, and gene conversion have been proposed to influence editing frequencies.

Results

We studied five mitochondrial genes in the monocot order Alismatales, all showing marked differences in editing frequencies among taxa. A general tendency to lose edited sites was observed in all taxa, but this tendency was particularly strong in two clades, with most of the edited sites lost in parallel in two different areas of the phylogeny. This pattern is observed in at least four of the five genes analyzed. Except in the groups that show an unusually low editing frequency, the rate of C-to-T changes in edited sites was not significantly higher that in non-edited 3^rd codon positions. This may indicate that selection is not actively removing edited sites in nine of the 12 families of the core Alismatales. In all genes but ccmB, a significant correlation was found between frequency of change in edited sites and synonymous substitution rate. In general, taxa with higher substitution rates tend to have fewer edited sites, as indicated by the phylogenetically independent correlation analyses. The elimination of edited sites in groups that lack or have reduced levels of editing could be a result of gene conversion involving a cDNA copy (retroprocessing). If so, this phenomenon could be relatively common in the Alismatales, and may have affected some groups recurrently. Indirect evidence of retroprocessing without a necessary correlation with substitution rate was found mostly in families Alismataceae and Hydrocharitaceae (e.g., groups that suffered a rapid elimination of all their edited sites, without a change in substitution rate).

Conclusions

The effects of substitution rate, selection, and/or gene conversion on the dynamics of edited sites in plant mitochondria remain poorly understood. Although we found an inverse correlation between substitution rate and editing frequency, this correlation is partially obscured by gene retroprocessing in lineages that have lost most of their edited sites. The presence of processed paralogs in plant mitochondria deserves further study, since most evidence of their occurrence is circumstantial.

Extensive loss of RNA editing sites in rapidly evolving Silene mitochondrial genomes: selection vs. retroprocessing as the driving force

Theoretical arguments suggest that mutation rates influence the proliferation and maintenance of RNA editing. We identified RNA editing sites in five species within the angiosperm genus Silene that exhibit highly divergent mitochondrial mutation rates. We found that mutational acceleration has been associated with rapid loss of mitochondrial editing sites. In contrast, we did not find a significant difference in the frequency of editing in chloroplast genes, which lack the mutation rate variation observed in the mitochondrial genome. As found in other angiosperms, the rate of substitution at RNA editing sites in Silene greatly exceeds the rate at synonymous sites, a pattern that has previously been interpreted as evidence for selection against RNA editing. Alternatively, we suggest that editing sites may experience higher rates of C-to-T mutation than other portions of the genome. Such a pattern could be caused by gene conversion with reverse-transcribed mRNA (i.e., retroprocessing). If so, the genomic distribution of RNA editing site losses in Silene suggests that such conversions must be occurring at a local scale such that only one or two editing sites are affected at a time. Because preferential substitution at editing sites appears to occur in angiosperms regardless of the mutation rate, we conclude that mitochondrial rate accelerations within Silene have "fast-forwarded" a preexisting pattern but have not fundamentally changed the evolutionary forces acting on RNA editing sites.

Large-scale detection and analysis of RNA editing in grape mtDNA by RNA deep-sequencing

RNA editing is a widespread post-transcriptional molecular phenomenon that can increase proteomic diversity, by modifying the sequence of completely or partially non-functional primary transcripts, through a variety of mechanistically and evolutionarily unrelated pathways. Editing by base substitution has been investigated in both animals and plants. However, conventional strategies based on directed Sanger sequencing are time-consuming and effectively preclude genome wide identification of RNA editing and assessment of partial and tissue-specific editing sites. In contrast, the high-throughput RNA-Seq approach allows the generation of a comprehensive landscape of RNA editing at the genome level. Short reads from Solexa/Illumina GA and ABI SOLiD platforms have been used to investigate the editing pattern in mitochondria of Vitis vinifera providing significant support for 401 C-to-U conversions in coding regions and an additional 44 modifications in non-coding RNAs. Moreover, 76% of all C-to-U conversions in coding genes represent partial RNA editing events and 28% of them were shown to be significantly tissue specific. Solexa/Illumina and SOLiD platforms showed different characteristics with respect to the specific issue of large-scale editing analysis, and the combined approach presented here reduces the false positive rate of discovery of editing events.

RNA editing of 10 Didymium iridis mitochondrial genes and comparison with the homologous genes in Physarum polycephalum

Regions of the Didymium iridis mitochondrial genome were identified with similarity to typical mitochondrial genes; however, these regions contained numerous stop codons. We used RT-PCR and DNA sequencing to determine whether, through RNA editing, these regions were transcribed into mRNAs that could encode functional proteins. Ten putative gene regions were examined: atp1, atp6, atp8, atp9, cox1, cox2, cytb, nad4L, nad6, and nad7. The cDNA sequences of each gene could encode a functional mitochondrial protein that was highly conserved compared with homologous genes. The type of editing events and editing sequence features were very similar to those observed in the homologous genes of Physarum polycephalum, though the actual editing locations showed a variable degree of conservation. Edited sites were compared with encoded sites in D. iridis and P. polycephalum for all 10 genes. Edited sequence for a portion of the cox1 gene was available for six myxomycetes, which, when compared, showed a high degree of conservation at the protein level. Different types of editing events showed varying degrees of site conservation with C-to-U base changes being the least conserved. Several aspects of single C insertion editing events led to the preferential creation of hydrophobic amino acid codons that may help to minimize adverse effects on the resulting protein structure.

RESOPS: a database for analyzing the correspondence of RNA editing sites to protein three-dimensional structures

Transcripts from mitochondrial and chloroplast DNA of land plants often undergo cytidine to uridine conversion-type RNA editing events. RESOPS is a newly built database that specializes in displaying RNA editing sites of land plant organelles on protein three-dimensional (3D) structures to help elucidate the mechanisms of RNA editing for gene expression regulation. RESOPS contains the following information: unedited and edited cDNA sequences with notes for the target nucleotides of RNA editing, conceptual translation from the edited cDNA sequence in pseudo-UniProt format, a list of proteins under the influence of RNA editing, multiple amino acid sequence alignments of edited proteins, the location of amino acid residues coded by codons under the influence of RNA editing in protein 3D structures and the statistics of biased distributions of the edited residues with respect to protein structures. Most of the data processing procedures are automated; hence, it is easy to keep abreast of updated genome and protein 3D structural data. In the RESOPS database, we clarified that the locations of residues switched by RNA editing are significantly biased to protein structural cores. The integration of different types of data in the database also help advance the understanding of RNA editing mechanisms. RESOPS is accessible at http://cib.cf.ocha.ac.jp/RNAEDITING/.

Rice OGR1 encodes a pentatricopeptide repeat-DYW protein and is essential for RNA editing in mitochondria

RNA editing is the alteration of RNA sequences via insertion, deletion and conversion of nucleotides. In flowering plants, specific cytidine residues of RNA transcribed from organellar genomes are converted into uridines. Approximately 35 editing sites are present in the chloroplasts of higher plants; six pentatricopeptide repeat genes involved in RNA editing have been identified in Arabidopsis. However, although approximately 500 editing sites are found in mitochondrial RNAs of flowering plants, only one gene in Arabidopsis has been reported to be involved in such editing. Here, we identified rice mutants that are defective in seven specific RNA editing sites on five mitochondrial transcripts. Their various phenotypes include delayed seed germination, retarded growth, dwarfism and sterility. Mutant seeds from heterozygous plants are opaque. This mutation, named opaque and growth retardation 1 (ogr1), was generated by T-DNA insertion into a gene that encodes a pentatricopeptide repeat protein containing the DYW motif. The OGR1-sGFP fusion protein is localized to mitochondria. Ectopic expression of OGR1 in the mutant complements the altered phenotypes. We conclude that OGR1 is essential for RNA editing in rice mitochondria and is required for normal growth and development.

A comparative genomics approach identifies a PPR-DYW protein that is essential for C-to-U editing of the Arabidopsis chloroplast accD transcript

Several nuclear-encoded proteins containing pentatricopeptide repeat (PPR) motifs have previously been identified to be trans-factors essential for particular chloroplast RNA editing events through analysis of mutants affected in chloroplast biogenesis or function. Other PPR genes are known to encode proteins involved in other aspects of organelle RNA metabolism. A function has not been assigned to most members of the large plant PPR gene family. Arabidopsis and rice each contain over 400 PPR genes, of which about a fifth exhibit a C-terminal DYW domain. We describe here a comparative genomics approach that will facilitate identification of the role of RNA-binding proteins in organelle RNA metabolism. We have implemented this strategy to identify an Arabidopsis nuclear-encoded gene RARE1 that is required for editing of the chloroplast accD transcript. RARE1 carries 15 PPR motifs, an E/E+ and a DYW domain, whereas previously reported editing factors CRR4, CRR21, and CLB19 lack a DYW domain. The accD gene encodes the beta carboxyltransferase subunit of acetyl coA carboxylase, which catalyzes the first step in fatty acid biosynthesis in chloroplasts. Despite a lack of accD C794 editing and lack of restoration of an evolutionarily conserved leucine residue in the beta carboxyltransferase protein, rare1 mutants are unexpectedly robust and reproduce under growth room conditions. Previously the serine-to-leucine alteration caused by editing was deemed essential in the light of the finding that a recombinantly expressed "unedited" form of the pea acetyl coA carboxylase was catalytically inactive.