Accurate transcriptome-wide identification and quantification of alternative polyadenylation from RNA-seq data with APAIQ.

Alternative polyadenylation (APA) enables a gene to generate multiple transcripts with different 3' ends, which is dynamic across different cell types or conditions. Many computational methods have been developed to characterize sample-specific APA using the corresponding RNA-seq data, but suffered from high error rate on both polyadenylation site (PAS) identification and quantification of PAS usage (PAU), and bias toward 3' untranslated regions. Here we developed a tool for APA identification and quantification (APAIQ) from RNA-seq data, which can accurately identify PAS and quantify PAU in a transcriptome-wide manner. Using 3' end-seq data as the benchmark, we showed that APAIQ outperforms current methods on PAS identification and PAU quantification, including DaPars2, Aptardi, mountainClimber, SANPolyA, and QAPA. Finally, applying APAIQ on 421 RNA-seq samples from liver cancer patients, we identified >540 tumor-associated APA events and experimentally validated two intronic polyadenylation candidates, demonstrating its capacity to unveil cancer-related APA with a large-scale RNA-seq data set.

3'UTR heterogeneity and cancer progression.

The majority of human mRNAs generate alternative 3' untranslated regions (UTRs) through various processes, including RNA modifications such as RNA editing, m6A methylation, and alternative polyadenylation (APA), with 3'UTR splicing as an emerging mechanism. Multiple factors, ranging from the genome to transcriptome level, regulate these processes and contribute to 3'UTR heterogeneity. Genomic variants in 3'UTR regions as well as aberrant 3'UTR processing alter the transcriptomic landscape and are associated with cancer. Increasing evidence, aided by high-resolution sequencing technologies and large-scale computational analyses, points towards potential crosstalk between these processes, whose deregulation may further contribute to cancer pathogenesis. In-depth characterization of these events will increase our appreciation of their significance and help to drive therapeutic development in this field.

Pan-cancer pervasive upregulation of 3' UTR splicing drives tumourigenesis.

Most mammalian genes generate messenger RNAs with variable untranslated regions (UTRs) that are important post-transcriptional regulators. In cancer, shortening at 3' UTR ends via alternative polyadenylation can activate oncogenes. However, internal 3' UTR splicing remains poorly understood as splicing studies have traditionally focused on protein-coding alterations. Here we systematically map the pan-cancer landscape of 3' UTR splicing and present this in SpUR ( http://www.cbrc.kaust.edu.sa/spur/home/ ). 3' UTR splicing is widespread, upregulated in cancers, correlated with poor prognosis and more prevalent in oncogenes. We show that antisense oligonucleotide-mediated inhibition of 3' UTR splicing efficiently reduces oncogene expression and impedes tumour progression. Notably, CTNNB1 3' UTR splicing is the most consistently dysregulated event across cancers. We validate its upregulation in hepatocellular carcinoma and colon adenocarcinoma, and show that the spliced 3' UTR variant is the predominant contributor to its oncogenic functions. Overall, our study highlights the importance of 3' UTR splicing in cancer and may launch new avenues for RNA-based anti-cancer therapeutics.

Global analysis of RNA-binding proteins identifies a positive feedback loop between LARP1 and MYC that promotes tumorigenesis.

In addition to genomic alterations, aberrant changes in post-transcriptional regulation can modify gene function and drive cancer development. RNA-binding proteins (RBPs) are a large class of post-transcriptional regulators that have been increasingly implicated in carcinogenesis. By integrating multi-omics data, we identify LARP1 as one of the most upregulated RBPs in colorectal cancer (CRC) and demonstrate its oncogenic properties. We perform LARP1:RNA interactome profiling and unveil a previously unexplored role for LARP1 in targeting the 3'UTR of oncogenes in CRC. Notably, we identify the proto-oncogenic transcription factor MYC as a key LARP1-regulated target. Our data show that LARP1 positively modulates MYC expression by associating with its 3'UTR. In addition, antisense oligonucleotide-mediated blocking of the interaction between LARP1 and the MYC 3'UTR reduces MYC expression and in vitro CRC growth. Furthermore, a systematic analysis of LARP1:protein interactions reveals IGF2BP3 and YBX1 as LARP1-interacting proteins that also regulate MYC expression and CRC development. Finally, we demonstrate that MYC reciprocally modulates LARP1 expression by targeting its enhancer. In summary, our data reveal a critical, previously uncharacterized role of LARP1 in promoting CRC tumorigenesis, validate its direct regulation of the proto-oncogene MYC and delineate a model of the positive feedback loop between MYC and LARP1 that promotes CRC growth and development.

