The effects of alternative splicing on miRNA binding sites in bladder cancer

MiR-6721-5p as a natural regulator of Meta-VCL is upregulated in the serum of patients with coronary artery disease

Background

Coronary artery disease (CAD), the leading cause of mortality globally, arises from atherosclerotic blockage of the coronary arteries. Meta-vinculin (meta-VCL), a large spliced isoform of VCL, co-localizes in muscular adhesive structures and plays significant roles in cardiac physiology and pathophysiology. This study aimed to identify microRNAs (miRNAs) regulating meta-VCL expression and investigate the expression alterations of the miRNAs of interest and meta-VCL as potential biomarkers in the serum of CAD patients.

Methods

Bioinformatics tools were employed to select miRNAs targeting meta-VCL. Cell-based ectopic expression analysis and a dual-luciferase assay were used to examine the interactions between miRNAs and meta-VCL. An ELISA assessed the concentrations of interleukin-6 (IL-6), IL-10, and tumor necrosis factor-α (TNF-α). MiRNA and meta-VCL expression patterns and biomarker suitability were evaluated in serum samples from CAD and non-CAD individuals using real-time PCR. A cardiac cell-line data set and CAD blood exosome samples were analyzed using bioinformatics and ROC curve analyses, respectively.

Results

miR-6721-5p directly interacted with the putative target sites at the 3'-UTR of meta-VCL and regulated its expression. IL-10 and TNF-α concentrations, which may act as anti-inflammatory factors, decreased following miR-6721-5p upregulation and meta-VCL downregulation. Bioinformatics and experimental expression analyses confirmed downregulated meta-VCL expression and upregulated miR-6721-5p expression in CAD samples. ROC curve analysis yielded an AUC score of 0.705 (P = 0.018), indicating the potential suitability of miR-6721-5p as a biomarker for CAD.

Conclusions

miR-6721-5p plays a regulatory role in meta-VCL expression and may contribute to CAD development by reducing anti-inflammatory factors. These findings suggest that miR-6721-5p could serve as a novel biomarker in the pathogenesis of CAD.

Impact of MICA 3'UTR allelic variability on miRNA binding prediction, a bioinformatic approach

