RNA editing in plants and its evolution

Fujian cytoplasmic male sterility and the fertility restorer gene OsRf19 provide a promising breeding system for hybrid rice

Although hybrid rice has been widely utilized for nearly half a century, tremendously improving rice productivity worldwide, the breeding of hybrids has been difficult because of genetic complications in male sterility and fertility-restoring systems currently available in rice. Here, we characterized Fujian Abortive cytoplasmic male sterility (CMS-FA) rice, which has shown stable male sterility controlled by the mitochondrial gene FA182; a single nuclear gene, OsRf19, completely restores fertility. This single-gene inheritance has greatly eased the breeding process. By converting CMS-WA hybrids with the CMS-FA system, we developed six hybrids that showed equivalent or better performance relative to their CMS-WA counterparts. CMS-FA/OsRf19 provides a promising system for future hybrid rice breeding.

Cytoplasmic male sterility (CMS) determined by mitochondrial genes and restorer of fertility (Rf) controlled by nuclear-encoded genes provide the breeding systems of many hybrid crops for the utilization of heterosis. Although several CMS/Rf systems have been widely exploited in rice, hybrid breeding using these systems has encountered difficulties due to either fertility instability or complications of two-locus inheritance or both. In this work, we characterized a type of CMS, Fujian Abortive cytoplasmic male sterility (CMS-FA), with stable sporophytic male sterility and a nuclear restorer gene that completely restores hybrid fertility. CMS is caused by the chimeric open reading frame FA182 that specifically occurs in the mitochondrial genome of CMS-FA rice. The restorer gene OsRf19 encodes a pentatricopeptide repeat (PPR) protein targeted to mitochondria, where it mediates the cleavage of FA182 transcripts, thus restoring male fertility. Comparative sequence analysis revealed that OsRf19 originated through a recent duplication in wild rice relatives, sharing a common ancestor with OsRf1a/OsRf5, a fertility restorer gene for Boro II and Hong-Lian CMS. We developed six restorer lines by introgressing OsRf19 into parental lines of elite CMS-WA hybrids; hybrids produced from these lines showed equivalent or better agronomic performance relative to their counterparts based on the CMS-WA system. These results demonstrate that CMS-FA/OsRf19 provides a highly promising system for future hybrid rice breeding.

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.

Genome-Wide Identification and Characterization of RNA/DNA Differences Associated with Fusarium graminearum Infection in Wheat

RNA/DNA difference (RDD) is a post-transcriptional modification playing a crucial role in regulating diverse biological processes in eukaryotes. Although it has been extensively studied in plant chloroplast and mitochondria genomes, RDDs in plant nuclear genomes are not well studied at present. Here, we investigated the RDDs associated with fusarium head blight (FHB) through a novel method by comparing the RNA-seq data between Fusarium-infected and control samples of four wheat genotypes. A total of 187 high-confidence unique RDDs in 36 genes were identified, representing the first landscape of the FHB-responsive RDD in wheat. The majority (26) of these 36 RDD genes were correlated either positively or negatively with FHB levels. Effects of these RDDs on RNA and protein sequences have been identified, their editing frequency and the expression level of the corresponding genes provided, and the prediction of the effect on the minimum folding free energy of mRNA, miRNA binding, and colocation of RDDs with conserved domains presented. RDDs were predicted to induce modifications in the mRNA and protein structures of the corresponding genes. In two genes, TraesCS1B02G294300 and TraesCS3A02G263900, editing was predicted to enhance their affinity with tae-miR9661-5p and tae-miR9664-3p, respectively. To our knowledge, this study is the first report of the association between RDD and FHB in wheat; this will contribute to a better understanding of the molecular basis underlying FHB resistance, and potentially lead to novel strategies to improve wheat FHB resistance through epigenetic methods.

Characterization of the complete chloroplast genome sequences of six Dalbergia species and its comparative analysis in the subfamily of Papilionoideae (Fabaceae)

Dalbergia spp. are numerous and widely distributed in pantropical areas in Asia, Africa and America, and most of the species have important economic and ecological value as precious timber. In this study, we determined and characterized six complete chloroplast genomes of Dalbergia species (Dalbergia obtusifolia, D. hupeana, D. mimosoides, D. sissoo, D. hancei, D. balansae), which displayed the typical quadripartite structure of angiosperms. The sizes of the genomes ranged from 155,698 bp (D. hancei) to 156,419 bp (D. obtusifolia). The complete chloroplast genomes of Dalbergia include 37 tRNA genes, eight rRNA genes and 84 protein-coding genes. We analysed the sequence diversity of Dalberigia chloroplast genomes coupled with previous reports. The results showed 12 noncoding regions (rps16-accD, trnR-UCU-trnG-UCC, ndhE-ndhG, trnG-UCC-psbZ, rps8-rpl14, trnP-UGG-psaJ, ndhH-rps15, trnQ-UUG-rps16, trnS-GCU-psbI, rps12-clpP, psbA-trnK-UUU, trnK-UUU-intron), and four coding regions (rps16, ycf1, rps15 and ndhF) showed many nucleotide variations that could be used as potential molecular markers. Based on a site-specific model, we analysed the selective pressure of chloroplast genes in Dalbergia species. Twenty-two genes with positively selected sites were detected, involving the photosynthetic system (ndhC, adhD, ndhF, petB, psaA, psaB, psbB, psbC, psbK and rbcL), self-replication category of genes (rpoA, rpoC2, rps3, rps12 and rps18) and others (accD, ccsA, cemA, clpP, matK, ycf1 and ycf2). Additionally, we identified potential RNA editing sites that were relatively conserved in the genus Dalbergia. Furthermore, the comparative analysis of cp genomes of Dalbergieae species indicated that the boundary of IRs/SSC was highly variable, which resulted in the size variation of cp genomes. Finally, phylogenetic analysis showed an inferred phylogenetic tree of Papilionoideae species with high bootstrap support and suggested that Amorpheae was the sister of the clade Dalbergieae. Moreover, three genera of the Pterocarpus clade showed a nested evolutionary relationship. These complete cp genomes provided valuable information for understanding the genetic variation and phylogenetic relationship of Dalbergia species with their relatives.

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

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

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.

Organization, Phylogenetic Marker Exploitation, and Gene Evolution in the Plastome of Thalictrum (Ranunculaceae)

