RNA editing machinery in plant organelles

Dissecting the molecular puzzle of the editosome core in Arabidopsis organelles

Over the last decade, the composition of the C-to-U RNA editing complex in embryophyte organelles has turned out to be much more complex than first expected. While PPR proteins were initially thought to act alone, significant evidences have clearly depicted a sophisticated mechanism with numerous protein-protein interaction involving PPR and non-PPR proteins. Moreover, the identification of specific functional partnership between PPRs also suggests that, in addition to the highly specific PPRs directly involved in the RNA target recognition, non-RNA-specific ones are required. Although some of them, such as DYW1 and DYW2, were shown to be the catalytic domains of the editing complex, the molecular function of others, such as NUWA, remains elusive. It was suggested that they might stabilize the complex by acting as a scaffold. We here performed functional complementation of the crr28-2 mutant with truncated CRR28 proteins mimicking PPR without the catalytic domain and show that they exhibit a specific dependency to one of the catalytic proteins DYW1 or DYW2. Moreover, we also characterized the role of the PPR NUWA in the editing reaction and show that it likely acts as a scaffolding factor. NUWA is no longer required for efficient editing of the CLB19 editing sites once this RNA specific PPR is fused to the DYW catalytic domain of its partner DYW2. Altogether, our results strongly support a flexible, evolutive and resilient editing complex in which RNA binding activity, editing activity and stabilization/scaffolding function can be provided by one or more PPRs.

Insights into structure, codon usage, repeats, and RNA editing of the complete mitochondrial genome of Perilla frutescens (Lamiaceae)

Perilla frutescens (L.) Britton, a member of the Lamiaceae family, stands out as a versatile plant highly valued for its unique aroma and medicinal properties. Additionally, P. frutescens seeds are rich in Îś-linolenic acid, holding substantial economic importance. While the nuclear and chloroplast genomes of P. frutescens have already been documented, the complete mitochondrial genome sequence remains unreported. To this end, the sequencing, annotation, and assembly of the entire Mitochondrial genome of P. frutescens were hereby conducted using a combination of Illumina and PacBio data. The assembled P. frutescens mitochondrial genome spanned 299,551 bp and exhibited a typical circular structure, involving a GC content of 45.23%. Within the genome, a total of 59 unique genes were identified, encompassing 37 protein-coding genes, 20 tRNA genes, and 2 rRNA genes. Additionally, 18 introns were observed in 8 protein-coding genes. Notably, the codons of the P. frutescens mitochondrial genome displayed a notable A/T bias. The analysis also revealed 293 dispersed repeat sequences, 77 simple sequence repeats (SSRs), and 6 tandem repeat sequences. Moreover, RNA editing sites preferentially produced leucine at amino acid editing sites. Furthermore, 70 sequence fragments (12,680 bp) having been transferred from the chloroplast to the mitochondrial genome were identified, accounting for 4.23% of the entire mitochondrial genome. Phylogenetic analysis indicated that among Lamiaceae plants, P. frutescens is most closely related to Salvia miltiorrhiza and Platostoma chinense. Meanwhile, inter-species Ka/Ks results suggested that Ka/Ks < 1 for 28 PCGs, indicating that these genes were evolving under purifying selection. Overall, this study enriches the mitochondrial genome data for P. frutescens and forges a theoretical foundation for future molecular breeding research.

Is RNA editing truly absent in the complex thalloid liverworts (Marchantiopsida)? Evidence of extensive RNA editing from Cyathodium cavernarum

RNA editing is a crucial modification in plants' organellar transcripts that converts cytidine to uridine (C-to-U; and sometimes uridine to cytidine) in RNA molecules. This post-transcriptional process is controlled by the PLS-class protein with a DYW domain, which belongs to the pentatricopeptide repeat (PPR) protein family. RNA editing is widespread in land plants; however, complex thalloid liverworts (Marchantiopsida) are the only group reported to lack both RNA editing and DYW-PPR protein. The liverwort Cyathodium cavernarum (Marchantiopsida, Cyathodiaceae), typically found in cave habitats, was newly found to have 129 C-to-U RNA editing sites in its chloroplast and 172 sites in its mitochondria. The Cyathodium genus, specifically C. cavernarum, has a large number of PPR editing factor genes, including 251 DYW-type PPR proteins. These DYW-type PPR proteins may be responsible for C-to-U RNA editing in C. cavernarum. Cyathodium cavernarum possesses both PPR DYW proteins and RNA editing. Our analysis suggests that the remarkable RNA editing capability of C. cavernarum may have been acquired alongside the emergence of DYW-type PPR editing factors. These findings provide insight into the evolutionary pattern of RNA editing in land plants.

A P-type pentatricopeptide repeat protein ZmRF5 promotes 5' region partial cleavages of atp6c transcripts to restore the fertility of CMS-C maize by recruiting a splicing factor

A fast evolution within mitochondria genome(s) often generates discords between nuclear and mitochondria, which is manifested as cytoplasmic male sterility (CMS) and fertility restoration (Rf) system. The maize CMS-C trait is regulated by the chimeric mitochondrial gene, atp6c, and can be recovered by the restorer gene ZmRf5. Through positional cloning in this study, we identified the nuclear restorer gene, ZmRf5, which encodes a P-type pentatricopeptide repeat (PPR) family protein. The over-expression of ZmRf5 brought back the fertility to CMS-C plants, whereas its genomic editing by CRISPR/Cas9 induced abortive pollens in the restorer line. ZmRF5 is sorted to mitochondria, and recruited RS31A, a splicing factor, through MORF8 to form a cleaving/restoring complex, which promoted the cleaving of the CMS-associated transcripts atp6c by shifting the major cleavage site from 480th nt to 344 th nt for fast degradation, and preserved just right amount of atp6c RNA for protein translation, providing adequate ATP6C to assembly complex V, thus restoring male fertility. Interestingly, ATP6C in the sterile line CMo17A, with similar cytology and physiology changes to YU87-1A, was accumulated much less than it in NMo17B, exhibiting a contrary trend in the YU87-1 nuclear genome previously reported, and was restored to normal level in the presence of ZmRF5. Collectively these findings unveil a new molecular mechanism underlying fertility restoration by which ZmRF5 cooperates with MORF8 and RS31A to restore CMS-C fertility in maize, complemented and perfected the sterility mechanism, and enrich the perspectives on communications between nucleus and mitochondria.

PlantC2U: deep learning of cross-species sequence landscapes predicts plastid C-to-U RNA editing in plants

In plants, C-to-U RNA editing mainly occurs in plastid and mitochondrial transcripts, which contributes to a complex transcriptional regulatory network. More evidence reveals that RNA editing plays critical roles in plant growth and development. However, accurate detection of RNA editing sites using transcriptome sequencing data alone is still challenging. In the present study, we develop PlantC2U, which is a convolutional neural network, to predict plastid C-to-U RNA editing based on the genomic sequence. PlantC2U achieves >95% sensitivity and 99% specificity, which outperforms the PREPACT tool, random forests, and support vector machines. PlantC2U not only further checks RNA editing sites from transcriptome data to reduce possible false positives, but also assesses the effect of different mutations on C-to-U RNA editing based on the flanking sequences. Moreover, we found the patterns of tissue-specific RNA editing in the mangrove plant Kandelia obovata, and observed reduced C-to-U RNA editing rates in the cold stress response of K. obovata, suggesting their potential regulatory roles in plant stress adaptation. In addition, we present RNAeditDB, available online at https://jasonxu.shinyapps.io/RNAeditDB/. Together, PlantC2U and RNAeditDB will help researchers explore the RNA editing events in plants and thus will be of broad utility for the plant research community.

Conservation of the moss RNA editing factor PPR78 despite the loss of its known cytidine-to-uridine editing sites is explained by a hidden extra target

