Massively parallel screen uncovers many rare 3' UTR variants regulating mRNA abundance of cancer driver genes.

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.

Multifaceted role of RNA editing in promoting loss-of-function of PODXL in cancer.

PODXL, a protein that is dysregulated in multiple cancers, plays an important role in promoting cancer metastasis. In this study, we report that RNA editing promotes the inclusion of a PODXL alternative exon. The resulting edited PODXL long isoform is more prone to protease digestion and has the strongest effects on reducing cell migration and cisplatin chemoresistance among the three PODXL isoforms (short, unedited long, and edited long isoforms). Importantly, the editing level of the PODXL recoding site and the inclusion level of the PODXL alternative exon are strongly associated with overall patient survival in Kidney Renal Clear Cell Carcinoma (KIRC). Supported by significant enrichment of exonic RNA editing sites in alternatively spliced exons, we hypothesize that exonic RNA editing sites may enhance proteomic diversity through alternative splicing, in addition to amino acid changes, a previously under-appreciated aspect of RNA editing function.

RNA editing in cancer impacts mRNA abundance in immune response pathways.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.

Allele-specific binding of RNA-binding proteins reveals functional genetic variants in the RNA.

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.

Regulation of RNA editing by RNA-binding proteins in human cells.

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.

Widespread RNA editing dysregulation in brains from autistic individuals.

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.

RNA editing in nascent RNA affects pre-mRNA splicing.

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.

Genomic analysis of ADAR1 binding and its involvement in multiple RNA processing pathways.

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.

Identification of allele-specific alternative mRNA processing via transcriptome sequencing.

Establishing the functional roles of genetic variants remains a significant challenge in the post-genomic era. Here, we present a method, allele-specific alternative mRNA processing (ASARP), to identify genetically influenced mRNA processing events using transcriptome sequencing (RNA-Seq) data. The method examines RNA-Seq data at both single-nucleotide and whole-gene/isoform levels to identify allele-specific expression (ASE) and existence of allele-specific regulation of mRNA processing. We applied the methods to data obtained from the human glioblastoma cell line U87MG and primary breast cancer tissues and found that 26-45% of all genes with sufficient read coverage demonstrated ASE, with significant overlap between the two cell types. Our methods predicted potential mechanisms underlying ASE due to regulations affecting either whole-gene-level expression or alternative mRNA processing, including alternative splicing, alternative polyadenylation and alternative transcriptional initiation. Allele-specific alternative splicing and alternative polyadenylation may explain ASE in hundreds of genes in each cell type. Reporter studies following these predictions identified the causal single nucleotide variants (SNVs) for several allele-specific alternative splicing events. Finally, many genes identified in our study were also reported as disease/phenotype-associated genes in genome-wide association studies. Future applications of our approach may provide ample insights for a better understanding of the genetic basis of gene regulation underlying phenotypic diversity and disease mechanisms.

Accurate identification of A-to-I RNA editing in human by transcriptome sequencing.

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.

Resveratrol-induced apoptosis is mediated by early growth response-1, Krüppel-like factor 4, and activating transcription factor 3.

Resveratrol, a dietary phytoalexin readily available in the diet, is reported to possess antitumorigenic properties in several cancers, including colorectal. However, the underlying mechanism(s) involved is not completely understood. In the present study, we investigated the effect of resveratrol treatment on gene modulation in human colorectal cancer cells and identified activating transcription factor 3 (ATF3) as the most highly induced gene after treatment. We confirmed that resveratrol upregulates ATF3 expression, both at the mRNA and protein level, and showed resveratrol involvement in ATF3 transcriptional regulation. Analysis of the ATF3 promoter revealed the importance of early growth response-1 (Egr-1; located at -245 to -236) and Krüppel-like factor 4 (KLF4; located at -178 to -174) putative binding sites in resveratrol-mediated ATF3 transactivation. Specificity of these sites to the Egr-1 and KLF4 protein was confirmed by electrophoretic mobility shift and chromatin immunoprecipitation assays. Resveratrol increased Egr-1 and KLF4 expression, which preceded ATF3 expression, and further suggests Egr-1 and KLF4 involvement in resveratrol-mediated activity. We provide evidence for Egr-1 and KLF4 interaction in the presence of resveratrol, which may facilitate ATF3 transcriptional regulation by this compound. Furthermore, we demonstrate that induction of apoptosis by resveratrol is mediated, in part, by increased ATF3 expression. Taken together, these results provide a novel mechanism by which resveratrol induces ATF3 expression and represent an additional explanation of how resveratrol exerts its antitumorigenic effects in human colorectal cancer cells.

Induction of cell growth arrest by atmospheric non-thermal plasma in colorectal cancer cells.

Plasma is generated by ionizing neutral gas molecules, resulting in a mixture of energy particles, including electrons and ions. Recent progress in the understanding of non-thermal atmospheric plasma has led to applications in biomedicine. However, the exact molecular mechanisms involved in plasma-induced cell growth arrest are unclear. In this study, we investigated the feasibility of non-thermal atmospheric plasma treatment for cancer therapy and examined the mechanism by which plasma induces anti-proliferative properties and cell death in human colorectal cancer cells. Non-thermal atmospheric plasma induced cell growth arrest and induced apoptosis. In addition, plasma reduced cell migration and invasion activities. As a result, we found that plasma treatment to the cells increases ß-catenin phosphorylation, suggesting that ß-catenin degradation plays a role at least in part in plasma-induced anti-proliferative activity. Therefore, non-thermal atmospheric plasma constitutes a new biologic tool with the potential for therapeutic applications that modulate cell signaling and function.