MiR-138 is a potent regulator of the heterogenous MYC transcript population in cancers.

3'UTR shortening in cancer has been shown to activate oncogenes, partly through the loss of microRNA-mediated repression. This suggests that many reported microRNA-oncogene target interactions may not be present in cancer cells. One of the most well-studied oncogenes is the transcription factor MYC, which is overexpressed in more than half of all cancers. MYC overexpression is not always accompanied by underlying genetic aberrations. In this study, we demonstrate that the MYC 3'UTR is shortened in colorectal cancer (CRC). Using unbiased computational and experimental approaches, we identify and validate microRNAs that target the MYC coding region. In particular, we show that miR-138 inhibits MYC expression and suppresses tumor growth of CRC and hepatocellular carcinoma (HCC) cell lines. Critically, the intravenous administration of miR-138 significantly impedes MYC-driven tumor growth in vivo. Taken together, our results highlight the previously uncharacterized shortening of the MYC 3'UTR in cancer, and identify miR-138 as a potent regulator of the heterogenous MYC transcript population.

Systematic Analysis of Intronic miRNAs Reveals Cooperativity within the Multicomponent FTX Locus to Promote Colon Cancer Development.

Approximately half of all miRNA reside within intronic regions and are often cotranscribed with their host genes. However, most studies of intronic miRNA focus on individual miRNA, while conversely most studies of protein-coding and noncoding genes frequently ignore any intron-derived miRNA. We hypothesize that the individual components of such multigenic loci may play cooperative or competing roles in driving disease progression and that examining the combinatorial effect of these components would uncover deeper insights into their functional importance. To address this, we performed systematic analyses of intronic miRNA:host loci in colon cancer. The FTX locus, comprising of a long noncoding RNA FTX and multiple intronic miRNA, was highly upregulated in cancer, and cooperativity within this multicomponent locus promoted cancer growth. FTX interacted with DHX9 and DICER and regulated A-to-I RNA editing and miRNA expression. These results show for the first time that a long noncoding RNA can regulate A-to-I RNA editing, further expanding the functional repertoire of long noncoding RNA. Intronic miR-374b and miR-545 inhibited tumor suppressors PTEN and RIG-I to enhance proto-oncogenic PI3K-AKT signaling. Furthermore, intronic miR-421 may exert an autoregulatory effect on miR-374b and miR-545. Taken together, our data unveil the intricate interplay between intronic miRNA and their host transcripts in the modulation of key signaling pathways and disease progression, adding new perspectives to the functional landscape of multigenic loci. SIGNIFICANCE: This study illustrates the functional relationships between individual components of multigenic loci in regulating cancer progression.See related commentary by Calin, p. 1212.

Noncoding RNA:RNA Regulatory Networks in Cancer.

Noncoding RNAs (ncRNAs) constitute the majority of the human transcribed genome. This largest class of RNA transcripts plays diverse roles in a multitude of cellular processes, and has been implicated in many pathological conditions, especially cancer. The different subclasses of ncRNAs include microRNAs, a class of short ncRNAs; and a variety of long ncRNAs (lncRNAs), such as lincRNAs, antisense RNAs, pseudogenes, and circular RNAs. Many studies have demonstrated the involvement of these ncRNAs in competitive regulatory interactions, known as competing endogenous RNA (ceRNA) networks, whereby lncRNAs can act as microRNA decoys to modulate gene expression. These interactions are often interconnected, thus aberrant expression of any network component could derail the complex regulatory circuitry, culminating in cancer development and progression. Recent integrative analyses have provided evidence that new computational platforms and experimental approaches can be harnessed together to distinguish key ceRNA interactions in specific cancers, which could facilitate the identification of robust biomarkers and therapeutic targets, and hence, more effective cancer therapies and better patient outcome and survival.

A FTH1 gene:pseudogene:microRNA network regulates tumorigenesis in prostate cancer.