Cytidine (C)-to-uridine (U) RNA editing in plant organelles relies on specific RNA-binding pentatricopeptide repeat (PPR) proteins. In the moss Physcomitrium patens, all such RNA editing factors feature a C-terminal DYW domain that acts as the cytidine deaminase for C-to-U conversion. PPR78 of Physcomitrium targets 2 mitochondrial editing sites, cox1eU755SL and rps14eU137SL. Remarkably, the latter is edited to highly variable degrees in different mosses. Here, we aimed to unravel the coevolution of PPR78 and its 2 target sites in mosses. Heterologous complementation in a Physcomitrium knockout line revealed that the variable editing of rps14eU137SL depends on the PPR arrays of different PPR78 orthologues but not their C-terminal domains. Intriguingly, PPR78 has remained conserved despite the simultaneous loss of editing at both known targets among Hypnales (feather mosses), suggesting it serves an additional function. Using a recently established RNA editing assay in Escherichia coli, we confirmed site-specific RNA editing by PPR78 in the bacterium and identified 4 additional off-targets in the bacterial transcriptome. Based on conservation profiles, we predicted ccmFNeU1465RC as a candidate editing target of PPR78 in moss mitochondrial transcriptomes. We confirmed editing at this site in several mosses and verified that PPR78 targets ccmFNeU1465RC in the bacterial editing system, explaining the conservation and functional adaptation of PPR78 during moss evolution.

The Complete Mitochondrial Genome of Paeonia lactiflora Pall. (Saxifragales: Paeoniaceae): Evidence of Gene Transfer from Chloroplast to Mitochondrial Genome

Paeonia lactiflora (P. lactiflora), a perennial plant renowned for its medicinal roots, provides a unique case for studying the phylogenetic relationships of species based on organelle genomes, as well as the transference of DNA across organelle genomes. In order to investigate this matter, we sequenced and characterized the mitochondrial genome (mitogenome) of P. lactiflora. Similar to the chloroplast genome (cpgenome), the mitogenome of P. lactiflora extends across 181,688 base pairs (bp). Its unique quadripartite structure results from a pair of extensive inverted repeats, each measuring 25,680 bp in length. The annotated mitogenome includes 27 protein-coding genes, 37 tRNAs, 8 rRNAs, and two pseudogenes (rpl5, rpl16). Phylogenetic analysis was performed to identify phylogenetic trees consistent with Paeonia species phylogeny in the APG Ⅳ system. Moreover, a total of 12 MTPT events were identified and 32 RNA editing sites were detected during mitogenome analysis of P. lactiflora. Our research successfully compiled and annotated the mitogenome of P. lactiflora. The study provides valuable insights regarding the taxonomic classification and molecular evolution within the Paeoniaceae family.

RNA editing events and expression profiles of mitochondrial protein-coding genes in the endemic and endangered medicinal plant, Corydalis saxicola

Corydalis saxicola, an endangered medicinal plant endemic to karst habitats, is widely used in Traditional Chinese Medicine to treat hepatitis, abdominal pain, bleeding hemorrhoids and other conditions. However, to date, the mitochondrial (mt) genome of C. saxicola has not been reported, which limits our understanding of the genetic and biological mechanisms of C. saxicola. Here, the mt genome of C. saxicola was assembled by combining the Nanopore and Illumina reads. The mt genome of C. saxicola is represented by a circular chromosome which is 587,939 bp in length, with an overall GC content of 46.50%. 40 unique protein-coding genes (PCGs), 22 tRNA genes and three rRNA genes were identified. Codon usage of the PCGs was investigated and 167 simple sequence repeats were identified. Twelve homologous fragments were identified between the mt and ct genomes of C. saxicola, accounting for 1.04% of the entire mt genome. Phylogenetic examination of the mt genomes of C. saxicola and 30 other taxa provided an understanding of their evolutionary relationships. We also predicted 779 RNA editing sites in 40 C. saxicola mt PCGs and successfully validated 506 (65%) of these using PCR amplification and Sanger sequencing. In addition, we transcriptionally profiled 24 core mt PCGs in C. saxicola roots treated with different concentrations of CaCl2, as well as in other organs. These investigations will be useful for effective utilization and molecular breeding, and will also provide a reference for further studies of the genus Corydalis.

Chloroplast gene expression: Recent advances and perspectives

Chloroplasts evolved from an ancient cyanobacterial endosymbiont more than 1.5 billion years ago. During subsequent coevolution with the nuclear genome, the chloroplast genome has remained independent, albeit strongly reduced, with its own transcriptional machinery and distinct features, such as chloroplast-specific innovations in gene expression and complicated post-transcriptional processing. Light activates the expression of chloroplast genes via mechanisms that optimize photosynthesis, minimize photodamage, and prioritize energy investments. Over the past few years, studies have moved from describing phases of chloroplast gene expression to exploring the underlying mechanisms. In this review, we focus on recent advances and emerging principles that govern chloroplast gene expression in land plants. We discuss engineering of pentatricopeptide repeat proteins and its biotechnological effects on chloroplast RNA research; new techniques for characterizing the molecular mechanisms of chloroplast gene expression; and important aspects of chloroplast gene expression for improving crop yield and stress tolerance. We also discuss biological and mechanistic questions that remain to be answered in the future.

Differential adaptive RNA editing signals between insects and plants revealed by a new measurement termed haplotype diversity

BACKGROUND

C-to-U RNA editing in plants is believed to confer its evolutionary adaptiveness by reversing unfavorable DNA mutations. This "restorative hypothesis" has not yet been tested genome-wide. In contrast, A-to-I RNA editing in insects like Drosophila and honeybee is already known to benefit the host by increasing proteomic diversity in a spatial-temporal manner (namely "diversifying hypothesis").

METHODS

We profiled the RNA editomes of multiple tissues of Arabidopsis thaliana, Drosophila melanogaster, and Apis melifera. We unprecedentedly defined the haplotype diversity (HD) of RNA molecules based on nonsynonymous editing events (recoding sites).

RESULTS

Signals of adaptation is confirmed in Arabidopsis by observing higher frequencies and levels at nonsynonymous editing sites over synonymous sites. Compared to A-to-I recoding sites in Drosophila, the C-to-U recoding sites in Arabidopsis show significantly lower HD, presumably due to the stronger linkage between C-to-U events.

CONCLUSIONS

C-to-U RNA editing in Arabidopsis is adaptive but it is not designed for diversifying the proteome like A-to-I editing in Drosophila. Instead, C-to-U recoding sites resemble DNA mutations. Our observation supports the restorative hypothesis of plant C-to-U editing which claims that editing is used for fixing unfavorable genomic sequences.

Deciphering the enigma of RNA editing in the ATP1_alpha subunit of ATP synthase in Triticum aestivum

There is evidence that RNA editing is related to plant cellular stress as well as electron transport organelles, such as mitochondria. The mitochondrial atp1 gene encodes the alpha-subunit of Atp synthase. Control as well as two periods of drought stress treatments were analyzed in the cDNAs generated from the mitochondrial atp1 gene of two cultivars of Triticum aestivum [Giza 168 (G168) and Gemmiza 10 (GM10)]. Following RNA-seq data assembly, atp1 cDNAs from the control (acc. no. OQ129415), 2-hour (acc. no. OQ129416), and 12-hour (acc. no. OQ129417) time points of the T. aestivum cultivar G168 were obtained. Control (acc. no. OQ129419), 2-hour (acc. no. OQ129420), and 12-hour (acc. no. OQ129421) samples all included reconstructed atp1 transcripts from Gemmiza 10. Atp1 transcripts were assembled using the wheat atp1 gene (acc. no. NC_036024). RNA-seq raw data was utilized to identify 11 RNA editing sites in atp1 in the tolerant cultivar Giza168 and 6 in the sensitive cultivar Gemmiza10. The significant difference in RNA editing observed between control and drought stress conditions in sites led to synonymous amino acids. This led to no change in tertiary structure between tolerant and sensitive cultivars. But the change was focused between produced protein and its correspondence sequence on DNA.

Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Two Species of Calla Lilies (Zantedeschia, Araceae)

In this study, the mitochondrial genomes of two calla species, Zantedeschia aethiopica Spreng. and Zantedeschia odorata Perry., were assembled and compared for the first time. The Z. aethiopica mt genome was assembled into a single circular chromosome, measuring 675,575 bp in length with a 45.85% GC content. In contrast, the Z. odorata mt genome consisted of bicyclic chromosomes (chromosomes 1 and 2), measuring 719,764 bp and exhibiting a 45.79% GC content. Both mitogenomes harbored similar gene compositions, with 56 and 58 genes identified in Z. aethiopica and Z. odorata, respectively. Analyses of codon usage, sequence repeats, gene migration from chloroplast to mitochondrial, and RNA editing were conducted for both Z. aethiopica and Z. odorata mt genomes. Phylogenetic examination based on the mt genomes of these two species and 30 other taxa provided insights into their evolutionary relationships. Additionally, the core genes in the gynoecium, stamens, and mature pollen grains of the Z. aethiopica mt genome were investigated, which revealed maternal mitochondrial inheritance in this species. In summary, this study offers valuable genomic resources for future research on mitogenome evolution and the molecular breeding of calla lily.

