Search

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.

Supplementary Figures S1-11 and Tables S1-2 are included in this file

RNA editing modifies single nucleotides of RNAs, regulating primary protein structure and protein abundance. In recent years, the diversity of proteins and complexity of gene regulation associated with RNA editing dysregulation has been increasingly appreciated in oncology. Large-scale shifts in editing have been observed in bulk tumors across various cancer types. However, RNA editing in single cells and individual cell types within tumors has not been explored. By profiling editing in single cells from lung adenocarcinoma biopsies, we found that the increased editing trend of bulk lung tumors was unique to cancer cells. Elevated editing levels were observed in cancer cells resistant to targeted therapy, and editing sites associated with drug response were enriched. Consistent with the regulation of antiviral pathways by RNA editing, higher editing levels in cancer cells were associated with reduced antitumor innate immune response, especially levels of natural killer cell infiltration. In addition, the level of RNA editing in cancer cells was positively associated with somatic point mutation burden. This observation motivated the definition of a new metric, RNA editing load, reflecting the amount of RNA mutations created by RNA editing. Importantly, in lung cancer, RNA editing load was a stronger predictor of patient survival than DNA mutations. This study provides the first single cell dissection of editing in cancer and highlights the significance of RNA editing load in cancer prognosis.Significance:RNA editing analysis in single lung adenocarcinoma cells uncovers RNA mutations that correlate with tumor mutation burden and cancer innate immunity and reveals the amount of RNA mutations that strongly predicts patient survival.See related commentary by Luo and Liang, p. 351

RNA splicing plays a critical role in post-transcriptional gene regulation. Exponential expansion of intron length poses a challenge for accurate splicing. Little is known about how cells prevent inadvertent and often deleterious expression of intronic elements due to cryptic splicing. In this study, we identify hnRNPM as an essential RNA binding protein that suppresses cryptic splicing through binding to deep introns, preserving transcriptome integrity. Long interspersed nuclear elements (LINEs) harbor large amounts of pseudo splice sites in introns. hnRNPM preferentially binds at intronic LINEs and represses LINE-containing pseudo splice site usage for cryptic splicing. Remarkably, a subgroup of the cryptic exons can form long dsRNAs through base-pairing of inverted Alu transposable elements scattered in between LINEs and trigger interferon immune response, a well-known antiviral defense mechanism. Notably, these interferon-associated pathways are found to be upregulated in hnRNPM-deficient tumors, which also exhibit elevated immune cell infiltration. These findings unveil hnRNPM as a guardian of transcriptome integrity. Targeting hnRNPM in tumors may be used to trigger an inflammatory immune response thereby boosting cancer surveillance.

RNA editing modifies single nucleotides of RNAs, regulating primary protein structure and protein abundance. In recent years, the diversity of proteins and complexity of gene regulation associated with RNA editing dysregulation has been increasingly appreciated in oncology. Large-scale shifts in editing have been observed in bulk tumors across various cancer types. However, RNA editing in single cells and individual cell types within tumors has not been explored. By profiling editing in single cells from lung adenocarcinoma biopsies, we found that the increased editing trend of bulk lung tumors was unique to cancer cells. Elevated editing levels were observed in cancer cells resistant to targeted therapy, and editing sites associated with drug response were enriched. Consistent with the regulation of antiviral pathways by RNA editing, higher editing levels in cancer cells were associated with reduced antitumor innate immune response, especially levels of natural killer cell infiltration. In addition, the level of RNA editing in cancer cells was positively associated with somatic point mutation burden. This observation motivated the definition of a new metric, RNA editing load, reflecting the amount of RNA mutations created by RNA editing. Importantly, in lung cancer, RNA editing load was a stronger predictor of patient survival than DNA mutations. This study provides the first single cell dissection of editing in cancer and highlights the significance of RNA editing load in cancer prognosis. Significance: RNA editing analysis in single lung adenocarcinoma cells uncovers RNA mutations that correlate with tumor mutation burden and cancer innate immunity and reveals the amount of RNA mutations that strongly predicts patient survival. See related commentary by Luo and Liang, p. 351

To identify genetic pathways for addiction, we analyzed intravenous self-administration of cocaine or saline in a panel of 84 inbred and recombinant inbred mouse strains over 10 days. We integrated the behavior data with RNA-Seq data from the medial frontal cortex and nucleus accumbens from 41 strains. The self-administration of cocaine and saline showed distinct genetic bases. We maximized power to map loci for cocaine intake by using a linear mixed model to account for this longitudinal phenotype while correcting for population structure. A total of 15 unique significant loci were identified in the genome-wide association study (GWAS). A transcriptome-wide association study (TWAS) highlighted theTrpv2ion channel as a key locus for cocaine self-administration from the GWAS. In addition, 17 genes supplementary to the GWAS were identified includingArhgef26, Slc18b1andSlco5a1. We found numerous instances where alternate splice site selection or RNA editing altered transcript abundance. Our work emphasizes the importance ofTrpv2, a known cannabinoid receptor, for the response to cocaine as well as identifying further relevant loci.

