Robust stratification of breast cancer subtypes using differential patterns of transcript isoform expression

Comprehensive Long-Read Sequencing Analysis Discloses the Transcriptome Features of Papillary Thyroid Microcarcinoma

CONTEXT

Papillary thyroid microcarcinoma (PTMC) is the most common type of thyroid cancer. It has been shown that lymph node metastasis is associated with poor prognosis in patients with PTMC.

OBJECTIVE

We aim to characterize the PTMC transcriptome landscape and identify the candidate transcripts that are associated with lateral neck lymph node metastasis of PTMC.

METHODS

We performed full-length transcriptome sequencing in 64 PTMC samples. Standard bioinformatic pipelines were applied to characterize and annotate the full-length expression profiles of 2 PTMC subtypes. Functional open reading frame (ORF) annotation of the known and novel transcripts were predicted by HMMER, DeepLoc, and DeepTMHMM tools. Candidate transcripts associated with the pN1b subtype were identified after transcript quantification and differential gene expression analyses.

RESULTS

We found that skipping exons accounted for the more than 27.82% of the alternative splicing events. At least 42.56% of the discovered transcripts were novel isoforms of annotated genes. A total of 39 193 ORFs in novel transcripts and 18 596 ORFs in known transcripts were identified. Distribution patterns of the characterized transcripts in functional domain, subcellular localization, and transmembrane structure were predicted. In total, 1033 and 1204 differentially expressed genes were identified in the pN0 and pN1b groups, respectively. Moreover, novel isoforms of FRMD3, NOD1, and SHROOM4 were highlighted for their association with pN1b subtype.

CONCLUSION

Our data provided the global transcriptome landscape of PTMC and also revealed the novel isoforms that associated with PTMC aggressiveness.

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.

Supervised Relevance-Redundancy assessments for feature selection in omics-based classification scenarios

BACKGROUND AND OBJECTIVE

Many classification tasks in translational bioinformatics and genomics are characterized by the high dimensionality of potential features and unbalanced sample distribution among classes. This can affect classifier robustness and increase the risk of overfitting, curse of dimensionality and generalization leaks; furthermore and most importantly, this can prevent obtaining adequate patient stratification required for precision medicine in facing complex diseases, like cancer. Setting up a feature selection strategy able to extract only proper predictive features by removing irrelevant, redundant, and noisy ones is crucial to achieving valuable results on the desired task.

METHODS

We propose a new feature selection approach, called ReRa, based on supervised Relevance-Redundancy assessments. ReRa consists of a customized step of relevance-based filtering, to identify a reduced subset of meaningful features, followed by a supervised similarity-based procedure to minimize redundancy. This latter step innovatively uses a combination of global and class-specific similarity assessments to remove redundant features while preserving those differentiated across classes, even when these classes are strongly unbalanced.

RESULTS

We compared ReRa with several existing feature selection methods to obtain feature spaces on which performing breast cancer patient subtyping using several classifiers: we considered two use cases based on gene or transcript isoform expression. In the vast majority of the assessed scenarios, when using ReRa-selected feature spaces, the performances were significantly increased compared to simple feature filtering, LASSO regularization, or even MRmr - another Relevance-Redundancy method. The two use cases represent an insightful example of translational application, taking advantage of ReRa capabilities to investigate and enhance a clinically-relevant patient stratification task, which could be easily applied also to other cancer types and diseases.

CONCLUSIONS

ReRa approach has the potential to improve the performance of machine learning models used in an unbalanced classification scenario. Compared to another Relevance-Redundancy approach like MRmr, ReRa does not require tuning the number of preserved features, ensures efficiency and scalability over huge initial dimensionalities and allows re-evaluation of all previously selected features at each iteration of the redundancy assessment, to ultimately preserve only the most relevant and class-differentiated features.

MAGOH is correlated with poor prognosis and is essential for cell proliferation in lower-grade glioma

OBJECTIVE

Mago-nashi homolog (MAGOH) has been shown to play a pivotal part in various tumors. However, its specific contribution in lower-grade glioma (LGG) is still unknown.

METHODS

Pan-cancer analysis was implemented to inspect the expression characteristics and prognostic significance of MAGOH in multiple tumors. The associations between MAGOH expression patterns and the pathological features of LGG were analyzed, as were the connections between MAGOH expression and the clinical traits, prognosis, biological activities, immune features, genomic variations, and responses to treatment in LGG. Additionally, in vitro studies were performed to detect the expression levels and biomedical functions of MAGOH in LGG.

RESULTS

Abnormally increased levels of MAGOH expression were connected with adverse prognosis in patients with several types of tumors, including LGG. Importantly, we found that levels of MAGOH expression were independent prognostic biomarker of patients with LGG. Increased MAGOH expression was also highly associated with several immune-related markers, immune cell infiltration, immune checkpoint genes (ICPGs), gene mutations, and responses to chemotherapy in patients with LGG. In vitro studies ascertained that abnormally increased MAGOH was essential for cell proliferation in LGG.

CONCLUSION

MAGOH is a valid predictive biomarker in LGG and may become a novel therapeutic target in these patients.

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.

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.

Multifaceted roles of MAGOH Proteins

MAGOH and MAGOHB are paralog proteins that can substitute each other in the exon junction complex (EJC). The EJC is formed of core components EIF4A3, RBM8A, and MAGOH/MAGOHB. As a part of the EJC, MAGOH proteins are required for mRNA splicing, export, translation and nonsense-mediated mRNA decay (NMD). MAGOH is also essential for embryonic development and normal cellular functioning. The haploinsufficiency of MAGOH results in disorders such as microcephaly and cancer. The present review discusses the discovery of MAGOH, its paralog MAGOHB, their roles in cellular function as part of the EJC, and other cellular roles that are not directly associated with mRNA processing. We also discuss how MAGOH haploinsufficiency in cancer cells can be exploited to develop a novel targeted cancer treatment.

Research progress on RNA−binding proteins in breast cancer

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

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.

A comprehensive long-read isoform analysis platform and sequencing resource for breast cancer

Tumors display widespread transcriptome alterations, but the full repertoire of isoform-level alternative splicing in cancer is unknown. We developed a long-read (LR) RNA sequencing and analytical platform that identifies and annotates full-length isoforms and infers tumor-specific splicing events. Application of this platform to breast cancer samples identifies thousands of previously unannotated isoforms; ~30% affect protein coding exons and are predicted to alter protein localization and function. We performed extensive cross-validation with -omics datasets to support transcription and translation of novel isoforms. We identified 3059 breast tumor–specific splicing events, including 35 that are significantly associated with patient survival. Of these, 21 are absent from GENCODE and 10 are enriched in specific breast cancer subtypes. Together, our results demonstrate the complexity, cancer subtype specificity, and clinical relevance of previously unidentified isoforms and splicing events in breast cancer that are only annotatable by LR-seq and provide a rich resource of immuno-oncology therapeutic targets.

