RNA Editing and Its Molecular Mechanism in Plant Organelles

Seeing the unseen in characterizing RNA editome during rice endosperm development

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

The pentatricopeptide repeat protein DG1 promotes the transition to bilateral symmetry during Arabidopsis embryogenesis through GUN1-mediated plastid signals

During Arabidopsis embryogenesis, the transition of the embryo's symmetry from radial to bilateral between the globular and heart stage is a crucial event, involving the formation of cotyledon primordia and concurrently the establishment of a shoot apical meristem (SAM). However, a coherent framework of how this transition is achieved remains to be elucidated. In this study, we investigated the function of DELAYED GREENING 1 (DG1) in Arabidopsis embryogenesis using a newly identified dg1-3 mutant. The absence of chloroplast-localized DG1 in the mutants led to embryos being arrested at the globular or heart stage, accompanied by an expansion of WUSCHEL (WUS) and SHOOT MERISTEMLESS (STM) expression. This finding pinpoints the essential role of DG1 in regulating the transition to bilateral symmetry. Furthermore, we showed that this regulation of DG1 may not depend on its role in plastid RNA editing. Nevertheless, we demonstrated that the DG1 function in establishing bilateral symmetry is genetically mediated by GENOMES UNCOUPLED 1 (GUN1), which represses the transition process in dg1-3 embryos. Collectively, our results reveal that DG1 functionally antagonizes GUN1 to promote the transition of the Arabidopsis embryo's symmetry from radial to bilateral and highlight the role of plastid signals in regulating pattern formation during plant embryogenesis.

Comparative mitochondrial genomics of Terniopsis yongtaiensis in Malpighiales: structural, sequential, and phylogenetic perspectives

BACKGROUND

Terniopsis yongtaiensis, a member of the Podostemaceae family, is an aquatic flowering plant displaying remarkable adaptive traits that enable survival in submerged, turbulent habitats. Despite the progressive expansion of chloroplast genomic information within this family, mitochondrial genome sequences have yet to be reported.

RESULTS

In current study, the mitochondrial genome of the T. yongtaiensis was characterized by a circular genome of 426,928 bp encoding 31 protein-coding genes (PCGs), 18 tRNAs, and 3 rRNA genes. Our comprehensive analysis focused on gene content, repeat sequences, RNA editing processes, intracellular gene transfer, phylogeny, and codon usage bias. Numerous repeat sequences were identified, including 130 simple sequence repeats, 22 tandem repeats, and 220 dispersed repeats. Phylogenetic analysis positioned T. yongtaiensis (Podostemaceae) within the Malpighiales order, showing a close relationship with the Calophyllaceae family, which was consistent with the APG IV classification. A comparative analysis with nine other Malpighiales species revealed both variable and conserved regions, providing insights into the genomic evolution within this order. Notably, the GC content of T. yongtaiensis was distinctively lower compared to other Malpighilales, primarily due to variations in non-coding regions and specific protein-coding genes, particularly the nad genes. Remarkably, the number of RNA editing sites was low (276), distributed unevenly across 27 PCGs. The dN/dS analysis showed only the ccmB gene of T. yongtaiensis was positively selected, which plays a crucial role in cytochrome c biosynthesis. Additionally, there were 13 gene-containing homologous regions between the mitochondrial and chloroplast genomes of T. yongtaiensis, suggesting the gene transfer events between these organellar genomes.

CONCLUSIONS

This study assembled and annotated the first mitochondrial genome of the Podostemaceae family. The comparison results of mitochondrial gene composition, GC content, and RNA editing sites provided novel insights into the adaptive traits and genetic reprogramming of this aquatic eudicot group and offered a foundation for future research on the genomic evolution and adaptive mechanisms of Podostemaceae and related plant families in the Malpighiales order.

Screening clusters of charged residues in plants' mitochondrial proteins and biological significance

Protein function is dependent on charge interactions and charge biased regions, which are involved in a wide range of cellular and biochemical processes. We report the development of a new algorithm implemented in Python and its use to identify charge clusters CC (NegativeCC: NCC, PositiveCC: PCC and MixedCC: MCC) and compare their presence in mitochondrial proteins of plant groups. To characterize the resulting CC, statistical, structural and functional analyses were conducted. The screening of 105 399 protein sequences showed that 2.6 %, 0.48 % and 0.03 % of the proteins contain NCC, PCC and MCC, respectively. Mitochondrial proteins encoded by the nuclear genome of green algae have the biggest proportion of both PCC (1.6 %) and MCC (0.4 %) and mitochondrial proteins coded by the nuclear genome of other plants group have the highest portion of NCC (7.5 %). The mapping of the identified CC showed that they are mainly located in the terminal regions of the protein. Annotation showed that proteins with CC are classified as binding proteins, are included in the transmembrane transport processes, and are mainly located in the membrane. The CC scanning revealed the presence of 2373 and 784 sites and 192 and 149 motif profiles within NCC and PCC, respectively. The investigation of CC within pentatricopeptide repeat-containing proteins revealed that they are involved in correct and specific RNA editing. CC were proven to play a key role in providing insightful structural and functional information of complex protein assemblies which could be useful in biotechnological applications.

Assembly and comparative analysis of the complete mitochondrial and chloroplast genome of Cyperus stoloniferus (Cyperaceae), a coastal plant possessing saline-alkali tolerance

BACKGROUND

Cyperus stoloniferus is an important species in coastal ecosystems and possesses economic and ecological value. To elucidate the structural characteristics, variation, and evolution of the organelle genome of C. stoloniferus, we sequenced, assembled, and compared its mitochondrial and chloroplast genomes.

RESULTS

We assembled the mitochondrial and chloroplast genomes of C. stoloniferus. The total length of the mitochondrial genome (mtDNA) was 927,413 bp, with a GC content of 40.59%. It consists of two circular DNAs, including 37 protein-coding genes (PCGs), 22 tRNAs, and five rRNAs. The length of the chloroplast genome (cpDNA) was 186,204 bp, containing 93 PCGs, 40 tRNAs, and 8 rRNAs. The mtDNA and cpDNA contained 81 and 129 tandem repeats, respectively, and 346 and 1,170 dispersed repeats, respectively, both of which have 270 simple sequence repeats. The third high-frequency codon (RSCU > 1) in the organellar genome tended to end at A or U, whereas the low-frequency codon (RSCU < 1) tended to end at G or C. The RNA editing sites of the PCGs were relatively few, with only 9 and 23 sites in the mtDNA and cpDNA, respectively. A total of 28 mitochondrial plastid DNAs (MTPTs) in the mtDNA were derived from cpDNA, including three complete trnT-GGU, trnH-GUG, and trnS-GCU. Phylogeny and collinearity indicated that the relationship between C. stoloniferus and C. rotundus are closest. The mitochondrial rns gene exhibited the greatest nucleotide variability, whereas the chloroplast gene with the greatest nucleotide variability was infA. Most PCGs in the organellar genome are negatively selected and highly evolutionarily conserved. Only six mitochondrial genes and two chloroplast genes exhibited Ka/Ks > 1; in particular, atp9, atp6, and rps7 may have undergone potential positive selection.

CONCLUSION

We assembled and validated the mtDNA of C. stoloniferus, which contains a 15,034 bp reverse complementary sequence. The organelle genome sequence of C. stoloniferus provides valuable genomic resources for species identification, evolution, and comparative genomic research in Cyperaceae.

De Novo Hybrid Assembly Unveils Multi-Chromosomal Mitochondrial Genomes in Ludwigia Species, Highlighting Genomic Recombination, Gene Transfer, and RNA Editing Events

