Widespread alternative exon usage in clinically distinct subtypes of Invasive Ductal Carcinoma

MYO6 contributes to tumor progression and enzalutamide resistance in castration-resistant prostate cancer by activating the focal adhesion signaling pathway

BACKGROUND

Enzalutamide (Enz) resistance is a poor prognostic factor for patients with castration-resistant prostate cancer (CRPC), which often involves aberrant expression of the androgen receptor (AR). Myosin VI (MYO6), one member of the myosin family, plays an important role in regulating cell survival and is highly expressed in prostate cancer (PCa). However, whether MYO6 is involved in Enz resistance in CRPC and its mechanism remain unclear.

METHODS

Multiple open-access databases were utilized to examine the relationship between MYO6 expression and PCa progression, and to screen differentially expressed genes (DEGs) and potential signaling pathways associated with the MYO6-regulated Enz resistance. Both in vitro and in vivo tumorigenesis assays were employed to examine the impact of MYO6 on the growth and Enz resistance of PCa cells. Human PCa tissues and related clinical biochemical data were utilized to identify the role of MYO6 in promoting PCa progression and Enz resistance. The molecular mechanisms underlying the regulation of gene expression, PCa progression, and Enz resistance in CRPC by MYO6 were investigated.

RESULTS

MYO6 expression increases in patients with PCa and is positively correlated with AR expression in PCa cell lines and tissues. Overexpression of AR increases MYO6 expression to promote PCa cell proliferation, migration and invasion, and to inhibit PCa cell apoptosis; whereas knockdown of MYO6 expression reverses these outcomes and enhances Enz function in suppressing the proliferation of the Enz- sensitive and resistant PCa cells both in vitro and in vivo. Mechanistically, AR binds directly to the promoter region (residues - 503 to - 283 base pairs) of MYO6 gene and promotes its transcription. Furthermore, MYO6 activates focal adhesion kinase (FAK) phosphorylation at tyrosine-397 through integrin beta 8 (ITGB8) modulation to promote PCa progression and Enz resistance. Notably, inhibition of FAK activity by Y15, an inhibitor of FAK, can resensitize CRPC cells to Enz treatment in cell lines and mouse xenograft models.

CONCLUSIONS

MYO6 has pro-tumor and Enz-resistant effects in CRPC, suggesting that targeting MYO6 may be beneficial for ENZ-resistant CRPC therapy through the AR/MYO6/FAK signaling pathway.

Identification of a novel alternative splicing isoform of the Hippo kinase STK3/MST2 with impaired kinase and cell growth suppressing activities

Mammalian Ste-20-like Kinases 1 and 2 (MST1/2) are core serine-threonine kinases of the Hippo pathway regulating several cellular processes, including cell cycle arrest and cell death. Here, we discovered a novel alternative splicing variant of the MST2 encoding gene, STK3, in malignant cells and tumor datasets. This variant, named STK3∆7 or MST2∆7 (for mRNA or protein, respectively), resulted from the skipping of exon 7. MST2∆7 exhibited increased ubiquitylation and interaction with the E3 ubiquitin-protein ligase CHIP compared to the full-length protein (MST2FL). Exon 7 in STK3 encodes a segment within the kinase domain, and its exclusion compromised MST2 interaction with and phosphorylation of MOB, a major MST1/2 substrate. Nevertheless, MST2∆7 was capable of interacting with MST1 and MST2FL. Unlike MST2FL, overexpression of MST2∆7 did not lead to increased cell death and growth arrest. Strikingly, we observed the exclusion of STK3 exon 7 in 3.2-15% of tumor samples from patients of several types of cancer, while STK3∆7 was seldomly found in healthy tissues. Our study identified a novel STK3 splicing variant with loss of function and the potential to disturb tissue homeostasis by impacting on MST2 activities in the regulation of cell death and quiescence.

Expression of Truncated Products at the 5'-Terminal Region of RIPK2 and Evolutive Aspects that Support Their Biological Importance

Alternative splicing is the process of generating different mRNAs from the same primary transcript, which contributes to increase the transcriptome and proteome diversity. Abnormal splicing has been associated with the development of several diseases including cancer. Given that mutations and abnormal levels of the RIPK2 transcript and RIP-2 protein are frequent in tumors, and that RIP-2 modulates immune and inflammatory responses, we investigated alternative splicing events that result in partial deletions of the kinase domain at the N-terminus of RIP-2. We also investigated the structure and expression of the RIPK2 truncated variants and isoforms in different environments. In addition, we searched data throughout Supraprimates evolution that could support the biological importance of RIPK2 alternatively spliced products. We observed that human variants and isoforms were differentially regulated following temperature stress, and that the truncated transcript was more expressed than the long transcript in tumor samples. The inverse was found for the longer protein isoform. The truncated variant was also detected in chimpanzee, gorilla, hare, pika, mouse, rat, and tree shrew. The fact that the same variant has been preserved in mammals with divergence times up to 70 million years raises the hypothesis that it may have a functional significance.

Subtype-specific alternative splicing events in breast cancer identified by large-scale data analysis

