An alternative splicing switch in FLNB promotes the mesenchymal cell state in human breast cancer

Characterization of Alternative Splicing in High-Risk Wilms' Tumors

The significant heterogeneity of Wilms' tumors between different patients is thought to arise from genetic and epigenetic distortions that occur during various stages of fetal kidney development in a way that is poorly understood. To address this, we characterized the heterogeneity of alternative mRNA splicing in Wilms' tumors using a publicly available RNAseq dataset of high-risk Wilms' tumors and normal kidney samples. Through Pareto task inference and cell deconvolution, we found that the tumors and normal kidney samples are organized according to progressive stages of kidney development within a triangle-shaped region in latent space, whose vertices, or "archetypes", resemble the cap mesenchyme, the nephrogenic stroma, and epithelial tubular structures of the fetal kidney. We identified a set of genes that are alternatively spliced between tumors located in different regions of latent space and found that many of these genes are associated with the epithelial-to-mesenchymal transition (EMT) and muscle development. Using motif enrichment analysis, we identified putative splicing regulators, some of which are associated with kidney development. Our findings provide new insights into the etiology of Wilms' tumors and suggest that specific splicing mechanisms in early stages of development may contribute to tumor development in different patients.

Alternative splicing in EMT and TGF-β signaling during cancer progression

Epithelial to mesenchymal transition (EMT) is a physiological process during development where epithelial cells transform to acquire mesenchymal characteristics, which allows them to migrate and colonize secondary tissues. Many cellular signaling pathways and master transcriptional factors exert a myriad of controls to fine tune this vital process to meet various developmental and physiological needs. Adding to the complexity of this network are post-transcriptional and post-translational regulations. Among them, alternative splicing has been shown to play important roles to drive EMT-associated phenotypic changes, including actin cytoskeleton remodeling, cell-cell junction changes, cell motility and invasiveness. In advanced cancers, transforming growth factor-β (TGF-β) is a major inducer of EMT and is associated with tumor cell metastasis, cancer stem cell self-renewal, and drug resistance. This review aims to provide an overview of recent discoveries regarding alternative splicing events and the involvement of splicing factors in the EMT and TGF-β signaling. It will emphasize the importance of various splicing factors involved in EMT and explore their regulatory mechanisms.

NONO promotes gallbladder cancer cell proliferation by enhancing oncogenic RNA splicing of DLG1 through interaction with IGF2BP3/RBM14

Gallbladder cancer (GBC) is a highly malignant and rapidly progressing tumor of the human biliary system, and there is an urgent need to develop new therapeutic targets and modalities. Non-POU domain-containing octamer-binding protein (NONO) is an RNA-binding protein involved in the regulation of transcription, mRNA splicing, and DNA repair. NONO expression is elevated in multiple tumors and can act as an oncogene to promote tumor progression. Here, we found that NONO was highly expressed in GBC and promoted tumor cells growth. The dysregulation of RNA splicing is a molecular feature of almost all tumor types. Accordingly, mRNA-seq and RIP-seq analysis showed that NONO promoted exon6 skipping in DLG1, forming two isomers (DLG1-FL and DLG1-S). Furthermore, lower Percent-Spliced-In (PSI) values of DLG1 were detected in tumor tissue relative to the paraneoplastic tissue, and were associated with poor patient prognosis. Moreover, DLG1-S and DLG1-FL act as tumor promoters and tumor suppressors, respectively, by regulating the YAP1/JUN pathway. N6-methyladenosine (m6A) is the most common and abundant RNA modification involved in alternative splicing processes. We identified an m6A reader, IGF2BP3, which synergizes with NONO to promote exon6 skipping in DLG1 in an m6A-dependent manner. Furthermore, IP/MS results showed that RBM14 was bound to NONO and interfered with NONO-mediated exon6 skipping of DLG1. In addition, IGF2BP3 disrupted the binding of RBM14 to NONO. Overall, our data elucidate the molecular mechanism by which NONO promotes DLG1 exon skipping, providing a basis for new therapeutic targets in GBC treatment.

Dual-Specificity Tyrosine Phosphorylation-Regulated Kinase 3 Expression and Its Correlation with Prognosis and Growth of Serous Ovarian Cancer: Correlation of DYRK3 with Ovarian Cancer Survival

Background

Epithelial ovarian cancer, primarily serous ovarian cancer (SOC), stands as a predominant cause of cancer-related mortality among women globally, emphasizing the urgent need for comprehensive research into its molecular underpinnings. Within this context, the dual-specificity tyrosine phosphorylation-regulated kinase 3 (DYRK3) has emerged as a potential key player with implications for prognosis and tumor progression.

Methods

This study conducted a meticulous retrospective analysis of 254 SOC cases from our medical center to unravel the prognostic significance of DYRK3. Survival analyses underscored DYRK3 as an independent adverse prognostic factor in SOC, with a hazard ratio of 2.60 (95% CI 1.67-4.07, P < 0.001). Experimental investigations involved DYRK3 knockdown in serous ovarian cancer cell lines (CAOV3 and OVCAR-3) through a shRNA strategy, revealing substantial decreases in cell growth and invasion capabilities. Bioinformatics analyses further hinted at DYRK3's involvement in modulating the tumor immune microenvironment. In vivo experiments with DYRK3-knockdown cell lines validated these findings, demonstrating a notable restriction in the growth of ovarian cancer xenografts.

Results

Our findings collectively illuminate DYRK3 as a pivotal tumor-promoting oncogene in SOC. Beyond its adverse prognostic implications, DYRK3 knockdown exhibited promising therapeutic potential by impeding cancer progression and potentially influencing the tumor immune microenvironment.

Conclusions

This study establishes a compelling foundation for further research into DYRK3's intricate role and therapeutic potential in ovarian cancer treatment. As we unravel the complexities surrounding DYRK3, our work not only contributes to the understanding of SOC pathogenesis but also unveils new prospects for targeted therapeutic interventions, holding promise for improved outcomes in ovarian cancer management.

Immunopeptidomics-based identification of naturally presented non-canonical circRNA-derived peptides

Circular RNAs (circRNAs) are covalently closed non-coding RNAs lacking the 5' cap and the poly-A tail. Nevertheless, it has been demonstrated that certain circRNAs can undergo active translation. Therefore, aberrantly expressed circRNAs in human cancers could be an unexplored source of tumor-specific antigens, potentially mediating anti-tumor T cell responses. This study presents an immunopeptidomics workflow with a specific focus on generating a circRNA-specific protein fasta reference. The main goal of this workflow is to streamline the process of identifying and validating human leukocyte antigen (HLA) bound peptides potentially originating from circRNAs. We increase the analytical stringency of our workflow by retaining peptides identified independently by two mass spectrometry search engines and/or by applying a group-specific FDR for canonical-derived and circRNA-derived peptides. A subset of circRNA-derived peptides specifically encoded by the region spanning the back-splice junction (BSJ) are validated with targeted MS, and with direct Sanger sequencing of the respective source transcripts. Our workflow identifies 54 unique BSJ-spanning circRNA-derived peptides in the immunopeptidome of melanoma and lung cancer samples. Our approach enlarges the catalog of source proteins that can be explored for immunotherapy.

Identification of somatic mutation-driven enhancers and their clinical utility in breast cancer

Somatic mutations contribute to cancer development by altering the activity of enhancers. In the study, a total of 135 mutation-driven enhancers, which displayed significant chromatin accessibility changes, were identified as candidate risk factors for breast cancer (BRCA). Furthermore, we identified four mutation-driven enhancers as independent prognostic factors for BRCA subtypes. In Her2 subtype, enhancer G > C mutation was associated with poorer prognosis through influencing its potential target genes FBXW9, TRIR, and WDR83. We identified aminoglutethimide and quinpirole as candidate drugs targeting the mutated enhancer. In normal subtype, enhancer G > A mutation was associated with poorer prognosis through influencing its target genes ALOX15B, LINC00324, and MPDU1. We identified eight candidate drugs such as erastin, colforsin, and STOCK1N-35874 targeting the mutated enhancer. Our findings suggest that somatic mutations contribute to breast cancer subtype progression by altering enhancer activity, which could be potential candidates for cancer therapy.

Epigenetic regulation of mRNA mediates the phenotypic plasticity of cancer cells during metastasis and therapeutic resistance (Review)