Integrative Modeling of Multiomics Data for Predicting Tumor Mutation Burden in Patients with Lung Cancer

Immunotherapy has been widely used in the treatment of lung cancer, and one of the most effective biomarkers for the prognosis of immunotherapy currently is tumor mutation burden (TMB). Although whole-exome sequencing (WES) could be utilized to assess TMB, several problems prevent its routine clinical application. To develop a simplified TMB prediction model, patients with lung adenocarcinoma (LUAD) in The Cancer Genome Atlas (TCGA) were randomly split into training and validation cohorts and categorized into the TMB-high (TMB-H) and TMB-low (TMB-L) groups, respectively. Based on the 610 differentially expressed genes, 50 differentially expressed miRNAs and 58 differentially methylated CpG sites between TMB-H and TMB-L patients, we constructed 4 predictive signatures and established TMB prediction model through machine learning methods that integrating the expression or methylation profiles of 7 genes, 7 miRNAs, and 6 CpG sites. The multiomics model exhibited excellent performance in predicting TMB with the area under curve (AUC) of 0.911 in the training cohort and 0.859 in the validation cohort. Besides, the significant correlation between the multiomics model score and TMB was observed. In summary, we developed a prognostic TMB prediction model by integrating multiomics data in patients with LUAD, which might facilitate the further development of quantitative real time-polymerase chain reaction- (qRT-PCR-) based TMB prediction assay.

Proteomic analysis of inhibitor of apoptosis protein‑like protein‑2 on breast cancer cell proliferation

Although inhibitor of apoptosis protein-like protein-2 (ILP-2) is considered to be a novel enhancer of breast cancer proliferation, its underlying mechanism of action remains unknown. Therefore, the present study aimed to investigate the expression profile of ILP-2-related proteins in MCF-7 cells to reveal their effect on promoting breast cancer cell proliferation. The isobaric tags for relative and absolute quantification (iTRAQ) method was used to analyse the expression profile of ILP-2-related proteins in MCF-7 breast cancer cells transfected with small interfering (si)RNA against ILP-2 (siRNA-5 group) and the negative control (NC) siRNA. The analysis of the iTRAQ data was carried out using western blotting and reverse transcription-quantitative PCR. A total of 4,065 proteins were identified in MCF-7 cells, including 241 differentially expressed proteins (DEPs; fold change ≥1.20 or ≤0.83; P<0.05). Among them, 156 proteins were upregulated and 85 were downregulated in the siRNA-5 group compared with in the NC group. The aforementioned DEPs were mainly enriched in ‘ECM-receptor interaction’. In addition, the top 10 biological processes related to these proteins were associated with signal transduction, cell proliferation and immune system processes. Furthermore, ILP-2 silencing upregulated N(4)-(β-N-acetylglucosaminyl)-L-asparaginase, metallothionein-1E and tryptophan 2,3-dioxygenase, whereas ILP-2 overexpression exerted the opposite effect. The results of the present study suggested that ILP-2 could promote breast cancer growth via regulating cell proliferation, signal transduction, immune system processes and other cellular physiological activities.

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.

The Estrogen Receptor α Signaling Pathway Controls Alternative Splicing in the Absence of Ligands in Breast Cancer Cells

Background: The transcriptional activity of estrogen receptor α (ERα) in breast cancer (BC) is extensively characterized. Our group has previously shown that ERα controls the expression of a number of genes in its unliganded form (apoERα), among which a large group of RNA-binding proteins (RBPs) encode genes, suggesting its role in the control of co- and post-transcriptional events. Methods: apoERα-mediated RNA processing events were characterized by the analysis of transcript usage and alternative splicing changes in an RNA-sequencing dataset from MCF-7 cells after siRNA-induced ERα downregulation. Results: ApoERα depletion induced an expression change of 681 RBPs, including 84 splicing factors involved in translation, ribonucleoprotein complex assembly, and 3′end processing. ApoERα depletion results in 758 isoform switching events with effects on 3′end length and the splicing of alternative cassette exons. The functional enrichment of these events shows that post-transcriptional regulation is part of the mechanisms by which apoERα controls epithelial-to-mesenchymal transition and BC cell proliferation. In primary BCs, the inclusion levels of the experimentally identified alternatively spliced exons are associated with overall and disease-free survival. Conclusion: Our data supports the role of apoERα in maintaining the luminal phenotype of BC cells by extensively regulating gene expression at the alternative splicing level.

Pseudogene RPL32P3 regulates the blood-tumor barrier permeability via the YBX2/HNF4G axis

The existence of the blood–tumor barrier (BTB) severely hinders the transport of anti-tumor drugs to brain tumor tissues. Selectively opening BTB is of great significance to improve the chemotherapy effect of glioma. Pseudogenes have been recognized as important regulators in various biologic processes. In this study, we identified that ribosomal protein L32 pseudogene 3 (RPL32P3) was highly expressed in glioma-exposed endothelial cells (GECs). Knockdown of RPL32P3 decreased the expression of tight junction-related proteins (TJPs) and increased BTB permeability. Subsequent analysis of the underlying mechanism indicated that RPL32P3 recruited lysine methyltransferase 2 A (KMT2A) to the Y-box binding protein 2 (YBX2) promoter region and mediated H3K4me3 to promote YBX2 transcription. Highly expressed YBX2 bound and stabilized hepatocyte nuclear factor 4 gamma (HNF4G) mRNA. Highly expressed HNF4G directly bound to the promoters of TJPs ZO-1, occludin and claudin-5 to promote their transcriptional activities and regulated BTB permeability. The simultaneous knockdown of RPL32P3, YBX2, and HNF4G combined with doxorubicin (DOX) increased the apoptosis of glioma cells. In conclusion, the current study indicated that RPL32P3 knockdown increased BTB permeability through the YBX2/HNF4G pathway. These findings may provide new targets for the comprehensive treatment of glioma.

Alternative splicing perturbation landscape identifies RNA binding proteins as potential therapeutic targets in cancer

Alternative splicing (AS) plays an important role in gene regulation, and AS perturbations are frequently observed in cancer. RNA binding protein (RBP) is one of the molecular determinants of AS, and perturbations in RBP-gene network activity are causally associated with cancer development. Here, we performed a systematic analysis to characterize the perturbations in AS events across 18 cancer types. We showed that AS alterations were prevalent in cancer and involved in cancer-related pathways. Given that the extent of AS perturbation was associated with disease severity, we proposed a computational pipeline to identify RBP regulators. Pan-cancer analysis identified a number of conserved RBP regulators, which play important roles in regulating AS of genes involved in cancer hallmark pathways. Our application analysis revealed that the expression of 68 RBP regulators helped in cancer subtyping. Specifically, we identified four subtypes of kidney cancer with differences in cancer hallmark pathway activities and prognosis. Finally, we identified the small molecules that can potentially target the RBP genes and suggested potential candidates for cancer therapy. In summary, our comprehensive AS perturbation landscape analysis identified RBPs as potential therapeutic targets in cancer and provided novel insights into the regulatory functions of RBPs in cancer.