MicroRNAs (miRNAs) are small non-coding RNAs that participate as powerful genetic regulators. MiRNAs can interfere with cellular processes by interacting with a broad spectrum of target genes under physiological and pathological states, including cancer development and progression. Major histocompatibility complex major histocompatibility complex class I-related chain A (MICA) belongs to a family of proteins that bind the natural-killer group 2, member D (NKG2D) receptor on Natural Killer cells and other cytotoxic lymphocytes. MICA plays a crucial role in the host's innate immune response to several disease settings, including cancer. MICA harbors various single nucleotide polymorphisms (SNPs) located in its 3'-untranslated region (3'UTR), a characteristic that increases the complexity of MICA regulation, favoring its post-transcriptional modulation by miRNAs under physiological and pathological conditions. Here, we conducted an in-depth analysis of MICA 3'UTR sequences according to each MICA allele described to date using NCBI database. We also systematically evaluated interactions between miRNAs and their putative targets on MICA 3'UTR containing SNPs using in silico analysis. Our in silico results showed that MICA SNPs rs9266829, rs 1880, and rs9266825, located in the target sequence of miRNAs hsa-miR-106a-5p, hsa-miR-17-5p, hsa-miR-20a-5p, hsa-miR-20b-5p, hsa-miR-93, hsa-miR-1207.5p, and hsa-miR-711 could modify the binding free energy between -8.62 and -18.14 kcal/mol, which may affect the regulation of MICA expression. We believe that our results may provide a starting point for further exploration of miRNA regulatory effects depending on MICA allelic variability; they may also be a guide to conduct miRNA in silico analysis for other highly polymorphic genes.

The Genome-Wide Characterization of Alternative Splicing and RNA Editing in the Development of Coprinopsis cinerea

Coprinopsis cinerea is one of the model species used in fungal developmental studies. This mushroom-forming Basidiomycetes fungus has several developmental destinies in response to changing environments, with dynamic developmental regulations of the organism. Although the gene expression in C. cinerea development has already been profiled broadly, previous studies have only focused on a specific stage or process of fungal development. A comprehensive perspective across different developmental paths is lacking, and a global view on the dynamic transcriptional regulations in the life cycle and the developmental paths is far from complete. In addition, knowledge on co- and post-transcriptional modifications in this fungus remains rare. In this study, we investigated the transcriptional changes and modifications in C. cinerea during the processes of spore germination, vegetative growth, oidiation, sclerotia formation, and fruiting body formation by inducing different developmental paths of the organism and profiling the transcriptomes using the high-throughput sequencing method. Transition in the identity and abundance of expressed genes drive the physiological and morphological alterations of the organism, including metabolism and multicellularity construction. Moreover, stage- and tissue-specific alternative splicing and RNA editing took place and functioned in C. cinerea. These modifications were negatively correlated to the conservation features of genes and could provide extra plasticity to the transcriptome during fungal development. We suggest that C. cinerea applies different molecular strategies in its developmental regulation, including shifts in expressed gene sets, diversifications of genetic information, and reversible diversifications of RNA molecules. Such features would increase the fungal adaptability in the rapidly changing environment, especially in the transition of developmental programs and the maintenance and balance of genetic and transcriptomic divergence. The multi-layer regulatory network of gene expression serves as the molecular basis of the functioning of developmental regulation.

Alternative splicing impacts microRNA regulation within coding regions

MicroRNAs (miRNAs) are small non-coding RNA molecules that bind to target sites in different gene regions and regulate post-transcriptional gene expression. Approximately 95% of human multi-exon genes can be spliced alternatively, which enables the production of functionally diverse transcripts and proteins from a single gene. Through alternative splicing, transcripts might lose the exon with the miRNA target site and become unresponsive to miRNA regulation. To check this hypothesis, we studied the role of miRNA target sites in both coding and non-coding regions using six cancer data sets from The Cancer Genome Atlas (TCGA) and Parkinson's disease data from PPMI. First, we predicted miRNA target sites on mRNAs from their sequence using TarPmiR. To check whether alternative splicing interferes with this regulation, we trained linear regression models to predict miRNA expression from transcript expression. Using nested models, we compared the predictive power of transcripts with miRNA target sites in the coding regions to that of transcripts without target sites. Models containing transcripts with target sites perform significantly better. We conclude that alternative splicing does interfere with miRNA regulation by skipping exons with miRNA target sites within the coding region.

Immunotherapeutic Significance of a Prognostic Alternative Splicing Signature in Bladder Cancer

Objectives: Bladder cancer is the fourth most common malignancy in men in the United States. Aberrant alternative splicing (AS) events are involved in the carcinogenesis, but the association between AS and bladder cancer remains unclear. This study aimed to construct an AS-based prognostic signature and elucidate the role of the tumor immune microenvironment (TIME) and the response to immunotherapy and chemotherapy in bladder cancer. Methods: Univariate Cox regression analysis was performed to detect prognosis-related AS events. The least absolute shrinkage and selection operator (LASSO) and multivariate Cox analyses were employed to build prognostic signatures. Kaplan-Meier survival analysis, multivariate Cox regression analysis, and receiver operating characteristic (ROC) curves were conducted to validate the prognostic signatures. Then, the Estimation of Stromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE) and tumor immune estimation resource (TIMER) databases were searched and the single-sample gene set enrichment analysis (ssGSEA) algorithm and CIBERSORT method were performed to uncover the context of TIME in bladder cancer. The Tumor Immune Dysfunction and Exclusion (TIDE) web tool and pRRophetic algorithm were used to predict the response to immunotherapy and chemotherapy. Finally, we constructed a correlation network between splicing factors (SFs) and survival-related AS events. Results: A total of 4684 AS events were significantly associated with overall survival in patients with bladder cancer. Eight prognostic signatures of bladder cancer were established, and a clinical survival prediction model was built. In addition, the consolidated prognostic signature was closely related to immune infiltration and the response to immunotherapy and chemotherapy. Furthermore, the correlation identified EIF3A, DDX21, SDE2, TNPO1, and RNF40 as hub SFs, and function analysis found ubiquitin-mediated proteolysis is correlated most significantly with survival-associated AS events. Conclusion: Our findings highlight the prognostic value of AS for patients with bladder cancer and reveal pivotal players of AS events in the context of TIME and the response to immunotherapy and chemotherapy, which may be important for patient management and treatment.

