RNA editing

Phylogenetic and structural analysis of Hydra ADAR

Adenosine deaminases acting on RNAs (ADARs) perform adenosine-to-inosine (A-to-I) RNA editing for essential biological functions. While studies of editing sites in diverse animals have revealed unique biological roles of ADAR editing including temperature adaptation and reproductive maturation, rigorous biochemical and structural studies of these ADARs are lacking. Here, we present a phylogenetic sequence analysis and AlphaFold computational structure prediction to reveal that medusozoan ADAR2s contain five dsRNA binding domains (dsRBDs) with several RNA binding residues in the dsRBDs and deaminase domain conserved. Additionally, we identified evolutionary divergence between the medusozoan (e.g. Hydra) and anthozoan cnidarian subphyla. The anthozoan ADAR deaminase domains more closely resemble human ADARs with longer 5' RNA binding loops, glutamate base-flipping residues, and a conserved TWDG dimerization motif. Conversely, medusozoan ADAR deaminase domains have short 5' binding loops, glutamine flipping residues, and non-conserved helix dimerization motif. We also report the direct detection of A-to-I RNA editing by an ADAR ortholog from the freshwater cnidarian Hydra vulgaris (hyADAR). We solved the crystal structure of the monomeric deaminase domain of hyADAR (hyADARd) to 2.0 Å resolution, showing conserved active site architecture and the presence of a buried inositol hexakisphosphate known to be required for ADAR activity. In addition, these data demonstrate that medusozoans have evolved novel ADAR structural features, however the physiological consequence of this remains unknown. In addition, these results provide a framework for biochemically and structurally characterizing ADARs from evolutionarily distant organisms to understand the diverse roles of ADAR editing amongst metazoans.

Assembly and comparative analysis of the complete mitochondrial genome of Daedaleopsissinensis (Polyporaceae, Basidiomycota), contributing to understanding fungal evolution and ecological functions

Daedaleopsissinensis is a crucial wood-decaying fungus with significant lignocellulose-degrading ability, which plays a vital role in the material cycle and energy flow of forest ecosystems. However, the mitochondrial genome of D.sinensis has not yet been revealed. In the present study, the complete mitochondrial genome of D.sinensis was assembled and compared with related species. The mitochondrial genome spans 69,155 bp and has a GC content of 25.0%. It comprises 15 protein-coding genes (PCGs), 26 transfer RNA genes, two ribosomal RNA genes and one DNA polymerase gene (dpo). Herein, we characterised and analysed the codon preferences, variation and evolution of PCGs, repeats, intron dynamics, as well as RNA editing events in the D.sinensis mitochondrial genome. Further, a phylogenetic analysis of D.sinensis and the other 86 Basidiomycota species was performed using mitochondrial genome data. The results revealed that four species, D.confragosa, D.sinensis, D.nitida and Fomesfomentarius, were grouped in a closely-related cluster with high support values, indicating that a close phylogenetic relationship existed between Daedaleopsis and Fomes. This study reported on the initial assembly and annotation of the mitochondrial genome of D.sinensis, which greatly improved the knowledge of the fungus. These results contribute to the limited understanding of the mitochondrial repository of wood-decaying fungi, thereby laying the foundation for subsequent research on fungal evolution and ecological functions.

mRNA Transcript Variants Expressed in Mammalian Cells

Gene expression or gene regulation studies often assume one gene expresses one mRNA. However, contrary to the conventional idea, a single gene in mammalian cells can express multiple transcript variants translated into several different proteins. The transcript variants are generated through transcription from alternative start sites and alternative post-transcriptional processing of the precursor mRNA (pre-mRNA). In addition, gene mutations and RNA editing further enhance the diversity of the transcript variants. The transcript variants can encode proteins with various domains, expanding the functional repertoire of a single gene. Some transcript variants may not encode proteins but function as non-coding RNAs and regulate gene expression. The expression level of the transcript variants may vary between cell types or within the same cells under different biological conditions. Transcript variants are characteristic of cell differentiation in a particular tissue, and the variants may play a key role in normal development and aging. Studies also reported that some transcript variants may have roles in disease pathogenesis. The biological significances urge studying the complexity of gene expression at the transcript level. This article updates the molecular basis of transcript variants in mammalian cells, including the formation mechanisms and potential roles in host biology. Gaining insight into the transcript variants will not only identify novel mechanisms of gene regulation but also unravel the role of the variants in health and disease.

Restoration of cytidine to uridine genetic code using an MS2-APOBEC1 artificial enzymatic approach

By employing site-directed RNA editing (SDRE) to restore point-mutated RNA molecules, it is possible to change gene-encoded information and synthesize proteins with different functionality from a single gene. Thymine (T) to cytosine (C) point mutations cause various genetic disorders, and when they occur in protein-coding regions, C-to-uridine (U) RNA changes can lead to non-synonymous alterations. By joining the deaminase domain of apolipoprotein B messenger RNA (mRNA) editing catalytic polypeptide 1 (APOBEC1) with a guide RNA (gRNA) complementary to a target mRNA, we created an artificial RNA editase. We used an mRNA encoding blue fluorescent protein (BFP), obtained from the green fluorescent protein (GFP) gene through the introduction of a T > C mutation, as our target RNA. In a proof of principle experiment, we reverted the T > C mutation at the RNA level using our APOBEC1 site-directed RNA editing system, recovering GFP signal. Sanger sequencing of cDNA from transfected cells and polymerase chain reaction-restriction length polymorphism analysis validated this result, indicating an editing of approximately 21 %. Our successful development of an artificial RNA editing system using the deaminase APOBEC1, in conjunction with the MS2 system, may lead to the development of treatments for genetic diseases based on the restoration of specific types of wild type sequences at the mRNA level.

REMR: Identification of RNA Editing-mediated MiRNA Regulation in Cancers

