Editing Sites Research Articles

Mitochondrial genome of Artemisia argyi L. suggested conserved mitochondrial protein-coding genes among genera Artemisia, Tanacetum and Chrysanthemum

BackgroundArtemisia argyi L., also known as mugwort, is a perennial herb whose leaves are commonly used as a source of traditional medicines. However, the evolution and structure of the mitochondrial genome (mitogenome) in A. argyi remain unclear. In this study, the mitogenome of A. argyi was assembled and characterized for the first time. ResultsThe mitogenome of A. argyi was a circular molecule of 229,354 bp. It encodes 56 genes, including 33 protein-coding genes (PCGs), 20 tRNA genes, and three rRNA genes, and three pseudogenes. Five trans-spliced introns were observed in three PCGs namely, nad1, nad2 and nad5. Repeat analysis identified 65 SSRs, 14 tandem repeats, and 167 dispersed repeats. The A. argyi mitogenome contains 12 plastid transfer sequences from 79 bp to 2552 bp. Five conserved MTPTs were identified in all 18 Asteraceae species. Comparison of mitogenome between A. argyi and one Artemisia specie and two Chrysanthemum species showed 14 conserved gene clusters. Phylogenetic analysis with organelle genomes of A. argyi and 18 other Anthemideae plants showed inconsistent phylogenetic trees, which implied that the evolutionary rates of PCGs and rrna genes derived from mitochondrion and plastid were incongruent. The Ka/Ks ratio of the 27 shared protein-coding genes in the 18 Anthemideae species are all less than 1 indicating that these genes were under the effect of purifying selection. Lastly, a total of 568 RNA editing sites in PCGs were further identified. The average editing frequency of non-synonymous changes was significantly higher than that of synonymous changes (one-sample Student’s t-test, p-values ≤ 0.05) in three tissues (root, leaf and stem). ConclusionsIn this study, the gene content, genome size, genome comparison, mitochondrial plastid sequences, dN/dS analysis of mitochondrial protein-coding genes, and RNA-editing events in A. argyi mitogenome were determined, providing insights into the phylogenetic relationships of Asteraceae plant.

Increased A-to-I RNA editing in atherosclerosis and cardiomyopathies.

Adenosine-to-inosine RNA editing is essential to prevent undesired immune activation. This diverse process alters the genetic content of the RNA and may recode proteins, change splice sites and miRNA targets, and mimic genomic mutations. Recent studies have associated or implicated aberrant editing with pathological conditions, including cancer, autoimmune diseases, and neurological and psychiatric conditions. RNA editing patterns in cardiovascular tissues have not been investigated systematically so far, and little is known about its potential role in cardiac diseases. Some hints suggest robust editing in this system, including the fact that ADARB1 (ADAR2), the main coding-sequence editor, is most highly expressed in these tissues. Here we characterized RNA editing in the heart and arteries and examined a contributory role to the development of atherosclerosis and two structural heart diseases -Ischemic and Dilated Cardiomyopathies. Analyzing hundreds of RNA-seq samples taken from the heart and arteries of cardiac patients and controls, we find that global editing, alongside inflammatory gene expression, is increased in patients with atherosclerosis, cardiomyopathies, and heart failure. We describe a single recoding editing site and suggest it as a target for focused research. This recoding editing site in the IGFBP7 gene is one of the only evolutionary conserved sites between mammals, and we found it exhibits consistently increased levels of editing in these patients. Our findings reveal that RNA editing is abundant in arteries and is elevated in several key cardiovascular conditions. They thus provide a roadmap for basic and translational research of RNA as a mediator of atherosclerosis and non-genetic cardiomyopathies.

Widespread RNA hypoediting in schizophrenia and its relevance to mitochondrial function.

RNA editing, the endogenous modification of nucleic acids, is known to be altered in genes with important neurological function in schizophrenia (SCZ). However, the global profile and molecular functions of disease-associated RNA editing remain unclear. Here, we analyzed RNA editing in postmortem brains of four SCZ cohorts and uncovered a significant and reproducible trend of hypoediting in patients of European descent. We report a set of SCZ-associated editing sites via WGCNA analysis, shared across cohorts. Using massively parallel reporter assays and bioinformatic analyses, we observed that differential 3' untranslated region (3'UTR) editing sites affecting host gene expression were enriched for mitochondrial processes. Furthermore, we characterized the impact of two recoding sites in the mitofusin 1 (MFN1) gene and showed their functional relevance to mitochondrial fusion and cellular apoptosis. Our study reveals a global reduction of editing in SCZ and a compelling link between editing and mitochondrial function in the disease.