The zebra finch (Taeniopygia guttata), a representative oscine songbird species, has been widely studied to investigate behavioral neuroscience, most notably the neurobiological basis of vocal learning, a rare trait shared in only a few animal groups including humans. In 2019, an updated zebra finch genome annotation (bTaeGut1_v1.p) was released from the Ensembl database and is substantially more comprehensive than the first version published in 2010. In this study, we utilized the publicly available RNA-seq data generated from Illumina-based short-reads and PacBio single-molecule real-time (SMRT) long-reads to assess the bird transcriptome. To analyze the high-throughput RNA-seq data, we adopted a hybrid bioinformatic approach combining short and long-read pipelines. From our analysis, we added 220 novel genes and 8,134 transcript variants to the Ensembl annotation, and predicted a new proteome based on the refined annotation. We further validated 18 different novel proteins by using mass-spectrometry data generated from zebra finch caudal telencephalon tissue. Our results provide additional resources for future studies of zebra finches utilizing this improved bird genome annotation and proteome.

Human stem-cell-derived models provide the promise of accelerating our understanding of brain disorders, but not knowing whether they possess the ability to mature beyond mid- to late-fetal stages potentially limits their utility. We leveraged a directed differentiation protocol to comprehensively assess maturation in vitro. Based on genome-wide analysis of the epigenetic clock and transcriptomics, as well as RNA editing, we observe that three-dimensional human cortical organoids reach postnatal stages between 250 and 300 days, a timeline paralleling in vivo development. We demonstrate the presence of several known developmental milestones, including switches in the histone deacetylase complex and NMDA receptor subunits, which we confirm at the protein and physiological levels. These results suggest that important components of an intrinsic in vivo developmental program persist in vitro. We further map neurodevelopmental and neurodegenerative disease risk genes onto in vitro gene expression trajectories to provide a resource and webtool (Gene Expression in Cortical Organoids, GECO) to guide disease modeling.

Alternative splicing is an RNA processing mechanism that affects most genes in human, contributing to disease mechanisms and phenotypic diversity. The regulation of splicing involves an intricate network of cis-regulatory elements and trans-acting factors. Due to their high sequence specificity, cis-regulation of splicing can be altered by genetic variants, significantly affecting splicing outcomes. Recently, multiple methods have been applied to understanding the regulatory effects of genetic variants on splicing. However, it is still challenging to go beyond apparent association to pinpoint functional variants. To fill in this gap, we utilized large-scale data sets of the Genotype-Tissue Expression (GTEx) project to study genetically modulated alternative splicing (GMAS) via identification of allele-specific splicing events. We demonstrate that GMAS events are shared across tissues and individuals more often than expected by chance, consistent with their genetically driven nature. Moreover, although the allelic bias of GMAS exons varies across samples, the degree of variation is similar across tissues versus individuals. Thus, genetic background drives the GMAS pattern to a similar degree as tissue-specific splicing mechanisms. Leveraging the genetically driven nature of GMAS, we developed a new method to predict functional splicing-altering variants, built upon a genotype-phenotype concordance model across samples. Complemented by experimental validations, this method predicted >1000 functional variants, many of which may alter RNA-protein interactions. Lastly, 72% of GMAS-associated SNPs were in linkage disequilibrium with GWAS-reported SNPs, and such association was enriched in tissues of relevance for specific traits/diseases. Our study enables a comprehensive view of genetically driven splicing variations in human tissues.

Single-cell RNA sequencing (scRNA-seq) data contain rich information at the gene, transcript, and nucleotide levels. Most analyses of scRNA-seq have focused on gene expression profiles, and it remains challenging to extract nucleotide variants and isoform-specific information. Here, we present scAllele, an integrative approach that detects single-nucleotide variants, insertions, deletions, and their allelic linkage with splicing patterns in scRNA-seq. We demonstrate that scAllele achieves better performance in identifying nucleotide variants than other commonly used tools. In addition, the read-specific variant calls by scAllele enables allele-specific splicing analysis, a unique feature not afforded by other methods. Applied to a lung cancer scRNA-seq dataset, scAllele identified variants with strong allelic linkage to alternative splicing, some of which are cancer specific and enriched in cancer-relevant pathways. scAllele represents a versatile tool to uncover multilayer information and previously unidentified biological insights from scRNA-seq data.

Using third-generation sequencers, long-read RNA-seq is increasingly applied in transcriptomic studies given its major advantage in characterizing full-length transcripts. A number of methods have been developed to analyze this new type of data for transcript isoforms and their abundance. Another application, which is significantly under-explored, is to identify and analyze single nucleotide variants (SNVs) in the RNA. Identification of SNVs, such as genetic mutations or RNA editing sites, is fundamental to many biomedical questions. In long-read RNA-seq, SNV analysis presents significant challenges, due to the well-known relatively high error rates of the third-generation sequencers. Here, we present the first study to detect and analyze RNA editing sites in long-read RNA-seq. Our new method, L-GIREMI, effectively handles sequencing errors and biases in the reads, and uses a model-based approach to score RNA editing sites. Applied to PacBio long-read RNA-seq data, L-GIREMI affords a high accuracy in RNA editing identification. In addition, the unique advantage of long reads allowed us to uncover novel insights about RNA editing occurrences in single molecules and double-stranded RNA (dsRNA) structures. L-GIREMI provides a valuable means to study RNA nucleotide variants in long-read RNA-seq.