Could MicroRNAs Be Useful Tools to Improve the Diagnosis and Treatment of Rare Gynecological Cancers? A Brief Overview

Gynecological cancers pose an important public health issue, with a high incidence among women of all ages. Gynecological cancers such as malignant germ-cell tumors, sex-cord-stromal tumors, uterine sarcomas and carcinosarcomas, gestational trophoblastic neoplasia, vulvar carcinoma and melanoma of the female genital tract, are defined as rare with an annual incidence of <6 per 100,000 women. Rare gynecological cancers (RGCs) are associated with poor prognosis, and given the low incidence of each entity, there is the risk of delayed diagnosis due to clinical inexperience and limited therapeutic options. There has been a growing interest in the field of microRNAs (miRNAs), a class of small non-coding RNAs of ∼22 nucleotides in length, because of their potential to regulate diverse biological processes. miRNAs usually induce mRNA degradation and translational repression by interacting with the 3' untranslated region (3'-UTR) of target mRNAs, as well as other regions and gene promoters, as well as activating translation or regulating transcription under certain conditions. Recent research has revealed the enormous promise of miRNAs for improving the diagnosis, therapy and prognosis of all major gynecological cancers. However, to date, only a few studies have been performed on RGCs. In this review, we summarize the data currently available regarding RGCs.

Bisphenol A Exposure Changes the Transcriptomic and Proteomic Dynamics of Human Retinoblastoma Y79 Cells

Bisphenol A (BPA) is a xenoestrogen chemical commonly used to manufacture polycarbonate plastics and epoxy resin and might affect various human organs. However, the cellular effects of BPA on the eyes have not been widely investigated. This study aimed to investigate the cellular cytotoxicity by BPA exposure on human retinoblastoma cells. BPA did not show cytotoxic effects, such as apoptosis, alterations to cell viability and cell cycle regulation. Comparative analysis of the transcriptome and proteome profiles were investigated after long-term exposure of Y79 cells to low doses of BPA. Transcriptome analysis using RNA-seq revealed that mRNA expression of the post-transcriptional regulation-associated gene sets was significantly upregulated in the BPA-treated group. Cell cycle regulation-associated gene sets were significantly downregulated by exposure to BPA. Interestingly, RNA-seq analysis at the transcript level indicated that alternative splicing events, particularly retained introns, were noticeably altered by low-dose BPA treatment. Additionally, proteome profiling using MALDI-TOF-MS identified a total of nine differentially expressed proteins. These results suggest that alternative splicing events and altered gene/protein expression patterns are critical phenomena affected by long-term low-dose BPA exposure. This represents a novel marker for the detection of various diseases associated with environmental pollutants such as BPA.

Alternative splicing of ceramide synthase 2 alters levels of specific ceramides and modulates cancer cell proliferation and migration in Luminal B breast cancer subtype

Global dysregulation of RNA splicing and imbalanced sphingolipid metabolism has emerged as promoters of cancer cell transformation. Here, we present specific signature of alternative splicing (AS) events of sphingolipid genes for each breast cancer subtype from the TCGA-BRCA dataset. We show that ceramide synthase 2 (CERS2) undergoes a unique cassette exon event specifically in Luminal B subtype tumors. We validated this exon 8 skipping event in Luminal B cancer cells compared to normal epithelial cells, and in patient-derived tumor tissues compared to matched normal tissues. Differential AS-based survival analysis shows that this AS event of CERS2 is a poor prognostic factor for Luminal B patients. As Exon 8 corresponds to catalytic Lag1p domain, overexpression of AS transcript of CERS2 in Luminal B cancer cells leads to a reduction in the level of very-long-chain ceramides compared to overexpression of protein-coding (PC) transcript of CERS2. We further demonstrate that this AS event-mediated decrease of very-long-chain ceramides leads to enhanced cancer cell proliferation and migration. Therefore, our results show subtype-specific AS of sphingolipid genes as a regulatory mechanism that deregulates sphingolipids like ceramides in breast tumors, and can be explored further as a suitable therapeutic target.

Identification of ribosomal protein family in triple-negative breast cancer by bioinformatics analysis

Triple-negative breast cancer (TNBC) accounts for ∼20% of all breast cancer (BC) cases. The management of TNBC represents a challenge due to its worse prognosis, heterogeneity and lack of targeted therapy. Moreover, its mechanisms are not fully clear. The aim of the study is to identify crucial genes between TNBC and non-TNBC for underlying targets for diagnostic and therapeutic methods of TNBC. The differentially expressed genes (DEGs) between TNBC and non-TNBC were selected from the Gene Expression Omnibus (GEO) database after the integrated analysis of two datasets (GSE65194 and GSE76124). Then Gene ontology (GO) and KEGG analysis were performed by DAVID database, protein-protein interaction (PPI) of DEGs was constructed by Search Tool for the Retrieval of Reciprocity Genes (STRING) database. Furthermore, centrality analysis and module analysis were carried out by Cytoscape to analyze the TNBC-related PPI. Subsequently, overall survival (OS) analysis was performed by GEPIA. Finally, the expressions of these key genes in TNBC and non-TNBC tissues were tested by qRT-PCR. The results showed that 955 DEGs were obtained, which were mainly enriched in ribosome, ribosomal subunit, and so on. Moreover, 19 candidate genes were focused on by centrality analysis and module analysis. Furthermore, we found the low expressions of ribosomal protein S9 (RPS9), ribosomal protein S14 (RPS14), ribosomal protein S27 (RPS27), ribosomal protein L11 (RPL11) and ribosomal protein L14 (RPL14) were related to a poor OS in BC patients. Additionally, qRT-PCR results suggested that these five genes were notably down-regulated in TNBC tissues. In summary, the present study suggests that ribosomal proteins are related to TNBC, and they may play an important role in the diagnosis, treatment and prognosis of TNBC.

MAGOH/MAGOHB Inhibits the Tumorigenesis of Gastric Cancer via Inactivation of b-RAF/MEK/ERK Signaling

Background

Gastric cancer is one of the most malignant tumors all over the world. It has been reported that proteins play key roles during the tumorigenesis of gastric cancer. To identify novel potential targets for gastric cancer, differential expressed proteins between gastric cancer and adjacent normal tissues were analyzed with proteomics and bioinformatics tool.

Methods

The differentially expressed proteins between gastric cancer and adjacent normal tissues were analyzed by Omicsbean (multi-omics data analysis tool). Cell viability was tested by CCK-8 assay. Flow cytometry was used to measure cell apoptosis and cycle. Transwell assay was used to test cell migration and invasion. Gene and protein expressions were detected by RT-qPCR, immunohistochemistry and Western blot, respectively.

Results