Genome analysis in cancer has focused mainly on elucidating the function and regulatory mechanisms of genes that exhibit differential expression or mutation in cancer samples compared to normal samples. Recently, transcriptome analysis revealed that abnormal splicing events in cancer samples could contribute to cancer pathogenesis. Moreover, splicing variants in cancer reportedly generate diverse cancer antigens. Although abnormal splicing events are expected to be potential targets in cancer immunotherapy, the exploration of such targets and their biological significance in cancer have not been fully understood. In this study, to explore subtype-specific alternative splicing events, we conducted a comprehensive analysis of splicing events for each breast cancer subtype using large-scale splicing data derived from The Cancer Genome Atlas and found subtype-specific alternative splicing patterns. Analyses indicated that genes that produce subtype-specific alternative splicing events are potential novel targets for immunotherapy against breast cancer. The subtype-specific alternative splicing events identified in this study, which were not identified by mutation or differential expression analysis, bring new significance to previously overlooked splicing events.

Aberrant spliceosome activity via elevated intron retention and upregulation and phosphorylation of SF3B1 in chronic lymphocytic leukemia

Splicing factor 3b subunit 1 (SF3B1) is the largest subunit and core component of the spliceosome. Inhibition of SF3B1 was associated with an increase in broad intron retention (IR) on most transcripts, suggesting that IR can be used as a marker of spliceosome inhibition in chronic lymphocytic leukemia (CLL) cells. Furthermore, we separately analyzed exonic and intronic mapped reads on annotated RNA-sequencing transcripts obtained from B cells (n = 98 CLL patients) and healthy volunteers (n = 9). We measured intron/exon ratio to use that as a surrogate for alternative RNA splicing (ARS) and found that 66% of CLL-B cell transcripts had significant IR elevation compared with normal B cells (NBCs) and that correlated with mRNA downregulation and low expression levels. Transcripts with the highest IR levels belonged to biological pathways associated with gene expression and RNA splicing. A >2-fold increase of active pSF3B1 was observed in CLL-B cells compared with NBCs. Additionally, when the CLL-B cells were treated with macrolides (pladienolide-B), a significant decrease in pSF3B1, but not total SF3B1 protein, was observed. These findings suggest that IR/ARS is increased in CLL, which is associated with SF3B1 phosphorylation and susceptibility to SF3B1 inhibitors. These data provide additional support to the relevance of ARS in carcinogenesis and evidence of pSF3B1 participation in this process.

Differential transcript usage analysis incorporating quantification uncertainty via compositional measurement error regression modeling

Differential transcript usage (DTU) occurs when the relative expression of multiple transcripts arising from the same gene changes between different conditions. Existing approaches to detect DTU often rely on computational procedures that can have speed and scalability issues as the number of samples increases. Here we propose a new method, CompDTU, that uses compositional regression to model the relative abundance proportions of each transcript that are of interest in DTU analyses. This procedure leverages fast matrix-based computations that make it ideally suited for DTU analysis with larger sample sizes. This method also allows for the testing of and adjustment for multiple categorical or continuous covariates. Additionally, many existing approaches for DTU ignore quantification uncertainty in the expression estimates for each transcript in RNA-seq data. We extend our CompDTU method to incorporate quantification uncertainty leveraging common output from RNA-seq expression quantification tool in a novel method CompDTUme. Through several power analyses, we show that CompDTU has excellent sensitivity and reduces false positive results relative to existing methods. Additionally, CompDTUme results in further improvements in performance over CompDTU with sufficient sample size for genes with high levels of quantification uncertainty, while also maintaining favorable speed and scalability. We motivate our methods using data from the Cancer Genome Atlas Breast Invasive Carcinoma data set, specifically using RNA-seq data from primary tumors for 740 patients with breast cancer. We show greatly reduced computation time from our new methods as well as the ability to detect several novel genes with significant DTU across different breast cancer subtypes.

SF3A2 promotes progression and cisplatin resistance in triple-negative breast cancer via alternative splicing of MKRN1

Triple-negative breast cancer (TNBC) is the deadliest subtype of breast cancer owing to the lack of effective therapeutic targets. Splicing factor 3a subunit 2 (SF3A2), a poorly defined splicing factor, was notably elevated in TNBC tissues and promoted TNBC progression, as confirmed by cell proliferation, colony formation, transwell migration, and invasion assays. Mechanistic investigations revealed that E3 ubiquitin-protein ligase UBR5 promoted the ubiquitination-dependent degradation of SF3A2, which in turn regulated UBR5, thus forming a feedback loop to balance these two oncoproteins. Moreover, SF3A2 accelerated TNBC progression by, at least in part, specifically regulating the alternative splicing of makorin ring finger protein 1 (MKRN1) and promoting the expression of the dominant and oncogenic isoform, MKRN1-T1. Furthermore, SF3A2 participated in the regulation of both extrinsic and intrinsic apoptosis, leading to cisplatin resistance in TNBC cells. Collectively, these findings reveal a previously unknown role of SF3A2 in TNBC progression and cisplatin resistance, highlighting SF3A2 as a potential therapeutic target for patients with TNBC.

Interrogations of single-cell RNA splicing landscapes with SCASL define new cell identities with physiological relevance

