Cellular and genetic drivers of RNA editing variation in the human brain

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk postmortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA editing and alternative polyadenylation (APA) within a cell-type-specific population of human neural progenitors and neurons. More RNA editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting that genetically mediated post-transcriptional regulation during brain development leads to differences in brain function.

Identifying cancer cells from calling single-nucleotide variants in scRNA-seq data

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) data are widely used to study cancer cell states and their heterogeneity. However, the tumour microenvironment is usually a mixture of healthy and cancerous cells and it can be difficult to fully separate these two populations based on transcriptomics alone. If available, somatic single-nucleotide variants (SNVs) observed in the scRNA-seq data could be used to identify the cancer population and match that information with the single cells' expression profile. However, calling somatic SNVs in scRNA-seq data is a challenging task, as most variants seen in the short-read data are not somatic, but can instead be germline variants, RNA edits or transcription, sequencing, or processing errors. In addition, only variants present in actively transcribed regions for each individual cell will be seen in the data.

RESULTS

To address these challenges, we develop CCLONE (Cancer Cell Labelling On Noisy Expression), an interpretable tool adapted to handle the uncertainty and sparsity of SNVs called from scRNA-seq data. CCLONE jointly identifies cancer clonal populations, and their associated variants. We apply CCLONE on two acute myeloid leukaemia datasets and one lung adenocarcinoma dataset and show that CCLONE captures both genetic clones and somatic events for multiple patients. These results show how CCLONE can be used to gather insight into the course of the disease and the origin of cancer cells in scRNA-seq data.

AVAILABILITY AND IMPLEMENTATION

Source code is available at github.com/HaghverdiLab/CCLONE.

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.

ADAR1 prevents ZBP1-dependent PANoptosis via A-to-I RNA editing in developmental sevoflurane neurotoxicity

It is well established that sevoflurane exposure leads to widespread neuronal cell death in the developing brain. Adenosine deaminase acting on RNA-1 (ADAR1) dependent adenosine-to-inosine (A-to-I) RNA editing is dynamically regulated throughout brain development. The current investigation is designed to interrogate the contributed role of ADAR1 in developmental sevoflurane neurotoxicity. Herein, we provide evidence to show that developmental sevoflurane priming triggers neuronal pyroptosis, apoptosis and necroptosis (PANoptosis), and elicits the release of inflammatory factors including IL-1β, IL-18, TNF-α and IFN-γ. Additionally, ADAR1-P150, but not ADAR1-P110, depresses cellular PANoptosis and inflammatory response by competing with Z-DNA/RNA binding protein 1 (ZBP1) for binding to Z-RNA in the presence of sevoflurane. Further investigation demonstrates that ADAR1-dependent A-to-I RNA editing mitigates developmental sevoflurane-induced neuronal PANoptosis. To restore RNA editing, we utilize adeno-associated virus (AAV) to deliver engineered circular ADAR-recruiting guide RNAs (cadRNAs) into cells, which is capable of recruiting endogenous adenosine deaminases to promote cellular A-to-I RNA editing. As anticipated, AAV-cadRNAs diminishes sevoflurane-induced cellular Z-RNA production and PANoptosis, which could be abolished by ADAR1-P150 shRNA transfection. Moreover, AAV-cadRNAs delivery ameliorates developmental sevoflurane-induced spatial and emotional cognitive deficits without influence on locomotor activity. Taken together, these results illustrate that ADAR1-P150 exhibits a prominent role in preventing ZBP1-dependent PANoptosis through A-to-I RNA editing in developmental sevoflurane neurotoxicity. Application of engineered cadRNAs to rectify the compromised ADAR1-dependent A-to-I RNA editing provides an inspiring direction for possible clinical preventions and therapeutics.

Divergent landscapes of A-to-I editing in postmortem and living human brain

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries offer more nuanced and accurate insights into the regulatory mechanisms of RNA editing in the human brain.

Inverted Alu repeats: friends or foes in the human transcriptome

Alu elements are highly abundant primate-specific short interspersed nuclear elements that account for ~10% of the human genome. Due to their preferential location in gene-rich regions, especially in introns and 3' UTRs, Alu elements can exert regulatory effects on the expression of both host and neighboring genes. When two Alu elements with inverse orientations are positioned in close proximity, their transcription results in the generation of distinct double-stranded RNAs (dsRNAs), known as inverted Alu repeats (IRAlus). IRAlus are key immunogenic self-dsRNAs and post-transcriptional cis-regulatory elements that play a role in circular RNA biogenesis, as well as RNA transport and stability. Recently, IRAlus dsRNAs have emerged as regulators of transcription and activators of Z-DNA-binding proteins. The formation and activity of IRAlus can be modulated through RNA editing and interactions with RNA-binding proteins, and misregulation of IRAlus has been implicated in several immune-associated disorders. In this review, we summarize the emerging functions of IRAlus dsRNAs, the regulatory mechanisms governing IRAlus activity, and their relevance in the pathogenesis of human diseases.