Biological invasions have been identified as the fifth cause of biodiversity loss, and their subsequent dispersal represents a major ecological challenge. The aquatic invasive species Ludwigia grandiflora subsp. hexapetala (Lgh) and Ludwigia peploides subsp. montevidensis (Lpm) are largely distributed in aquatic environments in North America and in Europe. However, they also present worrying terrestrial forms that are able to colonize wet meadows. To comprehend the mechanisms of the terrestrial adaptation of Lgh and Lpm, it is necessary to develop their genomic resources, which are currently poorly documented. We performed de novo assembly of the mitogenomes of Lgh and Lpm through hybrid assemblies, combining short reads (SR) and/or long reads (LR) before annotating both mitogenomes. We successfully assembled the mitogenomes of Lgh and Lpm into two circular molecules each, resulting in a combined total length of 711,578 bp and 722,518 bp, respectively. Notably, both the Lgh and Lpm molecules contained plastome-origin sequences, comprising 7.8% of the mitochondrial genome length. Additionally, we identified recombinations that were mediated by large repeats, suggesting the presence of multiple alternative conformations. In conclusion, our study presents the first high-quality mitogenomes of Lpm and Lgh, which are the only ones in the Myrtales order found as two circular molecules.

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.

Conquering new grounds: plant organellar C-to-U RNA editing factors can be functional in the plant cytosol

Plant mitochondrial and chloroplast transcripts are subject to numerous events of specific cytidine-to-uridine (C-to-U) RNA editing to correct genetic information. Key protein factors for this process are specific RNA-binding pentatricopeptide repeat (PPR) proteins, which are encoded in the nucleus and post-translationally imported into the two endosymbiotic organelles. Despite hundreds of C-to-U editing sites in the plant organelles, no comparable editing has been found for nucleo-cytosolic mRNAs raising the question why plant RNA editing is restricted to chloroplasts and mitochondria. Here, we addressed this issue in the model moss Physcomitrium patens, where all PPR-type RNA editing factors comprise specific RNA-binding and cytidine deamination functionalities in single proteins. To explore whether organelle-type RNA editing can principally also take place in the plant cytosol, we expressed PPR56, PPR65 and PPR78, three editing factors recently shown to also function in a bacterial setup, together with cytosolic co-transcribed native targets in Physcomitrium. While we obtained unsatisfying results upon their constitutive expression, we found strong cytosolic RNA editing under hormone-inducible expression. Moreover, RNA-Seq analyses revealed varying numbers of up to more than 900 off-targets in other cytosolic transcripts. We conclude that PPR-mediated C-to-U RNA editing is not per se incompatible with the plant cytosol but that its limited target specificity has restricted its occurrence to the much less complex transcriptomes of mitochondria and chloroplast in the course of evolution.

Piece and parcel of gymnosperm organellar genomes

MAIN CONCLUSION

Significant past, present, and potential future research into the organellar (plastid and mitochondrial) genomes of gymnosperms that can provide insight into the unknown origin and evolution of plants is highlighted. Gymnosperms are vascular seed plants that predominated the ancient world before their sister clade, angiosperms, took over during the Late Cretaceous. The divergence of gymnosperms and angiosperms took place around 300 Mya, with the latter evolving into the diverse group of flowering plants that dominate the plant kingdom today. Although gymnosperms have reportedly made some evolutionary innovations, the literature on their genome advances, particularly their organellar (plastid and mitochondrial) genomes, is relatively scattered and fragmented. While organellar genomes can shed light on plant origin and evolution, they are frequently overlooked, due in part to their limited contribution to gene expression and lack of evolutionary dynamics when compared to nuclear genomes. A better understanding of gymnosperm organellar genomes is critical because they reveal genetic changes that have contributed to their unique adaptations and ecological success, potentially aiding in plant survival, enhancement, and biodiversity conservation in the face of climate change. This review reveals significant information and gaps in the existing knowledge base of organellar genomes in gymnosperms, as well as the challenges and research needed to unravel their complexity.

Comprehensive identification, characterization, and expression analysis of the MORF gene family in Brassica napus

BACKGROUND

RNA editing in chloroplast and mitochondrion transcripts of plants is an important type of post-transcriptional RNA modification in which members of the multiple organellar RNA editing factor gene family (MORF) play a crucial role. However, a systematic identification and characterization of MORF members in Brassica napus is still lacking.

RESULTS

In this study, a total of 43 MORF genes were identified from the genome of the Brassica napus cultivar "Zhongshuang 11". The Brassica napus MORF (BnMORF) family members were divided into three groups through phylogenetic analysis. BnMORF genes distributed on 14 chromosomes and expanded due to segmental duplication and whole genome duplication repetitions. The majority of BnMORF proteins were predicted to be localized to mitochondria and chloroplasts. The promoter cis-regulatory element analysis, spatial-temporal expression profiling, and co-expression network of BnMORF genes indicated the involvement of BnMORF genes in stress and phytohormone responses, as well as growth and development.

CONCLUSION

This study provides a comprehensive analysis of BnMORF genes and lays a foundation for further exploring their physiological functions in Brassica napus.

Multiple factors interact in editing of PPR-E+-targeted sites in maize mitochondria and plastids

RNA cytidine-to-uridine editing is essential for plant organellar gene expression. Pentatricopeptide repeat (PPR)-E+ proteins have been proposed to bind to target sites and recruit the cytidine deaminase AtDYW2, facilitated by AtNUWA. Here we analyze the function of ZmNUWA, ZmDYW2A, and ZmDYW2B and their relationships with other editing factors in maize. The zmdyw2a and zmdyw2b single mutants are normal, but the zmdyw2a::zmdyw2b and zmnuwa mutants are severely arrested in seed development. ZmNUWA, ZmDYW2A, and ZmDYW2B are dual localized in mitochondria and plastids. Loss of ZmNUWA decreases the editing at 99 mitochondrial sites and 8 plastid sites. Surprisingly, loss of ZmDYW2A:ZmDYW2B affects almost the same set of sites targeted by PPR-E+ proteins. ZmNUWA interacts with ZmDYW2A and ZmDYW2B, suggesting that ZmNUWA recruits ZmDYW2A/2B in the editing of PPR-E+-targeted sites in maize. Further protein interaction analyses show that ZmNUWA and ZmDYW2A/2B interact with ZmMORF1, ZmMORF8, ZmMORF2, and ZmMORF9 and that ZmOZ1 interacts with ZmORRM1, ZmDYW2A, ZmDYW2B, ZmMORF8, and ZmMORF9. These results suggest that the maize mitochondrial PPR-E+ editosome contains PPR-E+, ZmDYW2A/2B, ZmNUWA, and ZmMORF1/8, whereas the plastid PPR-E+ editosome is composed of PPR-E+, ZmDYW2A/2B, ZmNUWA, ZmMORF2/8/9, ZmORRM1, and ZmOZ1.

Assembly and comparative analysis of the initial complete mitochondrial genome of Primulina hunanensis (Gesneriaceae): a cave-dwelling endangered plant

BACKGROUND

Primulina hunanensis, a troglobitic plant within the Primulina genus of Gesneriaceae family, exhibits robust resilience to arid conditions and holds great horticultural potential as an ornamental plant. The work of chloroplast genome (cpDNA) has been recently accomplished, however, the mitochondrial genome (mtDNA) that is crucial for plant evolution has not been reported.

RESULTS

