REDItools: high-throughput RNA editing detection made easy

Advances in RNA editing in hematopoiesis and associated malignancies

ABSTRACT

Adenosine-to-inosine (A-to-I) RNA editing is a prevalent RNA modification essential for cell survival. The process is catalyzed by the adenosine deaminase acting on RNA (ADAR) enzyme family that converts adenosines in double-stranded RNAs (dsRNAs) into inosines, which are read as guanosines during translation. Deep sequencing has helped to reveal that A-to-I editing occurs across various types of RNAs, affecting their functions. RNA editing detection is now so sophisticated that we can achieve a high level of accuracy and sensitivity to identify low-abundance edited events. Consequently, A-to-I editing has been implicated in various biological processes, including immune and stress responses, cancer progression, and stem cell fate determination. In particular, a crucial role for this process has been recently reported in hematopoietic cell development and hematologic malignancy progression. Results from genetic mouse models have demonstrated the impact of ADARs' catalytic activity on hematopoietic cells, complemented by insights from human cell studies. Meanwhile, clinical studies have implicated ADAR enzymes and RNA editing events in hematologic malignancies and highlighted their potential as prognostic indicators. In this review, we outline the regulatory mechanisms of RNA editing in both normal hematopoiesis and hematologic malignancies. We then speculate on how targeting ADAR expression and site-specific RNA substrates might serve as a therapeutic avenue for affected patients.

Unveiling the mitochondrial genome of Salvia splendens insights into the evolutionary traits within the genus Salvia

Previously, we resolved the complete sequences of the mitochondrial genomes (mitogenome) of two Salvia species (S. miltiorrhiza and S. officinalis). The major configurations of these two species were two circular chromosomes. In this study, we further studied the mitogenome of a floral species of Salvia (Salvia splendens) to understand the diversity and evolution of the Salvia mitogenomes. We sequenced the total DNAs of S. splendens using the Nanopore and Illumina platforms and assembled the mitogenome using a hybrid assembly strategy. The major configurations of the S. splendens were two circular chromosomes with lengths of 182,239 and 165,055 bp. There were 32 protein-coding genes (PCGs), three rRNA genes, and 18 tRNA genes annotated in the S. splendens mitogenome. We found 56 pairs of repetitive sequences in the S. splendens mitogenome. Three of them (R01, 04, and 07) could mediate recombination, whose products could be identified by the mapping of Nanopore reads, PCR amplifications, and Sanger sequencing of the PCR products. 457 RNA editing sites were identified in the S. splendens mitochondrial RNAs when comparing the RNA-seq data with their corresponding DNA templates. We showed that S. splendens was a sister taxon to S. miltiorrhiza based on the mitogenomes, consistent with the phylogeny determined with the plastome sequences. Crucially, we developed 12 mitochondrial markers sourced from mitochondrial intron regions to facilitate the identification of three Salvia species. Our study offers a comprehensive view of the structure of the Salvia mitogenomes and provides robust mitochondrial markers for Salvia species identification.

RNA editing of genomic neighbors controls antiviral response in fungi

Virus symptom expression involves complex interactions between viruses and their hosts, including antiviral defenses and counter-defenses, many of which are not well understood. This study utilizes Neurospora crassa as a model organism to investigate the role of RNA editing in the fungal antiviral response. We identify two adjacent genes in the genome: the A-to-I RNA-editing enzyme, OTT_1508-like deaminase (old), and its target, zinc fingers adjacent to old (zao). These genes regulate the transcriptional response to viral infection, with old modulating the expression of zao, which functions as a master transcription factor. This regulation contributes to asymptomatic infections by maintaining normal growth and development. However, in RNAi-deficient conditions, the overactivation of these genes leads to severe symptoms, akin to hypersensitive responses observed in plants. Additionally, homologs of zao-old are found as genomic neighbors in various filamentous ascomycetes, suggesting that this RNA-editing system may represent an evolutionarily conserved antiviral mechanism.

Porous Silicon Nanoneedles Efficiently Deliver Adenine Base Editor to Correct a Recurrent Pathogenic COL7A1 Variant in Recessive Dystrophic Epidermolysis Bullosa

Base editing, a CRISPR-based genome editing technology, enables precise correction of single-nucleotide variants, promising resolutive treatment for monogenic genetic disorders like recessive dystrophic epidermolysis bullosa (RDEB). However, the application of base editors in cell manufacturing is hindered by inconsistent efficiency and high costs, contributed by suboptimal delivery methods. Nanoneedles have emerged as an effective delivery approach, enabling highly efficient, non-perturbing gene therapies both in vitro and in vivo. Here we demonstrate that nanoneedle delivery of an adenine base editor corrects a heterozygous single-nucleotide pathogenic variant in COL7A1 in primary RDEB fibroblasts in vitro with 96.5% efficiency, without inducing off-target variants. The nanoneedle delivery maintains cell viability and displays modest phenotypical alterations unlike conventional cationic lipid transfection. The nanoneedle-mediated editing significantly increases the production and secretion of full-length type VII collagen protein, contributing to restore functional fibroblasts phenotype by improving cell adhesion. These findings underscore the suitability and safety of nanoneedles for gene editing in a clinically relevant context of cell manufacturing, establishing a foundation for their use in cell therapies.

Post-transcriptional modular synthetic receptors

Inspired by the power of transcriptional synthetic receptors and hoping to complement them to expand the toolbox for cell engineering, we establish LIDAR (Ligand-Induced Dimerization-Activating RNA editing), a modular post-transcriptional synthetic receptor platform that harnesses RNA editing by adenosine deaminases acting on RNA. LIDAR is compatible with various receptor architectures in different cellular contexts and enables the sensing of diverse ligands and the production of functional outputs. Furthermore, LIDAR can sense orthogonal signals in the same cell and produce synthetic spatial patterns, potentially enabling the programming of complex multicellular behaviors. Lastly, LIDAR is compatible with compact encoding and can be delivered as synthetic mRNA. Thus, LIDAR expands the family of synthetic receptors, holding the promise to empower basic research and therapeutic applications.

Assembly and comparative analysis of the complete mitochondrial of Spodiopogon sagittifolius, an endemic and protective species from Yunnan, China

BACKGROUND

Spodiopogon sagittifolius, a C4 plant closely related to cultivated crops, is an edible resource and a Class II nationally protected species in China. Endemic to Yunnan, its populations are declining due to habitat destruction, highlighting its resource and conservation importance. Despite its significance, the molecular phylogenetic relationships and genetic mechanisms of adaptive evolution in the genus Spodiopogon remain poorly understood.

RESULTS

We successfully assembled and annotated the first mitochondrial genome of S. sagittifolius using HiFi sequencing and the PMAT tool. The genome is 500,699 bp in length with a GC content of 43.15%. Synteny and dN/dS analyses revealed structural and functional conservation of mitochondrial genomes in closely related species, with most protein-coding genes under purifying selection (dN/dS < 1). Notably, nad2 exhibited signs of positive selection (dN/dS = 1.49), indicating potential adaptive evolution. Extensive RNA editing events were detected across 27 protein-coding genes, predominantly involving C-to-U conversions, with synonymous mutations accounting for 49.65% of the edits. Strong codon usage bias favoring A/U-ending codons and the identification of repeat sequences suggest enhanced mitochondrial efficiency and stress adaptation. Phylogenetic analyses confirmed the taxonomic placement of S. sagittifolius within the Andropogoneae tribe.

CONCLUSIONS

This study provides the first insights into the mitochondrial genome evolution of S. sagittifolius, highlighting key features linked to stress tolerance and adaptive evolution. These findings establish a foundation for its conservation and potential domestication, with implications for crop improvement and ecological resilience.

REDInet: a temporal convolutional network-based classifier for A-to-I RNA editing detection harnessing million known events

