miRNA-Based Regulation of Alternative RNA Splicing in Metazoans.

Alternative RNA splicing is an important regulatory process used by genes to increase their diversity. This process is mainly executed by specific classes of RNA binding proteins that act in a dosage-dependent manner to include or exclude selected exons in the final transcripts. While these processes are tightly regulated in cells and tissues, little is known on how the dosage of these factors is achieved and maintained. Several recent studies have suggested that alternative RNA splicing may be in part modulated by microRNAs (miRNAs), which are short, non-coding RNAs (~22 nt in length) that inhibit translation of specific mRNA transcripts. As evidenced in tissues and in diseases, such as cancer and neurological disorders, the dysregulation of miRNA pathways disrupts downstream alternative RNA splicing events by altering the dosage of splicing factors involved in RNA splicing. This attractive model suggests that miRNAs can not only influence the dosage of gene expression at the post-transcriptional level but also indirectly interfere in pre-mRNA splicing at the co-transcriptional level. The purpose of this review is to compile and analyze recent studies on miRNAs modulating alternative RNA splicing factors, and how these events contribute to transcript rearrangements in tissue development and disease.

1. Introduction

Multicellular organisms possess widely conserved post-transcriptional regulatory mechanisms used to develop and maintain cell and tissue identities. One of these mechanisms is alternative RNA splicing, which leads to the production of multiple protein isoforms from individual pre-mRNA molecules. Alternative RNA splicing is conserved in multicellular eukaryotes and occurs more frequently in complex organisms [1]. An estimated 95% of human multi-exonic genes are alternatively spliced [2].

Alternative RNA splicing is largely controlled by two well-characterized families of RNA binding proteins (RBPs): serine/arginine (SR)-rich proteins and heterogeneous nuclear ribonucleoproteins (hnRNPs). These two classes of proteins act in a dosage-dependent manner, competing for binding to specific sequences in pre-mRNAs to promote or inhibit exon inclusion, respectively [1]. SRs and hnRNPs are essential for constitutive as well as alternative RNA splicing. In most cases, SR proteins promote mRNA splicing activity and exon inclusion by binding to exonic (ESE) and intronic (ISE) splicing enhancer sequences on pre-mRNA molecules and recruit components of the spliceosome to the respective 5′ and 3′ splice sites [3,4,5]. However, in specific situations, SR proteins can also promote alternative exon skipping [2]. SR protein activity is restricted in regulating constitutive and alternative RNA splicing. Many of them have additional roles in modulating mRNA nuclear export, nonsense-mediated decay (NMD), and protein translation enhancement [6].

hnRNPs are instead typically described as inhibitors of splicing activity and exon inclusion. They bind to exon splicing silencers (ESSs) and intronic splicing silencers (ISSs) and prevent SR proteins and the spliceosome from binding to the pre-mRNA molecule [4,5]. However, similarly to SR proteins, the hnRNP family has highly diverse functions, which involve mRNA maturation, stabilizing mRNA during cellular transport, and activation or silencing of mRNA translation. hnRNPs can also promote or inhibit RNA splicing activity in specific circumstances [5,7,8].

While the SR and hnRNP proteins have been extensively studied and are fairly well characterized, numerous other RBPs are essential for the diversity seen in alternative RNA splicing. Examples of these RBPs include the RBFOX, Quaking, CUGBP and ETR-3-like factor (CELF), and RNA binding motif (RBM) protein families [9,10].

Despite the degree of specificity that can be achieved through several RNA splicing factors, other post-transcriptional regulatory mechanisms allow fine-tuning for tissue development and differentiation. Specifically, microRNAs (miRNAs), which are short non-coding (~22 nt) RNA molecules, exert a considerable influence by targeting motifs within the 3′ untranslated region (3′UTR) of genes, leading to inhibition of translation of target mRNAs [11] (Figure 1A). Although transcribing then degrading mRNA may seem counterproductive, it has been hypothesized that this mechanism evolved in metazoans to provide a way to “sculpt the transcriptome” to achieve optimal gene expression in various tissues after transcription occurs [11].

Figure 1

(A) miRNAs bind to target genes within their 3′UTR and repress their translation (Left). Since the pairing between the miRNAs and its target sequence is not perfectly complementary, a single miRNA may target and repress multiple mRNA transcripts (Right). (B) As shown in Tissue A and Tissue B, miRNAs may differentially regulate the abundance of RNA splicing factor proteins (shown as light green circles). Tissue 1 (Left) has a higher dosage of the RNA splicing factor due to a lack of miRNAs available to inhibit their translation. In this specific example, the abundance of splicing factors leads to an intron inclusion event (Isoform A). In contrast, in Tissue 2 (Right), miRNAs that target the RNA splicing factor lower its abundance, interfering with the splicing pattern of target genes. In this specific example, this leads to an intron exclusion event (Isoform B).

In metazoans, mature miRNAs are incorporated into the microRNA-induced silencing complex (miRISC), which recruits deadenylase complexes to the target 3′UTRs [11]. This pairing has been shown to require as little as six consecutive nucleotides in the 5′ end, or seed region, of the mature miRNA. Perfect complementarity within the seed region is considered the canonical indicator of miRNA targeting. However, many recent studies indicate that miRNAs can target non-canonical elements in target mRNAs [12,13]. Currently, bioinformatic approaches are the only tools for high-throughput miRNA target detection. Based on these algorithms, it is proposed that each miRNA controls hundreds of gene products [14] (Figure 1B).