In this study, we sequenced and assembled the P. hunanensis complete mtDNA, and elucidated its evolutionary and phylogenetic relationships. The assembled mtDNA spans 575,242 bp with 43.54% GC content, encompassing 60 genes, including 37 protein-coding genes (PCGs), 20 tRNA genes, and 3 rRNA genes. Notably, high number of repetitive sequences in the mtDNA and substantial sequence translocation from chloroplasts to mitochondria were observed. To determine the evolutionary and taxonomic positioning of P. hunanensis, a phylogenetic tree was constructed using mitochondrial PCGs from P. hunanensis and 32 other taxa. Furthermore, an exploration of PCGs relative synonymous codon usage, identification of RNA editing events, and an investigation of collinearity with closely related species were conducted.

CONCLUSIONS

This study reports the initial assembly and annotation of P. hunanensis mtDNA, contributing to the limited mtDNA repository for Gesneriaceae plants and advancing our understanding of their evolution for improved utilization and conservation.

The first mitochondrial genome of Calophyllum soulattri Burm.f

Calophyllum soulattri Burm.f. is traditionally used to treat skin infections and reduce rheumatic pain, yet genetic and genomic studies are still limited. Here, we present the first complete mitochondrial genome of C. soulattri. It is 378,262 bp long with 43.97% GC content, containing 55 genes (30 protein-coding, 5 rRNA, and 20 tRNA). Repeat analysis of the mitochondrial genome revealed 194 SSRs, mostly mononucleotides, and 266 pairs of dispersed repeats ( ≥ 30 bp) that were predominantly palindromic. There were 23 homologous fragments found between the mitochondrial and plastome genomes. We also predicted 345 C-to-U RNA editing sites from 30 protein-coding genes (PCGs) of the C. soulatrii mitochondrial genome. These RNA editing events created the start codon of nad1 and the stop codon of ccmFc. Most PCGs of the C. soulattri mitochondrial genome underwent negative selection, but atp4 and ccmB experienced positive selection. Phylogenetic analyses showed C. soulattri is a sister taxon of Garcinia mangostana. This study has shed light on C. soulattri's evolution and Malpighiales' phylogeny. As the first complete mitochondrial genome in Calophyllaceae, it can be used as a reference genome for other medicinal plant species within the family for future genetic studies.

DEFECTIVE KERNEL 56 functions in mitochondrial RNA editing and maize seed development

Proper seed development is essential for achieving grain production, successful seed germination, and seedling establishment in maize (Zea mays). In the past few decades, pentatricopeptide repeat (PPR) proteins have been proven to play an essential role in regulating the development of maize kernels through posttranscriptional RNA modification of mitochondrial genes. However, the underlying mechanisms remain largely unknown. Here, we characterized a mutant of DEFECTIVE KERNEL 56 (DEK56) with defective kernels that exhibited arrested development of both the embryo and endosperm. Accordingly, we isolated DEK56 through a map-based cloning strategy and found that it encoded an E subgroup PPR protein located in the mitochondria. Dysfunction of DEK56 resulted in altered cytidine (C)-to-uridine (U) editing efficiency at 48 editing sites across 21 mitochondrial transcripts. Notably, the editing efficiency of the maturase-related (matR)-1124 site was substantially reduced or abolished in the dek56 mutant. Furthermore, we found that the splicing efficiency of NADH dehydrogenase subunit 4 (nad4) Introns 1 and 3 was substantially reduced in dek56 kernels, which might be a consequence of the defective MatR function. Through a protein-protein interaction test, we hypothesized that DEK56 carries out its function by recruiting the PPR-DYW protein PPR motif, coiled-coil, and DYW domain-containing protein 1 (PCW1). This interaction is facilitated by Multiple Organellar RNA Editing Factors (ZmMORFs) and Glutamine-Rich Protein 23 (ZmGRP23). Based on these findings, we developed a working model of PPR-mediated mitochondrial processing that plays an essential role in the development of maize kernels. The present study will further broaden our understanding of PPR-mediated seed development and provide a theoretical basis for maize improvement.

Assembly and comparative analysis of the complete mitochondrial genome of Pinellia ternata

Pinellia ternata is an important natural medicinal herb in China. However, it is susceptible to withering when exposed to high temperatures during growth, which limits its tuber production. Mitochondria usually function in stress response. The P . ternata mitochondrial (mt) genome has yet to be explored. Therefore, we integrated PacBio and Illumina sequencing reads to assemble and annotate the mt genome of P . ternata . The circular mt genome of P . ternata is 876 608bp in length and contains 38 protein-coding genes (PCGs), 20 tRNA genes and three rRNA genes. Codon usage, sequence repeats, RNA editing and gene migration from chloroplast (cp) to mt were also examined. Phylogenetic analysis based on the mt genomes of P . ternata and 36 other taxa revealed the taxonomic and evolutionary status of P . ternata . Furthermore, we investigated the mt genome size and GC content by comparing P . ternata with the other 35 species. An evaluation of non-synonymous substitutions and synonymous substitutions indicated that most PCGs in the mt genome underwent negative selection. Our results provide comprehensive information on the P . ternata mt genome, which may facilitate future research on the high-temperature response of P . ternata and provide new molecular insights on the Araceae family.

Complete mitochondrial genome of Syzygium samarangense reveals genomic recombination, gene transfer, and RNA editing events

Wax apple (Syzygium samarangense) is a commercial fruit that belongs to one of the most species-rich tree genera in the world. We report here the first complete S. samarangense mitogenome obtained using a hybrid assembly strategy. The mitogenome was a 530,242 bp circular molecule encoding 61 unique genes accounting for 7.99% of the full-length genome. Additionally, 167 simple sequence repeats, 19 tandem repeats, and 529 pairs of interspersed repeats were identified. Long read mapping and Sanger sequencing revealed the involvement of two forward repeats (35,843 bp and 22,925 bp) in mediating recombination. Thirteen homologous fragments in the chloroplast genome were identified, accounting for 1.53% of the mitogenome, and the longest fragment was 2,432 bp. An evolutionary analysis showed that S. samarangense underwent multiple genomic reorganization events and lost at least four protein-coding genes (PCGs) (rps2, rps7, rps11, and rps19). A total of 591 RNA editing sites were predicted in 37 PCGs, of which nad1-2, nad4L-2, and rps10-2 led to the gain of new start codons, while atp6-1156, ccmFC-1315 and rps10-331 created new stop codons. This study reveals the genetic features of the S. samarangense mitogenome and provides a scientific basis for further studies of traits with an epistatic basis and for germplasm identification.

Mitochondrial RNA maturation

The vast majority of oxygen-utilizing eukaryotes need to express their own mitochondrial genome, mtDNA, to survive. In comparison to size of their nuclear genome, mtDNA is minimal, even in the most exceptional examples. Having evolved from bacteria in an endosymbiotic event, it might be expected that the process of mtDNA expression would be relatively simple. The aim of this short review is to illustrate just how wrong this assumption is. The production of functional mitochondrial RNA across species evolved in many directions. Organelles use a dizzying array of RNA processing, modifying, editing, splicing and maturation events that largely require the import of nuclear-encoded proteins from the cytosol. These processes are sometimes driven by the unusual behaviour of the mitochondrial genome from which the RNA is originally transcribed, but in many examples the complex processes that are essential for the production of functional RNA in the organelle, are fascinating and bewildering.

Codon usage bias in chloroplast genes implicate adaptive evolution of four ginger species