The ATP Synthase γ Subunit ATPC1 Regulates RNA Editing in Chloroplasts

RNA editing is the process of modifying RNA molecules by inserting, deleting, or substituting nucleotides. In flowering plants, RNA editing occurs predominantly in RNAs encoded by the organellar genomes of mitochondria and chloroplasts, and the main type of editing involves the substitution of cytidine with uridine at specific sites. Abnormal RNA editing in plants can affect gene expression, organelle function, plant growth, and reproduction. In this study, we report that ATPC1, the gamma subunit of ATP synthase in Arabidopsis chloroplasts, has an unexpected role in the regulation of editing at multiple sites of plastid RNAs. The loss of function of ATPC1 severely arrests chloroplast development, causing a pale-green phenotype and early seedling lethality. Disruption of ATPC1 increases the editing of matK-640, rps12-i-58, atpH-3'UTR-13210, and ycf2-as-91535 sites while decreasing the editing of rpl23-89, rpoA-200, rpoC1-488, and ndhD-2 sites. We further show that ATPC1 participates in RNA editing by interacting with known multiple-site chloroplast RNA editing factors, including MORFs, ORRM1, and OZ1. The transcriptome in the atpc1 mutant is profoundly affected, with a pattern of defective expression of chloroplast development-related genes. These results reveal that the ATP synthase γ subunit ATPC1 is involved in multiple-site RNA editing in Arabidopsis chloroplasts.

The PPR protein RARE1-mediated editing of chloroplast accD transcripts is required for fatty acid biosynthesis and heat tolerance in Arabidopsis

It has been reported that Arabidopsis chloroplast accD transcripts undergo RNA editing and that loss of accD-C794 RNA editing does not affect plant growth under normal conditions. To date, the exact biological role of accD-C794 editing has remained elusive. Here, we reveal an unexpected role for accD-C794 editing in response to heat stress. Loss of accD-C794 editing results in a yellow and dwarf phenotype with decreased chloroplast gene expression under heat stress, and artificial improvement of C794-edited accD gene expression enhances heat tolerance in Arabidopsis. These data suggest that accD-C794 editing confers heat tolerance in planta. We also found that treatment with the product of acetyl coenzyme A carboxylase (ACCase) could allay mutant phenotypic characteristics and showed that a mutation in the CAC3 gene for the α-subunit of ACCase was associated with dwarfism under heat stress. These observations indicate that defective accD-C794 editing may be intrinsic to reduced ACCase activity, thereby contributing to heat sensitivity. ACCase catalyzes the committed step of de novo fatty acid (FA) biosynthesis. FA content analysis revealed that unsaturated oleic (C18:1) and linoleic acids (C18:2) were low in the accD-C794 editing-defective mutant but high in the C794-edited accD-overexpressing plants compared with the wild type. Supplying exogenous C18:1 and C18:2 could rescue the mutant phenotype, suggesting that these FAs play an essential role in tolerance to heat stress. Transmission electron microscopy observations showed that heat stress seriously affected the membrane architecture in accD editing-defective mutants but not in accD-overexpressing plants. These results provide the first evidence that accD-C794 editing regulates FA biosynthesis for maintenance of membrane structural homeostasis under heat stress.

Complete sequence and comparative analysis of the mitochondrial genome of the rare and endangered Clematis acerifolia, the first clematis mitogenome to provide new insights into the phylogenetic evolutionary status of the genus

Clematis is one of the large worldwide genera of the Ranunculaceae Juss. Family, with high ornamental and medicinal value. China is the modern distribution centre of Clematis with abundant natural populations. Due to the complexity and high morphological diversity of Clematis, the genus is difficult to classify systematically, and in particular, the phylogenetic position of the endangered Clematis acerifolia is highly controversial. The use of the mitochondrial complete genome is a powerful molecular method that is frequently used for inferring plants phylogenies. However, studies on Clematis mitogenome are rare, thus limiting our full understanding of its phylogeny and genome evolution. Here, we sequenced and annotated the C. acerifolia mt genome using Illumina short- and Nanopore long-reads, characterized the species first complete mitogenome, and performed a comparative phylogenetic analysis with its close relatives. The total length of the C. acerifolia mitogenome is 698,247 bp and the main structure is multi-branched (linear molecule 1 and circular molecule 2). We annotated 55 genes, including 35 protein-coding, 17 tRNA, and 3 rRNA genes. The C. acerifolia mitogenome has extremely unconserved structurally, with extensive sequence transfer between the chloroplast and mitochondrial organelles, sequence repeats, and RNA editing. The phylogenetic position of C. acerifolia was determined by constructing the species mitogenome with 24 angiosperms. Further, our C. acerifolia mitogenome characteristics investigation included GC contents, codon usage, repeats and synteny analysis. Overall, our results are expected to provide fundamental information for C. acerifolia mitogenome evolution and confirm the validity of mitochondrial analysis in determining the phylogenetic positioning of Clematis plants.

A comparative study of the chloroplast genomes of five Lepidium species with high medicinal value

Plantgenomics is a rapidly developing field in medicinal plant research. This study analysed the relevant information of chloroplasts genome sequences of five medicinal plants from the genus Lepidium . We sequenced the complete chloroplast (cp) genomes of Lepidium apetalum Willd. and Lepidium perfoliatum Linnaeus., and assessed their genetic profiles against the reported profiles of Lepidium sativum Linnaeus., Lepidium meyenii Walp., and Lepidium virginicum Linn. We found that L. apetalum and L. perfoliatum possessed 130 distinct genes that included 85 protein-coding, 37 transfer RNA (tRNA), and eight ribosomal RNA (rRNA) genes. Our repeat analyses revealed that L. apetalum harboured 20 direct repeats, 16 palindrome repeats, 30 tandem repeats, and 87 simple sequence repeats, whereas, L. perfoliatum had 15 direct repeats, 20 palindrome repeats, four reverse repeats, 21 tandem repeats, and 98 simple sequence repeats. Using syntenic analysis, we also revealed a high degree of sequence similarity within the coding regions of Lepidium medicinal plant cp genomes, and a high degree of divergence among the intergenic spacers. Pairwise alignment and single-nucleotide polymorphism (SNP) examinations further revealed certain Lepidium -specific gene fragments. Codon usage analysis showed that codon 14 was the most frequently used codon in the Lepidium coding sequences. Further, correlation investigations suggest that L. apetalum and L. perfoliatum originate from similar genetic backgrounds. Analysis of codon usage bias of Lepidium cp genome was strongly influenced by mutation and natural selection. We showed that L. apetalum and L. perfoliatum will likely enhance breeding, species recognition, phylogenetic evolution, and cp genetic engineering of the Lepidium medicinal plants.

Biased gene expression reveals the contribution of subgenome to altitude adaptation in allopolyploid Isoetes sinensis

Allopolyploids are believed to inherit the genetic characteristics of its progenitors and exhibit stronger adaptability and vigor. The allotetraploid Isoetes sinensis was formed by the natural hybridization and polyploidization of two diploid progenitors, Isoetes taiwanensis and Isoetes yunguiensis, and was believed to have the potential to adapt to plateau environments. To explore the expression pattern of homoeologous genes and their contributions to altitude adaptation, we transplanted natural allotetraploid I. sinensis (TnTnYnYn) along the altitude gradient for a long-term, and harvested them in summer and winter, respectively. One year after transplanting, it still lived well, even in the extreme environment of the Qinghai-Tibet Plateau. Then, we performed high-throughput RNA sequencing to measure their gene expression level. A total of 7801 homoeologous genes were expressed, among which 5786 were identified as shared expression in different altitudes and seasons. We further found that altitude variations could change the subgenome bias trend of I. sinensis, but season could not. Moreover, the functions of uniquely expressed genes indicated that temperature might be an important restrictive factor during the adaptation process. Through the analysis of DEGs and uniquely expressed genes, we found that Y subgenome provided more contributions to high altitude adaptation than T subgenome. These adaptive traits to high altitude may be inherited from its plateau progenitor I. yunguiensis. Through weighted gene co-expression network analysis, pentatricopeptide repeats gene family and glycerophospholipid metabolism pathway were considered to play important roles in high-altitude adaptation. Totally, this study will enrich our understanding of allopolyploid in environmental adaptation.