Plasticity, the ability of cancer cells to transition between differentiation states without genomic alterations, has been recognized as a major source of intratumoral heterogeneity. It has a crucial role in cancer metastasis and treatment resistance. Thus, targeting plasticity holds tremendous promise. However, the molecular mechanisms of plasticity in cancer cells remain poorly understood. Several studies found that mRNA, which acts as a bridge linking the genetic information of DNA and protein, has an important role in translating genotypes into phenotypes. The present review provided an overview of the regulation of cancer cell plasticity occurring via changes in the transcription and editing of mRNAs. The role of the transcriptional regulation of mRNA in cancer cell plasticity was discussed, including DNA‑binding transcriptional factors, DNA methylation, histone modifications and enhancers. Furthermore, the role of mRNA editing in cancer cell plasticity was debated, including mRNA splicing and mRNA modification. In addition, the role of non‑coding (nc)RNAs in cancer plasticity was expounded, including microRNAs, long intergenic ncRNAs and circular RNAs. Finally, different strategies for targeting cancer cell plasticity to overcome metastasis and therapeutic resistance in cancer were discussed.

Quaking isoforms cooperate to promote the mesenchymal phenotype

The RNA-binding protein Quaking (QKI) has widespread effects on mRNA regulation including alternative splicing, stability, translation, and localization of target mRNAs. Recently, QKI was found to be induced during epithelial-mesenchymal transition (EMT), where it promotes a mesenchymal alternative splicing signature that contributes to the mesenchymal phenotype. QKI is itself alternatively spliced to produce three major isoforms, QKI-5, QKI-6, and QKI-7. While QKI-5 is primarily localized to the nucleus where it controls mesenchymal splicing during EMT, the functions of the two predominantly cytoplasmic isoforms, QKI-6 and QKI-7, in this context remain uncharacterized. Here we used CRISPR-mediated depletion of QKI in a human mammary epithelial cell model of EMT and studied the effects of expressing the QKI isoforms in isolation and in combination. QKI-5 was required to induce mesenchymal morphology, while combined expression of QKI-5 with either QKI-6 or QKI-7 further enhanced mesenchymal morphology and cell migration. In addition, we found that QKI-6 and QKI-7 can partially localize to the nucleus and contribute to alternative splicing of QKI target genes. These findings indicate that the QKI isoforms function in a dynamic and cooperative manner to promote the mesenchymal phenotype.

Genome-Related Mechanisms Contributing to Differences in Alzheimer's Disease Incidence Between White and Black Older US Adults

In this manuscript, we leverage a modified GWAS algorithm adapted for use with multidimensional Cox models and data from the Health and Retirement Study to explore how genetic variation influences the size of the disparity in Alzheimer's disease (AD) incidence between older Black and White US adults. We identified four loci that were associated with higher AD incidence levels in older Black adults: (1) rs11077034 (hazard ratio (HR), 4.98) from the RBFOX1 gene; (2) rs7144494 (HR, 1.68) from the HISLA gene; (3) rs7660552 (HR, 3.07) from the SLC25A4 gene; and (4) rs12599485 (HR, 3.181) from the NIP30 gene. The RBFOX1, HISLA, SLC25A4, and NIP30 genes are known to be associated with AD (RBFOX1, NIP30) directly, and also influence the risk of AD risk-related morbidities such as hypertension (RBFOX1, SLC25A4), depression (SLC25A4), and certain cancers (HISLA, SLC25A4). A likely disparity-generating mechanism may lie in endocytosis and abnormal tissue growing mechanisms, depending on inherited gene mutations and the structure of proxies as well as gene-environment interactions with other risk factors such as lifestyle, education level, and access to adequate medical services.

The landscape of alternative polyadenylation during EMT and its regulation by the RNA-binding protein Quaking

Epithelial-mesenchymal transition (EMT) plays important roles in tumour progression and is orchestrated by dynamic changes in gene expression. While it is well established that post-transcriptional regulation plays a significant role in EMT, the extent of alternative polyadenylation (APA) during EMT has not yet been explored. Using 3' end anchored RNA sequencing, we mapped the alternative polyadenylation (APA) landscape following Transforming Growth Factor (TGF)-β-mediated induction of EMT in human mammary epithelial cells and found APA generally causes 3'UTR lengthening during this cell state transition. Investigation of potential mediators of APA indicated the RNA-binding protein Quaking (QKI), a splicing factor induced during EMT, regulates a subset of events including the length of its own transcript. Analysis of QKI crosslinked immunoprecipitation (CLIP)-sequencing data identified the binding of QKI within 3' untranslated regions (UTRs) was enriched near cleavage and polyadenylation sites. Following QKI knockdown, APA of many transcripts is altered to produce predominantly shorter 3'UTRs associated with reduced gene expression. These findings reveal the changes in APA that occur during EMT and identify a potential role for QKI in this process.

RBFOX2 deregulation promotes pancreatic cancer progression and metastasis through alternative splicing

RNA splicing is an important biological process associated with cancer initiation and progression. However, the contribution of alternative splicing to pancreatic cancer (PDAC) development is not well understood. Here, we identify an enrichment of RNA binding proteins (RBPs) involved in splicing regulation linked to PDAC progression from a forward genetic screen using Sleeping Beauty insertional mutagenesis in a mouse model of pancreatic cancer. We demonstrate downregulation of RBFOX2, an RBP of the FOX family, promotes pancreatic cancer progression and liver metastasis. Specifically, we show RBFOX2 regulates exon splicing events in transcripts encoding proteins involved in cytoskeletal remodeling programs. These exons are differentially spliced in PDAC patients, with enhanced exon skipping in the classical subtype for several RBFOX2 targets. RBFOX2 mediated splicing of ABI1, encoding the Abelson-interactor 1 adapter protein, controls the abundance and localization of ABI1 protein isoforms in pancreatic cancer cells and promotes the relocalization of ABI1 from the cytoplasm to the periphery of migrating cells. Using splice-switching antisense oligonucleotides (AONs) we demonstrate the ABI1 ∆Ex9 isoform enhances cell migration. Together, our data identify a role for RBFOX2 in promoting PDAC progression through alternative splicing regulation.

RNA-binding proteins regulating the CD44 alternative splicing

Alternative splicing is often deregulated in cancer, and cancer-specific isoform switches are part of the oncogenic transformation of cells. Accumulating evidence indicates that isoforms of the multifunctional cell-surface glycoprotein CD44 play different roles in cancer cells as compared to normal cells. In particular, the shift of CD44 isoforms is required for epithelial to mesenchymal transition (EMT) and is crucial for the maintenance of pluripotency in normal human cells and the acquisition of cancer stem cells phenotype for malignant cells. The growing and seemingly promising use of splicing inhibitors for treating cancer and other pathologies gives hope for the prospect of using such an approach to regulate CD44 alternative splicing. This review integrates current knowledge about regulating CD44 alternative splicing by RNA-binding proteins.

Truncated SCRIB isoform promotes breast cancer metastasis through HNRNP A1 mediated exon 16 skipping

Breast cancer is one of the most common malignant tumors with high mortality due to metastases. SCRIB, a scaffold protein mainly distributed in the cell membrane, is a potential tumor suppressor. Mislocalization and aberrant expression of SCRIB stimulate the EMT pathway and promote tumor cell metastasis. SCRIB has two isoforms (with or without exon 16) produced by alternative splicing. In this study we investigated the function of SCRIB isoforms in breast cancer metastasis and their regulatory mechanisms. We showed that in contrast to the full-length isoform (SCRIB-L), the truncated SCRIB isoform (SCRIB-S) was overexpressed in highly metastatic MDA-MB-231 cells that promoted breast cancer metastasis through activation of the ERK pathway. The affinity of SCRIB-S for the catalytic phosphatase subunit PPP1CA was lower than that of SCRIB-L and such difference might contribute to the different function of the two isoforms in cancer metastasis. By conducting CLIP, RIP and MS2-GFP-based experiments, we revealed that the heterogeneous nuclear ribonucleoprotein A1 (hnRNP A1) promoted SCRIB exon 16 skipping by binding to the "AG"-rich sequence "caggauggaggccccccgugccgag" on intron 15 of SCRIB. Transfection of MDA-MB-231 cells with a SCRIB antisense oligodeoxynucleotide (ASO-SCRIB) designed on the basis of this binding sequence, not only effectively inhibited the binding of hnRNP A1 to SCRIB pre-mRNA and suppressed the production of SCRIB-S, but also reversed the activation of the ERK pathway by hnRNP A1 and inhibited the metastasis of breast cancer. This study provides a new potential target and a candidate drug for treating breast cancer.