Dysregulation of adenosine-to-inosine (A-to-I) RNA editing has been implicated in cancer progression. However, a comprehensive understanding of how A-to-I RNA editing is incorporated into miRNA regulation to modulate gene expression in cancer remains unclear, given the lack of effective identification methods. To this end, we introduced an information theory-based algorithm named REMR to systematically identify 12,006 A-to-I RNA editing-mediated miRNA regulatory triplets (RNA editing sites, miRNAs, and genes) across ten major cancer types based on multi-omics profiling data from The Cancer Genome Atlas (TCGA). Through analyses of functional enrichment, transcriptional regulatory networks, and protein-protein interaction (PPI) networks, we showed that RNA editing-mediated miRNA regulation potentially affects critical cancer-related functions, such as apoptosis, cell cycle, drug resistance, and immunity. Furthermore, triplets can serve as biomarkers for classifying cancer subtypes with distinct prognoses or drug responses, highlighting the clinical relevance of such regulation. In addition, an online resource (http://www.jianglab.cn/REMR/) was constructed to support the convenient retrieval of our findings. In summary, our study systematically dissected the RNA editing-mediated miRNA regulations, thereby providing a valuable resource for understanding the mechanism of RNA editing as an epitranscriptomic regulator in cancer.

The Biological Mechanisms and Clinical Roles of RNA-Binding Proteins in Cardiovascular Diseases

RNA-binding proteins (RBPs) have pivotal roles in cardiovascular biology, influencing various molecular mechanisms underlying cardiovascular diseases (CVDs). This review explores the significant roles of RBPs, focusing on their regulation of RNA alternative splicing, polyadenylation, and RNA editing, and their impact on CVD pathogenesis. For instance, RBPs are crucial in myocardial injury, contributing to disease progression and repair mechanisms. This review systematically analyzes the roles of RBPs in myocardial injury, arrhythmias, myocardial infarction, and heart failure, revealing intricate interactions that influence disease outcomes. Furthermore, the potential of RBPs as therapeutic targets for cardiovascular dysfunction is explored, highlighting the advances in drug development and clinical research. This review also discusses the emerging role of RBPs as biomarkers for cardiovascular diseases, offering insights into their diagnostic and prognostic potential. Despite significant progress, current research faces several limitations, which are critically examined. Finally, this review identifies the major challenges and outlines future research directions to advance the understanding and application of RBPs in cardiovascular medicine.

Unraveling the complex evolutionary features of the Cinnamomum camphora mitochondrial genome

KEY MESSAGE

We reported the mitochondrial genome of Cinnamomum camphora for the first time, revealing frequent rearrangement events in the non-coding regions of Magnoliids mitochondrial genomes. As one of the representative species in the Lauraceae family of Magnoliids, Cinnamomum camphora holds significant economic and ecological value. In this study, the mitochondrial genome (mitogenome) of C. camphora was complete assembled and annotated using PacBio HiFi sequencing. The C. camphora mitogenome is characterized by a branch structure, spans 900,894 bp, and contains 43 protein-coding genes (PCGs), 24 tRNAs, and 3 rRNAs. Most of these PCGs are under purifying selection, with only two (ccmFc and rps7) exhibiting signs of positive selection. The C. camphora mitogenome contains numerous repetitive sequences and intracellular gene transfers, with a total of 36 mitochondrial plastid DNAs, amounting to a combined length of 23,816 bp. Comparative analysis revealed that the non-coding regions of Magnoliids mitogenomes have undergone frequent rearrangements during evolution, but the coding sequences remain highly conserved (more than 98% similarity for protein-coding sequences). Furthermore, a maximum-likelihood phylogenetic tree was reconstructed based on 25 PCGs from 23 plant mitogenomes. The analysis supports the closest relationship between C. camphora and C. chekiangense, consistent with the APG IV classification system. This study elucidates the unique evolutionary features of the C. camphora mitogenome, which will provide valuable insights into the study of genetics and evolution of the family Lauraceae.

RNA editing landscape of adipose tissue in polycystic ovary syndrome provides insight into the obesity-related immune responses

Background

Polycystic ovary syndrome (PCOS) is the most common reproductive-endocrine disorder with wide-ranging metabolic implications, including obesity. RNA editing, a post-transcriptional modification, can fine-tune protein function and introduce heterogeneity. However, the role of RNA editing and its impact on adipose tissue function in PCOS remain poorly understood.

Methods

This study aimed to comprehensively analyze RNA-editing events in abdominal and subcutaneous adipose tissue of PCOS patients and healthy controls using high-throughput whole-genome sequencing (WGS) and RNA sequencing.

Results

Our results revealed that PCOS patients exhibited more RNA-editing sites, with adenosine-to-inosine (A-to-I) editing being prevalent. The expression of ADAR genes, responsible for A-to-I editing, was also higher in PCOS. Aberrant RNA-editing sites in PCOS adipose tissue was enriched in immune responses, and interleukin-12 biosynthetic process. Tumor necrosis factor (TNF) signaling, nuclear factor kappa B (NF-κB) signaling, Notch signaling, terminal uridylyl transferase 4 (TUT4), hook microtubule tethering protein 3 (HOOK3), and forkhead box O1 (FOXO1) were identified to be of significant differences. Differentially expressed genes (DEGs) in PCOS adipose tissue were enriched in immune responses compared with controls, and the DEGs between subcutaneous and abdominal adipose tissue were also enriched in immune responses suggesting the important role of subcutaneous adipose tissue. Furthermore, we identified the correlations between RNA editing levels and RNA expression levels of specific genes, such as ataxia-telangiectasia mutated (ATM) and mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1) in inflammation pathways and ATM, TUT4, and YTH N6-methyladenosine RNA-binding protein C2 (YTHDC2) in oocyte development pathway.

Conclusions

These findings suggest that RNA-editing dysregulation in PCOS adipose tissue may contribute to inflammatory dysregulations. Understanding the interplay between RNA editing and adipose tissue function may unveil potential therapeutic targets for PCOS management. However, further research and validation are required to fully elucidate the molecular mechanisms underlying these associations.

Epitranscriptomic Regulations in the Heart

RNA modifications affect key stages of the RNA life cycle, including splicing, export, decay, and translation. Epitranscriptomic regulations therefore significantly influence cellular physiology and pathophysiology. Here, we selected some of the most abundant modifications and reviewed their roles in the heart and in cardiovascular diseases: N6-methyladenosine (m6A), N6,2'-O-dimethyladenosine (m6Am), N1-methyladenosine (m1A), pseudouridine (?), 5 methylcytidine (m5C), and inosine (I). Dysregulation of epitranscriptomic machinery affecting these modifications vastly changes the cardiac phenotype and is linked with many cardiovascular diseases such as myocardial infarction, cardiomyopathies, or heart failure. Thus, a deeper understanding of these epitranscriptomic changes and their regulatory mechanisms can enhance our knowledge of the molecular underpinnings of prevalent cardiac diseases, potentially paving the way for novel therapeutic strategies. Keywords: Epitranscriptomics, RNA modifications, Epigenetics, m6A, RNA, Heart.

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

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

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 Complete Mitochondrial Genome of Paeonia lactiflora Pall. (Saxifragales: Paeoniaceae): Evidence of Gene Transfer from Chloroplast to Mitochondrial Genome

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

RNA Editing as a General Trait of Ebolaviruses

RNA editing has been discovered as an essential mechanism for the transcription of the glycoprotein (GP) gene of Ebola virus but not Marburg virus. We developed a rapid transcript quantification assay (RTQA) to analyze RNA transcripts generated through RNA editing and used immunoblotting with a pan-ebolavirus monoclonal antibody to confirm different GP gene-derived products. RTQA successfully quantified GP gene transcripts during infection with representative members of 5 ebolavirus species. Immunoblotting verified expression of the soluble GP and the transmembrane GP. Our results defined RNA editing as a general trait of ebolaviruses. The degree of editing, however, varies among ebolaviruses with Reston virus showing the lowest and Bundibugyo virus the highest degree of editing.