Characterization of the Complete Chloroplast Genome of Four Species in Callerya

BACKGROUND

Callerya reticulata (Bentham) Schot, Callerya dielsiana (Harms) P.K. Loc ex Z. Wei & Pedley, Callerya nitida var. hirsutissima (Z. Wei) X.Y. Zhu, and Callerya nitida (Bentham) R. Geesink, which belongs to the Leguminosae family, are important medicinal plants in China. The genus Callerya includes 26 species, 18 species are distributed in China, and the vine stems of some species are used as traditional medicinal herbs because they have important pharmacological activity. Due to the high similarity of appearance, it is difficult to identify them in the market by appearance alone. Therefore, circulating of Callerya-related materia medica on the market is confusing, sometimes even leading to drug safety problems. It is urgent to develop molecular methods for their identification.

OBJECTIVE

To sequence and analyze the complete chloroplast (cp) genomes of C. reticulata, C. dielsiana, C. nitida var. hirsutissima, and C. nitida and to analyze their cp genome differences as a basis for seeking easier DNA barcoding for their identification.

METHOD

After using Illumina high-throughput sequencing and nanopore sequencing to obtain the genome data, some bioinformatics software was used to assembly and analyze the molecular structure of cp genomes.

RESULTS

The complete cp genomes of the four species were circular molecules, which ranged from 130 435 to 132 546 bp, and GC contents ranged from 33.89% to 34.89%. Each of them includes a large single-copy region, a small single-copy region, and without large inverted repeat regions.

CONCLUSIONS

These results suggested that highly variable regions of the four cp genomes would provide useful plastid markers, which could be used as a potential genomic resource to resolve phylogenetic questions and provide a reference for mining specific DNA barcodes of these species.

HIGHLIGHTS

Our study provided highly effective molecular markers for subsequent phylogenetic analysis, species identification, and biogeographic analysis of Callerya.

Genome-wide identification and expression analysis of peach multiple organellar RNA editing factors reveals the roles of RNA editing in plant immunity

BACKGROUND

Multiple organellar RNA editing factor (MORF) genes play key roles in chloroplast developmental processes by mediating RNA editing of Cytosine-to-Uracil conversion. However, the function of MORF genes in peach (Prunus persica), a perennial horticultural crop species of Rosaceae, is still not well known, particularly the resistance to biotic and abiotic stresses that threaten peach yield seriously.

RESULTS

In this study, to reveal the regulatory roles of RNA editing in plant immunity, we implemented genome-wide analysis of peach MORF (PpMORF) genes in response to biotic and abiotic stresses. The chromosomal and subcellular location analysis showed that the identified seven PpMORF genes distributed on three peach chromosomes were mainly localized in the mitochondria and chloroplast. All the PpMORF genes were classified into six groups and one pair of PpMORF genes was tandemly duplicated. Based on the meta-analysis of two types of public RNA-seq data under different treatments (biotic and abiotic stresses), we observed down-regulated expression of PpMORF genes and reduced chloroplast RNA editing, especially the different response of PpMORF2 and PpMORF9 to pathogens infection between resistant and susceptible peach varieties, indicating the roles of MORF genes in stress response by modulating the RNA editing extent in plant immunity. Three upstream transcription factors (MYB3R-1, ZAT10, HSFB3) were identified under both stresses, they may regulate resistance adaption by modulating the PpMORF gene expression.

CONCLUSION

These results provided the foundation for further analyses of the functions of MORF genes, in particular the roles of RNA editing in plant immunity. In addition, our findings will be conducive to clarifying the resistance mechanisms in peaches and open up avenues for breeding new cultivars with high resistance.

RNA Editing in Chloroplast: Advancements and Opportunities

Many eukaryotic and prokaryotic organisms employ RNA editing (insertion, deletion, or conversion) as a post-transcriptional modification mechanism. RNA editing events are common in these organelles of plants and have gained particular attention due to their role in the development and growth of plants, as well as their ability to cope with abiotic stress. Owing to rapid developments in sequencing technologies and data analysis methods, such editing sites are being accurately predicted, and many factors that influence RNA editing are being discovered. The mechanism and role of the pentatricopeptide repeat protein family of proteins in RNA editing are being uncovered with the growing realization of accessory proteins that might help these proteins. This review will discuss the role and type of RNA editing events in plants with an emphasis on chloroplast RNA editing, involved factors, gaps in knowledge, and future outlooks.

Exploring the RNA Editing Events and Their Potential Regulatory Roles in Tea Plant (Camellia sinensis L.)

RNA editing is a post-transcriptional modification process that alters the RNA sequence relative to the genomic blueprint. In plant organelles (namely, mitochondria and chloroplasts), the most common type is C-to-U, and the absence of C-to-U RNA editing results in abnormal plant development, such as etiolation and albino leaves, aborted embryonic development and retarded seedling growth. Here, through PREP, RES-Scanner, PCR and RT-PCR analyses, 38 and 139 RNA editing sites were identified from the chloroplast and mitochondrial genomes of Camellia sinensis, respectively. Analysis of the base preference around the RNA editing sites showed that in the -1 position of the edited C had more frequent occurrences of T whereas rare occurrences of G. Three conserved motifs were identified at 25 bases upstream of the RNA editing site. Structural analyses indicated that the RNA secondary structure of 32 genes, protein secondary structure of 37 genes and the three-dimensional structure of 5 proteins were altered due to RNA editing. The editing level analysis of matK and ndhD in six tea cultivars indicated that matK-701 might be involved in the color change of tea leaves. Furthermore, 218 PLS-CsPPR proteins were predicted to interact with the identified RNA editing sites. In conclusion, this study provides comprehensive insight into RNA editing events, which will facilitate further study of the RNA editing phenomenon of the tea plant.

WAL3 encoding a PLS-type PPR protein regulates chloroplast development in rice

Chloroplast development is a complex process that is critical for the growth and development of plants. Pentapeptide repeat (PPR) proteins contain large members but only few of them have been characterized in rice. In this study, we identified a new PLS-type protein, WAL3 (Whole Albino Leaf on Chromosome 3), playing important roles in plant chloroplast development. Knockout of WAL3 gene in Nipponbare variety caused abnormal chloroplast development and showed an albino lethal phenotype. Expression analysis showed that WAL3 gene was constitutively expressed with the highest expression in leaves. The WAL3 protein localized in chloroplasts and affected the splicing of multiple group II introns. Transcriptome sequencing showed that WAL3 involved in multiple metabolic pathways including the chlorophyll synthesis and photosynthetic related metabolic pathways. The decreased abundance of photosynthesis-related proteins in wal3 mutants indicated WAL3 influence photosynthesis. In summary, our study revealed that WAL3 is essential for chloroplast development in rice.

Pentatricopeptide Repeat Gene-Mediated Mitochondrial RNA Editing Impacts on Rice Drought Tolerance

Mitochondrial RNA editing plays crucial roles in the plant development and environmental adaptation. Pentatricopeptide repeat (PPR) genes, which are involved in the regulating mitochondrial RNA editing, are potential gene resources in the improvement of rice drought tolerance. In this study, we investigated genome-wide mitochondrial RNA editing in response to drought between upland and lowland rice. Responses of mitochondrial RNA editing to drought exhibit site-specific and genotype-specific patterns. We detected 22 and 57 ecotype-differentiated editing sites under well-watered and drought-treated conditions, respectively. Interestingly, the RNA editing efficiency was positively correlated with many agronomic traits, while it was negatively correlated with drought tolerance. We further selected two mitochondrial-localized PPR proteins, PPR035 and PPR406, to validate their functions in drought tolerance. PPR035 regulated RNA editing at rps4-926 and orfX-406, while PPR406 regulated RNA editing at orfX-355. The defectiveness in RNA editing at these sites had no apparent penalties in rice respiration and vegetative growth. Meanwhile, the knockout mutants of ppr035 and ppr406 show enhanced drought- and salt tolerance. PPR035 and PPR406 were under the balancing selection in upland rice and highly differentiated between upland and lowland rice ecotypes. The upland-dominant haplotypes of PPR035 and PPR406 shall contribute to the better drought tolerance in upland rice. They have great prospective in the improvement of rice drought tolerance.