Colorectal cancer promoter methylation alteration affects the expression of glutamate ionotropic receptor AMPA type subunit 4 alternative isoforms potentially relevant in colon tissue

DNA methylation alterations are early events during tumourigenesis, affecting genes involved in the crosstalk between cells and surroundings in colorectal cancer (CRC). Among these genes, GRIA4, Glutamate Ionotropic Receptor AMPA Type Subunit 4, displays hypermethylation in the promoter region, and is an early diagnostic biomarker. It is well known that methylation can also affect alternative transcription. The purpose of this study is to evaluate the expression, at transcript and protein level, of GRIA4 main isoforms (the canonical one and a short variant) in 23 CRC and matched normal samples, of which we previously verified the methylation status. We further predicted miRNA/transcript target interactions as a possible post-transcriptional regulation using bioinformatics tools. As expected, downregulation of both variants has been observed in tumours. Interestingly, in contrast to what observed at transcriptional level, the GluR4 protein short isoform displayed higher expression than the canonical one either in normal or tumoural tissues. This may be explained by miRNA specifically targeting the canonical isoform. Our study is the first one that shows the expression of both isoforms in colon tissues. To note, the evident expression of the short isoform suggests a functional role in intestinal cell biology.

Epigenetic interplay between methylation and miRNA in bladder cancer: focus on isoform expression

BACKGROUND

Various epigenetic factors are responsible for the non-genetic regulation on gene expression. The epigenetically dysregulated oncogenes or tumor suppressors by miRNA and/or DNA methylation are often observed in cancer cells. Each of these epigenetic regulators has been studied well in cancer progressions; however, their mutual regulatory relationship in cancer still remains unclear. In this study, we propose an integrative framework to systematically investigate epigenetic interactions between miRNA and methylation at the alternatively spliced mRNA level in bladder cancer. Each of these epigenetic regulators has been studied well in cancer progressions; however, their mutual regulatory relationship in cancer still remains unclear.

RESULTS

The integrative analyses yielded 136 significant combinations (methylation, miRNA and isoform). Further, overall survival analysis on the 136 combinations based on methylation and miRNA, high and low expression groups resulted in 13 combinations associated with survival. Additionally, different interaction patterns were examined.

CONCLUSIONS

Our study provides a higher resolution of molecular insight into the crosstalk between two epigenetic factors, DNA methylation and miRNA. Given the importance of epigenetic interactions and alternative splicing in cancer, it is timely to identify and understand the underlying mechanisms based on epigenetic markers and their interactions in cancer, leading to alternative splicing with primary functional impact.

The microRNA target site landscape is a novel molecular feature associating alternative polyadenylation with immune evasion activity in breast cancer