RNA splicing shapes the gene regulatory programs that underlie various physiological and disease processes. Here, we present the SCASL (single-cell clustering based on alternative splicing landscapes) method for interrogating the heterogeneity of RNA splicing with single-cell RNA-seq data. SCASL resolves the issue of biased and sparse data coverage on single-cell RNA splicing and provides a new scheme for classifications of cell identities. With previously published datasets as examples, SCASL identifies new cell clusters indicating potentially precancerous and early-tumor stages in triple-negative breast cancer, illustrates cell lineages of embryonic liver development, and provides fine clusters of highly heterogeneous tumor-associated CD4 and CD8 T cells with functional and physiological relevance. Most of these findings are not readily available via conventional cell clustering based on single-cell gene expression data. Our study shows the potential of SCASL in revealing the intrinsic RNA splicing heterogeneity and generating biological insights into the dynamic and functional cell landscapes in complex tissues.

MYC-driven U2SURP regulates alternative splicing of SAT1 to promote triple-negative breast cancer progression

Triple-negative breast cancer (TNBC), although highly lethal, lacks validated therapeutic targets. Here, we report that U2 snRNP-associated SURP motif-containing protein (U2SURP), a poorly defined member of the serine/arginine rich protein family, was significantly upregulated in TNBC tissues, and its high expression was associated with poor prognosis of TNBC patients. MYC, a frequently amplified oncogene in TNBC tissues, enhanced U2SURP translation through an eIF3D (eukaryotic translation initiation factor 3 subunit D)-dependent mechanism, resulting in the accumulation of U2SURP in TNBC tissues. Functional assays revealed that U2SURP played an important role in facilitating tumorigenesis and metastasis of TNBC cells both in vitro and in vivo. Intriguingly, U2SURP had no significant effects on proliferative, migratory, and invasive potential of normal mammary epithelial cells. Furthermore, we found that U2SURP promoted alternative splicing of spermidine/spermine N1-acetyltransferase 1 (SAT1) pre-mRNA by removal of intron 3, resulting in an increase in the stability of SAT1 mRNA and subsequent protein expression levels. Importantly, spliced SAT1 promoted the oncogenic properties of TNBC cells, and re-expression of SAT1 in U2SURP-depleted cells partially rescued the impaired malignant phenotypes of TNBC cells caused by U2SURP knockdown both in vitro and in mice. Collectively, these findings reveal previously unknown functional and mechanism roles of the MYC-U2SURP-SAT1 signaling axis in TNBC progression and highlight U2SURP as a potential therapy target for TNBC.

Neoantigens: promising targets for cancer therapy

Recent advances in neoantigen research have accelerated the development and regulatory approval of tumor immunotherapies, including cancer vaccines, adoptive cell therapy and antibody-based therapies, especially for solid tumors. Neoantigens are newly formed antigens generated by tumor cells as a result of various tumor-specific alterations, such as genomic mutation, dysregulated RNA splicing, disordered post-translational modification, and integrated viral open reading frames. Neoantigens are recognized as non-self and trigger an immune response that is not subject to central and peripheral tolerance. The quick identification and prediction of tumor-specific neoantigens have been made possible by the advanced development of next-generation sequencing and bioinformatic technologies. Compared to tumor-associated antigens, the highly immunogenic and tumor-specific neoantigens provide emerging targets for personalized cancer immunotherapies, and serve as prospective predictors for tumor survival prognosis and immune checkpoint blockade responses. The development of cancer therapies will be aided by understanding the mechanism underlying neoantigen-induced anti-tumor immune response and by streamlining the process of neoantigen-based immunotherapies. This review provides an overview on the identification and characterization of neoantigens and outlines the clinical applications of prospective immunotherapeutic strategies based on neoantigens. We also explore their current status, inherent challenges, and clinical translation potential.

Sources of Cancer Neoantigens beyond Single-Nucleotide Variants

The success of checkpoint blockade therapy against cancer has unequivocally shown that cancer cells can be effectively recognized by the immune system and eliminated. However, the identity of the cancer antigens that elicit protective immunity remains to be fully explored. Over the last decade, most of the focus has been on somatic mutations derived from non-synonymous single-nucleotide variants (SNVs) and small insertion/deletion mutations (indels) that accumulate during cancer progression. Mutated peptides can be presented on MHC molecules and give rise to novel antigens or neoantigens, which have been shown to induce potent anti-tumor immune responses. A limitation with SNV-neoantigens is that they are patient-specific and their accurate prediction is critical for the development of effective immunotherapies. In addition, cancer types with low mutation burden may not display sufficient high-quality [SNV/small indels] neoantigens to alone stimulate effective T cell responses. Accumulating evidence suggests the existence of alternative sources of cancer neoantigens, such as gene fusions, alternative splicing variants, post-translational modifications, and transposable elements, which may be attractive novel targets for immunotherapy. In this review, we describe the recent technological advances in the identification of these novel sources of neoantigens, the experimental evidence for their presentation on MHC molecules and their immunogenicity, as well as the current clinical development stage of immunotherapy targeting these neoantigens.