Cell-type specific regulator RBPMS switches alternative splicing via higher-order oligomerization and heterotypic interactions with other splicing regulators

Alternative pre-mRNA splicing decisions are regulated by RNA binding proteins (RBPs) that can activate or repress regulated splice sites. Repressive RBPs typically harness multivalent interactions to bind stably to target RNAs. Multivalency can be achieved by homomeric oligomerization and heteromeric interactions with other RBPs, often mediated by intrinsically disordered regions (IDRs), and by possessing multiple RNA binding domains. Cell-specific splicing decisions often involve the action of widely expressed RBPs, which are able to bind multivalently around target exons, but without effect in the absence of a cell-specific regulator. To address how cell-specific regulators can collaborate with constitutive RBPs in alternative splicing regulation, we used the smooth-muscle specific regulator RBPMS. Recombinant RBPMS is sufficient to confer smooth muscle cell specific alternative splicing of Tpm1 exon 3 in cell-free assays by preventing assembly of ATP-dependent splicing complexes. This activity depends upon a C-terminal IDR that facilitates dynamic higher-order self-assembly, cooperative binding to multivalent RNA and interactions with widely expressed splicing co-regulators, including MBNL1 and RBFOX2, allowing cooperative assembly of stable cell-specific regulatory complexes.

Epithelial specific splicing regulator proteins as emerging oncogenes in aggressive prostate cancer

Prostate cancer progression is connected to the activity of conventional oncogenes and tumour suppressors and driven by circulating steroid hormones. A key issue has been how to identify and care for aggressively developing prostate tumours. Here we discuss how expression of the splicing regulators ESRP1 and ESRP2, and how their role as "masterminds" of epithelial splicing patterns, have been identified as markers of aggressively proliferating prostate primary tumours. We suggest that the origin of prostate cancer within epithelial cells, and the subsequent association of ESRP1 and ESRP2 expression with more aggressive disease progression, identify ESRP1 and ESRP2 as lineage survival oncogenes. To move this field on in the future it will be important to identify the gene expression targets controlled by ESRP1/2 that regulate prostate cancer proliferation. Potential future therapies could be designed to target ESRP1 and ESRP2 protein activity or their regulated splice isoforms in aggressive prostate tumours. Design of these therapies is potentially complicated by the risk of producing a more mesenchymal splicing environment that might promote tumour metastasis.

Highly multiplexed mRNA quantitation with CRISPR-Cas13

RNA quantitation tools are often either high-throughput or cost-effective, but rarely are they both. Existing methods can profile the transcriptome at great expense or are limited to quantifying a handful of genes by labor constraints. A technique that permits more throughput at a reduced cost could enable multi-gene kinetic studies, gene regulatory network analysis, and combinatorial genetic screens. Here, we introduce quantitative Combinatorial Arrayed Reactions for Multiplexed Evaluation of Nucleic acids (qCARMEN): an RNA quantitation technique which leverages the programmable RNA-targeting capabilities of CRISPR-Cas13 to address this challenge by quantifying over 4,500 gene-sample pairs in a single experiment. Using qCARMEN, we studied the response profiles of interferon-stimulated genes (ISGs) during interferon (IFN) stimulation and flavivirus infection. Additionally, we observed isoform switching kinetics during epithelial-mesenchymal transition. qCARMEN is a simple and inexpensive technique that greatly enhances the scalability of RNA quantitation for novel applications with performance similar to gold-standard methods.

MEN1 is a regulator of alternative splicing and prevents R-loop-induced genome instability through suppression of RNA polymerase II elongation

The fidelity of alternative splicing (AS) patterns is essential for growth development and cell fate determination. However, the scope of the molecular switches that regulate AS remains largely unexplored. Here we show that MEN1 is a previously unknown splicing regulatory factor. MEN1 deletion resulted in reprogramming of AS patterns in mouse lung tissue and human lung cancer cells, suggesting that MEN1 has a general function in regulating alternative precursor mRNA splicing. MEN1 altered exon skipping and the abundance of mRNA splicing isoforms of certain genes with suboptimal splice sites. Chromatin immunoprecipitation and chromosome walking assays revealed that MEN1 favored the accumulation of RNA polymerase II (Pol II) in regions encoding variant exons. Our data suggest that MEN1 regulates AS by slowing the Pol II elongation rate and that defects in these processes trigger R-loop formation, DNA damage accumulation and genome instability. Furthermore, we identified 28 MEN1-regulated exon-skipping events in lung cancer cells that were closely correlated with survival in patients with lung adenocarcinoma, and MEN1 deficiency sensitized lung cancer cells to splicing inhibitors. Collectively, these findings led to the identification of a novel biological role for menin in maintaining AS homeostasis and link this role to the regulation of cancer cell behavior.

circRNA-SFMBT2 orchestrates ERα activation to drive tamoxifen resistance in breast cancer cells

Dysregulated ERα signaling is responsible for endocrine resistance and eventual relapse in patients with estrogen receptor-positive (ER⁺) breast cancer. Thus, identifying novel ERα regulators is necessary to fully understand the mechanisms of endocrine resistance. Here, we identified circRNA-SFMBT2 to be highly expressed in ER⁺ breast cancer cells in comparison to ER^- cells and found that high circRNA-SFMBT2 levels were related to larger tumor size and poor prognosis in patients with ER⁺ breast cancer. In vitro and in vivo experiments confirmed that the circRNA-SFMBT2 level was positively correlated with the ERα protein level, implying a regulatory role for circRNA-SFMBT2 in ERα signaling. Moreover, we found that circRNA-SFMBT2 biogenesis could be facilitated via RNA-binding protein quaking (QKI), and biologically elevated circRNA-SFMBT2 expression promoted cell growth and tamoxifen resistance in ER⁺ breast cancer. Mechanistically, circRNA-SFMBT2 exhibits a specific tertiary structure that endows it with a high binding affinity for ERα and allows it to interact with the AF2 and DBD domains of ERα, enforcing recruitment of RNF181 to the AF1 domain of ERα. Furthermore, the circRNA-SFMBT2/RNF181 axis differentially regulated K48-linked and K63-linked ubiquitination of ERα to enhance ERα stability, resulting in increased expression of ERα target genes and tumor progression. In summary, circRNA-SFMBT2 is an important regulator of ERα signaling, and antagonizing circRNA-SFMBT2 expression may constitute a potential therapeutic strategy for breast cancer.

The possible molecular mechanism underlying the involvement of the variable shear factor QKI in the epithelial-mesenchymal transformation of oesophageal cancer

OBJECTIVE

Based on the GEO, TCGA and GTEx databases, we reveal the possible molecular mechanism of the variable shear factor QKI in epithelial mesenchymal transformation (EMT) of oesophageal cancer.

METHODS

Based on the TCGA and GTEx databases, the differential expression of the variable shear factor QKI in oesophageal cancer samples was analysed, and functional enrichment analysis of QKI was performed based on the TCGA-ESCA dataset. The percent-spliced in (PSI) data of oesophageal cancer samples were downloaded from the TCGASpliceSeq database, and the genes and variable splicing types that were significantly related to the expression of the variable splicing factor QKI were screened out. We further identified the significantly upregulated circRNAs and their corresponding coding genes in oesophageal cancer, screened the EMT-related genes that were significantly positively correlated with QKI expression, predicted the circRNA-miRNA binding relationship through the circBank database, predicted the miRNA-mRNA binding relationship through the TargetScan database, and finally obtained the circRNA-miRNA-mRNA network through which QKI promoted the EMT process.

RESULTS

Compared with normal control tissue, QKI expression was significantly upregulated in tumour tissue samples of oesophageal cancer patients. High expression of QKI may promote the EMT process in oesophageal cancer. QKI promotes hsa_circ_0006646 and hsa_circ_0061395 generation by regulating the variable shear of BACH1 and PTK2. In oesophageal cancer, QKI may promote the production of the above two circRNAs by regulating variable splicing, and these circRNAs further competitively bind miRNAs to relieve the targeted inhibition of IL-11, MFAP2, MMP10, and MMP1 and finally promote the EMT process.