Decoding RNA Editing Sites Through Transcriptome Analysis in Rice Under Alkaline Stress

Ribonucleic acid editing (RE) is a post-transcriptional process that altered the genetics of RNA which provide the extra level of gene expression through insertion, deletions, and substitutions. In animals, it converts nucleotide residues C-U. Similarly in plants, the role of RNA editing sites (RES) in rice under alkaline stress is not fully studied. Rice is a staple food for most of the world population. Alkaline stress cause reduction in yield. Here, we explored the effect of alkaline stress on RES in the whole mRNA from rice chloroplast and mitochondria. Ribonucleic acid editing sites in both genomes (3336 RESs) including chloroplast (345 RESs) and mitochondria (2991 RESs) with average RES efficiency ∼55% were predicted. Our findings showed that majority of editing events found in non-synonymous codon changes and change trend in amino acids was hydrophobic. Four types of RNA editing A-G (A-I), C-T (C-U), G-A, and T-C were identified in treated and untreated samples. Overall, RNA editing efficiency was increased in the treated samples. Analysis of Gene Ontology revealed that mapped genes were engaged in many biological functions and molecular processes. We also checked the expression of pentatricopeptide repeat (PPR), organelle zinc-finger (OZI), and multiple organellar RNA editing factors/RNA editing factor interacting proteins genes in control and treatment, results revealed upregulation of PPR and OZ1 genes in treated samples. This induction showed the role of these genes in RNA editing. The current findings report that RNA editing increased under alkaline stress which may contribute in adaptation for rice by changing amino acids in edited genes (88 genes). These findings will provide basis for identification of RES in other crops and also will be useful in alkaline tolerance development in rice.

A mitochondrial RNA processing protein mediates plant immunity to a broad spectrum of pathogens by modulating the mitochondrial oxidative burst

Mitochondrial function depends on the RNA processing of mitochondrial gene transcripts by nucleus-encoded proteins. This posttranscriptional processing involves the large group of nuclear-encoded pentatricopeptide repeat (PPR) proteins. Mitochondrial processes represent a crucial part in animal immunity, but whether mitochondria play similar roles in plants remains unclear. Here, we report the identification of RESISTANCE TO PHYTOPHTHORA PARASITICA 7 (AtRTP7), a P-type PPR protein, in Arabidopsis thaliana and its conserved function in immunity to diverse pathogens across distantly related plant species. RTP7 affects the levels of mitochondrial reactive oxygen species (mROS) by participating in RNA splicing of nad7, which encodes a critical subunit of the mitochondrial respiratory chain Complex I, the largest of the four major components of the mitochondrial oxidative phosphorylation system. The enhanced resistance of rtp7 plants to Phytophthora parasitica is dependent on an elevated mROS burst, but might be independent from the ROS burst associated with plasma membrane-localized NADPH oxidases. Our study reveals the immune function of RTP7 and the defective processing of Complex I subunits in rtp7 plants resulted in enhanced resistance to both biotrophic and necrotrophic pathogens without affecting overall plant development.

Genome-wide investigation and functional analysis of RNA editing sites in wheat

Wheat is an important cereal and half of the world population consumed it. Wheat faces environmental stresses and different techniques (CRISPR, gene silencing, GWAS, etc.) were used to enhance its production but RNA editing (RESs) is not fully explored in wheat. RNA editing has a special role in controlling environmental stresses. The genome-wide identification and functional characterization of RESs in different types of wheat genotypes was done. We employed six wheat genotypes by RNA-seq analyses to achieve RESs. The findings revealed that RNA editing events occurred on all chromosomes equally. RNA editing sites were distributed randomly and 10–12 types of RESs were detected in wheat genotypes. Higher number of RESs were detected in drought-tolerant genotypes. A-to-I RNA editing (2952, 2977, 1916, 2576, 3422, and 3459) sites were also identified in six wheat genotypes. Most of the genes were found to be engaged in molecular processes after a Gene Ontology analysis. PPR (pentatricopeptide repeat), OZ1 (organelle zinc-finger), and MORF/RIP gene expression levels in wheat were also examined. Normal growth conditions diverge gene expression of these three different gene families, implying that normal growth conditions for various genotypes can modify RNA editing events and have an impact on gene expression levels. While the expression of PPR genes was not change. We used Variant Effect Predictor (VEP) to annotate RNA editing sites, and Local White had the highest RESs in the CDS region of the protein. These findings will be useful for prediction of RESs in other crops and will be helpful in drought tolerance development in wheat.

Identification of the Complete Chloroplast Genome of Malus zhaojiaoensis Jiang and Its Comparison and Evolutionary Analysis with Other Malus Species

The genus Malus is rich in species and many of its plastid genomes have been released. However, limited resources and few markers are not conducive to the comparison of differences among species and resource identification and evaluation. In this study, the complete chloroplast genome of Malus zhaojiaoensis was studied by NGS sequencing, with a total length of 159998 bp. It consists of four regions, LSC (88,070 bp), IRB (26,359 bp), SSC (19,210 bp) and IRA (26,359 bp). M. zhaojiaoensis cp genome contained a total of 111 genes made up of three classes: 78 coding sequences, 29 tRNA genes, and four rRNA genes. In addition, a total of 91 SSRs and 43 INEs were found in the M. zhaojiaoensis cp genome, which was slightly different from M. baccata and M. hupehensis in number. The analysis of codon usage and RNA editing showed that high-frequency codons tended to end at A/U bases and RNA editing mainly occurred at the second codon. Comparative genome analysis suggested that the cp genomes of eight Malus species had higher overall similarity, but there were more variation hotspots (rps16_trnK-UUU, trnG-UCC_atpA, atpH_atpF, trnT-GGU_psbD, etc.) in the LSC region. By building evolutionary trees, it can be clearly observed that M. zhaojiaoensis formed a large group with eight species of Malus, but was relatively independent in differentiation. In conclusion, this study provides high-quality chloroplast genome resources of M. zhaojiaoensis and discusses the genetic variation characteristics of Malus genus. The findings of this study will provide a good reference for plastid genome assembly and interspecific comparison in the future.

Maize PPR278 Functions in Mitochondrial RNA Splicing and Editing

Pentatricopeptide repeat (PPR) proteins are a large protein family in higher plants and play important roles during seed development. Most reported PPR proteins function in mitochondria. However, some PPR proteins localize to more than one organelle; functional characterization of these proteins remains limited in maize (Zea mays L.). Here, we cloned and analyzed the function of a P-subfamily PPR protein, PPR278. Loss-function of PPR278 led to a lower germination rate and other defects at the seedling stage, as well as smaller kernels compared to the wild type. PPR278 was expressed in all investigated tissues. Furthermore, we determined that PPR278 is involved in the splicing of two mitochondrial transcripts (nad2 intron 4 and nad5 introns 1 and 4), as well as RNA editing of C-to-U sites in 10 mitochondrial transcripts. PPR278 localized to the nucleus, implying that it may function as a transcriptional regulator during seed development. Our data indicate that PPR278 is involved in maize seed development via intron splicing and RNA editing in mitochondria and has potential regulatory roles in the nucleus.

Genome-Wide Analysis of Multiple Organellar RNA Editing Factor (MORF) Family in Kiwifruit (Actinidia chinensis) Reveals Its Roles in Chloroplast RNA Editing and Pathogens Stress