Functional Interaction Between the Oncogenic Kinase NEK2 and Sam68 Promotes a Splicing Program Involved in Migration and Invasion in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) represents the most aggressive breast cancer subtype. Poor prognosis in TNBC is partly due to lack of efficacious targeted therapy and high propensity to metastasize. Dysregulation of alternative splicing has recently emerged as a trait of TNBC, suggesting that unveiling the molecular mechanisms underlying its regulation could uncover new druggable cancer vulnerabilities. The oncogenic kinase NEK2 is significantly upregulated in TNBC and contributes to shaping their unique splicing profile. Herein, we found that NEK2 interacts with the RNA binding protein Sam68 in TNBC cells and that NEK2-mediated phosphorylation of Sam68 enhances its splicing activity. Genome-wide transcriptome analyses identified the splicing targets of Sam68 in TNBC cells and revealed a common set of exons that are co-regulated by NEK2. Functional annotation of splicing-regulated genes highlighted cell migration and spreading as biological processes regulated by Sam68. Accordingly, Sam68 depletion reduces TNBC cell migration and invasion, and these effects are potentiated by the concomitant inhibition of NEK2 activity. Our findings indicate that Sam68 and NEK2 functionally cooperate in the regulation of a splicing program that sustains the pro-metastatic features of TNBC cells.

The oncogenic kinase NEK2 regulates an RBFOX2-dependent pro-mesenchymal splicing program in triple-negative breast cancer cells

BACKGROUND

Triple-negative breast cancer (TNBC) is the most heterogeneous and malignant subtype of breast cancer (BC). TNBC is defined by the absence of expression of estrogen, progesterone and HER2 receptors and lacks efficacious targeted therapies. NEK2 is an oncogenic kinase that is significantly upregulated in TNBC, thereby representing a promising therapeutic target. NEK2 localizes in the nucleus and promotes oncogenic splice variants in different cancer cells. Notably, alternative splicing (AS) dysregulation has recently emerged as a featuring trait of TNBC that contributes to its aggressive phenotype.

METHODS

To investigate whether NEK2 modulates TNBC transcriptome we performed RNA-sequencing analyses in a representative TNBC cell line (MDA-MB-231) and results were validated in multiple TNBC cell lines. Bioinformatics and functional analyses were carried out to elucidate the mechanism of splicing regulation by NEK2. Data from The Cancer Genome Atlas were mined to evaluate the potential of NEK2-sensitive exons as markers to identify the TNBC subtype and to assess their prognostic value.

RESULTS

Transcriptome analysis revealed a widespread impact of NEK2 on the transcriptome of TNBC cells, with 1830 AS events that are susceptible to its expression. NEK2 regulates the inclusion of cassette exons in splice variants that discriminate TNBC from other BC and that correlate with poor prognosis, suggesting that this kinase contributes to the TNBC-specific splicing program. NEK2 elicits its effects by modulating the expression of the splicing factor RBFOX2, a well-known regulator of epithelial to mesenchymal transition (EMT). Accordingly, NEK2 splicing-regulated genes are enriched in functional terms related to cell adhesion and contractile cytoskeleton and NEK2 depletion in mesenchymal TNBC cells induces phenotypic and molecular traits typical of epithelial cells. Remarkably, depletion of select NEK2-sensitive splice-variants that are prognostic in TNBC patients is sufficient to interfere with TNBC cell morphology and motility, suggesting that NEK2 orchestrates a pro-mesenchymal splicing program that modulates migratory and invasive properties of TNBC cells.

CONCLUSIONS

Our study uncovers an extensive splicing program modulated by NEK2 involving splice variants that confer an invasive phenotype to TNBCs and that might represent, together with NEK2 itself, valuable therapeutic targets for this disease.

Transcript-targeted analysis reveals isoform alterations and double-hop fusions in breast cancer

Although transcriptome alteration is an essential driver of carcinogenesis, the effects of chromosomal structural alterations on the cancer transcriptome are not yet fully understood. Short-read transcript sequencing has prevented researchers from directly exploring full-length transcripts, forcing them to focus on individual splice sites. Here, we develop a pipeline for Multi-Sample long-read Transcriptome Assembly (MuSTA), which enables construction of a transcriptome from long-read sequence data. Using the constructed transcriptome as a reference, we analyze RNA extracted from 22 clinical breast cancer specimens. We identify a comprehensive set of subtype-specific and differentially used isoforms, which extended our knowledge of isoform regulation to unannotated isoforms including a short form TNS3. We also find that the exon-intron structure of fusion transcripts depends on their genomic context, and we identify double-hop fusion transcripts that are transcribed from complex structural rearrangements. For example, a double-hop fusion results in aberrant expression of an endogenous retroviral gene, ERVFRD-1, which is normally expressed exclusively in placenta and is thought to protect fetus from maternal rejection; expression is elevated in several TCGA samples with ERVFRD-1 fusions. Our analyses provide direct evidence that full-length transcript sequencing of clinical samples can add to our understanding of cancer biology and genomics in general.

SRSF1-dependent inhibition of C9ORF72-repeat RNA nuclear export: genome-wide mechanisms for neuroprotection in amyotrophic lateral sclerosis

Background

Loss of motor neurons in amyotrophic lateral sclerosis (ALS) leads to progressive paralysis and death. Dysregulation of thousands of RNA molecules with roles in multiple cellular pathways hinders the identification of ALS-causing alterations over downstream changes secondary to the neurodegenerative process. How many and which of these pathological gene expression changes require therapeutic normalisation remains a fundamental question.

Methods

