TSVdb: a web-tool for TCGA splicing variants analysis

Differential MYC and PROM1 mRNA isoform expression in breast invasive carcinoma as biomarkers for subtyping and prognosis

BACKGROUND

Cancer stem cells are a subpopulation of tumor cells with the ability to self-renew; evidence suggests that cancer stem cells are responsible for recurrence, metastasis, and resistance to therapy. MYC and CD133 (PROM1 gene) are clinical biomarkers for cancer stem cells, and their dysregulation is involved in the progression of many cancers. Alternative splicing of these genes may contribute to cancer stem cell differentiation.

METHODS

Transcriptional and clinical data of PROM1 and MYC mRNA isoforms in breast cancer samples were downloaded from the TCGA Splicing Variants Database site, a web-tool to explore mRNA alternative-splicing based on TCGA samples. Data include RSEM isoform expression, clinical sample types, survival data, and clinical receptor expression. Breast cancer subtypes (luminal A, luminal B, Her2 positive, triple negative) were assigned on the basis of estrogen, progesterone, and HER2 expression.

RESULTS

Expression of MYC isoforms uc003ysh.1 and uc003ysi.3 was significantly greater in triple-negative breast cancer compared with all other breast cancer subtypes (P < .001). Isoform uc003ysi.3 was associated with greater 5-year survival in luminal A breast cancer (hazard ratio, 0.79; 95% confidence interval, 0.65-0.96; P = .02). PROM1 isoforms uc003gop.2, uc003goq.3, uc003gos.2, and uc003gou.2 were expressed greatest in triple-negative breast cancer (P < .001). PROM1 isoform uc003gou.2 was associated with better 5-year survival in luminal A breast cancer (hazard ratio, 0.79; 95% confidence interval, 0.65-0.97; P = .02).

CONCLUSIONS

MYC and PROM1 isoforms are differentially expressed in breast cancer subtypes. Certain isoforms confer better survival prognosis. Further work should be done to study alternative splicing in cancer stem cells.

CD4+FOXP3Exon2+ regulatory T cell frequency predicts breast cancer prognosis and survival

CD4+FOXP3+ regulatory T cells (Tregs) suppress immune responses to tumors, and their accumulation in the tumor microenvironment (TME) correlates with poor clinical outcome in several cancers, including breast cancer (BC). However, the properties of intratumoral Tregs remain largely unknown. Here, we found that a functionally distinct subpopulation of Tregs, expressing the FOXP3 Exon2 splicing variants, is prominent in patients with hormone receptor-positive BC with poor prognosis. Notably, a comprehensive examination of the TCGA validated FOXP3E2 as an independent prognostic marker in all other BC subtypes. We found that FOXP3E2 expression underlies BCs with defective mismatch repair and a stem-like signature and highlights pathways involved in tumor survival. Last, we found that the TME induces FOXP3E2 through the CXCL12/CXCR4 axis and confirmed the higher immunosuppressive capacity of FOXP3E2+ Tregs derived from patients with BC. Our study suggests that FOXP3E2+ Tregs might be used as an independent biomarker to predict BC prognosis and survival and to develop super-targeted immunotherapies.

A study of the role of androgen receptor and androgen receptor variant 7 in TNBC patients and the effect of their targeting by Enzalutamide and EPI-001 in MDA-MB-231

The lack of targeted therapy for triple-negative breast cancer (TNBC) is among the mainsprings of its poor prognosis. This study aimed to elucidate the role of the androgen receptor (AR) and its splice variant 7 (ARv7) in TNBC patients. Further, the molecular impact of their blockers, Enzalutamide and EPI-001, on the TNBC cell line MDA-MB-231 was investigated. Thereby, immunohistochemical expression of AR/ARv7 was assessed for TNBC Egyptian patients. Moreover, bioinformatics analysis of AR/ARv7 RNA status was carried out on TNBC patients from The Cancer Genome Atlas Breast Carcinoma project (TCGA-BRCA). Data from both groups was correlated with patients' clinicopathological features. Besides, scratch wound healing assay and ELISA were employed to assess the effect of AR/ARv7 blockers on several metastasis markers in MDA-MB-231 cell line. In the Egyptian-TNBC patients, AR expression was associated with worse 7-year DFS (40.6 ± 18.6 %). In addition, ARv7 showed cytoplasmic and nuclear patterns, and both cytoplasmic and nuclear ARv7+ patients demonstrated a worse 7-year DFS (22.7 ± 17.7 % and 20 ± 17.9 %) and overall survival (63.6 ± 14.5 % and 40 ± 21.8 %). Importantly, 80 % of the nuclear ARv7+ patients developed distant metastasis. The data of the TCGA-TNBC patients showed a tendency for poor outcomes in the high ARv7-expressing patients. Molecularly, in MDA-MB-231, both inhibitors modulated metastasis and epithelial to mesenchymal transition (EMT) markers ROCK1, ROCK2, c-Myc, E-cadherin and N-cadherin, with EPI-001 downregulating NF-ĸB level as well. We concluded that ARv7 indicated poor prognosis in the studied cohorts and that blocking of AR/ARv7 abated metastasis and key regulators of EMT in MDA-MB-231, at least in part by targeting ROCK/NF-ĸB/c-Myc axis.

Cat-E: A comprehensive web tool for exploring cancer targeting strategies

Identifying potential cancer-associated genes and drug targets from omics data is challenging due to its diverse sources and analyses, requiring advanced skills and large amounts of time. To facilitate such analysis, we developed Cat-E (Cancer Target Explorer), a novel R/Shiny web tool designed for comprehensive analysis with evaluation according to cancer-related omics data. Cat-E is accessible at https://cat-e.bioinfo-wuerz.eu/. Cat-E compiles information on oncolytic viruses, cell lines, gene markers, and clinical studies by integrating molecular datasets from key databases such as OvirusTB, TCGA, DrugBANK, and PubChem. Users can use all datasets and upload their data to perform multiple analyses, such as differential gene expression analysis, metabolic pathway exploration, metabolic flux analysis, GO and KEGG enrichment analysis, survival analysis, immune signature analysis, single nucleotide variation analysis, dynamic analysis of gene expression changes and gene regulatory network changes, and protein structure prediction. Cancer target evaluation by Cat-E is demonstrated here on lung adenocarcinoma (LUAD) datasets. By offering a user-friendly interface and detailed user manual, Cat-E eliminates the need for advanced computational expertise, making it accessible to experimental biologists, undergraduate and graduate students, and oncology clinicians. It serves as a valuable tool for investigating genetic variations across diverse cancer types, facilitating the identification of novel diagnostic markers and potential therapeutic targets.

Lower FGFR2 mRNA Expression and Higher Levels of FGFR2 IIIc in HER2-Positive Breast Cancer

Fibroblast growth factor receptor 2 (FGFR2) has been associated with breast cancer. We performed in silico analyses to investigate the FGFR2 mRNA expression and splice variants associated with breast cancer subtypes. Online databases, including cBioPortal and TCGA SpliceSeq, were used to examine the association between the FGFR2 expression and splice variants with breast cancer subtypes. A higher FGFR2 mRNA was significantly associated with luminal, oestrogen receptor (ER)-positive breast cancers, and invasive lobular carcinomas, whereas a lower FGFR2 was associated with human epidermal growth factor receptor 2 (HER2)-positive breast cancer and invasive ductal carcinomas. The epithelial alternatively spliced FGFR2 IIIb isoform was significantly enriched in ER+ breast cancer, while the mesenchymal FGFR2 IIIc isoform was significantly prevalent in HER2+ cancer. Increased levels of FGFR2 and IIIb splice isoforms are associated with less aggressive breast cancer phenotypes, while decreased levels of FGFR2 and increased IIIc splice isoform are associated with more aggressive phenotypes.

Overexpression of Osteopontin-a and Osteopontin-c Splice Variants Are Worse Prognostic Features in Colorectal Cancer