Codon usage bias (CUB) refers to different codons exhibiting varying frequencies of usage in the genome. Studying CUB is crucial for understanding genome structure, function, and evolutionary processes. Herein, we investigated the codon usage patterns and influencing factors of protein-coding genes in the chloroplast genomes of four sister genera (monophyletic Roscoea and Cautleya, and monophyletic Pommereschea and Rhynchanthus) from the Zingiberaceae family with contrasting habitats in southwestern China. These genera exhibit distinct habitats, providing a unique opportunity to explore the adaptive evolution of codon usage. We conducted a comprehensive analysis of nucleotide composition and codon usage on protein-coding genes in the chloroplast genomes. The study focused on understanding the relationship between codon usage and environmental adaptation, with a particular emphasis on genes associated with photosynthesis. Nucleotide composition analysis revealed that the overall G/C content of the coding genes was ˂ 48%, indicating an enrichment of A/T bases. Additionally, synonymous and optimal codons were biased toward ending with A/U bases. Natural selection is the primary factor influencing CUB characteristics, particularly photosynthesis-associated genes. We observed differential gene expressions related to light adaptation among sister genera inhabiting different environments. Certain codons were favored under specific conditions, possibly contributing to gene expression regulation in particular environments. This study provides insights into the adaptive evolution of these sister genera by analyzing CUB and offers theoretical assistance for understanding gene expression and regulation. In addition, the data support the relationship between RNA editing and CUB, and the findings shed light on potential research directions for investigating adaptive evolution.

Current technical advancements in plant epitranscriptomic studies

The growth and development of plants are the result of the interplay between the internal developmental programming and plant-environment interactions. Gene expression regulations in plants are made up of multi-level networks. In the past few years, many studies were carried out on co- and post-transcriptional RNA modifications, which, together with the RNA community, are collectively known as the "epitranscriptome." The epitranscriptomic machineries were identified and their functional impacts characterized in a broad range of physiological processes in diverse plant species. There is mounting evidence to suggest that the epitranscriptome provides an additional layer in the gene regulatory network for plant development and stress responses. In the present review, we summarized the epitranscriptomic modifications found so far in plants, including chemical modifications, RNA editing, and transcript isoforms. The various approaches to RNA modification detection were described, with special emphasis on the recent development and application potential of third-generation sequencing. The roles of epitranscriptomic changes in gene regulation during plant-environment interactions were discussed in case studies. This review aims to highlight the importance of epitranscriptomics in the study of gene regulatory networks in plants and to encourage multi-omics investigations using the recent technical advancements.

Assembly and analysis of the complete mitochondrial genome of Forsythia suspensa (Thunb.) Vahl

BACKGROUND

Forsythia suspensa (Thunb.) Vahl is a valuable ornamental and medicinal plant. Although the nuclear and chloroplast genomes of F. suspensa have been published, its complete mitochondrial genome sequence has yet to be reported. In this study, the genomic DNA of F. suspensa yellowish leaf material was extracted, sequenced by using a mixture of Illumina Novaseq6000 short reads and Oxford Nanopore PromethION long reads, and the sequencing data were assembled and annotated.

RESULT

The F. suspensa mitochondrial genome was obtained in the length of 535,692 bp with a circular structure, and the GC content was 44.90%. The genome contains 60 genes, including 36 protein-coding genes, 21 tRNA genes, and three rRNA genes. We further analyzed RNA editing of the protein-coding genes, relative synonymous codon usage, and sequence repeats based on the genomic data. There were 25 homologous sequences between F. suspensa mitochondria and chloroplast genome, which involved the transfer of 8 mitochondrial genes, and 9473 homologous sequences between mitochondrial and nuclear genomes. Analysis of the nucleic acid substitution rate, nucleic acid diversity, and collinearity of protein-coding genes of the F. suspensa mitochondrial genome revealed that the majority of genes may have undergone purifying selection, exhibiting a slower rate of evolution and a relatively conserved structure. Analysis of the phylogenetic relationships among different species revealed that F. suspensa was most closely related to Olea europaea subsp. Europaea.

CONCLUSION

In this study, we sequenced, assembled, and annotated a high-quality F. suspensa mitochondrial genome. The results of this study will enrich the mitochondrial genome data of Forsythia, lay a foundation for the phylogenetic development of Forsythia, and promote the evolutionary analysis of Oleaceae species.

The emerging role of epitranscriptome in shaping stress responses in plants

KEY MESSAGE

RNA modifications and editing changes constitute 'epitranscriptome' and are crucial in regulating the development and stress response in plants. Exploration of the epitranscriptome and associated machinery would facilitate the engineering of stress tolerance in crops. RNA editing and modifications post-transcriptionally decorate almost all classes of cellular RNAs, including tRNAs, rRNAs, snRNAs, lncRNAs and mRNAs, with more than 170 known modifications, among which m6A, Ψ, m5C, 8-OHG and C-to-U editing are the most abundant. Together, these modifications constitute the "epitranscriptome", and contribute to changes in several RNA attributes, thus providing an additional structural and functional diversification to the "cellular messages" and adding another layer of gene regulation in organisms, including plants. Numerous evidences suggest that RNA modifications have a widespread impact on plant development as well as in regulating the response of plants to abiotic and biotic stresses. High-throughput sequencing studies demonstrate that the landscapes of m6A, m5C, Am, Cm, C-to-U, U-to-G, and A-to-I editing are remarkably dynamic during stress conditions in plants. GO analysis of transcripts enriched in Ψ, m6A and m5C modifications have identified bonafide components of stress regulatory pathways. Furthermore, significant alterations in the expression pattern of genes encoding writers, readers, and erasers of certain modifications have been documented when plants are grown in challenging environments. Notably, manipulating the expression levels of a few components of RNA editing machinery markedly influenced the stress tolerance in plants. We provide updated information on the current understanding on the contribution of RNA modifications in shaping the stress responses in plants. Unraveling of the epitranscriptome has opened new avenues for designing crops with enhanced productivity and stress resilience in view of global climate change.

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.

The Genome-Wide Characterization of Alternative Splicing and RNA Editing in the Development of Coprinopsis cinerea

Coprinopsis cinerea is one of the model species used in fungal developmental studies. This mushroom-forming Basidiomycetes fungus has several developmental destinies in response to changing environments, with dynamic developmental regulations of the organism. Although the gene expression in C. cinerea development has already been profiled broadly, previous studies have only focused on a specific stage or process of fungal development. A comprehensive perspective across different developmental paths is lacking, and a global view on the dynamic transcriptional regulations in the life cycle and the developmental paths is far from complete. In addition, knowledge on co- and post-transcriptional modifications in this fungus remains rare. In this study, we investigated the transcriptional changes and modifications in C. cinerea during the processes of spore germination, vegetative growth, oidiation, sclerotia formation, and fruiting body formation by inducing different developmental paths of the organism and profiling the transcriptomes using the high-throughput sequencing method. Transition in the identity and abundance of expressed genes drive the physiological and morphological alterations of the organism, including metabolism and multicellularity construction. Moreover, stage- and tissue-specific alternative splicing and RNA editing took place and functioned in C. cinerea. These modifications were negatively correlated to the conservation features of genes and could provide extra plasticity to the transcriptome during fungal development. We suggest that C. cinerea applies different molecular strategies in its developmental regulation, including shifts in expressed gene sets, diversifications of genetic information, and reversible diversifications of RNA molecules. Such features would increase the fungal adaptability in the rapidly changing environment, especially in the transition of developmental programs and the maintenance and balance of genetic and transcriptomic divergence. The multi-layer regulatory network of gene expression serves as the molecular basis of the functioning of developmental regulation.

The first complete mitochondrial genome of Carex (C. breviculmis): a significantly expanded genome with highly structural variations

MAIN CONCLUSION