The human genome contains ~2000 distinct mature miRNAs located in intergenic regions, introns, and operon-like clusters. Metazoan miRNAs were first characterized in C. elegans as heterochronic regulators that temporally repress genes involved in developmental transitions [15,16]. miRNAs and their regulatory mechanisms are conserved from humans [17] to Cnidarians [18]. In the past years, several groups produced tissue-specific localization data for many miRNAs in worm, mouse, rat, and human somatic tissues [19,20,21] and in cancers [22]. These results unequivocally show distinct functional miRNA populations in tissues, which can reshape transcriptomes and contribute to cell identity acquisition and maintenance.

In recent years, several intriguing studies connected miRNA targeting to alternative RNA splicing. This novel form of regulation is achieved indirectly through miRNA targeting and lowering the dosage of several RNA splicing factors, which in turn change the splicing pattern of many genes in a tissue-specific manner. The degree of regulation that miRNA pathways exert in alternative RNA splicing and tissue differentiation has yet to be fully understood at an organismal level, but it is exciting. In this view, miRNAs do not only regulate gene expression at the post-transcriptional level, as canonically described, but also indirectly at the co-transcriptional level, throughout the modulation of these RNA splicing factors. In C. elegans, for example, virtually all SR and hnRNP orthologs are targeted by miRNAs, and the depletion of the miRISC leads to widespread changes in splice junction usage [19]. These broad splicing errors have also been observed in transcriptomes from Dicer knockout mice [23], suggesting that this form of gene regulation is not only present in vertebrates but may have evolved earlier in the tree of life.

In humans, it is very common for RNA splicing factors to contain predicted or validated miRNA targets (TargetScan) in their 3′UTRs. Still, the contribution of miRNAs to protein isoform production has not been fully characterized yet. This review highlights and discusses this new form of miRNA-based regulation and its potential implications in development and disease. Table 1 summarizes the RNA splicing factors with their respective miRNAs discussed below.

Table 1

Summary of the RNA splicing factors discussed in this review.

		miRNA(s)	Spliced Genes Affected	Development/Disease	Accession #(s)
SR Proteins	SRSF1	miR-10-a/bmiR-28miR-505miR-766-3p	Tau, Ich-1, Bcl-x	Various cancers	[24,25,26,27,28]
	SRSF2	miR-183-5pmiR-193a-3pmiR-200c-3p	-	Various cancers	[28,29,30,31,32,33]
	SRSF6	miR-146amiR-193a-5pmiR-506-3p	OGDHL, ECM1	Pancreatic cancer, Alzheimer’s disease, pleural fibrosis	[34,35,36,37]
	SRSF7	miR-30a-5pmiR-181a-5pmiR-188miR-216b-5p	SPP1	Renal cancer, kidney injuries	[38,39]
	SRSF9	miR-1/206miR-802	⎯	Bladder and cervical cancer	[40,41,42]
hnRNPs	hnRNP A1	miR-1-3pmiR-135a-5pmiR-137miR-149-5pmiR-206miR-424miR-451miR-503	-	Acute myeloid leukemia, colon cancer	[28,43]
	hnRNP A2	miR-124miR-340	PKM	-	[28]
PTBs	PTBP1	miR-1/206miR-124miR-133bmiR-137miR-194-5pmiR-340	PKM, APP	Various cancers	[44,45,46,47]
	PTBP2	miR-132miR-133	Tau, PKM	Glioblastoma	[48,49,50]
RBFOX	RBFOX1	miR-129-5pmiR-980	-	-	[51,52,53]
	RBFOX2	let-7gmiR-9miR-135a	Scn5a, Kcnd3	Myotonic dystrophy type 1 (DM1)	[54]
	RBFOX3	miR-129-5p	-	-	[53]
Quaking	QKI-5	miR-143-3pmiR-200cmiR-214miR-221miR-375	-	Various cancers	[55,56,57,58]
	QKI-6	miR-29amiR-214	-	Glioblastoma	[59,60]
	QKI-7	miR-214	-	-	[59]
CELF	CELF1	miR-23a/bmiR-322miR-503	-	-	[23,61]
	CELF2	miR-20amiR-23a/bmiR-95-3p	-	Glioblastoma	[23,62,63]
NOVA	NOVA1	miR-27a-3pmiR-181b-5pmiR-203a-3pmiR-339miR-592miR-146b-5p	-	Gastric and thyroid cancers	[64,65,66,67]
	NOVA2	miR-7-5p	-	Non-small cell lung cancer	[68]
RBM	RBM10	miR-133a/bmiR-335	Numb	Endometrial cancer	[69]
	RBM24	miR-222	Coro6, Fxr1, NACA	Myogenic differentiation	[70]
	ESRP1	miR-15a/bmiR-16miR-23a/b	CD44, FGFR2, EPB41L5, Rac1	Various cancers	[71,72]

2. RNA Splicing Factors Regulated by miRNAs

2.1. SR Proteins