Imbalance in Unc80 RNA Editing Disrupts Dynamic Neuronal Activity and Olfactory Perception

A-to-I RNA editing, catalyzed by the ADAR protein family, significantly contributes to the diversity and adaptability of mammalian RNA signatures, aligning with developmental and physiological needs. Yet, the functions of many editing sites are still to be defined. The Unc80 gene stands out in this context due to its brain-specific expression and the evolutionary conservation of its codon-altering editing event. The precise biological functions of Unc80 and its editing, however, are still largely undefined. In this study, we first demonstrated that Unc80 editing occurs in an ADAR2-dependent manner and is exclusive to the brain. By employing the CRISPR/Cas9 system to generate Unc80 knock-in mouse models that replicate the natural editing variations, our findings revealed that mice with the "gain-of-editing" variant (Unc80G/G) exhibit heightened basal neuronal activity in critical olfactory regions, compared to the "loss-of-editing" (Unc80S/S) counterparts. Moreover, an increase in glutamate levels was observed in the olfactory bulbs of Unc80G/G mice, indicating altered neurotransmitter dynamics. Behavioral analysis of odor detection revealed distinctive responses to novel odors-both Unc80 deficient (Unc80+/-) and Unc80S/S mice demonstrated prolonged exploration times and heightened dishabituation responses. Further elucidating the olfactory connection of Unc80 editing, transcriptomic analysis of the olfactory bulb identified significant alterations in gene expression that corroborate the behavioral and physiological findings. Collectively, our research advances the understanding of Unc80's neurophysiological functions and the impact of its editing on the olfactory sensory system, shedding light on the intricate molecular underpinnings of olfactory perception and neuronal activity.

Divergent landscapes of A-to-I editing in postmortem and living human brain

Adenosine-to-inosine (A-to-I) editing is a prevalent post-transcriptional RNA modification within the brain. Yet, most research has relied on postmortem samples, assuming it is an accurate representation of RNA biology in the living brain. We challenge this assumption by comparing A-to-I editing between postmortem and living prefrontal cortical tissues. Major differences were found, with over 70,000 A-to-I sites showing higher editing levels in postmortem tissues. Increased A-to-I editing in postmortem tissues is linked to higher ADAR1 and ADARB1 expression, is more pronounced in non-neuronal cells, and indicative of postmortem activation of inflammation and hypoxia. Higher A-to-I editing in living tissues marks sites that are evolutionarily preserved, synaptic, developmentally timed, and disrupted in neurological conditions. Common genetic variants were also found to differentially affect A-to-I editing levels in living versus postmortem tissues. Collectively, these discoveries illuminate the nuanced functions and intricate regulatory mechanisms of RNA editing within the human brain.

Temporal landscape and translational regulation of A-to-I RNA editing in mouse retina development

BACKGROUND

The significance of A-to-I RNA editing in nervous system development is widely recognized; however, its influence on retina development remains to be thoroughly understood.

RESULTS

In this study, we performed RNA sequencing and ribosome profiling experiments on developing mouse retinas to characterize the temporal landscape of A-to-I editing. Our findings revealed temporal changes in A-to-I editing, with distinct editing patterns observed across different developmental stages. Further analysis showed the interplay between A-to-I editing and alternative splicing, with A-to-I editing influencing splicing efficiency and the quantity of splicing events. A-to-I editing held the potential to enhance translation diversity, but this came at the expense of reduced translational efficiency. When coupled with splicing, it could produce a coordinated effect on gene translation.

CONCLUSIONS

Overall, this study presents a temporally resolved atlas of A-to-I editing, connecting its changes with the impact on alternative splicing and gene translation in retina development.

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.

In search of critical dsRNA targets of ADAR1

Recent studies have underscored the pivotal role of adenosine-to-inosine RNA editing, catalyzed by ADAR1, in suppressing innate immune interferon responses triggered by cellular double-stranded RNA (dsRNA). However, the specific ADAR1 editing targets crucial for this regulatory function remain elusive. We review analyses of transcriptome-wide ADAR1 editing patterns and their evolutionary dynamics, which offer valuable insights into this unresolved query. The growing appreciation of the significance of immunogenic dsRNAs and their editing in inflammatory and autoimmune diseases and cancer calls for a more comprehensive understanding of dsRNA immunogenicity, which may promote our understanding of these diseases and open doors to therapeutic avenues.