Here, we investigated genome-wide RNA changes in C9ORF72-ALS patient-derived neurons and Drosophila, as well as upon neuroprotection taking advantage of our gene therapy approach which specifically inhibits the SRSF1-dependent nuclear export of pathological C9ORF72-repeat transcripts. This is a critical study to evaluate (i) the overall safety and efficacy of the partial depletion of SRSF1, a member of a protein family involved itself in gene expression, and (ii) a unique opportunity to identify neuroprotective RNA changes.

Results

Our study shows that manipulation of 362 transcripts out of 2257 pathological changes, in addition to inhibiting the nuclear export of repeat transcripts, is sufficient to confer neuroprotection in C9ORF72-ALS patient-derived neurons. In particular, expression of 90 disease-altered transcripts is fully reverted upon neuroprotection leading to the characterisation of a human C9ORF72-ALS disease-modifying gene expression signature. These findings were further investigated in vivo in diseased and neuroprotected Drosophila transcriptomes, highlighting a list of 21 neuroprotective changes conserved with 16 human orthologues in patient-derived neurons. We also functionally validated the high neuroprotective potential of one of these disease-modifying transcripts, demonstrating that inhibition of ALS-upregulated human KCNN1–3 (Drosophila SK) voltage-gated potassium channel orthologs mitigates degeneration of human motor neurons and Drosophila motor deficits.

Conclusions

Strikingly, the partial depletion of SRSF1 leads to expression changes in only a small proportion of disease-altered transcripts, indicating that not all RNA alterations need normalization and that the gene therapeutic approach is safe in the above preclinical models as it does not disrupt globally gene expression. The efficacy of this intervention is also validated at genome-wide level with transcripts modulated in the vast majority of biological processes affected in C9ORF72-ALS. Finally, the identification of a characteristic signature with key RNA changes modified in both the disease state and upon neuroprotection also provides potential new therapeutic targets and biomarkers.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13024-021-00475-y.

Chromatin loop anchors predict transcript and exon usage

Epigenomics and transcriptomics data from high-throughput sequencing techniques such as RNA-seq and ChIP-seq have been successfully applied in predicting gene transcript expression. However, the locations of chromatin loops in the genome identified by techniques such as Chromatin Interaction Analysis with Paired End Tag sequencing (ChIA-PET) have never been used for prediction tasks. Here, we developed machine learning models to investigate if ChIA-PET could contribute to transcript and exon usage prediction. In doing so, we used a large set of transcription factors as well as ChIA-PET data. We developed different Gradient Boosting Trees models according to the different tasks with the integrated datasets from three cell lines, including GM12878, HeLaS3 and K562. We validated the models via 10-fold cross validation, chromosome-split validation and cross-cell validation. Our results show that both transcript and splicing-derived exon usage can be effectively predicted with at least 0.7512 and 0.7459 of accuracy, respectively, on all cell lines from all kinds of validations. Examining the predictive features, we found that RNA Polymerase II ChIA-PET was one of the most important features in both transcript and exon usage prediction, suggesting that chromatin loop anchors are predictive of both transcript and exon usage.

Development of alternative splicing signature in lung squamous cell carcinoma

Increasing evidence demonstrated that alternative splicing (AS) plays a vital role in tumorigenesis and clinical outcome of patient. However, systematical analysis of AS in lung squamous cell carcinoma (LUSC) is lacking and greatly necessary. Thus, this study was to systematically estimate the function of AS events served as prognostic indicators in LUSC. Among 31,345 mRNA AS events in 9633 genes, we detected 1996 AS in 1409 genes which have significant connection with overall survival (OS) of LUSC patients. Then, prognostic model based on seven types of AS events was established and we further constructed a combined prognostic model. The Kaplan–Meier curve results suggested that seven types of AS signatures and the combined prognostic model could exhibit robust performance in predicting prognosis. Patients in the high-risk group had significantly shorter OS than those in the low-risk group. The ROC showed all prognostic models had high accuracy and powerful predictive performance with different AUC ranging from 0.837 to 0.978. Moreover, the combined prognostic model had highest performance in risk stratification and predictive accuracy than single prognostic models and had higher accuracy than other mRNA model. Finally, a significant correlation network between survival-related AS genes and prognostic splicing factors (SFs) was established. In conclusion, our study provided several potential prognostic AS models and constructed splicing network between AS and SFs in LUSC, which could be used as potential indicators and treatment targets for LUSC patients.

Differential alternative RNA splicing and transcription events between tumors from African American and White patients in The Cancer Genome Atlas

Individuals of African ancestry suffer disproportionally from higher incidence, aggressiveness, and mortality for particular cancers. This disparity likely results from an interplay among differences in multiple determinants of health, including differences in tumor biology. We used The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA aggregate expression datasets and identified differential alternative RNA splicing and transcription events (ARS/T) in cancers between self-identified African American (AA) and White (W) patients. We found that retained intron events were enriched among race-related ARS/T. In addition, on average, 12% of the most highly ranked race-related ARS/T overlapped between any two analyzed cancers. Moreover, the genes undergoing race-related ARS/T functioned in cancer-promoting pathways, and a number of race-related ARS/T were associated with patient survival. We built a web-application, CanSplice, to mine genomic datasets by self-identified race. The race-related targets have the potential to aid in the development of new biomarkers and therapeutics to mitigate cancer disparity.