Kiwifruit (Actinidia chinensis) is well known for its high vitamin C content and good taste. Various diseases, especially bacterial canker, are a serious threat to the yield of kiwifruit. Multiple organellar RNA editing factor (MORF) genes are pivotal factors in the RNA editosome that mediates Cytosine-to-Uracil RNA editing, and they are also indispensable for the regulation of chloroplast development, plant growth, and response to stresses. Although the kiwifruit genome has been released, little is known about MORF genes in kiwifruit at the genome-wide level, especially those involved in the response to pathogens stress. In this study, we identified ten MORF genes in the kiwifruit genome. The genomic structures and chromosomal locations analysis indicated that all the MORF genes consisted of three conserved motifs, and they were distributed widely across the seven linkage groups and one contig of the kiwifruit genome. Based on the structural features of MORF proteins and the topology of the phylogenetic tree, the kiwifruit MORF gene family members were classified into six groups (Groups A-F). A synteny analysis indicated that two pairs of MORF genes were tandemly duplicated and five pairs of MORF genes were segmentally duplicated. Moreover, based on analysis of RNA-seq data from five tissues of kiwifruit, we found that both expressions of MORF genes and chloroplast RNA editing exhibited tissue-specific patterns. MORF2 and MORF9 were highly expressed in leaf and shoot, and may be responsible for chloroplast RNA editing, especially the ndhB genes. We also observed different MORF expression and chloroplast RNA editing profiles between resistant and susceptible kiwifruits after pathogen infection, indicating the roles of MORF genes in stress response by modulating the editing extend of mRNA. These results provide a solid foundation for further analyses of the functions and molecular evolution of MORF genes, in particular, for clarifying the resistance mechanisms in kiwifruits and breeding new cultivars with high resistance.

The Development of Mitochondrial Gene Editing Tools and Their Possible Roles in Crop Improvement for Future Agriculture

We are living in the era of genome editing. Nowadays, targeted editing of the plant nuclear DNA is prevalent in basic biological research and crop improvement since its first establishment a decade ago. However, achieving the same accomplishment for the plant mitochondrial genome has long been deemed impossible. Recently, the pioneer studies on editing plant mitogenome have been done using the mitochondria-targeted transcription activator-like effector nucleases (mitoTALENs) in rice, rapeseed, and Arabidopsis. It is well documented that mitochondria play essential roles in plant development and stress tolerance, particularly, in cytoplasmic male sterility widely used in production of hybrids. The success of mitochondrial genome editing enables studying the fundamentals of mitochondrial genome. Furthermore, mitochondrial RNA editing (mostly by nuclear-encoded pentatricopeptide repeat (PPR) proteins) in a sequence-specific manner can simultaneously change the production of translatable mitochondrial mRNA. Moreover, direct editing of the nuclear-encoding mitochondria-targeted factors required for plant mitochondrial genome dynamics and recombination may facilitate genetic manipulation of plant mitochondria. Here, the present state of knowledge on editing the plant mitochondrial genome is reviewed.

Molecular Analysis Uncovers the Mechanism of Fertility Restoration in Temperature-Sensitive Polima Cytoplasmic Male-Sterile Brassica napus

Temperature-sensitive male sterility is a heritable agronomic trait affected by genotype-environment interactions. In rapeseed (Brassica napus), Polima (pol) temperature-sensitive cytoplasmic male sterility (TCMS) is commonly used for two-line breeding, as the fertility of pol TCMS lines can be partially restored at certain temperatures. However, little is known about the underlying molecular mechanism that controls fertility restoration. Therefore, we aimed to investigate the fertility conversion mechanism of the pol TCMS line at two different ambient temperatures (16 °C and 25 °C). Our results showed that the anthers developed and produced vigorous pollen at 16 °C but not at 25 °C. In addition, we identified a novel co-transcript of orf224-atp6 in the mitochondria that might lead to fertility conversion of the pol TCMS line. RNA-seq analysis showed that 1637 genes were significantly differentially expressed in the fertile flowers of 596-L when compared to the sterile flower of 1318 and 596-H. Detailed analysis revealed that differentially expressed genes were involved in temperature response, ROS accumulation, anther development, and mitochondrial function. Single-molecule long-read isoform sequencing combined with RNA sequencing revealed numerous genes produce alternative splicing transcripts at high temperatures. Here, we also found that alternative oxidase, type II NAD(P)H dehydrogenases, and transcription factor Hsfs might play a crucial role in male fertility under the low-temperature condition. RNA sequencing and bulked segregant analysis coupled with whole-genome sequencing identified the candidate genes involved in the post-transcriptional modification of orf224. Overall, our study described a putative mechanism of fertility restoration in a pol TCMS line controlled by ambient temperature that might help utilise TCMS in the two-line breeding of Brassica crops.

RNA Editing and Its Roles in Plant Organelles

RNA editing, a vital supplement to the central dogma, yields genetic information on RNA products that are different from their DNA templates. The conversion of C-to-U in mitochondria and plastids is the main kind of RNA editing in plants. Various factors have been demonstrated to be involved in RNA editing. In this minireview, we summarized the factors and mechanisms involved in RNA editing in plant organelles. Recently, the rapid development of deep sequencing has revealed many RNA editing events in plant organelles, and we further reviewed these events identified through deep sequencing data. Numerous studies have shown that RNA editing plays essential roles in diverse processes, such as the biogenesis of chloroplasts and mitochondria, seed development, and stress and hormone responses. Finally, we discussed the functions of RNA editing in plant organelles.

Quantitative Analysis of RNA Editing at Specific Sites in Plant Mitochondria or Chloroplasts Using DNA Sequencing

Cytidine-to-uridine (C-to-U) RNA editing is one of the most important post-transcriptional RNA processing in plant mitochondria and chloroplasts. Several techniques have been developed to detect the RNA editing efficiency in plant mitochondria and chloroplasts, such as poisoned primer extension (PPE) assays, high-resolution melting (HRM) analysis, and DNA sequencing. Here, we describe a method for the quantitative detection of RNA editing at specific sites by sequencing cDNA from plant leaves to further evaluate the effect of different treatments or plant mutants on the C to U RNA editing in mitochondria and chloroplasts.

Genome Characterization, Comparison and Phylogenetic Analysis of Complete Mitochondrial Genome of Evolvulus alsinoides Reveals Highly Rearranged Gene Order in Solanales

Mitogenome sequencing provides an understanding of the evolutionary mechanism of mitogenome formation, mechanisms driving plant gene order, genome structure, and migration sequences. Data on the mitochondrial genome for family Convolvulaceae members is lacking. E. alsinoides, also known as shankhpushpi, is an important medicinal plant under the family Convolvulaceae, widely used in the Ayurvedic system of medicine. We identified the mitogenome of E. alsinoides using the Illumina mate-pair sequencing platform, and annotated using bioinformatics approaches in the present study. The mitogenome of E. alsinoides was 344184 bp in length and comprised 46 unique coding genes, including 31 protein-coding genes (PCGs), 12 tRNA genes, and 3 rRNA genes. The secondary structure of tRNAs shows that all the tRNAs can be folded into canonical clover-leaf secondary structures, except three trnW, trnG, and trnC. Measurement of the skewness of the nucleotide composition showed that the AT and GC skew is positive, indicating higher A's and G's in the mitogenome of E. alsinoides. The Ka/Ks ratios of 11 protein-coding genes (atp1, ccmC, cob, cox1, rps19, rps12, nad3, nad9, atp9, rpl5, nad4L) were <1, indicating that these genes were under purifying selection. Synteny and gene order analysis were performed to identify homologous genes among the related species. Synteny blocks representing nine genes (nad9, nad2, ccmFc, nad1, nad4, nad5, matR, cox1, nad7) were observed in all the species of Solanales. Gene order comparison showed that a high level of gene rearrangement has occurred among all the species of Solanales. The mitogenome data obtained in the present study could be used as the Convolvulaceae family representative for future studies, as there is no complex taxonomic history associated with this plant.

The Chloroplast Epitranscriptome: Factors, Sites, Regulation, and Detection Methods