MAGOH and MAGOHB were found to be notably upregulated in gastric cancer tissues compared with that in normal tissues. Knockdown of MAGOH significantly inhibited the proliferation of gastric cancer cells via inducing the cell apoptosis. In addition, MAGOH knockdown induced G2 phase arrest in gastric cancer cells. Moreover, MAGOH knockdown notably inhibited migration and invasion of gastric cancer cells. Importantly, double knockdown of MAGOH and MAGOHB exhibited much better anti-tumor effects on gastric cancer compared with alone treatment. Finally, double knockdown of MAGOH and MAGOHB mediated the tumorigenesis of gastric cancer via regulation of RAF/MEK/ERK signaling.

Conclusion

MAGOH knockdown inhibited the tumorigenesis of gastric cancer via mediation of b-RAF/MEK/ERK signaling, and double knockdown of MAGOH and MAGOHB exhibited much better anti-tumor effects. This finding might provide us a new strategy for the treatment of gastric cancer.

SSi6 promotes cell death by apoptosis through cell cycle arrest and inhibits migration and invasion of MDA-MB-231 human breast cancer cells

Triple-negative breast cancer subtype is the most aggressive type of breast cancer due to the lack of specific therapeutic targets, having limited treatment options, low survival prognosis and high recurrence rates. In this work, we describe the effects of a semisynthetic derivative of [6]-gingerol (6G) called SSi6, produced by the addition of a 2,4-dinitrophenylhydrazine reagent on several aspects of triple-negative breast cancer biology. Human breast cancer cell lines MDA-MB-231 and MCF-10A were used in the experiments. MTT assays were used to detect cell viability. Cell cycle and apoptosis assay were analyzed using flow cytometer Accuri C6 and analysis of proteins as retinoblastoma Rb and kinases Cdk4/6 were analyzed by western blotting. SSi6 induced cytotoxic effects on triple-negative breast cancer cells, with higher selectivity when compared to the non-tumor MCF-10A cells. In addition, SSi6 inhibited migration and invasion of triple-negative breast cancer cells and was able to arrest cell cycle at the G1-phase, mainly by decreasing Cdk4/6-Rb axis levels. Therefore, SSi6 provoked the induction of apoptosis in triple-negative breast cancer cells. SSi6 was more efficient in producing these effects, compared to the original 6G natural product. This study may contribute to a better understanding of the effects of natural and semisynthetic products on the in-vitro metastatic processes in the MDA-MB-231 triple-negative breast cancer cell line. Additional, it can be useful to understand the effects of chemical modifications on already effective natural compounds aiming at the improvement of their bioactive properties, such as in the increase of the cytotoxic selectivity against tumor cells, compared to non-tumor ones.

Wide-ranging analysis of survival-related alternative splicing events in invasive breast carcinoma

Invasive breast carcinoma (BRCA) is a serious disease that threatens the survival time of those affected. Alternative splicing (AS) involved in BRCA pathogenesis may be a potential therapeutic target. However, to the best of our knowledge, a systematic analysis of survival-related alternative splicing events (SREs) has not yet been reported. The aim of the present study was to identify SREs and analyze their potential biological functions as BRCA prognostic biomarkers. An UpSet plot demonstrated AS global characteristics. Cox's proportional hazards regression model quantitatively demonstrated the prognostic relevance of AS events. Functional enrichment analysis investigated the potential pathways through which AS events affect BRCA progression. The receiver operating characteristic curve model determined the clinical significance of AS events represented using percent-spliced-in (PSI) values. The regulatory network of splicing factors (SFs) and AS events laid the foundation for studying the role of SFs in BRCA. The present study identified 1,215 SREs and their distribution characteristics, suggesting that AS events in exon skipping (ES) primarily exerted normal physiological functions, while AS events in alternative terminator sites had the most significant prognostic effect. The present study demonstrated that survival-associated genes are involved primarily in certain biological processes of ribosomal proteins. In the diagnostic model, the alternative acceptor site, alternative donor site, alternative promoter site and ES performed well. ELAVL4 was the key gene associated with prognosis and SREs. In conclusion, a number of AS events affect BRCA initiation, progression and prognosis. The PSI value of AS events has the potential to diagnose BRCA and predict a prognosis; however, this must be confirmed in additional studies.

Discrete functional and mechanistic roles of chromodomain Y-like 2 (CDYL2) transcript variants in breast cancer growth and metastasis

Rationale: Chromodomain Y-like 2 (CDYL2) is a member of the CDY gene family involved in spermatogenesis, but its role in human cancer has not been reported. Analyses of publicly available databases demonstrate that CDYL2 is abundantly expressed in breast tumors. However, whether CDYL2 is involved in breast cancer progression remains unknown.

Methods: Quantitative real-time PCR and immunoblotting assays were used to determine the expression levels of CDYL2 transcript variants in breast cancer cell lines and primary breast tumors. The effect of CDYL2 transcript variants on the malignant phenotypes of breast cancer cells was examined through in vitro and in vivo assays. Immunofluorescent staining, RNA-seq, ATAC-seq, and ChIP-qPCR were used to investigate the underlying mechanisms behind the aforementioned observations.

Results: Here we show that CDYL2 generated four transcript variants, named CDYL2a-CDYL2d. CDYL2a and CDYL2b were the predominant variants expressed in breast cancer cell lines and breast tumors and exerted strikingly discrete functions in breast cancer growth and metastasis. CDYL2a was upregulated in the majority of the breast cancer cell lines and tumors, and promoted breast cancer cell proliferation, colony formation in vitro, and tumorigenesis in xenografts. In contrast, CDYL2b was mainly expressed in luminal- and HER2-positive types of breast cancer cell lines and tumors, and suppressed the migratory, invasive, and metastatic potential of breast cancer cells in vitro and in vivo. Mechanistically, CDYL2a partially localized to SC35-positive nuclear speckles and promoted alternative splicing of a subset of target genes, including FIP1L1, NKTR, and ADD3 by exon skipping. Elimination of full-length FIP1L1, NKTR, and ADD3 rescued the impaired cell proliferation through CDYL2a depletion. In contrast, CDYL2b localized to heterochromatin and transcriptionally repressed several metastasis-promoting genes, including HPSE, HLA-F, and SELL. Restoration of HPSE, HLA-F, or SELL expression in CDYL2b-overexpressing cells attenuated the ability of CDYL2b to suppress breast cancer cell migration and invasion.

Conclusions: Collectively, these findings establish an isoform-specific function of CDYL2 in breast cancer development and progression and highlight that pharmacological inhibition of the CDYL2a, but not the CDYL2b, isoform may be an effective strategy for breast cancer therapy.

dSreg: a Bayesian model to integrate changes in splicing and RNA-binding protein activity

MOTIVATION

Alternative splicing (AS) is an important mechanism in the generation of transcript diversity across mammals. AS patterns are dynamically regulated during development and in response to environmental changes. Defects or perturbations in its regulation may lead to cancer or neurological disorders, among other pathological conditions. The regulatory mechanisms controlling AS in a given biological context are typically inferred using a two-step framework: differential AS analysis followed by enrichment methods. These strategies require setting rather arbitrary thresholds and are prone to error propagation along the analysis.