Thalictrum is a phylogenetically and economically important genus in the family Ranunculaceae, but is also regarded as one of the most challengingly difficult in plants for resolving the taxonomical and phylogenetical relationships of constituent taxa within this genus. Here, we sequenced the complete plastid genomes of two Thalictrum species using Illumina sequencing technology via de novo assembly. The two Thalictrum plastomes exhibited circular and typical quadripartite structure that was rather conserved in overall structure and the synteny of gene order. By updating the previously reported plastome annotation of other nine Thalictrum species, we found that the expansion or contraction of the inverted repeat region affect the boundary of the single-copy regions in Thalictrum plastome. We identified eight highly variable noncoding regions-infA-rps8, ccsA-ndhD, trnSUGA-psbZ, trnHGUG-psbA, rpl16-rps3, ndhG-ndhI, ndhD-psaC, and ndhJ-ndhK-that can be further used for molecular identification, phylogenetic, and phylogeographic in different species. Selective pressure and codon usage bias of all the plastid coding genes were also analyzed for the 11 species. Phylogenetic relationships showed Thalictrum is monophyly and divided into two major clades based on 11 Thalictrum plastomes. The availability of these plastomes offers valuable genetic information for accurate identification of species and taxonomy, phylogenetic resolution, and evolutionary studies of Thalictrum, and should assist with exploration and utilization of Thalictrum plants.

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.

Complete Chloroplast Genome of Cnidium monnieri (Apiaceae) and Comparisons with Other Tribe Selineae Species

Cnidium monnieri is an economically important traditional Chinese medicinal plant. In this study, the complete chloroplast (cp) genome of C. monnieri was determined using the Illumina paired-end sequencing, the GetOrganelle de novo assembly strategy, as well as the GeSeq annotation method. Our results showed that the cp genome was 147,371 bp in length with 37.4% GC content and included a large single-copy region (94,361 bp) and a small single-copy region (17,552 bp) separated by a pair of inverted repeat regions (17,729 bp). A total of 129 genes were contained in the cp genome, including 85 protein-coding genes, 36 tRNA genes, and eight rRNA genes. We also investigated codon usage, RNA editing, repeat sequences, simple sequence repeats (SSRs), IR boundaries, and pairwise Ka/Ks ratios. Four hypervariable regions (trnD-trnY-trnE-trnT, ycf2, ndhF-rpl32-trnL, and ycf1) were identified as candidate molecular markers for species authentication. The phylogenetic analyses supported non-monophyly of Cnidium and C. monnieri located in tribe Selineae based on the cp genome sequences and internal transcribed spacer (ITS) sequences. The incongruence of the phylogenetic position of C. monnieri between ITS and cpDNA phylogenies suggested that C. monnieri might have experienced complex evolutions with hybrid and incomplete lineage sorting. All in all, the results presented herein will provide plentiful chloroplast genomic resources for studies of the taxonomy, phylogeny, and species authentication of C. monnieri. Our study is also conducive to elucidating the phylogenetic relationships and taxonomic position of Cnidium.

DEK48 Is Required for RNA Editing at Multiple Mitochondrial Sites and Seed Development in Maize

In flowering plants, C-to-U RNA editing can be critical to normal functions of mitochondrion-encoded proteins. Mitochondrial C-to-U RNA editing is facilitated by many factors from diverse protein families, of which the pentatricopeptide repeat (PPR) proteins play an important role. Owing to their large number and frequent embryo lethality in mutants, functions of many PPRs remain unknown. In this study, we characterized a mitochondrion-localized DYW-type PPR protein, DEK48, functioning in the C-to-U RNA editing at multiple mitochondrial transcripts in maize. Null mutation of Dek48 severely arrests embryo and endosperm development, causing a defective kernel (dek) phenotype, named dek48. DEK48 loss of function abolishes the C-to-U editing at nad3-185, -215, and nad4-376, -977 sites and decreases the editing at 11 other sites, resulting in the alteration of the corresponding amino acids. Consequently, the absence of editing caused reduced assembly and activity of complex I in dek48. Interestingly, we identified a point mutation in dek48-3 causing a deletion of the Tryptophan (W) residue in the DYW motif that abolishes the editing function. In sum, this study reveals the function of DEK48 in the C-to-U editing in mitochondrial transcripts and seed development in maize, and it demonstrates a critical role of the W residue in the DYW triplet motif of DEK48 for the C-to-U editing function in vivo.

Dek504 Encodes a Mitochondrion-Targeted E+-Type Pentatricopeptide Repeat Protein Essential for RNA Editing and Seed Development in Maize

In flowering plants, RNA editing is a post-transcriptional process that selectively deaminates cytidines (C) to uridines (U) in organellar transcripts. Pentatricopeptide repeat (PPR) proteins have been identified as site-specific recognition factors for RNA editing. Here, we report the map-based cloning and molecular characterization of the defective kernel mutant dek504 in maize. Loss of Dek504 function leads to delayed embryogenesis and endosperm development, which produce small and collapsed kernels. Dek504 encodes an E+-type PPR protein targeted to the mitochondria, which is required for RNA editing of mitochondrial NADH dehydrogenase 3 at the nad3-317 and nad3-44 sites. Biochemical analysis of mitochondrial protein complexes revealed a significant reduction in the mitochondrial NADH dehydrogenase complex I activity, indicating that the alteration of the amino acid sequence at nad3-44 and nad3-317 through RNA editing is essential for NAD3 function. Moreover, the amino acids are highly conserved in monocots and eudicots, whereas the events of C-to-U editing are not conserved in flowering plants. Thus, our results indicate that Dek504 is essential for RNA editing of nad3, which is critical for NAD3 function, mitochondrial complex I stability, and seed development in maize.

Physcomitrium patens pentatricopeptide repeat protein PpPPR_32 is involved in the accumulation of psaC mRNA encoding the iron sulfur protein of photosystem I

Pentatricopeptide repeat (PPR) proteins are involved in RNA metabolism and also play a role in posttranscriptional regulation during plant organellar gene expression. Although a hundred of PPR proteins exist in the moss Physcomitrium patens, their functions are not fully understood. Here, we report the function of P-class PPR protein PpPPR_32 in P. patens. A transient expression assay using green fluorescent protein demonstrated that the N-terminal region of PpPPR_32 functions as a chloroplast-targeting transit peptide, indicating that PpPPR_32 is localized in chloroplasts. PpPPR_32 knockout (KO) mutants grew autotrophically but with reduced protonema growth and the poor formation of photosystem I (PSI) complexes. Quantitative real-time reverse transcription-polymerase chain reaction and RNA gel blot hybridization analyses revealed a significant reduction in the transcript level of the psaC gene encoding the iron sulfur protein of PSI but no alteration to the transcript levels of other PSI genes. This suggests that PpPPR_32 is specifically involved in the expression level of the psaC gene. Our results indicate that PpPPR_32 is essential for the accumulation of psaC transcript and PSI complexes.

The plastome of Melocactus glaucescens Buining & Brederoo reveals unique evolutionary features and loss of essential tRNA genes