Alternative polyadenylation (APA) in breast tumor samples results in the removal/addition of cis-regulatory elements such as microRNA (miRNA) target sites in the 3'-untranslated region (3'-UTRs) of genes. Although previous computational APA studies focused on a subset of genes strongly affected by APA (APA genes), we identify miRNAs of which widespread APA events collectively increase or decrease the number of target sites [probabilistic inference of microRNA target site modification through APA (PRIMATA-APA)]. Using PRIMATA-APA on the cancer genome atlas (TCGA) breast cancer data, we found that the global APA events change the number of the target sites of particular microRNAs [target sites modified miRNA (tamoMiRNA)] enriched for cancer development and treatments. We also found that when knockdown (KD) of NUDT21 in HeLa cells induces a different set of widespread 3'-UTR shortening than TCGA breast cancer data, it changes the target sites of the common tamoMiRNAs. Since the NUDT21 KD experiment previously demonstrated the tumorigenic role of APA events in a miRNA dependent fashion, this result suggests that the APA-initiated tumorigenesis is attributable to the miRNA target site changes, not the APA events themselves. Further, we found that the miRNA target site changes identify tumor cell proliferation and immune cell infiltration to the tumor microenvironment better than the miRNA expression levels or the APA events themselves. Altogether, our computational analyses provide a proof-of-concept demonstration that the miRNA target site information indicates the effect of global APA events with a potential as predictive biomarker.

KH domain containing RNA-binding proteins coordinate with microRNAs to regulate Caenorhabditis elegans development

MicroRNAs (miRNAs) and RNA-binding proteins (RBPs) regulate gene expression at the post-transcriptional level, but the extent to which these key regulators of gene expression coordinate their activities and the precise mechanisms of this coordination are not well understood. RBPs often have recognizable RNA binding domains that correlate with specific protein function. Recently, several RBPs containing K homology (KH) RNA binding domains were shown to work with miRNAs to regulate gene expression, raising the possibility that KH domains may be important for coordinating with miRNA pathways in gene expression regulation. To ascertain whether additional KH domain proteins functionally interact with miRNAs during Caenorhabditis elegans development, we knocked down twenty-four genes encoding KH-domain proteins in several miRNA sensitized genetic backgrounds. Here, we report that a majority of the KH domain-containing genes genetically interact with multiple miRNAs and Argonaute alg-1. Interestingly, two KH domain genes, predicted splicing factors sfa-1 and asd-2, genetically interacted with all of the miRNA mutants tested, whereas other KH domain genes showed genetic interactions only with specific miRNAs. Our domain architecture and phylogenetic relationship analyses of the C. elegans KH domain-containing proteins revealed potential groups that may share both structure and function. Collectively, we show that many C. elegans KH domain RBPs functionally interact with miRNAs, suggesting direct or indirect coordination between these two classes of post-transcriptional gene expression regulators.

Frequency of MicroRNA Response Elements Identifies Pathologically Relevant Signaling Pathways in Triple-Negative Breast Cancer

Complex interactions between mRNAs and microRNAs influence cellular functions. The mRNA-microRNA interactions also determine the post-transcriptional availability of mRNAs and unbound microRNAs. MicroRNAs binds to one or more microRNA response elements (MREs) located on the 3′UTR of mRNAs. In this study, we leveraged MREs and their frequencies in cancer and matched normal tissues to obtain insights into disease-specific interactions between mRNAs and microRNAs. We developed a bioinformatics method “ReMIx” that utilizes RNA sequencing (RNA-Seq) data to quantify MRE frequencies across the transcriptome. We applied ReMIx to triple-negative (TN) breast cancer tumor-normal adjacent pairs and identified MREs specific to TN tumors. ReMIx identified candidate mRNAs and microRNAs in the MAPK signaling cascade. Further analysis of MAPK gene regulatory networks revealed microRNA partners that influence and modulate MAPK signaling. In conclusion, we demonstrate a novel method of using MREs in the identification of functionally relevant mRNA-microRNA interactions in TN breast cancer.

•

Bioinformatics method ReMIx identify differential microRNA response rlements (MRE)

•

Tumor-specific MREs frequency observed in triple-negative breast cancer (TNBC)

•

MRE analysis identify MAPK signaling genes as therapeutic target for TNBC

•

MREs frequency can be used to identify pathologically relevant pathways

WNT gene polymorphisms and predisposition to apical periodontitis