Modifications in nucleic acids are present in all three domains of life. More than 170 distinct chemical modifications have been reported in cellular RNAs to date. Collectively termed as epitranscriptome, these RNA modifications are often dynamic and involve distinct regulatory proteins that install, remove, and interpret these marks in a site-specific manner. Covalent nucleotide modifications-such as methylations at diverse positions in the bases, polyuridylation, and pseudouridylation and many others impact various events in the lifecycle of an RNA such as folding, localization, processing, stability, ribosome assembly, and translational processes and are thus crucial regulators of the RNA metabolism. In plants, the nuclear/cytoplasmic epitranscriptome plays important roles in a wide range of biological processes, such as organ development, viral infection, and physiological means. Notably, recent transcriptome-wide analyses have also revealed novel dynamic modifications not only in plant nuclear/cytoplasmic RNAs related to photosynthesis but especially in chloroplast mRNAs, suggesting important and hitherto undefined regulatory steps in plastid functions and gene expression. Here we report on the latest findings of known plastid RNA modifications and highlight their relevance for the post-transcriptional regulation of chloroplast gene expression and their role in controlling plant development, stress reactions, and acclimation processes.

OsMORF9 is necessary for chloroplast development and seedling survival in rice

Chloroplasts are closely associated with the growth and development of higher plants. Accumulating evidence has revealed that the multiple organellar RNA editing factors (MORF) family of proteins influences plastidic and mitochondrial development through post-transcriptional regulation. However, the role of MORFs in regulating the development of chloroplasts in rice is still unclear. The OsMORF9 gene belongs to a small family of 7 genes in rice and is highly expressed in young leaves. We used the CRISPR/Cas9 system to mutate OsMORF9. The resulting knockout lines osmorf9-1 and osmorf9-2 exhibited an albino seedling lethal phenotype. Besides, the expression of many plastid-encoded genes involved in photosynthesis, the biogenesis of plastidic ribosomes and the editing and splicing of specific plastidic RNA molecules were severely affected in these two OsMORF9 mutants. Furthermore, yeast two-hybrid analysis revealed that OsMORF9 could interact with OsSLA4 and DUA1 which are members of the pentatricopeptide repeat (PPR) family of proteins. Analysis of subcellular localization of OsMORF9 also suggested that it might function in chloroplasts. The findings from the present study demonstrated the critical role of OsMORF9 in the biogenesis of chloroplast ribosomes, chloroplast development and seedling survival. This therefore provides new insights on the function of MORF proteins in rice.

Plastid caseinolytic protease OsClpR1 regulates chloroplast development and chloroplast RNA editing in rice

BACKGROUND

Plant plastidic caseinolytic protease (Clp) is a central part of the plastid protease network and consists of multiple subunits. The molecular functions of many Clps in plants, especially in crops, are not well known.

RESULTS

In this study, we identified an albino lethal mutant al3 in rice, which produces albino leaves and dies at the seedling stage. Molecular cloning revealed that AL3 encodes a plastid caseinolytic protease, OsClpR1, homologous to Arabidopsis ClpR1 and is targeted to the chloroplast. Compared with the wild type, chloroplast structure in the al3 mutant was poorly developed. OsClpR1 was constitutively expressed in all rice tissues, especially in young leaves. The OsClpR1 mutation affected the transcript levels of chlorophyll biosynthesis and chloroplast development-related genes. The RNA editing efficiency of three chloroplast genes (rpl2, ndhB, ndhA) was remarkably reduced in al3. Using a yeast two-hybrid screen, we found that OsClpR1 interacted with OsClpP4, OsClpP5, OsClpP2, and OsClpS1.

CONCLUSIONS

Collectively, our results provide novel insights into the function of Clps in rice.

Characterization of the Complete Chloroplast Genome of Buddleja Lindleyana

BACKGROUND

Buddleja lindleyana Fort., which belongs to the Loganiaceae with a distribution throughout the tropics, is widely used as an ornamental plant in China. Buddleja contains several morphologically similar species, which need to be identified by molecular identification. But there is little molecular research on the genus Buddleja.

OBJECTIVE

Using molecular biology techniques to sequence and analyze the complete chloroplast (cp) genome of B. lindleyana.

METHODS

According to next-generation sequencing to sequence the genome data, a series of bioinformatics software were used to assembly and analysis the molecular structure of cp genome of B. lindleyana.

RESULTS

The complete cp genome of B. lindleyana is a circular 154,487-bp-long molecule with a GC content of 38.1%. It has a familiar quadripartite structure, including a large single-copy region (LSC; 85,489 bp), a small single-copy region (SSC; 17,898bp) and a pair of inverted repeats (IRs; 25,550 bp). A total of 133 genes were identified in the genome, including 86 protein-coding genes, 37 tRNA genes, 8 rRNA genes and 2 pseudogenes.

CONCLUSIONS

These results suggested that B. lindelyana cp genome could be used as a potential genomic resource to resolve the phylogenetic positions and relationships of Loganiaceae, and will offer valuable information for future research in the identification of Buddleja species and will conduce to genomic investigations of these species.

Identification of polycistronic transcriptional units and non-canonical introns in green algal chloroplasts based on long-read RNA sequencing data

BACKGROUND

Chloroplasts are important semi-autonomous organelles in plants and algae. Unlike higher plants, the chloroplast genomes of green algal linage have distinct features both in organization and expression. Despite the architecture of chloroplast genome having been extensively studied in higher plants and several model species of algae, little is known about the transcriptional features of green algal chloroplast-encoded genes.

RESULTS

Based on full-length cDNA (Iso-Seq) sequencing, we identified widely co-transcribed polycistronic transcriptional units (PTUs) in the green alga Caulerpa lentillifera. In addition to clusters of genes from the same pathway, we identified a series of PTUs of up to nine genes whose function in the plastid is not understood. The RNA data further allowed us to confirm widespread expression of fragmented genes and conserved open reading frames, which are both important features in green algal chloroplast genomes. In addition, a newly fragmented gene specific to C. lentillifera was discovered, which may represent a recent gene fragmentation event in the chloroplast genome. With the newly annotated exon-intron boundary information, gene structural annotation was greatly improved across the siphonous green algae lineages. Our data also revealed a type of non-canonical Group II introns, with a deviant secondary structure and intronic ORFs lacking known splicing or mobility domains. These widespread introns have conserved positions in their genes and are excised precisely despite lacking clear consensus intron boundaries.

CONCLUSION

Our study fills important knowledge gaps in chloroplast genome organization and transcription in green algae, and provides new insights into expression of polycistronic transcripts, freestanding ORFs and fragmented genes in algal chloroplast genomes. Moreover, we revealed an unusual type of Group II intron with distinct features and conserved positions in Bryopsidales. Our data represents interesting additions to knowledge of chloroplast intron structure and highlights clusters of uncharacterized genes that probably play important roles in plastids.

The FpPPR1 Gene Encodes a Pentatricopeptide Repeat Protein That Is Essential for Asexual Development, Sporulation, and Pathogenesis in Fusarium pseudograminearum

Fusarium crown rot (FCR) and Fusarium head blight (FHB) are caused by Fusarium pseudograminearum and are newly emerging diseases of wheat in China. In this study, we characterized FpPPR1, a gene that encodes a protein with 12 pentatricopeptide repeat (PPR) motifs. The radial growth rate of the ΔFpppr1 deletion mutant was significantly slower than the wild type strain WZ-8A on potato dextrose agar plates and exhibited significantly smaller colonies with sector mutations. The aerial mycelium of the mutant was almost absent in culture tubes. The ΔFpppr1 mutant was able to produce spores, but spores of abnormal size and altered conidium septum shape were produced with a significant reduction in sporulation compared to wild type. ΔFpppr1 failed to cause disease on wheat coleoptiles and barley leaves using mycelia plugs or spore suspensions. The mutant phenotypes were successfully restored to the wild type levels in complemented strains. FpPpr1-GFP signals in spores and mycelia predominantly overlapped with Mito-tracker signals, which substantiated the mitochondria targeting signal prediction of FpPpr1. RNAseq revealed significant transcriptional changes in the ΔFpppr1 mutant with 1,367 genes down-regulated and 1,333 genes up-regulated. NAD-binding proteins, thioredoxin, 2Fe-2S iron-sulfur cluster binding domain proteins, and cytochrome P450 genes were significantly down-regulated in ΔFpppr1, implying the dysfunction of mitochondria-mediated reductase redox stress in the mutant. The mating type idiomorphic alleles MAT1-1-1, MAT1-1-2, and MAT1-1-3 in F. pseudograminearum were also down-regulated after deletion of FpPPR1 and validated by real-time quantitative PCR. Additionally, 21 genes encoding putative heterokaryon incompatibility proteins were down-regulated. The yellow pigmentation of the mutant was correlated with reduced expression of PKS12 cluster genes. Taken together, our findings on FpPpr1 indicate that this PPR protein has multiple functions in fungal asexual development, regulation of heterokaryon formation, mating-type, and pathogenesis in F. pseudograminearum.

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.

