A large-scale analysis of alternative splicing reveals a key role of QKI in lung cancer

Qki5 safeguards spinal motor neuron function by defining the motor neuron-specific transcriptome via pre-mRNA processing

Many RNA-binding proteins (RBPs) are linked to the dysregulation of RNA metabolism in motor neuron diseases (MNDs). However, the molecular mechanisms underlying MN vulnerability have yet to be elucidated. Here, we found that such an RBP, Quaking5 (Qki5), contributes to formation of the MN-specific transcriptome profile, termed "MN-ness," through the posttranscriptional network and maintenance of the mature MNs. Immunohistochemical analysis and single-cell RNA sequencing (scRNA-seq) revealed that Qki5 is predominantly expressed in MNs, but not in other neuronal populations of the spinal cord. Furthermore, comprehensive RNA sequencing (RNA-seq) analyses revealed that Qki5-dependent RNA regulation plays a pivotal role in generating the MN-specific transcriptome through pre-messenger ribonucleic acid (mRNA) splicing for the synapse-related molecules and c-Jun N-terminal kinase/stress-activated protein kinase (JNK/SAPK) signaling pathways. Indeed, MN-specific ablation of the Qki5 caused neurodegeneration in postnatal mice and loss of Qki5 function resulted in the aberrant activation of stress-responsive JNK/SAPK pathway both in vitro and in vivo. These data suggested that Qki5 plays a crucial biological role in RNA regulation and safeguarding of MNs and might be associated with pathogenesis of MNDs.

Clinicopathologic Significance of Quaking Expression in Hepatocellular Carcinoma

BACKGROUND/AIM

The RNA binding protein quaking (QKI) is associated with the development and progression of tumor suppressors in various cancers. However, the clinical implications of QKI expression have not yet been fully elucidated. In this study, we aimed to investigate the clinicopathological and prognostic significance of QKI expression in hepatocellular carcinoma (HCC).

MATERIALS AND METHODS

We performed QKI, Zinc finger E-box-binding homeobox 1 (ZEB1), E-cadherin, and glutathione peroxidase 4 (GPX4) immunohistochemical staining on 166 HCC patient tissue samples. X-tile bioinformatics software was used to set the cut-off value for high QKI expression. Correlations between QKI expression and various clinicopathological parameters were assessed.

RESULTS

The best cut-off value for high QKI expression was 12.5. High QKI expression was observed in 28 of 166 patients (16.9%) and was an independent prognostic factor for inferior recurrence-free survival (RFS). In addition, high ZEB1 and GPX4 expression correlated with high QKI expression, but not with the loss of E-cadherin expression.

CONCLUSION

High QKI expression was identified in HCCs and associated with poor RFS. QKI might be a prognostic biomarker of HCCs associated with epithelial-to-mesenchymal transition and a potential candidate therapeutic target.

Development and validation of AI/ML derived splice-switching oligonucleotides

Splice-switching oligonucleotides (SSOs) are antisense compounds that act directly on pre-mRNA to modulate alternative splicing (AS). This study demonstrates the value that artificial intelligence/machine learning (AI/ML) provides for the identification of functional, verifiable, and therapeutic SSOs. We trained XGboost tree models using splicing factor (SF) pre-mRNA binding profiles and spliceosome assembly information to identify modulatory SSO binding sites on pre-mRNA. Using Shapley and out-of-bag analyses we also predicted the identity of specific SFs whose binding to pre-mRNA is blocked by SSOs. This step adds considerable transparency to AI/ML-driven drug discovery and informs biological insights useful in further validation steps. We applied this approach to previously established functional SSOs to retrospectively identify the SFs likely to regulate those events. We then took a prospective validation approach using a novel target in triple negative breast cancer (TNBC), NEDD4L exon 13 (NEDD4Le13). Targeting NEDD4Le13 with an AI/ML-designed SSO decreased the proliferative and migratory behavior of TNBC cells via downregulation of the TGFβ pathway. Overall, this study illustrates the ability of AI/ML to extract actionable insights from RNA-seq data.

Alternative splicing and related RNA binding proteins in human health and disease

Alternative splicing (AS) serves as a pivotal mechanism in transcriptional regulation, engendering transcript diversity, and modifications in protein structure and functionality. Across varying tissues, developmental stages, or under specific conditions, AS gives rise to distinct splice isoforms. This implies that these isoforms possess unique temporal and spatial roles, thereby associating AS with standard biological activities and diseases. Among these, AS-related RNA-binding proteins (RBPs) play an instrumental role in regulating alternative splicing events. Under physiological conditions, the diversity of proteins mediated by AS influences the structure, function, interaction, and localization of proteins, thereby participating in the differentiation and development of an array of tissues and organs. Under pathological conditions, alterations in AS are linked with various diseases, particularly cancer. These changes can lead to modifications in gene splicing patterns, culminating in changes or loss of protein functionality. For instance, in cancer, abnormalities in AS and RBPs may result in aberrant expression of cancer-associated genes, thereby promoting the onset and progression of tumors. AS and RBPs are also associated with numerous neurodegenerative diseases and autoimmune diseases. Consequently, the study of AS across different tissues holds significant value. This review provides a detailed account of the recent advancements in the study of alternative splicing and AS-related RNA-binding proteins in tissue development and diseases, which aids in deepening the understanding of gene expression complexity and offers new insights and methodologies for precision medicine.

Involvement of circRNA Regulators MBNL1 and QKI in the Progression of Esophageal Squamous Cell Carcinoma

OBJECTIVES

To investigate the role of circRNA regulators MBNL1 and QKI in the progression of esophageal squamous cell carcinoma.

BACKGROUND

MBNL1 and QKI are pivotal regulators of pre-mRNA alternative splicing, crucial for controlling circRNA production - an emerging biomarker and functional regulator of tumor progression. Despite their recognized roles, their involvement in ESCC progression remains unexplored.

METHODS

The expression levels of MBNL1 and QKI were examined in 28 tissue pairs from ESCC and adjacent normal tissues using data from the GEO database. Additionally, a total of 151 ESCC tissue samples, from stage T1 to T4, consisting of 13, 43, 87, and 8 cases per stage, respectively, were utilized for immunohistochemical (IHC) analysis. RNA sequencing was utilized to examine the expression profiles of circRNAs, lncRNAs, and mRNAs across 3 normal tissues, 3 ESCC tissues, and 3 pairs of KYSE150 cells in both wildtype (WT) and those with MBNL1 or QKI knockouts. Transwell, colony formation, and subcutaneous tumorigenesis assays assessed the impact of MBNL1 or QKI knockout on ESCC cell migration, invasion, and proliferation.

RESULTS

ESCC onset significantly altered MBNL1 and QKI expression levels, influencing diverse RNA species. Elevated MBNL1 or QKI expression correlated with patient age or tumor invasion depth, respectively. MBNL1 or QKI knockout markedly enhanced cancer cell migration, invasion, proliferation, and tumor growth. Moreover, the absence of either MBNL1 or QKI modulated the expression profiles of multiple circRNAs, causing extensive downstream alterations in the expression of numerous lncRNAs and mRNAs. While the functions of circRNA and lncRNA among the top 20 differentially expressed genes remain unclear, mRNAs like SLCO4C1, TMPRSS15, and MAGEB2 have reported associations with tumor progression.

CONCLUSIONS

This study underscores the tumor-suppressive roles of MBNL1 and QKI in ESCC, proposing them as potential biomarkers and therapeutic targets for ESCC diagnosis and treatment.