CONCLUSION

Variable shear factor QKI promotes hsa_circ_0006646 and hsa_circ_0061395 generation, and downstream related miRNAs can relieve the targeted inhibition of EMT-related genes (IL11, MFAP2, MMP10, MMP1) and promote the occurrence and development of oesophageal cancer, providing a new theoretical basis for screening prognostic markers of oesophageal cancer patients.

Disequilibrium between BRCA1 and BRCA2 Circular and Messenger RNAs Plays a Role in Breast Cancer

Breast cancer is a frequent disease for which the discovery of markers that enable early detection or prognostic assessment remains challenging. Circular RNAs (circRNAs) are single-stranded structures in closed loops that are produced by backsplicing. CircRNA and messenger RNA (mRNA) are generated co-transcriptionally, and backsplicing and linear splicing compete against each other. As mRNAs are key players in tumorigenesis, we hypothesize that a disruption of the balance between circRNAs and mRNAs could promote breast cancer. Hence, we developed an assay for a simultaneous study of circRNAs and mRNAs, which we have called splice and expression analyses by exon ligation and high-throughput sequencing (SEALigHTS). Following SEALigHTS validation for BRCA1 and BRCA2, our hypothesis was tested using an independent research set of 95 pairs from tumor and adjacent normal breast tissues. In this research set, ratios of BRCA1 and BRCA2 circRNAs/mRNAs were significantly lower in the tumor breast tissue compared to normal tissue (p = 1.6 × 10^-9 and p = 4.4 × 10^-5 for BRCA1 and BRCA2, respectively). Overall, we developed an innovative method to study linear splicing and backsplicing, described the repertoire of BRCA1 and BRCA2 circRNAs, including 15 novel ones, and showed for the first time that a disequilibrium between BRCA1 and BRCA2 circRNAs and mRNAs plays a role in breast cancer.

Role of identified proteins in the proteome profiles of CDK4/6 inhibitor-resistant breast cancer cell lines

Abemaciclib (Ab) and palbociclib (Pb) are CDK4/6 inhibitors used to cure advanced breast cancer (BC). However, acquired resistance is a major challenge. The molecular mechanisms and signature proteins of therapy resistance for Ab and Pb drugs need to be explored. Here we developed resistant cells for Ab and Pb drugs in MCF-7 cell lines and explored the mechanisms and signature proteins of therapy resistance in BC. Proteome profiling was performed using the label-free proteome-orbitrap-fusion-MS-MS technique. Gene ontology (GO)-terms, KEGG pathways and network analysis were performed for the proteome data. Drug-resistant cells showed increased drug tolerance, enhanced colony formation potential and an increased gap-healing tendency for the respective drug. Up-regulation of survival genes (BCL-2 and MCL-1) and down-regulation of apoptosis inducers were observed. Drug-resistance markers (MDR-1 and ABCG2 (BCRP)) along with ESR-1, CDK4, CDK6, and cyclin-D1 genes were up-regulated in resistant cells. A total of 237 and 239 proteins were found to be differentially expressed in the Ab and Pb-resistant cells, respectively. Down-regulated proteins induce apoptosis signalling and nucleotide metabolisms and restrict EGFR signalling; however, up-regulated proteins induce Erk, wnt-β-catenin, VEGFR-PI3K-AKT, glucose transportation, and hypoxia signalling pathways and regulate hydrogen peroxide signalling pathways. The panel of identified proteins associated with these pathways might have characteristics of molecular signature and new drug targets for overcoming drug resistance in breast cancer.

The Pattern of RNA Editing Changes in Pleural Mesothelioma upon Epithelial-Mesenchymal Transition

Pleural mesothelioma (PM) is a cancer where epithelioid, biphasic and sarcomatoid histotypes are observed. Sarcomatoid PM is characterized by mesenchymal features. Multi-omics have been used to characterize the epithelial-to-mesenchymal (EMT) phenotype at the molecular level. We contribute to this effort by including the analysis of RNA editing. We extracted samples with the highest vs. lowest Epithelial score from two PM cohorts and observed increased RNA editing in introns and decreased RNA editing in 3'UTR upon EMT. The same was observed in primary PM primary cultures stratified by transcriptomics analysis into two groups, one of them enriched with mesenchymal features. Our data demonstrate that, as has been observed in other cancer types, RNA editing associates to EMT phenotype in PM.

Heterogeneous RNA editing and influence of ADAR2 on mesothelioma chemoresistance and the tumor microenvironment

We previously observed increased levels of adenosine-deaminase-acting-on-dsRNA (Adar)-dependent RNA editing during mesothelioma development in mice exposed to asbestos. The aim of this study was to characterize and assess the role of ADAR-dependent RNA editing in mesothelioma. We found that tumors and mesothelioma primary cultures have higher ADAR-mediated RNA editing compared to mesothelial cells. Unsupervised clustering of editing in different genomic regions revealed heterogeneity between tumor samples as well as mesothelioma primary cultures. ADAR2 expression levels are higher in BRCA1-associated protein 1 wild-type tumors, with corresponding changes in RNA editing in transcripts and 3'UTR. ADAR2 knockdown and rescue models indicated a role in cell proliferation, altered cell cycle, increased sensitivity to antifolate treatment, and type-1 interferon signaling upregulation, leading to changes in the microenvironment in vivo. Our data indicate that RNA editing contributes to mesothelioma heterogeneity and highlights an important role of ADAR2 not only in growth regulation in mesothelioma but also in chemotherapy response, in addition to regulating inflammatory response downstream of sensing nucleic acid structures.

Dysregulation of splicing variants and spliceosome components in breast cancer

The dysregulation of the splicing process has emerged as a novel hallmark of metabolic and tumor pathologies. In breast cancer (BCa), which represents the most diagnosed cancer type among women worldwide, the generation and/or dysregulation of several oncogenic splicing variants have been described. This is the case of the splicing variants of HER2, ER, BRCA1, or the recently identified by our group, In1-ghrelin and SST5TMD4, which exhibit oncogenic roles, increasing the malignancy, poor prognosis, and resistance to treatment of BCa. This altered expression of oncogenic splicing variants has been closely linked with the dysregulation of the elements belonging to the macromolecular machinery that controls the splicing process (spliceosome components and the associated splicing factors). In this review, we compile the current knowledge demonstrating the altered expression of splicing variants and spliceosomal components in BCa, showing the existence of a growing body of evidence supporting the close implication of the alteration in the splicing process in mammary tumorigenesis.

A complex epigenome-splicing crosstalk governs epithelial-to-mesenchymal transition in metastasis and brain development

Epithelial-to-mesenchymal transition (EMT) renders epithelial cells migratory properties. While epigenetic and splicing changes have been implicated in EMT, the mechanisms governing their crosstalk remain poorly understood. Here we discovered that a C2H2 zinc finger protein, ZNF827, is strongly induced during various contexts of EMT, including in brain development and breast cancer metastasis, and is required for the molecular and phenotypic changes underlying EMT in these processes. Mechanistically, ZNF827 mediated these responses by orchestrating a large-scale remodelling of the splicing landscape by recruiting HDAC1 for epigenetic modulation of distinct genomic loci, thereby slowing RNA polymerase II progression and altering the splicing of genes encoding key EMT regulators in cis. Our findings reveal an unprecedented complexity of crosstalk between epigenetic landscape and splicing programme in governing EMT and identify ZNF827 as a master regulator coupling these processes during EMT in brain development and breast cancer metastasis.

RBFOX2-regulated TEAD1 alternative splicing plays a pivotal role in Hippo-YAP signaling

Alternative pre-mRNA splicing is key to proteome diversity; however, the biological roles of alternative splicing (AS) in signaling pathways remain elusive. Here, we focus on TEA domain transcription factor 1 (TEAD1), a YAP binding factor in the Hippo signaling pathway. Public database analyses showed that expression of YAP-TEAD target genes negatively correlated with the expression of a TEAD1 isoform lacking exon 6 (TEAD1ΔE6) but did not correlate with overall TEAD1 expression. We confirmed that the transcriptional activity and oncogenic properties of the full-length TEAD1 isoform were greater than those of TEAD1ΔE6, with the difference in transcription related to YAP interaction. Furthermore, we showed that RNA-binding Fox-1 homolog 2 (RBFOX2) promoted the inclusion of TEAD1 exon 6 via binding to the conserved GCAUG element in the downstream intron. These results suggest a regulatory mechanism of RBFOX2-mediated TEAD1 AS and provide insight into AS-specific modulation of signaling pathways.