A-to-I ribonucleic acid (RNA) editing detection is still a challenging task. Current bioinformatics tools rely on empirical filters and whole genome sequencing or whole exome sequencing data to remove background noise, sequencing errors, and artifacts. Sometimes they make use of cumbersome and time-consuming computational procedures. Here, we present REDInet, a temporal convolutional network-based deep learning algorithm, to profile RNA editing in human RNA sequencing (RNAseq) data. It has been trained on REDIportal RNA editing sites, the largest collection of human A-to-I changes from >8000 RNAseq data of the genotype-tissue expression project. REDInet can classify editing events with high accuracy harnessing RNAseq nucleotide frequencies of 101-base windows without the need for coupled genomic data.

Precise modelling of mitochondrial diseases using optimized mitoBEs

The development of animal models is crucial for studying and treating mitochondrial diseases. Here we optimized adenine and cytosine deaminases to reduce off-target effects on the transcriptome and the mitochondrial genome, improving the accuracy and efficiency of our newly developed mitochondrial base editors (mitoBEs)1. Using these upgraded mitoBEs (version 2 (v2)), we targeted 70 mouse mitochondrial DNA mutations analogous to human pathogenic variants2, establishing a foundation for mitochondrial disease mouse models. Circular RNA-encoded mitoBEs v2 achieved up to 82% editing efficiency in mice without detectable off-target effects in the nuclear genome. The edited mitochondrial DNA persisted across various tissues and was maternally inherited, resulting in F1 generation mice with mutation loads as high as 100% and some mice exhibiting editing only at the target site. By optimizing the transcription activator-like effector (TALE) binding site, we developed a single-base-editing mouse model for the mt-Nd5 A12784G mutation. Phenotypic evaluations led to the creation of mouse models for the mt-Atp6 T8591C and mt-Nd5 A12784G mutations, exhibiting phenotypes corresponding to the reduced heart rate seen in Leigh syndrome and the vision loss characteristic of Leber's hereditary optic neuropathy, respectively. Moreover, the mt-Atp6 T8591C mutation proved to be more deleterious than mt-Nd5 A12784G, affecting embryonic development and rapidly diminishing through successive generations. These upgraded mitoBEs offer a highly efficient and precise strategy for constructing mitochondrial disease models, laying a foundation for further research in this field.

Alternative splicing of transposable elements in human breast cancer

Transposable elements (TEs) drive genome evolution and can affect gene expression through diverse mechanisms. In breast cancer, disrupted regulation of TE sequences may facilitate tumor-specific transcriptomic alterations. We examine 142,514 full-length isoforms derived from long-read RNA sequencing (LR-seq) of 30 breast samples to investigate the effects of TEs on the breast cancer transcriptome. Approximately half of these isoforms contain TE sequences, and these contribute to half of the novel annotated splice junctions. We quantify splicing of these LR-seq derived isoforms in 1,135 breast tumors from The Cancer Genome Atlas (TCGA) and 1,329 healthy tissue samples from the Genotype-Tissue Expression (GTEx), and find 300 TE-overlapping tumor-specific splicing events. Some splicing events are enriched in specific breast cancer subtypes - for example, a TE-driven transcription start site upstream of ERBB2 in HER2 + tumors, and several TE-mediated splicing events are associated with patient survival and poor prognosis. The full-length sequences we capture with LR-seq reveal thousands of isoforms with signatures of RNA editing, including a novel isoform belonging to RHOA; a gene previously implicated in tumor progression. We utilize our full-length isoforms to discover polymorphic TE insertions that alter splicing and validate one of these events in breast cancer cell lines. Together, our results demonstrate the widespread effects of dysregulated TEs on breast cancer transcriptomes and highlight the advantages of long-read isoform sequencing for understanding TE biology. TE-derived isoforms may alter the expression of genes important in cancer and can potentially be used as novel, disease-specific therapeutic targets or biomarkers.One sentence summary: Transposable elements generate alternative isoforms and alter post-transcriptional regulation in human breast cancer.

The Taihangia mitogenome provides new insights into its adaptation and organelle genome evolution in Rosaceae

MAIN CONCLUSION

We present the first Taihangia mitogenome, uncovering frequent rearrangements and significant length variation in Rosaceae, likely driven by hybridization and repeat content, alongside widespread mito-chloroplast phylogenetic conflicts. Taihangia, an ancient and endangered monotypic genus within the subfamily Rosoideae of the family Rosaceae, is endemic to cliffs and serves as an ideal material for studying the adaptations of cliff-dwelling plants and the evolutionary processes of the Rosaceae family. In this study, the mitogenome and plastome of T. rupestris var. ciliata were assembled, with lengths of 265,633 bp and 155,467 bp, both exhibiting typical circular structures. Positive selection was detected in the nad4L and sdh4 genes, likely playing a role in adaptation to harsh environments. Comparative genomic analysis indicated that repetitive sequences are likely the main contributors to genome size variation in Rosaceae and also influence horizontal gene transfer between organelle genomes. In T. rupestris var. ciliata, 20 mitochondrial plastid DNA sequences were identified, including 16 complete plastid genes. Moreover, frequent rearrangements were observed in the non-coding regions of mitogenome within the subfamily Rosoideae, potentially linked to the complex evolutionary history and the presence of repetitive sequences. In contrast, coding regions remained highly conserved (over 83% similarity) to maintain essential mitochondrial functions. Phylogenomic analysis of the two organelle genomes revealed conflicts in the phylogenetic relationships within Rosaceae, potentially due to the inconsistent mutation rates and frequent hybridization events in the evolutionary history of the family. In conclusion, the organelle genome analysis of Taihangia provides crucial genomic resources for understanding the evolution and adaptation of Rosaceae species.

Features and evolutionary adaptations of the mitochondrial genome of Garuga forrestii W. W. Sm

Introduction

Garuga forrestii W. W. Sm. is a tree species of the Burseraceae family, endemic to China, found in hot/warm-dry valleys. This species plays a crucial role in maintaining biodiversity in these ecosystems.

Methods

We performed de novo assembly of the Garuga forrestii mitochondrial genome using PMAT (v.1.5.4), resulting in a typical circular molecule of 606,853 bp. The genome consists of 31 tRNA genes, 3 rRNA genes, 35 protein-coding genes, and 1 pseudogene. The study also investigates RNA editing sites and evolutionary patterns.

Results

The mitochondrial genome exhibits a low proportion of repetitive sequences (3.30%), suggesting a highly conserved structure. A high copy number of the trnM-CAT gene (4 copies) is noted, which may contribute to genomic rearrangement and adaptive evolution. Among the 476 RNA editing sites, hydrophilic-hydrophobic and hydrophobic-hydrophobic editing events are most common, accounting for 77.10%. Negative selection predominates among most genes (Ka/Ks < 1), while a few genes (e.g., matR, nad3, rps1, rps12, and rps4) show signs of positive selection (Ka/Ks > 1), potentially conferring evolutionary advantages. Additionally, a significant A/T bias is observed at the third codon position. Phylogenomic analysis supports the APG IV classification, with no evidence of horizontal gene transfer.

Discussion

This mitochondrial genome offers valuable insights into the adaptive mechanisms and evolutionary processes of Garuga forrestii. It enhances our understanding of the species' biogeography in tropical Southeast Asia and Southwest China, providing key information on the evolutionary history of this genus.

Mitochondrial Genome Insights into Evolution and Gene Regulation in Phragmites australis

As a globally distributed perennial Gramineae, Phragmites australis can adapt to harsh ecological environments and has significant economic and environmental values. Here, we performed a complete assembly and annotation of the mitogenome of P. australis using genomic data from the PacBio and BGI platforms. The P. australis mitogenome is a multibranched structure of 501,134 bp, divided into two circular chromosomes of 325,493 bp and 175,641 bp, respectively. A sequence-simplified succinate dehydrogenase 4 gene was identified in this mitogenome, which is often translocated to the nuclear genome in the mitogenomes of gramineous species. We also identified tissue-specific mitochondrial differentially expressed genes using RNAseq data, providing new insights into understanding energy allocation and gene regulatory strategies in the long-term adaptive evolution of P. australis mitochondria. In addition, we studied the mitogenome features of P. australis in more detail, including repetitive sequences, gene Ka/Ks analyses, codon preferences, intracellular gene transfer, RNA editing, and multispecies phylogenetic analyses. Our results provide an essential molecular resource for understanding the genetic characterisation of the mitogenome of P. australis and provide a research basis for population genetics and species evolution in Arundiaceae.