Complete Mitochondrial Genome Assembly of an Upland Wild Rice Species, Oryza granulata and Comparative Mitochondrial Genomic Analyses of the Genus Oryza

Wild upland rice species, including Oryza granulata, possess unique characteristics that distinguish them from other Oryza species. For instance, O. granulata characteristically has a GG genome and is accordingly classified as a basal lineage of the genus Oryza. Here, we deployed a versatile hybrid approach by integrating Illumina and PacBio sequencing data to generate a high-quality mitochondrial genome (mitogenome) assembly for O. granulata. The mitogenome of O. granulata was 509,311 base pairs (bp) with sixty-seven genes comprising two circular chromosomes, five ribosomal RNA (rRNA) coding genes, twenty-five transfer RNA (tRNA) coding genes, and thirty-seven genes coding for proteins. We identified a total of 378 simple sequence repeats (SSRs). The genome also contained 643 pairs of dispersed repeats comprising 340 palindromic and 303 forward. In the O. granulata mitogenome, the length of 57 homologous fragments in the chloroplast genome occupied 5.96% of the mitogenome length. Collinearity analysis of three Oryza mitogenomes revealed high structural variability and frequent rearrangements. Phylogenetic analysis showed that, compared to other related genera, O. granulata had the closest genetic relationship with mitogenomes reported for all members of Oryza, and occupies a position at the base of the Oryza phylogeny. Comparative analysis of complete mitochondrial genome assemblies for Oryza species revealed high levels of mitogenomic diversity, providing a foundation for future conservation and utilization of wild rice biodiversity.

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.

Identification of microRNA editing sites in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is a malignant tumor originating from the renal tubular epithelium. Although the microRNAs (miRNAs) transcriptome of ccRCC has been extensively studied, the role of miRNAs editing in ccRCC is largely unknown. By analyzing small RNA sequencing profiles of renal tissues of 154 ccRCC patients and 22 normal controls, we identified 1025 miRNA editing sites from 246 pre-miRNAs. There were 122 editing events with significantly different editing levels in ccRCC compared to normal samples, which include two A-to-I editing events in the seed regions of hsa-mir-376a-3p and hsa-mir-376c-3p, respectively, and one C-to-U editing event in the seed region of hsa-mir-29c-3p. After comparing the targets of the original and edited miRNAs, we found that hsa-mir-376a-1_49g, hsa-mir-376c_48g and hsa-mir-29c_59u had many new targets, respectively. Many of these new targets were deregulated in ccRCC, which might be related to the different editing levels of hsa-mir-376a-3p, hsa-mir-376c-3p, hsa-mir-29c-3p in ccRCC compared to normal controls. Our study sheds new light on miRNA editing events and their potential biological functions in ccRCC.

Benchmarking RNA Editing Detection Tools

RNA, like DNA and proteins, can undergo modifications. To date, over 170 RNA modifications have been identified, leading to the emergence of a new research area known as epitranscriptomics. RNA editing is the most frequent RNA modification in mammalian transcriptomes, and two types have been identified: (1) the most frequent, adenosine to inosine (A-to-I); and (2) the less frequent, cysteine to uracil (C-to-U) RNA editing. Unlike other epitranscriptomic marks, RNA editing can be readily detected from RNA sequencing (RNA-seq) data without any chemical conversions of RNA before sequencing library preparation. Furthermore, analyzing RNA editing patterns from transcriptomic data provides an additional layer of information about the epitranscriptome. As the significance of epitranscriptomics, particularly RNA editing, gains recognition in various fields of biology and medicine, there is a growing interest in detecting RNA editing sites (RES) by analyzing RNA-seq data. To cope with this increased interest, several bioinformatic tools are available. However, each tool has its advantages and disadvantages, which makes the choice of the most appropriate tool for bench scientists and clinicians difficult. Here, we have benchmarked bioinformatic tools to detect RES from RNA-seq data. We provide a comprehensive view of each tool and its performance using previously published RNA-seq data to suggest recommendations on the most appropriate for utilization in future studies.

Assembly of the Complete Mitochondrial Genome of Pereskia aculeata Revealed That Two Pairs of Repetitive Elements Mediated the Recombination of the Genome

Pereskia aculeata is a potential new crop species that has both food and medicinal (antinociceptive activity) properties. However, comprehensive genomic research on P. aculeata is still lacking, particularly concerning its organelle genome. In this study, P. aculeata was studied to sequence the mitochondrial genome (mitogenome) and to ascertain the assembly, informational content, and developmental expression of the mitogenome. The findings revealed that the mitogenome of P. aculeata is circular and measures 515,187 bp in length with a GC content of 44.05%. It contains 52 unique genes, including 33 protein-coding genes, 19 tRNA genes, and three rRNA genes. Additionally, the mitogenome analysis identified 165 SSRs, primarily consisting of tetra-nucleotides, and 421 pairs of dispersed repeats with lengths greater than or equal to 30, which were mainly forward repeats. Based on long reads and PCR experiments, we confirmed that two pairs of long-fragment repetitive elements were highly involved with the mitogenome recombination process. Furthermore, there were 38 homologous fragments detected between the mitogenome and chloroplast genome, and the longest fragment was 3962 bp. This is the first report on the mitogenome in the family Cactaceae. The decoding of the mitogenome of P. aculeata will provide important genetic materials for phylogenetic studies of Cactaceae and promote the utilization of species germplasm resources.

The mitochondrial genome of Heterosentis pseudobagri (Wang & Zhang, 1987) Pichelin & Cribb, 1999 reveals novel aspects of tRNA genes evolution in Acanthocephala

BACKGROUND

Acanthocephala is a clade of obligate endoparasites whose mitochondrial genomes (mitogenomes) and evolution remain relatively poorly understood. Previous studies reported that atp8 is lacking from acanthocephalan mitogenomes, and that tRNA genes often have nonstandard structures. Heterosentis pseudobagri (Arhythmacanthidae) is an acanthocephalan fish endoparasite for which no molecular data are currently available, and biological information is unavailable in the English language. Furthermore, there are currently no mitogenomes available for Arhythmacanthidae.

METHODS

We sequenced its mitogenome and transcriptome, and conducted comparative mitogenomic analyses with almost all available acanthocephalan mitogenomes.

RESULTS

The mitogenome had all genes encoded on the same strand and unique gene order in the dataset. Among the 12 protein-coding genes, several genes were highly divergent and annotated with difficulty. Moreover, several tRNA genes could not be identified automatically, so we had to identify them manually via a detailed comparison with orthologues. As common in acanthocephalans, some tRNAs lacked either the TWC arm or the DHU arm, but in several cases, we annotated tRNA genes only on the basis of the conserved narrow central segment comprising the anticodon, while the flanking 5' and 3' ends did not exhibit any resemblance to orthologues and they could not be folded into a tRNA secondary structure. We corroborated that these are not sequencing artefacts by assembling the mitogenome from transcriptomic data. Although this phenomenon was not observed in previous studies, our comparative analyses revealed the existence of highly divergent tRNAs in multiple acanthocephalan lineages.