SR proteins are regulators of both constitutive and alternative RNA splicing, and they are essential for proper tissue development and differentiation. The human genome contains at least 14 SR protein genes [73]. Dysregulation of many SR proteins has been characterized in several types of diseases, especially cancer and developmental disorders, where widespread aberrant RNA splicing events have been observed. As detected in several types of cancers, individual SR proteins may act as oncogenes or tumor suppressor genes, which is dependent on their tissue-specific expression [74,75]. Many SR proteins have miRNA targets in their 3′UTRs and have shown to be targeted and repressed by these miRNAs, leading to widespread RNA splicing errors. Several examples of miRNA-based SR protein regulation have been described for several specific SR proteins and are discussed below.

2.1.1. SRSF1 (SF2/ASF)

SRSF1 is a well-characterized general RNA splicing factor involved in both constitutive and alternative RNA splicing and influences the location of the splice site in a dosage-dependent manner. The RNA splicing function of SRSF1 is different depending on the cell subtype. When overexpressed in HeLa cells, SRSF1 promotes exon skipping in Ich-1, a member of the caspase family of proteases, to promote the longer pro-apoptotic Ich-1L isoform, which uses a more distal stop codon. In contrast, the shorter anti-apoptotic Ich-1s isoform retains the 61 bp exon (Figure 2). However, when SRSF1 is overexpressed in HEK 293 cells, the longer anti-apoptotic isoform of Bcl-xL, is favored over the shorter pro-apoptotic Bcl-x(s) isoform (Figure 2). The 3′UTR of SRSF1 possesses several miRNA targets, and coincidentally, the miRNAs that target SRSF1 are downregulated in many cancers. Verduci et al. demonstrated that the overexpression of miR-28 and miR-505 in mouse embryonic fibroblasts repressed SRSF1, which allowed for the ratios of Ich-1 and Bcl-x isoforms to become more balanced [24]. A previous study also supports these findings [76]. Taken together, these results suggest that miRNA-based regulation is crucial for the precise dosage of RNA splicing factors to produce the correct ratios of alternatively spliced isoforms [24] (Figure 2). However, it is worth noting that overexpressing miR-505 alone repressed SRSF1 expression in mouse embryonic fibroblasts but unexpectedly promoted the anti-apoptotic Ich-1s isoform [24], suggesting that other mechanisms are responsible for the proper expression of pro-apoptotic spliced gene isoforms in this mouse model [24].

Figure 2

In HeLa cells, SRSF1 promotes exon skipping in Ich-1 to produce the pro-apoptotic Ich-1L isoform. However, in HEK 293 cells, SRSF1 promotes the expression of Bcl-xL, the anti-apoptotic isoform of Bcl-x, by favoring a distal 5′ splice site. miR-28 and miR-505 directly target the 3′UTR of SRSF1, and the result induces a shift of the isoform ratios of Ich-1 and Bcl-x.

The relationship between miRNA regulation on SRSF1 expression and alternative RNA splicing has also been studied in neuroblastoma cells. When these cells are treated with retinoic acid, proliferation is inhibited, and apoptosis is promoted [77]. Further studies identified increased expression of miR-10a and miR-10b in these cells, which both target the 3′UTR of SRSF1. As predicted, SRSF1 levels decreased in neuroblastoma cells after retinoic acid treatment, which presumably altered the alternative RNA splicing of other mRNA transcripts. This was demonstrated by a switch to exon 10 inclusion in the tau gene [25], a key player in neurological disorders.

Interestingly, SRSF1 may also act as an oncogene in several cancers. When overexpressed, SRSF1 acts as a proto-oncogene by inducing alternative RNA splicing in multiple genes that operate in the apoptotic pathway, promoting cancer establishment or progression [78,79]. Predictably, an increase in SRSF1 expression in cancers is correlated with the appearance of new aberrant RNA isoforms [26] and poor prognosis [78,80].

SRSF1 is also targeted by miR-766-3p, a well-characterized tumor suppressor or inducer in several types of cancers [27,81,82,83], which is often markedly downregulated in cancers, leading to the overexpression of SRSF1 [27,82]. Unfortunately, while very important, the exact RNA splicing errors induced by the overexpression of SRSF1 in the absence of this miRNA are unclear and difficult to pinpoint due to the various roles beyond alternative RNA splicing in which SRSF1 is involved.

2.1.2. SRSF2

SRSF1 is not the only member of the SR protein family known to be targeted by miRNAs. SRSF2, a second member of the SR family of proteins, similarly to SRSF1, is also widely expressed and has a diverse role in many biological processes [84]. This gene is frequently upregulated in many human cancers, and this event is associated with poor prognosis in patients. Similar to SRSF1, SRSF2 overexpression leads to global alternative RNA splicing errors [29,30]. SRSF2 promotes apoptosis through alternative RNA splicing of specific cancer genes, such as c-flip, caspases-8 and -9, and Bcl-x, to their pro-apoptotic isoforms, and a decrease in its expression is associated with the inhibition of apoptosis [31,32]. Because of the significant role SRSF2 plays in alternative RNA splicing in apoptotic pathways, the dysregulation of miRNAs that target SRSF2 is a fingerprint in many types of cancers. In renal cell cancer, for example, miR-183-5p and miR-200c-3p are significantly upregulated and directly target SRSF2 to prevent the activation of apoptotic pathways [28,32,83]. Another miRNA, miR-193a-3p, which is upregulated in multiple other cancers, has been shown to directly target SRSF2 in nasopharyngeal cancer [33]. Interestingly, SRSF2 is not the only SR protein gene targeted by miR-193a-3p, as this miRNA has been shown to also directly target SRSF6 [34]. Unfortunately, the global RNA splicing defects caused by the change in SRSF2 expression levels by miRNAs have not yet been characterized in detail.