RESULTS

To overcome these limitations, we propose dSreg, a Bayesian model that integrates RNA-seq with data from regulatory features, e.g. binding sites of RNA-binding proteins. dSreg identifies the key underlying regulators controlling AS changes and quantifies their activity while simultaneously estimating the changes in exon inclusion rates. dSreg increased both the sensitivity and the specificity of the identified AS changes in simulated data, even at low read coverage. dSreg also showed improved performance when analyzing a collection of knock-down RNA-binding proteins' experiments from ENCODE, as opposed to traditional enrichment methods, such as over-representation analysis and gene set enrichment analysis. dSreg opens the possibility to integrate a large amount of readily available RNA-seq datasets at low coverage for AS analysis and allows more cost-effective RNA-seq experiments.

AVAILABILITY AND IMPLEMENTATION

dSreg was implemented in python using stan and is freely available to the community at https://bitbucket.org/cmartiga/dsreg.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

OncoOmics approaches to reveal essential genes in breast cancer: a panoramic view from pathogenesis to precision medicine

Breast cancer (BC) is the leading cause of cancer-related death among women and the most commonly diagnosed cancer worldwide. Although in recent years large-scale efforts have focused on identifying new therapeutic targets, a better understanding of BC molecular processes is required. Here we focused on elucidating the molecular hallmarks of BC heterogeneity and the oncogenic mutations involved in precision medicine that remains poorly defined. To fill this gap, we established an OncoOmics strategy that consists of analyzing genomic alterations, signaling pathways, protein-protein interactome network, protein expression, dependency maps in cell lines and patient-derived xenografts in 230 previously prioritized genes to reveal essential genes in breast cancer. As results, the OncoOmics BC essential genes were rationally filtered to 140. mRNA up-regulation was the most prevalent genomic alteration. The most altered signaling pathways were associated with basal-like and Her2-enriched molecular subtypes. RAC1, AKT1, CCND1, PIK3CA, ERBB2, CDH1, MAPK14, TP53, MAPK1, SRC, RAC3, BCL2, CTNNB1, EGFR, CDK2, GRB2, MED1 and GATA3 were essential genes in at least three OncoOmics approaches. Drugs with the highest amount of clinical trials in phases 3 and 4 were paclitaxel, docetaxel, trastuzumab, tamoxifen and doxorubicin. Lastly, we collected ~3,500 somatic and germline oncogenic variants associated with 50 essential genes, which in turn had therapeutic connectivity with 73 drugs. In conclusion, the OncoOmics strategy reveals essential genes capable of accelerating the development of targeted therapies for precision oncology.

Gene Expression Profiles Controlled by the Alternative Splicing Factor Nova2 in Endothelial Cells

Alternative splicing (AS) plays an important role in expanding the complexity of the human genome through the production of specialized proteins regulating organ development and physiological functions, as well as contributing to several pathological conditions. How AS programs impact on the signaling pathways controlling endothelial cell (EC) functions and vascular development is largely unknown. Here we identified, through RNA-seq, changes in mRNA steady-state levels in ECs caused by the neuro-oncological ventral antigen 2 (Nova2), a key AS regulator of the vascular morphogenesis. Bioinformatics analyses identified significant enrichment for genes regulated by peroxisome proliferator-activated receptor-gamma (Ppar-γ) and E2F1 transcription factors. We also showed that Nova2 in ECs controlled the AS profiles of Ppar-γ and E2F dimerization partner 2 (Tfdp2), thus generating different protein isoforms with distinct function (Ppar-γ) or subcellular localization (Tfdp2). Collectively, our results supported a mechanism whereby Nova2 integrated splicing decisions in order to regulate Ppar-γ and E2F1 activities. Our data added a layer to the sequential series of events controlled by Nova2 in ECs to orchestrate vascular biology.

Splicing Machinery is Dysregulated in Pituitary Neuroendocrine Tumors and is Associated with Aggressiveness Features

Pituitary neuroendocrine tumors (PitNETs) constitute approximately 15% of all brain tumors, and most have a sporadic origin. Recent studies suggest that altered alternative splicing and, consequently, appearance of aberrant splicing variants, is a common feature of most tumor pathologies. Moreover, spliceosome is considered an attractive therapeutic target in tumor pathologies, and the inhibition of SF3B1 (e.g., using pladienolide-B) has been shown to exert antitumor effects. Therefore, we aimed to analyze the expression levels of selected splicing-machinery components in 261 PitNETs (somatotropinomas/non-functioning PitNETS/corticotropinomas/prolactinomas) and evaluated the direct effects of pladienolide-B in cell proliferation/viability/hormone secretion in human PitNETs cell cultures and pituitary cell lines (AtT-20/GH3). Results revealed a severe dysregulation of splicing-machinery components in all the PitNET subtypes compared to normal pituitaries and a unique fingerprint of splicing-machinery components that accurately discriminate between normal and tumor tissue in each PitNET subtype. Moreover, expression of specific components was associated with key clinical parameters. Interestingly, certain components were commonly dysregulated throughout all PitNET subtypes. Finally, pladienolide-B reduced cell proliferation/viability/hormone secretion in PitNET cell cultures and cell lines. Altogether, our data demonstrate a drastic dysregulation of the splicing-machinery in PitNETs that might be associated to their tumorigenesis, paving the way to explore the use of specific splicing-machinery components as novel diagnostic/prognostic and therapeutic targets in PitNETs.

Determining the prognostic significance of alternative splicing events in soft tissue sarcoma using data from The Cancer Genome Atlas

BACKGROUND

Surgery, adjuvant chemotherapy, and radiotherapy are the primary treatment options for soft tissue sarcomas (STSs). However, identifying ways to improve the prognosis of patients with STS remains a considerable challenge. Evidence shows that the dysregulation of alternative splicing (AS) events is involved in tumor pathogenesis and progression. The present study objective was to identify survival-associated AS events that could serve as prognostic biomarkers and potentially serve as tumor-selective STS drug targets.

METHODS

STS-specific 'percent spliced in' (PSI) values for splicing events in 206 STS samples were downloaded from The Cancer Genome Atlas SpliceSeq^® database. Prognostic analyses were performed on seven types of AS events to determine their prognostic value in STS patients, for which prediction models were constructed with the risk score formula [Formula: see text]. Prediction models were also constructed to determine the prognostic value of AS events, and Spearman's rank correlation coefficients were calculated to determine the degree of correlation between splicing factor expression and the PSI values.

RESULTS