The first complete mitochondrial genome of Carex (C. breviculmis) was sequenced and assembled, and its genomic signature was analyzed and the possible conformations of its mitochondrial genome were validated. Carex breviculmis is a very adaptable grass that is highly resistant to environmental stresses such as drought and low light. It is also admired as a landscape plant with high development prospects and scientific research value. In this study, the mitochondrial genome of C. breviculmis was assembled using Pacbio and Illumina sequencing data. We detected 267 pairs of repeats and found that three pairs of repeats could mediate the recombination of its mitochondrial genome and formed four possible conformations, of which we verified the two conformations mediated by the shortest pair of repeats using PCR amplification and Sanger sequencing. The major conformation of the C. breviculmis mitochondrial genome is a 1,414,795 bp long circular molecule with 33 annotated protein-coding genes, 15 tRNA genes, and three rRNA genes. We detected a total of 25 homologous sequences between the chloroplast and mitochondrial genomes, corresponding to 0.40% of the mitochondrial genome. Combined with the low GC content (41.24%), we conclude that the reduction in RNA editing sites in the C. breviculmis mitochondrial genome may be due to an accumulation of point mutations in C-to-T or retroprocessing events within the genome. The relatively low number of RNA editing sites in its mitochondrial genome could serve as important material for subsequent studies on the selection pressure of RNA editing in angiosperms. A maximum likelihood analysis based on 23 conserved mitochondrial genes from 28 species reflects an accurate evolutionary and taxonomic position of C. breviculmis. This research provided us with a comprehensive understanding of the mitochondrial genome of Carex and also provided important information for its molecular breeding.

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.

Mitochondrial RNA editing sites affect the phylogenetic reconstruction of gymnosperms

•RNA editing sites may contain homoplasious signals that cause artifactual inferences in phylogenetic analyses.•Excluding RNA editing sites from gymnosperm mitochondrial genes restored the sister relationship of gnetophytes and Pinaceae.•Phylogenetic analysis based on mitochondrial genomic data should carefully evaluate the impact of RNA editing sites.

Comparative analysis of the chloroplast genomes of Rosa species and RNA editing analysis

BACKGROUND

The genus Rosa (Rosaceae) contains approximately 200 species, most of which have high ecological and economic values. Chloroplast genome sequences are important for studying species differentiation, phylogeny, and RNA editing.

RESULTS

In this study, the chloroplast genomes of three Rosa species, Rosa hybrida, Rosa acicularis, and Rosa rubiginosa, were assembled and compared with other reported Rosa chloroplast genomes. To investigate the RNA editing sites in R. hybrida (commercial rose cultivar), we mapped RNA-sequencing data to the chloroplast genome and analyzed their post-transcriptional features. Rosa chloroplast genomes presented a quadripartite structure and had highly conserved gene order and gene content. We identified four mutation hotspots (ycf3-trnS, trnT-trnL, psbE-petL, and ycf1) as candidate molecular markers for differentiation in the Rosa species. Additionally, 22 chloroplast genomic fragments with a total length of 6,192 bp and > 90% sequence similarity with their counterparts were identified in the mitochondrial genome, representing 3.96% of the chloroplast genome. Phylogenetic analysis including all sections and all subgenera revealed that the earliest divergence in the chloroplast phylogeny roughly distinguished species of sections Pimpinellifoliae and Rosa and subgenera Hulthemia. Moreover, DNA- and RNA-sequencing data revealed 19 RNA editing sites, including three synonymous and 16 nonsynonymous, in the chloroplast genome of R. hybrida that were distributed among 13 genes.

CONCLUSIONS

The genome structure and gene content of Rosa chloroplast genomes are similar across various species. Phylogenetic analysis based on the Rosa chloroplast genomes has high resolution. Additionally, a total of 19 RNA editing sites were validated by RNA-Seq mapping in R. hybrida. The results provide valuable information for RNA editing and evolutionary studies of Rosa and a basis for further studies on genomic breeding of Rosa species.

Structural insight into the activation of an Arabidopsis organellar C-to-U RNA editing enzyme by active site complementation

RNA-binding pentatricopeptide repeat (PPR) proteins catalyze hundreds of cytidine to uridine RNA editing events in plant organelles; these editing events are essential for proper gene expression. More than half of the PPR-type RNA editing factors, however, lack the DYW cytidine deaminase domain. Genetic analyses have suggested that their cytidine deaminase activity arises by association with a family of DYW1-like proteins that contain an N-terminally truncated DYW domain, but their molecular mechanism has been unclear. Here, we report the crystal structure of the Arabidopsis thaliana DYW1 deaminase domain at 1.8 Å resolution. DYW1 has a cytidine deaminase fold lacking the PG box. The internal insertion within the deaminase fold shows an α-helical fold instead of the β-finger reported for the gating domain of the A. thaliana ORGANELLE TRANSCRIPT PROCESSING 86. The substrate-binding pocket is incompletely formed and appears to be complemented in the complex by the E2 domain and the PG box of the interacting PPR protein. In vivo RNA editing assays corroborate the activation model for DYW1 deaminase. Our study demonstrates the common activation mechanism of the DYW1-like proteins by molecular complementation of the DYW domain and reconstitution of the substrate-binding pocket.

Plant organellar RNA maturation

Plant organellar RNA metabolism is run by a multitude of nucleus-encoded RNA-binding proteins (RBPs) that control RNA stability, processing, and degradation. In chloroplasts and mitochondria, these post-transcriptional processes are vital for the production of a small number of essential components of the photosynthetic and respiratory machinery-and consequently for organellar biogenesis and plant survival. Many organellar RBPs have been functionally assigned to individual steps in RNA maturation, often specific to selected transcripts. While the catalog of factors identified is ever-growing, our knowledge of how they achieve their functions mechanistically is far from complete. This review summarizes the current knowledge of plant organellar RNA metabolism taking an RBP-centric approach and focusing on mechanistic aspects of RBP functions and the kinetics of the processes they are involved in.

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 complete chloroplast genome sequence of Eugenia klotzschiana O. Berg unveils the evolutionary dynamics in plastomes of Myrteae DC. tribe (Myrtaceae)

Myrteae is the most diversified tribe in the Myrtaceae family and has great ecological and economic importance. Here, we performed the assembly and annotation of the chloroplast genome of Eugenia klotzschiana O. Berg and used this in a comparative analysis with other 13 species from the Myrteae tribe. The E. klotzschiana plastome exhibited a length of 158,977 bp and a very conserved structure and gene composition when compared with other Myrteae genomes. We identified 34 large repetitive sequences and 94 SSR repeats in E. klotzschiana plastome. The trnT-trnL, rpl32-trnL, ndhF-rpl32, psbE-petL, and ycf1 regions were identified as mutational hotspots. A negative selection signal was detected in 74 protein-coding genes while neutral selection was detected in two genes (rps12 and psaI). Furthermore, 222 RNA editing sites were identified in the E. klotzschiana plastome. We also obtained a plastome-based Myrtales phylogenetic tree, including E. klotzschiana for the first time in a molecular phylogeny, recovering its sister relationship for all other Eugenia species. Our results illuminate how evolution shaped the chloroplast genome structure and composition in the Myrteae tribe, especially in the E. klotzschiana plastome.

Characterisation of the complete mitochondrial genome of Taraxacum mongolicum revealed five repeat-mediated recombinations