Expression profiling of luminal B breast tumor in Indian women

PURPOSE

In this study, we aimed at profiling of luminal B breast cancer specific gene expression pattern in Indian women using mRNA-seq and validation based on TCGA expression data.

METHODS

RNA isolated from luminal B tumor and adjacent normal tissues was used for library construction and sequencing. Reference-based assemblies of these reads were used for differential gene expression analysis using DeSeq2. The DEGs were evaluated using TCGA expression data. Kaplan-Meier survival method was used to evaluate association between genes showing luminal B specific differential expression pattern and breast cancer prognosis and statistical significance was assessed using log-rank test. Alternate splicing analysis was done using rmats.

RESULTS

Differential expression analysis identified 2371 differentially expressed genes (DEGs) in luminal B breast tumors in comparison with adjacent normal tissues of Indian Women. Of them, 1692 DEGs were validated using TCGA luminal B paired samples. Integration of this data with the DEGs obtained by comparative analysis of unpaired luminal B with luminal A unpaired samples from TCGA resulted in 291 DEGs showing luminal B specific expression pattern. Further, 26 genes of prognostic value were identified. Differential splicing analysis between luminal B tumors and adjacent normal tissues in our cohort led to the identification of 687 genes showing significant differential alternate splicing events.

CONCLUSION

This study profiled gene expression pattern of luminal B tumors of Indian women and identified 26 key genes of prognostic value for luminal B breast cancer. This study also profiled differential alternate splicing and identified important alternate splicing events in luminal B breast cancer.

RNA‑binding protein quaking 5 inhibits the progression of non‑small cell lung cancer by upregulating netrin‑4 expression

It was recently reported that netrin‑4 (Ntn‑4), a component of the extracellular matrix, when downregulated, is involved in the progression of several types of cancer, including breast cancer, colorectal tumours, neuroblastoma and gastric cancer. In the present study, the level of Ntn‑4 was examined in a public non‑small cell lung cancer (NSCLC) dataset from the Netherlands Cancer Institute. This analysis revealed that the mRNA expression level of Ntn‑4 was lower in the samples of patients with NSCLC compared with that in the control samples. Consistent with the mRNA level, the protein level of Ntn‑4 was also found to be decreased in NSCLC cells. However, both the function of Ntn‑4 and the underlying mechanisms of Ntn‑4 downregulation in NSCLC have yet to be fully elucidated. As was anticipated, the overexpression of Ntn‑4 led to a marked decrease in the proliferation, migration and invasion of NSCLC cells. Notably, RNA‑binding protein quaking 5 (Qki‑5) was found to exhibit antitumor activity in lung cancer, not only by enhancing the level of Ntn‑4 by binding to Ntn‑4 mRNA, but also by suppressing the proliferation, invasion and migration of NSCLC cells. However, Qki‑5 is known to be frequently downregulated in NSCLC. Moreover, the knockdown of Ntn‑4 was found to reverse the suppressive effects of Qki‑5 on NSCLC progression both in vitro and in vivo. Taken together, the findings of the present study demonstrate that Ntn‑4 is able to suppress the progression of NSCLC, and that the level of Ntn‑4 can be regulated by Qki‑5. Therefore, Ntn‑4 may be a novel diagnostic and therapeutic target for the treatment of NSCLC.

Mapping the splicing landscape of the human immune system

Most human genes code for more than one transcript. Different ratios of transcripts of the same gene can be found in different cell types or states, indicating differential use of transcription start sites or differential splicing. Such differential transcript use (DTUs) events provide an additional layer of regulation and protein diversity. With the exceptions of PTPRC and CIITA, there are very few reported cases of DTU events in the immune system. To rigorously map DTUs between different human immune cell types, we leveraged four publicly available RNA sequencing datasets. We identified 282 DTU events between five human healthy immune cell types that appear in at least two datasets. The patterns of the DTU events were mostly cell-type-specific or lineage-specific, in the context of the five cell types tested. DTUs correlated with the expression pattern of potential regulators, namely, splicing factors and transcription factors. Of the several immune related conditions studied, only sepsis affected the splicing of more than a few genes and only in innate immune cells. Taken together, we map the DTUs landscape in human peripheral blood immune cell types, and present hundreds of genes whose transcript use changes between cell types or upon activation.

Role of the RNA-binding protein family in gynecologic cancers

Gynecological cancers pose a threat to women's health. Although early-stage gynecological cancers show good outcomes after standardized treatment, the prognosis of patients with advanced, met-astatic, and recurrent cancers is poor. RNA-binding proteins (RBPs) are important cellular proteins that interact with RNA through RNA-binding domains and participate extensively in post-transcriptional regulatory processes, such as mRNA alternative splicing, polyadenylation, intracellular localization and stability, and translation. Abnormal RBP expression affects the normal function of oncogenes and tumor suppressor genes in many malignancies, thus leading to the occurrence or progression of cancers. Similarly, RBPs play crucial roles in gynecological carcinogenesis. We summarize the role of RBPs in gynecological malignancies and explore their potential in the diagnosis and treatment of cancers. The findings summarized in this review may provide a guide for future research on the functions of RBPs.

Multi-omics data analysis reveals the biological implications of alternative splicing events in lung adenocarcinoma

Cancer is characterized by the dysregulation of alternative splicing (AS). However, the comprehensive regulatory mechanisms of AS in lung adenocarcinoma (LUAD) are poorly understood. Here, we displayed the AS landscape in LUAD based on the integrated analyses of LUAD's multi-omics data. We identified 13,995 AS events in 6309 genes as differentially expressed alternative splicing events (DEASEs) mainly covering protein-coding genes. These DEASEs were strongly linked to "cancer hallmarks", such as apoptosis, DNA repair, cell cycle, cell proliferation, angiogenesis, immune response, generation of precursor metabolites and energy, p53 signaling pathway and PI3K-AKT signaling pathway. We further built a regulatory network connecting splicing factors (SFs) and DEASEs. In addition, RNA-binding protein (RBP) mutations that can affect DEASEs were investigated to find some potential cancer drivers. Further association analysis demonstrated that DNA methylation levels were highly correlated with DEASEs. In summary, our results can bring new insight into understanding the mechanism of AS and provide novel biomarkers for personalized medicine of LUAD.

The role of alternative splicing in lung cancer

Aberrant alternative splicing (AS) events are frequently observed in lung cancer, which can be attributed to aberrant gene AS, alterations in splicing regulatory factors, or changes in splicing regulatory mechanisms. Consequently, the dysregulation of alternative RNA splicing is the fundamental cause of lung cancer. In this review, we have summarized the pivotal role of AS in the development, progression, invasion, metastasis, angiogenesis, and drug resistance of lung cancer. Ultimately, this review emphasizes the potential of AS as biomarkers in lung cancer prognosis and diagnosis, and introduces some applications of AS isoform in the treatment of lung cancer. The comprehension of the AS may provide a glimmer of hope for the eradication of lung cancer.

miR-31/QKI-5 axis facilitates cell cycle progression of non-small-cell lung cancer cells by interacting and regulating p21 and CDK4/6 expressions

BACKGROUND

RNA-binding protein Quaking-5 (QKI-5), a major isoform of QKIs, inhibits tumor progression in non-small cell lung cancer (NSCLC). However, the underlying molecular mechanisms of QKI-5 in the cell cycle of NSCLC are still largely unknown.

METHODS