The plastome of Melocactus glaucescens shows unique rearrangements, IR expansion, and unprecedented gene losses in Cactaceae. Our data indicate tRNA import from the cytosol to the plastids in this species. Cactaceae represents one of the richest families in keystone species of arid and semiarid biomes. This family shows various specific features comprehending morphology, anatomy, and metabolism, which allow them to grow under unfavorable environmental conditions. The subfamily Cactoideae contains the most divergence of species, which are highly variable in growth habit and morphology. This subfamily includes the endangered species Melocactus glaucescens (tribe Cereeae), which is a cactus endemic to the biome Caatinga in Brazil. Aiming to analyze the plastid evolution and develop molecular markers, we sequenced and analyzed in detail the plastome of M. glaucescens. Our analyses revealed that the M. glaucescens plastome is the most divergent among the species of the family Cactaceae sequenced so far. We characterized here unique rearrangements, expanded IRs containing an unusual set of genes, and several gene losses. Some genes related to the ndh complex were lost during the plastome evolution, while others have lost their functionality. Additionally, the loss of three tRNA genes (trnA-UGC, trnV-UAC, and trnV-GAC) suggests tRNA import from the cytosol to the plastids in M. glaucescens. Moreover, we identified high gene divergence, several putative positive signatures, and possible unique RNA-editing sites. Furthermore, we mapped 169 SSRs in the plastome of M. glaucescens, which are helpful to access the genetic diversity of natural populations and conservation strategies. Finally, our data provide new insights into the evolution of plastids in Cactaceae, which is an outstanding lineage adapted to extreme environmental conditions and a notorious example of the atypical evolution of plastomes.

A Structural Perspective on the RNA Editing of Plant Respiratory Complexes

The last steps of respiration, a core energy-harvesting process, are carried out by a chain of multi-subunit complexes in the inner mitochondrial membrane. Several essential subunits of the respiratory complexes are RNA-edited in plants, frequently leading to changes in the encoded amino acids. While the impact of RNA editing is clear at the sequence and phenotypic levels, the underlying biochemical explanations for these effects have remained obscure. Here, we used the structures of plant respiratory complex I, complex III₂ and complex IV to analyze the impact of the amino acid changes of RNA editing in terms of their location and biochemical features. Through specific examples, we demonstrate how the structural information can explain the phenotypes of RNA-editing mutants. This work shows how the structural perspective can bridge the gap between sequence and phenotype and provides a framework for the continued analysis of RNA-editing mutants in plant mitochondria and, by extension, in chloroplasts.

The Arabidopsis PeptideAtlas: Harnessing worldwide proteomics data to create a comprehensive community proteomics resource

We developed a resource, the Arabidopsis PeptideAtlas (www.peptideatlas.org/builds/arabidopsis/), to solve central questions about the Arabidopsis thaliana proteome, such as the significance of protein splice forms and post-translational modifications (PTMs), or simply to obtain reliable information about specific proteins. PeptideAtlas is based on published mass spectrometry (MS) data collected through ProteomeXchange and reanalyzed through a uniform processing and metadata annotation pipeline. All matched MS-derived peptide data are linked to spectral, technical, and biological metadata. Nearly 40 million out of ∼143 million MS/MS (tandem MS) spectra were matched to the reference genome Araport11, identifying ∼0.5 million unique peptides and 17,858 uniquely identified proteins (only isoform per gene) at the highest confidence level (false discovery rate 0.0004; 2 non-nested peptides ≥9 amino acid each), assigned canonical proteins, and 3,543 lower-confidence proteins. Physicochemical protein properties were evaluated for targeted identification of unobserved proteins. Additional proteins and isoforms currently not in Araport11 were identified that were generated from pseudogenes, alternative start, stops, and/or splice variants, and small Open Reading Frames; these features should be considered when updating the Arabidopsis genome. Phosphorylation can be inspected through a sophisticated PTM viewer. PeptideAtlas is integrated with community resources including TAIR, tracks in JBrowse, PPDB, and UniProtKB. Subsequent PeptideAtlas builds will incorporate millions more MS/MS data.

DYW domain structures imply an unusual regulation principle in plant organellar RNA editing catalysis

RNA editosomes selectively deaminate cytidines to uridines in plant organellar transcripts–mostly to restore protein functionality and consequently facilitate mitochondrial and chloroplast function. The RNA editosomal pentatricopeptide repeat proteins serve target RNA recognition, whereas the intensively studied DYW domain elicits catalysis. Here we present structures and functional data of a DYW domain in an inactive ground state and activated. DYW domains harbour a cytidine deaminase fold and a C-terminal DYW motif, with catalytic and structural zinc atoms, respectively. A conserved gating domain within the deaminase fold regulates the active site sterically and mechanistically in a process that we termed gated zinc shutter. Based on the structures, an autoinhibited ground state and its activation are cross-validated by RNA editing assays and differential scanning fluorimetry. We anticipate that, in vivo, the framework of an active plant RNA editosome triggers the release of DYW autoinhibition to ensure a controlled and coordinated cytidine deamination playing a key role in mitochondrial and chloroplast homeostasis.

Quantification of Mitochondrial RNA Editing Efficiency Using Sanger Sequencing Data

C-to-U RNA editing in mitochondria and plastids is widespread in almost all terrestrial plants, where it mainly changes codons to encode conserved amino acids in organelle mRNAs. In flowering plants, the number of RNA editing sites reaches 400-600 in mitochondria and about 40 in plastids, respectively. To date, more than 100 factors involved in RNA editing have been identified. Since target cytidines of each factor are often distributed across multiple transcripts, comprehensive monitoring of all RNA editing sites is necessary for their characterization. Comparing the signals of C and T in the Sanger sequencing chromatogram of RT-PCR products is the most frequently employed method for quantification of RNA editing efficiency, although several methods based on next-generation sequencing have been developed. I here describe a quick and easy method for quantification of RNA editing efficiency at several hundred sites using the Sanger sequencing chromatogram data.

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.

Precise plant genome editing using base editors and prime editors

The development of CRISPR-Cas systems has sparked a genome editing revolution in plant genetics and breeding. These sequence-specific RNA-guided nucleases can induce DNA double-stranded breaks, resulting in mutations by imprecise non-homologous end joining (NHEJ) repair or precise DNA sequence replacement by homology-directed repair (HDR). However, HDR is highly inefficient in many plant species, which has greatly limited precise genome editing in plants. To fill the vital gap in precision editing, base editing and prime editing technologies have recently been developed and demonstrated in numerous plant species. These technologies, which are mainly based on Cas9 nickases, can introduce precise changes into the target genome at a single-base resolution. This Review provides a timely overview of the current status of base editors and prime editors in plants, covering both technological developments and biological applications.

A photosynthesis operon in the chloroplast genome drives speciation in evening primroses

Genetic incompatibility between the cytoplasm and the nucleus is thought to be a major factor in species formation, but mechanistic understanding of this process is poor. In evening primroses (Oenothera spp.), a model plant for organelle genetics and population biology, hybrid offspring regularly display chloroplast-nuclear incompatibility. This usually manifests in bleached plants, more rarely in hybrid sterility or embryonic lethality. Hence, most of these incompatibilities affect photosynthetic capability, a trait that is under selection in changing environments. Here we show that light-dependent misregulation of the plastid psbB operon, which encodes core subunits of photosystem II and the cytochrome b6f complex, can lead to hybrid incompatibility, and this ultimately drives speciation. This misregulation causes an impaired light acclimation response in incompatible plants. Moreover, as a result of their different chloroplast genotypes, the parental lines differ in photosynthesis performance upon exposure to different light conditions. Significantly, the incompatible chloroplast genome is naturally found in xeric habitats with high light intensities, whereas the compatible one is limited to mesic habitats. Consequently, our data raise the possibility that the hybridization barrier evolved as a result of adaptation to specific climatic conditions.