Investigation of the effects of overexpression of jumping translocation breakpoint (JTB) protein in MCF7 cells for potential use as a biomarker in breast cancer

Jumping translocation breakpoint (JTB) gene acts as a tumor suppressor or an oncogene in different malignancies, including breast cancer (BC), where it was reported as overexpressed. However, the molecular functions, biological processes and underlying mechanisms through which JTB protein causes increased cell growth, proliferation and invasion is still not fully deciphered. Our goal is to identify the functions of JTB protein by cellular proteomics approaches. MCF7 breast cancer cells were transfected with sense orientation of hJTB cDNA in HA, His and FLAG tagged CMV expression vector to overexpress hJTB and the expression levels were confirmed by Western blotting (WB). Proteins extracted from transfected cells were separated by SDS-PAGE and the in-gel digested peptides were analyzed by nano-liquid chromatography tandem mass spectrometry (nanoLC-MS/MS). By comparing the proteome of cells with upregulated conditions of JTB vs control and identifying the protein dysregulation patterns, we aim to understand the function of this protein and its contribution to tumorigenesis. Gene Set Enrichment Analysis (GSEA) algorithm was performed to investigate the biological processes and pathways that are associated with the JTB protein upregulation. The results demonstrated four significantly enriched gene sets from the following significantly upregulated pathways: mitotic spindle assembly, estrogen response late, epithelial-to-mesenchymal transition (EMT) and estrogen response early. JTB protein itself is involved in mitotic spindle pathway by its role in cell division/cytokinesis, and within estrogen response early and late pathways, contributing to discrimination between luminal and mesenchymal breast cancer. Thus, the overexpressed JTB condition was significantly associated with an increased expression of ACTNs, FLNA, FLNB, EZR, MYOF, COL3A1, COL11A1, HSPA1A, HSP90A, WDR, EPPK1, FASN and FOXA1 proteins related to deregulation of cytoskeletal organization and biogenesis, mitotic spindle organization, ECM remodeling, cellular response to estrogen, proliferation, migration, metastasis, increased lipid biogenesis, endocrine therapy resistance, antiapoptosis and discrimination between different breast cancer subtypes. Other upregulated proteins for overexpressed JTB condition are involved in multiple cellular functions and pathways that become dysregulated, such as tumor microenvironment (TME) acidification, the transmembrane transport pathways, glycolytic flux, iron metabolism and oxidative stress, metabolic reprogramming, nucleocytosolic mRNA transport, transcriptional activation, chromatin remodeling, modulation of cell death pathways, stress responsive pathways, and cancer drug resistance. The downregulated proteins for overexpressed JTB condition are involved in adaptive communication between external and internal environment of cells and maintenance between pro-apoptotic and anti-apoptotic signaling pathways, vesicle trafficking and secretion, DNA lesions repair and suppression of genes involved in tumor progression, proteostasis, redox state regulation, biosynthesis of macromolecules, lipolytic pathway, carbohydrate metabolism, dysregulation of ubiquitin-mediated degradation system, cancer cell immune escape, cell-to-cell and cell-to-ECM interactions, and cytoskeletal behaviour. There were no significantly enriched downregulated pathways.

Transcriptional and post-transcriptional control of epithelial-mesenchymal plasticity: why so many regulators?

The dynamic transition between epithelial-like and mesenchymal-like cell states has been a focus for extensive investigation for decades, reflective of the importance of Epithelial-Mesenchymal Transition (EMT) through development, in the adult, and the contributing role EMT has to pathologies including metastasis and fibrosis. Not surprisingly, regulation of the complex genetic networks that underlie EMT have been attributed to multiple transcription factors and microRNAs. What is surprising, however, are the sheer number of different regulators (hundreds of transcription factors and microRNAs) for which critical roles have been described. This review seeks not to collate these studies, but to provide a perspective on the fundamental question of whether it is really feasible that so many regulators play important roles and if so, what does this tell us about EMT and more generally, the genetic machinery that controls complex biological processes.

Transcriptome profiling of ADAR1 targets in triple-negative breast cancer cells reveals mechanisms for regulating growth and invasion

ADARs catalyze Adenosine-to-Inosine (A-to-I) editing of double-stranded RNA and regulate global gene expression output through interactions with RNA and other proteins. ADARs play important roles in development and disease, and previous work has shown that ADAR1 is oncogenic in a growing list of cancer types. Here we show that ADAR1 is a critical gene for triple-negative breast cancer cells, as ADAR1 loss results in reduced growth (viability and cell cycle progression), invasion, and mammosphere formation. Whole transcriptome sequencing analyses demonstrate that ADAR1 regulates both coding and non-coding targets by altering gene expression level, A-to-I editing, and splicing. We determine that a recoding edit in filamin B (FLNB chr3:58156064) reduces the tumor suppressive activities of the protein to promote growth and invasion. We also show that several tumor suppressor microRNAs are upregulated upon ADAR1 loss and suppress cell cycle progression and invasion. Implications: This work describes several novel mechanisms of ADAR1-mediated oncogenesis in triple-negative breast cancer, providing support to strategies targeting ADAR1 in this aggressive cancer type that has few treatment options.

Characterization of distinct circular RNA signatures in solid tumors

Background

Circular RNAs (circRNAs) are differentially expressed between normal and cancerous tissues, contributing to tumor initiation and progression. However, comprehensive landscape of dysregulated circRNAs across cancer types remains unclear.

Methods

In this study, we conducted Ribo-Zero transcriptome sequencing on tumor tissues and their adjacent normal samples including glioblastoma, esophageal squamous cell carcinoma, lung adenocarcinoma, thyroid cancer, colorectal cancer, gastric cancer and hepatocellular carcinoma. CIRCexplorer2 was employed to identify circRNAs and dysregulated circRNAs and genes were determined by DESeq2 package. The expression of hsa_circ_0072309 (circLIFR) was measured by reverse transcription and quantitative real-time PCR, and its effect on cell migration was examined by Transwell and wound healing assays. The role of circLIFR in tumor metastasis was evaluated via mouse models of tail-vein injection and spleen injection for lung and liver metastasis, respectively.

Results

Distinct circRNA expression signatures were identified among seven types of solid tumors, and the dysregulated circRNAs exhibited cancer-specific expression or shared common expression signatures across cancers. Bioinformatics analyses indicated that aberrant expression of host genes and/or RNA-binding proteins contributed to circRNA dysregulation in cancer. Finally, circLIFR was experimentally validated to be downregulated in six solid tumors and to significantly inhibit cell migration in vitro and tumor metastasis in vivo.

Conclusions

Our results provide a comprehensive landscape of differentially expressed circRNAs in solid tumors and highlight that circRNAs are extensively involved in cancer pathogenesis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12943-022-01546-4.

Comprehensive profiling of mRNA splicing indicates that GC content signals altered cassette exon inclusion in Ewing sarcoma

Ewing sarcoma (EwS) is a small round blue cell tumor and is the second most frequent pediatric bone cancer. 85% of EwS tumors express the fusion oncoprotein EWS-FLI1, the product of a t(11;22) reciprocal translocation. Prior work has indicated that transcription regulation alone does not fully describe the oncogenic capacity of EWS-FLI1, nor does it provide an effective means to stratify patient tumors. Research using EwS cell lines and patient samples has suggested that EWS-FLI1 also disrupts mRNA biogenesis. In this work we both describe the underlying characteristics of mRNA that are aberrantly spliced in EwS tumor samples as well as catalogue mRNA splicing events across other pediatric tumor types. Here, we also use short- and long-read sequencing to identify cis-factors that contribute to splicing profiles we observe in Ewing sarcoma. Our analysis suggests that GC content upstream of cassette exons is a defining factor of mRNA splicing in EwS. We also describe specific splicing events that discriminate EwS tumor samples from the assumed cell of origin, human mesenchymal stem cells derived from bone marrow (hMSC-BM). Finally, we identify specific splicing factors PCBP2, RBMX, and SRSF9 by motif enrichment and confirm findings from tumor samples in EwS cell lines.

The Quaking RNA-binding proteins as regulators of cell differentiation