We reported the mitochondrial genome of Taraxacum mongolicum for the first time. Five pairs of repeats that can mediate recombination were validated, leading to multiple conformations of genome. Taraxacum mongolicum belongs to the Asteraceae family and has important pharmaceutical value. To explore the possible interaction between the organelle genomes, we assembled the complete mitochondrial genome (mitogenome) of T. mongolicum using Illumina and Oxford Nanopore sequencing data. This genome corresponded to a circular molecule 304,467 bp long. It encodes 52 unique genes including 31 protein-coding, 3 ribosomal RNA (rRNA) and 18 transfer RNA (tRNA) genes. In addition to the single circular conformation, the existence of alternative conformations mediated by five repetitive sequences in the mitogenome was identified and validated. Recombination mediated by the inverted repeats resulted in two conformations. Conversely, recombination mediated by the two direct repeats broke one large circular molecule into two subgenomic circular molecules. Furthermore, we identified 12 homologous fragments by comparing the sequences of mitogenome and plastome, including eight complete tRNA genes. Lastly, we identified a total of 278 RNA-editing sites in protein-coding sequences based on RNA-seq data. Among them, cox1 and nad5 gene has the most sites (21), followed by the nad2 gene with 19 sites. We successfully validated 213 predicted RNA-editing sites using PCR amplification and Sanger sequencing. This project reported the first mitogenome of T. mongolicum and demonstrated its multiple conformations generated by repeat-mediated recombination. This genome could provide critical information for the molecular breeding of T. mongolicum, and also be used as a reference genome for other species of the genus Taraxacum.

Poaceae Chloroplast Genome Sequencing: Great Leap Forward in Recent Ten Years

The first complete chloroplast genome of rice (Oryza sativa) was published in 1989, ushering in a new era of studies of chloroplast genomics in Poaceae. Progresses in Next-Generation Sequencing (NGS) and Third-Generation Sequencing (TGS) technologiesand in the development of genome assembly software, have significantly advanced chloroplast genomics research. Poaceae is one of the most targeted families in chloroplast genome research because of its agricultural, ecological, and economic importance. Over the last 30 years, 2,050 complete chloroplast genome sequences from 40 tribes and 282 genera have been generated, most (97%) of them in the recent ten years. The wealth of data provides the groundwork for studies on species evolution, phylogeny, genetic transformation, and other aspects of Poaceae chloroplast genomes. As a result, we have gained a deeper understanding of the properties of Poaceae chloroplast genomes. Here, we summarize the achievements of the studies of the Poaceae chloroplast genomes and envision the challenges for moving the area ahead.

SlRIP1b is a global organellar RNA editing factor, required for normal fruit development in tomato plants

RNA editing in plant organelles involves numerous C-U conversions, which often restore evolutionarily conserved codons and may generate new translation initiation and termination codons. These RNA maturation events rely on a subset of nuclear-encoded protein cofactors. Here, we provide evidence of the role of SlRIP1b on RNA editing of mitochondrial transcripts in tomato (Solanum lycopersicum) plants. SlRIP1b is a RIP/MORF protein that was originally identified as an interacting partner of the organellar editing factor SlORRM4. Mutants of SlRIP1b, obtained by CRISPR/Cas9 strategy, exhibited abnormal carpel development and grew into fruit with more locules. RNA-sequencing revealed that SlRIP1b affects the C-U editing of numerous mitochondrial pre-RNA transcripts and in particular altered RNA editing of various cytochrome c maturation (CCM)-related genes. The slrip1b mutants display increased H₂ O₂ and aberrant mitochondrial morphologies, which are associated with defects in cytochrome c biosynthesis and assembly of respiratory complex III. Taken together, our results indicate that SlRIP1b is a global editing factor that plays a key role in CCM and oxidative phosphorylation system biogenesis during fruit development in tomato plants. These data provide important insights into the molecular roles of organellar RNA editing factors during fruit development.

Effect of salinity on ccmfn gene RNA editing of mitochondria in wild barley and uncommon types of RNA editing

The primary function of mitochondria is cellular respiration and energy production. Cytochrome C complex is an essential complex that transports electrons in the respiratory chain between complex III and complex IV. One of this complex's main subunits is CcmFN, which is believed to be crucial for holocytochrome assembly. In wild-type plant Hordeum vulgare subsp. spontaneum, four ccmfn cDNAs are subjected to high salt stress (500 mM salinity), 0 h (or control) (GenBank accession no. ON764850), after 2 h (GenBank accession no. ON7648515), after 12 h (GenBank accession no. ON764852), and after 24 h (GenBank accession no. ON764853) and mtDNA of ccmfn gene (GenBank accession no. ON764854). Using raw data from RNA-seq, 47 sites with nucleotide and amino acid modifications were detected. There were ten different RNA editing types, with most of them are C to U. Unusual editing types in plants have also been found, such as A to C, C to A, A to G, A to U, T to A, T to C, C to G, G to C, and T to G. High levels of editing were observed in control as well as treatments of salinity stress. Amino acid changes were found in 43 sites; nearly all showed hydrophilic to hydrophilic alterations. Only C749 showed regulation under salinity stress.

Lactobacillus for ribosome peptide editing cancer

This study reviews newly discovered insect peptide point mutations as new possible cancer research targets. To interpret newly discovered peptide point mutations in insects as new possible cancer research targets, we focused on the numerous peptide changes found in the 'CSP' family on the sex pheromone gland of the female silkworm moth Bombyx mori. We predict that the Bombyx peptide modifications will have a significant effect on cancer CUP (cancers of unknown primary) therapy and that bacterial peptide editing techniques, specifically Lactobacillus combined to CRISPR, will be used to regulate ribosomes and treat cancer in humans.

Characterizing microRNA editing and mutation sites in Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a neurodevelopmental disorder whose pathogenesis is still unclear. MicroRNAs (miRNAs) are a kind of endogenous small non-coding RNAs that play important roles in the post-transcriptional regulation of genes. Recent researches show that miRNAs are edited in multiple ways especially in central nervous systems. A-to-I editing of RNA catalyzed by Adenosine deaminases acting on RNA (ADARs) happens intensively in brain and is also noticed in other organs and tissues. Although miRNAs are widely edited in human brain, miRNA editing in ASD is still largely unexplored. In order to reveal the editing events of miRNAs in ASD, we analyzed 131 miRNA-seq samples from 8 different brain regions of ASD patients and normal controls. We identified 834 editing sites with significant editing levels, of which 70 sites showed significantly different editing levels in the superior frontal gyrus samples of ASD patients (ASD-SFG) when compared with those of control samples. The editing level of an A-to-I editing site in hsa-mir-376a-1 (hsa-mir-376a-1_9_A_g) in ASD-SFG is higher than that of normal controls, and the difference is exaggerated in individuals under 10 years. The increased expression of ADAR1 is consistent with the increased editing level of hsa-mir-376a-1_9_A_g in ASD-SFG samples compared to normal SFG samples. Furthermore, we verify that A-to-I edited hsa-mir-376a-5p directly represses GPR85 and NAPB, which may contribute to the abnormal neuronal development of ASD patients. These results provide new insights into the mechanism of ASD.

Applications and Prospects of CRISPR/Cas9-Mediated Base Editing in Plant Breeding

The clustered regularly interspaced short palindromic repeats (CRISPR)/associated protein 9 system (Cas9) has been used at length to optimize multiple aspects of germplasm resources. However, large-scale genomic research has indicated that novel variations in crop plants are attributed to single-nucleotide polymorphisms (SNPs). Therefore, substituting single bases into a plant genome may produce desirable traits. Gene editing by CRISPR/Cas9 techniques frequently results in insertions-deletions (indels). Base editing allows precise single-nucleotide changes in the genome in the absence of double-strand breaks (DSBs) and donor repair templates (DRTs). Therefore, BEs have provided a new way of thinking about genome editing, and base editing techniques are currently being utilized to edit the genomes of many different organisms. As traditional breeding techniques and modern molecular breeding technologies complement each other, various genome editing technologies have emerged. How to realize the greater potential of BE applications is the question we need to consider. Here, we explain various base editings such as CBEs, ABEs, and CGBEs. In addition, the latest applications of base editing technologies in agriculture are summarized, including crop yield, quality, disease, and herbicide resistance. Finally, the challenges and future prospects of base editing technologies are presented. The aim is to provide a comprehensive overview of the application of BE in crop breeding to further improve BE and make the most of its value.