Quantification of AMPA receptor subunits and RNA editing-related proteins in the J20 mouse model of Alzheimer's disease by capillary western blotting

Introduction

Accurate modelling of molecular changes in Alzheimer's disease (AD) dementia is crucial for understanding the mechanisms driving neuronal pathology and for developing treatments. Synaptic dysfunction has long been implicated as a mechanism underpinning memory dysfunction in AD and may result in part from changes in adenosine deaminase acting on RNA (ADAR) mediated RNA editing of the GluA2 subunit of AMPA receptors and changes in AMPA receptor function at the post synaptic cleft. However, few studies have investigated changes in proteins which influence RNA editing and notably, AD studies that focus on studying changes in protein expression, rather than changes in mRNA, often use traditional western blotting.

Methods

Here, we demonstrate the value of automated capillary western blotting to investigate the protein expression of AMPA receptor subunits (GluA1-4), the ADAR RNA editing proteins (ADAR1-3), and proteins known to regulate RNA editing (PIN1, WWP2, FXR1P, and CREB1), in the J20 AD mouse model. We describe extensive optimisation and validation of the automated capillary western blotting method, demonstrating the use of total protein to normalise protein load, in addition to characterising the optimal protein/antibody concentrations to ensure accurate protein quantification. Following this, we assessed changes in proteins of interest in the hippocampus of 44-week-old J20 AD mice.

Results

We observed an increase in the expression of ADAR1 p110 and GluA3 and a decrease in ADAR2 in the hippocampus of 44-week-old J20 mice. These changes signify a shift in the balance of proteins that play a critical role at the synapse. Regression analysis revealed unique J20-specific correlations between changes in AMPA receptor subunits, ADAR enzymes, and proteins that regulate ADAR stability in J20 mice, highlighting potential mechanisms mediating RNA-editing changes found in AD.

Discussion

Our findings in J20 mice generally reflect changes seen in the human AD brain. This study underlines the importance of novel techniques, like automated capillary western blotting, to assess protein expression in AD. It also provides further evidence to support the hypothesis that a dysregulation in RNA editing-related proteins may play a role in the initiation and/or progression of AD.

Detection of Developmental Asexual Stage-Specific RNA Editing Events in Plasmodium falciparum 3D7 Malaria Parasite

Transcriptional variation has been studied but post-transcriptional modification due to RNA editing has not been investigated in Plasmodium. We investigated developmental stage-specific RNA editing in selected genes in Plasmodium falciparum 3D7. We detected extensive amination- and deamination-type RNA editing at 8, 16, 24, 32, 40, and 46 h in tightly synchronized Plasmodium. Most of the editing events were observed in 8 and 16 h ring-stage parasites. Extensive A-to-G deamination-type editing was detected more during the 16 h ring stage (25%) than the 8 h ring stage (20%). Extensive U-to-C amination-type editing was detected more during the 16 h ring stage (31%) than the 8 h ring stage (22%). In 28S, rRNA editing converted the loop structure to the stem structure. The hemoglobin binding activity of PF3D7_0216900 was also altered due to RNA editing. Among the expressed 28S rRNA genes, PF3D7_0532000 and PF3D7_0726000 expression was higher. Increased amounts of the transcripts of these two genes were found, particularly PF3D7_0726000 in the ring stage and PF3D7_0532000 in the trophozoite and schizont stages. Adenosine deaminase (ADA) expression did not correlate with the editing level. This first experimental report of RNA editing will help to identify the editing machinery that might be useful for antimalarial drug discovery and malaria control.

Changes in ADAR RNA editing patterns in CMV and ZIKV congenital infections

BACKGROUND

RNA editing is a process that increases transcriptome diversity, often through Adenosine Deaminases Acting on RNA (ADARs) that catalyze the deamination of adenosine to inosine. ADAR editing plays an important role in regulating brain function and immune activation, and is dynamically regulated during brain development. Additionally, the ADAR1 p150 isoform is induced by interferons in viral infection and plays a role in antiviral immune response. However, the question of how virus-induced ADAR expression affects host transcriptome editing remains largely unanswered. This question is particularly relevant in the context of congenital infections, given the dynamic regulation of ADAR editing during brain development, the importance of this editing for brain function, and subsequent neurological symptoms of such infections, including microcephaly, sensory issues, and other neurodevelopmental abnormalities. Here, we begin to address this question, examining ADAR expression in publicly available datasets of congenital infections of human cytomegalovirus (HCMV) microarray expression data, as well as mouse cytomegalovirus (MCMV) and mouse/ human induced pluripotent neuroprogenitor stem cell (hiNPC) Zika virus (ZIKV) RNA-seq data.