The RNA-binding protein Quaking (QKI) has emerged as a potent regulator of cellular differentiation in developmental and pathological processes. The QKI gene is itself alternatively spliced to produce three major isoforms, QKI-5, QKI-6, and QKI-7, that possess very distinct functions. Here, we highlight roles of the different QKI isoforms in neuronal, vascular, muscle, and monocyte cell differentiation, and during epithelial-mesenchymal transition in cancer progression. QKI isoforms control cell differentiation through regulating alternative splicing, mRNA stability and translation, with activities in gene transcription now also becoming evident. These diverse functions of the QKI isoforms contribute to their broad influences on RNA metabolism and cellular differentiation. This article is categorized under: RNA Interactions with Proteins and Other Molecules &gt; Protein-RNA Interactions: Functional Implications RNA Processing &gt; Splicing Regulation/Alternative Splicing RNA in Disease and Development &gt; RNA in Development.

Regulation of Alternative Splicing during Epithelial-Mesenchymal Transition

Alternative splicing is an essential mechanism of gene regulation, giving rise to remarkable protein diversity in higher eukaryotes. Epithelial-mesenchymal transition (EMT) is a developmental process that plays an essential role in metazoan embryogenesis. Recent studies have revealed that alternative splicing serves as a fundamental layer of regulation that governs cells to undergo EMT. In this review, we summarize recent findings on the functional impact of alternative splicing in EMT and EMT-associated activities. We then discuss the regulatory mechanisms that control alternative splicing changes during EMT.

A Multi-Omics Network of a Seven-Gene Prognostic Signature for Non-Small Cell Lung Cancer

There is an unmet clinical need to identify patients with early-stage non-small cell lung cancer (NSCLC) who are likely to develop recurrence and to predict their therapeutic responses. Our previous study developed a qRT-PCR-based seven-gene microfluidic assay to predict the recurrence risk and the clinical benefits of chemotherapy. This study showed it was feasible to apply this seven-gene panel in RNA sequencing profiles of The Cancer Genome Atlas (TCGA) NSCLC patients (n = 923) in randomly partitioned feasibility-training and validation sets (p < 0.05, Kaplan-Meier analysis). Using Boolean implication networks, DNA copy number variation-mediated transcriptional regulatory network of the seven-gene signature was identified in multiple NSCLC cohorts (n = 371). The multi-omics network genes, including PD-L1, were significantly correlated with immune infiltration and drug response to 10 commonly used drugs for treating NSCLC. ZNF71 protein expression was positively correlated with epithelial markers and was negatively correlated with mesenchymal markers in NSCLC cell lines in Western blots. PI3K was identified as a relevant pathway of proliferation networks involving ZNF71 and its isoforms formulated with CRISPR-Cas9 and RNA interference (RNAi) profiles. Based on the gene expression of the multi-omics network, repositioning drugs were identified for NSCLC treatment.

Emerging role of m6A methylation modification in ovarian cancer

m6A (N6-methyladenosine) methylation, a well-known modification in tumour epigenetics, dynamically and reversibly fine tunes the entire process of RNA metabolism. Aberrant levels of m6A and its regulators, which can predict the survival and outcomes of cancer patients, are involved in tumorigenesis, metastasis and resistance. Ovarian cancer (OC) ranks first among gynaecological tumours in the causes of death. At first diagnosis, patients with OC are usually at advanced stages owing to a lack of early biomarkers and effective targets. After treatment, patients with OC often develop drug resistance. This article reviews the recent experimental advances in understanding the role of m6A modification in OC, raising the possibility to treat m6A modification and its regulators as promising diagnostic markers and therapeutic targets for OC.

Alternative RNA splicing in stem cells and cancer stem cells: Importance of transcript-based expression analysis

Alternative ribonucleic acid (RNA) splicing can lead to the assembly of different protein isoforms with distinctive functions. The outcome of alternative splicing (AS) can result in a complete loss of function or the acquisition of new functions. There is a gap in knowledge of abnormal RNA splice variants promoting cancer stem cells (CSCs), and their prospective contribution in cancer progression. AS directly regulates the self-renewal features of stem cells (SCs) and stem-like cancer cells. Notably, octamer-binding transcription factor 4A spliced variant of octamer-binding transcription factor 4 contributes to maintaining stemness properties in both SCs and CSCs. The epithelial to mesenchymal transition pathway regulates the AS events in CSCs to maintain stemness. The alternative spliced variants of CSCs markers, including cluster of differentiation 44, aldehyde dehydrogenase, and doublecortin-like kinase, α6β1 integrin, have pivotal roles in increasing self-renewal properties and maintaining the pluripotency of CSCs. Various splicing analysis tools are considered in this study. LeafCutter software can be considered as the best tool for differential splicing analysis and identification of the type of splicing events. Additionally, LeafCutter can be used for efficient mapping splicing quantitative trait loci. Altogether, the accumulating evidence re-enforces the fact that gene and protein expression need to be investigated in parallel with alternative splice variants.

The RNA binding protein Quaking represses splicing of the Fibronectin EDA exon and downregulates the interferon response

Quaking (QKI) controls RNA metabolism in many biological processes including innate immunity, where its roles remain incompletely understood. To illuminate these roles, we performed genome scale transcriptome profiling in QKI knockout cells with or without poly(I:C) transfection, a double-stranded RNA analog that mimics viral infection. Analysis of RNA-sequencing data shows that QKI knockout upregulates genes induced by interferons, suggesting that QKI is an immune suppressor. Furthermore, differential splicing analysis shows that QKI primarily controls cassette exons, and among these events, we noted that QKI silences splicing of the extra domain A (EDA) exon in fibronectin (FN1) transcripts. QKI knockout results in elevated production and secretion of FN1-EDA protein, which is a known activator of interferons. Consistent with an upregulation of the interferon response in QKI knockout cells, our results show reduced production of dengue virus-2 and Japanese encephalitis virus in these cells. In conclusion, we demonstrate that QKI downregulates the interferon system and attenuates the antiviral state.

Therapeutically Targeting Cancers That Overexpress FOXC1: A Transcriptional Driver of Cell Plasticity, Partial EMT, and Cancer Metastasis

Metastasis accounts for more than 90% of cancer related mortality, thus the most pressing need in the field of oncology today is the ability to accurately predict future onset of metastatic disease, ideally at the time of initial diagnosis. As opposed to current practice, what would be desirable is that prognostic, biomarker-based detection of metastatic propensity and heightened risk of cancer recurrence be performed long before overt metastasis has set in. Without such timely information it will be impossible to formulate a rational therapeutic treatment plan to favorably alter the trajectory of disease progression. In order to help inform rational selection of targeted therapeutics, any recurrence/metastasis risk prediction strategy must occur with the paired identification of novel prognostic biomarkers and their underlying molecular regulatory mechanisms that help drive cancer recurrence/metastasis (i.e. recurrence biomarkers). Traditional clinical factors alone (such as TNM staging criteria) are no longer adequately prognostic for this purpose in the current molecular era. FOXC1 is a pivotal transcription factor that has been functionally implicated to drive cancer metastasis and has been demonstrated to be an independent predictor of heightened metastatic risk, at the time of initial diagnosis. In this review, we present our viewpoints on the master regulatory role that FOXC1 plays in mediating cancer stem cell traits that include cellular plasticity, partial EMT, treatment resistance, cancer invasion and cancer migration during cancer progression and metastasis. We also highlight potential therapeutic strategies to target cancers that are, or have evolved to become, “transcriptionally addicted” to FOXC1. The potential role of FOXC1 expression status in predicting the efficacy of these identified therapeutic approaches merits evaluation in clinical trials.

Single-cell-resolved differentiation of human induced pluripotent stem cells into pancreatic duct-like organoids on a microwell chip

Creating in vitro models of diseases of the pancreatic ductal compartment requires a comprehensive understanding of the developmental trajectories of pancreas-specific cell types. Here we report the single-cell characterization of the differentiation of pancreatic duct-like organoids (PDLOs) from human induced pluripotent stem cells (hiPSCs) on a microwell chip that facilitates the uniform aggregation and chemical induction of hiPSC-derived pancreatic progenitors. Using time-resolved single-cell transcriptional profiling and immunofluorescence imaging of the forming PDLOs, we identified differentiation routes from pancreatic progenitors through ductal intermediates to two types of mature duct-like cells and a few non-ductal cell types. PDLO subpopulations expressed either mucins or the cystic fibrosis transmembrane conductance regulator, and resembled human adult duct cells. We also used the chip to uncover ductal markers relevant to pancreatic carcinogenesis, and to establish PDLO co-cultures with stellate cells, which allowed for the study of epithelial-mesenchymal signalling. The PDLO microsystem could be used to establish patient-specific pancreatic duct models.