Post-transcriptional regulation plays a key role in transcriptome reprogramming during cardiac pathogenesis. In previous studies, we have identified that cardiac enriched RNA-binding protein, RBFox1 plays key role in cardiac hypertrophy through mRNA alternative splicing regulation in nuclei. However, RBFox1 gene also generates a cytosolic isoform (RBFox1c), suggesting additional functions of post-transcriptional regulation in heart. In adult heart, RBFox1c mRNA constituted ~ 40% of total RBFox1 level but was significantly repressed in pressure-overloaded failing mouse heart. Using CRISPR-Cas9 technology, we have established an isoform specific RBFox1c-cKO mouse. At baseline inactivation of RBFox1c led to decreased cardiac function along with induction of cardiac fibrosis. RBFox1c-cKO mice also showed macrophages infiltration into myocardium post 7days MI. In contrast, restoration of RBFox1c expression in adult intact hearts significantly reduced cardiac fibrosis post stress. RNA-seq analyses in RBFox1c expressing cardiomyocytes showed that RBFox1c specifically suppressed the expression of pro-inflammatory genes. Secondly, CLIP-Seq analysis and targeted RNA-IP showed that RBFox1c could directly interact with inflammatory pathway mRNAs. These results suggested the inflammatory mRNAs are direct downstream targets regulated by RBFox1c. Using both in vitro cultured cardiomyocytes and intact mouse hearts, we demonstrated that expression of RBFox1c reduces pro-inflammatory mRNA expression at baseline and upon hypertrophy stimulation. Lastly, we characterized the interactome of RBFox1c through proteomic analysis and found RBFox1c specifically interacted with a component of the RNA NMD machinery-Upf1. RBFox1c interaction with Upf1 in cardiomyocytes was diminished upon hypertrophic stress. Furthermore, by inactivation of Upf1 via siRNA, we demonstrated that RBFox1c mediated repression of proinflammatory genes was Upf1 dependent.RBFox1 regulates cardiac transcriptome reprogramming in two post-transcriptional processes via distinct isoforms. While the RBFox1n regulates RNA splicing, the RBFox1c functions through targeted mRNA repression of proinflammatory genes by recruitment of Upf1 mediated RNA degradation.

Lipin 1 regulates cellular lipid homeostasis through roles in glycerolipid synthesis (through phosphatidic acid phosphatase activity) and transcriptional coactivation. Lipin 1-deficient individuals exhibit episodic disease symptoms that are triggered by metabolic stress, such as stress caused by prolonged fasting. We sought to identify critical lipin 1 activities during fasting. We determined that lipin 1 deficiency induces widespread alternative mRNA splicing in liver during fasting, much of which is normalized by refeeding. The role of lipin 1 in mRNA splicing was largely independent of its enzymatic function. We identified interactions between lipin 1 and spliceosome proteins, as well as a requirement for lipin 1 to maintain homeostatic levels of spliceosome small nuclear RNAs and specific RNA splicing factors. In fasted Lpin1-/- liver, we identified a correspondence between alternative splicing of phospholipid biosynthetic enzymes and dysregulated phospholipid levels; splicing patterns and phospholipid levels were partly normalized by feeding. Thus, lipin 1 influences hepatic lipid metabolism through mRNA splicing, as well as through enzymatic and transcriptional activities, and fasting exacerbates the deleterious effects of lipin 1 deficiency on metabolic homeostasis.

Motivation Alternative polyadenylation (polyA) sites near the 3′ end of a pre-mRNA create multiple mRNA transcripts with different 3′ untranslated regions (3′ UTRs). The sequence elements of a 3′ UTR are essential for many biological activities such as mRNA stability, sub-cellular localization, protein translation, protein binding and translation efficiency. Moreover, numerous studies in the literature have reported the correlation between diseases and the shortening (or lengthening) of 3′ UTRs. As alternative polyA sites are common in mammalian genes, several machine learning tools have been published for predicting polyA sites from sequence data. These tools either consider limited sequence features or use relatively old algorithms for polyA site prediction. Moreover, none of the previous tools consider RNA secondary structures as a feature to predict polyA sites. Results In this paper, we propose a new deep learning model, called DeepPASTA, for predicting polyA sites from both sequence and RNA secondary structure data. The model is then extended to predict tissue-specific polyA sites. Moreover, the tool can predict the most dominant (i.e. frequently used) polyA site of a gene in a specific tissue and relative dominance when two polyA sites of the same gene are given. Our extensive experiments demonstrate that DeepPASTA signisficantly outperforms the existing tools for polyA site prediction and tissue-specific relative and absolute dominant polyA site prediction. Availability and implementation https://github.com/arefeen/DeepPASTA Supplementary information Supplementary data are available at Bioinformatics online.