The ADAR1 editome reveals drivers of editing-specificity for ADAR1-isoforms.

Adenosine deaminase acting on RNA ADAR1 promotes A-to-I conversion in double-stranded and structured RNAs. ADAR1 has two isoforms transcribed from different promoters: cytoplasmic ADAR1p150 is interferon-inducible while ADAR1p110 is constitutively expressed and primarily localized in the nucleus. Mutations in ADAR1 cause Aicardi - Goutières syndrome (AGS), a severe autoinflammatory disease associated with aberrant IFN production. In mice, deletion of ADAR1 or the p150 isoform leads to embryonic lethality driven by overexpression of interferon-stimulated genes. This phenotype is rescued by deletion of the cytoplasmic dsRNA-sensor MDA5 indicating that the p150 isoform is indispensable and cannot be rescued by ADAR1p110. Nevertheless, editing sites uniquely targeted by ADAR1p150 remain elusive. Here, by transfection of ADAR1 isoforms into ADAR-less mouse cells we detect isoform-specific editing patterns. Using mutated ADAR variants, we test how intracellular localization and the presence of a Z-DNA binding domain-α affect editing preferences. These data show that ZBDα only minimally contributes to p150 editing-specificity while isoform-specific editing is primarily directed by the intracellular localization of ADAR1 isoforms. Our study is complemented by RIP-seq on human cells ectopically expressing tagged-ADAR1 isoforms. Both datasets reveal enrichment of intronic editing and binding by ADAR1p110 while ADAR1p150 preferentially binds and edits 3'UTRs.

Hydrophobicity-Driven Increases in Editing in Mitochondrial mRNAs during the Evolution of Kinetoplastids.

Kinetoplastids are a diverse group of flagellates which exhibit editing by insertion/deletion of Us in the mitochondrial mRNAs. Some mRNAs require editing to build most of their coding sequences, a process known as pan-editing. Evidence suggests that pan-editing is an ancestral feature in kinetoplastids. Here, we investigate how the transition from nonedited to pan-edited states occurred. The mitochondrial mRNAs and protein sequences from nine kinetoplastids and related groups (diplonemids, euglenids, and jakobids) were analyzed. RNA editing increased protein hydrophobicity to extreme values by introducing Us in the second codon position, despite the absence of editing preferences related to codon position. In addition, hydrophobicity was maintained by purifying selection in species that lost editing by retroposition of the fully edited mRNA. Only a few hydrophobic to hydrophilic amino acid changes were inferred for such species. In the protein secondary structure, these changes occurred spatially close to other hydrophilic residues. The analysis of coevolving sites showed that multiple changes are required together for hydrophobicity to be lost, which suggest the proteins are locked into extended hydrophobicity. Finally, an analysis of the NAD7 protein-protein interactions showed they can also influence hydrophobicity increase in the protein and where editing can occur in the mRNA. In conclusion, our results suggest that protein hydrophobicity has influenced editing site selection and how editing expanded in mRNAs. In effect, the hydrophobicity increase was entrenched by a neutral ratchet moved by a mutational pressure to introduce Us, thus helping to explain both RNA editing increase and, possibly, persistence.

Host A-to-I RNA editing signatures in intracellular bacterial and single-strand RNA viral infections.

Microbial infection is accompanied by remodeling of the host transcriptome. Involvement of A-to-I RNA editing has been reported during viral infection but remains to be elucidated during intracellular bacterial infections. Herein we analyzed A-to-I RNA editing during intracellular bacterial infections based on 18 RNA-Seq datasets of 210 mouse samples involving 7 tissue types and 8 intracellular bacterial pathogens (IBPs), and identified a consensus signature of RNA editing for IBP infections, mainly involving neutrophil-mediated innate immunity and lipid metabolism. Further comparison of host RNA editing patterns revealed remarkable similarities between pneumonia caused by IBPs and single-strand RNA (ssRNA) viruses, such as altered editing enzyme expression, editing site numbers, and levels. In addition, functional enrichment analysis of genes with RNA editing highlighted that the Rab GTPase family played a common and vital role in the host immune response to IBP and ssRNA viral infections, which was indicated by the consistent up-regulated RNA editing of Ras-related protein Rab27a. Nevertheless, dramatic differences between IBP and viral infections were also observed, and clearly distinguished the two types of intracellular infections. Our study showed transcriptome-wide host A-to-I RNA editing alteration during IBP and ssRNA viral infections. By identifying and comparing consensus signatures of host A-to-I RNA editing, our analysis implicates the importance of host A-to-I RNA editing during these infections and provides new insights into the diagnosis and treatment of infectious diseases.