Matrix-screening reveals a vast potential for direct protein-protein interactions among RNA binding proteins

RNA-binding proteins (RBPs) are crucial factors of post-transcriptional gene regulation and their modes of action are intensely investigated. At the center of attention are RNA motifs that guide where RBPs bind. However, sequence motifs are often poor predictors of RBP-RNA interactions in vivo. It is hence believed that many RBPs recognize RNAs as complexes, to increase specificity and regulatory possibilities. To probe the potential for complex formation among RBPs, we assembled a library of 978 mammalian RBPs and used rec-Y2H matrix screening to detect direct interactions between RBPs, sampling > 600 K interactions. We discovered 1994 new interactions and demonstrate that interacting RBPs bind RNAs adjacently in vivo. We further find that the mRNA binding region and motif preferences of RBPs deviate, depending on their adjacently binding interaction partners. Finally, we reveal novel RBP interaction networks among major RNA processing steps and show that splicing impairing RBP mutations observed in cancer rewire spliceosomal interaction networks. The dataset we provide will be a valuable resource for understanding the combinatorial interactions of RBPs with RNAs and the resulting regulatory outcomes.

Systematic Profiling of mRNA Splicing Reveals the Prognostic Predictor and Potential Therapeutic Target for Glioblastoma Multiforme

Despite many changes in alternative splicing events (ASEs) are frequently involved in various cancers, prognosis-related ASEs and drug treatment targets in glioblastoma multiforme (GBM) have not been well explored. ASEs participate in many biological behaviors in the initiation and progression of tumors, the aberrant ASE has been considered another hallmark of cancer, and the systematic study of alternative splicing may provide potential biomarkers for malignancies. In this study, we carried out a systematic analysis to characterize the ASE signatures in GBM cohort. Through comparing GBM tissues and nontumor tissues, a total of 48,191 differently expressed ASEs from 10,727 genes were obtained, and these aberrant ASEs play an important role in the oncogenic process. Then, we identified 514 ASEs independently associated with patient survival in GBM by univariate and multivariate Cox regression, including exon skip in CD3D, alternate acceptor site in POLD2, and exon skip in DCN. Those prognostic models built on ASEs of each splice type can accurately predict the outcome of GBM patients, and values for the area under curve were 0.97 in the predictive model based on alternate acceptor site. In addition, the splicing-regulatory network revealed an interesting correlation between survival-associated splicing factors and prognostic ASE corresponding genes. Moreover, these three hub splicing factors in splicing regulation network are the potential targets of some drugs. In conclusion, a systematic analysis of ASE signatures in GBM could serve as an indicator for identifying novel prognostic biomarkers and guiding clinical treatment.

Microexon alternative splicing of small GTPase regulators: Implication in central nervous system diseases

Microexons are small sized (≤51 bp) exons which undergo extensive alternative splicing in neurons, microglia, embryonic stem cells, and cancer cells, giving rise to cell type specific protein isoforms. Due to their small sizes, microexons provide a unique challenge for the splicing machinery. They frequently lack exon splicer enhancers/repressors and require specialized neighboring trans-regulatory and cis-regulatory elements bound by RNA binding proteins (RBPs) for their inclusion. The functional consequences of including microexons within mRNAs have been extensively documented in the central nervous system (CNS) and aberrations in their inclusion have been observed to lead to abnormal processes. Despite the increasing evidence for microexons impacting cellular physiology within CNS, mechanistic details illustrating their functional importance in diseases of the CNS is still limited. In this review, we discuss the unique characteristics of microexons, and how RBPs participate in regulating their inclusion and exclusion during splicing. We consider recent findings of microexon alternative splicing and their implication for regulating the function of small GTPases in the context of the microglia, and we extrapolate these findings to what is known in neurons. We further discuss the emerging evidence for dysregulation of the Rho GTPase pathway in CNS diseases and the consequences contributed by the mis-splicing of microexons. This article is categorized under: RNA Processing > Splicing Mechanisms RNA Processing > Splicing Regulation/Alternative Splicing RNA in Disease and Development > RNA in Disease.

The Emerging Roles of the RNA Binding Protein QKI in Cardiovascular Development and Function

RNA binding proteins (RBPs) have a broad biological and physiological function and are critical in regulating pre-mRNA posttranscriptional processing, intracellular migration, and mRNA stability. QKI, also known as Quaking, is a member of the signal transduction and activation of RNA (STAR) family, which also belongs to the heterogeneous nuclear ribonucleoprotein K- (hnRNP K-) homology domain protein family. There are three major alternatively spliced isoforms, QKI-5, QKI-6, and QKI-7, differing in carboxy-terminal domains. They share a common RNA binding property, but each isoform can regulate pre-mRNA splicing, transportation or stability differently in a unique cell type-specific manner. Previously, QKI has been known for its important role in contributing to neurological disorders. A series of recent work has further demonstrated that QKI has important roles in much broader biological systems, such as cardiovascular development, monocyte to macrophage differentiation, bone metabolism, and cancer progression. In this mini-review, we will focus on discussing the emerging roles of QKI in regulating cardiac and vascular development and function and its potential link to cardiovascular pathophysiology.

A novel class of inhibitors that target SRSF10 and promote p53-mediated cytotoxicity on human colorectal cancer cells

The elevated expression of the splicing regulator SRSF10 in metastatic colorectal cancer (CRC) stimulates the production of the pro-tumorigenic BCLAF1-L splice variant. We discovered a group of small molecules with an aminothiazole carboxamide core (GPS167, GPS192 and others) that decrease production of BCLAF1-L. While additional alternative splicing events regulated by SRSF10 are affected by GPS167/192 in HCT116 cells (e.g. in MDM4, WTAP, SLK1 and CLK1), other events are shifted in a SRSF10-independent manner (e.g. in MDM2, NAB2 and TRA2A). GPS167/192 increased the interaction of SRSF10 with the CLK1 and CLK4 kinases, leading us to show that GPS167/192 can inhibit CLK kinases preferentially impacting the activity of SRSF10. Notably, GPS167 impairs the growth of CRC cell lines and organoids, inhibits anchorage-independent colony formation, cell migration, and promotes cytoxicity in a manner that requires SRSF10 and p53. In contrast, GPS167 only minimally affects normal colonocytes and normal colorectal organoids. Thus, GPS167 reprograms the tumorigenic activity of SRSF10 in CRC cells to elicit p53-dependent apoptosis.

Knockdown of Splicing Complex Protein PCBP2 Reduces Extravillous Trophoblast Differentiation Through Transcript Switching

Mutations in the LINC-HELLP non-coding RNA (HELLPAR) have been associated with familial forms of the pregnancy-specific HELLP syndrome. These mutations negatively affect extravillous trophoblast (EVT) differentiation from a proliferative to an invasive state and disturb the binding of RNA splicing complex proteins PCBP1, PCBP2, and YBX1 to LINC-HELLP. In this study, by using both in vitro and ex vivo experiments, we investigate if these proteins are involved in the regulation of EVT invasion during placentation. Additionally, we study if this regulation is due to alternative mRNA splicing. HTR-8/SVneo extravillous trophoblasts and human first trimester placental explants were used to investigate the effect of siRNA-mediated downregulation of PCBP1, PCBP2, and YBX1 genes on the differentiation of EVTs. Transwell invasion assays and proliferation assays indicated that upon knockdown of PCBP2 and, to a lesser extent, YBX1 and PCBP1, EVTs fail to differentiate toward an invasive phenotype. The same pattern was observed in placental explants where PCBP2 knockdown led to approximately 80% reduction in the number of explants showing any EVT outgrowth. Of the ones that still did show EVT outgrowth, the percentage of proliferating EVTs was significantly higher compared to explants transfected with non-targeting control siRNAs. To further investigate this effect of PCBP2 silencing on EVTs, we performed whole transcriptome sequencing (RNA-seq) on HTR-8/SVneo cells after PCBP2 knockdown. PCBP2 knockdown was found to have minimal effect on mRNA expression levels. In contrast, PCBP2 silencing led to a switch in splicing for a large number of genes with predominant functions in cellular assembly and organization, cellular function and maintenance, and cellular growth and proliferation and the cell cycle. EVTs, upon differentiation, alter their function to be able to invade the decidua of the mother by changing their cellular assembly and their proliferative activity by exiting the cell cycle. PCBP2 appears to be a paramount regulator of these differentiation mechanisms, where its disturbed binding to LINC-HELLP could contribute to dysfunctional placental development as seen in the HELLP syndrome.