CONCLUSIONS

These findings indicate either that multiple tRNA genes are non-functional or that (some) tRNA genes in (some) acanthocephalans might undergo extensive posttranscriptional tRNA processing which restores them to more conventional structures. It is necessary to sequence mitogenomes from yet unrepresented lineages and further explore the unusual patterns of tRNA evolution in Acanthocephala.

Complete mitochondrial genome assembly and comparison of Camellia sinensis var. Assamica cv. Duntsa

Camellia sinensis var. Assamica cv. Duntsa (C.duntsa), a valuable Theaceae from Hunan Province, has been looked at as a precious tea resource by local farmers because of its economic and ecological value. Genomics study on C.duntsa is essential for the domestication and enhancement of tea tree varieties. In the present study, we used a hybrid approach based on Illumina and PacBio data to sequence and assemble the mitochondrial genome of C.duntsa. The mitochondrial genome of C.duntsa was estimated to be 1,081,996 base pairs (bp) and eighty-one genes consisting of one pseudogene, three ribosomal RNA (rRNA) genes, thirty transfer RNA (tRNA) genes, and forty-seven protein-coding genes (PCGs). Tetramer repetitions made up 43.90% of simple sequence repeats (SSRs). The codon usage bias of the Theaceae mitochondrial gene atp9 was altered by mutation, but the codon usage of other genes was shaped by natural selection. Besides, there are eighteen gene-containing homologous regions between the chloroplast and mitochondrial genomes of C. duntsa.Some genomes including atp8, cox1, cox3, nad7, nad9, rpl16, rpl2, rps19, rps4, and sdh4 are absent in the mitochondrial genome of several Theaceae plant. However, C. duntsa maintains these genes integrity and functionality. Another gene, rps16, is either lacking from the mitochondrial genome of C. duntsa or is present as a pseudogene. C. duntsa and C. sinensis (OM809792) are very similar, as shown by a collinear match across four species of Theaceae; the most conservative genes are nad5, atp9, cox2, rps3, trnA-TGC, trnI-GAT, rrn18, trnV-GAC, and ccmFN. Similarly, the genome's phylogenetic trees revealed that C. duntsa was the sister species to C. sinensis. The results confirmed that the C. duntsa and C. sinensis (OM809792) mitochondrial genome underwent gene rearrangement.In general, our results shows that genomic information from organelles can help us understand plant phylogeny and can also be used to make molecular markers and study how genetic traits change over time. Our research will contribute to the population genetics and evolution of tea plant.

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

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

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.

RNA-Binding Proteins: Emerging Therapeutics for Vascular Dysfunction

Vascular diseases account for a significant number of deaths worldwide, with cardiovascular diseases remaining the leading cause of mortality. This ongoing, ever-increasing burden has made the need for an effective treatment strategy a global priority. Recent advances in regenerative medicine, largely the derivation and use of induced pluripotent stem cell (iPSC) technologies as disease models, have provided powerful tools to study the different cell types that comprise the vascular system, allowing for a greater understanding of the molecular mechanisms behind vascular health. iPSC disease models consequently offer an exciting strategy to deepen our understanding of disease as well as develop new therapeutic avenues with clinical translation. Both transcriptional and post-transcriptional mechanisms are widely accepted to have fundamental roles in orchestrating responses to vascular damage. Recently, iPSC technologies have increased our understanding of RNA-binding proteins (RBPs) in controlling gene expression and cellular functions, providing an insight into the onset and progression of vascular dysfunction. Revelations of such roles within vascular disease states have therefore allowed for a greater clarification of disease mechanisms, aiding the development of novel therapeutic interventions. Here, we discuss newly discovered roles of RBPs within the cardio-vasculature aided by iPSC technologies, as well as examine their therapeutic potential, with a particular focus on the Quaking family of isoforms.

Genetic and Functional Analyses of Archaeal ATP-Dependent RNA Ligase in C/D Box sRNA Circularization and Ribosomal RNA Processing

RNA ligases play important roles in repairing and circularizing RNAs post-transcriptionally. In this study, we generated an allelic knockout of ATP-dependent RNA ligase (Rnl) in the hyperthermophilic archaeon Thermococcus kodakarensis to identify its biological targets. A comparative analysis of circular RNA reveals that the Rnl-knockout strain represses circularization of C/D box sRNAs without affecting the circularization of tRNA and rRNA processing intermediates. Recombinant archaeal Rnl could circularize C/D box sRNAs with a mutation in the conserved C/D box sequence element but not when the terminal stem structures were disrupted, suggesting that proximity of the two ends could be critical for intramolecular ligation. Furthermore, T. kodakarensis accumulates aberrant RNA fragments derived from ribosomal RNA in the absence of Rnl. These results suggest that Rnl is responsible for C/D box sRNA circularization and may also play a role in ribosomal RNA processing.

Structural characterization of Platanthera ussuriensis chloroplast genome and comparative analyses with other species of Orchidaceae

BACKGROUND

The genus Tulotis has been classified into the genus Platanthera in the present taxonomic studies since the morphological characteristics of this genus is very similar to that of Platanthera. Platanthera ussuriensis, formerly named as Tulotis ussuriensis, is a small terrestrial orchid species and has been listed as wild plant under State protection (category II) in China. An improved understanding of the genomic information will enable future applications of conservation strategy as well as phylogenetic studies for this rare orchid species. The objective of this research was to characterize and compare the chloroplast genome of P. ussuriensis with other closely related species of Orchidaceae.

RESULTS

The chloroplast genome sequence of P. ussuriensis is 155,016 bp in length, which included a pair of inverted repeats (IRs) of 26,548 bp that separated a large single copy (LSC) region of 83,984 bp and a small single copy (SSC) region of 17,936 bp. The annotation contained a total of 132 genes, including 86 protein-coding genes, 38 tRNA genes and 8 rRNA genes. The simple sequence repeat (SSR) analysis showed that there were 104 SSRs in the chloroplast genome of P. ussuriensis. RNA editing sites recognition indicated 72 RNA editing events occurred, and all codon changes were C to T conversions. Comparative genomics showed that the chloroplast sequence of Platanthera related species were relatively conserved, while there were still some high variation regions that could be used as molecular markers. Moreover, Platanthera related species showed similar IR/SSC and IR/LSC borders. The phylogenetic analysis suggested that P. ussuriensis had a closer evolutionary relationship with P. japonica followed by the remaining Platanthera species.

CONCLUSION

Orchidaceae is a key group of biodiversity protection and also a hot spot group in the plant taxonomy and evolution studies due to their characteristics of high specialization and rapid evolution. This research determined the complete chloroplast genome of P. ussuriensis for the first time, and compared the sequence with other closely related orchid species. These results provide a foundation for future genomic and molecular evolution of the Orchidaceae species, and provide insights into the development of conservation strategy for Platanthera species.