Transcriptomic analyses of postmortem brains have begun to elucidate molecular abnormalities in autism spectrum disorder (ASD). However, a crucial pathway involved in synaptic development, RNA editing, has not yet been studied on a genome-wide scale. Here we profiled global patterns of adenosine-to-inosine (A-to-I) editing in a large cohort of postmortem brains of people with ASD. We observed a global bias for hypoediting in ASD brains, which was shared across brain regions and involved many synaptic genes. We show that the Fragile X proteins FMRP and FXR1P interact with RNA-editing enzymes (ADAR proteins) and modulate A-to-I editing. Furthermore, we observed convergent patterns of RNA-editing alterations in ASD and Fragile X syndrome, establishing this as a molecular link between these related diseases. Our findings, which are corroborated across multiple data sets, including dup15q (genomic duplication of 15q11.2-13.1) cases associated with intellectual disability, highlight RNA-editing dysregulation in ASD and reveal new mechanisms underlying this disorder.

BackgroundRNA editing generates modifications to the RNA sequences, thereby increasing protein diversity and shaping various layers of gene regulation. Recent studies have revealed global shifts in editing levels across many cancer types, as well as a few specific mechanisms implicating individual sites in tumorigenesis or metastasis. However, most tumor-associated sites, predominantly in noncoding regions, have unknown functional relevance.ResultsHere, we carry out integrative analysis of RNA editing profiles between epithelial and mesenchymal tumors, since epithelial-mesenchymal transition is a key paradigm for metastasis. We identify distinct editing patterns between epithelial and mesenchymal tumors in seven cancer types using TCGA data, an observation further supported by single-cell RNA sequencing data and ADAR perturbation experiments in cell culture. Through computational analyses and experimental validations, we show that differential editing sites between epithelial and mesenchymal phenotypes function by regulating mRNA abundance of their respective genes. Our analysis of RNA-binding proteins reveals ILF3 as a potential regulator of this process, supported by experimental validations. Consistent with the known roles of ILF3 in immune response, epithelial-mesenchymal differential editing sites are enriched in genes involved in immune and viral processes. The strongest target of editing-dependent ILF3 regulation is the transcript encoding PKR, a crucial player in immune and viral response.ConclusionsOur study reports widespread differences in RNA editing between epithelial and mesenchymal tumors and a novel mechanism of editing-dependent regulation of mRNA abundance. It reveals the broad impact of RNA editing in cancer and its relevance to cancer-related immune pathways.

In the mammalian heart, the left ventricle (LV) rapidly becomes more dominant in size and function over the right ventricle (RV) after birth. The molecular regulators responsible for this chamber-specific differential growth are largely unknown. We found that cardiomyocytes in the neonatal mouse RV had lower proliferation, more apoptosis, and a smaller average size compared with the LV. This chamber-specific growth pattern was associated with a selective activation of p38 mitogen-activated protein kinase (MAPK) activity in the RV and simultaneous inactivation in the LV. Cardiomyocyte-specific deletion of both the Mapk14 and Mapk11 genes in mice resulted in loss of p38 MAPK expression and activity in the neonatal heart. Inactivation of p38 activity led to a marked increase in cardiomyocyte proliferation and hypertrophy but diminished cardiomyocyte apoptosis, specifically in the RV. Consequently, the p38-inactivated hearts showed RV-specific enlargement postnatally, progressing to pulmonary hypertension and right heart failure at the adult stage. Chamber-specific p38 activity was associated with differential expression of dual-specific phosphatases (DUSPs) in neonatal hearts, including DUSP26. Unbiased transcriptome analysis revealed that IRE1α/XBP1-mediated gene regulation contributed to p38 MAPK-dependent regulation of neonatal cardiomyocyte proliferation and binucleation. These findings establish an obligatory role of DUSP/p38/IRE1α signaling in cardiomyocytes for chamber-specific growth in the postnatal heart.

Alternative splicing is an RNA processing mechanism that affects most genes in human, contributing to disease mechanisms and phenotypic diversity. The regulation of splicing involves an intricate network of cis-regulatory elements and trans-acting factors. Due to their high sequence specificity, cis-regulation of splicing can be altered by genetic variants, significantly affecting splicing outcomes. Recently, multiple methods have been applied to understanding the regulatory effects of genetic variants on splicing. However, it is still challenging to go beyond apparent association to pinpoint functional variants. To fill in this gap, we utilized large-scale datasets of the Genotype-Tissue Expression (GTEx) project to study genetically-modulated alternative splicing (GMAS) via identification of allele-specific splicing events. We demonstrate that GMAS events are shared across tissues and individuals more often than expected by chance, consistent with their genetically driven nature. Moreover, although the allelic bias of GMAS exons varies across samples, the degree of variation is similar across tissues vs. individuals. Thus, genetic background drives the GMAS pattern to a similar degree as tissue-specific splicing mechanisms. Leveraging the genetically driven nature of GMAS, we developed a new method to predict functional splicing-altering variants, built upon a genotype-phenotype concordance model across samples. Complemented by experimental validations, this method predicted >1000 functional variants, many of which may alter RNA-protein interactions. Lastly, 72% of GMAS-associated SNPs were in linkage disequilibrium with GWAS-reported SNPs, and such association was enriched in tissues of relevance for specific traits/diseases. Our study enables a comprehensive view of genetically driven splicing variations in human tissues.