The garden asparagus (Asparagus officinalis L.) mitochondrial genome revealed rich sequence variation throughout whole sequencing data

Garden asparagus (Asparagus officinalis L.) is a horticultural crop with high nutritional and medical value, considered an ideal plant for sex determination research among many dioecious plants, whose genomic information can support genetic analysis and breeding programs. In this research, the entire mitochondrial genome of A. officinalis was sequenced, annotated and assembled using a mixed Illumina and PacBio data. The garden asparagus circular mitochondrial genome measures 492,062 bp with a GC value of 45.9%. Thirty-six protein-coding genes, 17 tRNA and 6 rRNA genes were annotated, among which 8 protein-coding genes contained 16 introns. In addition, 254 SSRs with 10 complete tandem repeats and 293 non-tandem repeats were identified. It was found that the codons of edited sites located in the amino acids showed a leucine-formation trend, and RNA editing sites mainly caused the mutual transformation of amino acids with the same properties. Furthermore, 72 sequence fragments accounting for 20,240 bp, presentating 4.11% of the whole mitochondrial genome, were observed to migrate from chloroplast to mitochondrial genome of A. officinalis. The phylogenetic analysis showed that the closest genetic relationship between A. officinalis with onion (Allium cepa) inside the Liliaceae family. Our results demonstrated that high percentage of protein-coding genes had evolutionary conservative properties, with Ka/Ks values less than 1. Therefore, this study provides a high-quality garden asparagus mitochondrial genome, useful to promote better understanding of gene exchange between organelle genomes.

Selective Inhibition of ADAR1 Using 8-Azanebularine-Modified RNA Duplexes.

Adenosine deaminases acting on RNA (ADARs) are RNA editing enzymes that catalyze the hydrolytic deamination of adenosine (A) to inosine (I) in dsRNA. In humans, two catalytically active ADARs, ADAR1 and ADAR2, perform this A-to-I editing event. The growing field of nucleotide base editing has highlighted ADARs as promising therapeutic agents while multiple studies have also identified ADAR1's role in cancer progression. However, the potential for site-directed RNA editing as well as the rational design of inhibitors is being hindered by the lack of detailed molecular understanding of RNA recognition by ADAR1. Here, we designed short RNA duplexes containing the nucleoside analog, 8-azanebularine (8-azaN), to gain insight into molecular recognition by the human ADAR1 catalytic domain. From gel shift and in vitro deamination experiments, we validate ADAR1 catalytic domain's duplex secondary structure requirement and present a minimum duplex length for binding (14 bp, with 5 bp 5' and 8 bp 3' to editing site). These findings concur with predicted RNA-binding contacts from a previous structural model of the ADAR1 catalytic domain. Finally, we establish that neither 8-azaN as a free nucleoside nor a ssRNA bearing 8-azaN inhibits ADAR1 and demonstrate that the 8-azaN-modified RNA duplexes selectively inhibit ADAR1 and not the closely related ADAR2 enzyme.

Experimental evidence for the functional importance and adaptive advantage of A-to-I RNA editing in fungi.