Background: Osteopontin (OPN) is a glycoprotein involved in various physiological and pathological processes, and its aberrant expression in cancer cells is closely linked to tumor progression. In colorectal cancer (CRC), OPN is overexpressed, but the roles of its splice variants (OPN-SVs), OPNa, OPNb, and OPNc, are not well understood. This study aimed to characterize the expression patterns of OPN-SVs and their potential diagnostic and prognostic implications in CRC using transcriptomic data deposited in TSVdb and TCGA. Methods: The expression patterns of each OPN-SV were analyzed using transcriptomic data deposited in TSVdb and TCGA, which were correlated to patient data available at cBioPortal. Results: Bioinformatic analysis revealed that OPNa, OPNb, and OPNc are overexpressed in CRC samples compared to non-tumor samples. Notably, OPNa and OPNc are overexpressed in CRC stages (II, III, and IV) compared to stage I. Higher levels of OPNa and OPNc transcripts are associated with worse overall survival (OS) and shorter progression-free survival (PFS) in CRC patients. Additionally, the expression of OPNa, OPNb, and OPNc is correlated with BRAFV600E mutations in CRC samples. Conclusions: These findings suggest that OPNa and OPNc, in particular, have potential as diagnostic and prognostic biomarkers, paving the way for their further evaluation in CRC diagnosis and prognosis.

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Variant Impact Predictor database (VIPdb), version 2: Trends from 25 years of genetic variant impact predictors

Background

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

Results

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past 25 years, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 186 VIPs, resulting in a total of 403 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

Conclusions

VIPdb version 2 summarizes 403 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications.

Availability

VIPdb version 2 is available at https://genomeinterpretation.org/vipdb.

Isoform alterations in the ubiquitination machinery impacting gastrointestinal malignancies

The advancement of RNAseq and isoform-specific expression platforms has led to the understanding that isoform changes can alter molecular signaling to promote tumorigenesis. An active area in cancer research is uncovering the roles of ubiquitination on spliceosome assembly contributing to transcript diversity and expression of alternative isoforms. However, the effects of isoform changes on functionality of ubiquitination machineries (E1, E2, E3, E4, and deubiquitinating (DUB) enzymes) influencing onco- and tumor suppressor protein stabilities is currently understudied. Characterizing these changes could be instrumental in improving cancer outcomes via the identification of novel biomarkers and targetable signaling pathways. In this review, we focus on highlighting reported examples of direct, protein-coded isoform variation of ubiquitination enzymes influencing cancer development and progression in gastrointestinal (GI) malignancies. We have used a semi-automated system for identifying relevant literature and applied established systems for isoform categorization and functional classification to help structure literature findings. The results are a comprehensive snapshot of known isoform changes that are significant to GI cancers, and a framework for readers to use to address isoform variation in their own research. One of the key findings is the potential influence that isoforms of the ubiquitination machinery have on oncoprotein stability.

Tropomyosin1 isoforms underlie epithelial to mesenchymal plasticity, metastatic dissemination, and resistance to chemotherapy in high-grade serous ovarian cancer

Phenotypic plasticity, defined as the ability of individual cells with stable genotypes to exert different phenotypes upon exposure to specific environmental cues, represent the quintessential hallmark of the cancer cell en route from the primary lesion to distant organ sites where metastatic colonization will occur. Phenotypic plasticity is driven by a broad spectrum of epigenetic mechanisms that allow for the reversibility of epithelial-to-mesenchymal and mesenchymal-to-epithelial transitions (EMT/MET). By taking advantage of the co-existence of epithelial and quasi-mesenchymal cells within immortalized cancer cell lines, we have analyzed the role of EMT-related gene isoforms in the regulation of epithelial mesenchymal plasticity (EMP) in high grade serous ovarian cancer. When compared with colon cancer, a distinct spectrum of downstream targets characterizes quasi-mesenchymal ovarian cancer cells, likely to reflect the different modalities of metastasis formation between these two types of malignancy, i.e. hematogenous in colon and transcoelomic in ovarian cancer. Moreover, upstream RNA-binding proteins differentially expressed between epithelial and quasi-mesenchymal subpopulations of ovarian cancer cells were identified that underlie differential regulation of EMT-related isoforms. In particular, the up- and down-regulation of RBM24 and ESRP1, respectively, represent a main regulator of EMT in ovarian cancer cells. To validate the functional and clinical relevance of our approach, we selected and functionally analyzed the Tropomyosin 1 gene (TPM1), encoding for a protein that specifies the functional characteristics of individual actin filaments in contractile cells, among the ovarian-specific downstream AS targets. The low-molecular weight Tpm1.8/9 isoforms are specifically expressed in patient-derived ascites and promote invasion through activation of EMT and Wnt signaling, together with a broad spectrum of inflammation-related pathways. Moreover, Tpm1.8/9 expression confers resistance to taxane- and platinum-based chemotherapy. Small molecule inhibitors that target the Tpm1 isoforms support targeting Tpm1.8/9 as therapeutic targets for the development of future tailor-made clinical interventions.

Super enhancer loci of EGFR regulate EGFR variant 8 through enhancer RNA and strongly associate with survival in HNSCCs

Epidermal growth factor receptor (EGFR) has been shown to be overexpressed in human cancers due to mutation, amplification, and epigenetic hyperactivity, which leads to deregulated transcriptional mechanism. Among the eight different EGFR isoforms, the mechanism of regulation of full-length variant 1 is well-known, no studies have examined the function & factors regulating the expression of variant 8. This study aimed to understand the function of EGFR super-enhancer loci and its associated transcription factors regulating the expression of EGFR variant 8. Our study shows that overexpression of variant 8 and its transcription was more prevalent than variant 1 in many cancers and positively correlated with the EGFR-AS1 expression in oral cancer and HNSCC. Notably, individuals overexpressing variant 8 showed shorter overall survival and had a greater connection with other clinical traits than patients with overexpression of variant 1. In this study, TCGA enhancer RNA profiling on the constituent enhancer (CE1 and CE2) region revealed that the multiple enhancer RNAs formed from CE2 by employing CE1 as a promoter. Our bioinformatic analysis further supports the enrichment of enhancer RNA specific chromatin marks H3K27ac, H3K4me1, POL2 and H2AZ on CE2. GeneHancer and 3D chromatin capture analysis showed clustered interactions between CE1, CE2 loci and this interaction may regulates expression of both EGFR-eRNA and variant 8. Moreover, increased expression of SNAI2 and its close relationship to EGFR-AS1 and variant 8 suggest that SNAI2 could regulates variant 8 overexpression by building a MegaTrans complex with both EGFR-eRNA and EGFR-AS1. Our findings show that EGFR variant 8 and its transcriptional regulation & chromatin modification by eRNAs may provide a rationale for targeting RNA splicing in combination with targeted EGFR therapies in cancer.

MAJIQlopedia: an encyclopedia of RNA splicing variations in human tissues and cancer

Quantification of RNA splicing variations based on RNA-Sequencing can reveal tissue- and disease-specific splicing patterns. To study such splicing variations, we introduce MAJIQlopedia, an encyclopedia of splicing variations that encompasses 86 human tissues and 41 cancer datasets. MAJIQlopedia reports annotated and unannotated splicing events for a total of 486 175 alternative splice junctions in normal tissues and 338 317 alternative splice junctions in cancer. This database, available at https://majiq.biociphers.org/majiqlopedia/, includes a user-friendly interface that provides graphical representations of junction usage quantification for each junction across all tissue or cancer types. To demonstrate case usage of MAJIQlopedia, we review splicing variations in genes WT1, MAPT and BIN1, which all have known tissue or cancer-specific splicing variations. We also use MAJIQlopedia to highlight novel splicing variations in FDX1 and MEGF9 in normal tissues, and we uncover a novel exon inclusion event in RPS6KA6 that only occurs in two cancer types. Users can download the database, request the addition of data to the webtool, or install a MAJIQlopedia server to integrate proprietary data. MAJIQlopedia can serve as a reference database for researchers seeking to understand what splicing variations exist in genes of interest, and those looking to understand tissue- or cancer-specific splice isoform usage.

Pan-cancer transcriptomic data of ABI1 transcript variants and molecular constitutive elements identifies novel cancer metastatic and prognostic biomarkers

BACKGROUND

Abelson interactor 1 (ABI1) is associated with the metastasis and prognosis of many malignancies. The association between ABI1 transcript spliced variants, their molecular constitutive exons and exon-exon junctions (EEJs) in 14 cancer types and clinical outcomes remains unsolved.

OBJECTIVE

To identify novel cancer metastatic and prognostic biomarkers from ABI1 total mRNA, TSVs, and molecular constitutive elements.

METHODS

Using data from TCGA and TSVdb database, the standard median of ABI1 total mRNA, TSV, exon, and EEJ expression was used as a cut-off value. Kaplan-Meier analysis, Chi-squared test (X2) and Kendall's tau statistic were used to identify novel metastatic and prognostic biomarkers, and Cox regression analysis was performed to screen and identify independent prognostic factors.