A cell-to-patient machine learning transfer approach uncovers novel basal-like breast cancer prognostic markers amongst alternative splice variants

BACKGROUND

Breast cancer is amongst the 10 first causes of death in women worldwide. Around 20% of patients are misdiagnosed leading to early metastasis, resistance to treatment and relapse. Many clinical and gene expression profiles have been successfully used to classify breast tumours into 5 major types with different prognosis and sensitivity to specific treatments. Unfortunately, these profiles have failed to subclassify breast tumours into more subtypes to improve diagnostics and survival rate. Alternative splicing is emerging as a new source of highly specific biomarkers to classify tumours in different grades. Taking advantage of extensive public transcriptomics datasets in breast cancer cell lines (CCLE) and breast cancer tumours (TCGA), we have addressed the capacity of alternative splice variants to subclassify highly aggressive breast cancers.

RESULTS

Transcriptomics analysis of alternative splicing events between luminal, basal A and basal B breast cancer cell lines identified a unique splicing signature for a subtype of tumours, the basal B, whose classification is not in use in the clinic yet. Basal B cell lines, in contrast with luminal and basal A, are highly metastatic and express epithelial-to-mesenchymal (EMT) markers, which are hallmarks of cell invasion and resistance to drugs. By developing a semi-supervised machine learning approach, we transferred the molecular knowledge gained from these cell lines into patients to subclassify basal-like triple negative tumours into basal A- and basal B-like categories. Changes in splicing of 25 alternative exons, intimately related to EMT and cell invasion such as ENAH, CD44 and CTNND1, were sufficient to identify the basal-like patients with the worst prognosis. Moreover, patients expressing this basal B-specific splicing signature also expressed newly identified biomarkers of metastasis-initiating cells, like CD36, supporting a more invasive phenotype for this basal B-like breast cancer subtype.

CONCLUSIONS

Using a novel machine learning approach, we have identified an EMT-related splicing signature capable of subclassifying the most aggressive type of breast cancer, which are basal-like triple negative tumours. This proof-of-concept demonstrates that the biological knowledge acquired from cell lines can be transferred to patients data for further clinical investigation. More studies, particularly in 3D culture and organoids, will increase the accuracy of this transfer of knowledge, which will open new perspectives into the development of novel therapeutic strategies and the further identification of specific biomarkers for drug resistance and cancer relapse.

Molecular Analysis of ZNF71 KRAB in Non-Small-Cell Lung Cancer

Our previous study found that zinc finger protein 71 (ZNF71) mRNA expression was associated with chemosensitivity and its protein expression was prognostic of non-small-cell lung cancer (NSCLC). The Krüppel associated box (KRAB) transcriptional repression domain is commonly present in human zinc finger proteins, which are linked to imprinting, silencing of repetitive elements, proliferation, apoptosis, and cancer. This study revealed that ZNF71 KRAB had a significantly higher expression than the ZNF71 KRAB-less isoform in NSCLC tumors (n = 197) and cell lines (n = 117). Patients with higher ZNF71 KRAB expression had a significantly worse survival outcome than patients with lower ZNF71 KRAB expression (log-rank p = 0.04; hazard ratio (HR): 1.686 [1.026, 2.771]), whereas ZNF71 overall and KRAB-less expression levels were not prognostic in the same patient cohort. ZNF71 KRAB expression was associated with epithelial-to-mesenchymal transition (EMT) in both patient tumors and cell lines. ZNF71 KRAB was overexpressed in NSCLC cell lines resistant to docetaxel and paclitaxel treatment compared to chemo-sensitive cell lines, consistent with its association with poor prognosis in patients. Therefore, ZNF71 KRAB isoform is a more effective prognostic factor than ZNF71 overall and KRAB-less expression for NSCLC. Functional analysis using CRISPR-Cas9 and RNA interference (RNAi) screening data indicated that a knockdown/knockout of ZNF71 did not significantly affect NSCLC cell proliferation in vitro.

Quaking 5 suppresses TGF-β-induced EMT and cell invasion in lung adenocarcinoma

Quaking (QKI) proteins belong to the signal transduction and activation of RNA (STAR) family of RNA-binding proteins that have multiple functions in RNA biology. Here, we show that QKI-5 is dramatically decreased in metastatic lung adenocarcinoma (LUAD). QKI-5 overexpression inhibits TGF-β-induced epithelial-mesenchymal transition (EMT) and invasion, whereas QKI-5 knockdown has the opposite effect. QKI-5 overexpression and silencing suppresses and promotes TGF-β-stimulated metastasis in vivo, respectively. QKI-5 inhibits TGF-β-induced EMT and invasion in a TGFβR1-dependent manner. KLF6 knockdown increases TGFβR1 expression and promotes TGF-β-induced EMT, which is partly abrogated by QKI-5 overexpression. Mechanistically, QKI-5 directly interacts with the TGFβR1 3' UTR and causes post-transcriptional degradation of TGFβR1 mRNA, thereby inhibiting TGF-β-induced SMAD3 phosphorylation and TGF-β/SMAD signaling. QKI-5 is positively regulated by KLF6 at the transcriptional level. In LUAD tissues, KLF6 is lowly expressed and positively correlated with QKI-5 expression, while TGFβR1 expression is up-regulated and inversely correlated with QKI-5 expression. We reveal a novel mechanism by which KLF6 transcriptionally regulates QKI-5 and suggest that targeting the KLF6/QKI-5/TGFβR1 axis is a promising targeting strategy for metastatic LUAD.

QKI is a critical pre-mRNA alternative splicing regulator of cardiac myofibrillogenesis and contractile function

The RNA-binding protein QKI belongs to the hnRNP K-homology domain protein family, a well-known regulator of pre-mRNA alternative splicing and is associated with several neurodevelopmental disorders. Qki is found highly expressed in developing and adult hearts. By employing the human embryonic stem cell (hESC) to cardiomyocyte differentiation system and generating QKI-deficient hESCs (hESCs-QKI^del) using CRISPR/Cas9 gene editing technology, we analyze the physiological role of QKI in cardiomyocyte differentiation, maturation, and contractile function. hESCs-QKI^del largely maintain normal pluripotency and normal differentiation potential for the generation of early cardiogenic progenitors, but they fail to transition into functional cardiomyocytes. In this work, by using a series of transcriptomic, cell and biochemical analyses, and the Qki-deficient mouse model, we demonstrate that QKI is indispensable to cardiac sarcomerogenesis and cardiac function through its regulation of alternative splicing in genes involved in Z-disc formation and contractile physiology, suggesting that QKI is associated with the pathogenesis of certain forms of cardiomyopathies.

Splicing factors: Insights into their regulatory network in alternative splicing in cancer

More than 95% of all human genes are alternatively spliced after transcription, which enriches the diversity of proteins and regulates transcript and/or protein levels. The splicing isoforms produced from the same gene can manifest distinctly, even exerting opposite effects. Mounting evidence indicates that the alternative splicing (AS) mechanism is ubiquitous in various cancers and drives the generation and maintenance of various hallmarks of cancer, such as enhanced proliferation, inhibited apoptosis, invasion and metastasis, and angiogenesis. Splicing factors (SFs) play pivotal roles in the recognition of splice sites and the assembly of spliceosomes during AS. In this review, we mainly discuss the similarities and differences of SF domains, the details of SF function in AS, the effect of SF-driven pathological AS on different hallmarks of cancer, and the main drivers of SF expression level and subcellular localization. In addition, we briefly introduce the application prospects of targeted therapeutic strategies, including small-molecule inhibitors, siRNAs and splice-switching oligonucleotides (SSOs), from three perspectives (drivers, SFs and pathological AS). Finally, we share our insights into the potential direction of research on SF-centric AS-related regulatory networks.