MTT, flow cytometry, and colony formation assays were used to investigate cellular phenotypic changes. Mice xenograft model was used to evaluate the antitumor activities of QKI-5. Co-immunoprecipitation, RNA immunoprecipitation (RIP), and RIP sequencing were used to investigate protein-protein interaction and protein-mRNA interaction.

RESULTS

The QKI-5 expression was downregulated in NSCLC tissues compared with that in paired normal adjacent lung tissues. Overexpression of QKI-5 inhibited NSCLC cell proliferative and colony forming ability. In addition, QKI-5 induced cell cycle arrest at G0/G1 phase through upregulating p21^Waf1/Cip1 (p21) expression and downregulating cyclin D1, cyclin-dependent kinase 4 (CDK4), and CDK6 expressions. Further analyses showed that QKI-5 interacts with p21 protein and CDK4, CDK6 mRNAs, suggesting a critical function of QKI-5 in cell cycle regulation. In agreement with in vitro study, the mouse xenograft models validated tumor suppressive functions of QKI-5 in vivo through altering cell cycle G1-phase-associated proteins. Moreover, we demonstrated that QKI-5 is a direct target of miR-31. The QKI-5 expression was anticorrelated with the miR-31 expression in NSCLC patient samples.

CONCLUSION

Our results suggest that the miR-31/QKI-5/p21-CDK4-CDK6 axis might have critical functions in the progression of NSCLC, and targeting this axis could serve as a potential therapeutic strategy for NSCLC.

Kinetics of mRNA nuclear export regulate innate immune response gene expression

The abundance and stimulus-responsiveness of mature mRNA is thought to be determined by nuclear synthesis, processing, and cytoplasmic decay. However, the rate and efficiency of moving mRNA to the cytoplasm almost certainly contributes, but has rarely been measured. Here, we investigated mRNA export rates for innate immune genes. We generated high spatio-temporal resolution RNA-seq data from endotoxin-stimulated macrophages and parameterized a mathematical model to infer kinetic parameters with confidence intervals. We find that the effective chromatin-to-cytoplasm export rate is gene-specific, varying 100-fold: for some genes, less than 5% of synthesized transcripts arrive in the cytoplasm as mature mRNAs, while others show high export efficiency. Interestingly, effective export rates do not determine temporal gene responsiveness, but complement the wide range of mRNA decay rates; this ensures similar abundances of short- and long-lived mRNAs, which form successive innate immune response expression waves.

Research progress on RNA−binding proteins in breast cancer

Breast cancer is the most common malignancy in women and has a high incidence rate and mortality. Abnormal regulation of gene expression plays an important role in breast cancer occurrence and development. RNA-binding proteins (RBPs) are one kind of the key regulators for gene expression. By interacting with RNA, RBPs are widely involved in RNA cutting, transport, editing, intracellular localization, and translation regulation. RBPs are important during breast cancer occurrence and progression by engaging in many aspects, like proliferation, migration, invasion, and stemness. Therefore, comprehensively understanding the role of RBPs in breast cancer progression can facilitate early diagnosis, timely treatment, and long-term survival and quality of life of breast cancer patients.

Pathway analysis of smoking-induced changes in buccal mucosal gene expression

Background

Cigarette smoking is the leading preventable cause of death worldwide, and it is the most common cause of oral cancers. This study aims to provide a deeper understanding of the molecular pathways in the oral cavity that are altered by exposure to cigarette smoke.

Methods

The gene expression dataset (accession number GSE8987, GPL96) of buccal mucosa samples from smokers (n = 5) and never smokers (n = 5) was downloaded from The National Center for Biotechnology Information's (NCBI) Gene Expression Omnibus (GEO) repository. Differential expression was ascertained via NCBI's GEO2R software, and Ingenuity Pathway Analysis (IPA) software was used to perform a pathway analysis.

Results

A total of 459 genes were found to be significantly differentially expressed in smoker buccal mucosa (p  < 0.05). A total of 261 genes were over-expressed while 198 genes were under-expressed. The top canonical pathways predicted by IPA were nitric oxide and reactive oxygen production at macrophages, macrophages/fibroblasts and endothelial cells in rheumatoid arthritis, and thyroid cancer pathways. The IPA upstream analysis predicted that the TP53, APP, SMAD3, and TNF proteins as well as dexamethasone drug would be top transcriptional regulators.

Conclusions

IPA highlighted critical pathways of carcinogenesis, mainly nitric oxide and reactive oxygen production at macrophages, and confirmed widespread injury in the buccal mucosa due to exposure to cigarette smoke. Our findings suggest that cigarette smoking significantly impacts gene pathways in the buccal mucosa and may highlight potential targets for treating the effects of cigarette smoking.

Supplementary Information

The online version contains supplementary material available at 10.1186/s43042-022-00268-y.

Switched alternative splicing events as attractive features in lung squamous cell carcinoma

Background

Alternative splicing (AS) plays important roles in transcriptome and proteome diversity. Its dysregulation has a close affiliation with oncogenic processes. This study aimed to evaluate AS-based biomarkers by machine learning algorithms for lung squamous cell carcinoma (LUSC) patients.

Method

The Cancer Genome Atlas (TCGA) database and TCGA SpliceSeq database were utilized. After data composition balancing, Boruta feature selection and Spearman correlation analysis were used for differentially expressed AS events. Random forests and a nested fivefold cross-validation were applied for lymph node metastasis (LNM) classifier building. Random survival forest combined with Cox regression model was performed for a prognostic model, based on which a nomogram was developed. Functional enrichment analysis and Spearman correlation analysis were also conducted to explore underlying mechanisms. The expression of some switch-involved AS events along with parent genes was verified by qRT-PCR with 20 pairs of normal and LUSC tissues.

Results

We found 16 pairs of splicing events from same parent genes which were strongly related to the splicing switch (intrapair correlation coefficient = − 1). Next, we built a reliable LNM classifier based on 13 AS events as well as a nice prognostic model, in which switched AS events behaved prominently. The qRT-PCR presented consistent results with previous bioinformatics analysis, and some AS events like ITIH5-10715-AT and QKI-78404-AT showed remarkable detection efficiency for LUSC.

Conclusion

AS events, especially switched ones from the same parent genes, could provide new insights into the molecular diagnosis and therapeutic drug design of LUSC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02429-2.

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.

The RNA binding protein QKI5 suppresses ovarian cancer via downregulating transcriptional coactivator TAZ

RNA-binding proteins (RBPs) are a set of proteins involved in many steps of post-transcriptional regulation to maintain cellular homeostasis. Ovarian cancer (OC) is the most deadly gynecological cancer, but the roles of RBPs in OC are not fully understood. Here, we reported that the RBP QKI5 was significantly negatively correlated with aggressive tumor stage and worse prognosis in serous OC patients. QKI5 could suppress the growth and metastasis of OC cells both in vitro and in vivo. Transcriptome analysis showed that QKI5 negatively regulated the expression of the transcriptional coactivator TAZ and its downstream targets (e.g., CTGF and CYR61). Mechanistically, QKI5 bound to TAZ mRNA and recruited EDC4, thus decreasing the stability of TAZ mRNA. Functionally, TAZ was involved in the QKI5-mediated tumor suppression of OC cells, and QKI5 expression was inversely correlated with TAZ, CTGF, and CYR61 expression in OC patients. Together, our study indicates that QKI5 plays a tumor-suppressive role and negatively regulates TAZ expression in OC.