Single nucleotide polymorphisms (SNPs) in WNT genes may impact gene/protein function and contribute to individual predisposition to apical periodontitis (AP). Here, we investigated the association of SNPs in/nearby WNT3, WNT3A, WNT5A, WNT8A, WNT9B and WNT11 genes with AP using a case-control dataset. Cases were defined as individuals with deep caries and AP (n = 188); controls had deep caries and no AP (n = 230). Genotyping was performed using Taqman chemistry in real time PCR. Data analyses was performed using Fisher Exact tests assuming a Bonferroni correction threshold value of 0.005. Single-SNP association analysis revealed a trend for association with WNT3 rs9890413 genotypes (P = 0.009) under a dominant model and allelic association for WNT3A rs1745420 (P = 0.009). Haplotypes involving WNT3-WNT9B-WNT3A alleles were also significantly associated with AP (P ≤ 0.003). Luciferase reporter assays showed higher transcriptional activity (1.4-fold) with the alternate G allele in rs1745420. Expression of WNT3, WNT3A and WNT5A in AP tissues was significantly higher than in control tissues, and inversely correlated with the expression of SERPINB1, COL1A1 and TIMP1 (P < 0.05). Our results suggest that WNT genes have a role in modulating AP and polymorphisms in these genes may increase susceptibility to AP.

Aberrant DNA methylation defines isoform usage in cancer, with functional implications

Alternative transcript isoforms are common in tumors and act as potential drivers of cancer. Mechanisms determining altered isoform expression include somatic mutations in splice regulatory sites or altered splicing factors. However, since DNA methylation is known to regulate transcriptional isoform activity in normal cells, we predicted the highly dysregulated patterns of DNA methylation present in cancer also affect isoform activity. We analyzed DNA methylation and RNA-seq isoform data from 18 human cancer types and found frequent correlations specifically within 11 cancer types. Examining the top 25% of variable methylation sites revealed that the location of the methylated CpG site in a gene determined which isoform was used. In addition, the correlated methylation-isoform patterns classified tumors into known subtypes and predicted distinct protein functions between tumor subtypes. Finally, methylation-correlated isoforms were enriched for oncogenes, tumor suppressors, and cancer-related pathways. These findings provide new insights into the functional impact of dysregulated DNA methylation in cancer and highlight the relationship between the epigenome and transcriptome.

In eukaryotes, one gene can be transcribed into multiple RNA sequences (or isoforms) that are subsequently translated into proteins with different functions in response to specific cellular needs. Recent studies showed that cancer cells can obtain abnormal functions via expressing different isoforms. In normal cells, isoform expression can be regulated by DNA methylation–a molecular signature with attached methyl groups on DNA sequences. Given that dysregulation of DNA methylation is a cancer hallmark, we suspect the same regulation holds in cancer and contributes to cancer progression. In this study, we analyzed data from 18 human cancer types and found frequent correlations in 11 cancer types between specific isoform usage and DNA methylation depending on the location of the methylated site in a gene. These correlation patterns can classify heterogeneous tumors in a cancer type into homogeneous subtypes and are predicted to change protein functions via isoform switching between subtypes. Finally, we found cancer-related genes often harbored more DNA methylation-isoform correlations than genes not implicated in cancer. This finding could help us to better understand the functional impact of DNA methylation alterations via regulation of isoform expression in tumorigenesis and to further improve the cancer treatment.

Non-coding transcript variants of protein-coding genes - what are they good for?

The total number of protein-coding genes in the human genome is not significantly higher than those in much simpler eukaryotes, despite a general increase in genome size proportionate to the organismal complexity. The large non-coding transcriptome and extensive differential splicing, are increasingly being accepted as the factors contributing to the complex mammalian physiology and architecture. Recent studies reveal additional layers of functional complexity: some long non-coding RNAs have been re-defined as micropeptide or microprotein encoding transcripts, and in turn some protein-coding RNAs are bifunctional and display also non-coding functions. Moreover, several protein-coding genes express long non-coding RNA splice-forms and generate circular RNAs in addition to their canonical mRNA transcripts, revoking the strict definition of a gene as coding or non-coding. In this mini review, we discuss the current understanding of these hybrid genes and their possible roles and relevance.