Identification of RNA Editing Sites Reveals Functional Modifications with the Addition of Methionine to the Daily Rations of Yaks

Methionine is an amino acid necessary for the growth and development of all animals. Glutathione produced during methionine metabolism can reduce damage to cells caused by oxidative stress. Supplementing restricted amino acids in animals by scientific means will be beneficial to protein synthesis, which will affect the growth and development of animals and will bring huge economic benefits when applied to actual production and life. In this study, we collected three muscle tissues from 24 male Maiwa yaks, which were fattened for three months with different methionine concentrations in their diet. RNA-seq was performed to obtain expression reads. A total of 1116 editing sites were identified by at least two software; the editing site types were mainly T-to-C and A-to-G mutations. We found two significant RNA editing sites presenting high-risk editing types. One was located on the MSRA gene that regulates the reduction of methionine, and the other can make changes to the properties of encoded proteins. This provides further understanding of the mechanism of yak muscle tissue and regulation of gene expression after the addition of methionine to daily rations.

metsDB: a knowledgebase of cancer metastasis at bulk, single-cell and spatial levels

Cancer metastasis, the process by which tumour cells migrate and colonize distant organs from a primary site, is responsible for the majority of cancer-related deaths. Understanding the cellular and molecular mechanisms underlying this complex process is essential for developing effective metastasis prevention and therapy strategies. To this end, we systematically analysed 1786 bulk tissue samples from 13 cancer types, 988 463 single cells from 17 cancer types, and 40 252 spots from 45 spatial slides across 10 cancer types. The results of these analyses are compiled in the metsDB database, accessible at https://relab.xidian.edu.cn/metsDB/. This database provides insights into alterations in cell constitutions, cell relationships, biological pathways, molecular biomarkers, and drug responses during cancer metastasis at bulk, single-cell, and spatial levels. Users can perform cell or gene searches to obtain multi-view and multi-scale metastasis-related data. This comprehensive resource is invaluable for understanding the metastasis process and for designing molecular therapies.

REDIportal: toward an integrated view of the A-to-I editing

A-to-I RNA editing is the most common non-transient epitranscriptome modification. It plays several roles in human physiology and has been linked to several disorders. Large-scale deep transcriptome sequencing has fostered the characterization of A-to-I editing at the single nucleotide level and the development of dedicated computational resources. REDIportal is a unique and specialized database collecting ∼16 million of putative A-to-I editing sites designed to face the current challenges of epitranscriptomics. Its running version has been enriched with sites from the TCGA project (using data from 31 studies). REDIportal provides an accurate, sustainable and accessible tool enriched with interconnections with widespread ELIXIR core resources such as Ensembl, RNAcentral, UniProt and PRIDE. Additionally, REDIportal now includes information regarding RNA editing in putative double-stranded RNAs, relevant for the immune-related roles of editing, as well as an extended catalog of recoding events. Finally, we report a reliability score per site calculated using a deep learning model trained using a huge collection of positive and negative instances. REDIportal is available at http://srv00.recas.ba.infn.it/atlas/.

RPS 2.0: an updated database of RNAs involved in liquid-liquid phase separation