QKI-Regulated Alternative Splicing Events in Cervical Cancer: Pivotal Mechanism and Potential Therapeutic Strategy

QKI is a vital regulator in RNA splicing and maturation, but its role in cervical cancer (CC) is little known. In this study, we found that QKI is decreased in human CC, and overexpression of QKI inhibits HeLa cell proliferation and promotes the apoptosis of cancer cells. We identified hundreds of endogenous QKI-regulated alternative splicing events (ASEs) and differentially expressed genes (DEGs) in QKI-overexpressed HeLa cells by RNA-seq and selectively validated their expression by quantitative reverse-transcription polymerase chain reaction. The gene ontology and Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis showed that QKI-regulated ASEs and DEGs were closely related to cancer, apoptosis, and transcriptional regulatory functions. In short, QKI may affect the occurrence and development of CC by regulating gene expression through AS.

Biology of the mRNA Splicing Machinery and Its Dysregulation in Cancer Providing Therapeutic Opportunities

Dysregulation of messenger RNA (mRNA) processing—in particular mRNA splicing—is a hallmark of cancer. Compared to normal cells, cancer cells frequently present aberrant mRNA splicing, which promotes cancer progression and treatment resistance. This hallmark provides opportunities for developing new targeted cancer treatments. Splicing of precursor mRNA into mature mRNA is executed by a dynamic complex of proteins and small RNAs called the spliceosome. Spliceosomes are part of the supraspliceosome, a macromolecular structure where all co-transcriptional mRNA processing activities in the cell nucleus are coordinated. Here we review the biology of the mRNA splicing machinery in the context of other mRNA processing activities in the supraspliceosome and present current knowledge of its dysregulation in lung cancer. In addition, we review investigations to discover therapeutic targets in the spliceosome and give an overview of inhibitors and modulators of the mRNA splicing process identified so far. Together, this provides insight into the value of targeting the spliceosome as a possible new treatment for lung cancer.

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.

The RNA-binding protein Quaking regulates multiciliated and basal cell abundance in the developing lung

RNA-binding proteins (RBPs) form complexes with RNA, changing how the RNA is processed and thereby regulating gene expression. RBPs are important sources of gene regulation during organogenesis, including the development of lungs. The RBP called Quaking (QK) is critical for embryogenesis, yet it has not been studied in the developing lung. Here, we show that QK is widely expressed during rat lung development and into adulthood. The QK isoforms QK5 and QK7 colocalize to the nuclei of nearly all lung cells. QK6 is present in the nuclei and cytoplasm of mesenchymal cells and is only present in the epithelium during branching morphogenesis. QK knockdown in embryonic lung explants caused a greater number of multiciliated cells to appear in the airways, at the expense of basal cells. The mRNA of multiciliated cell genes and the abundance of FOXJ1/SOX2+ cells increased after knockdown, whereas P63/SOX2+ cells decreased. The cytokine IL-6, a known regulator of multiciliated cell differentiation, had increased mRNA levels after QK knockdown, although protein levels remained unchanged. Further studies are necessary to confirm whether QK acts as a blocker for the IL-6-induced differentiation of basal cells into multiciliated cells, and a conditional QK knockout would likely lead to additional discoveries on QK's role during lung development.

Profiling of specific long non-coding RNA signatures identifies ST8SIA6-AS1 AS a novel target for breast cancer

BACKGROUND

Breast cancer is the most commonly diagnosed cancer among women and is also the leading cause of cancer death for which the treatment and methods of diagnosis remain unsatisfied. Long non-coding RNA (lncRNA) plays an important role in the occurrence and development of tumors, including breast cancer. We aimed to seek new and efficient treatment targets by analyzing the lncRNA expression profiles of breast cancer.

METHODS

A competitive endogenous RNA microarray was used to investigate the profiles of differentially expressed lncRNAs. Quantitative real-time polymerase chain reaction analysis (qRT-PCR) validated the top differentially expressed lncRNAs in 107 pairs of breast cancer tissues and adjacent normal tissues. cis- and trans-regulation mRNAs of lncRNAs were used to perform enrichment analysis. Cell function assays were used to explore the functions of ST8SIA6-AS1.

RESULTS

Seven lncRNAs, comprising ST8SIA6-AS1, lnc-HIST1H2BJ-5:1, lnc-PRICKLE2-3:2, RP1-86C11.7, RP11-15F12.1, ZNF670-ZNF695 and lnc-STRN3-12:1, were shown to be significantly up-regulated in breast cancer. lncRNA ST8SIA6-AS1 was associated with TNM staging and Ki-67 index. The cell function assays showed that ST8SIA6-AS1 can promote the proliferation, migration and invasion of breast cancer cells. The functions of ST8SIA6-AS1 were explored and the competing endogenous RNA mode showed that miR-4252 was a potential candidate. Its target genes were further predicted. The lncRNA-protein mode showed three potential candidate RNA binding proteins: NONO, QKI and RBMX.

CONCLUSIONS

lncRNA ST8SIA6-AS1 can promote the proliferation, migration and invasion of breast cancer cells. By hypothesizing two different functional modes of ST8SIA6-AS1, we found lncRNA ST8SIA6-AS1 may contribute to breast cancer progression through miR-4252 or interacting with RNA binding proteins: NONO, QKI and RBMX.

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.

N6-Methyladenosine RNA Methylation Regulator-Related Alternative Splicing (AS) Gene Signature Predicts Non-Small Cell Lung Cancer Prognosis

Aberrant N6-methyladenosine (m6A) RNA methylation regulatory genes and related gene alternative splicing (AS) could be used to predict the prognosis of non-small cell lung carcinoma. This study focused on 13 m6A regulatory genes (METTL3, METTL14, WTAP, KIAA1429, RBM15, ZC3H13, YTHDC1, YTHDC2, YTHDF1, YTHDF2, HNRNPC, FTO, and ALKBH5) and expression profiles in TCGA-LUAD (n = 504) and TCGA-LUSC (n = 479) datasets from the Cancer Genome Atlas database. The data were downloaded and bioinformatically and statistically analyzed, including the gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analyses. There were 43,948 mRNA splicing events in lung adenocarcinoma (LUAD) and 46,020 in lung squamous cell carcinoma (LUSC), and the data suggested that m6A regulators could regulate mRNA splicing. Differential HNRNPC and RBM15 expression was associated with overall survival (OS) of LUAD and HNRNPC and METTL3 expression with the OS of LUSC patients. Furthermore, the non-small cell lung cancer prognosis-related AS events signature was constructed and divided patients into high- vs. low-risk groups using seven and 14 AS genes in LUAD and LUSC, respectively. The LUAD risk signature was associated with gender and T, N, and TNM stages, but the LUSC risk signature was not associated with any clinical features. In addition, the risk signature and TNM stage were independent prognostic predictors in LUAD and the risk signature and T stage were independent prognostic predictors in LUSC after the multivariate Cox regression and receiver operating characteristic analyses. In conclusion, this study revealed the AS prognostic signature in the prediction of LUAD and LUSC prognosis.

RNA-binding protein QKI suppresses breast cancer via RASA1/MAPK signaling pathway

Background

RNA-binding protein Quaking (QKI) has been linked with the pathogenesis and development of various human malignancies. Herein, we explored the particular role of QKI in breast cancer (BC) progression.

Methods