RESULTS

We found that in all three datasets, ADAR1 was overexpressed in infected samples compared to uninfected samples. In the RNA-seq datasets, editing rates were also analyzed. In all mouse infections cases, the number of editing sites was significantly increased in infected samples, albeit this was not the case for hiNPC ZIKV samples. Mouse ZIKV samples showed altered editing of well-established protein-recoding sites such as Gria3, Grik5, and Nova1, as well as editing sites that may impact miRNA binding.

CONCLUSIONS

Our findings provide evidence for changes in ADAR expression and subsequent dysregulation of ADAR editing of host transcriptomes in congenital infections. These changes in editing patterns of key neural genes have potential significance in the development of neurological symptoms, thus contributing to neurodevelopmental abnormalities. Further experiments should be performed to explore the full range of editing changes that occur in different congenital infections, and to confirm the specific functional consequences of these editing changes.

Identification of Bona Fide RNA Editing Sites: History, Challenges, and Opportunities

Adenosine-to-inosine (A-to-I) RNA editing, catalyzed by the adenosine deaminase acting on the RNA (ADAR) family of enzymes of which there are three members (ADAR1, ADAR2, and ADAR3), is a major gene regulatory mechanism that diversifies the transcriptome. It is widespread in many metazoans, including humans. As inosine is interpreted by cellular machineries mainly as guanosine, A-to-I editing effectively gives A-to-G nucleotide changes. Depending on its location, an editing event can generate new protein isoforms or influence other RNA processing pathways. Researchers have found that ADAR-mediated editing performs diverse functions. For example, it enables living organisms such as cephalopods to adapt rapidly to fluctuating environmental conditions such as water temperature. In development, the loss of ADAR1 is embryonically lethal partly because endogenous double-stranded RNAs (dsRNAs) are no longer marked by inosines, which signal "self", and thus cause the melanoma differentiation-associated protein 5 (MDA5) sensor to trigger a deleterious interferon response. Hence, ADAR1 plays a key role in preventing aberrant activation of the innate immune system. Furthermore, ADAR enzymes have been implicated in myriad human diseases. Intriguingly, some cancer cells are known to exploit ADAR1 activity to dodge immune responses. However, the exact identities of immunogenic RNAs in different biological contexts have remained elusive. Consequently, there is tremendous interest in identifying inosine-containing RNAs in the cell.The identification of A-to-I RNA editing sites is dependent on the sequencing of nucleic acids. Technological and algorithmic advancements over the past decades have revolutionized the way editing events are detected. At the beginning, the discovery of editing sites relies on Sanger sequencing, a first-generation technology. Both RNA, which is reverse transcribed into complementary DNA (cDNA), and genomic DNA (gDNA) from the same source are analyzed. After sequence alignment, one would require an adenosine to be present in the genome but a guanosine to be detected in the RNA sample for a position to be declared as an editing site. However, an issue with Sanger sequencing is its low throughput. Subsequently, Illumina sequencing, a second-generation technology, was invented. By permitting the simultaneous interrogation of millions of molecules, it enables many editing sites to be identified rapidly. However, a key challenge is that the Illumina platform produces short sequencing reads that can be difficult to map accurately. To tackle the challenge, we and others developed computational workflows with a series of filters to discard sites that are likely to be false positives. When Illumina sequencing data sets are properly analyzed, A-to-G variants should emerge as the most dominant mismatch type. Moreover, the quantitative nature of the data allows us to build a comprehensive atlas of editing-level measurements across different biological contexts, providing deep insights into the spatiotemporal dynamics of RNA editing. However, difficulties remain in identifying true A-to-I editing sites in short protein-coding exons or in organisms and diseases where DNA mutations and genomic polymorphisms are prevalent and mostly unknown. Nanopore sequencing, a third-generation technology, promises to address the difficulties, as it allows native RNAs to be sequenced without conversion to cDNA, preserving base modifications that can be directly detected through machine learning. We recently demonstrated that nanopore sequencing could be used to identify A-to-I editing sites in native RNA directly. Although further work is needed to enhance the detection accuracy in single molecules from fewer cells, the nanopore technology holds the potential to revolutionize epitranscriptomic studies.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

Genetics of cell-type-specific post-transcriptional gene regulation during human neurogenesis