Liquid-liquid phase separation (LLPS) is a crucial process for the formation of biomolecular condensates such as coacervate droplets, P-bodies and stress granules, which play critical roles in many physiological and pathological processes. Increasing studies have shown that not only proteins but also RNAs play a critical role in LLPS. To host LLPS-associated RNAs, we previously developed a database named 'RPS' in 2021. In this study, we present an updated version RPS 2.0 (https://rps.renlab.cn/) to incorporate the newly generated data and to host new LLPS-associated RNAs driven by post-transcriptional regulatory mechanisms. Currently, RPS 2.0 hosts 171 301 entries of LLPS-associated RNAs in 24 different biomolecular condensates with four evidence types, including 'Reviewed', 'High-throughput (LLPS enrichment)', 'High-throughput (LLPS perturbation)' and 'Predicted', and five event types, including 'Expression', 'APA', 'AS', 'A-to-I' and 'Modification'. Additionally, extensive annotations of LLPS-associated RNAs are provided in RPS 2.0, including RNA sequence and structure features, RNA-protein/RNA-RNA interactions, RNA modifications, as well as diseases related annotations. We expect that RPS 2.0 will further promote research of LLPS-associated RNAs and deepen our understanding of the biological functions and regulatory mechanisms of LLPS.

ncRNA Editing: Functional Characterization and Computational Resources

Non-coding RNAs (ncRNAs) play crucial roles in gene expression regulation, translation, and disease development, including cancer. They are classified by size in short and long non-coding RNAs. This chapter focuses on the functional implications of adenosine-to-inosine (A-to-I) RNA editing in both short (e.g., miRNAs) and long ncRNAs. RNA editing dynamically alters the sequence and structure of primary transcripts, impacting ncRNA biogenesis and function. Notable findings include the role of miRNA editing in promoting glioblastoma invasiveness, characterizing RNA editing hotspots across cancers, and its implications in thyroid cancer and ischemia. This chapter also highlights bioinformatics resources and next-generation sequencing (NGS) technologies that enable comprehensive ncRNAome studies and genome-wide RNA editing detection. Dysregulation of RNA editing machinery has been linked to various human diseases, emphasizing the potential of RNA editing as a biomarker and therapeutic target. This overview integrates current knowledge and computational tools for studying ncRNA editing, providing insights into its biological significance and clinical applications.

An improved transcriptome annotation reveals asymmetric expression and distinct regulation patterns in allotetraploid common carp

In allotetraploid common carp, protein-coding homoeologs presented divergent expression levels between the two subgenomes. However, whether subgenome dominance occurs in other transcriptional and post-transcriptional events remains unknown. Using Illumina RNA sequencing and PacBio full-length sequencing, we refined the common carp transcriptome annotation and explored differences in four transcriptional and post-transcriptional events between the two subgenomes. The results revealed that the B subgenome presented more alternative splicing events, as did lncRNAs and circRNAs. However, the expression levels, tissue specificity, sequence features, and functions of lncRNAs and circRNAs did not significantly differ between the two subgenomes, suggesting a common regulatory mechanism shared by the two subgenomes. Furthermore, both the number and base substitution frequency of RNA editing events were greater in the B subgenome. Functional analyses of these transcriptional events also revealed subgenome bias. Genes that undergo alternative splicing in the A subgenome participate in more biological processes, and lncRNA targets show a preference between subgenomes. CircRNA host genes in the B subgenome were associated with more biological functions, and RNA editing preferentially occurred in noncoding regions or led to nonsynonymous mutations in the B subgenome. Taken together, the refined transcriptome annotation revealed complicated and imbalanced expression strategies in allotetraploid common carp.

LoDEI: a robust and sensitive tool to detect transcriptome-wide differential A-to-I editing in RNA-seq data

RNA editing is a highly conserved process. Adenosine deaminase acting on RNA (ADAR) mediated deamination of adenosine (A-to-I editing) is associated with human disease and immune checkpoint control. Functional implications of A-to-I editing are currently of broad interest to academic and industrial research as underscored by the fast-growing number of clinical studies applying base editors as therapeutic tools. Analyzing the dynamics of A-to-I editing, in a biological or therapeutic context, requires the sensitive detection of differential A-to-I editing, a currently unmet need. We introduce the local differential editing index (LoDEI) to detect differential A-to-I editing in RNA-seq datasets using a sliding-window approach coupled with an empirical q value calculation that detects more A-to-I editing sites at the same false-discovery rate compared to existing methods. LoDEI is validated on known and novel datasets revealing that the oncogene MYCN increases and that a specific small non-coding RNA reduces A-to-I editing.

Repeat-mediated recombination results in Complex DNA structure of the mitochondrial genome of Trachelospermum jasminoides

BACKGROUND

Trachelospermum jasminoides has medicinal and ornamental value and is widely distributed in China. Although the chloroplast genome has been documented, the mitochondrial genome has not yet been studied.

RESULTS

The mitochondrial genome of T. jasminoides was assembled and functionally annotated using Illumina and nanopore reads. The mitochondrial genome comprises a master circular molecular structure of 605,764 bp and encodes 65 genes: 39 protein-coding genes, 23 transfer RNA (tRNA) genes and 3 ribosomal RNA genes. In addition to the single circular conformation, we found many alternative conformations of the T. jasminoides mitochondrial genome mediated by 42 repetitive sequences. Six repetitive sequences (DRS01-DRS06) were supported by nanopore long reads, polymerase chain reaction (PCR) amplifications, and Sanger sequencing of the PCR products. Eleven homologous fragments were identified by comparing the mitochondrial and chloroplast genome sequences, including three complete tRNA genes. Moreover, 531 edited RNA sites were identified in the protein-coding sequences based on RNA sequencing data, with nad4 having the highest number of sites (54).

CONCLUSION

To our knowledge, this is the first description of the mitochondrial genome of T. jasminoides. Our results demonstrate the existence of multiple conformations. These findings lay a foundation for understanding the genetics and evolutionary dynamics of Apocynaceae.

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.

ADAR1 promotes cisplatin resistance in intrahepatic cholangiocarcinoma by regulating BRCA2 expression through A-to-I editing manner

Aberrant A-to-I RNA editing, mediated by ADAR1 has been found to be associated with increased tumourigenesis and the development of chemotherapy resistance in various types of cancer. Intrahepatic cholangiocarcinoma (iCCA) is a highly aggressive malignancy with a poor prognosis, and overcoming chemotherapy resistance poses a significant clinical challenge. This study aimed to clarify the roles of ADAR1 in tumour resistance to cisplatin in iCCA. We discovered that ADAR1 expression is elevated in iCCA patients, particularly in those resistant to cisplatin, and associated with poor clinical outcomes. Downregulation of ADAR1 can increase the sensitivity of iCCA cells to cisplatin treatment, whereas its overexpression has the inverse effect. By integrating RNA sequencing and Sanger sequencing, we identified BRCA2, a critical DNA damage repair gene, as a downstream target of ADAR1 in iCCA. ADAR1 mediates the A-to-I editing in BRCA2 3'UTR, inhibiting miR-3157-5p binding, consequently increasing BRCA2 mRNA and protein levels. Furthermore, ADAR1 enhances cellular DNA damage repair ability and facilitates cisplatin resistance in iCCA cells. Combining ADAR1 targeting with cisplatin treatment markedly enhances the anticancer efficacy of cisplatin. In conclusion, ADAR1 promotes tumour progression and cisplatin resistance of iCCA. ADAR1 targeting could inform the development of innovative combination therapies for iCCA.

Complete mitochondrial genome assembly of Juglans regia unveiled its molecular characteristics, genome evolution, and phylogenetic implications

BACKGROUND

The Persian walnut (Juglans regia), an economically vital species within the Juglandaceae family, has seen its mitochondrial genome sequenced and assembled in the current study using advanced Illumina and Nanopore sequencing technology.

RESULTS

The 1,007,576 bp mitogenome of J. regia consisted of three circular chromosomes with a 44.52% GC content encoding 39 PCGs, 47 tRNA, and five rRNA genes. Extensive repetitive sequences, including 320 SSRs, 512 interspersed, and 83 tandem repeats, were identified, contributing to genomic complexity. The protein-coding sequences (PCGs) favored A/T-ending codons, and the codon usage bias was primarily shaped by selective pressure. Intracellular gene transfer occurred among the mitogenome, chloroplast, and nuclear genomes. Comparative genomic analysis unveiled abundant structure and sequence variation among J. regia and related species. The results of selective pressure analysis indicated that most PCGs underwent purifying selection, whereas the atp4 and ccmB genes had experienced positive selection between many species pairs. In addition, the phylogenetic examination, grounded in mitochondrial genome data, precisely delineated the evolutionary and taxonomic relationships of J. regia and its relatives. We identified a total of 539 RNA editing sites, among which 288 were corroborated by transcriptome sequencing data. Furthermore, expression profiling under temperature stress highlighted the complex regulation pattern of 28 differently expressed PCGs, wherein NADH dehydrogenase and ATP synthase genes might be critical in the mitochondria response to cold stress.

CONCLUSIONS

Our results provided valuable molecular resources for understanding the genetic characteristics of J. regia and offered novel perspectives for population genetics and evolutionary studies in Juglans and related woody species.

Complete organelle genomes of the threatened aquatic species Scheuchzeria palustris (Scheuchzeriaceae): Insights into adaptation and phylogenomic placement

Scheuchzeria palustris, the only species in the Scheuchzeriaceae family, plays a crucial role in methane production and transportation, influencing the global carbon cycle and maintaining ecosystem stability. However, it is now threatened by human activities and global warming. In this study, we generated new organelle genomes for S. palustris, with the plastome (pt) measuring 158,573 bp and the mitogenome (mt) measuring 420,724 bp. We predicted 296 RNA editing sites in mt protein-coding genes (PCGs) and 142 in pt-PCGs. Notably, abundant RNA editing sites in pt-PCGs likely originated from horizontal gene transfer between the plastome and mitogenome. Additionally, we identified positive selection signals in four mt-PCGs (atp4, ccmB, nad3, and sdh4) and one pt-PCG (rps7), which may contribute to the adaptation of S. palustris to low-temperature and high-altitude environments. Furthermore, we identified 35 mitochondrial plastid DNA (MTPT) segments totaling 58,479 bp, attributed to dispersed repeats near most MTPT. Phylogenetic trees reconstructed from mt- and pt-PCGs showed topologies consistent with the APG IV system. However, the conflicting position of S. palustris can be explained by significant differences in the substitution rates of its mt- and pt-PCGs (p < .001). In conclusion, our study provides vital genomic resources to support future conservation efforts and explores the adaptation mechanisms of S. palustris.

Prion diseases disrupt glutamate/glutamine metabolism in skeletal muscle

In prion diseases (PrDs), aggregates of misfolded prion protein (PrPSc) accumulate not only in the brain but also in extraneural organs. This raises the question whether prion-specific pathologies arise also extraneurally. Here we sequenced mRNA transcripts in skeletal muscle, spleen and blood of prion-inoculated mice at eight timepoints during disease progression. We detected gene-expression changes in all three organs, with skeletal muscle showing the most consistent alterations. The glutamate-ammonia ligase (GLUL) gene exhibited uniform upregulation in skeletal muscles of mice infected with three distinct scrapie prion strains (RML, ME7, and 22L) and in victims of human sporadic Creutzfeldt-Jakob disease. GLUL dysregulation was accompanied by changes in glutamate/glutamine metabolism, leading to reduced glutamate levels in skeletal muscle. None of these changes were observed in skeletal muscle of humans with amyotrophic lateral sclerosis, Alzheimer's disease, or dementia with Lewy bodies, suggesting that they are specific to prion diseases. These findings reveal an unexpected metabolic dimension of prion infections and point to a potential role for GLUL dysregulation in the glutamate/glutamine metabolism in prion-affected skeletal muscle.

Dissecting transposable elements and endogenous retroviruses upregulation by HDAC inhibitors in leiomyosarcoma cells: Implications for the interferon response

Transposable elements (TEs) are of interest as immunomodulators for cancer therapies. TEs can fold into dsRNAs that trigger the interferon response. Here, we investigated the effect of different HDAC inhibitors (HDACIs) on the expression of TEs in leiomyosarcoma cells. Our data show that endogenous retroviruses (ERVs), especially ERV1 elements, are upregulated after treatment with HDAC1/2/3-specific inhibitors. Surprisingly, the interferon response was not activated. We observed an increase in A-to-I editing of upregulated ERV1. This could have an impact on the stability of dsRNAs and the activation of the interferon response. We also found that H3K27ac levels are increased in the LTR12 subfamilies, which could be regulatory elements controlling the expression of proapoptotic genes such as TNFRSF10B. In summary, we provide a detailed characterization of TEs modulation in response to HDACIs and suggest the use of HDACIs in combination with ADAR inhibitors to induce cell death and support immunotherapy in cancer.

Impact of vaccination on SARS-CoV-2 evolution and immune escape variants

Vaccines and host genetic factors can influence the SARS-CoV-2 evolution and emergence of new variants. Even vaccinated cases get affected as virus continues to evolve, raising concerns about vaccine efficacy and the emergence of immune escape variants. Here, we have analyzed 2295 whole-genome sequences of SARS-CoV-2 collected from vaccinated and unvaccinated cases to evaluate the impact of vaccines on virus diversity within hosts. Our comparative analysis revealed a significant higher incidence of intra-host single nucleotides variants (iSNVs) in vaccinated cases compared to unvaccinated ones (p value<0.0001). Furthermore, we have found that specific mutational processes, including APOBEC (C > T) mediated and ADAR1 (A > G) mediated mutations, were found more prevalent in vaccinated cases. Vaccinated cases exhibited higher accumulation of nonsynonymous mutation than unvaccinated cases. Fixed iSNVs were predominantly located in the ORF1ab and spike genes, several key omicron defining immune escape variants S477N, Q493R, Q498R, Y505H, L452R, and N501Y were identified in the RBD domain of spike gene in vaccinated cases. Our findings suggest that vaccine plays an important role in the evolution of the virus genome. The virus genome acquires random mutations due to error-prone replication of the virus, host modification through APOBEC and ADAR1 mediated editing mechanism, and oxidative stress. These mutations become fixed in the viral population due to the selective pressure imposed by vaccination.

Comprehensive analysis of the Lycopodium japonicum mitogenome reveals abundant tRNA genes and cis-spliced introns in Lycopodiaceae species

Lycophytes and ferns represent one of the earliest-diverging lineages of vascular plants, with the Lycopodiaceae family constituting the basal clade among lycophytes. In this research, we successfully assembled and annotated the complete Lycopodium japonicum Thunb. (L. japonicum) mitochondrial genome (mitogenome) utilizing PacBio HiFi sequencing data, resulting in a single circular molecule with a size of 454,458 bp. 64 unique genes were annotated altogether, including 34 protein-coding genes, 27 tRNAs and 3 rRNAs. It also contains 32 group II introns, all of which undergo cis-splicing. We identified 195 simple sequence repeats, 1,948 dispersed repeats, and 92 tandem repeats in the L. japonicum mitogenome. Collinear analysis indicated that the mitogenomes of Lycopodiaceae are remarkably conserved compared to those of other vascular plants. We totally identified 326 RNA editing sites in 31 unique protein-coding genes with 299 sites converting cytosine to uracil and 27 sites the reverse. Notably, the L. japonicum mitogenome has small amounts foreign DNA from plastid or nuclear origin, accounting for only 2.81% of the mitogenome. The maximum likelihood phylogenetic analysis based on 23 diverse land plant mitogenomes and plastid genomes supports the basal position of lycophytes within vascular plants and they form a sister clade to all other vascular lineages, which is consistent with the PPG I classification system. As the first reported mitogenome of Lycopodioideae subfamily, this study enriches our understanding of Lycopodium mitogenomes, and sets the stage for future research on mitochondrial diversity and evolution within the lycophytes and ferns.

RNA Editing by ADAR Adenosine Deaminases in the Cell Models of CAG Repeat Expansion Diseases: Significant Effect of Differentiation from Stem Cells into Brain Organoids in the Absence of Substantial Influence of CAG Repeats on the Level of Editing

Expansion of CAG repeats in certain genes is a known cause of several neurodegenerative diseases, but exact mechanism behind this is not yet fully understood. It is believed that the double-stranded RNA regions formed by CAG repeats could be harmful to the cell. This study aimed to test the hypothesis that these RNA regions might potentially interfere with ADAR RNA editing enzymes, leading to the reduced A-to-I editing of RNA and activation of the interferon response. We studied induced pluripotent stem cells (iPSCs) derived from the patients with Huntington's disease or ataxia type 17, as well as midbrain organoids developed from these cells. A targeted panel for next-generation sequencing was used to assess editing in the specific RNA regions. Differentiation of iPSCs into brain organoids led to increase in the ADAR2 gene expression and decrease in the expression of protein inhibitors of RNA editing. As a result, there was increase in the editing of specific ADAR2 substrates, which allowed identification of differential substrates of ADAR isoforms. However, comparison of the pathology and control groups did not show differences in the editing levels among the iPSCs. Additionally, brain organoids with 42-46 CAG repeats did not exhibit global changes. On the other hand, brain organoids with the highest number of CAG repeats in the huntingtin gene (76) showed significant decrease in the level of RNA editing of specific transcripts, potentially involving ADAR1. Notably, editing of the long non-coding RNA PWAR5 was nearly absent in this sample. It could be stated in conclusion that in most cultures with repeat expansion, the hypothesized effect on RNA editing was not confirmed.

Insights into ADAR gene complement, expression patterns, and RNA editing landscape in Chlamys farreri

Adenosine Deaminases Acting on RNA (ADARs) are evolutionarily conserved enzymes known to convert adenosine to inosine in double-stranded RNAs and participate in host-virus interactions. Conducting a meta-analysis of available transcriptome data, we identified and characterised eight ADAR transcripts in Chlamys farreri, a farmed marine scallop susceptible to Acute viral necrosis virus (AVNV) infections and mortality outbreaks. Accordingly, we identified six ADAR genes in the Zhikong scallop genome, revised previous gene annotations, and traced alternative splicing variants. In detail, each ADAR gene encodes a unique combination of functional domains, always including the Adenosine deaminase domain, RNA binding domains and, in one case, two copies of a Z-DNA binding domain. After phylogenetic analysis, five C. farreri ADARs clustered in the ADAR1 clade along with sequences from diverse animal phyla. Gene expression analysis indicated CF051320 as the most expressed ADAR, especially in the eye and male gonad. The other four ADAR1 genes and one ADAR2 gene exhibited variable expression levels, with CF105370 and CF051320 significantly increasing during early scallop development. ADAR-mediated single-base editing, evaluated across adult C. farreri tissues and developmental stages, was mainly detectable in intergenic regions (83 % and 85 %, respectively). Overall, the expression patterns of the six ADAR genes together with the editing and hyper-editing values computed on scallops RNA-seq samples support the adaptive value of ADAR1-mediated editing, particularly in the pre-settling larval stages.

Programmable DNA pyrimidine base editing via engineered uracil-DNA glycosylase

DNA base editing technologies predominantly utilize engineered deaminases, limiting their ability to edit thymine and guanine directly. In this study, we successfully achieve base editing of both cytidine and thymine by leveraging the translesion DNA synthesis pathway through the engineering of uracil-DNA glycosylase (UNG). Employing structure-based rational design, exploration of homologous proteins, and mutation screening, we identify a Deinococcus radiodurans UNG mutant capable of effectively editing thymine. When fused with the nickase Cas9, the engineered DrUNG protein facilitates efficient thymine base editing at endogenous sites, achieving editing efficiencies up to 55% without enrichment and exhibiting minimal cellular toxicity. This thymine base editor (TBE) exhibits high editing specificity and significantly restores IDUA enzyme activity in cells derived from patients with Hurler syndrome. TBEs represent efficient, specific, and low-toxicity approaches to base editing with potential applications in treating relevant diseases.

Bioinformatics for Inosine: Tools and Approaches to Trace This Elusive RNA Modification

Inosine is a nucleotide resulting from the deamination of adenosine in RNA. This chemical modification process, known as RNA editing, is typically mediated by a family of double-stranded RNA binding proteins named Adenosine Deaminase Acting on dsRNA (ADAR). While the presence of ADAR orthologs has been traced throughout the evolution of metazoans, the existence and extension of RNA editing have been characterized in a more limited number of animals so far. Undoubtedly, ADAR-mediated RNA editing plays a vital role in physiology, organismal development and disease, making the understanding of the evolutionary conservation of this phenomenon pivotal to a deep characterization of relevant biological processes. However, the lack of direct high-throughput methods to reveal RNA modifications at single nucleotide resolution limited an extended investigation of RNA editing. Nowadays, these methods have been developed, and appropriate bioinformatic pipelines are required to fully exploit this data, which can complement existing approaches to detect ADAR editing. Here, we review the current literature on the "bioinformatics for inosine" subject and we discuss future research avenues in the field.

Highly active repeat-mediated recombination in the mitogenome of the aquatic grass Hygroryza aristata

BACKGROUND

Floating bamboo (Hygroryza aristata) is an endangered species with a narrow native distribution and is renowned for its unique aesthetic qualities, which holds significant ecological and ornamental value. However, the lack of genetic information research, with only one complete plastome available, significantly hampers conservation efforts and further research for this species.

RESULTS

In this research, we sequenced and assembled the organelle genomes of floating bamboo, including the mitogenome (587,847 bp) and plastome (135,675 bp). The mitogenome can recombine into various configurations, which are mediated by 25 repeat pairs (13 SRs, 6 MRs, 1 LR, and 5 CRs). LR1 and SR5 are particularly notable as they have the ability to combine with other contigs, forming complex repeat units that facilitate further homologous recombination. The rate of homologous recombination varies significantly among species, yet there is still a pronounced positive correlation observed between the length of these repeat pairs and the rate of recombination they mediate. The mitogenome integrates seven intact protein-coding genes from the chloroplast. The codon usage patterns in both organelles are similar, with a noticeable bias towards C and T on the third codon. The gene map of Poales shows the entire loss of rpl6, succinate dehydrogenase subunits (sdh3 and sdh4). Additionally, the BOP clade retained more variable genes compared to the PACMAD clade.

CONCLUSIONS

We provided a high-quality and well-annotated mitogenome for floating bamboo and demonstrated the presence of diverse configurations. Our study has revealed the correlation between repeat length and their corresponding recombination rate despite variations among species. Although the mitogenome can potentially exist in the form of a unicircular in vivo, this occurrence is rare and may not be stable.

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.

A-to-I Editing Is Subtype-Specific in Non-Hodgkin Lymphomas

Cancer is a complex and heterogeneous disease, in which a number of genetic and epigenetic changes occur in tumor onset and progression. Recent studies indicate that changes at the RNA level are also involved in tumorigenesis, such as adenosine-to-inosine (A-to-I) RNA editing. Here, we systematically investigate transcriptome-wide A-to-I editing events in a large number of samples from Non-Hodgkin lymphomas (NHLs). Using a computational pipeline that determines significant differences in editing level between NHL and normal samples at known A-to-I editing sites, we identify a number of differentially edited editing sites between NHL subtypes and normal samples. Most of the differentially edited sites are located in non-coding regions, and many such sites show a strong correlation between gene expression level and editing efficiency, indicating that RNA editing might have direct consequences for the cancer cell's aberrant gene regulation status in these cases. Moreover, we establish a strong link between RNA editing and NHL by demonstrating that NHL and normal samples and even NHL subtypes can be distinguished based on genome-wide RNA editing profiles alone. Our study establishes a strong link between RNA editing, cancer and aberrant gene regulation in NHL.

Assembly and comparative analysis of the first complete mitochondrial genome of the invasive water hyacinth, Eichhornia crassipes

Eichhornia crassipes is an aquatic plant in tropical and subtropical regions, renowned for its notorious invasive tendencies. In this study, we assembled the complete mitogenome of E. crassipes into a single circle molecule of 397,361 bp. The mitogenome has 58 unique genes, including 37 protein-coding genes (PCGs), 18 tRNA genes, three rRNA genes, and 47 % GC content. Sixteen (6.93 %) homologous fragments, ranging from 31 bp to 8548 bp, were identified, indicating the transfer of genetic material from chloroplasts to mitochondria. In addition, we detected positive selection in six PCGs (ccmB, ccmC, ccmFC, nad3, nad4 and sdh4), along with the identification of 782 RNA editing sites across 37 mt-PCGs. These findings suggest a potential contribution to the robust adaptation of this invasive plant to the stressful environment. Lastly, we inferred that phylogenetic conflicts of E. crassipes between the plastome and mitogenome may be attributed to the difference in nucleotide substitution rates between the two organelle genomes. In conclusion, our study provided vital genomic resources for further understanding the invasive mechanism of this species and exploring the dynamic evolution of mitogenomes within the monocot clade.

The complete mitochondrial genome of Pontederia crassipes: using HiFi reads to investigate genome recombination and gene transfer from chloroplast genome

Water hyacinth (Pontederia crassipes Mart.) is a monocotyledonous aquatic plant renowned for its rapid growth, extensive proliferation, biological invasiveness, and ecological resilience to variations in pH, nutrients, and temperature. The International Union for Conservation of Nature (IUCN) has listed P. crassipes among the top 100 invasive species. However, comprehensive genomic information, particularly concerning its mitochondrial genome (mitogenome), remains surprisingly limited. In this study, the complete mitogenome of P. crassipes was analyzed using bioinformatics approaches. The mitogenome is 399,263 bp long and contains 38 protein-coding genes (PCGs), 24 tRNA genes, and 3 rRNA genes. Sequence analysis revealed that the complete mitogenome of the species contains 3,289 dispersed repeats, and 765 RNA editing sites in protein-coding genes. The P. crassipes mitogenome possessed un-conserved structures, including extensive sequence transfer between its chloroplasts and mitochondria. Our study on the mitogenome of P. crassipes offers critical insights into its evolutionary patterns and phylogenetic relationships with related taxa. This research enhances our understanding of this invasive species, known for its significant biomass and rapid overgrowth in aquatic environments.

Hyphal editing of the conserved premature stop codon in CHE1 is stimulated by oxidative stress in Fusarium graminearum

Although genome-wide A-to-I editing mediated by adenosine-deaminase-acting-on-tRNA (ADAT) occurs during sexual reproduction in the presence of stage-specific cofactors, RNA editing is not known to occur during vegetative growth in filamentous fungi. Here we identified 33 A-to-I RNA editing events in vegetative hyphae of Fusarium graminearum and functionally characterized one conserved hyphal-editing site. Similar to ADAT-mediated editing during sexual reproduction, majority of hyphal-editing sites are in coding sequences and nonsynonymous, and have strong preference for U at -1 position and hairpin loops. Editing at TA437G, one of the hyphal-specific editing sites, is a premature stop codon correction (PSC) event that enables CHE1 gene to encode a full-length zinc fingertranscription factor. Manual annotations showed that this PSC site is conserved in CHE1 orthologs from closely-related Fusarium species. Whereas the che1 deletion and CHE1TAA (G438 to A) mutants had no detectable phenotype, the CHE1TGG (A437 to G) mutant was defective in hyphal growth, conidiation, sexual reproduction, and plant infection. However, the CHE1TGG mutant was increased in tolerance against oxidative stress and editing of TA437G in CHE1 was stimulated by H2O2 treatment in F. graminearum. These results indicate that fixation of the premature stop codon in CHE1 has a fitness cost on normal hyphal growth and reproduction but provides a benefit to tolerance against oxidative stress. Taken together, A-to-I editing events, although rare (not genome-wide), occur during vegetative growth and editing in CHE1 plays a role in response to oxidative stress in F. graminearum and likely in other fungal pathogens.

Elucidation of the Epitranscriptomic RNA Modification Landscape of Chikungunya Virus

The genomes of positive-sense (+) single-stranded RNA (ssRNA) viruses are believed to be subjected to a wide range of RNA modifications. In this study, we focused on the chikungunya virus (CHIKV) as a model (+) ssRNA virus to study the landscape of viral RNA modification in infected human cells. Among the 32 distinct RNA modifications analysed by mass spectrometry, inosine was found enriched in the genomic CHIKV RNA. However, orthogonal validation by Illumina RNA-seq analyses did not identify any inosine modification along the CHIKV RNA genome. Moreover, CHIKV infection did not alter the expression of ADAR1 isoforms, the enzymes that catalyse the adenosine to inosine conversion. Together, this study highlights the importance of a multidisciplinary approach to assess the presence of RNA modifications in viral RNA genomes.

Comprehensive analysis of the mitochondrial genome of Rehmannia glutinosa: insights into repeat-mediated recombinations and RNA editing-induced stop codon acquisition

Rehmannia glutinosa is an economically significant medicinal plant. Yet, the structure and sequence of its mitochondrial genome has not been published, which plays a crucial role in evolutionary analysis and regulating respiratory-related macromolecule synthesis. In this study, the R. glutinosa mitogenome was sequenced employing a combination of Illumina short reads and Nanopore long reads, with subsequent assembly using a hybrid strategy. We found that the predominant configuration of the R. glutinosa mitogenome comprises two circular chromosomes. The primary structure of the mitogenome encompasses two mitochondrial chromosomes corresponding to the two major configurations, Mac1-1 and Mac1-2. The R. glutinosa mitogenome encoded an angiosperm-typical set of 24 core genes, nine variable genes, three rRNA genes, and 15 tRNA genes. A phylogenetic analysis using the 16 shared protein-coding genes (PCG) yielded a tree consistent with the phylogeny of Lamiales species and two outgroup taxa. Mapping RNA-seq data to the coding sequences (CDS) of the PCGs revealed 507 C-to-U RNA editing sites across 31 PCGs of the R. glutinosa mitogenome. Furthermore, one start codon (nad4L) and two stop codons (rpl10 and atp6) were identified as products of RNA editing events in the R. glutinosa mitogenome.

Unveiling the A-to-I mRNA editing machinery and its regulation and evolution in fungi

A-to-I mRNA editing in animals is mediated by ADARs, but the mechanism underlying sexual stage-specific A-to-I mRNA editing in fungi remains unknown. Here, we show that the eukaryotic tRNA-specific heterodimeric deaminase FgTad2-FgTad3 is responsible for A-to-I mRNA editing in Fusarium graminearum. This editing capacity relies on the interaction between FgTad3 and a sexual stage-specific protein called Ame1. Although Ame1 orthologs are widely distributed in fungi, the interaction originates in Sordariomycetes. We have identified key residues responsible for the FgTad3-Ame1 interaction. The expression and activity of FgTad2-FgTad3 are regulated through alternative promoters, alternative translation initiation, and post-translational modifications. Our study demonstrates that the FgTad2-FgTad3-Ame1 complex can efficiently edit mRNA in yeasts, bacteria, and human cells, with important implications for the development of base editors in therapy and agriculture. Overall, this study uncovers mechanisms, regulation, and evolution of RNA editing in fungi, highlighting the role of protein-protein interactions in modulating deaminase function.

txtools: an R package facilitating analysis of RNA modifications, structures, and interactions

We present txtools, an R package that enables the processing, analysis, and visualization of RNA-seq data at the nucleotide-level resolution, seamlessly integrating alignments to the genome with transcriptomic representation. txtools' main inputs are BAM files and a transcriptome annotation, and the main output is a table, capturing mismatches, deletions, and the number of reads beginning and ending at each nucleotide in the transcriptomic space. txtools further facilitates downstream visualization and analyses. We showcase, using examples from the epitranscriptomic field, how a few calls to txtools functions can yield insightful and ready-to-publish results. txtools is of broad utility also in the context of structural mapping and RNA:protein interaction mapping. By providing a simple and intuitive framework, we believe that txtools will be a useful and convenient tool and pave the path for future discovery. txtools is available for installation from its GitHub repository at https://github.com/AngelCampos/txtools.

Characterization and Function Analysis of miRNA Editing during Fat Deposition in Chinese Indigenous Ningxiang Pigs

This study aimed to identify active miRNA editing sites during adipose development in Ningxiang pigs and analyze their characteristics and functions. Based on small RNA-seq data from the subcutaneous adipose tissues of Ningxiang pigs at four stages-30 days (piglet), 90 days (nursery), 150 days (early fattening), and 210 days (late fattening)-we constructed a developmental map of miRNA editing in the adipose tissues of Ningxiang pigs. A total of 505 miRNA editing sites were identified using the revised pipeline, with C-to-U editing types being the most prevalent, followed by U-to-C, A-to-G, and G-to-U. Importantly, these four types of miRNA editing exhibited base preferences. The number of editing sites showed obvious differences among age groups, with the highest occurrence of miRNA editing events observed at 90 days of age and the lowest at 150 days of age. A total of nine miRNA editing sites were identified in the miRNA seed region, with significant differences in editing levels (p < 0.05) located in ssc-miR-23a, ssc-miR-27a, ssc-miR-30b-5p, ssc-miR-15a, ssc-miR-497, ssc-miR-15b, and ssc-miR-425-5p, respectively. Target gene prediction and KEGG enrichment analyses indicated that the editing of miR-497 might potentially regulate fat deposition by inhibiting adipose synthesis via influencing target binding. These results provide new insights into the regulatory mechanism of pig fat deposition.

Restoration of virulence in the attenuated Candid#1 vaccine virus requires reversion at both positions 168 and 427 in the envelope glycoprotein GPC

Live-attenuated virus vaccines provide long-lived protection against viral disease but carry inherent risks of residual pathogenicity and genetic reversion. The live-attenuated Candid#1 vaccine was developed to protect Argentines against lethal infection by the Argentine hemorrhagic fever arenavirus, Junín virus. Despite its safety and efficacy in Phase III clinical study, the vaccine is not licensed in the US, in part due to concerns regarding the genetic stability of attenuation. Previous studies had identified a single F427I mutation in the transmembrane domain of the Candid#1 envelope glycoprotein GPC as the key determinant of attenuation, as well as the propensity of this mutation to revert upon passage in cell culture and neonatal mice. To ascertain the consequences of this reversion event, we introduced the I427F mutation into recombinant Candid#1 (I427F rCan) and investigated the effects in two validated small-animal models: in mice expressing the essential virus receptor (human transferrin receptor 1; huTfR1) and in the conventional guinea pig model. We report that I427F rCan displays only modest virulence in huTfR1 mice and appears attenuated in guinea pigs. Reversion at another attenuating locus in Candid#1 GPC (T168A) was also examined, and a similar pattern was observed. By contrast, virus bearing both revertant mutations (A168T+I427F rCan) approached the lethal virulence of the pathogenic Romero strain in huTfR1 mice. Virulence was less extreme in guinea pigs. Our findings suggest that genetic stabilization at both positions is required to minimize the likelihood of reversion to virulence in a second-generation Candid#1 vaccine.IMPORTANCELive-attenuated virus vaccines, such as measles/mumps/rubella and oral poliovirus, provide robust protection against disease but carry with them the risk of genetic reversion to the virulent form. Here, we analyze the genetics of reversion in the live-attenuated Candid#1 vaccine that is used to protect against Argentine hemorrhagic fever, an often-lethal disease caused by the Junín arenavirus. In two validated small-animal models, we find that restoration of virulence in recombinant Candid#1 viruses requires back-mutation at two positions specific to the Candid#1 envelope glycoprotein GPC, at positions 168 and 427. Viruses bearing only a single change showed only modest virulence. We discuss strategies to genetically harden Candid#1 GPC against these two reversion events in order to develop a safer second-generation Candid#1 vaccine virus.

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.

Assessing the complementary information from an increased number of biologically relevant features in liquid biopsy-derived RNA-Seq data

Liquid biopsy-derived RNA sequencing (lbRNA-seq) exhibits significant promise for clinic-oriented cancer diagnostics due to its non-invasiveness and ease of repeatability. Despite substantial advancements, obstacles like technical artefacts and process standardisation impede seamless clinical integration. Alongside addressing technical aspects such as normalising fluctuating low-input material and establishing a standardised clinical workflow, the lack of result validation using independent datasets remains a critical factor contributing to the often low reproducibility of liquid biopsy-detected biomarkers. Considering the outlined drawbacks, our objective was to establish a workflow/methodology characterised by: 1. Harness the rich diversity of biological features accessible through lbRNA-seq data, encompassing a holistic range of molecular and functional attributes. These components are seamlessly integrated via a Machine Learning-based Ensemble Classification framework, enabling a unified and comprehensive analysis of the intricate information encoded within the data. 2. Implementing and rigorously benchmarking intra-sample normalisation methods to heighten their relevance within clinical settings. 3. Thoroughly assessing its efficacy across independent test sets to ascertain its robustness and potential utility. Using ten datasets from several studies comprising three different sources of biological material, we first show that while the best-performing normalisation methods depend strongly on the dataset and coupled Machine Learning method, the rather simple Counts Per Million method is generally very robust, showing comparable performance to cross-sample methods. Subsequently, we demonstrate that the innovative biofeature types introduced in this study, such as the Fraction of Canonical Transcript, harbour complementary information. Consequently, their inclusion consistently enhances prediction power compared to models relying solely on gene expression-based biofeatures. Finally, we demonstrate that the workflow is robust on completely independent datasets, generally from different labs and/or different protocols. Taken together, the workflow presented here outperforms generally employed methods in prediction accuracy and may hold potential for clinical diagnostics application due to its specific design.

Geographical Expansion of Avian Metapneumovirus Subtype B: First Detection and Molecular Characterization of Avian Metapneumovirus Subtype B in US Poultry

Avian metapneumovirus (aMPV), classified within the Pneumoviridae family, wreaks havoc on poultry health. It typically causes upper respiratory tract and reproductive tract infections, mainly in turkeys, chickens, and ducks. Four subtypes of AMPV (A, B, C, D) and two unclassified subtypes have been identified, of which subtypes A and B are widely distributed across the world. In January 2024, an outbreak of severe respiratory disease occurred on turkey and chicken farms across different states in the US. Metagenomics sequencing of selected tissue and swab samples confirmed the presence of aMPV subtype B. Subsequently, all samples were screened using an aMPV subtype A and B multiplex real-time RT-PCR kit. Of the 221 farms, 124 (56%) were found to be positive for aMPV-B. All samples were negative for subtype A. Six whole genomes were assembled, five from turkeys and one from chickens; all six assembled genomes showed 99.29 to 99.98% nucleotide identity, indicating a clonal expansion event for aMPV-B within the country. In addition, all six sequences showed 97.74 to 98.58% nucleotide identity with previously reported subtype B sequences, e.g., VCO3/60616, Hungary/657/4, and BR/1890/E1/19. In comparison to these two reference strains, the study sequences showed unique 49-62 amino acid changes across the genome, with maximum changes in glycoprotein (G). One unique AA change from T (Threonine) to I (Isoleucine) at position 153 in G protein was reported only in the chicken aMPV sequence, which differentiated it from turkey sequences. The twelve unique AA changes along with change in polarity of the G protein may indicate that these unique changes played a role in the adaptation of this virus in the US poultry. This is the first documented report of aMPV subtype B in US poultry, highlighting the need for further investigations into its genotypic characterization, pathogenesis, and evolutionary dynamics.

Sexual stage-specific A-to-I mRNA editing is mediated by tRNA-editing enzymes in fungi

A-to-I RNA editing catalyzed by adenosine-deaminase-acting-on-RNA (ADARs) was assumed to be unique to metazoans because fungi and plants lack ADAR homologs. However, genome-wide messenger RNA (mRNA) editing was found to occur specifically during sexual reproduction in filamentous ascomycetes. Because systematic characterization of adenosine/cytosine deaminase genes has implicated the involvement of TAD2 and TAD3 orthologs in A-to-I editing, in this study, we used genetic and biochemical approaches to characterize the role of FgTAD2, an essential adenosine-deaminase-acting-on-tRNA (ADAT) gene, in mRNA editing in Fusarium graminearum. FgTAD2 had a sexual-stage-specific isoform and formed heterodimers with enzymatically inactive FgTAD3. Using a repeat-induced point (RIP) mutation approach, we identified 17 mutations in FgTAD2 that affected mRNA editing during sexual reproduction but had no effect on transfer RNA (tRNA) editing and vegetative growth. The functional importance of the H352Y and Q375*(nonsense) mutations in sexual reproduction and mRNA editing were confirmed by introducing specific point mutations into the endogenous FgTAD2 allele in the wild type. An in vitro assay was developed to show that FgTad2-His proteins purified from perithecia, but not from vegetative hyphae, had mRNA editing activities. Moreover, the H352Y mutation affected the enzymatic activity of FgTad2 to edit mRNA but had no effect on its ADAT activity. We also identified proteins co-purified with FgTad2-His by mass spectrometry analysis and found that two of them have the RNA recognition motif. Taken together, genetic and biochemical data from this study demonstrated that FgTad2, an ADAT, catalyzes A-to-I mRNA editing with the stage-specific isoform and cofactors during sexual reproduction in fungi.

Recommendations for detection, validation, and evaluation of RNA editing events in cardiovascular and neurological/neurodegenerative diseases

RNA editing, a common and potentially highly functional form of RNA modification, encompasses two different RNA modifications, namely adenosine to inosine (A-to-I) and cytidine to uridine (C-to-U) editing. As inosines are interpreted as guanosines by the cellular machinery, both A-to-I and C-to-U editing change the nucleotide sequence of the RNA. Editing events in coding sequences have the potential to change the amino acid sequence of proteins, whereas editing events in noncoding RNAs can, for example, affect microRNA target binding. With advancing RNA sequencing technology, more RNA editing events are being discovered, studied, and reported. However, RNA editing events are still often overlooked or discarded as sequence read quality defects. With this position paper, we aim to provide guidelines and recommendations for the detection, validation, and follow-up experiments to study RNA editing, taking examples from the fields of cardiovascular and brain disease. We discuss all steps, from sample collection, storage, and preparation, to different strategies for RNA sequencing and editing-sensitive data analysis strategies, to validation and follow-up experiments, as well as potential pitfalls and gaps in the available technologies. This paper may be used as an experimental guideline for RNA editing studies in any disease context.

Global characterization of RNA editing in genetic regulation of multiple ovarian cancer subtypes

RNA editing plays an extensive role in the initiation and progression of cancer. However, the overall profile and molecular functions of RNA editing in different ovarian cancer subtypes have not been fully characterized and elucidated. Here, we conducted a study on RNA editing in four cohorts of ovarian cancer subtypes through large-scale parallel reporting and bioinformatics analysis. Our findings revealed that RNA editing patterns exhibit subtype-specific characteristics within cancer subtypes. The expression pattern of ADAR and the number of differential editing sites varied under different conditions. CCOC and EOC exhibited significant editing deficiency, whereas HGSC and MOC displayed significant editing excess. The sites within the turquoise module of the coedited network also revealed their correlation with ovarian cancer. In addition, we identified an average of over 40,000 cis-edQTLs in the four subtypes. Finally, we explored the association between RNA editing and drug response, uncovering several potentially effective editing-drug pairs (EDP) and suggesting the conceivable utility of RNA editing sites as therapeutic targets for cancer treatment. Overall, our comprehensive study has identified and characterized RNA editing events in various subtypes of ovarian cancer, providing a new perspective for ovarian cancer research and facilitating the development of medical interventions and treatments.