Functional and conformational impact of cancer-associated SF3B1 mutations depends on the position and the charge of amino acid substitution

The hotspot mutations of SF3B1, the most frequently mutated splicing gene in cancers, contribute to oncogenesis by corrupting the mRNA splicing. Further SF3B1 mutations have been reported in cancers but their consequences remain unclear. Here, we screened for SF3B1 mutations in the vicinity of the hotspot region in tumors. We then performed in-silico prediction of the functional outcome followed by in-cellulo modelling of different SF3B1 mutants. We show that cancer-associated SF3B1 mutations present varying functional consequences that are loosely predicted by the in-silico algorithms. Analysis of the tertiary structure of SF3B1 mutants revealed that the resulting splicing errors may be due to a conformational change in SF3B1 N-terminal region, which mediates binding with other splicing factors. Our study demonstrates a varying functional impact of SF3B1 mutations according to the mutated codon and the amino acid substitution, implying unequal pathogenic and prognostic potentials of SF3B1 mutations in cancers.

iMOKA: k-mer based software to analyze large collections of sequencing data

iMOKA (interactive multi-objective k-mer analysis) is a software that enables comprehensive analysis of sequencing data from large cohorts to generate robust classification models or explore specific genetic elements associated with disease etiology. iMOKA uses a fast and accurate feature reduction step that combines a Naïve Bayes classifier augmented by an adaptive entropy filter and a graph-based filter to rapidly reduce the search space. By using a flexible file format and distributed indexing, iMOKA can easily integrate data from multiple experiments and also reduces disk space requirements and identifies changes in transcript levels and single nucleotide variants. iMOKA is available at https://github.com/RitchieLabIGH/iMOKA and Zenodo https://doi.org/10.5281/zenodo.4008947 .

The Possible Influence of Mediterranean Diet on Extracellular Vesicle miRNA Expression in Breast Cancer Survivors

The Mediterranean diet (MD) has been reported to have beneficial effects on breast cancer and cardiovascular diseases. Recently, microRNAs (miRNAs) have been suggested as biomarkers for the diagnosis and disease prognosis in cancer and cardiovascular diseases. We evaluated the influence of the MD on the plasma-derived extracellular vesicle miRNA signature of overweight breast cancer survivors. Sixteen participants instructed to adhere to the MD for eight weeks were included in this study. To curate differentially expressed miRNAs after MD intervention, we employed two methods: significance analysis of microarrays and DESeq2. The selected miRNAs were analyzed using ingenuity pathway analysis. After an eight-week intervention, body mass index, waist circumference, fasting glucose, fasting insulin, and homeostatic model assessment for insulin resistance were significantly improved. Expression levels of 798 miRNAs were comprehensively analyzed, and 42 extracellular vesicle miRNAs were significantly differentially regulated after the eight-week MD (36 were up and 6 were down-regulated). We also identified enriched pathways in genes regulated by differentially expressed 42 miRNAs, which include signaling associated with breast cancer, energy metabolism, glucose metabolism, and insulin. Our study indicates that extracellular vesicle miRNAs differentially expressed as a result of the MD might be involved in the mechanisms that relate to cardiometabolic risk factors in overweight breast cancer survivors.

RNA Splicing and Cancer

RNA splicing is an essential process that governs many aspects of cellular proliferation, survival, and differentiation. Considering the importance of RNA splicing in gene regulation, alterations in this pathway have been implicated in many human cancers. Large-scale genomic studies have uncovered a spectrum of splicing machinery mutations that contribute to tumorigenesis. Moreover, cancer cells are capable of hijacking the expression of RNA-binding proteins (RBPs), leading to dysfunctional gene splicing and tumor-specific dependencies. Advances in next-generation RNA sequencing have revealed tumor-specific isoforms associated with these alterations, including the presence of neoantigens, which serve as potential immunotherapeutic targets. In this review, we discuss the various mechanisms by which cancer cells exploit RNA splicing to promote tumor growth and the current therapeutic landscape for splicing-based therapies.

Epigenetic regulation of the ITGB4 gene in prostate cancer

BACKGROUND

Examination of epigenetic changes at the ITGB4 gene promoter reveals altered methylation at different stages of prostate tumour progression and these changes may, in part, explain the complex patterns of gene expression of this integrin observed. Transcriptional re-programming perturbs expression of cell adhesion molecules and underpins metastatic tumour cell behaviour. Decreasing expression of the cell adhesion molecule ITGB4, which encodes the beta subunit of the integrin, alpha6 beta4 (α6β4), has been correlated with increased tumour aggressiveness and metastasis in multiple tumour types including prostate cancer. Paradoxically, in vitro studies in tumour cell models demonstrate that ITGB4 mediates cell mobility and invasion. Herein we examined whether transcriptional re-programming by methylation influenced ITGB4 gene expression at different stages of prostate cancer progression. Bisulphite sequencing of a large CpG island in the ITGB4 gene promoter identified differentially methylated regions in prostate cancer cell lines representing a localised tumour (22Rv1), lymph node metastasis (LNCaP), and a bone metastasis (PC-3). The highest levels of methylation were observed in the CpG island surrounding the ITGB4 transcription start site in PC-3 cells, and this observation also correlated with higher gene expression of ITGB4 in these cells. Furthermore, PC-3 cells expressed two distinct transcripts, using an alternate transcription start site, which was not detected in other cell lines. In prostate tumour biopsy samples, patterns of methylation across the ITGB4 promoter were similar overall in matched primary and metastatic samples (n = 4 pairs), with a trend toward loss of methylation at specific sites in metastatic lesions.