The function of some genetic variants associated with brain-relevant traits has been explained through colocalization with expression quantitative trait loci (eQTL) conducted in bulk post-mortem adult brain tissue. However, many brain-trait associated loci have unknown cellular or molecular function. These genetic variants may exert context-specific function on different molecular phenotypes including post-transcriptional changes. Here, we identified genetic regulation of RNA-editing and alternative polyadenylation (APA), within a cell-type-specific population of human neural progenitors and neurons. More RNA-editing and isoforms utilizing longer polyadenylation sequences were observed in neurons, likely due to higher expression of genes encoding the proteins mediating these post-transcriptional events. We also detected hundreds of cell-type-specific editing quantitative trait loci (edQTLs) and alternative polyadenylation QTLs (apaQTLs). We found colocalizations of a neuron edQTL in CCDC88A with educational attainment and a progenitor apaQTL in EP300 with schizophrenia, suggesting genetically mediated post-transcriptional regulation during brain development lead to differences in brain function.

Transcriptome-wide profiling and quantification of N6-methyladenosine by enzyme-assisted adenosine deamination

N6-methyladenosine (m6A), the most abundant internal messenger RNA modification in higher eukaryotes, serves myriad roles in regulating cellular processes. Functional dissection of m6A is, however, hampered in part by the lack of high-resolution and quantitative detection methods. Here we present evolved TadA-assisted N6-methyladenosine sequencing (eTAM-seq), an enzyme-assisted sequencing technology that detects and quantifies m6A by global adenosine deamination. With eTAM-seq, we analyze the transcriptome-wide distribution of m6A in HeLa and mouse embryonic stem cells. The enzymatic deamination route employed by eTAM-seq preserves RNA integrity, facilitating m6A detection from limited input samples. In addition to transcriptome-wide m6A profiling, we demonstrate site-specific, deep-sequencing-free m6A quantification with as few as ten cells, an input demand orders of magnitude lower than existing quantitative profiling methods. We envision that eTAM-seq will enable researchers to not only survey the m6A landscape at unprecedented resolution, but also detect m6A at user-specified loci with a simple workflow.

A novel computational method enables RNA editome profiling during human hematopoiesis from scRNA-seq data

RNA editing is a post-transcriptional modification with a cell-specific manner and important biological implications. Although single-cell RNA-seq (scRNA-seq) is an effective method for studying cellular heterogeneity, it is difficult to detect and study RNA editing events from scRNA-seq data because of the low sequencing coverage. To overcome this, we develop a computational method to systematically identify RNA editing sites of cell types from scRNA-seq data. To demonstrate its effectiveness, we apply it to scRNA-seq data of human hematopoietic stem/progenitor cells (HSPCs) with an annotated lineage differentiation relationship according to previous research and study the impacts of RNA editing on hematopoiesis. The dynamic editing patterns reveal the relevance of RNA editing on different HSPCs. For example, four microRNA (miRNA) target sites on 3' UTR of EIF2AK2 are edited across all HSPC populations, which may abolish the miRNA-mediated inhibition of EIF2AK2. Elevated EIF2AK2 may thus activate the integrated stress response (ISR) pathway to initiate global translational attenuation as a protective mechanism to maintain cellular homeostasis during HSPCs' differentiation. Besides, our findings also indicate that RNA editing plays an essential role in the coordination of lineage commitment and self-renewal of hematopoietic stem cells (HSCs). Taken together, we demonstrate the capacity of scRNA-seq data to exploit RNA editing events of cell types, and find that RNA editing may exert multiple modules of regulation in hematopoietic processes.

Spatiotemporal and genetic regulation of A-to-I editing throughout human brain development

Posttranscriptional RNA modifications by adenosine-to-inosine (A-to-I) editing are abundant in the brain, yet elucidating functional sites remains challenging. To bridge this gap, we investigate spatiotemporal and genetically regulated A-to-I editing sites across prenatal and postnatal stages of human brain development. More than 10,000 spatiotemporally regulated A-to-I sites were identified that occur predominately in 3' UTRs and introns, as well as 37 sites that recode amino acids in protein coding regions with precise changes in editing levels across development. Hyper-edited transcripts are also enriched in the aging brain and stabilize RNA secondary structures. These features are conserved in murine and non-human primate models of neurodevelopment. Finally, thousands of cis-editing quantitative trait loci (edQTLs) were identified with unique regulatory effects during prenatal and postnatal development. Collectively, this work offers a resolved atlas linking spatiotemporal variation in editing levels to genetic regulatory effects throughout distinct stages of brain maturation.