A total 10,439 events were found to significantly correlate with patient survival rates. The area under the time-dependent receiver operating characteristic curve for the prognostic predictor of STS overall survival was 0.826. Notably, the splicing events of certain STS key genes were significantly associated with STS 2-year overall survival in the present study, including exon skip (ES) events in MDM2 and EWSR1, alternate terminator events in CDKN2A and HMGA2 for dedifferentiated liposarcoma, ES in MDM2 and alternate promoter events in CDKN2A for leiomyosarcoma, and ES in EWSR1 for undifferentiated pleomorphic sarcoma. Moreover, splicing correlation networks between AS events and splicing factors revealed that almost all of the AS events showed negatively correlations with the expression of splicing factors.

CONCLUSION

An in-depth analysis of alternative RNA splicing could provide new insights into the mechanisms of STS oncogenesis and the potential for novel approaches to this type of cancer therapy.

Whole Genome Analysis and Prognostic Model Construction Based on Alternative Splicing Events in Endometrial Cancer

Objectives

A growing body of evidence has shown that aberrant alternative splicing (AS) is closely related to the occurrence and development of cancer. However, prior studies mainly have concentrated on a few genes that exhibit aberrant AS. This study aimed to determine AS events through whole genome analysis and construct a prognostic model of endometrial cancer (EC).

Methods

We downloaded gene expression RNAseq data from UCSC Xena, and seven types of AS events from TCGA SpliceSeq. Univariate Cox regression was employed to analyze the prognostic-related alternative splicing events (PASEs) and splicing factors; multivariate Cox regression was conducted to analyze the effect of risk score (All) and clinicopathological parameters on EC prognosis. An underlying interaction network of PASEs of EC was constructed by Cytoscape Reactome FI, GO, and KEGG pathway enrichment was performed by DAVID. ROC curves and Kaplan-Meir analysis were used to assess the diagnostic value of prognostic model. The correlation between PASEs and splicing factors was analyzed by GraphPad Prism; then a network was constructed using Cytoscape.

Results

In total, 28,281 AS events in EC were identified, which consisted of 1166 PASEs. RNPS1, NEK2, and CTNNB1 were the hub genes in the network of the top 600 PASEs. The area under the curve (AUC) of risk score (All) reached 0.819. Risk score (All) together with FIGO stage, cancer status, and primary therapy outcome success was risk factors of the prognosis of EC patients. Splicing factors YBX1, HNRNPDL, and HNRNPA1 were significantly related to the overall survival (OS). The splicing network indicated that the expression of splicing factors was significantly correlated with percent-splice-in (PSI) value of PASEs.

Conclusion

We constructed a model for predicting the prognosis of EC patients based on PASEs using whole genome analysis of AS events and thereby provided a reliable theoretical basis for EC clinical prognosis evaluation.

LINC00958 regulated miR-627-5p/YBX2 axis to facilitate cell proliferation and migration in oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC), the subtype of head and neck cancers, is notorious for its high incidence and death rate. The role of long non-coding RNAs (lncRNAs) is discovered to be significant for the canceration and cancer progression. Long intergenic non-protein coding RNA 958 (LINC00958) is discovered as a carcinogene in multiple cancers, such as gastric cancer, pancreatic cancer, and glioma, but there has been no report about how LINC00958 functions in OSCC. The objective of our study is to unfold function and mechanism investigation on LINC00958 in OSCC. First, TCGA database showed the upregulation and prognostic significance of LINC00958 in head and neck squamous carcinoma. Then, we discovered in OSCC clinical samples that LINC00958 presented high expression and predicted poor prognosis. Also, LINC00958 was elevated in OSCC cells. In vitro gain- and loss-function experiments proved that LINC00958 facilitated cell growth, retarded apoptosis, accelerated migration, and epithelial-to-mesenchymal transition (EMT) in OSCC. Mechanistically, we confirmed the cytoplasmic expression of LINC00958 in OSCC cells, and revealed that LINC00958 sequestered miR-627-5p to upregulate YBX2 expression. Rescue assays indicated that LINC00958 regulated OSCC cell proliferation, motility and EMT through YBX2. Together, we showed that LINC00958 promoted OSCC progression through miR-627-5p/YBX2 axis, indicating LINC00958 as a new prognostic marker, and provided new perspectives for molecular targeted treatment for OSCC.

Transcriptional consequences of impaired immune cell responses induced by cystic fibrosis plasma characterized via dual RNA sequencing

Background

In cystic fibrosis (CF), impaired immune cell responses, driven by the dysfunctional CF transmembrane conductance regulator (CFTR) gene, may determine the disease severity but clinical heterogeneity remains a major therapeutic challenge. The characterization of molecular mechanisms underlying impaired immune responses in CF may reveal novel targets with therapeutic potential. Therefore, we utilized simultaneous RNA sequencing targeted at identifying differentially expressed genes, transcripts, and miRNAs that characterize impaired immune responses triggered by CF and its phenotypes.

Methods

Peripheral blood mononuclear cells (PBMCs) extracted from a healthy donor were stimulated with plasma from CF patients (n = 9) and healthy controls (n = 3). The PBMCs were cultured (1 × 10⁵ cells/well) for 9 h at 37 ° C in 5% CO₂. After culture, total RNA was extracted from each sample and used for simultaneous total RNA and miRNA sequencing.

Results

Analysis of expression signatures from peripheral blood mononuclear cells induced by plasma of CF patients and healthy controls identified 151 genes, 154 individual transcripts, and 41 miRNAs differentially expressed in CF compared to HC while the expression signatures of 285 genes, 241 individual transcripts, and seven miRNAs differed due to CF phenotypes. Top immune pathways influenced by CF included agranulocyte adhesion, diapedesis signaling, and IL17 signaling, while those influenced by CF phenotypes included natural killer cell signaling and PI3K signaling in B lymphocytes. Upstream regulator analysis indicated dysregulation of CCL5, NF-κB and IL1A due to CF while dysregulation of TREM1 and TP53 regulators were associated with CF phenotype. Five miRNAs showed inverse expression patterns with three target genes relevant in CF-associated impaired immune pathways while two miRNAs showed inverse expression patterns with two target genes relevant to a dysregulated immune pathway associated with CF phenotypes.

Conclusions

Our results indicate that miRNAs and individual transcript variants are relevant molecular targets contributing to impaired immune cell responses in CF.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0529-0) contains supplementary material, which is available to authorized users.

Metastasis inhibition in breast cancer by targeting cancer cell extravasation

The spread of cells from primary tumors toward distant tissues and organs, also known as metastasis, is responsible for most cancer-associated deaths. The metastasis cascade comprises a series of events, characterized by the displacement of tumor cells (TCs) from the primary tumor to distant organs by traveling through the bloodstream, and their subsequent colonization. The first step in metastasis involves loss of cell-cell and cell-matrix adhesions, increased invasiveness and migratory abilities, leading to intravasation of TCs into the blood or lymphatic vessels. Stationary TCs must undergo the process of epithelial-mesenchymal transition in order to achieve this migratory and invasive phenotype. Circulating tumor cells that have survived in the circulation and left the blood or lymphatic vessels will reach distant sites where they may stay dormant for many years or grow to form secondary tumors. To do this, cells need to go through the mesenchymal-epithelial transition to revert the phenotype in order to regain epithelial cell-to-cell junctions, grow and become a clinically relevant and detectable tumor mass. This work will review the main steps of the metastatic cascade and describe some strategies to inhibit metastasis by reducing cancer cell extravasation presenting recent studies in the context of breast cancer.