Motivation RNA-sequencing (RNA-seq) enables global identification of RNA-editing sites in biological systems and disease. A salient step in many studies is to identify editing sites that statistically associate with treatment (e.g. case versus control) or covary with biological factors, such as age. However, RNA-seq has technical features that incumbent tests (e.g. t-test and linear regression) do not consider, which can lead to false positives and false negatives. Results In this study, we demonstrate the limitations of currently used tests and introduce the method, RNA-editing tests (REDITs), a suite of tests that employ beta-binomial models to identify differential RNA editing. The tests in REDITs have higher sensitivity than other tests, while also maintaining the type I error (false positive) rate at the nominal level. Applied to the GTEx dataset, we unveil RNA-editing changes associated with age and gender, and differential recoding profiles between brain regions. Availability and implementation REDITs are implemented as functions in R and freely available for download at https://github.com/gxiaolab/REDITs. The repository also provides a code example for leveraging parallelization using multiple cores.

microRNAs (miRNAs) are small non-coding RNAs that play critical roles in gene regulation. The presence of miRNAs in extracellular biofluids is increasingly recognized. However, most previous characterization of extracellular miRNAs focused on their overall expression levels. Alternative sequence isoforms and modifications of miRNAs were rarely considered in the extracellular space. Here, we developed a highly accurate bioinformatic method, called miNTA, to identify 3’ non-templated additions (NTAs) of miRNAs using small RNA-sequencing data. Using miNTA, we conducted an in-depth analysis of miRNA 3’ NTA profiles in 1047 extracellular RNA-sequencing data sets of 4 types of biofluids. This analysis identified hundreds of miRNAs with 3’ uridylation or adenylation, with the former being more prevalent. Among these miRNAs, up to 53% (22%) had an average 3’ uridylation (adenylation) level of at least 10% in a specific biofluid. Strikingly, we found that 3’ uridylation levels enabled segregation of different types of biofluids, more effectively than overall miRNA expression levels. This observation suggests that 3’ NTA levels possess fluid-specific information relatively robust to batch effects. In addition, we observed that extracellular miRNAs with 3’ uridylations are enriched in processes related to angiogenesis, apoptosis, and inflammatory response, and this type of modification may stabilize base-pairing between miRNAs and their target genes. Together, our study provides a comprehensive landscape of miRNA NTAs in human biofluids, which paves way for further biomarker discoveries. The insights generated in our work built a foundation for future functional, mechanistic, and translational discoveries.

Recent studies revealed global shifts in RNA editing, the modification of RNA sequences, across many cancers. Besides a few sites implicated in tumorigenesis or metastasis, most tumor-associated sites, predominantly in noncoding regions, have unknown function. Here, we characterize editing profiles between epithelial (E) and mesenchymal (M) phenotypes in seven cancer types, as epithelial-mesenchymal transition (EMT) is a key paradigm for metastasis. We observe distinct editing patterns between E and M tumors and EMT induction upon loss of ADAR enzymes in cultured cells. E-M differential sites are highly enriched in genes involved in immune and viral processes, some of which regulate mRNA abundance of their respective genes. We identify a novel mechanism in which ILF3 preferentially stabilizes edited transcripts. Among editing-dependent ILF3 targets is the transcript encoding PKR, a crucial player in immune response. Our study demonstrates the broad impact of RNA editing in cancer and relevance of editing to cancer-related immune pathways.

The essential biological cofactor heme is synthesized in cells in the Fe(II) form. Oxidized Fe(III) heme is specifically required for processing primary transcripts of microRNAs (pri-miRNAs) by the RNA-binding protein DGCR8, a core component of the Microprocessor complex. It is unknown how readily available Fe(III) heme is in the largely reducing environment in human cells and how changes in cellular Fe(III) heme availability alter microRNA (miRNA) expression. Here we address the first question by characterizing DGCR8 mutants with various degrees of deficiency in heme-binding. We observed a strikingly simple correlation between Fe(III) heme affinityin vitroand the Microprocessor activity in HeLa cells, with the heme affinity threshold for activation estimated to be between 0.6-5 pM under typical cell culture conditions. The threshold is strongly influenced by cellular heme synthesis and uptake. We suggest that the threshold reflects a labile Fe(III) heme pool in cells. Based on our understanding of DGCR8 mutants, we reanalyzed recently reported miRNA sequencing data and conclude that heme is generally required for processing canonical pri-miRNAs, that heme modulates the specificity of Microprocessor, and that cellular heme level and differential DGCR8 heme occupancy alter the expression of distinct groups of miRNAs in a hierarchical fashion. Overall, our study provides the first glimpse of a labile Fe(III) heme pool important for a fundamental physiological function and reveal principles governing how Fe(III) heme modulates miRNA maturation at a genomic scale. We also discuss potential states and biological significance of the labile Fe(III) heme pool.