Regulation of MicroRNAs

MicroRNAs are RNAs of about 18-24 nucleotides in lengths, which are found in the small noncoding RNA class and have a crucial role in the posttranscriptional regulation of gene expression, cellular metabolic pathways, and developmental events. These small but essential molecules are first processed by Drosha and DGCR8 in the nucleus and then released into the cytoplasm, where they cleaved by Dicer to form the miRNA duplex. These duplexes are bound by the Argonaute (AGO) protein to form the RNA-induced silencing complex (RISC) in a process called RISC loading. Transcription of miRNAs, processing with Drosha and DGCR8 in the nucleus, cleavage by Dicer, binding to AGO proteins and forming RISC are the most critical steps in miRNA biogenesis. Additional molecules involved in biogenesis at these stages can enhance or inhibit these processes, which can radically change the fate of the cell. Biogenesis is regulated by many checkpoints at every step, primarily at the transcriptional level, in the nucleus, cytoplasm, with RNA regulation, RISC loading, miRNA strand selection, RNA methylation/uridylation, and turnover rate. Moreover, in recent years, different regulation mechanisms have been discovered in noncanonical Drosha or Dicer-independent pathways. This chapter seeks answers to how miRNA biogenesis and function are regulated through both canonical and non-canonical pathways.

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

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

Artificial RNA Editing with ADAR for Gene Therapy

Editing mutated genes is a potential way for the treatment of genetic diseases. G-to-A mutations are common in mammals and can be treated by adenosine-to-inosine (A-to-I) editing, a type of substitutional RNA editing. The molecular mechanism of A-to-I editing involves the hydrolytic deamination of adenosine to an inosine base; this reaction is mediated by RNA-specific deaminases, adenosine deaminases acting on RNA (ADARs), family protein. Here, we review recent findings regarding the application of ADARs to restoring the genetic code along with different approaches involved in the process of artificial RNA editing by ADAR. We have also addressed comparative studies of various isoforms of ADARs. Therefore, we will try to provide a detailed overview of the artificial RNA editing and the role of ADAR with a focus on the enzymatic site directed A-to-I editing.

Genetic variation and microRNA targeting of A-to-I RNA editing fine tune human tissue transcriptomes

BACKGROUND

A-to-I RNA editing diversifies the transcriptome and has multiple downstream functional effects. Genetic variation contributes to RNA editing variability between individuals and has the potential to impact phenotypic variability.

RESULTS

We analyze matched genetic and transcriptomic data in 49 tissues across 437 individuals to identify RNA editing events that are associated with genetic variation. Using an RNA editing quantitative trait loci (edQTL) mapping approach, we identify 3117 unique RNA editing events associated with a cis genetic polymorphism. Fourteen percent of these edQTL events are also associated with genetic variation in their gene expression. A subset of these events are associated with genome-wide association study signals of complex traits or diseases. We determine that tissue-specific levels of ADAR and ADARB1 are able to explain a subset of tissue-specific edQTL events. We find that certain microRNAs are able to differentiate between the edited and unedited isoforms of their targets. Furthermore, microRNAs can generate an expression quantitative trait loci (eQTL) signal from an edQTL locus by microRNA-mediated transcript degradation in an editing-specific manner. By integrative analyses of edQTL, eQTL, and microRNA expression profiles, we computationally discover and experimentally validate edQTL-microRNA pairs for which the microRNA may generate an eQTL signal from an edQTL locus in a tissue-specific manner.

CONCLUSIONS

Our work suggests a mechanism in which RNA editing variability can influence the phenotypes of complex traits and diseases by altering the stability and steady-state level of critical RNA molecules.

Assembly and comparative analysis of the complete mitochondrial genome of Suaeda glauca

Background

Suaeda glauca (S. glauca) is a halophyte widely distributed in saline and sandy beaches, with strong saline-alkali tolerance. It is also admired as a landscape plant with high development prospects and scientific research value. The S. glauca chloroplast (cp) genome has recently been reported; however, the mitochondria (mt) genome is still unexplored.

Results

The mt genome of S. glauca were assembled based on the reads from Pacbio and Illumina sequencing platforms. The circular mt genome of S. glauca has a length of 474,330 bp. The base composition of the S. glauca mt genome showed A (28.00%), T (27.93%), C (21.62%), and G (22.45%). S. glauca mt genome contains 61 genes, including 27 protein-coding genes, 29 tRNA genes, and 5 rRNA genes. The sequence repeats, RNA editing, and gene migration from cp to mt were observed in S. glauca mt genome. Phylogenetic analysis based on the mt genomes of S. glauca and other 28 taxa reflects an exact evolutionary and taxonomic status of S. glauca. Furthermore, the investigation on mt genome characteristics, including genome size, GC contents, genome organization, and gene repeats of S. gulaca genome, was investigated compared to other land plants, indicating the variation of the mt genome in plants. However, the subsequently Ka/Ks analysis revealed that most of the protein-coding genes in mt genome had undergone negative selections, reflecting the importance of those genes in the mt genomes.

Conclusions

In this study, we reported the mt genome assembly and annotation of a halophytic model plant S. glauca. The subsequent analysis provided us a comprehensive understanding of the S. glauca mt genome, which might facilitate the research on the salt-tolerant plant species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07490-9.

Evolutionary history of cotranscriptional editing in the paramyxoviral phosphoprotein gene

The phosphoprotein gene of the paramyxoviruses encodes multiple protein products. The P, V, and W proteins are generated by transcriptional slippage. This process results in the insertion of non-templated guanosine nucleosides into the mRNA at a conserved edit site. The P protein is an essential component of the viral RNA polymerase and is encoded by a faithful copy of the gene in the majority of paramyxoviruses. However, in some cases, the non-essential V protein is encoded by default and guanosines must be inserted into the mRNA in order to encode P. The number of guanosines inserted into the P gene can be described by a probability distribution, which varies between viruses. In this article, we review the nature of these distributions, which can be inferred from mRNA sequencing data, and reconstruct the evolutionary history of cotranscriptional editing in the paramyxovirus family. Our model suggests that, throughout known history of the family, the system has switched from a P default to a V default mode four times; complete loss of the editing system has occurred twice, the canonical zinc finger domain of the V protein has been deleted or heavily mutated a further two times, and the W protein has independently evolved a novel function three times. Finally, we review the physical mechanisms of cotranscriptional editing via slippage of the viral RNA polymerase.

SMK6 mediates the C-to-U editing at multiple sites in maize mitochondria