Effects of atmospheric nonthermal plasma on invasion of colorectal cancer cells.

The effect that the gas content and plasma power of atmospheric, nonthermal plasma has on the invasion activity in colorectal cancer cells has been studied. Helium and helium plus oxygen plasmas were induced through a nozzle and operated with an ac power of less than 10 kV which exhibited a length of 2.5 cm and a diameter of 3-4 mm in ambient air. Treatment of cancer cells with the plasma jet resulted in a decrease in cell migrationinvasion with higher plasma intensity and the addition of oxygen to the He flow gas.

ESE-1/EGR-1 pathway plays a role in tolfenamic acid-induced apoptosis in colorectal cancer cells.

Nonsteroidal anti-inflammatory drugs (NSAIDs) are known to prevent colorectal tumorigenesis. Although antitumor effects of NSAIDs are mainly due to inhibition of cyclooxygenase activity, there is increasing evidence that cyclooxygenase-independent mechanisms may also play an important role. The early growth response-1 (EGR-1) gene is a member of the immediate-early gene family and has been identified as a tumor suppressor gene. Tolfenamic acid is a NSAID that exhibits anticancer activity in a pancreatic cancer model. In the present study, we investigated the anticancer activity of tolfenamic acid in human colorectal cancer cells. Tolfenamic acid treatment inhibited cell growth and induced apoptosis as measured by caspase activity and bioelectric impedance. Tolfenamic acid induced EGR-1 expression at the transcription level, and analysis of the EGR-1 promoter showed that a putative ETS-binding site, located at -400 and -394 bp, was required for activation by tolfenamic acid. The electrophoretic mobility shift assay and chromatin immunoprecipitation assay confirmed that this sequence specifically bound to the ETS family protein epithelial-specific ETS-1 (ESE-1) transcription factor. Tolfenamic acid also facilitated translocation of endogenous and exogenous ESE-1 to the nucleus in colorectal cancer cells, and gene silencing using ESE-1 small interfering RNA attenuated tolfenamic acid-induced EGR-1 expression and apoptosis. Overexpression of EGR-1 increased apoptosis and decreased bioelectrical impedance, and silencing of endogenous EGR-1 prevented tolfenamic acid-induced apoptosis. These results show that activation of ESE-1 via enhanced nuclear translocation mediates tolfenamic acid-induced EGR-1 expression, which plays a critical role in the activation of apoptosis.

Changes in pyridoxal kinase immunoreactivity in the gerbil hippocampus following spontaneous seizure.

To identify the roles of pyridoxal kinase (PLK) in epileptogenesis and the recovery mechanisms in spontaneous seizure, a chronological and comparative analysis of PLK expression in the gerbil hippocampus was conducted. PLK immunoreactivity in a pre-seizure group of seizure sensitive (SS) gerbils was more strongly detected than that in a seizure resistant (SR) group. The density of PLK immunoreactivity in a 30-min postictal group was significantly lower than that of a pre-seizure group. In a 12 h postictal group, PLK immunodensity recovered to pre-seizure level. The over-expression of PLK in the hippocampus of pre-seizure SS gerbils suggests that PLP play an important role in the modulation of GAD activity and GABA reuptake as mediated by membrane transporter via neurons.

The decreases in calcium binding proteins and neurofilament immunoreactivities in the Purkinje cell of the seizure sensitive gerbils.

In previous studies, it has been reported that Purkinje cell degeneration during seizure is evoked by excitotoxicity due to an increase in the intracellular Ca(2+) level, though calbindin D-28k (CB) and parvalbumin (PV), intracellular free calcium buffers, are abundantly colocalized in these cells. In the present study, we investigated the expressions of CB, PV, neurofilament (NF) 68, 150, 200, and polyphosphorylated epitope in NF (RT 97), in the cerebellum of gerbils to identify the mechanism of Purkinje cell damages induced by seizure. In seizure resistant gerbils, nearly all the Purkinje cells showed CB, PA, NF 150, NF 200 and RT 97 immunoreactivity. In SS gerbils, however, a clear decrease in the number of CB(+) and PV(+) Purkinje cells was observed. The NF and RT 97 immunoreactivities, in the Purkinje cells was also lower (except NF 68), but not absent. These results suggest several points. First, the decrease in the concentrations of CB and PV may render the Purkinje cells more susceptible to intermittent Ca(2+) fluctuations and more prone to accumulating intolerable quantities of Ca(2+). Second, during the Ca(2+)-PV interaction PV plays an important role in facilitating donations of Mg(2+), which is a potent enzyme activator in phosphorylation. Thus the decline in PV concentration also implicated the defects of phosphorylation in the NF. Third, increases in both the intracellular Ca(2+) level and dephosphorylation trigger the degradation of the NF, particularly NF 200. Finally, these degradations in the NF induce the functional defects in Purkinje cell, which then cause Purkinje cell degeneration.