Adenosine-to-inosine (A-to-I) editing is the most prevalent type of RNA editing in animals, and it occurs in fungi specifically during sexual reproduction. However, it is debatable whether A-to-I RNA editing is adaptive. Deciphering the functional importance of individual editing sites is essential for the mechanistic understanding of the adaptive advantages of RNA editing. Here, by performing gene deletion for 17 genes with conserved missense editing (CME) sites and engineering underedited (ue) and overedited (oe) mutants for 10 CME sites using site-specific mutagenesis at the native locus in Fusarium graminearum, we demonstrated that two CME sites in CME5 and CME11 genes are functionally important for sexual reproduction. Although the overedited mutant was normal in sexual reproduction, the underedited mutant of CME5 had severe defects in ascus and ascospore formation like the deletion mutant, suggesting that the CME site of CME5 is co-opted for sexual development. The preediting residue of Cme5 is evolutionarily conserved across diverse classes of Ascomycota, while the postediting one is rarely hardwired into the genome, implying that editing at this site leads to higher fitness than a genomic A-to-G mutation. More importantly, mutants expressing only the underedited or the overedited allele of CME11 are defective in ascosporogenesis, while those expressing both alleles displayed normal phenotypes, indicating that concurrently expressing edited and unedited versions of Cme11 is more advantageous than either. Our study provides convincing experimental evidence for the long-suspected adaptive advantages of RNA editing in fungi and likely in animals.

Rate accelerations in plastid and mitochondrial genomes of Cyperaceae occur in the same clades

Cyperaceae, the second largest family in the monocot order Poales, comprises >5500 species and includes the genus Eleocharis with ∼ 250 species. A previous study of complete plastomes of two Eleocharis species documented extensive structural heteroplasmy, gene order changes, high frequency of dispersed repeats along with gene losses and duplications. To better understand the phylogenetic distribution of gene and intron content as well as rates and patterns of sequence evolution within and between mitochondrial and plastid genomes of Eleocharis and Cyperaceae, an additional 29 Eleocharis organelle genomes were sequenced and analyzed. Eleocharis experienced extensive gene loss in both genomes while loss of introns was mitochondria-specific. Eleocharis has higher rates of synonymous (dS) and nonsynonymous (dN) substitutions in the plastid and mitochondrion than most sampled angiosperms, and the pattern was distinct from other eudicot lineages with accelerated rates. Several clades showed higher dS and dN in mitochondrial genes than in plastid genes. Furthermore, nucleotide substitution rates of mitochondrial genes were significantly accelerated on the branch leading to Cyperaceae compared to most angiosperms. Mitochondrial genes of Cyperaceae exhibited dramatic loss of RNA editing sites and a negative correlation between RNA editing and dS values was detected among angiosperms. Mutagenic retroprocessing and dysfunction of DNA replication, repair and recombination genes are the most likely cause of striking rate accelerations and loss of edit sites and introns in Eleocharis and Cyperaceae organelle genomes.

Genomic Landscape and Potential Regulation of RNA Editing in Drug Resistance.

Adenosine-to-inosine RNA editing critically affects the response of cancer therapies. However, comprehensive identification of drug resistance-related RNA editing events and systematic understanding of how RNA editing mediates anticancer drug resistance remain unclear. Here, 7157 differential editing sites (DESs) are identified from 98127 informative RNA editing sites in tumor tissues, many of which are validated in cancer cell lines. Diverse editing patterns of DESs are discovered in resistant samples, which could not be fully explained by adenosine deaminase acting on RNA enzymes. Some RNA-binding proteins are identified that potentially regulate these editing events. Notably, the DESs are significantly enriched in 3'-untranslated regions (3'-UTRs). The impact of DESs in 3'-UTR on the microRNA (miRNA) regulations is explored, and some triplets (DES, miRNA, and gene) that may contribute to drug resistance are identified. In addition, it is determined that the functions of genes enriched with DESs are associated with drug resistance, such as apoptosis, drug metabolism, and DNA synthesis involved in DNA repair. An online resource (http://www.jianglab.cn/REDR/) to support convenient retrieval of DESs is also built. The findings reveal the landscape and potential regulatory mechanism of RNA editing in drug resistance, providing new therapeutic targets for reversing drug resistance.

Genome-wide profiling of rice Double-stranded RNA-Binding Protein 1-associated RNAs by targeted RNA editing.