The Role of Chloroplast Gene Expression in Plant Responses to Environmental Stress

Chloroplasts are plant organelles that carry out photosynthesis, produce various metabolites, and sense changes in the external environment. Given their endosymbiotic origin, chloroplasts have retained independent genomes and gene-expression machinery. Most genes from the prokaryotic ancestors of chloroplasts were transferred into the nucleus over the course of evolution. However, the importance of chloroplast gene expression in environmental stress responses have recently become more apparent. Here, we discuss the emerging roles of the distinct chloroplast gene expression processes in plant responses to environmental stresses. For example, the transcription and translation of psbA play an important role in high-light stress responses. A better understanding of the connection between chloroplast gene expression and environmental stress responses is crucial for breeding stress-tolerant crops better able to cope with the rapidly changing environment.

Plant editosome database: a curated database of RNA editosome in plants

RNA editing plays an important role in plant development and growth, enlisting a number of editing factors in the editing process and accordingly revealing the diversity of plant editosomes for RNA editing. However, there is no resource available thus far that integrates editosome data for a variety of plants. Here, we present Plant Editosome Database (PED; http://bigd.big.ac.cn/ped), a curated database of RNA editosome in plants that is dedicated to the curation, integration and standardization of plant editosome data. Unlike extant relevant databases, PED incorporates high-quality editosome data manually curated from related publications and organelle genome annotations. In the current version, PED integrates a complete collection of 98 RNA editing factors and 20 836 RNA editing events, covering 203 organelle genes and 1621 associated species. In addition, it contains functional effects of editing factors in regulating plant phenotypes and includes detailed experimental evidence. Together, PED serves as an important resource to help researchers investigate the RNA editing process across a wide range of plants and thus would be of broad utility for the global plant research community.

Dynamic response of RNA editing to temperature in grape by RNA deep sequencing

RNA editing is a post-transcriptional process of modifying genetic information on RNA molecules, which provides cells an additional level of gene expression regulation. Unlike mammals, in land plants, RNA editing converts C-to-U residues in organelles. However, its potential roles in response to different stressors (heat, salt, and so on) remains unclear. Grape is one of the most popular and economically important fruits in the world, and its production, like other crops, must deal with abiotic and biotic stresses, which cause reductions in yield and fruit quality. In our study, we tested the influence of the environmental factor temperature on RNA editing process in the whole mRNA from grape organelle. In total, we identified 122 and 627 RNA editing sites in chloroplast and mitochondria respectively with the average editing efficiency nearly ~ 60%. The analyses revealed that number of non-synonymous editing sites were higher than that of synonymous editing sites, and the amino acid substitution type tends to be hydrophobic. Additionally, the overall editing level decreased with the temperature rises, especially for several gene transcripts in chloroplast and mitochondria (matK, ndhB, etc.). We also found that the expression level of most PPR genes decreased with the temperature rises, which may contribute to the decline of RNA editing efficiency at high temperature. Our findings suggested that the RNA editing events were very sensitive to heat stress; the changes of amino acid in RNA editing genes may contribute to the stress adaption for grape.

Characterization, Expression, and Interaction Analyses of OsMORF Gene Family in Rice

The multiple organellar RNA editing factors (MORF) gene family plays a key role in organelle RNA editing in flowering plants. MORF genes expressions are also affected by abiotic stress. Although seven OsMORF genes have been identified in rice, few reports have been published on their expression patterns in different tissues and under abiotic stress, and OsMORF–OsMORF interactions. In this study, we analyzed the gene structure of OsMORF family genes. The MORF family members were divided into six subgroups in different plants based on phylogenetic analysis. Seven OsMORF genes were highly expressed in leaves. Six and seven OsMORF genes expressions were affected by cold and salt stresses, respectively. OsMORF–OsMORF interaction analysis indicated that OsMORF1, OsMORF8a, and OsMORF8b could each interact with themselves to form homomers. Moreover, five OsMORF proteins were shown to be able to interact with each other, such as OsMORF8a and OsMORF8b interacting with OsMORF1 and OsMORF2b, respectively, to form heteromers. These results provide information for further study of OsMORF gene function.

GhYGL1d, a pentatricopeptide repeat protein, is required for chloroplast development in cotton

BACKGROUND

The pentatricopeptide repeat (PPR) gene family, which contains multiple 35-amino acid repeats, constitutes one of the largest gene families in plants. PPR proteins function in organelles to target specific transcripts and are involved in plant development and growth. However, the function of PPR proteins in cotton is still unknown.

RESULTS

In this study, we characterized a PPR gene YELLOW-GREEN LEAF (GhYGL1d) that is required for cotton plastid development. The GhYGL1d gene has a DYW domain in C-terminal and is highly express in leaves, localized to the chloroplast fractions. GhYGL1d share high amino acid-sequence homology with AtECB2. In atecb2 mutant, overexpression of GhYGL1d rescued the seedling lethal phenotype and restored the editing of accD and ndhF transcripts. Silencing of GhYGL1d led to the reduction of chlorophyll and phenotypically yellow-green leaves in cotton. Compared with wild type, GhYGL1d-silenced cotton showed significant deformations of thylakoid structures. Furthermore, the transcription levels of plastid-encoded polymerase (PEP) and nuclear-encoded polymerase (NEP) dependent genes were decreased in GhYGL1d-silenced cotton.

CONCLUSIONS

Our data indicate that GhYGL1d not only contributes to the editing of accD and ndhF genes, but also affects the expression of NEP- and PEP-dependent genes to regulate the development of thylakoids, and therefore regulates leaf variegation in cotton.

Genome-Wide Analysis of Multiple Organellar RNA Editing Factor Family in Poplar Reveals Evolution and Roles in Drought Stress

Poplar (Populus) is one of the most important woody plants worldwide. Drought, a primary abiotic stress, seriously affects poplar growth and development. Multiple organellar RNA editing factor (MORF) genes-pivotal factors in the RNA editosome in Arabidopsis thaliana-are indispensable for the regulation of various physiological processes, including organelle C-to-U RNA editing and plasmid development, as well as in the response to stresses. Although the poplar genome sequence has been released, little is known about MORF genes in poplar, especially those involved in the response to drought stress at the genome-wide level. In this study, we identified nine MORF genes in the Populus genome. Based on the structural features of MORF proteins and the topology of the phylogenetic tree, the P. trichocarpa (Ptr) MORF family members were classified into six groups (Groups I⁻VI). A microsynteny analysis indicated that two (22.2%) PtrMORF genes were tandemly duplicated and seven genes (77.8%) were segmentally duplicated. Based on the d_N/d_S ratios, purifying selection likely played a major role in the evolution of this family and contributed to functional divergence among PtrMORF genes. Moreover, analysis of qRT-PCR data revealed that PtrMORFs exhibited tissue- and treatment-specific expression patterns. PtrMORF genes in all group were involved in the stress response. These results provide a solid foundation for further analyses of the functions and molecular evolution of MORF genes in poplar, and, in particular, for improving the drought resistance of poplar by genetics manipulation.

High Nighttime Temperature Induces Antioxidant Molecule Perturbations in Heat-Sensitive and Heat-Tolerant Coisogenic Rice ( Oryza sativa) Strains

Global warming-associated increases in temperature, particularly at nighttime, are detrimental to rice yield and quality. Metabolomic profiling was used to examine and compare the short-term extreme high nighttime temperature-induced molecular perturbations in rice ( Oryza sativa) coisogenic strains with contrasting heat-tolerances at the first stage of seed ripening. Compared to the heat-sensitive strain, antioxidant molecules were higher in abundance in the heat-tolerant strain, whereas the abundances of molecules involved in photosynthesis, nucleotide catabolism, and the S-adenosylmethionine (SAM) cycle varied only slightly. Thus, we proposed that the high abundance of antioxidant molecules in the heat-tolerant strain alleviated cellular oxidative stress, which protected photosynthesis, nucleotide catabolism, and the SAM cycle, leading to good grain filling.