Competing endogenous RNA networks: tying the essential knots for cancer biology and therapeutics
Journal of Hematology & Oncology volume 8, Article number: 30 (2015)
A recently discovered dimension of post-transcriptional gene regulation involves co-regulatory crosstalk between RNA transcripts, which compete for common pools of microRNA (miRNA) molecules. These competing endogenous RNAs (ceRNAs), or natural miRNA sponges, have an active role in regulating miRNA availability within the cell and form intertwined regulatory networks. Recent reports have implicated diverse RNA species including protein-coding messenger RNAs and non-coding RNAs as ceRNAs in human development and diseases including human cancer. In this review, we discuss the most recent discoveries that implicate natural miRNA decoys in human cancer biology, as well as exciting advances in the study of ceRNA networks and dynamics. The structure and topology of intricate genome-scale ceRNA networks can be predicted computationally, and their dynamic response to fluctuations in ceRNA and miRNA levels can be studied via mathematical modeling. Additionally, the development of new methods to quantitatively determine absolute expression levels of miRNA and ceRNA molecules have expanded the capacity to accurately study the efficiency of ceRNA crosstalk in diverse biological models. These major milestones are of critical importance to identify key components of ceRNA regulatory networks that could aid the development of new approaches to cancer diagnostics and oligonucleotide-based therapeutics.
MicroRNAs (miRNAs) comprise a class of small non-coding RNAs implicated in post-transcriptional RNA regulation. These RNA molecules are approximately 22 nt in length and bind to miRNA response elements (MREs) on target transcripts through sequence complementarity, usually inducing repression or destabilization of the target transcript [1,2]. A single miRNA can regulate up to thousands of target transcripts , and miRNAs can act in a combinatorial manner by binding to separate MREs on a single RNA transcript. MiRNA-mediated regulation is one of the most widespread post-transcriptional regulatory mechanism in eukaryotes, and is estimated to affect the majority of human transcripts [4,5]. In fact, accumulating evidence highlights the role of miRNA-mediated regulation in cell growth, differentiation, proliferation, and apoptosis . Furthermore, miRNAs have been implicated in many diseases including cancer [7,8]. Alterations in the miRNA balance in the cell can lead to dysregulation of tumor suppressor genes and/or oncogenes controlled by aberrantly expressed miRNAs, leading to cancer [9,10].
Although the miRNA field is relatively young, a miRNA-based therapeutic has already entered phase 2 clinical trials [11,12]. In cancer, therapeutic strategies targeting miRNAs can be divided into two opposing approaches: (1) MiRNA mimics, which restore the function of tumor suppressive miRNAs; and (2) MiRNA inhibitors, which antagonize the function of oncogenic miRNAs. This can be achieved either by using chemically-modified RNA oligos which mimic or inhibit miRNA function or by artificially generated ‘miRNA sponges’ [13-15]. These exogenously expressed transcripts contain tandem repeats of MRE sites that allow them to specifically bind a distinct miRNA or combination of miRNAs. Once the miRNA(s) of interest has been decoyed, it will be unavailable to bind to its targets, leading to the effective de-repression of these transcripts [13,16].
In a similar fashion to how artificial miRNA sponges regulate the expression of endogenous miRNAs, it has been recently demonstrated that endogenous transcripts which contain MREs can interact and co-regulate each other through competition for common miRNA pools in the cell, acting as endogenous miRNA sponges or competing endogenous RNAs (ceRNAs) [17,18]. This intricate interplay between diverse RNA species, whereby one transcript can reciprocally modulate the expression of another transcript by sequestering shared miRNAs, has been referred to as ceRNA crosstalk . CeRNA activity has been attributed to both protein-coding and non-coding RNA transcripts. Examples of this competing endogenous crosstalk have been described in plants, such as Arabidopsis thaliana , the primate virus Herpesvirus saimiri , zebrafish , and mice as well as humans [23,24].
The implication of ceRNAs in humans has been extended from normal cell differentiation to tissue development and pathology, since relevant ceRNA functionality has been associated with the self-renewal capacity of embryonic stem cells , muscle differentiation , idiopathic pulmonary fibrosis , and liver cancer  among other biological processes. Although ceRNA research is in its infancy, accumulating evidence suggests that this additional dimension of post-transcriptional gene regulation, in which RNA transcripts titrate miRNA availability, represents a biologically relevant, well-conserved and widespread mechanism of regulation in eukaryotes.
Competing endogenous RNAs in cancer
Several examples of MRE-containing non-coding RNA transcripts acting as ceRNAs have been described. This suggests that many long non-coding RNAs may potentially function as ceRNAs, and that the analysis of ceRNA interactions may facilitate the functional characterization of such non-coding transcripts. Interestingly, there is a widespread prevalence of somatic mutations, gross genetic alterations, and chromosome rearrangements affecting the non-coding genome in cancer [28,29]. Additionally, cancer cells often express alternative isoforms of mRNAs with shorter 3′UTR regions . These events may lead to the mutation or loss of MREs, which will not only affect miRNA regulation in cis but also the associated ceRNA interactions. The dysregulation of ceRNA interactions may, in turn, have profound implications for cancer initiation, maintenance, or progression. Thus, it is perhaps unsurprising that significant research efforts thus far on ceRNA-mediated post-transcriptional regulation have focused on cancer initiation, maintenance, or progression (Figure 1). The role of ceRNAs in cancer biology has been reviewed in detail elsewhere [19,31]. As such, in this review, we will focus mainly on the most recent developments in the field.
The most extensively characterized ceRNA network (ceRNET) is perhaps that of the tumor suppressor gene PTEN. The processed pseudogene PTENP1, which is highly homologous to PTEN and shares many of its 3′UTR MREs, was the first non-coding transcript shown to post-transcriptionally regulate PTEN expression in a miRNA-dependent manner . In a recent study, PTENP1 was reported to be downregulated in clear-cell renal cell carcinoma and shown to regulate PTEN expression. PTENP1 tumor suppressive properties were due in part to sequestration of miR-21, providing additional support for the functional relevance of PTENP1 as a tumor suppressive ceRNA for PTEN . Additional studies have identified various mRNA transcripts that function as PTEN ceRNAs and regulate PTEN expression in a 3′UTR and miRNA-dependent manner in brain, breast, colon, prostate, and skin cancers [33-36].
Recent reports have provided increasing support for the ceRNA activity of long non-coding RNAs (lncRNAs), suggesting that the characterization of ceRNA networks may enable researchers to systematically study the function of these largely uncharacterized RNA transcripts. In addition, lncRNAs often display specific expression patterns in different cell types, tissue types, developmental stages, and disease , making them ideal candidates for participating in finely tuned post-transcriptional regulation.
One such example is linc-RoR, which had already been shown to control self-renewal and maintain pluripotency of human embryonic stem cells by acting as a miR-145 ‘sponge’ and thus controlling OCT4, NANOG, and SOX2 expression . Similarly, linc-RoR was found to regulate expression of the same transcription factors in endometrial cancer stem cells, and inhibit their differentiation, in a miR-145 dependent manner .
Another example of lncRNA-mediated ceRNA regulation involves the tumor suppressor gene BARD1. The newly discovered lncRNA BARD1 9′L is transcribed from an alternative intronic promoter of the BARD1 gene, and was found to share both miR-203 and miR-101 MREs with BARD1 mRNA in their homologous 3′UTRs . Moreover, full-length BARD1 and cancer-associated BARD1 mRNAs were downregulated by miR-203 and miR-101, and lncRNA BARD1 9′L counteracted the effect of these miRNAs. These data support a role for BARD1 9′L as a tumor suppressor transcript through its ceRNA activity.
Conversely, the ceRNA activity of lncRNAs has also been shown to have an oncogenic effect: The lncRNA HOTAIR was shown to display ceRNA activity in gastric cancer cells, in which it was found to specifically bind the tumor suppressor miR-331-3p, modulating HER2 de-repression . In addition, the authors observed an association of HOTAIR expression with cancer progression and malignancy in clinical samples from patients with advanced gastric cancer.
A very recent report explored the role of the BRAF pseudogene (BRAFP1) in regulating the expression of BRAF in a miRNA-dependent manner . This study demonstrated that BRAFP1 functions as a ceRNA for BRAF in part by sponging miR-30a, miR-182, and miR-876, and that the oncogenic effect of BRAFP1 overexpression in vitro was at least partly mediated by the upregulation of BRAF expression and subsequent activation of MAPK. Moreover, the authors presented a mouse model overexpressing full-length murine B-Raf pseudogene that developed an aggressive malignancy resembling human diffuse large B cell lymphoma. This study is the first one to demonstrate the oncogenicity of a pseudogene in an engineered mouse model, and provides further in vivo support for the role of ceRNA-mediated miRNA sequestration on cancer development.
Molecular considerations for ceRNA interactions
Although there are increasing reports describing bona fide ceRNA interactions in diverse biological settings, several open questions remain about the physiological relevance of these interactions . Factors such as miRNA- and ceRNA-expression levels and subcellular localization, number of shared MREs, and miRNA/ceRNA binding affinity have all been suggested to modulate the effectiveness of ceRNA crosstalk [18,43,44]. Indeed, it has been shown that molecular interactions that are dependent on titration mechanisms, such as ceRNA crosstalk, usually display a threshold-like behavior determined by the relative abundance of interacting molecules .
Several recent reports have described various computational models that have shed light on the specific molecular requirements for effective ceRNA interactions both at steady-state conditions, and in dynamic conditions when the system transiently responds to perturbations around the threshold level [46-50]. Using an in silico mathematical mass-action model to determine the optimal conditions for ceRNA interactions, Ala et al. identified several molecular determinants for effective ceRNA crosstalk . These determinants include factors such as miRNA and ceRNA expression levels, the number of MREs shared between ceRNAs, and the total number of MREs for the miRNAs of interest. In an alternative approach, Figliuzzi et al. utilized a minimal rate-equation-based model of post-transcriptional regulation at steady-state levels to characterize the nature of ceRNA interactions . Interestingly, they suggest that ceRNA interactions can be either symmetrical, whereby two ceRNAs co-regulate each other, or asymmetrical, whereby one ceRNA regulates expression of the other in a unidirectional manner. In addition to steady-state system assumptions, researchers have harnessed various mathematical models to investigate the dynamics of ceRNA crosstalk. Studies have described the equilibrium phenomenon that the system displays in the threshold proximity that is found at near equimolar conditions , as well as the transient response of the system to perturbations . In addition, a recent report presented a computational model to quantitatively describe a minimum ceRNA network, and experimentally validated the predictions in cultured human cells using synthetic gene circuits .
Despite the differences in the mathematical approaches, interaction parameters, and assumptions for these studies, several common conclusions can be drawn. For instance, optimal ceRNA-mediated regulation was observed to depend heavily on the miRNA/ceRNA ratio [46,47,50], consistent with quantitative experimental data demonstrating that the optimal conditions for miRNA functionality depended largely on miRNA/target relative abundance [50,51]. Optimal ceRNA crosstalk, whereby one ceRNA has the most striking effect on its interacting ceRNA partners, was predicted to occur when miRNA and ceRNA levels were near equimolarity [47,51] (Figure 2A). Another common finding was that the number of shared miRNAs among RNA transcripts played an important role in determining the effectiveness of ceRNA crosstalk [46-48].
Interestingly, ceRNA crosstalk appeared to represent a viable mechanism to achieve fast, positive changes in transcript levels . Additionally, mathematical models predicted that ceRNA crosstalk could represent a powerful noise-processing mechanism in gene expression that can lead to noise reduction or amplification [47,48]. These bioinformatic analyses suggest that ceRNAs may potentially play a central role in post-transcriptional regulation.
These mathematical models also led to some unexpected findings. For example, ceRNAs may interact with other ceRNAs directly and/or indirectly . In the hypothetical network depicted in Figure 2B, ceRNA 1 interacts directly with ceRNA 2 via miRNA 1, and directly with ceRNA 3 via miRNA 2. In a situation where ceRNA 2 also interacts with ceRNA 3 via an additional miRNA 3 (Figure 2B, bottom panel), then ceRNA 1 would interact indirectly with ceRNA 2 via miRNA 2-ceRNA 3-miRNA 3. Ala et al. found that such indirect interactions were surprisingly effective at enhancing ceRNA crosstalk, and that ceRNAs were able to respond to changes in other ceRNAs even in the absence of shared miRNAs [46,47]. The relevance of indirect interactions to overall ceRNET architecture was corroborated by an independent study analyzing the interactions between distant ceRNET nodes in HEK293 cells . Additionally, in optimal molecular conditions, ceRNETs may be altered by changes in the transcription rate of a single key ceRNA , and this effect is predicted to be greater for highly expressed ceRNAs . These findings have enormous implications not only for the molecular study of cancer gene dysregulation but also for the exploration of miRNA sponges as new therapeutic tools.
As an essential complement to the above-mentioned analyses, which were largely computational, a number of recent studies have sought to experimentally quantify miRNA and target cellular abundance as well as miRNA target affinity. In a landmark study, a quantitative assessment of genome-wide miRNA regulation in the context of ceRNA crosstalk was performed in embryonic and mesenchymal stem cells, using Ago2 individual-nucleotide resolution UV crosslinking and immunoprecipitation, together with absolute miRNA and mRNA quantification . They found that miRNA/ceRNA ratios determined not only Ago accumulation to its different affinity sites but also the susceptibility and threshold of miRNA-dependent ceRNA regulation. Single-cell miRNA reporter assays confirmed that low miRNA/ceRNA pool ratios, such as those found endogenously for the miR-92/25 family, possess an increased susceptibility for ceRNA regulation. On the other hand, highly expressed miRNA families such as miR-294 and let-7 were found to be insensitive to changes in relative ceRNA abundance in this experimental model. These data establish an experimental framework for the biochemical analysis of ceRNA interactions and provide critical experimental support for the physiological relevance of ceRNA crosstalk.
As discussed above, several computational and experimental studies have suggested that high abundance miRNAs may not be susceptible to effective ceRNA regulation. However, exceptions to this hypothesis may exist. The HMGA2 mRNA was shown to decoy the highly abundant let-7 miRNA family, thus regulating TGFBR3 expression and enhancing TGF-β signaling . This ceRNA activity of the HMGA2 transcript was proposed to promote lung cancer progression. Critically, the authors found that HMGA2 and TGFBR3 transcripts were expressed at similar levels, and that total let-7 family expression was within an order of magnitude of these transcripts in the lung cancer cell line tested. These data suggest that in certain settings, even highly abundant miRNAs such as let-7 may be effectively regulated by ceRNA interactions. Furthermore, the precise topology of ceRNETs may result in bona fide ceRNA regulation which may seem improbable based purely on specific miRNA/mRNA ratios. For example, mathematical modeling has predicted that effective interactions between two ceRNAs can be mediated by a large number of shared miRNAs, which may each only have a weak effect on ceRNA levels . This suggests that future studies focused on high-throughput techniques to identify and characterize interconnected ceRNA networks will lead to important insights into the functional relevance of ceRNA interactions in complex physiological systems.
Indeed, the effectiveness of ceRNA crosstalk has been shown to vary in different physiological settings. For example, two recent studies examined the effect of sequestering the liver-specific miR-122. Denzler et al. analyzed the stoichiometric relationship between miR-122 and its targets in hepatocytes in both physiological and pathological conditions . Due to the high expression levels of miR-122, they found that effective ceRNA crosstalk in this system could only be mediated by non-physiological increases in ceRNA levels and that the system was insensitive to up to threefold changes in the miRNA level. However, even though ceRNA crosstalk was lacking in physiological conditions for highly abundant miR-122, the system showed effective ceRNA interaction near the threshold where miRNA and target levels were near equimolarity, in agreement with the general consensus in the field. In contrast, a second study by Luna et al. described a quantitative mathematical model of miR-122 sequestration by the hepatitis C virus (HCV) RNA, and utilized HITS-CLIP and single-cell measurements to demonstrate that functional sequestration of miR-122 by the HCV RNA reduced miR-122 binding to its endogenous target mRNAs .
These studies present accumulating in silico and in vitro evidence indicating that the overall stoichiometry of the interacting molecules and miRNA/ceRNA binding affinity are critical determinants of effective crosstalk between ceRNAs in vivo. These factors should be considered when investigating novel ceRNA transcripts and networks, especially those of relevance to specific human disease conditions, to facilitate the successful translation of research results to clinical settings.
Prediction of ceRNA networks in human cancer
Most ceRNA interactions reported to date have been between a pair of interacting transcripts. However, the ability of a single miRNA to bind and regulate multiple transcripts and the regulation of a single transcript by multiple different miRNA molecules allow the effects of relative changes in miRNA/ceRNA ratios to percolate through intertwined ceRNETs and affect the expression of distant ceRNAs. It is thus critical to begin to unravel the complexities of this regulatory crosstalk in genome-scale ceRNETs.
The availability of genome-wide miRNA-, mRNA-, and other transcript expression data in large cancer patient datasets, together with the development of novel algorithms and bioinformatic tools, enable the in silico prediction of ceRNA crosstalk. These predictions can be expanded, with the appropriate statistical treatment, to potentially generate genome-scale ceRNA network data for different biological processes, including cancer. Using this approach, breast and thyroid cancer-specific mRNA ceRNETs were recently described by two independent groups [57,58]. Although these results remain to be confirmed experimentally, they have proven useful to predict the risk of metastasis in breast cancer patients , and to highlight specific biological processes that may be central to certain cancers and regulated by ceRNA crosstalk, such as immune response regulation in papillary thyroid carcinoma .
Recently, a lncRNA-miRNA-mRNA ceRNET was predicted computationally in breast cancer patients . This report found that ceRNETs may be significantly altered between normal and pathological breast tissues, and identified the lncRNA PVT1 as a key ceRNA on breast cancer. A second ceRNET analysis involving lncRNAs was described using genome-wide expression data from gastric cancer patients . The authors integrated microarray data, bioinformatics, and a miRNA target database to generate a lncRNA-miRNA-mRNA network. Interestingly, this network contained a high number of mRNA targets implicated in cancer. Although relevant experimental validation will be needed to determine the implications of these findings, the authors showed co-regulation of the same ceRNAs in independent datasets from six other cancer types.
Recently, an intriguing new facet of ceRNET dysregulation in cancers was reported. Exome sequencing data from a large set of prostate cancer patients was analyzed to examine the effect of alternative cleavage on 3′UTRs and polyadenylation dynamics on putative ceRNET interactions . Cancer patients could be stratified into groups with differing risks according to their alternative polyadenylation profiles, and a ceRNA subnetwork enriched with prostate cancer genes was found to be specifically dysregulated only in high-risk cancers. As 3′UTRs harbor multiple MREs which mediate miRNA binding, it is perhaps unsurprising that 3′UTR shortening can affect the associated ceRNETs. However, the observation that an unusually large number of MREs are shared between these significantly altered 3′UTRs and the affected ceRNA interactions suggests that ceRNETs may indeed be selectively dysregulated in high-risk prostate cancer.
The potential of studying ceRNA crosstalk in the context of larger interconnected networks, rather than isolated ceRNA pair interactions, opens exciting possibilities for the study of ceRNA-mediated post-transcriptional regulation in cancer in a setting that is closer to physiological and pathological conditions. However, there are still several considerations limiting the practical applicability of ceRNET study. One main limitation is the dependence of many ceRNA prediction tools on miRNA target prediction algorithms. Despite improvements in the field, these algorithms are not fool-proof due to our incomplete understanding of the rules for miRNA target recognition. This is particularly relevant for the prediction of mammalian miRNA and subsequent ceRNA networks, as mammalian miRNAs often pair with imperfect complementarity to their target transcripts. Although ceRNA prediction analysis has been successfully performed using validated miRNA interactors in the context of PTEN , miRNA target prediction tools were still required to identify binding of these validated PTEN-targeting miRNAs to potential ceRNAs of interest.
Additionally, given the scale of high-throughput data and the number of interactions taken into account, careful statistical treatment of the data needs to be performed for global ceRNET study. Moreover, the discovery of novel ceRNA interactions rely, in many cases, on the significant co-regulation of transcripts in cancer patient expression data. However, it has been reported that significant expression correlation among ceRNAs can emerge in experimental readouts due to transcriptional fluctuations even in the absence of miRNA-mediated crosstalk .
The development of high-throughput experimental validation of ceRNA interactions would thus be central to the expansion of our knowledge about ceRNETs and applicability for clinical practice in the future. Recent advances in techniques enabling the high-throughput sequencing of immunoprecipitated RNAs after crosslinking (such as CLIP-Seq, HITS-CLIP, and PAR-CLIP) provide a biochemical method to identify relevant miRNA target interactions . For instance, the integration of HITS-CLIP data with computational approaches was shown to improve miRNA target prediction by over 20-fold . An integrated database termed starBase v2.0 has been developed to identify experimentally supported regulatory RNA-RNA and protein-RNA interaction networks from multiple CLIP-Seq datasets . starBase identified about 10,000 ceRNA pairs from CLIP-supported MREs, including miRNA-mRNA and miRNA-lncRNA interactions. As miRNAs, lncRNAs and even mRNAs often display tissue-, temporal-, and disease-specific expression patterns, such high-throughput biochemical analyses will be invaluable to characterize ceRNA interactions in specific cellular contexts.
In recent years, diverse RNA species including mRNAs, pseudogenes, small non-coding RNAs, lncRNAs, and circular RNAs have all been shown to possess ceRNA activity in defined biological settings. Given the vast amount of genome-wide transcript expression data, and advances in the computational management of this publicly available data, we anticipate that reports of the ceRNA function of various RNA transcripts will continue to increase. Since miRNA decoying capacity is independent of protein-coding function, it has been proposed that the ceRNA activity of coding transcripts may allow them to exert independent and even opposing roles to their encoded protein in initiating or maintaining tumorigenic properties . This in turn has important implications in the study, and potentially the clinical management, of human cancers. The existence of new ceRNA transcript classes also represents an exciting area for future study. Small RNA species such as rRNA and tRNA have been shown to interact with Argonaute proteins , raising the question of whether they could possess ceRNA activity.
In addition, ceRNA-mediated post-transcriptional regulation has opened the exciting possibility of methodically studying and functionalizing the non-coding transcriptome. In this respect, the identification of previously uncharacterized lncRNAs as ceRNA counterparts for mRNAs is flourishing. Interestingly, protein-coding genes that produce lncRNAs are over-represented in cancer association studies, being twice as likely to be associated to cancer as human protein-coding genes overall . Moreover, specific regulatory patterns for several lncRNA transcripts have been found for a number of biological processes . The ceRNA activity of non-coding transcripts may be enhanced by their escape from any form of interference from the translational machinery. All of this make lncRNAs great candidates for participating in finely tuned post-transcriptional regulation.
Most of the work done so far has focused on the study and quantification of key interacting RNA species such as miRNAs, mRNAs, lncRNAs, and other ceRNA transcripts. Future work investigating the potential interplay of ceRNAs with RNA binding proteins (RBPs) will provide further insights into the complexities of ceRNA crosstalk. Interestingly, the dynamic recruitment of Ago-interacting RNA-biding proteins has been shown to alter miRNA/target affinity [66,67]. A recent report has detailed the development of a database to identify potential RNA:RNA and RNA:protein regulatory interactions in multiple CLIP-Seq datasets encompassing miRNA, mRNA, lncRNA, circRNA, pseudogene, and RBP associations . Experimental characterization of these RNA:RBP interactions will be needed to fully understand their impact on the dynamics of ceRNA interactions.
One of the most exciting advances in the ceRNA field in recent times has been the potential ability to computationally determine the structure and topology of miRNA and ceRNA interactions in interconnected ceRNETs, and harness mathematical modeling to study the dynamic behavior of these ceRNETs in response to fluctuations in the expression of ceRNAs and miRNAs in the network. Furthermore, the quantitative determination of absolute expression levels of miRNA and ceRNA molecules expands the capacity to accurately study the efficiency of ceRNA crosstalk in diverse physiological and pathological conditions. We anticipate that the systematic analysis of ceRNET dynamics would allow for the discovery of new oligonucleotide gene therapies, targeting highly connected nodes, which, unlike non-essential players in the network, may cause severe cellular damage when disrupted .
A better understanding of the molecular conditions in which ceRNA interactions occur is of critical importance to facilitate the study of this dimension of post-transcriptional regulation. Fluctuations in both miRNA and ceRNA levels in vivo need to be carefully examined to determine whether this dysregulation will lead to significant perturbation of ceRNA crosstalk in specific biological settings. The integration of this information with genome-scale ceRNET dynamics will be important to determine which ceRNAs are highly connected target hubs, less likely to be susceptible to compensation by alternative signaling pathways.
Multiple mechanisms may contribute to the dysregulation of miRNA and associated ceRNA interactions in cancer. For example, single nucleotide polymorphisms (SNPs) located in specific miRNA binding sites have been shown to modulate the binding of miRNAs to these target sites. A panel of 30 proto-oncogene associated SNPs that may impair miRNA target interactions was recently identified, and a subset of these was further associated with therapeutic outcome in cancer patients . In addition to genomic alterations, sequence alterations at the transcript level such as those introduced through RNA editing have also been shown to effectively regulate miRNA target specification . The contribution of such events to the dysregulation of ceRNA networks in cancer and potential subsequent impact on tumor progression and metastasis should be investigated in future studies.
In conclusion, the study of miRNA and ceRNA networks will open up new avenues for basic cancer research, as well as facilitate the development of novel diagnostic and therapeutic tools. For example, the silencing of aberrantly expressed miRNAs in cancer has been achieved through antisense oligomers, and recent advances in targeted delivery to tumor cells have shown promise in mouse models in vivo . The identification of key miRNAs and ceRNAs in human cancer may thus represent promising new therapeutic targets. ceRNA network components could be targeted alone, to reduce oncogenic capacity of the cells, or in combination with traditional therapies to impair acquired drug resistance through compensatory pathways.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.
Miranda KC, Huynh T, Tay Y, Ang YS, Tam WL, Thomson AM, et al. A pattern-based method for the identification of MicroRNA binding sites and their corresponding heteroduplexes. Cell. 2006;126:1203–17.
Friedman RC, Farh KK, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009;19:92–105.
Thomas M, Lieberman J, Lal A. Desperately seeking microRNA targets. Nat Struct Mol Biol. 2010;17:1169–74.
Esquela-Kerscher A, Slack FJ. Oncomirs - microRNAs with a role in cancer. Nat Rev Cancer. 2006;6:259–69.
Calin GA, Dumitru CD, Shimizu M, Bichi R, Zupo S, Noch E, et al. Frequent deletions and down-regulation of micro- RNA genes miR15 and miR16 at 13q14 in chronic lymphocytic leukemia. Proc Natl Acad Sci U S A. 2002;99:15524–9.
Lu M, Zhang Q, Deng M, Miao J, Guo Y, Gao W, et al. An analysis of human microRNA and disease associations. PLoS One. 2008;3:e3420.
Calin GA, Croce CM. MicroRNA signatures in human cancers. Nat Rev Cancer. 2006;6:857–66.
Plaisier CL, Pan M, Baliga NS. A miRNA-regulatory network explains how dysregulated miRNAs perturb oncogenic processes across diverse cancers. Genome Res. 2012;22:2302–14.
Bader AG. miR-34 - a microRNA replacement therapy is headed to the clinic. Front Genet. 2012;3(120):1–9.
Beg M, Brenner A, Sachdev J, Borad M, Cortes J, Tibes R, et al. 4LBA A phase 1 study of first-in-class microRNA-34 mimic, MRX34, in patients with hepatocellular carcinoma or advanced cancer with liver metastasis. Eur J Can. 2014;50:196.
Ebert MS, Neilson JR, Sharp PA. MicroRNA sponges: competitive inhibitors of small RNAs in mammalian cells. Nat Methods. 2007;4:721–6.
Brown BD, Gentner B, Cantore A, Colleoni S, Amendola M, Zingale A, et al. Endogenous microRNA can be broadly exploited to regulate transgene expression according to tissue, lineage and differentiation state. Nat Biotechnol. 2007;25:1457–67.
Gentner B, Schira G, Giustacchini A, Amendola M, Brown BD, Ponzoni M, et al. Stable knockdown of microRNA in vivo by lentiviral vectors. Nat Methods. 2009;6:63–6.
Ebert MS, Sharp PA. MicroRNA sponges: progress and possibilities. RNA. 2010;16:2043–50.
Poliseno L, Salmena L, Zhang J, Carver B, Haveman WJ, Pandolfi PP. A coding-independent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010;465:1033–8.
Salmena L, Poliseno L, Tay Y, Kats L, Pandolfi PP. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146:353–8.
Tay Y, Rinn J, Pandolfi PP. The multilayered complexity of ceRNA crosstalk and competition. Nature. 2014;505:344–52.
Franco-Zorrilla JM, Valli A, Todesco M, Mateos I, Puga MI, Rubio-Somoza I, et al. Target mimicry provides a new mechanism for regulation of microRNA activity. Nat Genet. 2007;39:1033–7.
Cazalla D, Yario T, Steitz JA. Down-regulation of a host microRNA by a Herpesvirus saimiri noncoding RNA. Science. 2010;328:1563–6.
Memczak S, Jens M, Elefsinioti A, Torti F, Krueger J, Rybak A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333–8.
Cesana M, Cacchiarelli D, Legnini I, Santini T, Sthandier O, Chinappi M, et al. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. 2001;147:358–69.
Marques AC, Tan J, Lee S, Kong L, Heger A, Ponting CP. Evidence for conserved post-transcriptional roles of unitary pseudogenes and for frequent bifunctionality of mRNAs. Genome Biol. 2012;13:R102.
Wang Y, Xu Z, Jiang J, Xu C, Kang J, Xiao L, et al. Endogenous miRNA sponge lincRNA-RoR regulates Oct4, Nanog, and Sox2 in human embryonic stem cell self-renewal. Dev Cell. 2013;25:69–80.
Song X, Cao G, Jing L, Lin S, Wang X, Zhang J, et al. Analysing the relationship between lncRNA and protein-coding gene and the role of lncRNA as ceRNA in pulmonary fibrosis. J Cell Mol Med. 2014;18:991–1003.
Wang J, Liu X, Wu H, Ni P, Gu Z, Qiao Y, et al. CREB up-regulates long non-coding RNA, HULC expression through interaction with microRNA-372 in liver cancer. Nucleic Acids Res. 2010;38:5366–83.
Weinhold N, Jacobsen A, Schultz N, Sander C, Lee W. Genome-wide analysis of noncoding regulatory mutations in cancer. Nat Genet. 2014;46:1160–5.
Forbes SA, Beare D, Gunasekaran P, Leung K, Bindal N, Boutselakis H, et al. COSMIC: exploring the world’s knowledge of somatic mutations in human cancer. Nucleic Acids Res. 2015;43:D805–11.
Mayr C, Bartel DP. Widespread shortening of 3′UTRs by alternative cleavage and polyadenylation activates oncogenes in cancer cells. Cell. 2009;138:673–84.
Karreth FA, Pandolfi PP. ceRNA cross-talk in cancer: when ce-bling rivalries go awry. Cancer Discov. 2013;3:1113–21.
Yu G, Yao W, Gumireddy K, Li A, Wang J, Xiao W, et al. Pseudogene PTENP1 functions as a competing endogenous RNA to suppress clear cell renal cell carcinoma progression. Mol Cancer Ther. 2014; In press.
Lee DY, Jeyapalan Z, Fang L, Yang J, Zhang Y, Yee AY, et al. Expression of versican 3′-untranslated region modulates endogenous microRNA functions. PLoS One. 2010;5:e13599.
Tay Y, Kats L, Salmena L, Weiss D, Tan SM, Ala U, et al. Coding-independent regulation of the tumor suppressor PTEN by competing endogenous mRNAs. Cell. 2011;147:344–57.
Karreth FA, Tay Y, Perna D, Ala U, Tan SM, Rust AG, et al. In vivo identification of tumor- suppressive PTEN ceRNAs in an oncogenic BRAF-induced mouse model of melanoma. Cell. 2011;147:382–95.
Sumazin P, Yang X, Chiu HS, Chung WJ, Iyer A, Llobet-Navas D, et al. An extensive microRNA-mediated network of RNA–RNA interactions regulates established oncogenic pathways in glioblastoma. Cell. 2011;147:370–81.
Guttman M, Rinn JL. Modular regulatory principles of large non-coding RNAs. Nature. 2012;482:339–46.
Zhou X, Gao Q, Wang J, Zhang X, Liu K, Duan Z. Linc-RNA-RoR acts as a “sponge” against mediation of the differentiation of endometrial cancer stem cells by microRNA-145. Gynecol Oncol. 2014;133:333–9.
Pilyugin M, Irminger-Finger I. Long non-coding RNA and microRNAs might act in regulating the expression of BARD1 mRNAs. Int J Biochem Cell Biol. 2014;54:356–67.
Liu XH, Sun M, Nie FQ, Ge YB, Zhang EB, Yin DD, et al. Lnc RNA HOTAIR functions as a competing endogenous RNA to regulate HER2 expression by sponging miR-331-3p in gastric cancer. Mol Cancer. 2014;13:92.
Karreth FK, Reschke M, Ruocco A, Ng C, Chapuy B, Léopold V, et al. The BRAF pseudogene functions as a competitive endogenous RNA and induces lymphoma in vivo. Cell. 2015; In press.
Jens M, Rajewsky N. Competition between target sites of regulators shapes post-transcriptional gene regulation. Nat Rev Genet. 2015;16:113–26.
Mukherji S, Ebert MS, Zheng GX, Tsang JS, Sharp PA, van Oudenaarden A. MicroRNAs can generate thresholds in target gene expression. Nat Genet. 2011;43:854–9.
Wee LM, Flores-Jasso CF, Salomon WE, Zamore PD. Argonaute divides its RNA guide into domains with distinct functions and RNA-binding properties. Cell. 2012;151:1055–67.
Buchler N, Louis M. Molecular titration and ultrasensitivity in regulatory networks. J Mol Biol. 2008;384:1106–19.
Ala U, Karreth FA, Bosia C, Pagnani A, Taulli R, Léopold V, et al. Integrated transcriptional and competitive endogenous RNA networks are cross-regulated in permissive molecular environments. Proc Natl Acad Sci U S A. 2013;110:7154–9.
Bosia C, Pagnani A, Zecchina R. Modelling competing endogenous RNA networks. PLoS One. 2013;8:e66609.
Figliuzzi M, Marinari E, De Martino A. MicroRNAs as a selective channel of communication between competing RNAs: a steady-state theory. Biophys J. 2013;104:1203–13.
Figliuzzi M, De Martino A, Marinari E. RNA-based regulation: dynamics and response to perturbations of competing RNAs. Biophys J. 2014;107:1011–22.
Yuan Y, Liu B, Xie P, Zhang MQ, Li Y, Xie Z, et al. Model-guided quantitative analysis of microRNA-mediated regulation on competing endogenous RNAs using a synthetic gene circuit. Proc Natl Acad Sci U S A. 2015;112:3158–63.
Mullokandov G, Baccarini A, Ruzo A, Jayaprakash AD, Tung N, Israelow B, et al. High-throughput assessment of microRNA activity and function using microRNA sensor and decoy libraries. Nat Methods. 2012;9:840–6.
Nitzan M, Steiman-Shimony A, Altuvia Y, Biham O, Margalit H. Interactions between distant ceRNAs in regulatory networks. Biophys J. 2014;106:2254–66.
Bosson AD, Zamudio JR, Sharp PA. Endogenous miRNA and target concentrations determine susceptibility to potential ceRNA competition. Mol Cell. 2014; In press.
Kumar MS, Armenteros-Monterroso E, East P, Chakravorty P, Matthews N, Winslow MM, et al. HMGA2 functions as a competing endogenous RNA to promote lung cancer progression. Nature. 2014;505:212–7.
Denzler R, Agarwal V, Stefano J, Bartel DP, Stoffel M. Assessing the ceRNA hypothesis with quantitative measurements of miRNA and target abundance. Mol Cell. 2014;54:766–76.
Luna JM, Scheel TK, Danino T, Shaw KS, Mele A, Fak JJ, et al. Hepatitis C virus RNA functionally sequesters miR-122. Cell. 2015;160:1099–110.
Zhou X, Liu J, Wang W. Construction and investigation of breast-cancer-specific ceRNA network based on the mRNA and miRNA expression data. IET Syst Biol. 2014;8:96–103.
Huang CT, Oyang YJ, Huang HC, Juan HF. MicroRNA-mediated networks underlie immune response regulation in papillary thyroid carcinoma. Sci Rep. 2014;4:6495.
Paci P, Colombo T, Farina L. Computational analysis identifies a sponge interaction network between long non-coding RNAs and messenger RNAs in human breast cancer. BMC Syst Biol. 2014;8:83.
Xia T, Liao Q, Jiang X, Shao Y, Xiao B, Xi Y, et al. Long noncoding RNA associated -competing endogenous RNAs in gastric cancer. Sci Rep. 2014;4:6088.
Li L, Wang D, Xue M, Mi X, Liang Y, Wang P. 3′UTR shortening identifies high-risk cancers with targeted dysregulation of the ceRNA network. Sci Rep. 2014;4:5406.
Thomson DW, Bracken CP, Goodall GJ. Experimental strategies for microRNA target identification. Nucleic Acids Res. 2011;39:6845–53.
Li JH, Liu S, Zhou H, Qu LH, Yang JH. starBase v2.0: decoding miRNA-ceRNA, miRNA-ncRNA and protein-RNA interaction networks from large-scale CLIP-Seq data. Nucleic Acids Res. 2014;42(Database issue):D92–7.
Tollervey D. Molecular biology: RNA lost in translation. Nature. 2006;440:425–6.
Khachane AN, Harrison PM. Mining mammalian transcript data for functional long non-coding RNAs. PLoS One. 2010;5:e10316.
Kim HH, Kuwano Y, Srikantan S, Lee EK, Martindale JL, Gorospe M. HuR recruits let-7/RISC to repress c-Myc expression. Genes Dev. 2009;23:1743–8.
van Kouwenhove M, Kedde M, Agami R. MicroRNA regulation by RNA-binding proteins and its implications for cancer. Nat Rev Cancer. 2011;11:644–56.
Giza DE, Vasilescu C, Calin GA. MicroRNAs and ceRNAs: therapeutic implications of RNA networks. Expert Opin Biol Ther. 2014;14:1285–93.
Manikandan M, Munirajan AK. Single nucleotide polymorphisms in microRNA binding sites of oncogenes: implications in cancer and pharmacogenomics. OMICS. 2014;18:142–54.
Kawahara Y, Zinshteyn B, Sethupathy P, Iizasa H, Hatzigeorgiou AG, Nishikura K. Redirection of silencing targets by adenosine-to-inosine editing of miRNAs. Science. 2007;315:1137–40.
Cheng CJ, Bahal R, Babar IA, Pincus Z, Barrera F, Liu C, et al. MicroRNA silencing for cancer therapy targeted to the tumour microenvironment. Nature. 2015;518:107–10.
Rutnam ZJ, Yang BB. The non-coding 3ʹUTR of CD44 induces metastasis by regulating extracellular matrix functions. J Cell Sci. 2012;125:2075–85.
Yang J, Li T, Gao C, Lv X, Liu K, Song H, et al. FOXO1 3′UTR functions as a ceRNA in repressing the metastases of breast cancer cells via regulating miRNA activity. FEBS Lett. 2014;588:3218–324.
Fang L, Du WW, Yang X, Chen K, Ghanekar A, Levy G, et al. Versican 3′-untranslated region (3′-UTR) functions as a ceRNA in inducing the development of hepatocellular carcinoma by regulating miRNA activity. FASEB J. 2013;27:907–19.
Li C, Wang Y, Wang S, Wu B, Hao J, Fan H, et al. Hepatitis B virus mRNA-mediated miR-122 inhibition upregulates PTTG1-binding protein, which promotes hepatocellular carcinoma tumor growth and cell invasion. J Virol. 2013;87:2193–205.
Shen K, Mao R, Ma L, Li Y, Qiu Y, Cui D, et al. Post-transcriptional regulation of the tumor suppressor miR-139-5p and a network of miR-139-5p-mediated mRNA interactions in colorectal cancer. FEBS J. 2014;281:3609–24.
We apologize to all colleagues whose work could not be cited due to space constraints. We thank Tay lab members and Shen Mynn Tan for critical reading of the manuscript. Y.T. was supported by a Singapore National Research Foundation Fellowship and a National University of Singapore President’s Assistant Professorship.
Authors declare that they have no competing interests.
ASM and YT helped draft and edit the manuscript. Both authors read and approved the manuscript.
About this article
Cite this article
Sanchez-Mejias, A., Tay, Y. Competing endogenous RNA networks: tying the essential knots for cancer biology and therapeutics. J Hematol Oncol 8, 30 (2015). https://doi.org/10.1186/s13045-015-0129-1