RESULTS

A total of 35 ABI1-related factors were found to be closely related to the prognosis of eight candidate cancer types. A total of 14 ABI1 TSVs and molecular constitutive elements were identified as novel metastatic and prognostic biomarkers in four cancer types. A total of 13 ABI1 molecular constitutive elements were identified as independent prognostic biomarkers in six cancer types.

CONCLUSIONS

In this study, we identified 14 ABI1-related novel metastatic and prognostic markers and 21 independent prognostic factors in total 8 candidate cancer types.

Predicting the Structural Impact of Human Alternative Splicing

Protein structure prediction with neural networks is a powerful new method for linking protein sequence, structure, and function, but structures have generally been predicted for only a single isoform of each gene, neglecting splice variants. To investigate the structural implications of alternative splicing, we used AlphaFold2 to predict the structures of more than 11,000 human isoforms. We employed multiple metrics to identify splicing-induced structural alterations, including template matching score, secondary structure composition, surface charge distribution, radius of gyration, accessibility of post-translational modification sites, and structure-based function prediction. We identified examples of how alternative splicing induced clear changes in each of these properties. Structural similarity between isoforms largely correlated with degree of sequence identity, but we identified a subset of isoforms with low structural similarity despite high sequence similarity. Exon skipping and alternative last exons tended to increase the surface charge and radius of gyration. Splicing also buried or exposed numerous post-translational modification sites, most notably among the isoforms of BAX. Functional prediction nominated numerous functional differences among isoforms of the same gene, with loss of function compared to the reference predominating. Finally, we used single-cell RNA-seq data from the Tabula Sapiens to determine the cell types in which each structure is expressed. Our work represents an important resource for studying the structure and function of splice isoforms across the cell types of the human body.

FilGAP regulates tumor growth in Glioma through the regulation of mTORC1 and mTORC2

The mechanistic target of rapamycin (mTOR) is a serine/threonine protein kinase that forms the two different protein complexes, known as mTORC1 and mTORC2. mTOR signaling is activated in a variety of tumors, including glioma that is one of the malignant brain tumors. FilGAP (ARHGAP24) is a negative regulator of Rac, a member of Rho family small GTPases. In this study, we found that FilGAP interacts with mTORC1/2 and is involved in tumor formation in glioma. FilGAP interacted with mTORC1 via Raptor and with mTORC2 via Rictor and Sin1. Depletion of FilGAP in KINGS-1 glioma cells decreased phosphorylation of S6K and AKT. Furthermore, overexpression of FilGAP increased phosphorylation of S6K and AKT, suggesting that FilGAP activates mTORC1/2. U-87MG, glioblastoma cells, showed higher mTOR activity than KINGS-1, and phosphorylation of S6K and AKT was not affected by suppression of FilGAP expression. However, in the presence of PI3K inhibitors, phosphorylation of S6K and AKT was also decreased in U-87MG by depletion of FilGAP, suggesting that FilGAP may also regulate mTORC2 in U-87MG. Finally, we showed that depletion of FilGAP in KINGS-1 and U-87MG cells significantly reduced spheroid growth. These results suggest that FilGAP may contribute to tumor growth in glioma by regulating mTORC1/2 activities.

Alternative splicing of its 5'-UTR limits CD20 mRNA translation and enables resistance to CD20-directed immunotherapies

Aberrant skipping of coding exons in CD19 and CD22 compromises the response to immunotherapy in B-cell malignancies. Here, we showed that the MS4A1 gene encoding human CD20 also produces several messenger RNA (mRNA) isoforms with distinct 5' untranslated regions. Four variants (V1-4) were detected using RNA sequencing (RNA-seq) at distinct stages of normal B-cell differentiation and B-lymphoid malignancies, with V1 and V3 being the most abundant. During B-cell activation and Epstein-Barr virus infection, redirection of splicing from V1 to V3 coincided with increased CD20 positivity. Similarly, in diffuse large B-cell lymphoma, only V3, but not V1, correlated with CD20 protein levels, suggesting that V1 might be translation-deficient. Indeed, the longer V1 isoform contained upstream open reading frames and a stem-loop structure, which cooperatively inhibited polysome recruitment. By modulating CD20 isoforms with splice-switching morpholino oligomers, we enhanced CD20 expression and anti-CD20 antibody rituximab-mediated cytotoxicity in a panel of B-cell lines. Furthermore, reconstitution of CD20-knockout cells with V3 mRNA led to the recovery of CD20 positivity, whereas V1-reconstituted cells had undetectable levels of CD20 protein. Surprisingly, in vitro CD20-directed chimeric antigen receptor T cells were able to kill both V3- and V1-expressing cells, but the bispecific T-cell engager mosunetuzumab was only effective against V3-expressing cells. To determine whether CD20 splicing is involved in immunotherapy resistance, we performed RNA-seq on 4 postmosunetuzumab follicular lymphoma relapses and discovered that in 2 of them, the downregulation of CD20 was accompanied by a V3-to-V1 shift. Thus, splicing-mediated mechanisms of epitope loss extend to CD20-directed immunotherapies.

Formation of an invasion-permissive matrix requires TGFβ/SNAIL1-regulated alternative splicing of fibronectin

BACKGROUND

As in most solid cancers, the emergence of cells with oncogenic mutations in the mammary epithelium alters the tissue homeostasis. Some soluble factors, such as TGFβ, potently modify the behavior of healthy stromal cells. A subpopulation of cancer-associated fibroblasts expressing a TGFβ target, the SNAIL1 transcription factor, display myofibroblastic abilities that rearrange the stromal architecture. Breast tumors with the presence of SNAIL1 in the stromal compartment, and with aligned extracellular fiber, are associated with poor survival prognoses.

METHODS

We used deep RNA sequencing and biochemical techniques to study alternative splicing and human tumor databases to test for associations (correlation t-test) between SNAIL1 and fibronectin isoforms. Three-dimensional extracellular matrices generated from fibroblasts were used to study the mechanical properties and actions of the extracellular matrices on tumor cell and fibroblast behaviors. A metastatic mouse model of breast cancer was used to test the action of fibronectin isoforms on lung metastasis.

RESULTS

In silico studies showed that SNAIL1 correlates with the expression of the extra domain A (EDA)-containing (EDA+) fibronectin in advanced human breast cancer and other types of epithelial cancers. In TGFβ-activated fibroblasts, alternative splicing of fibronectin as well as of 500 other genes was modified by eliminating SNAIL1. Biochemical analyses demonstrated that SNAIL1 favors the inclusion of the EDA exon by modulating the activity of the SRSF1 splicing factor. Similar to Snai1 knockout fibroblasts, EDA- fibronectin fibroblasts produce an extracellular matrix  that does not sustain TGFβ-induced fiber organization, rigidity, fibroblast activation, or tumor cell invasion. The presence of EDA+ fibronectin changes the action of metalloproteinases on fibronectin fibers. Critically, in an mouse orthotopic breast cancer model, the absence of the fibronectin EDA domain completely prevents lung metastasis.

CONCLUSIONS

Our results support the requirement of EDA+ fibronectin in the generation of a metastasis permissive stromal architecture in breast cancers and its molecular control by SNAIL1. From a pharmacological point of view, specifically blocking EDA+ fibronectin deposition could be included in studies to reduce the formation of a pro-metastatic environment.

MDM2-Targeting Reassembly Peptide (TRAP) Nanoparticles for p53-Based Cancer Therapy

Gene mutations and functional inhibition are the major obstacles for p53-mediated oncotherapy. For p53-wild-type tumors, the underlying mechanisms of functional inhibition of p53 during oncogenesis are unknown. The results reveal that the expression of the MDM2 inhibitor ARF is inhibited in p53-wild-type tumors, indicating that the restoration of ARF could be a potential oncotherapy strategy for p53-wild-type tumors. Therefore, ARF-mimetic MDM2-targeting reassembly peptide nanoparticles (MtrapNPs) for p53-based tumor therapy is developed. The results elucidated that the MtrapNPs respond to and form a nanofiber structure with MDM2. By trapping MDM2, the MtrapNPs stabilize and activate p53 for the inhibition of p53-wild-type tumors. In most cases, reactivated mutant p53 is inhibited and degraded by MDM2. In the present study, MtrapNPs are used to load and deliver arsenic trioxide, a p53 mutation rescuer, for p53-mutated tumor treatment in both orthotopic and metastatic models, and they exhibit significant therapeutic effects. Therefore, the study provides evidence supporting a link between decreased ARF expression and tumor development in patients with p53-wild-type tumors. Thus, the MDM2-trap strategy, which addresses both the inhibition and mutations of p53, is an efficient strategy for the treatment of p53-mutated tumors.