The methods employed in the study included public dataset analysis, western blot, quantitative real-time PCR (qRT-PCR), cell count kit-8 (CCK8) assay, colony formation assay, flow cytometric analysis, RNA immunoprecipitation (RIP), messenger RNA (mRNA) stability assay, QKI overexpression and knockdown, and Ras p21 protein activator 1 (RASA1) knockdown.

Results

Aberrant expression levels of QKI and RASA1 were detected in BC and compared with those in noncancerous tissues. A moderately positive correlation between QKI and RASA1 was verified within BC tissues. Low expression of QKI was associated with positive estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER2) status, non-triple-negative breast cancer (TNBC), non-basal-like BC, and poor clinical outcomes in BC patients. QKI overexpression suppressed BC cell proliferation and colony formation, and arrested cell cycle at G1 phase. RIP assay and mRNA stability assay confirmed that QKI directly bound to RASA1 transcript and increased its stability, thus inactivating the MAPK pathway and inhibiting BC progression. RASA1 knockdown could partly attenuate the inhibitory effect of QKI on BC cell proliferation via activating the mitogen-activated protein kinase (MAPK) pathway.

Conclusions

QKI, which was frequently downregulated in BC, could significantly inhibit cell proliferation and arrest cell cycle at G1 phase by binding and enhancing RASA1 mRNA expression. Low expression of QKI was prominently associated with unfavorable clinical outcomes in BC patients, indicating the prognostic value of QKI in BC.

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.

Exploration of predictive and prognostic alternative splicing signatures in lung adenocarcinoma using machine learning methods

Background

Alternative splicing (AS) plays critical roles in generating protein diversity and complexity. Dysregulation of AS underlies the initiation and progression of tumors. Machine learning approaches have emerged as efficient tools to identify promising biomarkers. It is meaningful to explore pivotal AS events (ASEs) to deepen understanding and improve prognostic assessments of lung adenocarcinoma (LUAD) via machine learning algorithms.

Method

RNA sequencing data and AS data were extracted from The Cancer Genome Atlas (TCGA) database and TCGA SpliceSeq database. Using several machine learning methods, we identified 24 pairs of LUAD-related ASEs implicated in splicing switches and a random forest-based classifiers for identifying lymph node metastasis (LNM) consisting of 12 ASEs. Furthermore, we identified key prognosis-related ASEs and established a 16-ASE-based prognostic model to predict overall survival for LUAD patients using Cox regression model, random survival forest analysis, and forward selection model. Bioinformatics analyses were also applied to identify underlying mechanisms and associated upstream splicing factors (SFs).

Results

Each pair of ASEs was spliced from the same parent gene, and exhibited perfect inverse intrapair correlation (correlation coefficient = − 1). The 12-ASE-based classifier showed robust ability to evaluate LNM status of LUAD patients with the area under the receiver operating characteristic (ROC) curve (AUC) more than 0.7 in fivefold cross-validation. The prognostic model performed well at 1, 3, 5, and 10 years in both the training cohort and internal test cohort. Univariate and multivariate Cox regression indicated the prognostic model could be used as an independent prognostic factor for patients with LUAD. Further analysis revealed correlations between the prognostic model and American Joint Committee on Cancer stage, T stage, N stage, and living status. The splicing network constructed of survival-related SFs and ASEs depicts regulatory relationships between them.

Conclusion

In summary, our study provides insight into LUAD researches and managements based on these AS biomarkers.

QKI-5 regulates the alternative splicing of cytoskeletal gene ADD3 in lung cancer

Accumulating evidence indicates that the alternative splicing program undergoes extensive changes during cancer development and progression. The RNA-binding protein QKI-5 is frequently downregulated and exhibits anti-tumor activity in lung cancer. Howeve-r, little is known about the functional targets and regulatory mechanism of QKI-5. Here, we report that upregulation of exon 14 inclusion of cytoskeletal gene Adducin 3 (ADD3) significantly correlates with a poor prognosis in lung cancer. QKI-5 inhibits cell proliferation and migration in part through suppressing the splicing of ADD3 exon 14. Through genome-wide mapping of QKI-5 binding sites in vivo at nucleotide resolution by iCLIP-seq analysis, we found that QKI-5 regulates alternative splicing of its target mRNAs in a binding position-dependent manner. By binding to multiple sites in an upstream intron region, QKI-5 represses the splicing of ADD3 exon 14. We also identified several QKI mutations in tumors, which cause dysregulation of the splicing of QKI targets ADD3 and NUMB. Taken together, our results reveal that QKI-mediated alternative splicing of ADD3 is a key lung cancer-associated splicing event, which underlies in part the tumor suppressor function of QKI.

Transcriptome profile analysis reveals a silica-induced immune response and fibrosis in a silicosis rat model

Silicosis is a type of pneumoconiosis caused by the inhalation of silica dust. It is characterized by inflammation and fibrosis of the lung. Although many studies have reported that crystalline silica-inhalation into the lung initiates the immune response, activating effector cells and triggering the inflammatory cascade with subsequent elaboration of the extracellular matrix and fibrosis, the mechanism of silicosis pathogenesis remains unclear. In the present study, we established a silica inhalation-induced silicosis rat model validated by histological and cytokine analyses. RNA-seq and bioinformatic analyses showed that 600 genes were upregulated and 537 genes were downregulated in the silica-treated group. GO enrichment analysis indicates that these differentially expressed genes are enriched in several biological processes including immune response and organism remodeling. KEGG enrichment analysis showed that 53 enriched pathways were mainly associated with human diseases, immune response, signal transduction, and fibrosis process. Since alternative splicing of pre-mRNAs is also essential for the regulation of gene expression, we identified several alternative pre-mRNA splicing events in the fibrotic process. This study will provide a foundation to understand the molecular mechanism of the pulmonary fibrosis caused by silica.

RNA-binding proteins in tumor progression

RNA-binding protein (RBP) has a highly dynamic spatiotemporal regulation process and important biological functions. They are critical to maintain the transcriptome through post-transcriptionally controlling the processing and transportation of RNA, including regulating RNA splicing, polyadenylation, mRNA stability, mRNA localization, and translation. Alteration of each process will affect the RNA life cycle, produce abnormal protein phenotypes, and thus lead to the occurrence and development of tumors. Here, we summarize RBPs involved in tumor progression and the underlying molecular mechanisms whereby they are regulated and exert their effects. This analysis is an important step towards the comprehensive characterization of post-transcriptional gene regulation involved in tumor progression.

Upstream analysis of alternative splicing: a review of computational approaches to predict context-dependent splicing factors

Alternative splicing (AS) has shown to play a pivotal role in the development of diseases, including cancer. Specifically, all the hallmarks of cancer (angiogenesis, cell immortality, avoiding immune system response, etc.) are found to have a counterpart in aberrant splicing of key genes. Identifying the context-specific regulators of splicing provides valuable information to find new biomarkers, as well as to define alternative therapeutic strategies. The computational models to identify these regulators are not trivial and require three conceptual steps: the detection of AS events, the identification of splicing factors that potentially regulate these events and the contextualization of these pieces of information for a specific experiment. In this work, we review the different algorithmic methodologies developed for each of these tasks. Main weaknesses and strengths of the different steps of the pipeline are discussed. Finally, a case study is detailed to help the reader be aware of the potential and limitations of this computational approach.

SWATH-MS based proteomic profiling of pancreatic ductal adenocarcinoma tumours reveals the interplay between the extracellular matrix and related intracellular pathways