2.1.3. SRSF6

SRSF6 is another SR protein member also known to be targeted by miRNAs. SRSF6 is unfortunately less well characterized when compared to other SR proteins, but it has been shown to be dysregulated in a variety of human diseases. SRSF6 overexpression has been observed in lung and colon cancers [85] and Huntington’s disease [35], but surprisingly, it is downregulated in pancreatic cancer [34].

Among other genes, SRSF6 modulates alternative RNA splicing of the genes OGDHL (oxoglutarate dehydrogenase-like protein) and ECM1 (extracellular matrix protein 1), both known to play roles in cancer metastasis. Upregulation of miR-193a-5p (also a regulator of SRSF2 described previously) directly targets and inhibits SRSF6, and as a result, the splice variants of OGDHL and ECM1 conducive to tumorigenicity are expressed in pancreatic cancer [34].

In a mouse model studying miRNA-based treatments for Alzheimer’s disease (AD), increased activity of miR-146a caused the depletion of SRSF6 expression [36]. The exact mechanisms involving miRNA regulation and SRSF6 in neurological disorders have yet to be fully understood, but unbalanced SRSF6 expression alters the splicing pattern of several neuronal target mRNAs, including the tau gene [86], a well-known player in AD also expressed with aberrant RNA splicing in other neurological diseases.

Furthermore, SRSF6 is a direct target of miR-506-3p, and its overexpression caused by the repression of miR-506-3p was observed in pleural fibrosis, a thickening of the membranes that cover the lungs [37]. Unfortunately, the downstream effects have yet to be directly identified. As discussed, SRSF6 has also been shown to be up- or downregulated in a variety of diseases, and this regulation was often found to be dependent on SRSF6-specific miRNAs.

2.1.4. SRSF7

SRSF7 is an essential splicing regulator, and its aberrant expression has also been implicated in several types of cancers. SRSF7 is known to direct alternative RNA splicing of several genes, including BRCA1 [87] and SPP1 [38], which are well-known regulators of tumorigenesis.

The SRSF7 3′UTR contains at least seven miRNA targets (TargetScan), an indicator of the need to precisely dose the abundance of this RNA splicing factor in cells. Boguslawska et al. demonstrated specifically in renal cancer that miR-30a-5p, miR-181a-5p, and miR-216b directly target the 3′UTR of SRSF7 and repress its expression, and in turn, affect the alternative RNA splicing of SPP1. Interestingly, miR-216b repression of SRSF7 led to change in the alternative RNA splicing pattern of SPP1, favoring the SPP1-c isoform, while decreasing the SPP1-b isoform, which is instead boosted in the presence of SRSF7 [38].

These three miRNAs are not the only regulators of SRSF7. MiR-188 was also found to target the SRSF7 3′UTR in an acute kidney injury model. This interaction led to a decrease in cell viability, potentially induced by the alteration of RNA splicing patterns of anti-apoptotic gene isoforms caused by the lack of SRSF7 [39].

2.1.5. SRSF9

SRSF9 is an important myogenic regulator. Its depletion, together with other RNA splicing factors, is an essential step in the myogenesis of muscle development. Specifically, miR-1 and miR-206 directly target the SRSF9 3′UTR to allow for the transition from proliferation to differentiation in myogenesis [40]. In developed tissues, SRSF9 acts as an oncogene as it represses apoptosis through alternative RNA splicing of genes, such as Bcl-x [88]. Continued repression of SRSF9 is dependent on miR-1, which is downregulated in several types of cancers, including bladder cancer [41]. In cervical cancer cells, miR-802, which targets the 3′UTR of SRSF9, was found to be downregulated while SRSF9 was upregulated. Upon transfection of miR-802 in cervical cancer cells, SRSF9 expression was inhibited, leading to inhibited cell proliferation, cell cycle arrest, and apoptosis [42]. While a few players have been identified [88], the widespread changes in RNA splicing patterns caused by miRNAs’ deregulation of SRSF9 function in various cancers have yet to be identified.

2.2. hnRNP Proteins

Working as antagonists to SR proteins, hnRNP proteins typically promote exon exclusion in a dosage-dependent manner. Dysregulation of hnRNPs is described in several cancers and is also characterized in other diseases, such as amyotrophic lateral sclerosis and neurodegenerative disorders [89]. Additionally, the hnRNP family is comprised of about 20 proteins and is the most abundantly expressed group of RNA binding proteins in mammalian cells [89]. Similar to SR proteins, changes in their abundances lead to widespread RNA splicing errors. A few examples of these errors induced by miRNA deregulation are discussed below.

2.2.1. hnRNP A1

hnRNP A1 is one of the main components in the spliceosome, and is essential for constitutive RNA splicing [90] but is also known to modulate alternative RNA splicing in a dosage-dependent manner [89]. hnRNP A1 is ubiquitously expressed, and dysregulation of different miRNAs that target this gene has been observed in diseases or disorders in distinct tissues. For example, a study by Sokół et al. (2018) found hnRNP A1 to be targeted by miR-135a, miR-1-3p, miR-206, and miR-149-5p in renal cancer cells and other types of cancer [28]; however, in acute myeloid leukemia, hnRNP A1 is overexpressed due to downregulation of miR-451, which directly targets hnRNP A1 [43].