A planar polarized MYO6-DOCK7-RAC1 axis promotes tissue fluidification in mammary epithelia

Tissue fluidification and collective motility are pivotal in regulating embryonic morphogenesis, wound healing, and tumor metastasis. These processes frequently require that each cell constituent of a tissue coordinates its migration activity and directed motion through the oriented extension of lamellipodium cell protrusions, promoted by RAC1 activity. While the upstream RAC1 regulators in individual migratory cells or leader cells during invasion or wound healing are well characterized, how RAC1 is controlled in follower cells remains unknown. Here, we identify a MYO6-DOCK7 axis essential for spatially restricting RAC1 activity in a planar polarized fashion in model tissue monolayers. The MYO6-DOCK7 axis specifically controls the extension of cryptic lamellipodia required to drive tissue fluidification and cooperative-mode motion in otherwise solid and static carcinoma cell collectives.

Alternative splicing of its 5'-UTR limits CD20 mRNA translation and enables resistance to CD20-directed immunotherapies

Aberrant skipping of coding exons in CD19 and CD22 compromises responses to immunotherapy for B-cell malignancies. Here, we show that the MS4A1 gene encoding human CD20 also produces several mRNA isoforms with distinct 5' untranslated regions (5'-UTR). Four variants (V1-4) were detectable by RNA-seq in distinct stages of normal B-cell differentiation and B-lymphoid malignancies, with V1 and V3 being the most abundant by far. During B-cell activation and Epstein-Barr virus infection, redirection of splicing from V1 to V3 coincided with increased CD20 positivity. Similarly, in diffuse large B-cell lymphoma only V3, but not V1, correlated with CD20 protein levels, suggesting that V1 might be translation-deficient. Indeed, the longer V1 isoform was found to contain upstream open reading frames (uORFs) and a stem-loop structure, which cooperatively inhibited polysome recruitment. By modulating CD20 isoforms with splice-switching Morpholino oligomers, we enhanced CD20 expression and anti-CD20 antibody rituximab-mediated cytotoxicity in a panel of B-cell lines. Furthermore, reconstitution of CD20-knockout cells with V3 mRNA led to the recovery of CD20 positivity, while V1-reconstituted cells had undetectable levels of CD20 protein. Surprisingly, in vitro CD20-directed CAR T cells were able to kill both V3- and V1-expressing cells, but the bispecific T cell engager mosunetuzumab was only effective against V3-expressing cells. To determine whether CD20 splicing is involved in immunotherapy resistance, we performed RNA-seq on four post-mosunetuzumab follicular lymphoma relapses and discovered that in two of them downregulation of CD20 was accompanied by the V3-to-V1 shift. Thus, splicing-mediated mechanisms of epitope loss extend to CD20-directed immunotherapies.

Key Points

In normal & malignant human B cells, CD20 mRNA is alternatively spliced into four 5'-UTR isoforms, some of which are translation-deficient.The balance between translation-deficient and -competent isoforms modulates CD20 protein levels & responses to CD20-directed immunotherapies.

Explanation of Novelty

We discovered that in normal and malignant B-cells, CD20 mRNA is alternatively spliced to generate four distinct 5'-UTRs, including the longer translation-deficient V1 variant. Cells predominantly expressing V1 were still sensitive to CD20-targeting chimeric antigen receptor T-cells. However, they were resistant to the bispecific anti-CD3/CD20 antibody mosunetuzumab, and the shift to V1 were observed in CD20-negative post-mosunetuzumab relapses of follicular lymphoma.

Dynamics of TERT regulation via alternative splicing in stem cells and cancer cells

Part of the regulation of telomerase activity includes the alternative splicing (AS) of the catalytic subunit telomerase reverse transcriptase (TERT). Although a therapeutic window for telomerase/TERT inhibition exists between cancer cells and somatic cells, stem cells express TERT and rely on telomerase activity for physiological replacement of cells. Therefore, identifying differences in TERT regulation between stem cells and cancer cells is essential for developing telomerase inhibition-based cancer therapies that reduce damage to stem cells. In this study, we measured TERT splice variant expression and telomerase activity in induced pluripotent stem cells (iPSCs), neural progenitor cells (NPCs), and non-small cell lung cancer cells (NSCLC, Calu-6 cells). We observed that a NOVA1-PTBP1-PTBP2 axis regulates TERT alternative splicing (AS) in iPSCs and their differentiation into NPCs. We also found that splice-switching of TERT, which regulates telomerase activity, is induced by different cell densities in stem cells but not cancer cells. Lastly, we identified cell type-specific splicing factors that regulate TERT AS. Overall, our findings represent an important step forward in understanding the regulation of TERT AS in stem cells and cancer cells.

Targeting CD44 Variant 5 with an Antibody-Drug Conjugate Is an Effective Therapeutic Strategy for Intrahepatic Cholangiocarcinoma

Intrahepatic cholangiocarcinoma (ICC) is the second most frequent type of primary liver cancer. ICC is among the deadliest malignancies, highlighting that novel treatments are urgently needed. Studies have shown that CD44 variant isoforms, rather than the CD44 standard isoform, are selectively expressed in ICC cells, providing an opportunity for the development of an antibody-drug conjugate (ADC)-based targeted therapeutic strategy. In this study, we observed the specific expression of CD44 variant 5 (CD44v5) in ICC tumors. CD44v5 protein was expressed on the surface of most ICC tumors (103 of 155). A CD44v5-targeted ADC, H1D8-DC (H1D8-drug conjugate), was developed that comprises a humanized anti-CD44v5 mAb conjugated to the microtubule inhibitor monomethyl auristatin E (MMAE) via a cleavable valine-citrulline-based linker. H1D8-DC exhibited efficient antigen binding and internalization in cells expressing CD44v5 on the cell surface. Because of the high expression of cathepsin B in ICC cells, the drug was preferentially released in cancer cells but not in normal cells, thus inducing potent cytotoxicity at picomolar concentrations. In vivo studies showed that H1D8-DC was effective against CD44v5-positive ICC cells and induced tumor regression in patient-derived xenograft models, whereas no significant adverse toxicities were observed. These data demonstrate that CD44v5 is a bona fide target in ICC and provide a rationale for the clinical investigation of a CD44v5-targeted ADC-based approach.

SIGNIFICANCE

Elevated expression of CD44 variant 5 in intrahepatic cholangiocarcinoma confers a targetable vulnerability using the newly developed antibody-drug conjugate H1D8-DC, which induces potent growth suppressive effects without significant toxicity.

NFX1-123: A potential therapeutic target in cervical cancer