Roles and mechanisms of alternative splicing in cancer - implications for care

Removal of introns from messenger RNA precursors (pre-mRNA splicing) is an essential step for the expression of most eukaryotic genes. Alternative splicing enables the regulated generation of multiple mRNA and protein products from a single gene. Cancer cells have general as well as cancer type-specific and subtype-specific alterations in the splicing process that can have prognostic value and contribute to every hallmark of cancer progression, including cancer immune responses. These splicing alterations are often linked to the occurrence of cancer driver mutations in genes encoding either core components or regulators of the splicing machinery. Of therapeutic relevance, the transcriptomic landscape of cancer cells makes them particularly vulnerable to pharmacological inhibition of splicing. Small-molecule splicing modulators are currently in clinical trials and, in addition to splice site-switching antisense oligonucleotides, offer the promise of novel and personalized approaches to cancer treatment.

Putatively cancer-specific exon-exon junctions are shared across patients and present in developmental and other non-cancer cells

This study probes the distribution of putatively cancer-specific junctions across a broad set of publicly available non-cancer human RNA sequencing (RNA-seq) datasets. We compared cancer and non-cancer RNA-seq data from The Cancer Genome Atlas (TCGA), the Genotype-Tissue Expression (GTEx) Project and the Sequence Read Archive. We found that (i) averaging across cancer types, 80.6% of exon-exon junctions thought to be cancer-specific based on comparison with tissue-matched samples (σ = 13.0%) are in fact present in other adult non-cancer tissues throughout the body; (ii) 30.8% of junctions not present in any GTEx or TCGA normal tissues are shared by multiple samples within at least one cancer type cohort, and 87.4% of these distinguish between different cancer types; and (iii) many of these junctions not found in GTEx or TCGA normal tissues (15.4% on average, σ = 2.4%) are also found in embryological and other developmentally associated cells. These findings refine the meaning of RNA splicing event novelty, particularly with respect to the human neoepitope repertoire. Ultimately, cancer-specific exon-exon junctions may have a substantial causal relationship with the biology of disease.

Systemic Analysis of RNA Alternative Splicing Signals Related to the Prognosis for Head and Neck Squamous Cell Carcinoma

Alternative splicing (AS) is an important mechanism that is responsible for the production of protein diversity. An increasing body of evidence has suggested that out-of-control AS is closely related to the genesis and development of cancer. Systematic analysis of genome-wide AS in head and neck squamous cell carcinoma (HNSCC) has not yet been carried out, and consideration of this topic remains at the preliminary stage and requires further investigation. In this study, systemic bioinformatic analysis was carried out on the genome-wide AS events of 555 clinical HNSCC samples from the TCGA database. Firstly, we statistically analyzed the distributions of seven AS event types in HNSCC samples. Then, through univariate survival analysis, we observed the relationship between AS and the prognosis of the disease and found that 437 intersections of AS events were significantly related to overall survival. Among them, 335 cross-genes showed a high degree of consistency in the genes associated with overall survival and recurrence. The overall survival was significantly related to AS events. Besides, the frequency of overall survival-related ES events was evidently reduced, while the AP and the AT events were increased. In addition, AT events accounted for the largest proportion. Further, multiple regression model analysis proved that AS could become a new classification method for HNSCC, and KEGG enrichment analysis proved that most genes and proteins interacting with AS events had different biological functions and were associated with a variety of diseases. Finally, through the selection of characteristic HNSCC genes and the construction of a prognostic model, seven cross-genes related to survival and recurrence were screened out, and these characteristic genes were verified by multivariate survival model analysis so as to classify the prognosis at different splicing times and gene expression levels. These results have laid a solid foundation for our further research and play a decisive role in showing the correlation of AS with the prognosis of HNSCC.

Genome-wide analysis reveals the association between alternative splicing and DNA methylation across human solid tumors

Background

Dysregulation of alternative splicing (AS) is a critical signature of cancer. However, the regulatory mechanisms of cancer-specific AS events, especially the impact of DNA methylation, are poorly understood.

Methods

By using The Cancer Genome Atlas (TCGA) SpliceSeq and TCGA data for ten solid tumor types, association analysis was performed to characterize the potential link between cancer-specific AS and DNA methylation. Functional and pathway enrichment analyses were performed, and the protein-protein interaction (PPI) network was constructed with the String website. The prognostic analysis was carried out with multivariate Cox regressions models.

Results

15,818 AS events in 3955 annotated genes were identified across ten solid tumor types. The different DNA methylation patterns between tumor and normal tissues at the corresponding alternative spliced exon boundaries were shown, and 51.3% of CpG sites (CpGs) revealed hypomethylated in tumors. Notably, 607 CpGs were found to be highly correlated with 369 cancer-specific AS events after permutation tests. Among them, the hypomethylated CpGs account for 52.7%, and the number of down-regulated exons was 173. Furthermore, we found 38 AS events in 35 genes could serve as new molecular biomarkers to predict patient survival.

Conclusions