A recent study showed that in the presence of miR-424 and miR-503, these two miRNAs were able to directly repress hnRNP A1 and potentially alter the RNA splicing pattern of hnRNP A1-regulated genes, but more extensive studies beyond bioinformatics approaches have yet to support this hypothesis [91].

These various examples involving hnRNP A1 suggest that the abundance of this gene is tightly controlled by miRNA networks, and the deregulation of these networks may lead to dramatic rearrangements in RNA splicing patterns in many transcripts controlled by this RNA splicing regulator, leading to disease states.

2.2.2. hnRNP A2

The RNA splicing factor hnRNP A2 is closely related to hnRNP A1 and is hypothesized to have originated from a duplication event [92]. As previously mentioned, several hnRNP members, including both hnRNP A1 and hnRNP A2, alternatively splice the PKM gene to exclude exon 9 and include exon 10 to produce the PKM2 isoform [93,94,95]. A recent study using exon-tiling microarrays showed that knockdown of hnRNP A2 with RNAi promoted at least two exon inclusion and four exon skipping events in the genes MTA3, MAP9, EPB41L4A, TP53INP2, RXFP1, and PHF14 [96], and potentially more, suggesting that hnRNP A2 is also able to induce alternative RNA splicing in a large number of genes.

Similar to hnRNP A1, the hnRNP A2 3′UTR contains several predicted and validated miRNA targets, suggesting that the dosage of this gene is also tightly regulated in cells. hnRNP A2 is directly targeted by at least three miRNAs. A 2012 study by Sun et al. found that hnRNP A2 is repressed by miR-124 and miR-340, while hnRNP A1 is repressed by miR-137 in a luciferase assay [94], which suggests that despite their similarities, hnRNP A1 and hnRNP A2 are regulated by different miRNAs.

Although miRNAs are typically repressors of gene expression, the expression of hnRNP A2, as well as the SFPQ protein (another RNA splicing factor), increased when miR-369 was transfected in mouse adipose-derived mesenchymal stem cells [97]. In this study, hnRNP A1 only showed a slight increase when compared to hnRNP A2 and SFPQ. While an exact mechanism was not established, it is important to note that the 3′UTR of hnRNP A2 but not hnRNP A1 contained a predicted miR-369 binding site, which seems to be able to stabilize hnRNP A2 mRNAs [97,98].

2.2.3. Polypyrimidine Tract Binding Proteins

Polypyrimidine tract binding (PTB) proteins are named because of their ability to bind intronic polypyrimidine tracts during the pre-mRNA splicing event influencing pre-mRNA processing [99]. Although their primary role is to modulate pre-mRNA splicing, they have been shown to also play a role in mRNA metabolism and mRNA transport. The PTB proteins have well-characterized roles in the biogenesis and processing of specific miRNAs, which often form regulatory feedback loops and allow for the transition of specialized cells [100,101,102]. PTBP1 and PTBP2 are the two most well-characterized members of the PTB family. PTBP1 is a major repressive regulator of alternative RNA splicing, causing exon skipping in many alternatively spliced precursor mRNAs [5,103]. PTBP1 is highly expressed in most of the tissue types but is repressed in neuronal tissues, allowing for the inclusion of neuron-specific exons in many genes [44]. Both PTBP1 and PTBP2 have several predicted (TargetScan) and validated miRNA binding sites in their 3′UTRs. MiRNAs have been shown to exert considerable influence in neuron-specific alternative RNA splicing by targeting PTBP1 transcripts, which in turn upregulate neuron-specific PTBP2 expression. Coincidentally, several of these miRNAs, including miR-124 and miR-137 [45,104], are highly abundant in neuronal tissues, which supports the idea that an miRNA-based form of regulation is required for the correct tissue-specific expression of PTB isoforms.

Dysregulation of PTBP1-specific miRNAs is implicated in neurological diseases, including AD and brain cancers. Abnormal RNA splicing of the amyloid precursor protein (APP) in the neurons of AD patients was observed with downregulation of miR-124, a direct target of PTBP1, which instead is upregulated [46].

Several non-neuronal-specific miRNAs are also regulators of PTBP1 in various cancers. PTBP1 promotes the expression of the pyruvate kinase M1/2 isoform (PKM2) through alternative RNA splicing, which is expressed in proliferating cells [105] and is crucial for establishing the “Warburg effect” in many types of cancers [106]. Overexpression of PTBP1 in colorectal tumors was attributed to the downregulation of muscle-specific miR-1 and miR-133b, which target the PTBP1 3′UTR [47]. Additionally, muscle-specific miR-206 directly targets the PTBP1 3′UTR [45] and its dysregulation is implicated in sarcomas [107].

PTBP2, also known as nPTB, is the neuron-specific PTB isoform downregulated in differentiated muscle tissue. Similar to PTBP1, the muscle-specific miR-133 targets PTBP2 transcripts during myoblast differentiation. Changes in miR-133 expression and its regulation of both PTB isoforms resulted in mRNA splicing errors in several genes including those involved in muscle cell differentiation [48]. Dysregulation of miR-132, a brain-specific miRNA that targets the PTBP2 3′UTR, has been implicated in several disorders. In glioblastoma cells, miR-132 promotes proliferation and self-renewal potential by inhibiting PTBP2, which is implicated in cell differentiation [49]. Proper expression of PTBP2 in neuronal cells is crucial for alternative RNA splicing of neuronal genes. For example, in the brain, the tau proteins are equally expressed as two isoforms: 4R-tau and 3R-tau, which either include or exclude exon 10 from the primary tau mRNA. PTBP2 was shown to promote the 4R-tau isoform. The upregulation of miR-132, which targets the 3′UTR of PTBP2, causes a decrease in PTBP2 expression, leading to a tau isoform imbalance, which is observed in progressive supranuclear palsy [50].