The recently identified PPR-E+/NVWA/DYW2 RNA editing complex provides insights into the mechanism of RNA editing in higher plant organelles. However, whether the complex works together with the previously identified editing factors RIPs/MORFs is unclear. In this paper, we identified a maize Smk6 gene, which encodes a mitochondrion-targeted PPR-E+protein with E1 and E2 domains at the C terminus. Loss of Smk6 function affects the C-to-U editing at nad1-740, nad4L-110, nad7-739, and mttB-138,139 sites, impairs mitochondrial activity and blocks embryogenesis and endosperm development. Genetic and molecular analysis indicated that SMK6 is the maize ortholog of the Arabidopsis SLO2, which is a component of the PPR-E+/NVWA/DYW2 editing complex. However, yeast two-hybrid analyses did not detect any interaction between SMK6 and any of the mitochondrion-targeted RIPs/MORFs, suggesting that RIPs/MORFs may not be a component of PPR-E+/NVWA/DYW2 RNA editing complex. Further analyses are required to provide evidence that RIP/MORFs and SMK6 do not physically interact in vivo.

FEDRO: a software tool for the automatic discovery of candidate ORFs in plants with c →u RNA editing

BACKGROUND

RNA editing is an important mechanism for gene expression in plants organelles. It alters the direct transfer of genetic information from DNA to proteins, due to the introduction of differences between RNAs and the corresponding coding DNA sequences. Software tools successful for the search of genes in other organisms not always are able to correctly perform this task in plants organellar genomes. Moreover, the available software tools predicting RNA editing events utilise algorithms that do not account for events which may generate a novel start codon.

RESULTS

We present FEDRO, a Java software tool implementing a novel strategy to generate candidate Open Reading Frames (ORFs) resulting from Cytidine to Uridine (c→u) editing substitutions which occur in the mitochondrial genome (mtDNA) of a given input plant. The goal is to predict putative proteins of plants mitochondria that have not been yet annotated. In order to validate the generated ORFs, a screening is performed by checking for sequence similarity or presence in active transcripts of the same or similar organisms. We illustrate the functionalities of our framework on a model organism.

CONCLUSIONS

The proposed tool may be used also on other organisms and genomes. FEDRO is publicly available at http://math.unipa.it/rombo/FEDRO .

Genome-wide identification of RNA editing in seven porcine tissues by matched DNA and RNA high-throughput sequencing

BACKGROUND

RNA editing is a co/posttranscriptional modification mechanism that increases the diversity of transcripts, with potential functional consequences. The advent of next-generation sequencing technologies has enabled the identification of RNA edits at unprecedented throughput and resolution. However, our knowledge of RNA editing in swine is still limited.

RESULTS

Here, we utilized RES-Scanner to identify RNA editing sites in the brain, subcutaneous fat, heart, liver, muscle, lung and ovary in three 180-day-old Large White gilts based on matched strand-specific RNA sequencing and whole-genome resequencing datasets. In total, we identified 74863 editing sites, and 92.1% of these sites caused adenosine-to-guanosine (A-to-G) conversion. Most A-to-G sites were located in noncoding regions and generally had low editing levels. In total, 151 A-to-G sites were detected in coding regions (CDS), including 94 sites that could lead to nonsynonymous amino acid changes. We provide further evidence supporting a previous observation that pig transcriptomes are highly editable at PRE-1 elements. The number of A-to-G editing sites ranged from 4155 (muscle) to 25001 (brain) across the seven tissues. The expression levels of the ADAR enzymes could explain some but not all of this variation across tissues. The functional analysis of the genes with tissue-specific editing sites in each tissue revealed that RNA editing might play important roles in tissue function. Specifically, more pathways showed significant enrichment in the fat and liver than in other tissues, while no pathway was enriched in the muscle.

CONCLUSIONS

This study identified a total of 74863 nonredundant RNA editing sites in seven tissues and revealed the potential importance of RNA editing in tissue function. Our findings largely extend the porcine editome and enhance our understanding of RNA editing in swine.

Single-nucleotide variants in human RNA: RNA editing and beyond

Through analysis of paired high-throughput DNA-Seq and RNA-Seq data, researchers quickly recognized that RNA-Seq can be used for more than just gene expression quantification. The alternative applications of RNA-Seq data are abundant, and we are particularly interested in its usefulness for detecting single-nucleotide variants, which arise from RNA editing, genomic variants and other RNA modifications. A stunning discovery made from RNA-Seq analyses is the unexpectedly high prevalence of RNA-editing events, many of which cannot be explained by known RNA-editing mechanisms. Over the past 6-7 years, substantial efforts have been made to maximize the potential of RNA-Seq data. In this review we describe the controversial history of mining RNA-editing events from RNA-Seq data and the corresponding development of methodologies to identify, predict, assess the quality of and catalog RNA-editing events as well as genomic variants.

Regulation of RNA editing by RNA-binding proteins in human cells

Adenosine-to-inosine (A-to-I) editing, mediated by the ADAR enzymes, diversifies the transcriptome by altering RNA sequences. Recent studies reported global changes in RNA editing in disease and development. Such widespread editing variations necessitate an improved understanding of the regulatory mechanisms of RNA editing. Here, we study the roles of >200 RNA-binding proteins (RBPs) in mediating RNA editing in two human cell lines. Using RNA-sequencing and global protein-RNA binding data, we identify a number of RBPs as key regulators of A-to-I editing. These RBPs, such as TDP-43, DROSHA, NF45/90 and Ro60, mediate editing through various mechanisms including regulation of ADAR1 expression, interaction with ADAR1, and binding to Alu elements. We highlight that editing regulation by Ro60 is consistent with the global up-regulation of RNA editing in systemic lupus erythematosus. Additionally, most key editing regulators act in a cell type-specific manner. Together, our work provides insights for the regulatory mechanisms of RNA editing.

Molecular and Functional Diversity of RNA Editing in Plant Mitochondria

RNA editing is a fundamental biochemical process relating to the modification of nucleotides in messenger RNAs of functional genes in cells. RNA editing leads to re-establishment of conserved amino acid residues for functional proteins in nuclei, chloroplasts, and mitochondria. Identification of RNA editing factors that contributes to target site recognition increases our understanding of RNA editing mechanisms. Significant progress has been made in recent years in RNA editing studies for both animal and plant cells. RNA editing in nuclei and mitochondria of animal cells and in chloroplast of plant cells has been extensively documented and reviewed. RNA editing has been also extensively documented on plant mitochondria. However, functional diversity of RNA editing factors in plant mitochondria is not overviewed. Here, we review the biological significance of RNA editing, recent progress on the molecular mechanisms of RNA editing process, and function diversity of editing factors in plant mitochondrial research. We will focus on: (1) pentatricopeptide repeat proteins in Arabidopsis and in crop plants; (2) the progress of RNA editing process in plant mitochondria; (3) RNA editing-related RNA splicing; (4) RNA editing associated flower development; (5) RNA editing modulated male sterile; (6) RNA editing-regulated cell signaling; and (7) RNA editing involving abiotic stress. Advances described in this review will be valuable in expanding our understanding in RNA editing. The diverse functions of RNA editing in plant mitochondria will shed light on the investigation of molecular mechanisms that underlies plant development and abiotic stress tolerance.