Non-coding RNAs play a vital role in diverse cellular processes. Pseudogenes, which are non-coding homologs of protein-coding genes, were once considered non-functional evolutional relics. However, recent studies have shown that pseudogene transcripts can regulate their parental transcripts by sequestering shared microRNAs (miRNAs), thus acting as competing endogenous RNAs (ceRNAs). In this study, we utilize an unbiased screen to identify the ferritin heavy chain 1 (FTH1) transcript and multiple FTH1 pseudogenes as targets of several oncogenic miRNAs in prostate cancer (PCa). We characterize the critical role of this FTH1 gene:pseudogene:miRNA network in regulating tumorigenesis in PCa, whereby oncogenic miRNAs downregulate the expression of FTH1 and its pseudogenes to drive oncogenesis. We further show that impairing miRNA binding and subsequent ceRNA crosstalk completely rescues the slow growth phenotype in vitro and in vivo. Our results also demonstrate the reciprocal regulation between the pseudogenes and intracellular iron levels, which are crucial for multiple physiological and pathophysiological processes. In summary, we describe an extensive gene:pseudogene network comprising multiple miRNAs and multiple pseudogenes derived from a single parental gene. The network could be regulated through multiple mechanisms to modulate iron storage in various signaling pathways, the deregulation of which results in PCa development and progression.

Screening for protein-protein interactions using Förster resonance energy transfer (FRET) and fluorescence lifetime imaging microscopy (FLIM).

We present a high content multiwell plate cell-based assay approach to quantify protein interactions directly in cells using Förster resonance energy transfer (FRET) read out by automated fluorescence lifetime imaging (FLIM). Automated FLIM is implemented using wide-field time-gated detection, typically requiring only 10 s per field of view (FOV). Averaging over biological, thermal and shot noise with 100's to 1000's of FOV enables unbiased quantitative analysis with high statistical power. Plotting average donor lifetime vs. acceptor/donor intensity ratio clearly identifies protein interactions and fitting to double exponential donor decay models provides estimates of interacting population fractions that, with calibrated donor and acceptor fluorescence intensities, can yield dissociation constants. We demonstrate the application to identify binding partners of MST1 kinase and estimate interaction strength among the members of the RASSF protein family, which have important roles in apoptosis via the Hippo signalling pathway. KD values broadly agree with published biochemical measurements.

PLCɛ and the RASSF family in tumour suppression and other functions.

Not all proteins implicated in direct binding to Ras appear to have a positive role in the generation and progression of tumours; examples include Phospholipase C epsilon (PLCɛ) and some members of the Ras-association domain family (RASSF). The RASSF family comprises of ten members, known as RASSF1 to RASSF10. PLCɛ and RASSF members carry a common Ras-association domain (RA) that can potentially bind Ras oncoproteins and mediate Ras-regulated functions. RASSF1 to RASSF6 also share a common SARAH domain that facilitates protein-protein interactions with other SARAH domain proteins. The majority of the family are frequently downregulated by epigenetic silencing in cancers. They are implicated in various important biological processes including apoptosis, microtubule stabilisation and cell cycle regulation. Recent studies have reinforced the tumour suppressive properties of the RASSF family, with new evidence of emerging pathways and novel functions that suggest a wider role for these proteins. This review will first describe an emerging role of PLCɛ in tumour suppression and then focus on and summarise the new findings on the RASSF family in the last five years to consolidate their well-established functions, and highlight the new regulatory roles of specific RASSF members.

Comparative analysis of interactions of RASSF1-10.

Members of the RASSF family (RASSF1-10) have been identified as candidate tumour suppressors that are frequently downregulated by promoter hypermethylation in cancers. These proteins carry a common Ras-association (RA) and SARAH domain (RASSF1-6) that can potentially bind Ras oncoproteins and mediate protein-protein interactions with other SARAH domain proteins. However, there is a notable lack of comparative characterisation of the RASSF family, as well as molecular and structural information that facilitate their tumour suppressive functions. As part of our comparative analysis, we modelled the RA and SARAH domains of the RASSF members based on existing structures and predicted their potential interactions. These in silico predictions were compared to in vitro interaction studies with Ras and MST kinase (a SARAH domain-containing protein). Our data shows a diversity of interaction within the RASSF family RA domain, whereas the SARAH domain-mediated interactions for RASSF1-6 are consistent with the predictions. This suggests that different members, despite shared general architecture, could have distinct functional properties. Additionally, we identify a new interacting partner for MST kinase in the form of RASSF7. Current data supports an interaction model where RASSF serves as an adaptor for the assembly of multiple protein complexes and further functional interactions, involving MST kinases and other SARAH domain proteins, which could be regulated by Ras.