Characterization and phylogenetic analysis of the complete mitochondrial genome sequence of Photinia serratifolia

Plant mitochondrial genomes (mitogenomes) are a valuable source of genetic information for a better understanding of phylogenetic relationships. However, no mitogenome of any species in the genus of Photinia has been reported. In this study, using NGS sequencing, we reported the mitogenome assembly and annotation of Photinia serratifolia, which is 473,579 bp in length, contains 38 protein-coding genes, 23 tRNAs, and 6 rRNAs, with 61 genes have no introns. The rps2 and rps11 genes are missing in the P. serratifolia mitogenome. Although there are more editing sites (488) in the P. serratifolia mitogenome than in most angiosperms, fewer editing types were found in the P. serratifolia mitogenome, showing a clear bias in RNA-editing. Phylogenetic analysis based on the mitogenomes of P. serratifolia and 8 other taxa of the Rosaceae family reflected the exact evolutionary and taxonomic status of P. serratifolia. However, Ka/Ks analysis revealed that 72.69% of the protein-coding genes in the P. serratifolia mitogenome had undergone negative selections, reflecting the importance of those genes in the P. serratifolia mitogenome. Collectively, these results will provide valuable information for the evolution of P. serratifolia and provide insight into the evolutionary relationships within Photinia and the Rosaceae family.

Maize cytolines as models to study the impact of different cytoplasms on gene expression under heat stress conditions

BACKGROUND

Crops are under constant pressure due to global warming, which unfolds at a much faster pace than their ability to adapt through evolution. Agronomic traits are linked to cytoplasmic-nuclear genome interactions. It thus becomes important to understand the influence exerted by the organelles on gene expression under heat stress conditions and profit from the available genetic diversity. Maize (Zea mays) cytolines allow us to investigate how the gene expression changes under heat stress conditions in three different cytoplasmic environments, but each having the same nucleus. Analyzing retrograde signaling in such an experimental set-up has never been done before. Here, we quantified the response of three cytolines to heat stress as differentially expressed genes (DEGs), and studied gene expression patterns in the context of existing polymorphism in their organellar genomes.

RESULTS

Our study unveils a plethora of new genes and GO terms that are differentially expressed or enriched, respectively, in response to heat stress. We report 19,600 DEGs as responding to heat stress (out of 30,331 analyzed), which significantly enrich 164 GO biological processes, 30 GO molecular functions, and 83 GO cell components. Our approach allowed for the discovery of a significant number of DEGs and GO terms that are not common in the three cytolines and could therefore be linked to retrograde signaling. Filtering for DEGs with a fold regulation > 2 (absolute values) that are exclusive to just one of the cytolines, we find a total of 391 up- and down-DEGs. Similarly, there are 19 GO terms with a fold enrichment > 2 that are cytoline-specific. Using GBS data we report contrasting differences in the number of DEGs and GO terms in each cytoline, which correlate with the genetic distances between the mitochondrial genomes (but not chloroplast) and the original nuclei of the cytolines, respectively.

CONCLUSIONS

The experimental design used here adds a new facet to the paradigm used to explain how gene expression changes in response to heat stress, capturing the influence exerted by different organelles upon one nucleus rather than investigating the response of several nuclei in their innate cytoplasmic environments.

A Bibliometric Study for Plant RNA Editing Research: Trends and Future Challenges

RNA editing is a post-transcriptional process that introduces changes in RNA sequences encoded by nuclear, mitochondrial, or plastid genomes. To understand the research progress of plant RNA editing, we comprehensively analyze the articles on plant RNA editing from 2001 to 2022 through bibliometric methods. Nucleic Acids Research, Plant Journal and Plant cell are the journals that deserve attention with their high production, total local citation scores (TLCS), and h-indexes. The USA, China, and Germany are the top three countries with highly productive publications. Ulm University, Cornell University, and Chinese Acad Sci are excellent cooperative institutions with a high level of influence in the field, and KNOOP V and TAKENAKA M are good partnership. Plant RNA editing researches concentrate on the subject categories of Biochemistry & Molecular Biology, Plant Sciences, Genetics & Heredity, etc. Plant mitochondria, genome editing and messenger-RNA may be the research hotspots in the future. The main plant RNA editing research tools are JACUSA, SPRINT, and REDO, and the main databases are REDIdb, PED, and dbRES. At present, the research streams are (1) RNA editing sites; (2) Pentapeptide repeat protein (PPR) involved in RNA editing; (3) RNA editing factors. Overall, this article summarizes the research overview of plant RNA editing until 2022 and provides theoretical implications for its possible future directions.

The Mitogenome of Sedum plumbizincicola (Crassulaceae): Insights into RNA Editing, Lateral Gene Transfer, and Phylogenetic Implications

As the largest family within the order Saxifragales, Crassulaceae contains about 34 genera with 1400 species. Mitochondria play a critical role in cellular energy production. Since the first land plant mitogenome was reported in Arabidopsis, more than 400 mitogenomic sequences have been deposited in a public database. However, no entire mitogenome data have been available for species of Crassulaceae to date. To better understand the evolutionary history of the organelles of Crassulaceae, we sequenced and performed comprehensive analyses on the mitogenome of Sedum plumbizincicola. The master mitogenomic circle is 212,159 bp in length, including 31 protein-coding genes (PCGs), 14 tRNA genes, and 3 rRNA genes. We further identified totally 508 RNA editing sites in PCGs, and demonstrated that the second codon positions of mitochondrial genes are most prone to RNA editing events. Notably, by neutrality plot analyses, we observed that the mitochondrial RNA editing events have large effects on the driving forces of plant evolution. Additionally, 4 MTPTs and 686 NUMTs were detected in the mitochondrial and nuclear genomes of S. plumbizincicola, respectively. Additionally, we conducted further analyses on gene transfer, secondary structures of mitochondrial RNAs, and phylogenetic implications. Therefore, the findings presented here will be helpful for future investigations on plant mitogenomes.

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.

Sequential action of a tRNA base editor in conversion of cytidine to pseudouridine

Post-transcriptional RNA editing modulates gene expression in a condition-dependent fashion. We recently discovered C-to-Ψ editing in Vibrio cholerae tRNA. Here, we characterize the biogenesis, regulation, and functions of this previously undescribed RNA editing process. We show that an enzyme, TrcP, mediates the editing of C-to-U followed by the conversion of U to Ψ, consecutively. AlphaFold-2 predicts that TrcP consists of two globular domains (cytidine deaminase and pseudouridylase) and a long helical domain. The latter domain tethers tRNA substrates during both the C-to-U editing and pseudouridylation, likely enabling a substrate channeling mechanism for efficient catalysis all the way to the terminal product. C-to-Ψ editing both requires and suppresses other modifications, creating an interdependent network of modifications in the tRNA anticodon loop that facilitates coupling of tRNA modification states to iron availability. Our findings provide mechanistic insights into an RNA editing process that likely promotes environmental adaptation.

Identification and expression analysis of the grape pentatricopeptide repeat (PPR) gene family in abiotic stress