RNA-binding proteins (RBPs) play essential roles in regulating gene expression. However, the RNA ligands of RBPs are poorly understood in plants, not least due to the lack of efficient tools for genome-wide identification of RBP-bound RNAs. An RBP-fused adenosine deaminase acting on RNA (ADAR) can edit RBP-bound RNAs, which allows efficient identification of RNA ligands of RBPs in vivo. Here, we report the RNA editing activities of the ADAR deaminase domain (ADARdd) in plants. Protoplast experiments indicated that RBP-ADARdd fusions efficiently edited adenosines within 41 nucleotides (nt) of their binding sites. We then engineered ADARdd to profile the RNA ligands of rice (Oryza sativa) Double-stranded RNA-Binding Protein 1 (OsDRB1). Overexpressing the OsDRB1-ADARdd fusion protein in rice introduced thousands of A-to-G and T-to-C RNA‒DNA variants (RDVs). We developed a stringent bioinformatic approach to identify A-to-I RNA edits from RDVs, which removed 99.7% to 100% of background single-nucleotide variants in RNA-seq data. This pipeline identified a total of 1,798 high-confidence RNA editing (HiCE) sites, which marked 799 transcripts as OsDRB1-binding RNAs, from the leaf and root samples of OsDRB1-ADARdd-overexpressing plants. These HiCE sites were predominantly located in repetitive elements, 3'-UTRs, and introns. Small RNA sequencing also identified 191 A-to-I RNA edits in miRNAs and other sRNAs, confirming that OsDRB1 is involved in sRNA biogenesis or function. Our study presents a valuable tool for genome-wide profiling of RNA ligands of RBPs in plants and provides a global view of OsDRB1-binding RNAs.

Autocatalytic base editing for RNA-responsive translational control

Genetic circuits that control transgene expression in response to pre-defined transcriptional cues would enable the development of smart therapeutics. To this end, here we engineer programmable single-transcript RNA sensors in which adenosine deaminases acting on RNA (ADARs) autocatalytically convert target hybridization into a translational output. Dubbed DART VADAR (Detection and Amplification of RNA Triggers via ADAR), our system amplifies the signal from editing by endogenous ADAR through a positive feedback loop. Amplification is mediated by the expression of a hyperactive, minimal ADAR variant and its recruitment to the edit site via an orthogonal RNA targeting mechanism. This topology confers high dynamic range, low background, minimal off-target effects, and a small genetic footprint. We leverage DART VADAR to detect single nucleotide polymorphisms and modulate translation in response to endogenous transcript levels in mammalian cells.

Mitochondrial Genome Sequence of Salvia officinalis (Lamiales: Lamiaceae) Suggests Diverse Genome Structures in Cogeneric Species and Finds the Stop Gain of Genes through RNA Editing Events.

Our previous study was the first to confirm that the predominant conformation of mitochondrial genome (mitogenome) sequence of Salvia species contains two circular chromosomes. To further understand the organization, variation, and evolution of Salvia mitogenomes, we characterized the mitogenome of Salvia officinalis. The mitogenome of S. officinalis was sequenced using Illumina short reads and Nanopore long reads and assembled using a hybrid assembly strategy. We found that the predominant conformation of the S. officinalis mitogenome also had two circular chromosomes that were 268,341 bp (MC1) and 39,827 bp (MC2) in length. The S. officinalis mitogenome encoded an angiosperm-typical set of 24 core genes, 9 variable genes, 3 rRNA genes, and 16 tRNA genes. We found many rearrangements of the Salvia mitogenome through inter- and intra-specific comparisons. A phylogenetic analysis of the coding sequences (CDs) of 26 common protein-coding genes (PCGs) of 11 Lamiales species and 2 outgroup taxa strongly indicated that the S. officinalis was a sister taxon to S. miltiorrhiza, consistent with the results obtained using concatenated CDs of common plastid genes. The mapping of RNA-seq data to the CDs of PCGs led to the identification of 451 C-to-U RNA editing sites from 31 PCGs of the S. officinalis mitogenome. Using PCR amplification and Sanger sequencing methods, we successfully validated 113 of the 126 RNA editing sites from 11 PCGs. The results of this study suggest that the predominant conformation of the S. officinalis mitogenome are two circular chromosomes, and the stop gain of rpl5 was found through RNA editing events of the Salvia mitogenome.

Pan-cancer investigation of C-to-U editing reveals its important role in cancer development and new targets for cancer treatment.