The 27 kDa Trypanosoma brucei Pentatricopeptide Repeat Protein is a G-tract Specific RNA Binding Protein

Pentatricopeptide repeat (PPR) proteins, a helical repeat family of organellar RNA binding proteins, play essential roles in post-transcriptional RNA processing. In Trypanosoma brucei, an expanded family of PPR proteins localize to the parasite’s single mitochondrion, where they are believed to perform important roles in both RNA processing and translation. We studied the RNA binding specificity of the simplest T. brucei PPR protein (KRIPP11) using electrophoretic mobility shift assays, fluorescence anisotropy, circular dichroism spectroscopy, and in vitro selection. We found KRIPP11 to be an RNA binding protein with specificity for sequences of four or more consecutive guanosine residues (G-tracts). Such G-tracts are dramatically enriched in T. brucei mitochondrial transcripts that are destined for extensive uridine insertion/deletion editing but are not present in mRNAs following editing. We further found that the quadruplex oligoguanosine RNA conformation is preferentially recognized by KRIPP11 over other conformational forms, and is bound without disruption of the quadruplex structure. In combination with prior data demonstrating association of KRIPP11 with the small ribosomal subunit, these results suggest possible roles for KRIPP11 in bridging mRNA maturation and translation or in facilitating translation of unusual dual-coded open reading frames.

Advances in genomic hepatocellular carcinoma research

Background

Hepatocellular carcinoma (HCC) is the cancer with the second highest mortality in the world due to its late presentation and limited treatment options. As such, there is an urgent need to identify novel biomarkers for early diagnosis and to develop novel therapies. The availability of next-generation sequencing (NGS) data from tumors of liver cancer patients has provided us with invaluable resources to better understand HCC through the integration of data from different sources to facilitate the identification of promising biomarkers or therapeutic targets.

Findings

Here, we review key insights gleaned from more than 20 NGS studies of HCC tumor samples, comprising approximately 582 whole genomes and 1,211 whole exomes mainly from the East Asian population. Through consolidation of reported somatic mutations from multiple studies, we identified genes with different types of somatic mutations, including single nucleotide variations, insertion/deletions, structural variations, and copy number alterations as well as genes with multiple frequent viral integration. Pathway analysis showed that this curated list of somatic mutations is critically involved in cancer-related pathways, viral carcinogenesis, and signaling pathways. Lastly, we addressed the future directions of HCC research as more NGS datasets become available.

Conclusions

Our review is a comprehensive resource for the current NGS research in HCC, consolidating published articles, potential gene candidates, and their related biological pathways.

RNA variant identification discrepancy among splice-aware alignment algorithms

Next-generation sequencing (NGS) techniques have been generating various molecular maps, including transcriptomes via RNA-seq. Although the primary purpose of RNA-seq is to quantify the expression level of known genes, RNA variants are also identifiable. However, care must be taken to account for RNA's dynamic nature. In this study, we evaluated the following popular splice-aware alignment algorithms in the context of RNA variant-calling analysis: HISAT2, STAR, STAR (two-pass mode), Subread, and Subjunc. For this, we performed RNA-seq with ten pieces of invasive ductal carcinoma from breast tissue and three pieces of adjacent normal tissue from a single patient. These RNA-seq data were used to evaluate the performance of splice-aware aligners. Surprisingly, the number of common potential RNA editing sites (pRESs) identified by all alignment algorithms was less than 2% of the total. The main cause of this difference was the mapped reads on the splice junctions. In addition, the RNA quality significantly affected the outcome. Therefore, researchers must consider these experimental and bioinformatic features during RNA variant analysis. Further investigations of common pRESs discovered that BDH1, CCDC137, and TBC1D10A transcripts contained a single non-synonymous RNA variant that was unique to breast cancer tissue compared to adjacent normal tissue; thus, further clinical validation is required.

Evidence of non-tandemly repeated rDNAs and their intragenomic heterogeneity in Rhizophagus irregularis

Arbuscular mycorrhizal fungus (AMF) species are some of the most widespread symbionts of land plants. Our much improved reference genome assembly of a model AMF, Rhizophagus irregularis DAOM-181602 (total contigs = 210), facilitated a discovery of repetitive elements with unusual characteristics. R. irregularis has only ten or 11 copies of complete 45S rDNAs, whereas the general eukaryotic genome has tens to thousands of rDNA copies. R. irregularis rDNAs are highly heterogeneous and lack a tandem repeat structure. These findings provide evidence for the hypothesis that rDNA heterogeneity depends on the lack of tandem repeat structures. RNA-Seq analysis confirmed that all rDNA variants are actively transcribed. Observed rDNA/rRNA polymorphisms may modulate translation by using different ribosomes depending on biotic and abiotic interactions. The non-tandem repeat structure and intragenomic heterogeneity of AMF rDNA/rRNA may facilitate successful adaptation to various environmental conditions, increasing host compatibility of these symbiotic fungi.

Transcription-Associated Mutation Promotes RNA Complexity in Highly Expressed Genes-A Major New Source of Selectable Variation

Alternatively spliced transcript isoforms are thought to play a critical role for functional diversity. However, the mechanism generating the enormous diversity of spliced transcript isoforms remains unknown, and its biological significance remains unclear. We analyzed transcriptomes in saker falcons, chickens, and mice to show that alternative splicing occurs more frequently, yielding more isoforms, in highly expressed genes. We focused on hemoglobin in the falcon, the most abundantly expressed genes in blood, finding that alternative splicing produces 10-fold more isoforms than expected from the number of splice junctions in the genome. These isoforms were produced mainly by alternative use of de novo splice sites generated by transcription-associated mutation (TAM), not by the RNA editing mechanism normally invoked. We found that high expression of globin genes increases mutation frequencies during transcription, especially on nontranscribed DNA strands. After DNA replication, transcribed strands inherit these somatic mutations, creating de novo splice sites, and generating multiple distinct isoforms in the cell clone. Bisulfate sequencing revealed that DNA methylation may counteract this process by suppressing TAM, suggesting DNA methylation can spatially regulate RNA complexity. RNA profiling showed that falcons living on the high Qinghai-Tibetan Plateau possess greater global gene expression levels and higher diversity of mean to high abundance isoforms (reads per kilobases per million mapped reads ≥18) than their low-altitude counterparts, and we speculate that this may enhance their oxygen transport capacity under low-oxygen environments. Thus, TAM-induced RNA diversity may be physiologically significant, providing an alternative strategy in lifestyle evolution.

REDO: RNA Editing Detection in Plant Organelles Based on Variant Calling Results