A novel L1CAM isoform with angiogenic activity generated by NOVA2-mediated alternative splicing

The biological players involved in angiogenesis are only partially defined. Here, we report that endothelial cells (ECs) express a novel isoform of the cell-surface adhesion molecule L1CAM, termed L1-ΔTM. The splicing factor NOVA2, which binds directly to L1CAM pre-mRNA, is necessary and sufficient for the skipping of L1CAM transmembrane domain in ECs, leading to the release of soluble L1-ΔTM. The latter exerts high angiogenic function through both autocrine and paracrine activities. Mechanistically, L1-ΔTM-induced angiogenesis requires fibroblast growth factor receptor-1 signaling, implying a crosstalk between the two molecules. NOVA2 and L1-ΔTM are overexpressed in the vasculature of ovarian cancer, where L1-ΔTM levels correlate with tumor vascularization, supporting the involvement of NOVA2-mediated L1-ΔTM production in tumor angiogenesis. Finally, high NOVA2 expression is associated with poor outcome in ovarian cancer patients. Our results point to L1-ΔTM as a novel, EC-derived angiogenic factor which may represent a target for innovative antiangiogenic therapies.

Growing tumors stimulate the formation of new blood vessels to supply the oxygen and nutrients the cancerous cells need to stay alive. Stopping tumors from forming the blood vessels could therefore help us to treat cancer. To do so, we need to understand how different proteins control when and how blood vessels develop.

Cells make proteins by first ‘transcribing’ genes to form RNA molecules. In many cases, the RNA then goes through a process called alternative splicing. Proteins known as splicing factors cut out different segments of the RNA molecule and stick together the remaining segments to form templates for protein production. This enables a single gene to produce many different variants of a protein.

Angiolini, Belloni, Giordano et al. have now studied mouse and human versions of the cells that line the blood vessels grown by tumors. This revealed that a splicing factor called NOVA2 targets a protein called L1CAM, which is normally responsible for gluing adjacent cells together. Angiolini et al. found that NOVA2 splices L1CAM into a form not seen before. Instead of remaining anchored to cell surfaces, the newly identified form of L1CAM is released into the blood circulation, where it stimulates new blood vessels to grow.

Samples taken from the blood vessels of human ovarian tumors showed high levels of both NOVA2 and the modified form of L1CAM, while blood vessels in healthy tissue contain no, or very low levels of both proteins. Therefore, if the new form of L1CAM can be detected in the blood, it could be used to help cancer diagnosis, and to indicate which patients would benefit from treatments that restrict the growth of blood vessels in tumors. Further work is now needed to explore these possibilities.

Gliosarcoma Is Driven by Alterations in PI3K/Akt, RAS/MAPK Pathways and Characterized by Collagen Gene Expression Signature

Gliosarcoma is a very rare brain tumor reported to be a variant of glioblastoma (GBM), IDH-wildtype. While differences in molecular and histological features between gliosarcoma and GBM were reported, detailed information on the genetic background of this tumor is lacking. We intend to fill in this knowledge gap by the complex analysis of somatic mutations, indels, copy number variations, translocations and gene expression patterns in gliosarcomas. Using next generation sequencing, we determined somatic mutations, copy number variations (CNVs) and translocations in 10 gliosarcomas. Six tumors have been further subjected to RNA sequencing analysis and gene expression patterns have been compared to those of GBMs. We demonstrate that gliosarcoma bears somatic alterations in gene coding for PI3K/Akt (PTEN, PI3K) and RAS/MAPK (NF1, BRAF) signaling pathways that are crucial for tumor growth. Interestingly, the frequency of PTEN alterations in gliosarcomas was much higher than in GBMs. Aberrations of PTEN were the most frequent and occurred in 70% of samples. We identified genes differentially expressed in gliosarcoma compared to GBM (including collagen signature) and confirmed a difference in the protein level by immunohistochemistry. We found several novel translocations (including translocations in the RABGEF1 gene) creating potentially unfavorable combinations. Collected results on genetic alterations and transcriptomic profiles offer new insights into gliosarcoma pathobiology, highlight differences in gliosarcoma and GBM genetic backgrounds and point out to distinct molecular cues for targeted treatment.

Transcriptome Profiling and Molecular Therapeutic Advances in Cystic Fibrosis: Recent Insights

In cystic fibrosis (CF), mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene disrupt the capacity of the encoded protein to function as a channel to transport chloride ions and water across cell membranes. The consequences are deleterious, system-wide, and immensely variable, even among patients with the same CFTR genotype. This underscores the need to characterize the mechanisms contributing to CF pathophysiology. Gene replacement and gene editing therapies have been pursued intensively and are expected to provide a one-time treatment for CF. However, gene replacement therapy is limited by the lack of efficient vectors to deliver functional copies of CFTR to cells without immunological complications, while gene editing technologies such as CRISPR/Cas9 are still in their infancy, mainly useful in somatic cells and limited by off-target insertions. Small molecule treatments targeted at potentiating or correcting CFTR have shown clinical benefits, but they are limited to a few CFTR mutations and insufficient to overcome challenges related to clinical heterogeneity. Transcriptome profiling approaches have emerged as robust tools capable of characterizing phenotypic variability and revealing novel molecular targets with therapeutic potential for CF. We summarize current insights gained through transcriptome profiling approaches in CF studies and recent advances in molecular therapeutics.

Distinct co-expression networks using multi-omic data reveal novel interventional targets in HPV-positive and negative head-and-neck squamous cell cancer

The human papillomavirus (HPV) is present in a significant fraction of head-and-neck squamous cell cancer (HNSCC). The main goal of this study was to identify distinct co-expression patterns between HPV+ and HPV- HNSCC and to provide insights into potential regulatory mechanisms/effects within the analyzed networks. We selected cases deposited in The Cancer Genome Atlas database comprising data of gene expression, methylation profiles and mutational patterns, in addition to clinical information. The intersection among differentially expressed and differentially methylated genes showed the negative correlations between the levels of methylation and expression, suggesting that these genes have their expression levels regulated by methylation alteration patterns in their promoter. Weighted correlation network analysis was used to identify co-expression modules and a systematic approach was applied to refine them and identify key regulatory elements integrating results from the other omics. Three distinct co-expression modules were associated with HPV status and molecular signatures. Validation using independent studies reporting biological experimental data converged for the most significant genes in all modules. This study provides insights into complex genetic and epigenetic particularities in the development and progression of HNSCC according to HPV status, and contribute to unveiling specific genes/pathways as novel therapeutic targets in HNSCC.