Pancreatic cancer accounts for 2.8% of new cancer cases worldwide and is projected to become the second leading cause of cancer-related deaths by 2030. Patients of African ancestry appear to be at an increased risk for pancreatic ductal adenocarcinoma (PDAC), with more severe disease and outcomes. The purpose of this study was to map the proteomic and genomic landscape of a cohort of PDAC patients of African ancestry. Thirty tissues (15 tumours and 15 normal adjacent tissues) were obtained from consenting South African PDAC patients. Optimisation of the sample preparation method allowed for the simultaneous extraction of high-purity protein and DNA for SWATH-MS and OncoArray SNV analyses. We quantified 3402 proteins with 49 upregulated and 35 downregulated proteins at a minimum 2.1 fold change and FDR adjusted p-value (q-value) ≤ 0.01 when comparing tumour to normal adjacent tissue. Many of the upregulated proteins in the tumour samples are involved in extracellular matrix formation (ECM) and related intracellular pathways. In addition, proteins such as EMIL1, KBTB2, and ZCCHV involved in the regulation of ECM proteins were observed to be dysregulated in pancreatic tumours. Downregulation of pathways involved in oxygen and carbon dioxide transport were observed. Genotype data showed missense mutations in some upregulated proteins, such as MYPN, ESTY2 and SERPINB8. Approximately 11% of the dysregulated proteins, including ISLR, BP1, PTK7 and OLFL3, were predicted to be secretory proteins. These findings help in further elucidating the biology of PDAC and may aid in identifying future plausible markers for the disease.

The RNA binding protein Quaking represses host interferon response by downregulating MAVS

Quaking (QKI) is an RNA-binding protein (RBP) involved in multiple aspects of RNA metabolism and many biological processes. Despite a known immune function in regulating monocyte differentiation and inflammatory responses, the degree to which QKI regulates the host interferon (IFN) response remains poorly characterized. Here we show that QKI ablation enhances poly(I:C) and viral infection-induced IFNβ transcription. Characterization of IFN-related signalling cascades reveals that QKI knockout results in higher levels of IRF3 phosphorylation. Interestingly, complementation with QKI-5 isoform alone is sufficient to rescue this phenotype and reduce IRF3 phosphorylation. Further analysis shows that MAVS, but not RIG-I or MDA5, is robustly upregulated in the absence of QKI, suggesting that QKI downregulates MAVS and thus represses the host IFN response. As expected, MAVS depletion reduces IFNβ activation and knockout of MAVS in the QKI knockout cells completely abolishes IFNβ induction. Consistently, ectopic expression of RIG-I activates stronger IFNβ induction via MAVS-IRF3 pathway in the absence of QKI. Collectively, these findings demonstrate a novel role for QKI in negatively regulating host IFN response by reducing MAVS levels.

QKI-6 inhibits bladder cancer malignant behaviours through down-regulating E2F3 and NF-κB signalling

Quaking homolog (QKI) is a member of the RNA-binding signal transduction and activator of proteins family. Previous studies showed that QKI possesses the tumour suppressor activity in human cancers by interacting with the 3'-untraslated region (3'-UTR) of various gene transcripts via the STAR domain. This study first assessed the association of QKI-6 expression with clinicopathological and survival data from bladder cancer patients and then investigated the underlying molecular mechanisms. Bladder cancer tissues (n = 223) were subjected to immunohistochemistry, and tumour cell lines and nude mice were used for different in vitro and in vivo assays following QKI-6 overexpression or knockdown. QKI-6 down-regulation was associated with advanced tumour TNM stages and poor patient overall survival. QKI-6 overexpression inhibited bladder cancer cell growth and invasion capacity, but induced tumour cell apoptosis and cell cycle arrest. Furthermore, ectopic expression of QKI-6 reduced tumour xenograft growth and expression of proliferation markers, Ki67 and PCNA. However, knockdown of QKI-6 expression had opposite effects in vitro and in vivo. QKI-6 inhibited expression of E2 transcription factor 3 (E2F3) by directly binding to the E2F3 3'-UTR, whereas E2F3 induced QKI-6 transcription by binding to the QKI-6 promoter in negative feedback mechanism. QKI-6 expression also suppressed activity and expression of nuclear factor-κB (NF-κB) signalling proteins in vitro, implying a novel multilevel regulatory network downstream of QKI-6. In conclusion, QKI-6 down-regulation contributes to bladder cancer development and progression.

psichomics: graphical application for alternative splicing quantification and analysis

Alternative pre-mRNA splicing generates functionally distinct transcripts from the same gene and is involved in the control of multiple cellular processes, with its dysregulation being associated with a variety of pathologies. The advent of next-generation sequencing has enabled global studies of alternative splicing in different physiological and disease contexts. However, current bioinformatics tools for alternative splicing analysis from RNA-seq data are not user-friendly, disregard available exon-exon junction quantification or have limited downstream analysis features. To overcome such limitations, we have developed psichomics, an R package with an intuitive graphical interface for alternative splicing quantification and downstream dimensionality reduction, differential splicing and gene expression and survival analyses based on The Cancer Genome Atlas, the Genotype-Tissue Expression project, the Sequence Read Archive project and user-provided data. These integrative analyses can also incorporate clinical and molecular sample-associated features. We successfully used psichomics in a laptop to reveal alternative splicing signatures specific to stage I breast cancer and associated novel putative prognostic factors.

Gene Expression Analyses in Non Muscle Invasive Bladder Cancer Reveals a Role for Alternative Splicing and Tp53 Status

Non-muscle invasive bladder cancer (NMIBC) represents a crucial problem for the national health care systems due to its high rates of recurrence and the consequent need of frequent follow-ups. Here, gene expression analyses in patients diagnosed as NMIBC were performed to determine those molecular pathways involved in tumor initiation, finding that both MYC and E2F are up regulated and helps to tumor initiation and progression. Our results also support an important involvement of alternative splicing events, modifying key pathways to favour bladder tumor evolution. Finally, since MDM2 showed differential exon usage, mutations in TP53 and its protein expression have been also studied in the same patients. Our data support that recurrence is epigenetically mediated and favoured by an increase protein expression of TP53, which appears more frequently mutated in advanced stages and grades, being associated to a worse prognosis. Therefore, TP53 mutational status could be used as a potential biomarker in the first stages of NMIBC to predict recurrence and prognosis.

Overexpression of miR-20a promotes the progression of osteosarcoma by directly targeting QKI2

Osteosarcoma (OS) is the most common type of malignant primary bone neoplasm. Although the application of neoadjuvant chemotherapy has improved the 5-year survival rate of patients suffering from OS, prognosis remains poor. Therefore, it is important to elucidate the molecular mechanisms underlying the occurrence, progression and metastasis of OS. The RNA-binding protein Quaking (QKI) is a member of the STAR family of proteins, and can function as a tumor suppressor gene to suppress the occurrence and progression of a variety of tumors; however, the role of QKI in OS remains to be fully elucidated. In the present study, it was identified that the expression of QKI2 was downregulated in OS using western blot analysis. In addition, subsequent functional investigations, including MTT, Transwell invasion and migration assays, revealed that QKI2 inhibited the proliferation, invasion and migration of an OS cell line in vitro. By implementing a series of experimental techniques in molecular biology, including reverse transcription-quantitative polymerase chain reaction and a double fluorescence reporter assay, it was demonstrated that the expression of miR-20a was high and inhibited the expression of QKI2 in OS. In conclusion, it was revealed that aberrantly upregulated miR-20a inhibited the expression of QKI2 in OS by targeting QKI2 mRNA, subsequently promoting the proliferation, migration and invasion of OS cells.