The first report of RNA U to C or G editing in the mitochondrial NADH dehydrogenase subunit 5 (Nad5) transcript of wild barley

BACKGROUND

Nad dehydrogenase complex in mtDNA has a significant role in cellular respiration. One of the largest subunits in the complex is subunit 5 (Nad5).

METHODS AND RESULTS

Four cDNAs of the Hordeum vulgare subsp. spontaneum nad5 gene have been characterized and subjected to four phases of 0.5 M salinity, at 0 h (control, accession no. MT235236), after 2 h (acc. no. MT235237), after 12 h (acc. no. MT235238) and after 24 h (acc. no. MT235239). Utilizing raw data from RNA-seq, ten RNA editing sites were reported. Seven sites have common editing from C to U in positions (C1490, C1859, C1895, C1900, C1901, C1916, C1918). A rare editing event U to C was detected in two positions (U1650 and U1652) and a novel editing event U to G was for the first time in positions nad5-U231. The highest editing level was shown in 2 and 12 h after salinity exposure. After 24 h, these edits were disrupted, possibly due to the launch of the programed cell death mechanism. However, the RNA editing in positions U1650, U1652 and U231 was fixed at all exposure times.

CONCLUSIONS

Although study clarified the role of salinity stress in nad5 RNA editing sites, the main achievements are first report of U to G RNA editing in plants at position U231 and first report of U to C editing in the nad5 gene at U1650 and U1652.

The U-to-C RNA editing affects the mRNA stability of nuclear genes in Arabidopsis thaliana