Splicing dysregulation as a driver of breast cancer

Breast cancer is known to be a heterogeneous disease driven by a large repertoire of molecular abnormalities, which contribute to its diverse clinical behaviour. Despite the success of targeted therapy approaches for breast cancer patient management, there is still a lack of the molecular understanding of aggressive forms of the disease and clinical management of these patients remains difficult. The advent of high-throughput sequencing technologies has paved the way for a more complete understanding of the molecular make-up of the breast cancer genome. As such, it is becoming apparent that disruption of canonical splicing within breast cancer governs its clinical progression. In this review, we discuss the role of dysregulation of spliceosomal component genes and associated factors in the progression of breast cancer, their role in therapy resistance and the use of quantitative isoform expression as potential prognostic and predictive biomarkers with a particular focus on oestrogen receptor-positive breast cancer.

NF1 deficiency correlates with estrogen receptor signaling and diminished survival in breast cancer

The key negative regulatory gene of the RAS pathway, NF1, is mutated or deleted in numerous cancer types and is associated with increased cancer risk and drug resistance. Even though women with neurofibromatosis (germline NF1 mutations) have a substantially increased breast cancer risk at a young age and NF1 is commonly mutated in sporadic breast cancers, we have a limited understanding of the role of NF1 in breast cancer. We utilized CRISPR-Cas9 gene editing to create Nf1 rat models to evaluate the effect of Nf1 deficiency on tumorigenesis. The resulting Nf1 indels induced highly penetrant, aggressive mammary adenocarcinomas that express estrogen receptor (ER) and progesterone receptor (PR). We identified distinct Nf1 mRNA and protein isoforms that were altered during tumorigenesis. To evaluate NF1 in human breast cancer, we analyzed genomic changes in a data set of 2000 clinically annotated breast cancers. We found NF1 shallow deletions in 25% of sporadic breast cancers, which correlated with poor clinical outcome. To identify biological networks impacted by NF1 deficiency, we constructed gene co-expression networks using weighted gene correlation network analysis (WGCNA) and identified a network connected to ESR1 (estrogen receptor). Moreover, NF1-deficient cancers correlated with established RAS activation signatures. Estrogen-dependence was verified by estrogen-ablation in Nf1 rats where rapid tumor regression was observed. Additionally, Nf1 deficiency correlated with increased estrogen receptor phosphorylation in mammary adenocarcinomas. These results demonstrate a significant role for NF1 in both NF1-related breast cancer and sporadic breast cancer, and highlight a potential functional link between neurofibromin and the estrogen receptor.

Perspective on Oncogenic Processes at the End of the Beginning of Cancer Genomics

The Cancer Genome Atlas (TCGA) has catalyzed systematic characterization of diverse genomic alterations underlying human cancers. At this historic junction marking the completion of genomic characterization of over 11,000 tumors from 33 cancer types, we present our current understanding of the molecular processes governing oncogenesis. We illustrate our insights into cancer through synthesis of the findings of the TCGA PanCancer Atlas project on three facets of oncogenesis: (1) somatic driver mutations, germline pathogenic variants, and their interactions in the tumor; (2) the influence of the tumor genome and epigenome on transcriptome and proteome; and (3) the relationship between tumor and the microenvironment, including implications for drugs targeting driver events and immunotherapies. These results will anchor future characterization of rare and common tumor types, primary and relapsed tumors, and cancers across ancestry groups and will guide the deployment of clinical genomic sequencing.

hnRNPM guides an alternative splicing program in response to inhibition of the PI3K/AKT/mTOR pathway in Ewing sarcoma cells

Ewing sarcomas (ES) are biologically aggressive tumors of bone and soft tissues for which no cure is currently available. Most ES patients do not respond to chemotherapeutic treatments or acquire resistance. Since the PI3K/AKT/mTOR axis is often deregulated in ES, its inhibition offers therapeutic perspective for these aggressive tumors. Herein, by using splicing sensitive arrays, we have uncovered an extensive splicing program activated upon inhibition of the PI3K/AKT/mTOR signaling pathway by BEZ235. Bioinformatics analyses identified hnRNPM as a key factor in this response. HnRNPM motifs were significantly enriched in introns flanking the regulated exons and proximity of binding represented a key determinant for hnRNPM-dependent splicing regulation. Knockdown of hnRNPM expression abolished a subset of BEZ235-induced splicing changes that contained hnRNPM binding sites, enhanced BEZ235 cytotoxicity and limited the clonogenicity of ES cells. Importantly, hnRNPM up-regulation correlates with poor outcome in sarcoma patients. These findings uncover an hnRNPM-dependent alternative splicing program set in motion by inhibition of the mTOR/AKT/PI3K pathway in ES cells that limits therapeutic efficacy of pharmacologic inhibitors, suggesting that combined inhibition of the PI3K/AKT/mTOR pathway and hnRNPM activity may represent a novel approach for ES treatment.

The clinicopathological significance of UBE2C in breast cancer: a study based on immunohistochemistry, microarray and RNA-sequencing data

BACKGROUND

Ubiquitin-conjugating enzyme E2C (UBE2C) has been previously reported to correlate with the malignant progression of various human cancers, however, the exact molecular function of UBE2C in breast carcinoma (BRCA) remained elusive. We aimed to investigate UBE2C expression in BRCA and its clinical significance.

METHODS

The expression of UBE2C in 209 BRCA tissue samples and 53 adjacent normal tissue samples was detected using immunohistochemistry. The clinical role of UBE2C was analyzed. Public databases including the human protein atlas and Oncomine were used to assess UBE2C expression in BRCA. Moreover, the cancer genome atlas (TCGA) database was employed to investigate the prognostic value of UBE2C in BRCA.

RESULTS

The positive expression rate of UBE2C in BRCA was 70.8% (148/209), and UBE2C expression in the adjacent breast tissue was negative. The expression of UBE2C was positively correlated with tumor size (r = 0.32, P < 0.001), histological grade (r = 0.237, P = 0.001), clinical stage (r = 0.198, P = 0.004), lymph node metastasis (r = 0.155, P = 0.026), HER2 expression level (r = 0.356, P < 0.001), Ki-67 expression level (r = 0.504, P < 0.001), and P53 expression level (r = 0.32, P = 0.001). Negative correlations were found between UBE2C expression and the ER (r = - 0.403, P < 0.001) and PR (r = - 0.468, P < 0.001) status. UBE2C gene expression data from the public databases all proved that UBE2C was overexpressed in BRCA. According to the TCGA data analysis, a higher positive expression of UBE2C was associated with worse survival of BRCA patients (P = 0.0428), and data from cBioPortal indicated that 11% of all sequenced BRCA patients possessed a gene alteration of UBE2C, predominately gene amplification and mRNA regulation.

CONCLUSION

Ubiquitin-conjugating enzyme E2C might pose an oncogenic effect on the progression of BRCA.