2.3. Other RNA Splicing Factors

Several other RNA splicing factor protein families are also essential for alternative RNA splicing and are modulated through miRNA regulation. We will now focus on five of the most well-characterized families: RBFOX, CELF, NOVA, Quaking, and RBM protein families.

2.3.1. RBFOX

The RBFOX protein family is highly conserved across metazoans and consists of three paralogs in mammals: RFBOX1, RBFOX2, and RBFOX3, which are expressed in a tissue-specific pattern [108]. All RBFOX paralogs have been found to synergize or antagonize the activities of other RNA splicing factors, which would allow for a higher level of tissue-specific alternative RNA splicing [108].

The RNA splicing factor RBFOX1 is implicated in mammalian neuronal development, myoblast fusion, and skeletal muscle and heart function [108]. The C. elegans RBFOX-1 homolog (fox-1), which has been shown to control male or hermaphrodite development by splicing several genes including xol-1 pre-mRNA [109], is expressed in body muscle tissue and is targeted by miRNAs in intestinal tissue [19]. While the mammalian homologs of RBFOX1 may have distinct or unrelated functions, RBFOX1 is also expressed in muscle tissue and neurons [110], and is needed for proper neuronal development [111], and its absence or alteration has been implicated in a variety of neurological processes and disorders [112]. The Rbfox1 3′UTR contains several predicted (PicTar) and validated miRNA targets. A 2016 study identified Rbfox1 in Drosophila melanogaster as a memory-promoting gene and conversely, one of its miRNA partners, miR-980, as a memory suppressor [51]. However, later research clarified that only Rbfox1 transcripts with extended 3′UTRs were targeted by miR-980 [52], leaving shorter transcripts intact. The ability of genes to express different 3′UTR lengths is possible through the mechanism of alternative polyadenylation (APA). APA has been shown by several groups as a potent mechanism to bypass miRNA regulation [19,113,114,115]. In addition, another miRNA, miR-129-5p, has been shown to directly inhibit the expression of Rbfox1 and Rbfox3 transcripts through targeting their 3′UTRs, which both have multiple miR-129-5p binding sites [53].

Another member of the RBFOX family, Rbfox2, is directly targeted by let-7g, miR-9, and miR-135a, and is expressed with two tissue-specific isoforms, Rbfox2 and Rbfox2, which differ by an alternatively spliced exon (either 40 or 43 nt long). Although these two isoforms seem to have redundant traits, they are functionally diverse, and their dysregulation results in severe complications [54]. The non-muscle-specific isoform, RBFOX240, is overexpressed in the hearts in myotonic dystrophy type 1 (DM1) patients. A combination of repressed expression of let-7g, miR-9, and miR-135a and elevated expression of the RNA splicing factor CELF1, discussed below, were identified to drive Rbfox2 expression in a mouse heart model of DM1 by Misra et al. 2020. RNA mis-splicing of Scn5a and Kcnd3 ion channel transcripts and potentially other genes are consequential of abundant Rbfox2 expression [54], leading to ion channels that function more slowly and are more prone to arrhythmias [116]. Taken together, all these studies indicate the RBFOX family is strongly regulated by miRNAs, which suggest that alternative RNA splicing is also subject to this type of regulation.

2.3.2. CELF Protein Family

CELF (CUGBP and ETR-3-like factor) family protein homologs are widely conserved and are involved in post-transcriptional mechanisms, including RNA editing and translation, in addition to alternative RNA splicing. CELF1 and CELF2 regulate several developmentally related RNA splicing events in the heart and are highly expressed in the brain [23,117]. As previously mentioned, overexpression of CELF1 promotes the expression of the Rbfox2 isoform observed in DM1 patients [54]. The regulation of CELF1 and CELF2 through miRNAs is essential for proper heart development. Kalsotra et al. (2008) found that the change in CELF1 and CELF2 protein levels decreased by 10- and 18-fold, respectively, from embryonic to adult hearts from mice with no change in mRNA levels. To investigate if these changes in protein levels were regulated by miRNA, Kalsotra et al. (2010) performed a knockout of Dicer in mouse heart tissue, which significantly increased CELF protein expression. Both of the CELF1 and CELF2 3′UTRs contain predicted targets for miRNA-23a/b, which were confirmed by luciferase assays [23]. Furthermore, administration of an miR-23a antagomir in the adult mouse heart sequestered miR-23a/b, which caused a significant increase of CELF1 and CELF2 protein levels and shifted the alternative RNA splicing of pre-mRNAs to embryonic-specific isoforms. In addition to miR-23a/b, luciferase assays indicated that miR-322/-503 also targeted CELF1 as shown by the luciferase activity decreasing by more than 50% in the presence of these miRNAs with the CELF1 3′UTR [61], which suggests CELF1 and CELF2 have distinct patterns of expression or function.