Our study described the relationship between DNA methylation and AS events across ten human solid tumor types and provided new insights into intragenic DNA methylation and exon usage during the AS process.

Targeting mRNA processing as an anticancer strategy

Discoveries in the past decade have highlighted the potential of mRNA as a therapeutic target for cancer. Specifically, RNA sequencing revealed that, in addition to gene mutations, alterations in mRNA can contribute to the initiation and progression of cancer. Indeed, precursor mRNA processing, which includes the removal of introns by splicing and the formation of 3' ends by cleavage and polyadenylation, is frequently altered in tumours. These alterations result in numerous cancer-specific mRNAs that generate altered levels of normal proteins or proteins with new functions, leading to the activation of oncogenes or the inactivation of tumour-suppressor genes. Abnormally spliced and polyadenylated mRNAs are also associated with resistance to cancer treatment and, unexpectedly, certain cancers are highly sensitive to the pharmacological inhibition of splicing. This Review summarizes recent progress in our understanding of how splicing and polyadenylation are altered in cancer and highlights how this knowledge has been translated for drug discovery, resulting in the production of small molecules and oligonucleotides that modulate the spliceosome and are in clinical trials for the treatment of cancer.

Aberrant RNA Splicing in Cancer and Drug Resistance

More than 95% of the 20,000 to 25,000 transcribed human genes undergo alternative RNA splicing, which increases the diversity of the proteome. Isoforms derived from the same gene can have distinct and, in some cases, opposing functions. Accumulating evidence suggests that aberrant RNA splicing is a common and driving event in cancer development and progression. Moreover, aberrant splicing events conferring drug/therapy resistance in cancer is far more common than previously envisioned. In this review, aberrant splicing events in cancer-associated genes, namely BCL2L1, FAS, HRAS, CD44, Cyclin D1, CASP2, TMPRSS2-ERG, FGFR2, VEGF, AR and KLF6, will be discussed. Also highlighted are the functional consequences of aberrant splice variants (BCR-Abl35INS, BIM-γ, IK6, p61 BRAF V600E, CD19-∆2, AR-V7 and PIK3CD-S) in promoting resistance to cancer targeted therapy or immunotherapy. To overcome drug resistance, we discuss opportunities for developing novel strategies to specifically target the aberrant splice variants or splicing machinery that generates the splice variants. Therapeutic approaches include the development of splice variant-specific siRNAs, splice switching antisense oligonucleotides, and small molecule inhibitors targeting splicing factors, splicing factor kinases or the aberrant oncogenic protein isoforms.

Comprehensive Analysis of Alternative Splicing Across Tumors from 8,705 Patients

Our comprehensive analysis of alternative splicing across 32 The Cancer Genome Atlas cancer types from 8,705 patients detects alternative splicing events and tumor variants by reanalyzing RNA and whole-exome sequencing data. Tumors have up to 30% more alternative splicing events than normal samples. Association analysis of somatic variants with alternative splicing events confirmed known trans associations with variants in SF3B1 and U2AF1 and identified additional trans-acting variants (e.g., TADA1, PPP2R1A). Many tumors have thousands of alternative splicing events not detectable in normal samples; on average, we identified ≈930 exon-exon junctions ("neojunctions") in tumors not typically found in GTEx normals. From Clinical Proteomic Tumor Analysis Consortium data available for breast and ovarian tumor samples, we confirmed ≈1.7 neojunction- and ≈0.6 single nucleotide variant-derived peptides per tumor sample that are also predicted major histocompatibility complex-I binders ("putative neoantigens").

Profiling plasma extracellular vesicle by pluronic block-copolymer based enrichment method unveils features associated with breast cancer aggression, metastasis and invasion

Extracellular vesicle (EV)-based liquid biopsies have been proposed to be a readily obtainable biological substrate recently for both profiling and diagnostics purposes. Development of a fast and reliable preparation protocol to enrich such small particles could accelerate the discovery of informative, disease-related biomarkers. Though multiple EV enrichment protocols are available, in terms of efficiency, reproducibility and simplicity, precipitation-based methods are most amenable to studies with large numbers of subjects. However, the selectivity of the precipitation becomes critical. Here, we present a simple plasma EV enrichment protocol based on pluronic block copolymer. The enriched plasma EV was able to be verified by multiple platforms. Our results showed that the particles enriched from plasma by the copolymer were EV size vesicles with membrane structure; proteomic profiling showed that EV-related proteins were significantly enriched, while high-abundant plasma proteins were significantly reduced in comparison to other precipitation-based enrichment methods. Next-generation sequencing confirmed the existence of various RNA species that have been observed in EVs from previous studies. Small RNA sequencing showed enriched species compared to the corresponding plasma. Moreover, plasma EVs enriched from 20 advanced breast cancer patients and 20 age-matched non-cancer controls were profiled by semi-quantitative mass spectrometry. Protein features were further screened by EV proteomic profiles generated from four breast cancer cell lines, and then selected in cross-validation models. A total of 60 protein features that highly contributed in model prediction were identified. Interestingly, a large portion of these features were associated with breast cancer aggression, metastasis as well as invasion, consistent with the advanced clinical stage of the patients. In summary, we have developed a plasma EV enrichment method with improved precipitation selectivity and it might be suitable for larger-scale discovery studies.