Cytidine-to-uridine (C-to-U) RNA editing has been generally observed in land plants; however, reverse (U-to-C) RNA editing is a rare phenomenon. In this study, we investigated the U-to-C RNA editing-related genes in Arabidopsis tissues and the effects on mRNA stability, with a special focus on PPR proteins. A previous study showed the extensive occurrence of U-to-C RNA editing in 12-day and 20-dayold Arabidopsis seedlings. Here, we have demonstrated the effects of this "reverse" RNA editing on the mRNA stability for all seven edited genes. We also identified U-to-C RNA editing in the nuclear PPR gene (AT2G19280) in 12-day-old seedlings of Arabidopsis thaliana. The U-to-C RNA editing sites were found in the untranslated region (3' UTR) of the mature mRNA and may affect its secondary structure. We also examined the correlation between U-to-C RNA editing-related genes and their mRNA abundance. Furthermore, we investigated the effects of U-to-C RNA editing in Arabidopsis using the transcription inhibitor actinomycin D (Act D). The addition of Act D to the seedlings of transgenic Arabidopsis generated by Agrobacterium-mediated transformation showed that single nucleotide base conversion adversely affected the mRNA secondary structure and stability.

Plant super-barcode: a case study on genome-based identification for closely related species of Fritillaria

Background

Although molecular analysis offers a wide range of options for species identification, a universal methodology for classifying and distinguishing closely related species remains elusive. This study validated the effectiveness of utilizing the entire chloroplast (cp) genome as a super-barcode to help identify and classify closely related species.

Methods

We here compared 26 complete cp genomes of ten Fritillaria species including 18 new sequences sequenced in this study. Each species had repeats and the cp genomes were used as a whole DNA barcode to test whether they can distinguish Fritillaria species.

Results

The cp genomes of Fritillaria medicinal plants were conserved in genome structure, gene type, and gene content. Comparison analysis of the Fritillaria cp genomes revealed that the intergenic spacer regions were highly divergent compared with other regions. By constructing the phylogenetic tree by the maximum likelihood and maximum parsimony methods, we found that the entire cp genome showed a high discrimination power for Fritillaria species with individuals of each species in a monophyletic clade. These results indicate that cp genome can be used to effectively differentiate medicinal plants from the genus Fritillaria at the species level.

Conclusions

This study implies that cp genome can provide distinguishing differences to help identify closely related Fritillaria species, and has the potential to be served as a universal super-barcode for plant identification.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13020-021-00460-z.

RNA Editing Detection in HPC Infrastructures

RNA editing by A-to-I deamination is a relevant co/posttranscriptional modification carried out by ADAR enzymes. In humans, it has pivotal cellular effects and its deregulation has been linked to a variety of human disorders including neurological and neurodegenerative diseases and cancer. Despite its biological relevance, the detection of RNA editing variants in large transcriptome sequencing experiments (RNAseq) is yet a challenging computational task. To drastically reduce computing times we have developed a novel REDItools version able to identify A-to-I events in huge amount of RNAseq data employing High Performance Computing (HPC) infrastructures.Here we show how to use REDItools v2 in HPC systems.

Transcriptome and Comparative Chloroplast Genome Analysis of Vincetoxicum versicolor: Insights Into Molecular Evolution and Phylogenetic Implication

Vincetoxicum versicolor (Bunge) Decne is the original plant species of the Chinese herbal medicine Cynanchi Atrati Radix et Rhizoma. The lack of information on the transcriptome and chloroplast genome of V. versicolor hinders its evolutionary and taxonomic studies. Here, the V. versicolor transcriptome and chloroplast genome were assembled and functionally annotated. In addition, the comparative chloroplast genome analysis was conducted between the genera Vincetoxicum and Cynanchum. A total of 49,801 transcripts were generated, and 20,943 unigenes were obtained from V. versicolor. One thousand thirty-two unigenes from V. versicolor were classified into 73 functional transcription factor families. The transcription factors bHLH and AP2/ERF were the most significantly abundant, indicating that they should be analyzed carefully in the V. versicolor ecological adaptation studies. The chloroplast genomes of Vincetoxicum and Cynanchum exhibited a typical quadripartite structure with highly conserved gene order and gene content. They shared an analogous codon bias pattern in which the codons of protein-coding genes had a preference for A/U endings. The natural selection pressure predominantly influenced the chloroplast genes. A total of 35 RNA editing sites were detected in the V. versicolor chloroplast genome by RNA sequencing (RNA-Seq) data, and one of them restored the start codon in the chloroplast ndhD of V. versicolor. Phylogenetic trees constructed with protein-coding genes supported the view that Vincetoxicum and Cynanchum were two distinct genera.

Correlation between Organelle Genetic Variation and RNA Editing in Dinoflagellates Associated with the Coral Acropora digitifera

In order to develop successful strategies for coral reef preservation, it is critical that the biology of both host corals and symbiotic algae are investigated. In the Ryukyu Archipelago, which encompasses many islands spread over ∼500 km of the Pacific Ocean, four major populations of the coral Acropora digitifera have been studied using whole-genome shotgun (WGS) sequence analysis (Shinzato C, Mungpakdee S, Arakaki N, Satoh N. 2015. Genome-wide single-nucleotide polymorphism (SNP) analysis explains coral diversity and recovery in the Ryukyu Archipelago. Sci Rep. 5:18211.). In contrast, the diversity of the symbiotic dinoflagellates associated with these A. digitifera populations is unknown. It is therefore unclear if these two core components of the coral holobiont share a common evolutionary history. This issue can be addressed for the symbiotic algal populations by studying the organelle genomes of their mitochondria and plastids. Here, we analyzed WGS data from ∼150 adult A. digitifera, and by mapping reads to the available reference genome sequences, we extracted 2,250 sequences representing 15 organelle genes of Symbiodiniaceae. Molecular phylogenetic analyses of these mitochondrial and plastid gene sets revealed that A. digitifera from the southern Yaeyama islands harbor a different Symbiodiniaceae population than the islands of Okinawa and Kerama in the north, indicating that the distribution of symbiont populations partially matches that of the four host populations. Interestingly, we found that numerous SNPs correspond to known RNA-edited sites in 14 of the Symbiodiniaceae organelle genes, with mitochondrial genes showing a stronger correspondence than plastid genes. These results suggest a possible correlation between RNA editing and SNPs in the two organelle genomes of symbiotic dinoflagellates.

Atomic structures of respiratory complex III2, complex IV, and supercomplex III2-IV from vascular plants

Mitochondrial complex III (CIII2) and complex IV (CIV), which can associate into a higher-order supercomplex (SC III2+IV), play key roles in respiration. However, structures of these plant complexes remain unknown. We present atomic models of CIII2, CIV, and SC III2+IV from Vigna radiata determined by single-particle cryoEM. The structures reveal plant-specific differences in the MPP domain of CIII2 and define the subunit composition of CIV. Conformational heterogeneity analysis of CIII2 revealed long-range, coordinated movements across the complex, as well as the motion of CIII2's iron-sulfur head domain. The CIV structure suggests that, in plants, proton translocation does not occur via the H channel. The supercomplex interface differs significantly from that in yeast and bacteria in its interacting subunits, angle of approach and limited interactions in the mitochondrial matrix. These structures challenge long-standing assumptions about the plant complexes and generate new mechanistic hypotheses.

Uncovering dynamic evolution in the plastid genome of seven Ligusticum species provides insights into species discrimination and phylogenetic implications

Ligusticum L., one of the largest members in Apiaceae, encompasses medicinally important plants, the taxonomic statuses of which have been proved to be difficult to resolve. In the current study, the complete chloroplast genomes of seven crucial plants of the best-known herbs in Ligusticum were presented. The seven genomes ranged from 148,275 to 148,564 bp in length with a highly conserved gene content, gene order and genomic arrangement. A shared dramatic decrease in genome size resulted from a lineage-specific inverted repeat (IR) contraction, which could potentially be a promising diagnostic character for taxonomic investigation of Ligusticum, was discovered, without affecting the synonymous rate. Although a higher variability was uncovered in hotspot divergence regions that were unevenly distributed across the chloroplast genome, a concatenated strategy for rapid species identification was proposed because separate fragments inadequately provided variation for fine resolution. Phylogenetic inference using plastid genome-scale data produced a concordant topology receiving a robust support value, which revealed that L. chuanxiong had a closer relationship with L. jeholense than L. sinense, and L. sinense cv. Fuxiong had a closer relationship to L. sinense than L. chuanxiong, for the first time. Our results not only furnish concrete evidence for clarifying Ligusticum taxonomy but also provide a solid foundation for further pharmaphylogenetic investigation.

The novel E-subgroup pentatricopeptide repeat protein DEK55 is responsible for RNA editing at multiple sites and for the splicing of nad1 and nad4 in maize

BACKGROUND

Pentatricopeptide repeat (PPR) proteins compose a large protein family whose members are involved in both RNA processing in organelles and plant growth. Previous reports have shown that E-subgroup PPR proteins are involved in RNA editing. However, the additional functions and roles of the E-subgroup PPR proteins are unknown.

RESULTS

In this study, we developed and identified a new maize kernel mutant with arrested embryo and endosperm development, i.e., defective kernel (dek) 55 (dek55). Genetic and molecular evidence suggested that the defective kernels resulted from a mononucleotide alteration (C to T) at + 449 bp within the open reading frame (ORF) of Zm00001d014471 (hereafter referred to as DEK55). DEK55 encodes an E-subgroup PPR protein within the mitochondria. Molecular analyses showed that the editing percentage of 24 RNA editing sites decreased and that of seven RNA editing sites increased in dek55 kernels, the sites of which were distributed across 14 mitochondrial gene transcripts. Moreover, the splicing efficiency of nad1 introns 1 and 4 and nad4 intron 1 significantly decreased in dek55 compared with the wild type (WT). These results indicate that DEK55 plays a crucial role in RNA editing at multiple sites as well as in the splicing of nad1 and nad4 introns. Mutation in the DEK55 gene led to the dysfunction of mitochondrial complex I. Moreover, yeast two-hybrid assays showed that DEK55 interacts with two multiple organellar RNA-editing factors (MORFs), i.e., ZmMORF1 (Zm00001d049043) and ZmMORF8 (Zm00001d048291).

CONCLUSIONS

Our results demonstrated that a mutation in the DEK55 gene affects the mitochondrial function essential for maize kernel development. Our results also provide novel insight into the molecular functions of E-subgroup PPR proteins involved in plant organellar RNA processing.

The ADAR Family in Amphioxus: RNA Editing and Conserved Orthologous Site Predictions

RNA editing is a relatively unexplored process in which transcribed RNA is modified at specific nucleotides before translation, adding another level of regulation of gene expression. Cephalopods use it extensively to increase the regulatory complexity of their nervous systems, and mammals use it too, but less prominently. Nevertheless, little is known about the specifics of RNA editing in most of the other clades and the relevance of RNA editing from an evolutionary perspective remains unknown. Here we analyze a key element of the editing machinery, the ADAR (adenosine deaminase acting on RNA) gene family, in an animal with a key phylogenetic position at the root of chordates: the cephalochordate amphioxus. We show, that as in cephalopods, ADAR genes in amphioxus are predominantly expressed in the nervous system; we identify a number of RNA editing events in amphioxus; and we provide a newly developed method to identify RNA editing events in highly polymorphic genomes using orthology as a guide. Overall, our work lays the foundations for future comparative analysis of RNA-editing events across the metazoan tree.

Polyploidization is accompanied by synonymous codon usage bias in the chloroplast genomes of both cotton and wheat

Synonymous codon usage bias (SCUB) of both nuclear and organellar genes can mirror the evolutionary specialization of plants. The polyploidization process exposes the nucleus to genomic shock, a syndrome which promotes, among other genetic variants, SCUB. Its effect on organellar genes has not, however, been widely addressed. The present analysis targeted the chloroplast genomes of two leading polyploid crop species, namely cotton and bread wheat. The frequency of codons in the chloroplast genomes ending in either adenosine (NNA) or thymine (NNT) proved to be higher than those ending in either guanidine or cytosine (NNG or NNC), and this difference was conserved when comparisons were made between polyploid and diploid forms in both the cotton and wheat taxa. Preference for NNA/T codons was heterogeneous among genes with various numbers of introns and was also differential among the exons. SCUB patterns distinguished tetraploid cotton from its diploid progenitor species, as well as bread wheat from its diploid/tetraploid progenitor species, indicating that SCUB in the chloroplast genome partially mirrors the formation of polyploidies.

Plastomes of eight Ligusticum species: characterization, genome evolution, and phylogenetic relationships

BACKGROUND

The genus Ligusticum consists of approximately 60 species distributed in the Northern Hemisphere. It is one of the most taxonomically difficult taxa within Apiaceae, largely due to the varied morphological characteristics. To investigate the plastome evolution and phylogenetic relationships of Ligusticum, we determined the complete plastome sequences of eight Ligusticum species using a de novo assembly approach.

RESULTS

Through a comprehensive comparative analysis, we found that the eight plastomes were similar in terms of repeat sequence, SSR, codon usage, and RNA editing site. However, compared with the other seven species, L. delavayi exhibited striking differences in genome size, gene number, IR/SC borders, and sequence identity. Most of the genes remained under the purifying selection, whereas four genes showed relaxed selection, namely ccsA, rpoA, ycf1, and ycf2. Non-monophyly of Ligusticum species was inferred from the plastomes and internal transcribed spacer (ITS) sequences phylogenetic analyses.

CONCLUSION

The plastome tree and ITS tree produced incongruent tree topologies, which may be attributed to the hybridization and incomplete lineage sorting. Our study highlighted the advantage of plastome with mass informative sites in resolving phylogenetic relationships. Moreover, combined with the previous studies, we considered that the current taxonomy system of Ligusticum needs to be improved and revised. In summary, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of Ligusticum species.

Deciphering the organelle genomes and transcriptomes of a common ornamental plant Ligustrum quihoui reveals multiple fragments of transposable elements in the mitogenome

Ligustrum quihoui (L. quihoui) is an important hedge material for landscaping and also possesses medicinal value. To generate genomic resources for better understanding the evolutionary history of this important plant, the organelle genomes of L. quihoui are de novo assembled and functionally annotated. Compared with other Oleaceae species, the 163,069 bp chloroplast genome of L. quihoui exhibits a typical quadripartite structure with highly conserved gene content and gene order, while the 848,451 bp mitochondrial genome of L. quihoui exhibits highly divergent genome size and gene content. Codon usage analyses show that genes related with photosynthesis and mitochondrial respiratory chain show inconsistent codon biases. A total of 48,760 bp transposable elements (TEs) fragments and 41,887 bp chloroplast-like sequences are found in the L. quihoui mitochondrial genome. A striking discrepancy of RNA editing between the two organelle genomes is found in L. quihoui, in which 146 mitochondrial editing sites coexist with only 43 such sites in chloroplast. Based on DNA and RNA-Seq data, we propose that GTG may act as the start codon of mitochondrial rpl16 in Oleaceae species. Phylogenetic analysis based on chloroplast genome shows that L. quihoui and L. japonicum form a sister clade within the genus Ligustrum.

Fine mapping and candidate gene prediction of a major quantitative trait locus for tassel branch number in maize

Tassel branch number (TBN) is the principal component of tassel inflorescence architecture in the maize plant. TBN is believed to be controlled by a set of quantitative trait loci (QTLs). However, it is necessary to identify and genetically evaluate these QTLs before the TBN can be improved upon using a molecular breeding approach. Therefore, in this study, we developed the chromosome segment introgression line (CSIL) TBN1 with the Zong31 (Z31) background and a higher TBN, and then we utilized the CSIL-TBN1-derived populations and identified a major QTL, qTBN6a, by linkage analysis. Fine mapping of the qTBN6a QTL was validated using a set of sub-CSILs and located in a 240-kb genomic region (Bin6.07) in B73RefGen_v4. One allele included in the introgression fragment had a positive effect, noticeably increasing the TBN and demonstrating the potential to improve the TBN of Z31. Afterward, in the qTBN6a interval, gene expression, sequence alignment, functional analysis, and the analysis of motifs in the 5' UTR suggested that candidate genes of qTBN6a are important functional genes at the early stage of immature infected tassel development. Among these candidate genes, a long W22::Mu-insertion/deletion in exon one and an 11-bp insertion/deletion in the promoter region may affect the variation of the qTBN6a QTL observed between Z31 and TBN1. In addition, the candidate genes of qTBN6a were found to encode a pentatricopeptide repeat (PPR)-containing protein and a histone deacetylase (HDA), which are known to be closely associated with RNA editing and stability and chromatin state activity for the transcription of gene expression, respectively. Finally, a model of qTBN6a based on the synergistic regulation of PPR and HDA for the maintenance of inflorescence meristem (IM) identity and its differentiation to the branch meristem (BM) in TBN1 was suggested. Collectively, our results provide an available locus for the molecular improvement of TBN and the isolation of functional genes underlying this QTL.

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 complete plastome sequences of five Aponogeton species (Aponogetonaceae): Insights into the structural organization and mutational hotspots

Members of the aquatic plant genus Aponogeton are widely used commercially in aquariums because of their variable leaf shape and unique inflorescences. However, due to extensive similarity between species in this genus, morphological characters are generally inadequate for taxonomic classification. Currently, molecular makers available for taxonomic and phylogenetic studies of Aponogeton are limited. One approach to clarifying relationships between species in these complex groups is to use divergence hotspot regions within the genome. Here, we sequenced and analyzed the plastomes of five Aponogeton species collected from China, Zambia, and Kenya, and subsequently screened these plastomes for divergent DNA hotspots. The five plastomes are circular structures with sizes ranging from 154,167 bp to 154,860 bp. The Large and the Small Single Copies are separated by two Inverted Repeats. One hundred and thirteen unique genes were identified including 79 protein-coding, 30 tRNA, and four rRNA genes. We found that the most abundant repeats in all but one species were mononucleotide repeats (A/T) and that there were 23 potential RNA ending sites. Interestingly, a ~3 kb inversion, which includes the accD gene, was detected within the Asian species of Aponogeton. The inversion may be related to more frequent exchanges between this region and the nuclear genome. Furthermore, we detected mutational hotspot sites among the five Aponogeton species. Three of these hotspots are intergenic spacer regions (accD-psaI, rbcL-accD and trnH-GUG-psbA) that might be suitable for use as barcodes to resolve intra-generic relationships. We also identified four highly variable protein-coding genes (ccsA, rpl22, rps16 and ycf1) may be used as barcodes to resolve the higher-level phylogenies. Our study will provide valuable molecular resources for the taxonomic and phylogenomic study of the complex genus Aponogeton.

Extreme plastid RNA editing may confound phylogenetic reconstruction: A case study of Selaginella (lycophytes)

Cytidine-to-uridine (C-to-U) RNA editing is common in coding regions of organellar genomes throughout land plants. In most cases RNA editing alters translated amino acids or creates new start codons, potentially confounds phylogenetic reconstructions. In this study, we used the spike moss genus Selaginella (lycophytes), which has the highest frequency of RNA editing, as a model to test the effects of extreme RNA editing on phylogenetic reconstruction. We predicted the C-to-U RNA editing sites in coding regions of 18 Selaginella plastomes, and reconstructed the phylogenetic relationships within Selaginella based on three data set pairs consisted of plastome or RNA-edited coding sequences, first and second codon positions, and translated amino acid sequences, respectively. We predicted between 400 and 3100 RNA editing sites of 18 Selaginella plastomes. The numbers of RNA editing sites in plastomes were highly correlated with the GC content of first and second codon positions, but not correlated with the GC content of plastomes as a whole. Contrast phylogenetic analyses showed that there were substantial differences (e.g., the placement of clade B in Selaginella) between the phylogenies generated by the plastome and RNA-edited data sets. This empirical study provides evidence that extreme C-to-U RNA editing in the coding regions of organellar genomes alters the sequences used for phylogenetic reconstruction, and might even confound phylogenetic reconstruction. Therefore, RNA editing sites should be corrected when plastid or mitochondrial genes are used for phylogenetic studies, particularly in those lineages with abundant organellar RNA editing sites, such as hornworts, quillworts, spike mosses, and some seed plants.

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.

Emerging Roles of RNA-Binding Proteins in Seed Development and Performance

Seed development, dormancy, and germination are key physiological events that are not only important for seed generation, survival, and dispersal, but also contribute to agricultural production. RNA-binding proteins (RBPs) directly interact with target mRNAs and fine-tune mRNA metabolism by governing post-transcriptional regulation, including RNA processing, intron splicing, nuclear export, trafficking, stability/decay, and translational control. Recent studies have functionally characterized increasing numbers of diverse RBPs and shown that they participate in seed development and performance, providing significant insight into the role of RBP-mRNA interactions in seed processes. In this review, we discuss recent research progress on newly defined RBPs that have crucial roles in RNA metabolism and affect seed development, dormancy, and germination.

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.

Substitutional RNA Editing in Plant Organelles

RNA editing by cytidine (C) to uridine (U) conversions frequently occurs in land plant mitochondria and plastids. Target cytidines are specifically recognized by nuclear-encoded pentatricopeptide repeat (PPR) proteins in a sequence-specific manner. In the moss Physcomitrella patens, all PPR editing factors possess the DYW-deaminase domain at the C-terminus. Here, we describe methods for the direct sequencing of cDNA to detect RNA editing events and the RNA electrophoresis mobility shift assay (REMSA) to analyze the specific binding of PPR editing factors to their target RNA.

High level of intraspecific divergence and low frequency of RNA editing in the chloroplast genome sequence of Tagetes erecta

Tagetes erecta L. is an important commercial and medicinal plant. In this study, we reported the complete chloroplast genome sequence of T. erecta. The genome has a circular structure of 152,076 bp containing a large single-copy region (LSC) of 83,914 bp, a small copy region (SSC) of 18,064 bp, and two inverted repeats (IR) of 25,049 bp by each. It harbors 111 unique genes, including 79 protein-coding genes, 4 ribosomal RNA genes, and 28 transfer RNA genes. A total of 41 microsatellite, 20 tandem, and 37 interspersed repeats were detected in the genome. The phylogenomic analysis shows that T. erecta is a single phylogenetic cluster. The complete chloroplast genome of T. erecta lays the foundation for the phylogenetic, evolutionary, and conservation studies of the genus Tagetes. Furthermore, the intergenic region of atpB-rbcL was variable among the species T. erecta. This suggests that this region might be a mutation hotspot and will be useful for phylogenetic study and the development of molecular markers. At last, we systematically identified the RNA editing sites in the chloroplast genome of T. erecta based on the transcriptome downloaded from the SRA database. This study identified the characteristics of the T. erecta chloroplast genome, SNPs, and RNA editing sites, which will facilitate species identification and phylogenetic analysis within T. erecta.

Disruption of a Rice Chloroplast-Targeted Gene OsHMBPP Causes a Seedling-Lethal Albino Phenotype

BACKGROUND

Chloroplast development is coordinately regulated by plastid- and nuclear-encoding genes. Although many regulators have been reported to be involved in chloroplast development, new factors remain to be identified, given the complexity of this process.

RESULTS

In this study, we characterized a rice mutant lethal albinic seedling 1(las1)form of a 4-hydroxy-3-methylbut-2-enyl diphosphate reductase (OsHMBPP) that was targeted to the chloroplasts. The LAS1 mutation caused the albino lethal phenotype in seedlings. Transmission electron microscopy indicated that las1 were defective in early chloroplast development. LAS1 is preferentially expressed in leaves, implying its role in controlling chloroplast development. The expression levels of many chloroplast-encoded genes were altered significantly in las1. The expression levels of nuclear-encoded gene involved in Chl biosynthesis were also decreased in las1. We further investigated plastidic RNA editing in las1 and found that the edit efficiency of four chloroplast genes were markly altered. Compared with WT, las1 exhibited defective in biogenesis of chloroplast ribosomes.

CONCLUSIONS

Our results show that LAS1/OsHMBPP plays an essential role in the early chloroplast development in rice.

Crystal structure of Arabidopsis thaliana cytidine deaminase

Cytidine deaminase (CDA) catalyzes the (deoxy)cytidine deamination to (deoxy)uridine, which involves in the catabolic and salvage pathways of pyrimidine nucleotides in plants. CDA serves as a prototype of the cytidine deaminase superfamily that contains a number of RNA editing enzymes. Arabidopsis thaliana has only one functional CDA, AtCDA1. We solved the crystal structures of AtCDA1, which is a dimeric zinc-containing enzyme and each protomer consists of an N-terminal zinc-binding catalytic domain and a C-terminal non-catalytic domain. Both domains adopt a typical α/β/α sandwich fold. In vitro biochemical assays showed that the ribose moiety of cytidine is required for ligand binding, and structural analyses revealed a conserved catalytic mechanism is adopted by AtCDA1.

Phylogenetic and evolutionary features of the plastome of Tropaeolum pentaphyllum Lam. (Tropaeolaceae)

Complete plastome sequence of Tropaeolum pentaphyllum revealed molecular markers, hotspots of nucleotide polymorphism, RNA editing sites and phylogenetic aspects Tropaeolaceae Juss. ex DC. comprises approximately 95 species across North and South Americas. Tropaeolum pentaphyllum Lam. is an unconventional and endangered species with occurrence in some countries of South America. Although this species presents nutritional, medicinal and ornamental uses, genetic studies involving natural populations or promising genotypes are practically non-existent. Here, we report the nucleotide sequence of T. pentaphyllum plastome. It represents the first complete plastome sequence of the family Tropaeolaceae to be fully sequenced and analyzed in detail. The sequencing data revealed that the T. pentaphyllum plastome is highly similar to the plastomes of other Brassicales. Notwithstanding, our analyses detected some specific features concerning events of IR expansion and structural changes in some genes such as matK, rpoA, and rpoC2. We also detected 251 SSR loci, nine hotspots of nucleotide polymorphism, and two specific RNA editing sites in the plastome of T. pentaphyllum. Moreover, plastid phylogenomic inference indicated a closed relationship between the families Tropaeolaceae and Akaniaceae, which formed a sister group to Moringaceae-Caricaceae. Finally, our data bring new molecular markers and evolutionary features to be applied in the natural population, germplasm collection, and genotype selection aiming conservation, genetic diversity evaluation, and exploitation of this endangered species.

The Mitochondrial Genome of Eleusine indica and Characterization of Gene Content within Poaceae

Plant mitochondrial (mt) genome assembly provides baseline data on size, structure, and gene content, but resolving the sequence of these large and complex organelle genomes remains challenging due to fragmentation, frequent recombination, and transfers of DNA from neighboring plastids. The mt genome for Eleusine indica (Poaceae: goosegrass) is comprehensibly analyzed here, providing key reference data for an economically significant invasive species that is also the maternal parent of the allotetraploid crop Finger millet (Eleusine coracana). The assembled E. indica genome contains 33 protein coding genes, 6 rRNA subunits, 24 tRNA, 8 large repetitive regions 15 kb of transposable elements across a total of 520,691 bp. Evidence of RNA editing and loss of rpl2, rpl5, rps14, rps11, sdh4, and sdh3 genes is evaluated in the context of an updated survey of mt genomic gene content across the grasses through an analysis of publicly available data. Hypothesized patterns of Poaceae mt gene loss are examined in a phylogenetic context to clarify timing, showing that rpl2 was transferred to the nucleus from the mitochondrion prior to the origin of the PACMAD clade.

Mutation of YL Results in a Yellow Leaf with Chloroplast RNA Editing Defect in Soybean

RNA editing plays a key role in organelle gene expression. Little is known about how RNA editing factors influence soybean plant development. Here, we report the isolation and characterization of a soybean yl (yellow leaf) mutant. The yl plants showed decreased chlorophyll accumulation, lower PS II activity, an impaired net photosynthesis rate, and an altered chloroplast ultrastructure. Fine mapping of YL uncovered a point mutation in Glyma.20G187000, which encodes a chloroplast-localized protein homologous to Arabidopsis thaliana (Arabidopsis) ORRM1. YL is mainly expressed in trifoliate leaves, and its deficiency affects the editing of multiple chloroplast RNA sites, leading to inferior photosynthesis in soybean. Taken together, these results demonstrate the importance of the soybean YL protein in chloroplast RNA editing and photosynthesis.

Gene editing: an instrument for practical application of gene biology to plant breeding

Plant breeding is well recognized as one of the most important means to meet food security challenges caused by the ever-increasing world population. During the past three decades, plant breeding has been empowered by both new knowledge on trait development and regulation (e.g., functional genomics) and new technologies (e.g., biotechnologies and phenomics). Gene editing, particularly by clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein (Cas) and its variants, has become a powerful technology in plant research and may become a game-changer in plant breeding. Traits are conferred by coding and non-coding genes. From this perspective, we propose different editing strategies for these two types of genes. The activity of an encoded enzyme and its quantity are regulated at transcriptional and post-transcriptional, as well as translational and post-translational, levels. Different strategies are proposed to intervene to generate gene functional variations and consequently phenotype changes. For non-coding genes, trait modification could be achieved by regulating transcription of their own or target genes via gene editing. Also included is a scheme of protoplast editing to make gene editing more applicable in plant breeding. In summary, this review provides breeders with a host of options to translate gene biology into practical breeding strategies, i.e., to use gene editing as a mechanism to commercialize gene biology in plant breeding.

SMALL KERNEL4 is required for mitochondrial cox1 transcript editing and seed development in maize

In land plants, cytidine-to-uridine (C-to-U) editing of organellar transcripts is an important post-transcriptional process, which is considered to remediate DNA genetic mutations to restore the coding of functional proteins. Pentatricopeptide repeat (PPR) proteins have key roles in C-to-U editing. Owing to its large number, however, the biological functions of many PPR proteins remain to be identified. Through characterizing a small kernel4 (smk4) mutant, here we report the function of Smk4 and its role in maize growth and development. Null mutation of Smk4 slows plant growth and development, causing small plants, delayed flowering time, and small kernels. Cloning revealed that Smk4 encodes a new E-subclass PPR protein, and localization indicated that SMK4 is exclusively localized in mitochondria. Loss of Smk4 function abolishes C-to-U editing at position 1489 of the cytochrome c oxidase1 (cox1) transcript, causing an amino acid change from serine to proline at 497 in Cox1. Cox1 is a core component of mitochondrial complex IV. Indeed, complex IV activity is reduced in the smk4, along with drastically elevated expression of alternative oxidases (AOX). These results indicate that SMK4 functions in the C-to-U editing of cox1-1489, and this editing is crucial for mitochondrial complex IV activity, plant growth, and kernel development in maize.

MEF46 and MEF47 are novel specificity factors for RNA editing sites in mitochondrial nad transcripts

Terrestrial plants have C-to-U RNA editing in the transcripts of plastids and mitochondria. Target specificity for more than several hundred editing sites are governed by PLS (PPR, Long and Short) class Pentatricopeptide repeat (PPR) proteins with additional C-terminal domains. Half of these PPR proteins have DYW (Aspartate (D), Tyrosine (Y) and Tryptophan (W)) domains, which most likely harbour cytidine deaminase activity. The other half of them, E subclass and E+ subclass proteins, contain no or only a part of the DYW domain. Missing DYW domains in the E and E+ subclass PPR proteins are likely to be complemented by other DYW containing proteins. All target sites of so far characterized E+ subclass PPR proteins show defects in dyw2 mutants, suggesting that the DYW2 protein complements the missing DYW domains in the E+ subclass PPR proteins. Here we report two novel RNA editing factors, MEF46 and MEF47, which belong to E+ and E subclass, respectively. The defective editing site in mef46, nad5-1958, overlaps with the affected sites in dyw2 mutants, while that in mef47, nad3-64 and ccmC-614 do not, further supporting the specific functional connection between E+ subclass PPR proteins and DYW2.