Integration of CLIP experiments of RNA-binding proteins: a novel approach to predict context-dependent splicing factors from transcriptomic data

BACKGROUND

Splicing is a genetic process that has important implications in several diseases including cancer. Deciphering the complex rules of splicing regulation is crucial to understand and treat splicing-related diseases. Splicing factors and other RNA-binding proteins (RBPs) play a key role in the regulation of splicing. The specific binding sites of an RBP can be measured using CLIP experiments. However, to unveil which RBPs regulate a condition, it is necessary to have a priori hypotheses, as a single CLIP experiment targets a single protein.

RESULTS

In this work, we present a novel methodology to predict context-specific splicing factors from transcriptomic data. For this, we systematically collect, integrate and analyze more than 900 CLIP experiments stored in four CLIP databases: POSTAR2, CLIPdb, DoRiNA and StarBase. The analysis of these experiments shows the strong coherence between the binding sites of RBPs of similar families. Augmenting this information with expression changes, we are able to correctly predict the splicing factors that regulate splicing in two gold-standard experiments in which specific splicing factors are knocked-down.

CONCLUSIONS

The methodology presented in this study allows the prediction of active splicing factors in either cancer or any other condition by only using the information of transcript expression. This approach opens a wide range of possible studies to understand the splicing regulation of different conditions. A tutorial with the source code and databases is available at https://gitlab.com/fcarazo.m/sfprediction .

Alternative splicing in lung cancer

Lung cancer has the highest mortality rate of all cancers worldwide. Lung cancer is a very heterogeneous disease that is often diagnosed at later stages which have a poor prognosis. Aberrant alternative splicing patterns found in lung cancer contribute to important cell functions. These include changes in splicing for the BCL2L1, MDM2, MDM4, NUMB and MET genes during lung tumourigenesis, to affect pathways involved in apoptosis, cell proliferation and cellular cohesion. Global analyses of RNASeq datasets suggest there may be many more potentially influential aberrant splicing events that need to be investigated in lung cancer. Changes in expression of the splicing factors that regulate alternative splicing events have also been identified in lung cancer. Of these, changes in expression of QKI, RBM4, RBM5, RBM6, RBM10 and SRSF1 proteins regulate many of the most frequently referenced aberrant splicing events in lung cancer. The expanding list of genes known to be aberrantly spliced in lung cancer along with the altered expression of splicing factors that regulate them are providing new clues as to how lung cancer develops, and how these events can be exploited for better treatment. This article is part of a Special Issue entitled: RNA structure and splicing regulation edited by Francisco Baralle, Ravindra Singh and Stefan Stamm.

YY1 Complex Promotes Quaking Expression via Super-Enhancer Binding during EMT of Hepatocellular Carcinoma

Quaking (QKI) is an alternative splicing factor that can regulate circRNA formation in the progression of epithelial-mesenchymal transition, but the mechanism remains unclear. High expression of QKI is correlated with short survival time, metastasis, and high clinical stage and pathology grade in hepatocellular carcinoma (HCC). Here we report that transcription of the QKI gene was activated by the Yin-Yang 1 (YY1)/p65/p300 complex, in which YY1 bound to the super-enhancer and promoter of QKI, p65 combined with the promoter, and p300 served as a mediator to maintain the stability of the complex. This YY1/p65/p300 complex increased QKI expression to promote the malignancy of HCC as well as an increased circRNA formation in vitro and in vivo. Hyperoside is one of several plant-derived flavonol glycoside compounds. Through virtual screening and antitumor activity analysis, we found that hyperoside inhibited QKI expression by targeting the YY1/p65/p300 complex. Overall, our study suggests that the regulatory mechanism of QKI depends on the YY1/p65/p300 complex and that it may serve as a potential target for treatment of HCC. SIGNIFICANCE: These findings identify the YY1/p65/p300 complex as a regulator of QKI expression, identifying several potential therapeutic targets for the treatment of HCC.

5 protein-based signature for resectable lung squamous cell carcinoma improves the prognostic performance of the TNM staging

INTRODUCTION

Prognostic biomarkers have been very elusive in the lung squamous cell carcinoma (SCC) and none is currently being used in the clinical setting. We aimed to identify and validate the clinical utility of a protein-based prognostic signature to stratify patients with early lung SCC according to their risk of recurrence or death.

METHODS

Patients were staged following the new International Association for the Study of Lung Cancer (IASLC) staging criteria (eighth edition, 2018). Three independent retrospective cohorts of 117, 96 and 105 patients with lung SCC were analysed to develop and validate a prognostic signature based on immunohistochemistry for five proteins.

RESULTS

We identified a five protein-based signature whose prognostic index (PI) was an independent and significant predictor of disease-free survival (DFS) (p<0.001; HR=4.06, 95% CI 2.18 to 7.56) and overall survival (OS) (p=0.004; HR=2.38, 95% CI 1.32 to 4.31). The prognostic capability of PI was confirmed in an external multi-institutional cohort for DFS (p=0.042; HR=2.01, 95% CI 1.03 to 3.94) and for OS (p=0.031; HR=2.29, 95% CI 1.08 to 4.86). Moreover, PI added complementary information to the newly established IASLC TNM 8th edition staging system. A combined prognostic model including both molecular and anatomical (TNM) criteria improved the risk stratification in both cohorts (p<0.05).

CONCLUSION

We have identified and validated a clinically feasible protein-based prognostic model that complements the updated TNM system allowing more accurate risk stratification. This signature may be used as an advantageous tool to improve the clinical management of the patients, allowing the reduction of lung SCC mortality through a more accurate knowledge of the patient's potential outcome.

Comparison of RNA-seq and microarray platforms for splice event detection using a cross-platform algorithm

Background

RNA-seq is a reference technology for determining alternative splicing at genome-wide level. Exon arrays remain widely used for the analysis of gene expression, but show poor validation rate with regard to splicing events. Commercial arrays that include probes within exon junctions have been developed in order to overcome this problem.

We compare the performance of RNA-seq (Illumina HiSeq) and junction arrays (Affymetrix Human Transcriptome array) for the analysis of transcript splicing events. Three different breast cancer cell lines were treated with CX-4945, a drug that severely affects splicing. To enable a direct comparison of the two platforms, we adapted EventPointer, an algorithm that detects and labels alternative splicing events using junction arrays, to work also on RNA-seq data. Common results and discrepancies between the technologies were validated and/or resolved by over 200 PCR experiments.

Results

As might be expected, RNA-seq appears superior in cases where the technologies disagree and is able to discover novel splicing events beyond the limitations of physical probe-sets. We observe a high degree of coherence between the two technologies, however, with correlation of EventPointer results over 0.90. Through decimation, the detection power of the junction arrays is equivalent to RNA-seq with up to 60 million reads.

Conclusions

Our results suggest, therefore, that exon-junction arrays are a viable alternative to RNA-seq for detection of alternative splicing events when focusing on well-described transcriptional regions.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5082-2) contains supplementary material, which is available to authorized users.

Dysregulated transcriptional networks in KMT2A- and MLLT10-rearranged T-ALL