RNA editing is a post-transcriptional or cotranscriptional process that changes the sequence of the precursor transcript by substitutions, insertions, or deletions. Almost all of the land plants undergo RNA editing in organelles (plastids and mitochondria). Although several software tools have been developed to identify RNA editing events, there has been a great challenge to distinguish true RNA editing events from genome variation, sequencing errors, and other factors. Here we introduce REDO, a comprehensive application tool for identifying RNA editing events in plant organelles based on variant call format files from RNA-sequencing data. REDO is a suite of Perl scripts that illustrate a bunch of attributes of RNA editing events in figures and tables. REDO can also detect RNA editing events in multiple samples simultaneously and identify the significant differential proportion of RNA editing loci. Comparing with similar tools, such as REDItools, REDO runs faster with higher accuracy, and more specificity at the cost of slightly lower sensitivity. Moreover, REDO annotates each RNA editing site in RNAs, whereas REDItools reports only possible RNA editing sites in genome, which need additional steps to obtain RNA editing profiles for RNAs. Overall, REDO can identify potential RNA editing sites easily and provide several functions such as detailed annotations, statistics, figures, and significantly differential proportion of RNA editing sites among different samples.

R/G editing in GluA2Rflop modulates the functional difference between GluA1 flip and flop variants in GluA1/2R heteromeric channels

In α-amino-3-hydroxy-5-methyl-4-isoxazole-propionate (AMPA) receptors, RNA editing and alternative splicing generate sequence variants, and those variants, as in GluA2-4 AMPA receptor subunits, generally show different properties. Yet, earlier studies have shown that the alternatively spliced, flip and flop variants of GluA1 AMPA receptor subunit exhibit no functional difference in homomeric channel form. Using a laser-pulse photolysis technique, combined with whole-cell recording, we measured the rate of channel opening, among other kinetic properties, for a series of AMPA channels with different arginine/glycine (R/G) editing and flip/flop status. We find that R/G editing in the GluA2 subunit modulates the channel properties in both homomeric (GluA2Q) and complex (GluA2Q/2R and GluA1/2R) channel forms. However, R/G editing is only effective in flop channels. Specifically, editing at the R/G site on the GluA2R flop isoform accelerates the rate of channel opening and desensitization for GluA1/2R channels more pronouncedly with the GluA1 being in the flop form than in the flip form; yet R/G editing has no effect on either channel-closing rate or EC₅₀. Our results suggest R/G editing via GluA2R serve as a regulatory mechanism to modulate the function of GluA2R-containing, native receptors involved in fast excitatory synaptic transmission.

Highly diversified expansions shaped the evolution of membrane bound proteins in metazoans

The dramatic increase in membrane proteome complexity is arguably one of the most pivotal evolutionary events that underpins the origin of multicellular animals. However, the origin of a significant number of membrane families involved in metazoan development has not been clarified. In this study, we have manually curated the membrane proteomes of 22 metazoan and 2 unicellular holozoan species. We identify 123,014 membrane proteins in these 24 eukaryotic species and classify 86% of the dataset. We determine 604 functional clusters that are present from the last holozoan common ancestor (LHCA) through many metazoan species. Intriguingly, we show that more than 70% of the metazoan membrane protein families have a premetazoan origin. The data show that enzymes are more highly represented in the LHCA and expand less than threefold throughout metazoan species; in contrast to receptors that are relatively few in the LHCA but expand nearly eight fold within metazoans. Expansions related to cell adhesion, communication, immune defence, and developmental processes are shown in conjunction with emerging biological systems, such as neuronal development, cytoskeleton organization, and the adaptive immune response. This study defines the possible LHCA membrane proteome and describes the fundamental functional clusters that underlie metazoan diversity and innovation.

Non-coding RNAs as a new dawn in tumor diagnosis

The current knowledge about non-coding RNAs (ncRNAs) as important regulators of gene expression in both physiological and pathological conditions, has been the main engine for the design of innovative platforms to finalize the pharmacological application of ncRNAs as either therapeutic tools or as molecular biomarkers in cancer. Biochemical alterations of cancer cells are, in fact, largely supported by ncRNA disregulation in the tumor site, which, in turn, reflects the cancer-associated specific modification of circulating ncRNA expression pattern. The aim of this review is to describe the state of the art of pre-clinical and clinical studies that analyze the involvement of miRNAs and lncRNAs in cancer-related processes, such as proliferation, invasion and metastases, giving emphasis to their functional role. A central node of our work has been also the examination of advantages and criticisms correlated with the clinical use of ncRNAs, taking into account the pressing need to refine the profiling methods aimed at identify novel diagnostic and prognostic markers and the request to optimize the delivery of such nucleic acids for a therapeutic use in an imminent future.

Whole-transcriptome RNA-seq, gene set enrichment pathway analysis, and exon coverage analysis of two plastid RNA editing mutants

In land plants, plastid and mitochondrial RNAs are subject to post-transcriptional C-to-U RNA editing. T-DNA insertions in the ORGANELLE RNA RECOGNITION MOTIF PROTEIN6 gene resulted in reduced photosystem II (PSII) activity and smaller plant and leaf sizes. Exon coverage analysis of the ORRM6 gene showed that orrm6-1 and orrm6-2 are loss-of-function mutants. Compared to other ORRM proteins, ORRM6 affects a relative small number of RNA editing sites. Sanger sequencing of reverse transcription-PCR products of plastid transcripts revealed 2 plastid RNA editing sites that are substantially affected in the orrm6 mutants: psbF-C77 and accD-C794. The psbF gene encodes the β subunit of cytochrome b₅₅₉, an essential component of PSII. The accD gene encodes the β subunit of acetyl-CoA carboxylase, a protein required in plastid fatty acid biosynthesis. Whole-transcriptome RNA-seq demonstrated that editing at psbF-C77 is nearly absent and the editing extent at accD-C794 was significantly reduced. Gene set enrichment pathway analysis showed that expression of multiple gene sets involved in photosynthesis, especially photosynthetic electron transport, is significantly upregulated in both orrm6 mutants. The upregulation could be a mechanism to compensate for the reduced PSII electron transport rate in the orrm6 mutants. These results further demonstrated that Organelle RNA Recognition Motif protein ORRM6 is required in editing of specific RNAs in the Arabidopsis (Arabidopsis thaliana) plastid.

The Evolving Definition of the Term "Gene"

This paper presents a history of the changing meanings of the term "gene," over more than a century, and a discussion of why this word, so crucial to genetics, needs redefinition today. In this account, the first two phases of 20th century genetics are designated the "classical" and the "neoclassical" periods, and the current molecular-genetic era the "modern period." While the first two stages generated increasing clarity about the nature of the gene, the present period features complexity and confusion. Initially, the term "gene" was coined to denote an abstract "unit of inheritance," to which no specific material attributes were assigned. As the classical and neoclassical periods unfolded, the term became more concrete, first as a dimensionless point on a chromosome, then as a linear segment within a chromosome, and finally as a linear segment in the DNA molecule that encodes a polypeptide chain. This last definition, from the early 1960s, remains the one employed today, but developments since the 1970s have undermined its generality. Indeed, they raise questions about both the utility of the concept of a basic "unit of inheritance" and the long implicit belief that genes are autonomous agents. Here, we review findings that have made the classic molecular definition obsolete and propose a new one based on contemporary knowledge.