The pentatricopeptide repeat (PPR) is one of the largest gene family in plants, and play important role in regulating plant growth, development and abiotic stress response. However, PPR genes have been poorly studied in grapes. In this study, based on the grape genome database, bioinformatics methods and quantitative real-time PCR (qRT-PCR) were used to identify the VvPPR family and the response to abiotic stress. A total of 181 PPR genes were identified in grape and divided into two subfamilies. Subcellular localization predicted that this gene family mainly functions in chloroplasts, nucleus, and mitochondria. Protein-protein interaction prediction indicated that there may be interaction between VvPPR44,53 and VvPPR44. The promoter region of VvPPR gene family contained various cis-acting elements, which were related to light and hormone. Expression pattern analysis showed that the VvPPR gene family was highly expressed in grape leaves, buds and carpel tissues. qRT-PCR results showed that the expression of VvPPR genes in roots was higher than stems and leaves under NAA, SA, ABA, MeJA and GA3 treatments. VvPPR8 was significantly up-regulated after GA3 and MeJA treatment for 24 h, VvPPR53 was significantly up-regulated after SA, NAA, ABA and MeJA treatment. In addition, In grape leaves, VvPPR53 was up-regulated under PEG, Nacl and 4 ℃ treatments. These data indicate that VvPPR gene family members are responsive to hormones and abiotic stresses, and that there are some differences in the degree of response and expression in different grape tissues. This study provides a certain theoretical basis for grape resistance breeding.

The complete mitochondrial genome of Isochrysis galbana harbors a unique repeat structure and a specific trans-spliced cox1 gene

The haptophyte Isochrysis galbana is considered as a promising source for food supplements due to its rich fucoxanthin and polyunsaturated fatty acids content. Here, the I. galbana mitochondrial genome (mitogenome) was sequenced using a combination of Illumina and PacBio sequencing platforms. This 39,258 bp circular mitogenome has a total of 46 genes, including 20 protein-coding genes, 24 tRNA genes and two rRNA genes. A large block of repeats (~12.7 kb) was segregated in one region of the mitogenome, accounting for almost one third of the total size. A trans-spliced gene cox1 was first identified in I. galbana mitogenome and was verified by RNA-seq and DNA-seq data. The massive expansion of tandem repeat size and cis- to trans-splicing shift could be explained by the high mitogenome rearrangement rates in haptophytes. Strict SNP calling based on deep transcriptome sequencing data suggested the lack of RNA editing in both organelles in this species, consistent with previous studies in other algal lineages. To gain insight into haptophyte mitogenome evolution, a comparative analysis of mitogenomes within haptophytes and among eight main algal lineages was performed. A core gene set of 15 energy and metabolism genes is present in haptophyte mitogenomes, consisting of 1 cob, 3 cox, 7 nad, 2 atp and 2 ribosomal genes. Gene content and order was poorly conserved in this lineage. Haptophyte mitogenomes have lost many functional genes found in many other eukaryotes including rps/rpl, sdh, tat, secY genes, which make it contain the smallest gene set among all algal taxa. All these implied the rapid-evolving and more recently evolved mitogenomes of haptophytes compared to other algal lineages. The phylogenetic tree constructed by cox1 genes of 204 algal mitogenomes yielded well-resolved internal relationships, providing new evidence for red-lineages that contained plastids of red algal secondary endosymbiotic origin. This newly assembled mitogenome will add to our knowledge of general trends in algal mitogenome evolution within haptophytes and among different algal taxa.

Assembly and comparative analysis of the complete mitochondrial genome of Bupleurum chinense DC

BACKGROUND

Bupleurum chinense(B. chinense) is a plant that is widely distributed globally and has strong pharmacological effects. Though the chloroplast(cp) genome of B. chinense has been studied, no reports regarding the mitochondrial(mt) genome of B. chinense have been published yet.

RESULTS

The mt genome of B.chinense was assembled and functionally annotated. The circular mt genome of B. chinense was 435,023 bp in length, and 78 genes, including 39 protein-coding genes, 35 tRNA genes, and 4 rRNA genes, were annotated. Repeat sequences were analyzed and sites at which RNA editing would occur were predicted. Gene migration was observed to occur between the mt and cp genomes of B. chinense via the detection of homologous gene fragments. In addition, the sizes of plant mt genomes and their GC content were analyzed and compared. The sizes of mt genomes of plants varied greatly, but their GC content was conserved to a greater extent during evolution. Ka/Ks analysis was based on code substitutions, and the results showed that most of the coding genes were negatively selected. This indicates that mt genes were conserved during evolution.

CONCLUSION

In this study, we assembled and annotated the mt genome of the medicinal plant B. chinense. Our findings provide extensive information regarding the mt genome of B. chinense, and help lay the foundation for future studies on the genetic variations, phylogeny, and breeding of B. chinense via an analysis of the mt genome.

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

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

A chloroplast-localized pentatricopeptide repeat protein involved in RNA editing and splicing and its effects on chloroplast development in rice

BACKGROUND

The chloroplast is the organelle responsible for photosynthesis in higher plants. The generation of functional chloroplasts depends on the precise coordination of gene expression in the nucleus and chloroplasts and is essential for the development of plants. However, little is known about nuclear-plastid regulatory mechanisms at the early stage of chloroplast generation in rice.

RESULTS

In this study, we identified a rice (Oryza sativa) mutant that exhibited albino and seedling-lethal phenotypes and named it ssa1(seedling stage albino1). Transmission electron microscopy (TEM) analysis indicated that the chloroplasts of ssa1 did not have organized thylakoid lamellae and that the chloroplast structure was destroyed. Genetic analysis revealed that the albino phenotypes of ssa1 were controlled by a pair of recessive nuclear genes. Map-based cloning experiments found that SSA1 encoded a pentapeptide repeat (PPR) protein that was allelic to OSOTP51,which was previously reported to participate in Photosystem I (PSI) assembly. The albino phenotype was reversed to the wild type (WT) phenotype when the normal SSA1 sequence was expressed in ssa1 under the drive of the actin promoter. Knockout experiments further created mutants ssa1-2/1-9, which had a phenotype similar to that of ssa1. SSA1 consisted of 7 pentatricopeptide repeat domains and two C-terminal LAGLIDADG tandem sequence motifs and was located in the chloroplast. GUS staining and qRT-PCR analysis showed that SSA1 was mainly expressed in young leaves and stems. In the ssa1 mutants, plastid genes transcribed by plastid-encoded RNA polymerase decreased, while those transcribed by nuclear-encoded RNA polymerase increased at the mRNA level. Loss-of-function SSA1 destroys RNA editing of ndhB-737 and intron splicing of atpF and ycf3-2 in the plastid genome. Yeast two-hybrid and BiFC assays revealed that SSA1 physically interacted with two new RNA editing partners, OsMORF8 and OsTRXz, which have potential functions in RNA editing and chloroplast biogenesis.

CONCLUSIONS

Rice SSA1 encodes a pentatricopeptide repeat protein, which is targeted to the chloroplast. SSA1 regulates early chloroplast development and plays a critical role in RNA editing and intron splicing in rice. These data will facilitate efforts to further elucidate the molecular mechanism of chloroplast biogenesis.

C-to-U RNA Editing: A Site Directed RNA Editing Tool for Restoration of Genetic Code

The restoration of genetic code by editing mutated genes is a potential method for the treatment of genetic diseases/disorders. Genetic disorders are caused by the point mutations of thymine (T) to cytidine (C) or guanosine (G) to adenine (A), for which gene editing (editing of mutated genes) is a promising therapeutic technique. In C-to-Uridine (U) RNA editing, it converts the base C-to-U in RNA molecules and leads to nonsynonymous changes when occurring in coding regions; however, for G-to-A mutations, A-to-I editing occurs. Editing of C-to-U is not as physiologically common as that of A-to-I editing. Although hundreds to thousands of coding sites have been found to be C-to-U edited or editable in humans, the biological significance of this phenomenon remains elusive. In this review, we have tried to provide detailed information on physiological and artificial approaches for C-to-U RNA editing.