For children and young adults with T-lineage acute lymphoblastic leukemia (T-ALL), event free survival following relapse is < 10%. We recently showed that rearrangements of the mixed lineage leukemia gene (KMT2A-R) are associated with induction failure and an inferior survival in T-ALL. Because there are currently no molecular features that inform treatment strategies in T-ALL, we hypothesized that transcriptional alterations related to KMT2A-R and MLLT10-R T-ALL could identify biologically relevant genes and signaling pathways for the development of targeted therapies for these groups of patients. We analyzed microarray data from a retrospective cohort of 100 T-ALL patients to identify novel targets for KMT2A (n = 12) or MLLT10 (n = 9) chimeras. We identified 330 probe sets that could discriminate between these groups, including novel targets, like RUNX2, TCF4 or MYO6. The results were further validated in two independent data sets and the functional networks were analyzed to identify pathways that may be of pathogenic or therapeutic relevance.

Alternative-splicing defects in cancer: Splicing regulators and their downstream targets, guiding the way to novel cancer therapeutics

Defects in alternative splicing are frequently found in human tumors and result either from mutations in splicing-regulatory elements of specific cancer genes or from changes in the regulatory splicing machinery. RNA splicing regulators have emerged as a new class of oncoproteins and tumor suppressors, and contribute to disease progression by modulating RNA isoforms involved in the hallmark cancer pathways. Thus, dysregulation of alternative RNA splicing is fundamental to cancer and provides a potentially rich source of novel therapeutic targets. Here, we review the alterations in splicing regulatory factors detected in human tumors, as well as the resulting alternatively spliced isoforms that impact cancer hallmarks, and discuss how they contribute to disease pathogenesis. RNA splicing is a highly regulated process and, as such, the regulators are themselves tightly regulated. Differential transcriptional and posttranscriptional regulation of splicing factors modulates their levels and activities in tumor cells. Furthermore, the composition of the tumor microenvironment can also influence which isoforms are expressed in a given cell type and impact drug responses. Finally, we summarize current efforts in targeting alternative splicing, including global splicing inhibition using small molecules blocking the spliceosome or splicing-factor-modifying enzymes, as well as splice-switching RNA-based therapeutics to modulate cancer-specific splicing isoforms. This article is categorized under: RNA in Disease and Development > RNA in Disease RNA Processing > Splicing Regulation/Alternative Splicing.

A novel protein-based prognostic signature improves risk stratification to guide clinical management in early-stage lung adenocarcinoma patients

Each of the pathological stages (I-IIIa) of surgically resected non-small-cell lung cancer has hidden biological heterogeneity, manifested as heterogeneous outcomes within each stage. Thus, the finding of robust and precise molecular classifiers with which to assess individual patient risk is an unmet medical need. Here, we identified and validated the clinical utility of a new prognostic signature based on three proteins (BRCA1, QKI, and SLC2A1) to stratify early-stage lung adenocarcinoma patients according to their risk of recurrence or death. Patients were staged according to the new International Association for the Study of Lung Cancer (IASLC) staging criteria (8th edition, 2018). A test cohort (n = 239) was used to assess the value of this new prognostic index (PI) based on the three proteins. The prognostic signature was developed by Cox regression with the use of stringent statistical criteria (TRIPOD: Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis). The model resulted in a highly significant predictor of 5-year outcome for disease-free survival (p < 0.001) and overall survival (p < 0.001). The prognostic ability of the model was externally validated in an independent multi-institutional cohort of patients (n = 114, p = 0.021). We also demonstrated that this molecular classifier adds relevant information to the gold standard TNM-based pathological staging, with a highly significant improvement of the likelihood ratio. We subsequently developed a combined PI including both the molecular and the pathological data that improved the risk stratification in both cohorts (p ≤ 0.001). Moreover, the signature may help to select stage I-IIA patients who might benefit from adjuvant chemotherapy. In summary, this protein-based signature accurately identifies those patients with a high risk of recurrence and death, and adds further prognostic information to the TNM-based clinical staging, even when the new IASLC 8th edition staging criteria are applied. More importantly, it may be a valuable tool for selecting patients for adjuvant therapy. Copyright © 2018 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

miR-200/375 control epithelial plasticity-associated alternative splicing by repressing the RNA-binding protein Quaking

Members of the miR‐200 family are critical gatekeepers of the epithelial state, restraining expression of pro‐mesenchymal genes that drive epithelial–mesenchymal transition (EMT) and contribute to metastatic cancer progression. Here, we show that miR‐200c and another epithelial‐enriched miRNA, miR‐375, exert widespread control of alternative splicing in cancer cells by suppressing the RNA‐binding protein Quaking (QKI). During EMT, QKI‐5 directly binds to and regulates hundreds of alternative splicing targets and exerts pleiotropic effects, such as increasing cell migration and invasion and restraining tumour growth, without appreciably affecting mRNA levels. QKI‐5 is both necessary and sufficient to direct EMT‐associated alternative splicing changes, and this splicing signature is broadly conserved across many epithelial‐derived cancer types. Importantly, several actin cytoskeleton‐associated genes are directly targeted by both QKI and miR‐200c, revealing coordinated control of alternative splicing and mRNA abundance during EMT. These findings demonstrate the existence of a miR‐200/miR‐375/QKI axis that impacts cancer‐associated epithelial cell plasticity through widespread control of alternative splicing.

Identification and characterization of host proteins bound to dengue virus 3' UTR reveal an antiviral role for quaking proteins

The four dengue viruses (DENV1-4) are rapidly reemerging infectious RNA viruses. These positive-strand viral genomes contain structured 3' untranslated regions (UTRs) that interact with various host RNA binding proteins (RBPs). These RBPs are functionally important in viral replication, pathogenesis, and defense against host immune mechanisms. Here, we combined RNA chromatography and quantitative mass spectrometry to identify proteins interacting with DENV1-4 3' UTRs. As expected, RBPs displayed distinct binding specificity. Among them, we focused on quaking (QKI) because of its preference for the DENV4 3' UTR (DENV-4/SG/06K2270DK1/2005). RNA immunoprecipitation experiments demonstrated that QKI interacted with DENV4 genomes in infected cells. Moreover, QKI depletion enhanced infectious particle production of DENV4. On the contrary, QKI did not interact with DENV2 3' UTR, and DENV2 replication was not affected consistently by QKI depletion. Next, we mapped the QKI interaction site and identified a QKI response element (QRE) in DENV4 3' UTR. Interestingly, removal of QRE from DENV4 3' UTR abolished this interaction and increased DENV4 viral particle production. Introduction of the QRE to DENV2 3' UTR led to QKI binding and reduced DENV2 infectious particle production. Finally, reporter assays suggest that QKI reduced translation efficiency of viral RNA. Our work describes a novel function of QKI in restricting viral replication.

Alternative Splicing as a Target for Cancer Treatment

Alternative splicing is a key mechanism determinant for gene expression in metazoan. During alternative splicing, non-coding sequences are removed to generate different mature messenger RNAs due to a combination of sequence elements and cellular factors that contribute to splicing regulation. A different combination of splicing sites, exonic or intronic sequences, mutually exclusive exons or retained introns could be selected during alternative splicing to generate different mature mRNAs that could in turn produce distinct protein products. Alternative splicing is the main source of protein diversity responsible for 90% of human gene expression, and it has recently become a hallmark for cancer with a full potential as a prognostic and therapeutic tool. Currently, more than 15,000 alternative splicing events have been associated to different aspects of cancer biology, including cell proliferation and invasion, apoptosis resistance and susceptibility to different chemotherapeutic drugs. Here, we present well established and newly discovered splicing events that occur in different cancer-related genes, their modification by several approaches and the current status of key tools developed to target alternative splicing with diagnostic and therapeutic purposes.