Dysregulation of CELF2 has been described in several cancers, particularly glioblastomas. Profiling of miRNA populations from human glioma tissues and glioblastoma cell lines identified the upregulation of several miRNAs that are predicted to target CELF2. Subsequent research utilizing luciferase assays with miR-95-3p and miR-20a provided substantial evidence of their role in CELF2 inhibition, promotion of tumor development, and progression in glioblastomas [61,62,63]. The distinct expression patterns of CELF1 and CELF2 could possibly be reliant on proper miRNA regulation.

2.3.3. NOVA

The NOVA (neuro-oncological ventral antigen) protein family members, NOVA1 and NOVA2, are both RNA splicing factors that were first identified for their roles in neurological diseases. The roles of NOVA1 in cancer are not yet fully understood, but correlations between Nova1-targeting miRNAs and NOVA1 expression have been described in multiple types of cancers. Nova1 is directly targeted by miR-181-5p and miR-203a-3p in neurons [64,65], and is also expressed in various other tissues and is often implicated in cancers. MiR-592 is considerably underexpressed in thyroid cancer while its target, Nova1, is overexpressed [66]. NOVA1 is also overexpressed in gastric cancer with the downregulation of miR-339, another miRNA that directly targets the Nova1 3′UTR [67]. In addition, miR-146b-5p directly targets Nova1, and its upregulation was observed in tissues surrounding the sites of tumors in gastric cancer after gastrectomy [118]. Contradictory to other studies, repression of NOVA1 by the overexpression of miR-27a-3p, which directly targets the Nova1 3′UTR, has been identified in gastric cancer [119,120,121]. A study examining the role of miR-7-5p in non-small cell lung cancer (NSCLC) identified Nova2 as a direct target [68]. In NSCLC cell lines, decreased expression of miR-7-5p and subsequent overexpression of NOVA2 was identified as the cause of tumor cell proliferation, migration, and invasion [68].

2.3.4. Quaking Protein Family

The Quaking (QKI) protein is a member of the STAR (signal transduction and activation of RNA) family of RBPs and was originally named due to the persistent tremors seen in QKI mutant mice [122]. In humans, each of the three QKI isoforms has distinct roles in cell differentiation [55,56,123] despite differing in their C-termini and 3′UTRs [57]. QKI-5 is expressed in the nucleus during embryogenesis and regulates pre-mRNAs for myelination proteins by alternative RNA splicing, and then is downregulated after birth [58]. In contrast, QKI-6 and QKI-7 are expressed later in embryogenesis in the cytoplasm [56], and their mRNAs possess different 3′UTRs than QKI-5 [124]. Several studies indicate that QKI proteins are potent mediators of differentiation in several cell types, for example, neurogenesis in neuronal progenitor cells (NPCs) [124] and endothelial cells that enter endothelial-to-mesenchymal transition (EMT) [125]. EMT is a process where epithelial cells lose cell polarity and adhesion and gain invasive properties typical of mesenchymal cells. EMT is a common hallmark of malignant cancers.

Coincidentally, several miRNAs that specifically target QKI proteins are upregulated during cell proliferation and development. The timing of neuronal differentiation depends on appropriate expression and regulation of miRNAs, such as miR-214-3p, which is upregulated in NPCs during neurogenesis and dendritic development, inhibiting the expression of all three QKI isoforms [59,124]. Similar patterns of QKI inhibition by miR-214-3p have also been identified in epithelial cells during angiogenesis [125]. Additionally, analysis of the 2.3 kb 3′UTR of Qki-5 has shown multiple predicted binding sites for miR-200c and miR-375; consequently, overexpressing these specific miRNAs significantly decreased QKI-5 expression in epithelial cells, which resulted in a decrease of their target mesenchymal-specific RNA alternatively spliced genes [126]. Additionally, upregulation of miRNA-200c was associated with multiple types of cancers, presumably by decreasing QKI-5 expression, resulting in dysregulation of alternative RNA splicing events [126]. Similarly, in colorectal cancer, Qki-5 was identified as a target for miR-221, which is a necessary factor for in vivo tumor growth [127]. Furthermore, QKI-5 expression is associated with tumorigenic properties in esophageal squamous cell carcinoma, such as proliferation, tumor invasion, and mitigation. These tumorigenic traits are inhibited, and apoptosis is induced by miR-143-3p, which directly targets Qki-5 [128]. Similarly, QKI-6 was found to be upregulated in glioblastomas and its targeting miRNA, miR-29a, to be downregulated. The underlying mechanism could be attributed to QKI-6 regulation of WTAP, a protein that stimulates epidermal growth factor (EGF) signaling [60], but more experiments need to be performed to validate this possibility.

2.3.5. RNA Binding Motif (RBM) Proteins

RBM proteins are loosely defined by the presence of at least one RNA-recognition motif (RRM), which they use to target specific mRNAs [129]. Several RBM proteins are involved in cardiomyocyte differentiation, myofibrillogenesis, and muscle-specific alternative RNA splicing.

One of its members, RBM24, is instrumental in cardiac and skeletal muscle development, by promoting exon inclusion to produce muscle-specific protein isoforms of several genes. A recent study by Cardinali et al. (2016) found that the Rbm24 3′UTR is directly targeted by miR-222 during the differentiation of skeletal muscle cells. Consequently, the overexpression of miR-222 led to the repression of Rbm24, resulting in defective exon inclusion of Coro6, Fxr1, and NACA muscle-specific gene isoforms [70]. This research supports the role of dysregulated miR-222 in several muscular disorders [21].