RNA editing is prevalent in the transcriptome and is important for multiple cellular processes. C-to-U RNA editing sites (RES) are relatively rare and understudied in humans, compared to A-to-I editing. However, the functional impact of C-to-U editing in human cancers also remains elusive. Here, we conducted the first comprehensive survey of pan-cancer C-to-U RESs. Surprisingly, we found that the same subset of RESs were associated with multiple features, including patient survival, cancer stemness, tumor mutation burden (TMB), and tumor-infiltrated immune cell compositions (ICC), suggesting an RES-mediated close relationship between these features. For example, editing sites for GALM or IFI6 that led to higher expression were linked to lower survival and more cancer stemness. Also, TMB was found to be lower in prostate cancer cases with ICC-associated RESs in CAVIN1 or VWA8 or higher in prostate cancer cases with thymoma. With experimental support, we also found RESs in CST3, TPI1, or TNC that are linked to immune checkpoint blockade by anti-PD1. We also confirmed through experiments that two C-to-U RESs in CSNK2B or RPS14 had different effects on colon cancer cells. Patients with CSNK2B editing, which increased the expression of the oncogene CLDN18, had a lower response to drugs. On the other hand, drugs worked better on people who had RPS14 editing, which greatly increased ribosome production. In summary, our study demonstrated the important roles of C-to-U RESs across cancers and shed light on personalized cancer therapy.

A Database of Lung Cancer-Related Genes for the Identification of Subtype-Specific Prognostic Biomarkers

The molecular subtype is critical for accurate treatment and follow-up in patients with lung cancer; however, information regarding subtype-associated genes is dispersed among thousands of published studies. Systematic curation and cross-validation of the scientific literature would provide a solid foundation for comparative genetic studies of the major molecular subtypes of lung cancer. Here, we constructed a literature-based lung cancer gene database (LCGene). In the current release, we collected and curated 2507 unique human genes, including 2267 protein-coding and 240 non-coding genes from comprehensive manual examination of 10,960 PubMed article abstracts. Extensive annotations were added to aid identification of differentially expressed genes, potential gene editing sites, and non-coding gene regulation. For instance, we prepared 607 curated genes with CRISPR knockout information in 43 lung cancer cell lines. Further comparison of these implicated genes among different subtypes identified several subtype-specific genes with high mutational frequencies. Common tumor suppressors and oncogenes shared by lung adenocarcinoma and lung squamous cell carcinoma, for example, exhibited different mutational frequencies and prognostic features, suggesting the presence of subtype-specific biomarkers. Our retrospective analysis revealed 43 small cell lung cancer-specific genes. Moreover, 52 tumor suppressors and oncogenes shared by lung adenocarcinoma and squamous cell carcinoma confirmed the different molecular mechanisms of these two cancer subtypes. The subtype-based genetic differences, when combined, may provide insight into subtype-specific biomarkers for genetic testing.

Genome-wide identification of A-to-I RNA editing events provides the functional implications in PDAC.

RNA editing, a wide-acknowledged post-transcriptional mechanism, has been reported to be involved in the occurrence and development of cancer, especially the abnormal alteration of adenosine to inosine. However, fewer studies focus on pancreaticcancer. Therefore, we aimed to explore the possible linkages between altered RNA editing events and the development of PDAC. We characterized the global A-to-I RNA editing spectrum from RNA and matched whole-genome sequencing data of 41 primary PDAC and adjacent normal tissues. The following analyses were performed: different editing level and RNA expression analysis,pathway analysis, motif analysis, RNA secondary structure analysis, alternative splicing events analysis, and survival analysis.The RNA editing of single-cell RNA public sequencing data was also characterized. A large number of adaptive RNA editing events with significant differences in editing levels were identified, which are mainly regulated by ADAR1. Moreover, RNA editing in tumors has a higher editing level and more abundant editing sites in general. 140genes were screened out since they were identified with significantly different RNA editing events and were significantly different in expression level between tumor and matched normal samples. Further analysis showed a preference that in the tumor-specific group, they are mainly enriched in cancer-related signal pathways, while in the normal tissue-specific group, they are mainly enriched in pancreatic secretion. At the same time, we also found positively selected differentially edited sites in a series of cancer immune genes, including EGF, IGF1R, and PIK3CD. RNA editing might participate in pathogenisis of PDAC through regulating the alternative splicing and RNA secondary structure of important genesto further regulate gene expression and protein synthesis, including RAB27B and CERS4. Furthermore, single cell sequencing results showed that type2 ductal cells contributed the most to RNA editing events in tumors. RNA editing is an epigenetic mechanism involved in the occurrence and development of pancreatic cancer, which has the potential to diagnose of PDAC and is closely related to the prognosis.

De Novo Assembly and Comparative Analysis of the Complete Mitochondrial Genome of Chaenomeles speciosa (Sweet) Nakai Revealed the Existence of Two Structural Isomers.