Alcohol (ethanol, EtOH) consumption during pregnancy can result in fetal alcohol spectrum disorders (FASDs), which are characterized by prenatal and postnatal growth restriction and craniofacial dysmorphology. Recently, cell-derived extracellular vesicles, including exosomes and microvesicles containing several species of RNAs (exRNAs), have emerged as a mechanism of cell-to-cell communication. However, EtOH’s effects on the biogenesis and function of non-coding exRNAs during fetal development have not been explored. Therefore, we studied the effects of maternal EtOH exposure on the composition of exosomal RNAs in the amniotic fluid (AF) using rat fetal alcohol exposure (FAE) model. Through RNA-Seq analysis we identified and verified AF exosomal miRNAs with differential expression levels specifically associated with maternal EtOH exposure. Uptake of purified FAE AF exosomes by rBMSCs resulted in significant alteration of molecular markers associated with osteogenic differentiation of rBMSCs. We also determined putative functional roles for AF exosomal miRNAs (miR-199a-3p, miR-214-3p and let-7g) that are dysregulated by FAE in osteogenic differentiation of rBMSCs. Our results demonstrate that FAE alters AF exosomal miRNAs and that exosomal transfer of dysregulated miRNAs has significant molecular effects on stem cell regulation and differentiation. Our results further suggest the usefulness of assessing molecular alterations in AF exRNAs to study the mechanisms of FAE teratogenesis that should be further investigated by using an in vivo model.

BACKGROUND It was recently discovered that abundant and stable extracellular RNA (exRNA) species exist in bodily fluids. Saliva is an emerging biofluid for biomarker development for noninvasive detection and screening of local and systemic diseases. Use of RNA-Sequencing (RNA-Seq) to profile exRNA is rapidly growing; however, no single preparation and analysis protocol can be used for all biofluids. Specifically, RNA-Seq of saliva is particularly challenging owing to high abundance of bacterial contents and low abundance of salivary exRNA. Given the laborious procedures needed for RNA-Seq library construction, sequencing, data storage, and data analysis, saliva-specific and optimized protocols are essential. METHODS We compared different RNA isolation methods and library construction kits for long and small RNA sequencing. The role of ribosomal RNA (rRNA) depletion also was evaluated. RESULTS The miRNeasy Micro Kit (Qiagen) showed the highest total RNA yield (70.8 ng/mL cell-free saliva) and best small RNA recovery, and the NEBNext library preparation kits resulted in the highest number of detected human genes [5649–6813 at 1 reads per kilobase RNA per million mapped (RPKM)] and small RNAs [482–696 microRNAs (miRNAs) and 190–214 other small RNAs]. The proportion of human RNA-Seq reads was much higher in rRNA-depleted saliva samples (41%) than in samples without rRNA depletion (14%). In addition, the transfer RNA (tRNA)-derived RNA fragments (tRFs), a novel class of small RNAs, were highly abundant in human saliva, specifically tRF-4 (4%) and tRF-5 (15.25%). CONCLUSIONS Our results may help in selection of the best adapted methods of RNA isolation and small and long RNA library constructions for salivary exRNA studies.

Motivation The length of the 3′ untranslated region (3′ UTR) of an mRNA is essential for many biological activities such as mRNA stability, sub-cellular localization, protein translation, protein binding and translation efficiency. Moreover, correlation between diseases and the shortening (or lengthening) of 3′ UTRs has been reported in the literature. This length is largely determined by the polyadenylation cleavage site in the mRNA. As alternative polyadenylation (APA) sites are common in mammalian genes, several tools have been published recently for detecting APA sites from RNA-Seq data or performing shortening/lengthening analysis. These tools consider either up to only two APA sites in a gene or only APA sites that occur in the last exon of a gene, although a gene may generally have more than two APA sites and an APA site may sometimes occur before the last exon. Furthermore, the tools are unable to integrate the analysis of shortening/lengthening events with APA site detection. Results We propose a new tool, called TAPAS, for detecting novel APA sites from RNA-Seq data. It can deal with more than two APA sites in a gene as well as APA sites that occur before the last exon. The tool is based on an existing method for finding change points in time series data, but some filtration techniques are also adopted to remove change points that are likely false APA sites. It is then extended to identify APA sites that are expressed differently between two biological samples and genes that contain 3′ UTRs with shortening/lengthening events. Our extensive experiments on simulated and real RNA-Seq data demonstrate that TAPAS outperforms the existing tools for APA site detection or shortening/lengthening analysis significantly. Availability and implementation https://github.com/arefeen/TAPAS Supplementary information Supplementary data are available at Bioinformatics online.

RNA editing introduces single nucleotide changes to RNA, thus potentially diversifying gene expression. Recent studies have reported significant changes in RNA editing profiles in disease and development. The functional consequences of these widespread alterations remain elusive because of the unknown function of most RNA editing sites. Here, we carry out a comprehensive analysis of A-to-I editomes in human populations. Surprisingly, we observe highly similar editing profiles across populations despite striking differences in the expression levels of ADAR genes. Striving to explain this discrepancy, we uncover a functional mechanism of A-to-I editing in regulating mRNA abundance. We show that A-to-I editing stabilizes RNA secondary structures and reduces the accessibility of AGO2-miRNA to target sites in mRNAs. The editing-dependent stabilization of mRNAs in turn alters the observed editing levels in the stable RNA repertoire. Our study provides valuable insights into the functional impact of RNA editing in human cells., RNA A-to-I editing introduces single nucleotide changes to RNA, but its role in cells remains unclear. Here, the authors analyse A-to-I editomes in human samples and find that A-to-I editing stabilizes RNA secondary structures and reduces the accessibility of AGO2-miRNA to target sites in mRNAs.