In addition, another RBM protein, RBM10, is involved in alternative RNA splicing of the conserved gene Numb, which possesses both a short and a long gene isoform, Numb-S and Numb-L, respectively. Expression of Numb-L is associated with cell proliferation, and its upregulation leads to various types of cancers [69]. Rbm10 is targeted by multiple miRNAs, including miR-133a, miR-133b, and miR-335 [69]. Interestingly, in tumors, only miR-335 is increased and inhibits RBM10 expression, leading to the expression of the oncogenic Numb-L isoform, specifically in endometrial cancer [69], and perhaps interfering with the splicing pattern of other genes.

Endothelial splicing regulatory proteins 1 and 2 (ESRP1 and ESRP2) are two members of the RBM family that regulate EMT, an event critical for proper tissue differentiation, through their downregulation [130,131]. Dysregulation of EMT events through the overexpression of ESRP1 in particular, has dire consequences for cell differentiation and consequently is typically observed in many types of cancers [132,133]. An exception to this pattern has been observed in pancreatic cancer, where ESRP1 acts as a tumor suppressor while miR-23a, which directly targets the 3′UTR of ESRP1, promotes invasion and metastasis of cancer cells [71].

Considering these examples, miRNA-based regulation has clear implications involving RBM proteins and their alternative RNA splicing activities.

2.4. Indirect Modulators of miRNA-Based RNA Alternative Splicing Events

This review has mainly focused on the canonical mechanism of miRNA-based regulation, where a given miRNA directly targets a 3′UTR of an RNA splicing factor, dosing its expression, and leading to alternative RNA splicing events to target genes. Recently, several exciting new studies highlighted a few additional mechanisms, where the change in dosage of a given RNA splicing factor, is either due to the depletion of its target miRNAs through RNA sponges, or the depletion of transcription factors involved in transcribing the RNA splicing factors by miRNAs.

2.4.1. circRNAs

Circular RNAs (circRNAs) are a recently discovered class of long ncRNAs in which the 5′ and 3′ termini are covalently linked, forming a circular structure. CircRNA are broadly expressed in mammalian cells and are involved in many cellular processes [134]. CircRNAs have been recently described as robust modulators of miRNA regulation of alternative RNA splicing factors, by acting as miRNA sponges. As expected, circRNA misregulation is correlated with specific forms of cancer, for example, ovarian cancer. A 2020 study identified that pervasive overexpression of ESRP1 in ovarian cancer is dependent on the expression of circ-0005585, which acts as an miRNA sponge for the ESRP1-specific miRNAs miR-15a/15b/16 and miR-23a/b, by expressing mutual miRNA targets as the 3′UTR of ESRP1 [72]. The level of competition between circ-0005585 and the 3′UTR of ESRP1 for mutual miRNAs was significant enough to allow ESRP1 overexpression, which suggest that circRNAs have a considerable influence on development and disease by modulating miRNA regulation.

In addition, it is important to note that circRNAs may modulate miRNA activity upstream of specific alternative RNA splicing factors, which nonetheless may significantly alter the downstream expression of a myriad of genes. This potent form of miRNA regulation is demonstrated by a 2018 study, which identified circANKS1B, a circRNA, to be distinctly correlated with breast cancer metastasis. Specifically, circANKS1B upregulated the expression of an ESRP1-specific transcription factor, USF1, by competitively binding the USF1 3′UTR-specific miRNAs. Thus, the increased expression of USF1 allowed for the increased expression of ESRP1 [135].

2.4.2. Transcription Factors

CircRNAs are not the only modulators of miRNA-induced RNA alternative splicing events. Recently, an interesting study linked miRNA-controlled transcription factors with RNA alternative splicing events. The sequencing of RNAs in human airway epithelial cells undergoing EMT showed that miR-133a is significantly overexpressed in these cells [136], forcing these cells into spontaneous EMT, through the downregulation of the grainyhead-like 2 (GRHL2) gene, an epithelial transcription factor involved in several disorders, such as the ectodermal dysplasia/short stature syndrome [137] and deafness autosomal dominant 28 [138]. The loss of GRHL2 expression led to the downregulation of the ESRP1 [136]. One of its characterized targets is p120-catenin. Loss of ESRP1 lead to an isoform switch that caused the loss of E-cadherin, which is one of the hallmarks of mesenchymal cells [136]. Although this is only one example, it is relatively common for transcription factors to possess many predicted and validated miRNA targets in their 3′UTRs, suggesting that perhaps this novel regulatory mechanism is more widespread.

3. Conclusions

It is clear that miRNAs are important modulators of alternative RNA splicing networks. Disruption of specific miRNA targeting pathways leads to RNA splicing errors that have been described in many types of cancers, neurological disorders, and developmental dysregulation. Unfortunately, apart from some of the examples discussed in this review, relatively little is known about splicing changes associated with micro-RNA-mediated repression. Still, the field is growing, with more interactions identified in recent years. In this review, we covered several examples of alternative RNA splicing induced by the dysregulation of RNA splicing factors. It is very challenging to pinpoint the direct effects that miRNAs have on alternative RNA splicing, but the widespread consequences of these changes contribute significantly to disease and developmental disorders.