As a valuable Chinese traditional medicinal species, Chaenomeles speciosa (Sweet) Nakai (C. speciosa) is a natural resource with significant economic and ornamental value. However, its genetic information is not well understood. In this study, the complete mitochondrial genome of C. speciosa was assembled and characterized to explore the repeat sequences, recombination events, rearrangements, and IGT, to predict RNA editing sites, and to clarify the phylogenetic and evolutionary relationship. The C. speciosa mitochondrial genome was found to have two circular chromosomes as its major conformation, with a total length of 436,464 bp and 45.2% GC content. The mitochondrial genome contained 54 genes, including 33 unique protein-coding genes, 18 tRNAs, and 3 rRNA genes. Seven pairs of repeat sequences involving recombination events were analyzed. Both the repeat pairs, R1 and R2, played significant roles in mediating the major and minor conformations. In total, 18 MTPTs were identified, 6 of which were complete tRNA genes. There were 454 RNA editing sites in the 33 protein-coding sequences predicted by the PREPACT3 program. A phylogenetic analysis based on 22 species of mitochondrial genomes was constructed and indicated highly conserved PCG sequences. Synteny analyses showed extensive genomic rearrangements in the mitochondrial genome of C. speciosa and closely related species. This work is the first to report the C. speciosa mitochondrial genome, which is of great significance for conducting additional genetic studies on this organism.

Extensive Recoding of the Neural Proteome in Cephalopods by RNA Editing.

The coleoid cephalopods have the largest brains, and display the most complex behaviors, of all invertebrates. The molecular and cellular mechanisms that underlie these remarkable advancements remain largely unexplored. Early molecular cloning studies of squid ion channel transcripts uncovered an unusually large number of A→I RNA editing sites that recoded codons. Further cloning of other neural transcripts showed a similar pattern. The advent of deep-sequencing technologies and the associated bioinformatics allowed the mapping of RNA editing events across the entire neural transcriptomes of various cephalopods. The results were remarkable: They contained orders of magnitude more recoding editing sites than any other taxon. Although RNA editing sites are abundant in most multicellular metazoans, they rarely recode. In cephalopods, the majority of neural transcripts are recoded. Recent studies have focused on whether these events are adaptive, as well as other noncanonical aspects of cephalopod RNA editing.

Prognostic RNA-editing signature predicts immune functions and therapy responses in gliomas.

Background: RNA-editing refers to post-transcriptional transcript alterations that lead to the formation of protein isoforms and the progression of various tumors. However, little is known about its roles in gliomas. Aim: The aim of this study is to identify prognosis-related RNA-editing sites (PREs) in glioma, and to explore their specific effects on glioma and potential mechanisms of action. Methods: Glioma genomic and clinical data were obtained from TCGA database and SYNAPSE platform. The PREs was identified with regression analyses and the corresponding prognostic model was evaluated with survival analysis and receiver operating characteristic curve. Functional enrichment of differentially expressed genes between risk groups was performed to explore action mechanisms. The CIBERSORT, ssGSEA, gene set variation analysis, and ESTIMATE algorithms were employed to assess the association between PREs risk score and variations of tumor microenvironment, immune cell infiltration, immune checkpoints, and immune responses. The maftools and pRRophetic packages were used to evaluate tumor mutation burden and predict drug sensitivity. Results: A total of thirty-five RNA-editing sites were identified as prognosis-related in glioma. Functional enrichment implied variation of immune-related pathways between groups. Notably, glioma samples with higher PREs risk score exhibited higher immune score, lower tumor purity, increased infiltration of macrophage and regulatory T cells, suppressed NK cell activation, elevated immune function score, upregulated immune checkpoint gene expression, and higher tumor mutation burden, all of which implied worse response to immune therapy. Finally, high-risk glioma samples are more sensitive to Z-LLNle-CHO and temozolomide, while the low-risk ones respond better to Lisitinib. Conclusion: We identified a PREs signature of thirty-five RNA editing sites and calculated their corresponding risk coefficients. Higher total signature risk score indicates worse prognosis and worse immune response and lower sensitivity to immune therapy. The novel PREs signature could help risk stratification, immunotherapy response prediction, individualized treatment strategy-making for glioma patients, and development of novel therapeutic approaches.