Alternative transcription start (ATS) and alternative polyadenylation (APA) create alternative RNA isoforms and modulate many aspects of RNA expression and protein production. However, ATS and APA remain difficult to detect in RNA sequencing (RNA-seq). Here, we developed mountainClimber, a de novo cumulative-sum-based approach to identify ATS and APA as change points. Unlike many existing methods, mountainClimber runs on a single sample and identifies multiple ATS or APA sites anywhere in the transcript. We analyzed 2,342 GTEx samples (36 tissues, 215 individuals) and found that tissue type is the predominant driver of transcript end variations. 75% and 65% of genes exhibited differential APA and ATS across tissues, respectively. In particular, testis displayed longer 5' untranslated regions (UTRs) and shorter 3' UTRs, often in genes related to testis-specific biology. Overall, we report the largest study of transcript ends across human tissues to our knowledge. mountainClimber is available at github.com/gxiaolab/mountainClimber.

© 2019 Elsevier Inc. The Extracellular RNA Communication Consortium (ERCC) was launched to accelerate progress in the new field of extracellular RNA (exRNA) biology and to establish whether exRNAs and their carriers, including extracellular vesicles (EVs), can mediate intercellular communication and be utilized for clinical applications. Phase 1 of the ERCC focused on exRNA/EV biogenesis and function, discovery of exRNA biomarkers, development of exRNA/EV-based therapeutics, and construction of a robust set of reference exRNA profiles for a variety of biofluids. Here, we present progress by ERCC investigators in these areas, and we discuss collaborative projects directed at development of robust methods for EV/exRNA isolation and analysis and tools for sharing and computational analysis of exRNA profiling data. The Extracellular RNA Communication Consortium (ERCC) presents progress toward understanding the biology of extracellular RNAs and their use as biomarkers and therapeutics.

Allele-specific protein-RNA binding is an essential aspect that may reveal functional genetic variants (GVs) mediating post-transcriptional regulation. Recently, genome-wide detection of in vivo binding of RNA-binding proteins is greatly facilitated by the enhanced crosslinking and immunoprecipitation (eCLIP) method. We developed a new computational approach, called BEAPR, to identify allele-specific binding (ASB) events in eCLIP-Seq data. BEAPR takes into account crosslinking-induced sequence propensity and variations between replicated experiments. Using simulated and actual data, we show that BEAPR largely outperforms often-used count analysis methods. Importantly, BEAPR overcomes the inherent overdispersion problem of these methods. Complemented by experimental validations, we demonstrate that the application of BEAPR to ENCODE eCLIP-Seq data of 154 proteins helps to predict functional GVs that alter splicing or mRNA abundance. Moreover, many GVs with ASB patterns have known disease relevance. Overall, BEAPR is an effective method that helps to address the outstanding challenge of functional interpretation of GVs., Differential binding of RNA-binding proteins mediated by genetic variants (GVs) can influence posttranscriptional regulation. Here, the authors develop BEAPR, a computational approach to identify allele-specific binding events in eCLIP-Seq data.

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

Identification of functional genetic variants and elucidation of their regulatory mechanisms represent significant challenges of the post-genomic era. A poorly understood topic is the involvement of genetic variants in mediating post-transcriptional RNA processing, including alternative splicing. Thus far, little is known about the genomic, evolutionary, and regulatory features of genetically modulated alternative splicing (GMAS). Here, we systematically identified intronic tag variants for genetic modulation of alternative splicing using RNA-seq data specific to cellular compartments. Combined with our previous method that identifies exonic tags for GMAS, this study yielded 622 GMAS exons. We observed that GMAS events are highly cell type independent, indicating that splicing-altering genetic variants could have widespread function across cell types. Interestingly, GMAS genes, exons, and single-nucleotide variants (SNVs) all demonstrated positive selection or accelerated evolution in primates. We predicted that GMAS SNVs often alter binding of splicing factors, with SRSF1 affecting the most GMAS events and demonstrating global allelic binding bias. However, in contrast to their GMAS targets, the predicted splicing factors are more conserved than expected, suggesting that cis-regulatory variation is the major driving force of splicing evolution. Moreover, GMAS-related splicing factors had stronger consensus motifs than expected, consistent with their susceptibility to SNV disruption. Intriguingly, GMAS SNVs in general do not alter the strongest consensus position of the splicing factor motif, except the more than 100 GMAS SNVs in linkage disequilibrium with polymorphisms reported by genome-wide association studies. Our study reports many GMAS events and enables a better understanding of the evolutionary and regulatory features of this phenomenon.

