Extracellular RNAs (exRNAs) in human body fluids are emerging as effective biomarkers for detection of diseases. Saliva, as the most accessible and noninvasive body fluid, has been shown to harbor ...exRNA biomarkers for several human diseases. However, the entire spectrum of exRNA from saliva has not been fully characterized.
Using high-throughput RNA sequencing (RNA-Seq), we conducted an in-depth bioinformatic analysis of noncoding RNAs (ncRNAs) in human cell-free saliva (CFS) from healthy individuals, with a focus on microRNAs (miRNAs), piwi-interacting RNAs (piRNAs), and circular RNAs (circRNAs).
Our data demonstrated robust reproducibility of miRNA and piRNA profiles across individuals. Furthermore, individual variability of these salivary RNA species was highly similar to those in other body fluids or cellular samples, despite the direct exposure of saliva to environmental impacts. By comparative analysis of >90 RNA-Seq data sets of different origins, we observed that piRNAs were surprisingly abundant in CFS compared with other body fluid or intracellular samples, with expression levels in CFS comparable to those found in embryonic stem cells and skin cells. Conversely, miRNA expression profiles in CFS were highly similar to those in serum and cerebrospinal fluid. Using a customized bioinformatics method, we identified >400 circRNAs in CFS. These data represent the first global characterization and experimental validation of circRNAs in any type of extracellular body fluid.
Our study provides a comprehensive landscape of ncRNA species in human saliva that will facilitate further biomarker discoveries and lay a foundation for future studies related to ncRNAs in human saliva.
RNA editing enhances the diversity of gene products at the post-transcriptional level. Approaches for genome-wide identification of RNA editing face two main challenges: separating true editing sites ...from false discoveries and accurate estimation of editing levels. We developed an approach to analyze transcriptome sequencing data (RNA-seq) for global identification of RNA editing in cells for which whole-genome sequencing data are available. We applied the method to analyze RNA-seq data of a human glioblastoma cell line, U87MG. Around 10,000 DNA-RNA differences were identified, the majority being putative A-to-I editing sites. These predicted A-to-I events were associated with a low false-discovery rate (∼5%). Moreover, the estimated editing levels from RNA-seq correlated well with those based on traditional clonal sequencing. Our results further facilitated unbiased characterization of the sequence and evolutionary features flanking predicted A-to-I editing sites and discovery of a conserved RNA structural motif that may be functionally relevant to editing. Genes with predicted A-to-I editing were significantly enriched with those known to be involved in cancer, supporting the potential importance of cancer-specific RNA editing. A similar profile of DNA-RNA differences as in U87MG was predicted for another RNA-seq data set obtained from primary breast cancer samples. Remarkably, significant overlap exists between the putative editing sites of the two transcriptomes despite their difference in cell type, cancer type, and genomic backgrounds. Our approach enabled de novo identification of the RNA editome, which sets the stage for further mechanistic studies of this important step of post-transcriptional regulation.
Adenosine deaminases acting on RNA (ADARs) are the primary factors underlying adenosine to inosine (A-to-I) editing in metazoans. Here we report the first global study of ADAR1-RNA interaction in ...human cells using CLIP-seq. A large number of CLIP sites are observed in Alu repeats, consistent with ADAR1's function in RNA editing. Surprisingly, thousands of other CLIP sites are located in non-Alu regions, revealing functional and biophysical targets of ADAR1 in the regulation of alternative 3' UTR usage and miRNA biogenesis. We observe that binding of ADAR1 to 3' UTRs precludes binding by other factors, causing 3' UTR lengthening. Similarly, ADAR1 interacts with DROSHA and DGCR8 in the nucleus and possibly out-competes DGCR8 in primary miRNA binding, which enhances mature miRNA expression. These functions are dependent on ADAR1's editing activity, at least for a subset of targets. Our study unfolds a broad landscape of the functional roles of ADAR1.
In eukaryotes, nascent RNA transcripts undergo an intricate series of RNA processing steps to achieve mRNA maturation. RNA editing and alternative splicing are two major RNA processing steps that can ...introduce significant modifications to the final gene products. By tackling these processes in isolation, recent studies have enabled substantial progress in understanding their global RNA targets and regulatory pathways. However, the interplay between individual steps of RNA processing, an essential aspect of gene regulation, remains poorly understood. By sequencing the RNA of different subcellular fractions, we examined the timing of adenosine-to-inosine (A-to-I) RNA editing and its impact on alternative splicing. We observed that >95% A-to-I RNA editing events occurred in the chromatin-associated RNA prior to polyadenylation. We report about 500 editing sites in the 3' acceptor sequences that can alter splicing of the associated exons. These exons are highly conserved during evolution and reside in genes with important cellular function. Furthermore, we identified a second class of exons whose splicing is likely modulated by RNA secondary structures that are recognized by the RNA editing machinery. The genome-wide analyses, supported by experimental validations, revealed remarkable interplay between RNA editing and splicing and expanded the repertoire of functional RNA editing sites.
BACE1 is the rate-limiting enzyme for amyloid-β peptides (Aβ) generation, a key event in the pathogenesis of Alzheimer’s disease (AD). By an unknown mechanism, levels of BACE1 and a BACE1 ...mRNA-stabilizing antisense RNA (BACE1-AS) are elevated in the brains of AD patients, implicating that dysregulation of BACE1 expression plays an important role in AD pathogenesis. We found that nuclear factor erythroid-derived 2-related factor 2 (NRF2/NFE2L2) represses the expression of BACE1 and BACE1-AS through binding to antioxidant response elements (AREs) in their promoters of mouse and human. NRF2-mediated inhibition of BACE1 and BACE1-AS expression is independent of redox regulation. NRF2 activation decreases production of BACE1 and BACE1-AS transcripts and Aβ production and ameliorates cognitive deficits in animal models of AD. Depletion of NRF2 increases BACE1 and BACE1-AS expression and Aβ production and worsens cognitive deficits. Our findings suggest that activation of NRF2 can prevent a key early pathogenic process in AD.
RNA 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.
Here, 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.
Our 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.
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.
Understanding the function of rare non-coding variants represents a significant challenge. Using MapUTR, a screening method, we studied the function of rare 3' UTR variants affecting mRNA abundance ...post-transcriptionally. Among 17,301 rare gnomAD variants, an average of 24.5% were functional, with 70% in cancer-related genes, many in critical cancer pathways. This observation motivated an interrogation of 11,929 somatic mutations, uncovering 3928 (33%) functional mutations in 155 cancer driver genes. Functional MapUTR variants were enriched in microRNA- or protein-binding sites and may underlie outlier gene expression in tumors. Further, we introduce untranslated tumor mutational burden (uTMB), a metric reflecting the amount of somatic functional MapUTR variants of a tumor and show its potential in predicting patient survival. Through prime editing, we characterized three variants in cancer-relevant genes (MFN2, FOSL2, and IRAK1), demonstrating their cancer-driving potential. Our study elucidates the function of tens of thousands of non-coding variants, nominates non-coding cancer driver mutations, and demonstrates their potential contributions to cancer.
Autism spectrum disorder (ASD) is a heterogeneously pervasive developmental disorder in which various genetic and environmental factors are believed to underlie its development. Recently, epigenetics ...has been suggested as a novel concept for ASD aetiology with a proposition that epigenetic marks can be transgenerationally inherited. Based on this assumption of epigenetics, we investigated the transgenerational inheritance of ASD-like behaviours and their related synaptic changes in the VPA animal model of ASD. The first generation (F1) VPA-exposed offspring exhibited autistic-like impaired sociability and increased marble burying. They also showed increased seizure susceptibility, hyperactivity and decreased anxiety. We mated the VPA-exposed F1 male offspring with naïve females to produce the second generation (F2), and then similarly mated the F2 to deliver the third generation (F3). Remarkably, the autism-like behavioural phenotypes found in F1 persisted to the F2 and F3. Additionally, the frontal cortices of F1 and F3 showed some imbalanced expressions of excitatory/inhibitory synaptic markers, suggesting a transgenerational epigenetic inheritance. These results open the idea that E/I imbalance and ASD-like behavioural changes induced by environmental insults in mice can be epigenetically transmitted, at least, to the third generation. This study could help explain the unprecedented increase in ASD prevalence.
Drosophila neuropeptide F (NPF), a homolog of vertebrate neuropeptide Y, functions in feeding and coordination of behavioral changes in larvae and in modulation of alcohol sensitivity in adults, ...suggesting diverse roles for this peptide. To gain more insight into adult-specific NPF neuronal functions, we studied how npf expression is regulated in the adult brain. Here, we report that npf expression is regulated in both sex-nonspecific and male-specific manners. Our data show that male-specific npf (ms-npf) expression is controlled by the transformer (tra)-dependent sex-determination pathway. Furthermore, fruitless, one of the major genes functioning downstream of tra, is apparently an upstream regulator of ms-npf transcription. Males lacking ms-npf expression (through tra(F)-mediated feminization) or npf-ablated male flies display significantly reduced male courtship activity, suggesting that one function of ms-npf neurons is to modulate fruitless-regulated sexual behavior. Interestingly, one of the ms-npf neuronal groups belongs to the previously defined clock-controlling dorsolateral neurons. Such ms-npf expression in the dorsolateral neurons is absent in arrhythmic Clock(Jrk) and cycle⁰² mutants, suggesting that npf is under dual regulation by circadian and sex-determining factors. Based on these data, we propose that NPF also plays a role in clock-controlled sexual dimorphism in adult Drosophila.