NFX1-123 is a splice variant isoform of the NFX1 gene. It is highly expressed in cervical cancers caused by HPV, and NFX1-123 is a protein partner with the HPV oncoprotein E6. Together, NFX1-123 and E6 affect cellular growth, longevity, and differentiation. The expression status of NFX1-123 in cancers beyond cervical and head and neck cancers, and its potential as therapeutic target, have not been investigated. TSVdb of TCGA was used to quantify NFX1-123 expression in 24 cancers compared with normal tissues. The NFX1-123 protein structure was predicted and then submitted to retrieve suitable drug molecules. The top four compounds, found to bind in silico to NFX1-123, were tested experimentally to determine their effects on NFX1-123-related cellular growth, survival, and migration. 46% of cancers (11 of 24 had significant differences in NFX1-123 expression, with nine having had greater NFX1-123 expression, when compared with adjacent normal tissues. Bioinformatics and proteomic predictive analysis modeled the three-dimensional structure of NFX1-123, and drug libraries were screened for high-binding affinity compounds using this modeled structure. Seventeen drugs with binding energies ranging from -1.3 to -10 Kcal/mol were identified. The top four compounds were used to treat HPV- and HPV+ cervical cancer cell lines, three of which (Ropitoin, R428 and Ketoconazole) reduced NFX1-123 protein levels, inhibited cellular growth, survival, and migration, and enhanced the cytotoxicity of Cisplatin. These findings highlight cancers expressing high levels of NFX1-123, and drugs that target it, may reduce cellular growth, survival, and migration, making NFX1-123 a potential novel therapeutic target.

Alternative Splicing Changes Promoted by NOVA2 Upregulation in Endothelial Cells and Relevance for Gastric Cancer

Angiogenesis is crucial for cancer progression. While several anti-angiogenic drugs are in use for cancer treatment, their clinical benefits are unsatisfactory. Thus, a deeper understanding of the mechanisms sustaining cancer vessel growth is fundamental to identify novel biomarkers and therapeutic targets. Alternative splicing (AS) is an essential modifier of human proteome diversity. Nevertheless, AS contribution to tumor vasculature development is poorly known. The Neuro-Oncological Ventral Antigen 2 (NOVA2) is a critical AS regulator of angiogenesis and vascular development. NOVA2 is upregulated in tumor endothelial cells (ECs) of different cancers, thus representing a potential driver of tumor blood vessel aberrancies. Here, we identified novel AS transcripts generated upon NOVA2 upregulation in ECs, suggesting a pervasive role of NOVA2 in vascular biology. In addition, we report that NOVA2 is also upregulated in ECs of gastric cancer (GC), and its expression correlates with poor overall survival of GC patients. Finally, we found that the AS of the Rap Guanine Nucleotide Exchange Factor 6 (RapGEF6), a newly identified NOVA2 target, is altered in GC patients and associated with NOVA2 expression, tumor angiogenesis, and poor patient outcome. Our findings provide a better understanding of GC biology and suggest that AS might be exploited to identify novel biomarkers and therapeutics for anti-angiogenic GC treatments.

Interactions of circadian clock genes with the hallmarks of cancer

The molecular machinery of the circadian clock regulates the expression of many genes and processes in the organism, allowing the adaptation of cellular activities to the daily light-dark cycles. Disruption of the circadian rhythm can lead to various pathologies, including cancer. Thus, disturbance of the normal circadian clock at both genetic and environmental levels has been described as an independent risk factor for cancer. In addition, researchers have proposed that circadian genes may have a tissue-dependent and/or context-dependent role in tumorigenesis and may function both as tumor suppressors and oncogenes. Finally, circadian clock core genes may trigger or at least be involved in different hallmarks of cancer. Hence, expanding the knowledge of the molecular basis of the circadian clock would be helpful to identify new prognostic markers of tumorigenesis and potential therapeutic targets.

Matrix stiffness regulates tumor cell intravasation through expression and ESRP1-mediated alternative splicing of MENA

During intravasation, cancer cells cross the endothelial barrier and enter the circulation. Extracellular matrix stiffening has been correlated with tumor metastatic potential; however, little is known about the effects of matrix stiffness on intravasation. Here, we utilize in vitro systems, a mouse model, specimens from patients with breast cancer, and RNA expression profiles from The Cancer Genome Atlas Program (TCGA) to investigate the molecular mechanism by which matrix stiffening promotes tumor cell intravasation. Our data show that heightened matrix stiffness increases MENA expression, which promotes contractility and intravasation through focal adhesion kinase activity. Further, matrix stiffening decreases epithelial splicing regulatory protein 1 (ESRP1) expression, which triggers alternative splicing of MENA, decreases the expression of MENA11a, and enhances contractility and intravasation. Altogether, our data indicate that matrix stiffness regulates tumor cell intravasation through enhanced expression and ESRP1-mediated alternative splicing of MENA, providing a mechanism by which matrix stiffness regulates tumor cell intravasation.

Meta-analysis of Osteopontin splice variants in cancer

BACKGROUND

The cytokine Osteopontin is a mediator of tumor progression and cancer metastasis. In 2006, we reported that (in addition to the full-length form -a) splice variants of Osteopontin (forms -b and -c) are produced selectively by transformed cells. Through June 2021, 36 PubMed-indexed journal articles have studied Osteopontin splice variants in various cancer patients.

METHODS

Applying a categorical approach previously developed by us, here we conduct a meta-analysis of the pertinent literature. We supplement this with evaluation of the relevant entries in the TSVdb database, which focusses on splice variant expression, thus including the additional variants -4 and -5. The analysis covers 5886 patients across 15 tumors from the literature and 10,446 patients across 33 tumors from TSVdb.

RESULTS

The database yields positive results more frequently than the categorical meta-analysis. The two sources are in agreement on the elevation of OPN-a, OPN-b, and OPN-c in lung cancer and the elevation of OPN-c in breast cancer as compared to healthy tissue. Specific splice variants are associated with grade, stage, or patient survival pertaining to various cancers.

CONCLUSIONS

There are cases of persisting discrepancies, which require further investigation to clarify the Osteopontin splice variant utilization, so that their diagnostic, prognostic and potentially predictive potential can be brought to fruition.

The small GTPase ARF3 controls invasion modality and metastasis by regulating N-cadherin levels

ARF GTPases are central regulators of membrane trafficking that control local membrane identity and remodeling facilitating vesicle formation. Unraveling their function is complicated by the overlapping association of ARFs with guanine nucleotide exchange factors (GEFs), GTPase-activating proteins (GAPs), and numerous interactors. Through a functional genomic screen of three-dimensional (3D) prostate cancer cell behavior, we explore the contribution of ARF GTPases, GEFs, GAPs, and interactors to collective invasion. This revealed that ARF3 GTPase regulates the modality of invasion, acting as a switch between leader cell-led chains of invasion or collective sheet movement. Functionally, the ability of ARF3 to control invasion modality is dependent on association and subsequent control of turnover of N-cadherin. In vivo, ARF3 levels acted as a rheostat for metastasis from intraprostatic tumor transplants and ARF3/N-cadherin expression can be used to identify prostate cancer patients with metastatic, poor-outcome disease. Our analysis defines a unique function for the ARF3 GTPase in controlling how cells collectively organize during invasion and metastasis.

PredDSMC: A predictor for driver synonymous mutations in human cancers

Introduction: Driver mutations play a critical role in the occurrence and development of human cancers. Most studies have focused on missense mutations that function as drivers in cancer. However, accumulating experimental evidence indicates that synonymous mutations can also act as driver mutations.

Methods: Here, we proposed a computational method called PredDSMC to accurately predict driver synonymous mutations in human cancers. We first systematically explored four categories of multimodal features, including sequence features, splicing features, conservation scores, and functional scores. Further feature selection was carried out to remove redundant features and improve the model performance. Finally, we utilized the random forest classifier to build PredDSMC.

Results: The results of two independent test sets indicated that PredDSMC outperformed the state-of-the-art methods in differentiating driver synonymous mutations from passenger mutations.

Discussion: In conclusion, we expect that PredDSMC, as a driver synonymous mutation prediction method, will be a valuable method for gaining a deeper understanding of synonymous mutations in human cancers.

RAC1B function is essential for breast cancer stem cell maintenance and chemoresistance of breast tumor cells

Breast cancer stem cells (BCSC) are presumed to be responsible for treatment resistance, tumor recurrence and metastasis of breast tumors. However, development of BCSC-targeting therapies has been held back by their heterogeneity and the lack of BCSC-selective molecular targets. Here, we demonstrate that RAC1B, the only known alternatively spliced variant of the small GTPase RAC1, is expressed in a subset of BCSCs in vivo and its function is required for the maintenance of BCSCs and their chemoresistance to doxorubicin. In human breast cancer cell line MCF7, RAC1B is required for BCSC plasticity and chemoresistance to doxorubicin in vitro and for tumor-initiating abilities in vivo. Unlike Rac1, Rac1b function is dispensable for normal mammary gland development and mammary epithelial stem cell (MaSC) activity. In contrast, loss of Rac1b function in a mouse model of breast cancer hampers the BCSC activity and increases their chemosensitivity to doxorubicin treatment. Collectively, our data suggest that RAC1B is a clinically relevant molecular target for the development of BCSC-targeting therapies that may improve the effectiveness of doxorubicin-mediated chemotherapy.

Labyrinthin Expression Is Associated with Poor Prognosis in Patients with Non-Small-Cell Lung Cancer

To determine Labyrinthin (LAB) expression in non-small-cell lung cancer (NSCLC), we immunostained and scored for LAB immunohistochemistry (IHC) expression on sections of tissue microarrays (TMAs) prepared from 256 archival tissue blocks of NSCLC. Propensity-score-weighted Kaplan-Meier curves and weighted Cox models were used to associate LAB expression with overall survival. LAB mRNA expression was assessed in The Cancer Genome Atlas (TCGA) and correlated with clinical phenotype and outcome. Positive LAB IHC expression (>5% of tumor cells) was detected in 208/256 (81.3%) of NSCLC samples, and found in both lung adenocarcinomas (LUAD) and lung squamous cell cancer (LUSC). LAB positivity was associated with poor overall survival (HR = 3.56, 95% CI: 2.3-5.4; p < 0.0001) and high tumor differentiation grade or metastasis compared with negative LAB expression. Univariant and multivariate survival analyses demonstrated LAB expression as an independent prognostic factor for NSCLC patients. LAB RNA expression in TCGA-LUAD was higher in primary and advanced-stage tumors than in normal tissue, and was associated with poorer overall survival. No significant differences or associations were found with LAB RNA expression in TCGA-LUSC. The LAB IHC assay is being used to identify candidate cancer patients for the first-in-human phase I trial evaluating the LAB vaccines (UCDCC#296, NCT051013560).

FOXP3 Isoforms Expression in Cervical Cancer: Evidence about the Cancer-Related Properties of FOXP3Δ2Δ7 in Keratinocytes

Cervical cancer (CC) is the fourth most common type of cancer among women; the main predisposing factor is persistent infection by high-risk human papillomavirus (hr-HPV), mainly the 16 or 18 genotypes. Both hr-HPVs are known to manipulate the cellular machinery and the immune system to favor cell transformation. FOXP3, a critical transcription factor involved in the biology of regulatory T cells, has been detected as highly expressed in the tumor cells of CC patients. However, its biological role in CC, particularly in the keratinocytes, remained unclarified. Therefore, this work aimed to uncover the effect of FOXP3 on the biology of the tumoral cells. First, public databases were analyzed to identify the FOXP3 expression levels and the transcribed isoforms in CC and normal tissue samples. The study's findings demonstrated an increased expression of FOXP3 in HPV16+ CC samples. Additionally, the FOXP3Δ2 variant was detected as the most frequent splicing isoform in tumoral cells, with a high differential expression level in metastatic samples. However, the analysis of FOXP3 expression in different CC cell lines, HPV+ and HPV-, suggests no relationship between the presence of HPV and FOXP3 expression. Since the variant FOXP3Δ2Δ7 was found highly expressed in the HPV16+ SiHa cell line, a model with constitutive expression of FOXP3Δ2Δ7 was established to evaluate its role in proliferation, migration, and cell division. Finally, RNAseq was performed to identify differentially expressed genes and enriched pathways modulated by FOXP3Δ2Δ7. The exogenous expression of FOXP3Δ2Δ7 promotes cell division, proliferation, and migration. The transcriptomic analyses highlight the upregulation of multiple genes with protumor activities. Moreover, immunological and oncogenic pathways were detected as highly enriched. These data support the hypothesis that FOXP3Δ2Δ7 in epithelial cells induces cancer-related hallmarks and provides information about the molecular events triggered by this isoform, which could be important for developing CC.

N6-Methyladenosine-modified lncRNA LINREP promotes Glioblastoma progression by recruiting the PTBP1/HuR complex

Glioblastoma multiforme (GBM) is acknowledged as the most aggressive primary brain tumor in adults. It is typically characterized by the high heterogeneity which corresponds to extensive genetic mutations and complex alternative splicing (AS) profiles. Known as a major repressive splicing factor in AS, polypyrimidine tract-binding protein 1 (PTBP1) is involved in the exon skipping events of multiple precursor mRNAs (pre-mRNAs) in GBM. However, precise mechanisms that modulate the expression and activity of PTBP1 remain to be elucidated. In present study, we provided evidences for the role of a long intergenic noncoding RNA (LINREP) implicated in the regulation of PTBP1-induced AS. LINREP interacted with PTBP1 and human antigen R (HuR, ELAVL1) protein complex and protected PTBP1 from the ubiquitin-proteasome degradation. Consequently, a broad spectrum of PTBP1-induced spliced variants was generated by exon skipping, especially for the skipping of reticulon 4 (RTN4) exon 3. Interestingly, LINREP also promoted the dissociation of nuclear UPF1 from PTBP1, which increased the binding of PTBP1 to RTN4 transcripts, thus enhancing the skipping of RTN4 exon 3 to some extent. Besides, HuR recruitment was essential for the stabilization of LINREP via a manner dependent on N6-methyladenosine (m6A) formation and identification. Taken together, our results demonstrated the functional significance of LINREP in human GBM for its dual regulation of PTBP1-induced AS and its m6A modification modality, implicating that HuR/LINREP/PTBP1 axis might serve as a potential therapeutic target for GBM.

Establishment of quantitative nested-PCR of Abelson interactor 1 transcript variant-11

Abelson interactor 1 (ABI1), which presents 18 Transcript Variants (TSV), plays an important role in CRC metastasis. Different ABI1-TSVs play synergistic or antagonistic roles in the same pathophysiological events. ABI1 Transcript Variant-11 (ABI1-TSV-11) functionally promotes lymph node metastasis of left-sided colorectal cancer (LsCC) and is an independent molecular marker to evaluate the prognosis of patients with LsCC. However, there is still lack of a quick and accurate method to detect the expression of ABI-TSV-11, distinguishing ABI1-TSV-11 from other 17 TSVs. To establish a rapid method specific for ABI1-TSV-11detection, we developed a quantitative nested-PCR method composed of pre-amplification regular PCR using ABI1 universal primer pair and the followed Real Time (RT)-qPCR using ABI1-TSV-11 specific primer pair spanning exon-exon junction. ABI1-TSV-11-overexpressed SW480 and LoVo cell lines were used to verify the quantitative nested-PCR assay, and the sequencing data was used to evaluate the accuracy of ABI1-TSV-11 quantitative nested-PCR assay. The detection limit was 5.24×10⁴ copies/ml. ABI1-TSV-11 quantitative nested-PCR provides a new technical means for the detection of ABI1-TSV-11.

Prognostic value of PNN in prostate cancer and its correlation with therapeutic significance

Prostate cancer (PCa) is the most common malignancy. New biomarkers are in demand to facilitate the management. The role of the pinin protein (encoded by PNN gene) in PCa has not been thoroughly explored yet. Using The Cancer Genome Atlas (TCGA-PCa) dataset validated with Gene Expression Omnibus (GEO) and protein expression data retrieved from the Human Protein Atlas, the prognostic and diagnostic values of PNN were studied. Highly co-expressed genes with PNN (HCEG) were constructed for pathway enrichment analysis and drug prediction. A prognostic signature based on methylation status using HCEG was constructed. Gene set enrichment analysis (GSEA) and the TISIDB database were utilised to analyse the associations between PNN and tumour-infiltrating immune cells. The upregulated PNN expression in PCa at both transcription and protein levels suggests its potential as an independent prognostic factor of PCa. Analyses of the PNN's co-expression network indicated that PNN plays a role in RNA splicing and spliceosomes. The prognostic methylation signature demonstrated good performance for progression-free survival. Finally, our results showed that the PNN gene was involved in splicing-related pathways in PCa and identified as a potential biomarker for PCa.

Single-molecule array assay reveals the prognostic impact of plasma LRIG1 in ovarian carcinoma

BACKGROUND

Ovarian carcinoma is the eighth most common cause of cancer death in women worldwide. The disease is predominantly diagnosed at a late stage. This contributes to high recurrence rates, eventually leading to the development of treatment-resistant disease. Leucine-rich repeats and immunoglobulin-like domains protein 1 (LRIG1) is a transmembrane protein that functions as a tumor suppressor and regulator of growth factor signaling. LRIG1 levels have not been investigated in human plasma previously.

MATERIALS AND METHODS

A quantitative LRIG1-specific single molecule array assay was developed and validated. LRIG1 levels were quantified in plasma samples from 486 patients with suspicious ovarian masses.

RESULTS

Among women with ovarian carcinoma, LRIG1 levels were significantly elevated compared to women with benign or borderline type tumors. High LRIG1 plasma levels were associated with worse overall survival and shorter disease-free survival both in the group of all malignant cases and among the stage 3 cases only. LRIG1 was an independent prognostic factor in patients with stage 3 ovarian carcinoma.

CONCLUSION

LRIG1 plasma levels were elevated in patients with ovarian carcinoma, and high levels were associated with poor prognosis, suggesting that LRIG1 might be an etiologic factor and a potentially useful biomarker in ovarian carcinoma.

Intronic Cis-Element DR8 in hTERT Is Bound by Splicing Factor SF3B4 and Regulates hTERT Splicing in Non-Small Cell Lung Cancer

Splicing of the hTERT gene to produce the full-length (FL) transcript is necessary for telomerase enzyme activity and telomere-dependent cellular immortality in the majority of human tumors, including non-small cell lung cancer (NSCLC) cells. The molecular machinery to splice hTERT to the FL isoform remains mostly unknown. Previously, we reported that an intron 8 cis-element termed "direct repeat 8" (DR8) promotes FL hTERT splicing, telomerase, and telomere length maintenance when bound by NOVA1 and PTBP1 in NSCLC cells. However, some NSCLC cells and patient tumor samples lack NOVA1 expression. This leaves a gap in knowledge about the splicing factors and cis-elements that promote telomerase in the NOVA1-negative context. We report that DR8 regulates FL hTERT splicing in the NOVA1-negative and -positive lung cancer contexts. We identified splicing factor 3b subunit 4 (SF3B4) as an RNA trans-factor whose expression is increased in lung adenocarcinoma (LUAD) tumors compared with adjacent normal tissue and predicts poor LUAD patient survival. In contrast to normal lung epithelial cells, which continued to grow with partial reductions of SF3B4 protein, SF3B4 knockdown reduced hTERT splicing, telomerase activity, telomere length, and cell growth in lung cancer cells. SF3B4 was also demonstrated to bind the DR8 region of hTERT pre-mRNA in both NOVA1-negative and -positive NSCLC cells. These findings provide evidence that DR8 is a critical binding hub for trans-factors to regulate FL hTERT splicing in NSCLC cells. These studies help define mechanisms of gene regulation important to the generation of telomerase activity during carcinogenesis.

IMPLICATIONS

Manipulation of a core spliceosome protein reduces telomerase/hTERT splicing in lung cancer cells and results in slowed cancer cell growth and cell death, revealing a potential therapeutic strategy.

Proteogenomic, Epigenetic, and Clinical Implications of Recurrent Aberrant Splice Variants in Clear Cell Renal Cell Carcinoma

BACKGROUND

Alternative mRNA splicing can be dysregulated in cancer, resulting in the generation of aberrant splice variants (SVs). Given the paucity of actionable genomic mutations in clear cell renal cell carcinoma (ccRCC), aberrant SVs may be an avenue to novel mechanisms of pathogenesis.

OBJECTIVE

To identify and characterize aberrant SVs enriched in ccRCC.

DESIGN, SETTING, AND PARTICIPANTS

Using RNA-seq data from the Cancer Cell Line Encyclopedia, we identified neojunctions uniquely expressed in ccRCC. Candidate SVs were then checked for expression across normal tissue in the Genotype-Tissue Expression Project and primary tumor tissue from The Cancer Genome Atlas (TCGA), Clinical Proteomic Tumor Analysis Consortium (CPTAC), and our institutional Total Cancer Care database.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

Clinicopathologic, genomic, and survival data were available for all cohorts. Epigenetic data were available for the TCGA and CPTAC cohorts. Proteomic data were available for the CPTAC cohort. The association of aberrant SV expression with these variables was examined using the Kruskal-Wallis test, pairwise t test, Spearman correlation test, and Cox regression analysis.

RESULTS AND LIMITATIONS

Our pipeline identified 16 ccRCC-enriched SVs. EGFR, HPCAL1-SV and RNASET2-SV expression was negatively correlated with gene-specific CpG methylation. We derived a survival risk score based primarily on the expression of five SVs (RNASET2, FGD1, PDZD2, COBLL1, and PTPN14), which was consistent and applicable across multiple cohorts on multivariate analysis. The splicing factor RBM4, which modulates splicing of HIF-1α, exhibited significantly lower expression at the protein level in the high-risk group, as defined by our SV-based score.

CONCLUSIONS

We describe 16 aberrant SVs enriched in ccRCC, many of which are associated with disease biology and/or clinical outcomes. This study provides a novel strategy for identifying and characterizing disease-specific aberrant SVs.

PATIENT SUMMARY

We describe a method to identify disease targets and biomarkers using transcriptomic analysis beyond somatic mutations or gene expression. Kidney tumors express unique splice variants that may provide additional prognostic information following surgery.

Paracrine signalling between intestinal epithelial and tumour cells induces a regenerative programme

Tumours are complex ecosystems composed of different types of cells that communicate and influence each other. While the critical role of stromal cells in affecting tumour growth is well established, the impact of mutant cancer cells on healthy surrounding tissues remains poorly defined. Here, using mouse intestinal organoids, we uncover a paracrine mechanism by which intestinal cancer cells reactivate foetal and regenerative YAP-associated transcriptional programmes in neighbouring wildtype epithelial cells, rendering them adapted to thrive in the tumour context. We identify the glycoprotein thrombospondin-1 (THBS1) as the essential factor that mediates non-cell-autonomous morphological and transcriptional responses. Importantly, Thbs1 is associated with bad prognosis in several human cancers. This study reveals the THBS1-YAP axis as the mechanistic link mediating paracrine interactions between epithelial cells in intestinal tumours.

Mitochondrial-nuclear epistasis underlying phenotypic variation in breast cancer pathology

The interplay between genes harboring single nucleotide polymorphisms (SNPs) is vital to better understand underlying contributions to the etiology of breast cancer. Much attention has been paid to epistasis between nuclear genes or mutations in the mitochondrial genome. However, there is limited understanding about the epistatic effects of genetic variants in the nuclear and mitochondrial genomes jointly on breast cancer. We tested the interaction of germline SNPs in the mitochondrial (mtSNPs) and nuclear (nuSNPs) genomes of female breast cancer patients in The Cancer Genome Atlas (TCGA) for association with morphological features extracted from hematoxylin and eosin (H&E)-stained pathology images. We identified 115 significant (q-value < 0.05) mito-nuclear interactions that increased nuclei size by as much as 12%. One interaction between nuSNP rs17320521 in an intron of the WSC Domain Containing 2 (WSCD2) gene and mtSNP rs869096886, a synonymous variant mapped to the mitochondrially-encoded NADH dehydrogenase 4 (MT-ND4) gene, was confirmed in an independent breast cancer data set from the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC). None of the 10 mito-nuclear interactions identified from non-diseased female breast tissues from the Genotype-Expression (GTEx) project resulted in an increase in nuclei size. Comparisons of gene expression data from the TCGA breast cancer patients with the genotype homozygous for the minor alleles of the SNPs in WSCD2 and MT-ND4 versus the other genotypes revealed core transcriptional regulator interactions and an association with insulin. Finally, a Cox proportional hazards ratio = 1.7 (C.I. 0.98-2.9, p-value = 0.042) and Kaplan-Meier plot suggest that the TCGA female breast cancer patients with low gene expression of WSCD2 coupled with large nuclei have an increased risk of mortality. The intergenomic dependency between the two variants may constitute an inherent susceptibility of a more severe form of breast cancer and points to genetic targets for further investigation of additional determinants of the disease.

GFAP splice variants fine-tune glioma cell invasion and tumour dynamics by modulating migration persistence

Glioma is the most common form of malignant primary brain tumours in adults. Their highly invasive nature makes the disease incurable to date, emphasizing the importance of better understanding the mechanisms driving glioma invasion. Glial fibrillary acidic protein (GFAP) is an intermediate filament protein that is characteristic for astrocyte- and neural stem cell-derived gliomas. Glioma malignancy is associated with changes in GFAP alternative splicing, as the canonical isoform GFAPα is downregulated in higher-grade tumours, leading to increased dominance of the GFAPδ isoform in the network. In this study, we used intravital imaging and an ex vivo brain slice invasion model. We show that the GFAPδ and GFAPα isoforms differentially regulate the tumour dynamics of glioma cells. Depletion of either isoform increases the migratory capacity of glioma cells. Remarkably, GFAPδ-depleted cells migrate randomly through the brain tissue, whereas GFAPα-depleted cells show a directionally persistent invasion into the brain parenchyma. This study shows that distinct compositions of the GFAPnetwork lead to specific migratory dynamics and behaviours of gliomas.

OUP accepted manuscript

NDST1 (glucosaminyl N-deacetylase/N-sulfotransferase) is a key enzyme in heparan sulfate (HS) biosynthesis, where it is responsible for HS N-deacetylation and N-sulfation. In addition to the full length human enzyme of 882 amino acids, here designated NDST1A, a shorter form containing 825 amino acids (NDST1B) is synthesized after alternative splicing of the NDST1 mRNA. NDST1B is mostly expressed at a low level, but increased amounts are seen in several types of cancer where it is associated with shorter survival. In this study, we aimed at characterizing the enzymatic properties of NDST1B and its effect on HS biosynthesis. Purified recombinant NDST1B lacked both N-deacetylase and N-sulfotransferase activities. Interestingly, HEK293 cells overexpressing NDST1B synthesized HS with reduced sulfation and altered domain structure. Fluorescence resonance energy transfer-microscopy demonstrated that both NDST1A and NDST1B had the capacity to interact with the HS copolymerase subunits EXT1 and EXT2 and also to form NDST1A/NDST1B dimers. Since lysates from cells overexpressing NDST1B contained less NDST enzyme activity than control cells, we suggest that NDST1B works in a dominant negative manner, tentatively by replacing the active endogenous NDST1 in the enzyme complexes taking part in biosynthesis.

High expression of NFX1-123 in HPV positive head and neck squamous cell carcinomas

BACKGROUND

High-risk human papillomaviruses (HR HPV) cause nearly all cervical cancers and, in the United States, the majority of head and neck cancers (HNSCCs). NFX1-123 is overexpressed in cervical cancers, and NFX1-123 partners with the HR HPV type 16 E6 oncoprotein to affect multiple growth, differentiation, and immune response genes. However, neither the expression of NFX1-123 nor the levels of these genes have been investigated in HPV positive (HPV+) or negative (HPV-) HNSCCs.

METHODS

The Cancer Genome Atlas Splicing Variants Database and HNSCC cell lines were used to quantify expression of NFX1-123 and cellular genes increased in cervical cancers.

RESULTS

NFX1-123 was increased in HPV+ HNSCCs compared to HPV- HNSCCs. LCE1B, KRT16, SPRR2G, and FBN2 were highly expressed in HNSCCs compared to normal tissues. Notch1 and CCNB1IP1 had greater expression in HPV+ HNSCCs compared to HPV- HNSCCs.

CONCLUSION

NFX1-123 and a subset of its known targets were increased in HPV+ HNSCCs.

Systematic Genome-Wide Profiles Reveal Alternative Splicing Landscape and Implications of Splicing Regulator DExD-Box Helicase 21 in Aggressive Progression of Adrenocortical Carcinoma

Alternative splicing (AS) in the tumor biological process has provided a novel perspective on carcinogenesis. However, the clinical significance of individual AS patterns of adrenocortical carcinoma (ACC) has been underestimated, and in-depth investigations are lacking. We selected 76 ACC samples from the Cancer Genome Atlas (TCGA) SpliceSeq and SpliceAid2 databases, and 39 ACC samples from Fudan University Shanghai Cancer Center (FUSCC). Prognosis-related AS events (PASEs) and survival analysis were evaluated based on prediction models constructed by machine-learning algorithm. In total, 23,984 AS events and 3,614 PASEs were detected in the patients with ACC. The predicted risk score of each patient suggested that eight PASEs groups were significantly correlated with the clinical outcomes of these patients (p < 0.001). Prognostic models produced AUC values of 0.907 in all PASEs' groups. Eight splicing factors (SFs), including BAG2, CXorf56, DExD-Box Helicase 21 (DDX21), HSPB1, MBNL3, MSI1, RBMXL2, and SEC31B, were identified in regulatory networks of ACC. DDX21 was identified and validated as a novel clinical promoter and therapeutic target in 115 patients with ACC from TCGA and FUSCC cohorts. In conclusion, the strict standards used in this study ensured the systematic discovery of profiles of AS events using genome-wide cohorts. Our findings contribute to a comprehensive understanding of the landscape and underlying mechanism of AS, providing valuable insights into the potential usages of DDX21 for predicting prognosis for patients with ACC.

Supplementary Information

The online version contains supplementary material available at 10.1007/s43657-021-00026-x.

Cyclin E1 in Murine and Human Liver Cancer: A Promising Target for Therapeutic Intervention during Tumour Progression

Cyclin E1 (CCNE1) is a regulatory subunit of Cyclin-dependent kinase 2 (CDK2) and is thought to control the transition of quiescent cells into cell cycle progression. Recently, we identified CCNE1 and CDK2 as key factors for the initiation of hepatocellular carcinoma (HCC). In the present study, we dissected the contributions of CCNE1 and CDK2 for HCC progression in mice and patients. Therefore, we generated genetically modified mice allowing inducible deletion of Ccne1 or Cdk2. After initiation of HCC, using the hepatocarcinogen diethylnitrosamine (DEN), we deleted Ccne1 or Cdk2 and subsequently analysed HCC progression. The relevance of CCNE1 or CDK2 for human HCC progression was investigated by in silico database analysis. Interventional deletion of Ccne1, but not of Cdk2, substantially reduced the HCC burden in mice. Ccne1-deficient HCCs were characterised by attenuated proliferation, impaired DNA damage response and downregulation of markers for stemness and microinvasion. Additionally, the tumour microenvironment of Ccne1-deficient mice showed a reduction in immune mediators, myeloid cells and cancer-associated fibroblasts. In sharp contrast, Cdk2 was dispensable for HCC progression in mice. In agreement with our mouse data, CCNE1 was overexpressed in HCC patients independent of risk factors, and associated with reduced disease-free survival, a common signature for enhanced chromosomal instability, proliferation, dedifferentiation and invasion. However, CDK2 lacked diagnostic or prognostic value in HCC patients. In summary, CCNE1 drives HCC progression in a CDK2-independent manner in mice and man. Therefore, interventional inactivation of CCNE1 represents a promising strategy the treatment of liver cancer.

Characterization and prognostic value of LXR splice variants in triple-negative breast cancer

Activity of liver x receptor (LXR), the homeostatic regulator of cholesterol metabolism, is elevated in triple-negative breast cancer (BCa) relative to other BCa subtypes, driving drug resistance and metastatic gene signatures. The loci encoding LXRα and LXRβ produce multiple alternatively spliced proteins, but the true range of variants and their relevance to cancer remain poorly defined. Here, we report seven LXR splice variants, three of which have not previously been reported and five that were prognostic for disease-free survival. Expression of full-length LXRα splice variants was associated with poor prognosis, consistent with a role as an oncogenic driver of triple-negative tumor pathophysiology. Contrary to this was the observation that high expression of truncated LXRα splice variants or any LXRβ splice variant was associated with longer survival. These findings indicate that LXR isoform abundance is an important aspect of understanding the link between dysregulated cholesterol metabolism and cancer pathophysiology.

The role of alternative splicing in human cancer progression

In eukaryotes, alternative splicing refers to a process via which a single precursor RNA (pre-RNA) is transcribed into different mature RNAs. Thus, alternative splicing enables the translation of a limited number of coding genes into a large number of proteins with different functions. Although, alternative splicing is common in normal cells, it also plays an important role in cancer development. Alteration in splicing mechanisms and even the participation of non-coding RNAs may cause changes in the splicing patterns of cancer-related genes. This article reviews the latest research on alternative splicing in cancer, with a view to presenting new strategies and guiding future studies related to pathological mechanisms associated with cancer.

A ligand-insensitive UNC5B splicing isoform regulates angiogenesis by promoting apoptosis

The Netrin-1 receptor UNC5B is an axon guidance regulator that is also expressed in endothelial cells (ECs), where it finely controls developmental and tumor angiogenesis. In the absence of Netrin-1, UNC5B induces apoptosis that is blocked upon Netrin-1 binding. Here, we identify an UNC5B splicing isoform (called UNC5B-Δ8) expressed exclusively by ECs and generated through exon skipping by NOVA2, an alternative splicing factor regulating vascular development. We show that UNC5B-Δ8 is a constitutively pro-apoptotic splicing isoform insensitive to Netrin-1 and required for specific blood vessel development in an apoptosis-dependent manner. Like NOVA2, UNC5B-Δ8 is aberrantly expressed in colon cancer vasculature where its expression correlates with tumor angiogenesis and poor patient outcome. Collectively, our data identify a mechanism controlling UNC5B’s necessary apoptotic function in ECs and suggest that the NOVA2/UNC5B circuit represents a post-transcriptional pathway regulating angiogenesis.

UNC5B is a Netrin-1 receptor expressed in endothelial cells that in the absence of ligand induces apoptosis. Here the authors identify an UNC5B splicing isoform that is insensitive to the pro-survival ligand Netrin-1 and is required for apoptosis-dependent blood vessel development.