Male vertebrate social displays vary from physically simple to complex, with the latter involving exquisite motor command of the body and appendages. Studies of these displays have, in turn, provided substantial insight into neuromotor mechanisms. The neotropical golden-collared manakin (Manacus vitellinus) has been used previously as a model to investigate intricate motor skills because adult males of this species perform an acrobatic and androgen-dependent courtship display. To support this behavior, these birds express elevated levels of androgen receptors (AR) in their skeletal muscles. Here we use RNA sequencing to explore how testosterone (T) modulates the muscular transcriptome to support male manakin courtship displays. In addition, we explore how androgens influence gene expression in the muscles of the zebra finch (Taenopygia guttata), a model passerine bird with a limited courtship display and minimal muscle AR. We identify androgen-dependent, muscle-specific gene regulation in both species. In addition, we identify manakin-specific effects that are linked to muscle use during the manakin display, including androgenic regulation of genes associated with muscle fiber contractility, cellular homeostasis, and energetic efficiency. Overall, our results point to numerous genes and gene networks impacted by androgens in male birds, including some that underlie optimal muscle function necessary for performing acrobatic display routines. Manakins are excellent models to explore gene regulation promoting athletic ability.

Autism spectrum disorder (ASD) is a genetically complex, clinically heterogeneous neurodevelopmental disease. Recently, our understanding of the molecular abnormalities in ASD has been expanded through transcriptomic analyses of postmortem brains. However, a crucial molecular pathway involved in synaptic development, RNA editing, has not yet been studied on a genome-wide scale. Here, we profiled the global patterns of adenosine-to-inosine (A-to-I) editing in a large cohort of post-mortem ASD brains. Strikingly, we observed a global bias of hypo-editing in ASD brains, common to different brain regions and involving many genes with known neurobiological functions. Through genome-wide protein-RNA binding analyses and detailed molecular assays, we show that the Fragile X proteins, FMRP and FXR1P, interact with ADAR proteins and modulate A-to-I editing. Furthermore, we observed convergent patterns of RNA editing alterations in ASD and Fragile X syndrome, thus establishing RNA editing as a molecular link underlying these two highly related diseases. Our findings support a role for RNA editing dysregulation in ASD and highlight novel mechanisms for RNA editing regulation.

© International & American Associations for Dental Research 2018. Coevolution of the human host and its associated microbiota has led to sophisticated interactions to maintain a delicate homeostasis. Emerging evidence suggests that in addition to small molecules, peptides, and proteins, small regulatory noncoding RNAs (sRNAs) might play an important role in cross-domain interactions. In this study, we revealed the presence of diverse host transfer RNA–derived small RNAs (tsRNAs) among human salivary sRNAs. We selected 2 tsRNAs (tsRNA-000794 and tsRNA-020498) for further study based on their high sequence similarity to specific tRNAs from a group of Gram-negative oral bacteria, including Fusobacterium nucleatum, a key oral commensal and opportunistic pathogen. We showed that the presence of F. nucleatum triggers exosome-mediated release of tsRNA-000794 and tsRNA-020498 by human normal oral keratinocyte cells. Furthermore, both tsRNA candidates exerted a growth inhibition effect on F. nucleatum, likely through interference with bacterial protein biosynthesis, but did not affect the growth of Streptococcus mitis, a health-associated oral Gram-positive bacterium whose genome does not carry sequences bearing high similarity to either tsRNA. Our data provide the first line of evidence for the modulatory role of host-derived tsRNAs in the microbial-host interaction.

Background: RBFox1 is known to be an RNA splicing regulator with enriched expression in cardiac muscle. Loss of RBFox1 expression is a molecular hallmark associated with heart failure in mouse and human. Genetic manipulation of RBFox1 reveals a major function of RBFox1 in the pathogenesis of cardiac hypertrophy, fibrosis and dysfunction under pathological stresses. However, much of our current knowledge about RBFox1 focuses on the nuclear RBFox1 with a major impact on global alternative splicing changes. Yet, RBFox1 gene also generates a cytosolic isoform-RBFox1c through alternative splicing, and the specific function of RBFox1c has not been characterized. Goal and Methods: This study investigated the functional impact and the underlying mechanism of the RBFox1c in cardiac pathological remodeling and targeted gene regulation. Results: RBFox1c expression is significantly repressed in stressed cardiomyocytes in vivo and in vitro. RNA-Seq combined with IPA analyses from cardiomyocytes revealed RBFox1c but not the nucleus RBFox1 specifically suppressed pro-inflammatory genes. In the cardiac specific RBFox1 knockout mice, enhanced cardiac fibrosis is observed following I/R injury associated with elevated expression of pro-inflammatory genes. In contrast, cardiac specific expression of RBFox1c reduced cardiac fibrosis and inflammatory gene expression following pressure-overload and myocardial infarct injury associated with improved ejection fraction. Motif enrichment analysis identified significant enrichment of the RBFox1 binding motif in the 3’UTR of the RBFox1c regulated genes. We performed CLIP analysis followed by RT-PCR and observed RBFox1c interacted with inflammatory gene 3’UTR. Lastly, we explored the interactome of RBFox1c and found RBFox1c specifically interacted with a component of RNA decay machinery-Upf1. Indeed, while expression of RBFox1c repressed inflammatory gene expression, inactivation of Upf1 abolished this effect in cardiomyocytes. Conclusion: RBFox1 regulates cardiac transcriptome reprogramming at two post-transcriptional steps. The RBFox1 nuclei isoform regulates RNA splicing reprogramming, and the RBFox1c represses proinflammatory genes via recruitment of Upf1 mediated RNA degradation.