Alternative splicing and biological heterogeneity in prostate cancer

Aberrant PI3Kδ splice isoform as a potential biomarker and novel therapeutic target for endocrine cancers

Introduction

PI3K/AKT signaling pathway is upregulated in a broad spectrum of cancers. Among the class I PI3Ks (PI3Kδ/β/δ isoforms), PI3Kδ has been implicated in hematologic cancers and solid tumors. Alternative splicing is a post-transcriptional process for acquiring proteomic diversity in eukaryotic cells. Emerging evidence has highlighted the involvement of aberrant mRNA splicing in cancer development/progression.

Methods

Our previous studies revealed that PIK3CD-S is an oncogenic splice variant that promotes tumor aggressiveness and drug resistance in prostate cancer (PCa). To further evaluate the potential of utilizing PI3Kδ-S (encoded from PIK3CD-S) as a cancer biomarker and/or drug target, comprehensive analyses were performed in a series of patient samples and cell lines derived from endocrine/solid tumors. Specifically, IHC, immunofluorescence, western blot and RT-PCR assay results have demonstrated that PI3Kδ isoforms were highly expressed in endocrine/solid tumor patient specimens and cell lines.

Results

Differential PIK3CD-S/PIK3CD-L expression profiles were identified in a panel of endocrine/solid tumor cells. SiRNA knockdown of PIK3CD-L or PIK3CD-S differentially inhibits AKT/mTOR signaling in PCa, breast, colon and lung cancer cell lines. Moreover, siRNA knockdown of PTEN increased PI3Kδ levels and activated AKT/mTOR signaling, while overexpression of PTEN reduced PI3Kδ levels and inhibited AKT/mTOR signaling in cancer cells. Intriguingly, PI3Kδ-S levels remained unchanged upon either siRNA knockdown or overexpression of PTEN. Taken together, these results suggested that PTEN negatively regulates PI3Kδ-L and its downstream AKT/mTOR signaling, while PI3Kδ-S promotes AKT/mTOR signaling without regulation by PTEN. Lastly, PI3Kδ inhibitor Idelalisib and SRPK1/2 inhibitor SRPIN340 were employed to assess their efficacies on inhibiting the PI3Kδ-expressing endocrine/solid tumors. Our results have shown that Idelalisib effectively inhibited PI3Kδ-L (but not PI3Kδ-S) mediated AKT/mTOR signaling. In contrast, SRPIN340 reversed the aberrant mRNA splicing, thereby inhibiting AKT/mTOR signaling. In-vitro functional assays have further demonstrated that a combination of Idelalisib and SRPIN340 achieved a synergistic drug effect (with drastically reduced cell viabilities/growths of tumor spheroids) in inhibiting the advanced tumor cells.

Conclusion

In summary, our study has suggested a promising potential of utilizing PI3Kδ-S (an oncogenic isoform conferring drug resistance and exempt from PTEN regulation) as a prognostic biomarker and drug target in advanced endocrine cancers.

A Splicing Transcriptome-Wide Association Study Identifies Candidate Altered Splicing for Prostate Cancer Risk

Prostate cancer (PCa) represents a huge public health burden among men. Many susceptibility genetic factors for PCa still remain unknown. In this study, we performed a large splicing transcriptome-wide association study (spTWAS) using three modeling strategies to develop alternative splicing genetic prediction models for identifying novel susceptibility loci and splicing introns for PCa risk by assessing 79,194 cases and 61,112 controls of European ancestry in the PRACTICAL, CRUK, CAPS, BPC3, and PEGASUS consortia. We identified 120 splicing introns of 97 genes showing an association with PCa risk at false discovery rate (FDR)-corrected threshold (FDR <0.05). Of them, 33 genes were enriched in PCa-related diseases and function categories. Fine-mapping analysis suggested that 21 splicing introns of 19 genes were likely causally associated with PCa risk. Thirty-five splicing introns of 34 novel genes were identified to be related to PCa susceptibility for the first time, and 11 of the genes were enriched in a cancer-related network. Our study identified novel loci and splicing introns associated with PCa risk, which can improve our understanding of the etiology of this common malignancy.

Identification and characterization of novel ETV4 splice variants in prostate cancer

ETV4, one of ETS proteins overexpressed in prostate cancer, promotes migration, invasion, and proliferation in prostate cells. This study identifies a series of previously unknown ETV4 alternatively spliced transcripts in human prostate cell lines. Their expression has been validated using several unbiased techniques, including Nanopore sequencing. Most of these transcripts originate from an in-frame exon skipping and, thus, are expected to be translated into ETV4 protein isoforms. Functional analysis of the most abundant among these isoforms shows that they still bear an activity, namely a reduced ability to promote proliferation and a residual ability to regulate the transcription of ETV4 target genes. Alternatively spliced genes are common in cancer cells: an analysis of the TCGA dataset confirms the abundance of these novel ETV4 transcripts in prostate tumors, in contrast to peritumoral tissues. Since none of their translated isoforms have acquired a higher oncogenic potential, such abundance is likely to reflect the tumor deranged splicing machinery. However, it is also possible that their interaction with the canonical variants may contribute to the biology and the clinics of prostate cancer. Further investigations are needed to elucidate the biological role of these ETV4 transcripts and of their putative isoforms.

Molecular Insight into Drug Resistance Mechanism Conferred by Aberrant PIK3CD Splice Variant in African American Prostate Cancer

Targeting PI3Kδ has emerged as a promising therapy for hematologic and non-hematologic malignancies. Previously, we identified an oncogenic splice variant, PIK3CD-S, conferring Idelalisib resistance in African American (AA) prostate cancer (PCa). In the current study, we employed a comprehensive analysis combining molecular biology, biochemistry, histology, in silico simulation, and in vitro functional assays to investigate the PIK3CD-S expression profiles in PCa samples and to elucidate the drug resistance mechanism mediated by PI3Kδ-S (encoded by PIK3CD-S). The immunohistochemistry, RT-PCR, and Western blot assays first confirmed that PI3Kδ-S is highly expressed in AA PCa. Compared with PCa expressing the full-length PI3Kδ-L, PCa expressing PI3Kδ-S exhibits enhanced drug resistance properties, including a higher cell viability, more antiapoptotic and invasive capacities, and constitutively activated PI3K/AKT signaling, in the presence of PI3Kδ/PI3K inhibitors (Idelalisib, Seletalisib, Wortmannin, and Dactolisib). Molecular docking, ATP-competitive assays, and PI3 kinase assays have further indicated a drastically reduced affinity of PI3Kδ inhibitors with PI3Kδ-S vs. PI3Kδ-L, attributed to the lack of core binding residues in the PI3Kδ-S catalytic domain. Additionally, SRSF2 has been identified as a critical splicing factor mediating exon 20 skipping in PIK3CD pre-mRNA. The inhibition of the SRSF2 activity by SRPIN340 successfully sensitizes AA PCa cells to PI3Kδ inhibitors, suggesting a novel therapeutic option for Idelalisib-resistant tumors.

FOXA1 regulates alternative splicing in prostate cancer

Dysregulation of alternative splicing in prostate cancer is linked to transcriptional programs activated by AR, ERG, FOXA1, and MYC. Here, we show that FOXA1 functions as the primary orchestrator of alternative splicing dysregulation across 500 primary and metastatic prostate cancer transcriptomes. We demonstrate that FOXA1 binds to the regulatory regions of splicing-related genes, including HNRNPK and SRSF1. By controlling trans-acting factor expression, FOXA1 exploits an "exon definition" mechanism calibrating alternative splicing toward dominant isoform production. This regulation especially impacts splicing factors themselves and leads to a reduction of nonsense-mediated decay (NMD)-targeted isoforms. Inclusion of the NMD-determinant FLNA exon 30 by FOXA1-controlled oncogene SRSF1 promotes cell growth in vitro and predicts disease recurrence. Overall, we report a role for FOXA1 in rewiring the alternative splicing landscape in prostate cancer through a cascade of events from chromatin access, to splicing factor regulation, and, finally, to alternative splicing of exons influencing patient survival.

DeepIII: Predicting Isoform-Isoform Interactions by Deep Neural Networks and Data Fusion

Alternative splicing enables a gene translating into different isoforms and into the corresponding proteoforms, which actually accomplish various biological functions of a living body. Isoform-isoform interactions (IIIs) provide a higher resolution interactome to explore the cellular processes and disease mechanisms than the canonically studied protein-protein interactions (PPIs), which are often recorded at the coarse gene level. The knowledge of IIIs is critical to map pathways, understand protein complexity and functional diversity, but the known IIIs are very scanty. In this paper, we propose a deep learning based method called DeepIII to systematically predict genome-wide IIIs by integrating diverse data sources, including RNA-seq datasets of different human tissues, exon array data, domain-domain interactions (DDIs) of proteins, nucleotide sequences and amino acid sequences. Particularly, DeepIII fuses these data to learn the representation of isoform pairs with a four-layer deep neural networks, and then performs binary classification on the learnt representation to achieve the prediction of IIIs. Experimental results show that DeepIII achieves a superior prediction performance to the state-of-the-art solutions and the III network constructed by DeepIII gives more accurate isoform function prediction. Case studies further confirm that DeepIII can differentiate the individual interaction partners of different isoforms spliced from the same gene. The code and datasets of DeepIII are available at http://mlda.swu.edu.cn/codes.php?name=DeepIII.

Evolution of the Human Brain Can Help Determine Pathophysiology of Neurodevelopmental Disorders

The evolution of humans brought about a co-occurring evolution of the human brain, which is far larger and more complex than that of many other organisms. The brain has evolved characteristically in humans in many respects, including macro-and micro-anatomical changes in the brain structure, changes in gene expression, and cell populations and ratios. These characteristics are essential for the execution of higher functions, such as sociality, language, and cognition, which express humanity, and are thought to have been acquired over evolutionary time. However, with the acquisition of higher functions also comes the risk of the disease in which they fail. This review focuses on human brain evolution and neurodevelopmental disorders (NDDs) and discusses brain development, molecular evolution, and human brain evolution. Discussing the potential for the development and pathophysiology of NDDs acquired by human brain evolution will provide insights into the acquisition and breakdown of higher functions from a new perspective.

The aberrant upregulation of exon 10-inclusive SREK1 through SRSF10 acts as an oncogenic driver in human hepatocellular carcinoma

Deregulation of alternative splicing is implicated as a relevant source of molecular heterogeneity in cancer. However, the targets and intrinsic mechanisms of splicing in hepatocarcinogenesis are largely unknown. Here, we report a functional impact of a Splicing Regulatory Glutamine/Lysine-Rich Protein 1 (SREK1) variant and its regulator, Serine/arginine-rich splicing factor 10 (SRSF10). HCC patients with poor prognosis express higher levels of exon 10-inclusive SREK1 (SREK1^L). SREK1^L can sustain BLOC1S5-TXNDC5 (B-T) expression, a targeted gene of nonsense-mediated mRNA decay through inhibiting exon-exon junction complex binding with B-T to exert its oncogenic role. B-T plays its competing endogenous RNA role by inhibiting miR-30c-5p and miR-30e-5p, and further promoting the expression of downstream oncogenic targets SRSF10 and TXNDC5. Interestingly, SRSF10 can act as a splicing regulator for SREK1^L to promote hepatocarcinogenesis via the formation of a SRSF10-associated complex. In summary, we demonstrate a SRSF10/SREK1^L/B-T signalling loop to accelerate the hepatocarcinogenesis.

Alternative splicing is dysregulated in hepatocellular carcinoma. Here, the authors investigate the role of the splice variant of Splicing Regulatory Glutamic Acid and Lysine Rich Protein 1 (SREK1) and its upstream regulator, Serine/arginine-rich splicing factor 10 (SRSF10) in sustaining the oncogenic signal.

Summit: Automated Analysis of Arrayed Single-Cell Gel Electrophoresis

New pipelines are required to automate the quantitation of emerging high-throughput electrophoretic (EP) assessment of DNA damage, or proteoform expression in single cells. EP cytometry consists of thousands of Western blots performed on a microscope slide-sized gel microwell array for single cells. Thus, EP cytometry images pose an analysis challenge that blends requirements for accurate and reproducible analysis encountered for both standard Western blots and protein microarrays. Here, we introduce the Summit algorithm to automate array segmentation, peak background subtraction, and Gaussian fitting for EP cytometry. The data structure storage of parameters allows users to perform quality control on identically processed data, yielding a ~6.5% difference in coefficient of quartile variation (CQV) of protein peak area under the curve (AUC) distributions measured by four users. Further, inspired by investigations of background subtraction methods to reduce technical variation in protein microarray measurements, we aimed to understand the trade-offs between EP cytometry analysis throughput and variation. We found an 11%-50% increase in protein peaks that passed quality control with a subtraction method similar to microarray "average on-boundary" versus an axial subtraction method. The background subtraction method only mildly influences AUC CQV, which varies between 1% and 4.5%. Finally, we determined that the narrow confidence interval for peak location and peak width parameters from Gaussian fitting yield minimal uncertainty in protein sizing. The AUC CQV differed by only ~1%-2% when summed over the peak width bounds versus the 95% peak width confidence interval. We expect Summit to be broadly applicable to other arrayed EP separations, or traditional Western blot analysis.

Androgen receptor and its splice variant, AR-V7, differentially induce mRNA splicing in prostate cancer cells

Prostate cancer (PCa) is dependent on the androgen receptor (AR). Advanced PCa is treated with an androgen deprivation therapy-based regimen; tumors develop resistance, although they typically remain AR-dependent. Expression of constitutively active AR variants lacking the ligand-binding domain including the variant AR-V7 contributes to this resistance. AR and AR-V7, as transcription factors, regulate many of the same genes, but also have unique activities. In this study, the capacity of the two AR isoforms to regulate splicing was examined. RNA-seq data from models that endogenously express AR and express AR-V7 in response to doxycycline were used. Both AR isoforms induced multiple changes in splicing and many changes were isoform-specific. Analyses of two endogenous genes, PGAP2 and TPD52, were performed to examine differential splicing. A novel exon that appears to be a novel transcription start site was preferentially induced by AR-V7 in PGAP2 although it is induced to a lesser extent by AR. The previously described AR induced promoter 2 usage that results in a novel protein derived from TPD52 (PrLZ) was not induced by AR-V7. AR, but not AR-V7, bound to a site proximal to promoter 2, and induction was found to depend on FOXA1.

Age-related injury responses of human oligodendrocytes to metabolic insults: link to BCL-2 and autophagy pathways

Myelin destruction and oligodendrocyte (OL) death consequent to metabolic stress is a feature of CNS disorders across the age spectrum. Using cells derived from surgically resected tissue, we demonstrate that young (<age 5) pediatric-aged sample OLs are more resistant to in-vitro metabolic injury than fetal O4+ progenitor cells, but more susceptible to cell death and apoptosis than adult-derived OLs. Pediatric but not adult OLs show measurable levels of TUNEL+ cells, a feature of the fetal cell response. The ratio of anti- vs pro-apoptotic BCL-2 family genes are increased in adult vs pediatric (<age 5) mature OLs and in more mature OL lineage cells. Lysosomal gene expression was increased in adult and pediatric compared to fetal OL lineage cells. Cell death of OLs was increased by inhibiting pro-apoptotic BCL-2 gene and autophagy activity. These distinct age-related injury responses should be considered in designing therapies aimed at reducing myelin injury.

A Roadmap Toward the Definition of Actionable Tumor-Specific Antigens

The search for tumor-specific antigens (TSAs) has considerably accelerated during the past decade due to the improvement of proteogenomic detection methods. This provides new opportunities for the development of novel antitumoral immunotherapies to mount an efficient T cell response against one or multiple types of tumors. While the identification of mutated antigens originating from coding exons has provided relatively few TSA candidates, the possibility of enlarging the repertoire of targetable TSAs by looking at antigens arising from non-canonical open reading frames opens up interesting avenues for cancer immunotherapy. In this review, we outline the potential sources of TSAs and the mechanisms responsible for their expression strictly in cancer cells. In line with the heterogeneity of cancer, we propose that discrete families of TSAs may be enriched in specific cancer types.

Elevated KLF7 levels may serve as a prognostic signature and might contribute to progression of squamous carcinoma

Global efforts have been undertaken to define the genome-wide distribution of epigenetic markers in cancerous tissues, which provide an invaluable opportunity to understand cancer biology and identify predictive signatures. Several studies have focused on the gene expression patterns of squamous carcinoma to identify tumor subtypes and find prognostic and therapeutic targets because squamous carcinoma genomes showed high instability. However, the number of reliable reports referring prognostic significance of genes and their role in squamous carcinoma is still quite limited. Krüppel-like factor 7 (KLF7) is a transcription factor that is widely expressed in numerous human tissues at low levels. Members of the KLF family have established roles in tumor cell fate, stress response, cell survival and the tumor-initiating properties of cancer stem-like cells. Hence to investigate whether KFL7 expression from cancer tissue holds promise as a prognostic and/or therapeutic target, we analyzed gene expression profiles from squamous carcinoma and surgical margin tissues in The Cancer Genome Atlas. We identified significant up-regulation of KLF7 in squamous carcinoma, which was confirmed by immunohistochemical staining. Elevated KLF7 expression was associated with poor squamous carcinoma prognosis before and after correcting for confounding factors by multivariate Cox regression analysis. Several pathways, such as Neurotrophin and GnRH pathways, were activated in KLF7-up-regulated squamous carcinoma samples through Gene Set Enrichment Analysis. In conclusion, we consolidate the potential role(s) of KLF7 in squamous carcinoma carcinogenesis from The Cancer Genome Atlas surgical margin tissue, offering insights into expression signatures that are potentially useful for prognosis modalities.

DSCAM-AS1-Driven Proliferation of Breast Cancer Cells Involves Regulation of Alternative Exon Splicing and 3'-End Usage

DSCAM-AS1 is a cancer-related long noncoding RNA with higher expression levels in Luminal A, B, and HER2-positive Breast Carcinoma (BC), where its expression is strongly dependent on Estrogen Receptor Alpha (ERα). DSCAM-AS1 expression is analyzed in 30 public datasets and, additionally, by qRT-PCR in tumors from 93 BC patients, to uncover correlations with clinical data. Moreover, the effect of DSCAM-AS1 knockdown on gene expression and alternative splicing is studied by RNA-Seq in MCF-7 cells. We confirm DSCAM-AS1 overexpression in high grade Luminal A, B, and HER2+ BCs and find a significant correlation with disease relapse. In total, 908 genes are regulated by DSCAM-AS1-silencing, primarily involved in the cell cycle and inflammatory response. Noteworthily, the analysis of alternative splicing and isoform regulation reveals 2085 splicing events regulated by DSCAM-AS1, enriched in alternative polyadenylation sites, 3′UTR (untranslated region) shortening and exon skipping events. Finally, the DSCAM-AS1-interacting splicing factor heterogeneous nuclear ribonucleoprotein L (hnRNPL) is predicted as the most enriched RBP for exon skipping and 3′UTR events. The relevance of DSCAM-AS1 overexpression in BC is confirmed by clinical data and further enhanced by its possible involvement in the regulation of RNA processing, which is emerging as one of the most important dysfunctions in cancer.

SpliceHetero: An information theoretic approach for measuring spliceomic intratumor heterogeneity from bulk tumor RNA-seq

Motivation

Intratumor heterogeneity (ITH) represents the diversity of cell populations that make up cancer tissue. The level of ITH in a tumor is usually measured by a genomic variation profile, such as copy number variation and somatic mutation. However, a recent study has identified ITH at the transcriptome level and suggested that ITH at gene expression levels is useful for predicting prognosis. Measuring ITH levels at the spliceome level is a natural extension. There are serious technical challenges in measuring spliceomic ITH (sITH) from bulk tumor RNA sequencing (RNA-seq) due to the complex splicing patterns.

Results

We propose an information-theoretic method to measure the sITH of bulk tumors to overcome the above challenges. This method has been extensively tested in experiments using synthetic data, xenograft tumor data, and TCGA pan-cancer data. As a result, we showed that sITH is closely related to cancer progression and clonal heterogeneity, along with clinically significant features such as cancer stage, survival outcome and PAM50 subtype. As far as we know, it is the first study to define ITH at the spliceome level. This method can greatly improve the understanding of cancer spliceome and has great potential as a diagnostic and prognostic tool.

Comprehensive investigation of alternative splicing and development of a prognostic risk score for prostate cancer based on six-gene signatures

Purpose: To systematically document alternative splicing profiles of prostate cancer in relatively large populations in order to construct a prognostic predictors model for prostate cancer.

Methods: Splicing data and clinical information of 495 prostate cancer patients were obtained from The Cancer Genome Atlas (TCGA). The SpliceSeq database was used to extract information regarding splicing events. Multiple bioinformatic tools were used for functional and pathway enrichment analysis as well as for construction of gene interaction networks. Candidate gene expression profiles were verified with clinical samples using QRT-PCR.

Results: We detected a total of 44070 alternative splicing events of 10381 genes in prostate cancer. 7 and 14 KEGG pathways were enriched and were associated with overall and recurrence-free survival, respectively. The expression of 396 genes among the 1526 overall survival genes associated alternative splicing events were associated with overall survival. The expression of 483 genes among the 1916 recurrence-free survival genes associated alternative splicing events were associated with recurrence-free survival. Lastly, we constructed the prognosis risk score system based on the expression profiles of six-gene signatures which in combination had an AUC of 0.941 for overall survival associated alternative splicing events, followed by overall survival associated gene expressions with an AUC of 0.794, a recurrence-free survival associated gene expression with an AUC of 0.752 and recurrence-free survival associated alternative splicing events with an AUC of 0.735, indicating its strong ability to predict patient outcome. The expression profile of the six genes was also confirmed in different prostate cell lines and clinic samples.

Conclusion: Our comprehensive investigation of alternative splicing not only provided insight into the biological pathways of alternative splicing involved in the development of prostate cancer but also revealed new potential biomarkers for prognosticating as well as novel therapeutic targets for development of prostate cancer treatment.

A Novel FGFR3 Splice Variant Preferentially Expressed in African American Prostate Cancer Drives Aggressive Phenotypes and Docetaxel Resistance

Alternative splicing (AS) has been shown to participate in prostate cancer development and progression; however, a link between AS and prostate cancer health disparities has been largely unexplored. Here we report on the cloning of a novel splice variant of FGFR3 that is preferentially expressed in African American (AA) prostate cancer. This novel variant (FGFR3-S) omits exon 14, comprising 123 nucleotides that encode the activation loop in the intracellular split kinase domain. Ectopic overexpression of FGFR3-S in European American (EA) prostate cancer cell lines (PC-3 and LNCaP) led to enhanced receptor autophosphorylation and increased activation of the downstream signaling effectors AKT, STAT3, and ribosomal S6 compared with FGFR3-L (retains exon 14). The increased oncogenic signaling imparted by FGFR3-S was associated with a substantial gain in proliferative and antiapoptotic activities, as well as a modest but significant gain in cell motility. Moreover, the FGFR3-S-conferred proliferative and motility gains were highly resistant to the pan-FGFR small-molecule inhibitor dovitinib and the antiapoptotic gain was insensitive to the cytotoxic drug docetaxel, which stands in marked contrast with dovitinib- and docetaxel-sensitive FGFR3-L. In an in vivo xenograft model, mice injected with PC-3 cells overexpressing FGFR3-S exhibited significantly increased tumor growth and resistance to dovitinib treatment compared with cells overexpressing FGFR3-L. In agreement with our in vitro and in vivo findings, a high FGFR3-S/FGFR3-L expression ratio in prostate cancer specimens was associated with poor patient prognosis. IMPLICATIONS: This work identifies a novel FGFR3 splice variant and supports the hypothesis that differential AS participates in prostate cancer health disparities.

Dysregulated Transcriptional Control in Prostate Cancer

Recent advances in whole-genome and transcriptome sequencing of prostate cancer at different stages indicate that a large number of mutations found in tumors are present in non-protein coding regions of the genome and lead to dysregulated gene expression. Single nucleotide variations and small mutations affecting the recruitment of transcription factor complexes to DNA regulatory elements are observed in an increasing number of cases. Genomic rearrangements may position coding regions under the novel control of regulatory elements, as exemplified by the TMPRSS2-ERG fusion and the amplified enhancer identified upstream of the androgen receptor (AR) gene. Super-enhancers are increasingly found to play important roles in aberrant oncogenic transcription. Several players involved in these processes are currently being evaluated as drug targets and may represent new vulnerabilities that can be exploited for prostate cancer treatment. They include factors involved in enhancer and super-enhancer function such as bromodomain proteins and cyclin-dependent kinases. In addition, non-coding RNAs with an important gene regulatory role are being explored. The rapid progress made in understanding the influence of the non-coding part of the genome and of transcription dysregulation in prostate cancer could pave the way for the identification of novel treatment paradigms for the benefit of patients.

Role of Alternative Splicing in Prostate Cancer Aggressiveness and Drug Resistance in African Americans

Alternative splicing, the process of removing introns and joining exons of pre-mRNA, is critical for growth, development, tissue homeostasis, and species diversity. Dysregulation of alternative splicing can initiate and drive disease. Aberrant alternative splicing has been shown to promote the "hallmarks of cancer" in both hematological and solid cancers. Of interest, recent work has focused on the role of alternative splicing in prostate cancer and prostate cancer health disparities. We will provide a review of prostate cancer health disparities involving the African American population, alternative RNA splicing, and alternative splicing in prostate cancer. Lastly, we will summarize our work on differential alternative splicing in prostate cancer disparities and its implications for disparate health outcomes and therapeutic targets.

Androgen-regulated transcription of ESRP2 drives alternative splicing patterns in prostate cancer

Prostate is the most frequent cancer in men. Prostate cancer progression is driven by androgen steroid hormones, and delayed by androgen deprivation therapy (ADT). Androgens control transcription by stimulating androgen receptor (AR) activity, yet also control pre-mRNA splicing through less clear mechanisms. Here we find androgens regulate splicing through AR-mediated transcriptional control of the epithelial-specific splicing regulator ESRP2. Both ESRP2 and its close paralog ESRP1 are highly expressed in primary prostate cancer. Androgen stimulation induces splicing switches in many endogenous ESRP2-controlled mRNA isoforms, including splicing switches correlating with disease progression. ESRP2 expression in clinical prostate cancer is repressed by ADT, which may thus inadvertently dampen epithelial splice programmes. Supporting this, treatment with the AR antagonist bicalutamide (Casodex) induced mesenchymal splicing patterns of genes including FLNB and CTNND1. Our data reveals a new mechanism of splicing control in prostate cancer with important implications for disease progression.

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.

Androgen-dependent alternative mRNA isoform expression in prostate cancer cells

Background: Androgen steroid hormones are key drivers of prostate cancer. Previous work has shown that androgens can drive the expression of alternative mRNA isoforms as well as transcriptional changes in prostate cancer cells. Yet to what extent androgens control alternative mRNA isoforms and how these are expressed and differentially regulated in prostate tumours is unknown. Methods: Here we have used RNA-Seq data to globally identify alternative mRNA isoform expression under androgen control in prostate cancer cells, and profiled the expression of these mRNA isoforms in clinical tissue. Results: Our data indicate androgens primarily switch mRNA isoforms through alternative promoter selection. We detected 73 androgen regulated alternative transcription events, including utilisation of 56 androgen-dependent alternative promoters, 13 androgen-regulated alternative splicing events, and selection of 4 androgen-regulated alternative 3' mRNA ends. 64 of these events are novel to this study, and 26 involve previously unannotated isoforms. We validated androgen dependent regulation of 17 alternative isoforms by quantitative PCR in an independent sample set. Some of the identified mRNA isoforms are in genes already implicated in prostate cancer (including LIG4, FDFT1 and RELAXIN), or in genes important in other cancers (e.g. NUP93 and MAT2A). Importantly, analysis of transcriptome data from 497 tumour samples in the TGCA prostate adenocarcinoma (PRAD) cohort identified 13 mRNA isoforms (including TPD52, TACC2 and NDUFV3) that are differentially regulated in localised prostate cancer relative to normal tissue, and 3 ( OSBPL1A, CLK3 and TSC22D3) which change significantly with Gleason grade and  tumour stage. Conclusions: Our findings dramatically increase the number of known androgen regulated isoforms in prostate cancer, and indicate a highly complex response to androgens in prostate cancer cells that could be clinically important.

Discovery and expression analysis of novel transcripts of the human SR-related CTD-associated factor 1 (SCAF1) gene in human cancer cells using Next-Generation Sequencing

The human SR-related CTD associated factor 1 (SCAF1) gene is a new member of the human SR (Ser/Arg-rich) superfamily of pre-mRNA splicing factors, which has been discovered and cloned by members of our lab. SCAF1 interacts with the CTD domain of the RNA polymerase II polypeptide A and is firmly involved in pre-mRNA splicing. Although it was found to be expressed widely in multiple human tissues, its mRNA levels vary a lot. The significant relation of SCAF1 with cancer has been confirmed by many studies, since SCAF1 mRNA transcript was found to be overexpressed in breast and ovarian tumors, confirming its significant prognostic value as a cancer biomarker in both these human malignancies. In this study, we describe the discovery and cloning of fifteen novel transcripts of the human SCAF1 gene (SCAF1 v.2 - v.16), using nested PCR and NGS technology. In detail, extensive bioinformatic analysis revealed that these novel SCAF1 splice variants comprise a total of nine novel alternative splicing events between the annotated exons of the gene, thus producing seven novel SCAF1 transcripts with open-reading frames, which are predicted to encode novel SCAF1 isoforms and eight novel SCAF1 transcripts with premature termination codons that are likely long non-coding RNAs. Additionally, a novel 3' UTR was discovered and cloned using nested 3' RACE and was validated with Sanger sequencing. In order to validate the NGS findings as well as to investigate the expression profile of each novel transcript, RT-PCR experiments were carried out with the use of variant-specific primers. Since SCAF1 is implicated in many human malignancies, qualifying as a potential biomarker, the quantification of the presented novel transcripts in human samples may have clinical applications in different types of cancer.

Profiles of miRNA Isoforms and tRNA Fragments in Prostate Cancer

MicroRNA (miRNA) isoforms (“isomiRs”) and tRNA-derived fragments (“tRFs”) are powerful regulatory non-coding RNAs (ncRNAs). In human tissues, both types of molecules are abundant, with expression patterns that depend on a person’s race, sex and population origin. Here, we present our analyses of the Prostate Cancer (PRAD) datasets of The Cancer Genome Atlas (TCGA) from the standpoint of isomiRs and tRFs. This study represents the first simultaneous examination of isomiRs and tRFs in a large cohort of PRAD patients. We find that isomiRs and tRFs have extensive correlations with messenger RNAs (mRNAs). These correlations are disrupted in PRAD, which suggests disruptions of the regulatory network in the disease state. Notably, we find that the profiles of isomiRs and tRFs differ in patients belonging to different races. We hope that the presented findings can lay the groundwork for future research efforts aimed at elucidating the functional roles of the numerous and distinct members of these two categories of ncRNAs that are present in PRAD.

Estradiol shows anti-skin cancer activities through decreasing MDM2 expression

Estradiol plays important roles in many biological responses inducing tumor genesis and cancer treatment. However, the effects of estradiol on tumors were inconsistent among a lot of researches and the mechanism is not fully understood. Our previous study indicated that splicing factor hnRNPA1 could bind to the human homologue of mouse double minute (MDM2), an oncogene which has been observed to be over-expressed in numerous types of cancers. In this research, we investigated whether and how estradiol correlate to cancer cell behaviors through heterogeneous nuclear ribonucleoprotein (hnRNPA1) and MDM2. Results showed that 10×10^-13Mestradiol elevated the expression of hnRNPA1 regardless ER expression in cells, and then down-regulated the expression of MDM2. At the same time, estradiol inhibited cell proliferation, migration and epithelial-mesenchymal transition progression of A375 and GLL19 cells. While, knocking down hnRNPA1 through the transfection of hnRNPA1 siRNA led to the increase of MDM2 at both protein level and gene level In vivo experiment, subcutaneous injection with estradiol every two days near the tumor at doses of 2.5mg/kg/d suppressed tumor growth and reduced MDM2 expression. In a word, via increasing hnRNPA1 level and then reducing the expression of MDM2, estradiol prevented carcinogenesis in melanomas. We confirmed therapeutic effect of estradiol, as well as a new way for estradiol to resist skin cancer.

Effect of miR‑146a‑5p on tumor growth in NSCLC using chick chorioallantoic membrane assay and bioinformatics investigation

Our previous study demonstrated that the expression of miR-146a-5p was downregulated in non-small cell lung cancer (NSCLC) tissue, which affected the progression and prognosis of patients with NSCLC. Thus, the present study was conducted to investigate the functional mechanism of miR-146a-5p in tumorigenesis and angiogenesis in NSCLC. Following the construction of a H460 NSCLC cell line in which miR-146a-5p was overexpressed via lentivirus transduction, the NSCLC chick embryo chorioallantoic membrane (CAM) model was established by transplanting miR-146a-5p-overexpressing NSCLC cells into the CAM. Then, the size of the neoplasms within the CAM was measured, the vessel ratio was calculated, and the cellular morphology, metastasis and inflammation of tumor cell was observed using hematoxylin and eosin staining. The target genes of miR-146a-5p were predicted by 12 online software programs; these genes were then subjected to Gene Ontology enrichment analysis and Kyoto Encyclopedia of Genes and Genomes pathway annotations using the Database for Annotation, Visualization and Integrated Discovery 6.7 as well as constructed into a protein interaction network using protein-protein interaction from Search Tool for the Retrieval of Interacting Genes/Proteins. The xenograft tumor size and angiogenesis conditions of the miR-146a-5p-overexpressing group (volume 6.340±0.066 mm³, vessel ratio 9.326±0.083) was obviously restricted (P<0.001) when compared with the low expression group (volume 30.13±0.06 mm³, vessel ratio 16.94±0.11). In addition, marked necrosis along with inflammatory cell infiltration was observed with the HE-stained slices from the miR-146a-5p low expression group. Regarding the results of the target gene prediction, cancer and toll-like receptor signaling were the two most significant pathways represented among the target genes, while JUN, EGFR and RAC1 were the most relevant proteins among the selected potential targets of miR-146a-5p. In a CAM xenograft tumor model, overexpression of miR-146a-5p inhibited the tumorigenesis and angiogenesis of an NSCLC cell line. miR-146a-5p may act as a tumor suppressor gene in NSCLC and have moderate prognostic value in lung cancer.

RNA splicing and splicing regulator changes in prostate cancer pathology

Changes in mRNA splice patterns have been associated with key pathological mechanisms in prostate cancer progression. The androgen receptor (abbreviated AR) transcription factor is a major driver of prostate cancer pathology and activated by androgen steroid hormones. Selection of alternative promoters by the activated AR can critically alter gene function by switching mRNA isoform production, including creating a pro-oncogenic isoform of the normally tumour suppressor gene TSC2. A number of androgen-regulated genes generate alternatively spliced mRNA isoforms, including a prostate-specific splice isoform of ST6GALNAC1 mRNA. ST6GALNAC1 encodes a sialyltransferase that catalyses the synthesis of the cancer-associated sTn antigen important for cell mobility. Genetic rearrangements occurring early in prostate cancer development place ERG oncogene expression under the control of the androgen-regulated TMPRSS2 promoter to hijack cell behaviour. This TMPRSS2-ERG fusion gene shows different patterns of alternative splicing in invasive versus localised prostate cancer. Alternative AR mRNA isoforms play a key role in the generation of prostate cancer drug resistance, by providing a mechanism through which prostate cancer cells can grow in limited serum androgen concentrations. A number of splicing regulator proteins change expression patterns in prostate cancer and may help drive key stages of disease progression. Up-regulation of SRRM4 establishes neuronal splicing patterns in neuroendocrine prostate cancer. The splicing regulators Sam68 and Tra2β increase expression in prostate cancer. The SR protein kinase SRPK1 that modulates the activity of SR proteins is up-regulated in prostate cancer and has already given encouraging results as a potential therapeutic target in mouse models.

STARs in the CNS

STAR (signal transduction and activation of RNA) proteins regulate splicing of target genes that have roles in neural connectivity, survival and myelination in the vertebrate nervous system. These regulated splicing targets include mRNAs such as the Neurexins (Nrxn), SMN2 (survival of motor neuron) and MAG (myelin-associated glycoprotein). Recent work has made it possible to identify and validate STAR protein splicing targets in vivo by using genetically modified mouse models. In this review, we will discuss the importance of STAR protein splicing targets in the CNS (central nervous system).

The oncogenic role of the spliced somatostatin receptor sst5TMD4 variant in prostate cancer

sst5TMD4, a splice variant of the sst5 gene, is overexpressed and associated with aggressiveness in various endocrine-related tumors, but its presence, functional role, and mechanisms of actions in prostate cancer (PCa)-the most common cancer type in males-is completely unexplored. In this study, formalin-fixed, paraffin-embedded prostate pieces from patients with localized PCa, which included tumoral and nontumoral adjacent regions (n = 45), fresh biopsies from patients with high-risk PCa (n = 52), and healthy fresh prostates from cystoprostatectomies (n = 14) were examined. In addition, PCa cell lines and xenograft models were used to determine the presence and functional role of sst5TMD4. Results demonstrated that sst5TMD4 is overexpressed (mRNA/protein) in PCa samples, and this is especially drastic in metastatic and/or high Gleason score tumor samples. Remarkably, sst5TMD4 expression was associated with an altered frequency of 2 single-nucleotide polymorphisms: rs197055 and rs12599155. In addition, PCa cell lines and xenograft models were used to demonstrate that sst5TMD4 overexpression increases cell proliferation and migration in PCa cells and induces larger tumors in nude mice, whereas its silencing decreased proliferation and migration. Remarkably, sst5TMD4 overexpression activated multiple intracellular pathways (ERK/JNK, MYC/MAX, WNT, retinoblastoma), altered oncogenes and tumor suppressor gene expression, and disrupted the normal response to somatostatin analogs in PCa cells. Altogether, we demonstrate that sst5TMD4 is overexpressed in PCa, especially in those patients with a worse prognosis, and plays an important pathophysiologic role in PCa, which suggesting its potential as a biomarker and/or therapeutic target.-Hormaechea-Agulla, D., Jiménez-Vacas, J. M., Gómez-Gómez, E., L.-López, F., Carrasco-Valiente, J., Valero-Rosa, J., Moreno, M. M., Sánchez-Sánchez, R., Ortega-Salas, R., Gracia-Navarro, F., Culler, M. D., Ibáñez-Costa, A., Gahete, M. D., Requena, M. J., Castaño, J. P., Luque, R. M. The oncogenic role of the spliced somatostatin receptor sst5TMD4 variant in prostate cancer.

Alternative splicing promotes tumour aggressiveness and drug resistance in African American prostate cancer

Clinical challenges exist in reducing prostate cancer (PCa) disparities. The RNA splicing landscape of PCa across racial populations has not been fully explored as a potential molecular mechanism contributing to race-related tumour aggressiveness. Here, we identify novel genome-wide, race-specific RNA splicing events as critical drivers of PCa aggressiveness and therapeutic resistance in African American (AA) men. AA-enriched splice variants of PIK3CD, FGFR3, TSC2 and RASGRP2 contribute to greater oncogenic potential compared with corresponding European American (EA)-expressing variants. Ectopic overexpression of the newly cloned AA-enriched variant, PIK3CD-S, in EA PCa cell lines enhances AKT/mTOR signalling and increases proliferative and invasive capacity in vitro and confers resistance to selective PI3Kδ inhibitor, CAL-101 (idelalisib), in mouse xenograft models. High PIK3CD-S expression in PCa specimens associates with poor survival. These results highlight the potential of RNA splice variants to serve as novel biomarkers and molecular targets for developmental therapeutics in aggressive PCa.

The Landscape of Isoform Switches in Human Cancers

Alternative usage of transcript isoforms from the same gene has been hypothesized as an important feature in cancers. However, differential usage of gene transcripts between conditions (isoform switching) has not been comprehensively characterized in and across cancer types. To this end, we developed methods for identification and visualization of isoform switches with predicted functional consequences. Using these methods, we characterized isoform switching in RNA-seq data from >5,500 cancer patients covering 12 solid cancer types. Isoform switches with potential functional consequences were common, affecting approximately 19% of multiple transcript genes. Among these, isoform switches leading to loss of DNA sequence encoding protein domains were more frequent than expected, particularly in pancancer switches. We identified several isoform switches as powerful biomarkers: 31 switches were highly predictive of patient survival independent of cancer types. Our data constitute an important resource for cancer researchers, available through interactive web tools. Moreover, our methods, available as an R package, enable systematic analysis of isoform switches from other RNA-seq datasets.Implications: This study indicates that isoform switches with predicted functional consequences are common and important in dysfunctional cells, which in turn means that gene expression should be analyzed at the isoform level. Visual Overview: http://mcr.aacrjournals.org/content/molcanres/15/9/1206/F1.large.jpg.Mol Cancer Res; 15(9); 1206-20. ©2017 AACR.

Integrative gene set enrichment analysis utilizing isoform-specific expression

Gene set enrichment analysis (GSEA) aims at identifying essential pathways, or more generally, sets of biologically related genes that are involved in complex human diseases. In the past, many studies have shown that GSEA is a very useful bioinformatics tool that plays critical roles in the innovation of disease prevention and intervention strategies. Despite its tremendous success, it is striking that conclusions of GSEA drawn from isolated studies are often sparse, and different studies may lead to inconsistent and sometimes contradictory results. Further, in the wake of next generation sequencing technologies, it has been made possible to measure genome-wide isoform-specific expression levels, calling for innovations that can utilize the unprecedented resolution. Currently, enormous amounts of data have been created from various RNA-seq experiments. All these give rise to a pressing need for developing integrative methods that allow for explicit utilization of isoform-specific expression, to combine multiple enrichment studies, in order to enhance the power, reproducibility, and interpretability of the analysis. We develop and evaluate integrative GSEA methods, based on two-stage procedures, which, for the first time, allow statistically efficient use of isoform-specific expression from multiple RNA-seq experiments. Through simulation and real data analysis, we show that our methods can greatly improve the performance in identifying essential gene sets compared to existing methods that can only use gene-level expression.

Protease Expression Levels in Prostate Cancer Tissue Can Explain Prostate Cancer-Associated Seminal Biomarkers-An Explorative Concept Study

Previously, we described prostate cancer (PCa) detection (83% sensitivity; 67% specificity) in seminal plasma by CE-MS/MS. Moreover, advanced disease was distinguished from organ-confined tumors with 80% sensitivity and 82% specificity. The discovered biomarkers were naturally occurring fragments of larger seminal proteins, predominantly semenogelin 1 and 2, representing endpoints of the ejaculate liquefaction. Here we identified proteases putatively involved in PCa specific protein cleavage, and examined gene expression and tissue protein levels, jointly with cell localization in normal prostate (nP), benign prostate hyperplasia (BPH), seminal vesicles and PCa using qPCR, Western blotting and confocal laser scanning microscopy. We found differential gene expression of chymase (CMA1), matrix metalloproteinases (MMP3, MMP7), and upregulation of MMP14 and tissue inhibitors (TIMP1 and TIMP2) in BPH. In contrast tissue protein levels of MMP14 were downregulated in PCa. MMP3/TIMP1 and MMP7/TIMP1 ratios were decreased in BPH. In seminal vesicles, we found low-level expression of most proteases and, interestingly, we also detected TIMP1 and low levels of TIMP2. We conclude that MMP3 and MMP7 activity is different in PCa compared to BPH due to fine regulation by their inhibitor TIMP1. Our findings support the concept of seminal plasma biomarkers as non-invasive tool for PCa detection and risk stratification.

Identification of Gene Transcription Start Sites and Enhancers Responding to Pulmonary Carbon Nanotube Exposure in Vivo

Increased use of nanomaterials in industry, medicine, and consumer products has raised concerns over their toxicity. To ensure safe use of nanomaterials, understanding their biological effects at the molecular level is crucial. In particular, the regulatory mechanisms responsible for the cascade of genes activated by nanomaterial exposure are not well-characterized. To this end, we profiled the genome-wide usage of gene transcription start sites and linked active enhancer regions in lungs of C57BL/6 mice 24 h after intratracheal instillation of a single dose of the multiwalled carbon nanotube (MWCNT) Mitsui-7. Our results revealed a massive gene regulatory response, where expression of key inflammatory genes (e.g., Csf3, Il24, and Fgf23) was increased >100-fold 24 h after Mitsui-7 exposure. Many of the Mitsui-7-responsive transcription start sites were alternative transcription start sites for known genes, and the number of alternative transcription start sites used in a given gene was correlated with overall Mitsui-7 response. Strikingly, genes that were up-regulated after Mitsui-7 exposure only through their main annotated transcription start site were linked to inflammatory and defense responses, while genes up-regulated only through alternative transcription start sites were functionally heterogeneous and not inflammation-associated. Furthermore, we identified almost 12 000 active enhancers, many of which were Mitsui-7-responsive, and we identified similarly responding putative target genes. Overall, our study provides the location and activity of Mitsui-7-induced enhancers and transcription start sites, providing a useful resource for targeted experiments elucidating the biological effects of nanomaterials and the identification of biomarkers for early detection of MWCNT-induced inflammation.

The Involvement of Splicing Factor hnRNP A1 in UVB-induced Alternative Splicing of hdm2

Human homolog double minute 2 (hdm2), an oncoprotein, which binds to tumor suppressor p53 to facilitate its degradation, has been known to contribute to tumorigenesis. Its splicing variants are reported to be highly expressed in many cancers and can be induced by ultraviolet B light (UVB). However, the mechanisms of how UVB radiation induces hdm2 alternative splicing still remain unclear. In this study, we investigated the roles of two common splicing factors, heterogeneous nuclear ribonucleoproteins (hnRNP) A1 and serine/arginine-rich splicing factor 1 (SRSF1), in regulating UVB-induced hdm2 splicing. Our study indicated that while the expression of both hnRNP A1 and SRSF1 are induced, only hnRNP A1 is involved in hdm2 alternative splicing upon UVB irradiation. Overexpression of hnRNP A1 resulted in decrease of full-length hdm2 (hdm2-FL) and increase of hdm2B, one of hdm2 alternate-splicing forms; while down-regulated hnRNP A1 expression led to the decrease of the hdm2-FL and hdm2B in HaCaT cells. Protein-mRNA binding assay confirmed that UVB irradiation could increase the binding of hnRNP A1 to hdm2 pre-mRNA. In conclusion, we elucidated that UVB induces alternative splicing of hdm2 by increasing the expression and the binding of hnRNP A1 to hdm2 full-length mRNA.

The RNA-binding protein Sam68 regulates expression and transcription function of the androgen receptor splice variant AR-V7

Castration-resistant (CR) prostate cancer (PCa) partly arises due to persistence of androgen receptor (AR) transcriptional activity in the absence of cognate ligand. An emerging mechanism underlying the CRPCa phenotype and predicting response to therapy is the expression of the constitutively-active AR-V7 splice variant generated by AR cryptic exon 3b inclusion. Here, we explore the role of the RNA-binding protein (RBP) Sam68 (encoded by KHDRBS1), which is over-expressed in clinical PCa, on AR-V7 expression and transcription function. Using a minigene reporter, we show that Sam68 controls expression of exon 3b resulting in an increase in endogenous AR-V7 mRNA and protein expression in RNA-binding-dependent manner. We identify a novel protein-protein interaction between Sam68 and AR-V7 mediated by a common domain shared with full-length AR, and observe these proteins in the cell nucleoplasm. Using a luciferase reporter, we demonstrate that Sam68 co-activates ligand-independent AR-V7 transcriptional activity in an RNA-binding-independent manner, and controls expression of the endogenous AR-V7-specific gene target UBE2C. Our data suggest that Sam68 has separable effects on the regulation of AR-V7 expression and transcriptional activity, through its RNA-binding capacity. Sam68 and other RBPs may control expression of AR-V7 and other splice variants as well as their downstream functions in CRPCa.

Unravelling personalized dysfunctional gene network of complex diseases based on differential network model

In the conventional analysis of complex diseases, the control and case samples are assumed to be of great purity. However, due to the heterogeneity of disease samples, many disease genes are even not always consistently up-/down-regulated, leading to be under-estimated. This problem will seriously influence effective personalized diagnosis or treatment. The expression variance and expression covariance can address such a problem in a network manner. But, these analyses always require multiple samples rather than one sample, which is generally not available in clinical practice for each individual. To extract the common and specific network characteristics for individual patients in this paper, a novel differential network model, e.g. personalized dysfunctional gene network, is proposed to integrate those genes with different features, such as genes with the differential gene expression (DEG), genes with the differential expression variance (DEVG) and gene-pairs with the differential expression covariance (DECG) simultaneously, to construct personalized dysfunctional networks. This model uses a new statistic-like measurement on differential information, i.e., a differential score (DEVC), to reconstruct the differential expression network between groups of normal and diseased samples; and further quantitatively evaluate different feature genes in the patient-specific network for each individual. This DEVC-based differential expression network (DEVC-net) has been applied to the study of complex diseases for prostate cancer and diabetes. (1) Characterizing the global expression change between normal and diseased samples, the differential gene networks of those diseases were found to have a new bi-coloured topological structure, where their non hub-centred sub-networks are mainly composed of genes/proteins controlling various biological processes. (2) The differential expression variance/covariance rather than differential expression is new informative sources, and can be used to identify genes or gene-pairs with discriminative power, which are ignored by traditional methods. (3) More importantly, DEVC-net is effective to measure the expression state or activity of different feature genes and their network or modules in one sample for an individual. All of these results support that DEVC-net indeed has a clear advantage to effectively extract discriminatively interpretable features of gene/protein network of one sample (i.e. personalized dysfunctional network) even when disease samples are heterogeneous, and thus can provide new features like gene-pairs, in addition to the conventional individual genes, to the analysis of the personalized diagnosis and prognosis, and a better understanding on the underlying biological mechanisms.

IIIDB: a database for isoform-isoform interactions and isoform network modules

BACKGROUND

Protein-protein interactions (PPIs) are key to understanding diverse cellular processes and disease mechanisms. However, current PPI databases only provide low-resolution knowledge of PPIs, in the sense that "proteins" of currently known PPIs generally refer to "genes." It is known that alternative splicing often impacts PPI by either directly affecting protein interacting domains, or by indirectly impacting other domains, which, in turn, impacts the PPI binding. Thus, proteins translated from different isoforms of the same gene can have different interaction partners.

RESULTS

Due to the limitations of current experimental capacities, little data is available for PPIs at the resolution of isoforms, although such high-resolution data is crucial to map pathways and to understand protein functions. In fact, alternative splicing can often change the internal structure of a pathway by rearranging specific PPIs. To fill the gap, we systematically predicted genome-wide isoform-isoform interactions (IIIs) using RNA-seq datasets, domain-domain interaction and PPIs. Furthermore, we constructed an III database (IIIDB) that is a resource for studying PPIs at isoform resolution. To discover functional modules in the III network, we performed III network clustering, and then obtained 1025 isoform modules. To evaluate the module functionality, we performed the GO/pathway enrichment analysis for each isoform module.

CONCLUSIONS

The IIIDB provides predictions of human protein-protein interactions at the high resolution of transcript isoforms that can facilitate detailed understanding of protein functions and biological pathways. The web interface allows users to search for IIIs or III network modules. The IIIDB is freely available at http://syslab.nchu.edu.tw/IIIDB.

Identification of a candidate prognostic gene signature by transcriptome analysis of matched pre- and post-treatment prostatic biopsies from patients with advanced prostate cancer

BACKGROUND

Although chemotherapy for prostate cancer (PCa) can improve patient survival, some tumours are chemo-resistant. Tumour molecular profiles may help identify the mechanisms of drug action and identify potential prognostic biomarkers. We performed in vivo transcriptome profiling of pre- and post-treatment prostatic biopsies from patients with advanced hormone-naive prostate cancer treated with docetaxel chemotherapy and androgen deprivation therapy (ADT) with an aim to identify the mechanisms of drug action and identify prognostic biomarkers.

METHODS

RNA sequencing (RNA-Seq) was performed on biopsies from four patients before and ~22 weeks after docetaxel and ADT initiation. Gene fusion products and differentially-regulated genes between treatment pairs were identified using TopHat and pathway enrichment analyses undertaken. Publically available datasets were interrogated to perform survival analyses on the gene signatures identified using cBioportal.

RESULTS

A number of genomic rearrangements were identified including the TMPRSS2/ERG fusion and 3 novel gene fusions involving the ETS family of transcription factors in patients, both pre and post chemotherapy. In total, gene expression analyses showed differential expression of at least 2 fold in 575 genes in post-chemotherapy biopsies. Of these, pathway analyses identified a panel of 7 genes (ADAM7, FAM72B, BUB1B, CCNB1, CCNB2, TTK, CDK1), including a cell cycle-related geneset, that were differentially-regulated following treatment with docetaxel and ADT. Using cBioportal to interrogate the MSKCC-Prostate Oncogenome Project dataset we observed a statistically-significant reduction in disease-free survival of patients with tumours exhibiting alterations in gene expression of the above panel of 7 genes (p = 0.015).

CONCLUSIONS

Here we report on the first "real-time" in vivo RNA-Seq-based transcriptome analysis of clinical PCa from pre- and post-treatment TRUSS-guided biopsies of patients treated with docetaxel chemotherapy plus ADT. We identify a chemotherapy-driven PCa transcriptome profile which includes the down-regulation of important positive regulators of cell cycle progression. A 7 gene signature biomarker panel has also been identified in high-risk prostate cancer patients to be of prognostic value. Future prospective study is warranted to evaluate the clinical value of this panel.

Whole transcriptome sequencing reveals extensive unspliced mRNA in metastatic castration-resistant prostate cancer

Men with metastatic prostate cancer who are treated with androgen deprivation therapies (ADT) usually relapse within 2 to 3 years with disease that is termed castration-resistant prostate cancer (CRPC). To identify the mechanism that drives these advanced tumors, paired-end RNA-sequencing (RNA-seq) was performed on a panel of CRPC bone marrow biopsy specimens. From this genome-wide approach, mutations were found in a series of genes with prostate cancer relevance, including AR, NCOR1, KDM3A, KDM4A, CHD1, SETD5, SETD7, INPP4B, RASGRP3, RASA1, TP53BP1, and CDH1, and a novel SND1:BRAF gene fusion. Among the most highly expressed transcripts were 10 noncoding RNAs (ncRNAs), including MALAT1 and PABPC1, which are involved in RNA processing. Notably, a high percentage of sequence reads mapped to introns, which were determined to be the result of incomplete splicing at canonical splice junctions. Using quantitative PCR (qPCR), a series of genes (AR, KLK2, KLK3, STEAP2, CPSF6, and CDK19) were confirmed to have a greater proportion of unspliced RNA in CRPC specimens than in normal prostate epithelium, untreated primary prostate cancer, and cultured prostate cancer cells. This inefficient coupling of transcription and mRNA splicing suggests an overall increase in transcription or defect in splicing.

IMPLICATIONS

Inefficient splicing in advanced prostate cancer provides a selective advantage through effects on microRNA networks but may render tumors vulnerable to agents that suppress rate-limiting steps in splicing.

Next-generation sequencing of advanced prostate cancer treated with androgen-deprivation therapy

BACKGROUND

Androgen-deprivation therapy (ADT) is standard treatment for locally advanced or metastatic prostate cancer (PCa). Many patients develop castration resistance (castration-resistant PCa [CRPC]) after approximately 2-3 yr, with a poor prognosis. The molecular mechanisms underlying CRPC progression are unclear.

OBJECTIVE

To undertake quantitative tumour transcriptome profiling prior to and following ADT to identify functionally important androgen-regulated pathways or genes that may be reactivated in CRPC.

DESIGN, SETTING, AND PARTICIPANTS

RNA sequencing (RNA-seq) was performed on tumour-rich, targeted prostatic biopsies from seven patients with locally advanced or metastatic PCa before and approximately 22 wk after ADT initiation. Differentially regulated genes were identified in treatment pairs and further investigated by quantitative reverse transcription-polymerase chain reaction (qRT-PCR) on cell lines and immunohistochemistry on a separate CRPC patient cohort. Functional assays were used to determine the effect of pathway modulation on cell phenotypes.

OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS

We searched for gene expression changes affecting key cell signalling pathways that may be targeted as proof of principle in a CRPC in vitro cell line model.

RESULTS AND LIMITATIONS

We identified ADT-regulated signalling pathways, including the Wnt/β-catenin signalling pathway, and observed overexpression of β-catenin in a subset of CRPC by immunohistochemistry. We validated 6 of 12 (50%) pathway members by qRT-PCR on LNCaP/LNCaP-AI cell RNAs, of which 4 (67%) demonstrated expression changes consistent with RNA-seq data. We show that the tankyrase inhibitor XAV939 (which promotes β-catenin degradation) reduced androgen-independent LNCaP-AI cell line growth compared with androgen-responsive LNCaP cells via an accumulation of cell proportions in the G0/G1 phase and reduction in the S and G2/M phases. Our biopsy protocol did not account for tumour heterogeneity, and pathway inhibition was limited to pharmacologic approaches.

CONCLUSIONS

RNA-seq of paired PCa samples revealed ADT-regulated signalling pathways. Proof-of-principle inhibition of the Wnt/β-catenin signalling pathway specifically delays androgen-independent PCa cell cycle progression and proliferation and warrants further investigation as a potential target for therapy for CRPC.

The RNA-binding protein hnRNPA2 regulates β-catenin protein expression and is overexpressed in prostate cancer

INTRODUCTION

The RNA-binding protein hnRNPA2 (HNRNPA2B1) is upregulated in cancer, where it controls alternative pre-mRNA splicing of cancer-relevant genes. Cytoplasmic hnRNPA2 is reported in aggressive cancers, but is functionally uncharacterized. We explored the role of hnRNPA2 in prostate cancer (PCa).

METHODS

hnRNPA2 function/localization/expression in PCa was determined using biochemical approaches (colony forming/proliferation/luciferase reporter assays/flow cytometry/immunohistocytochemistry). Binding of hnRNPA2 within cancer-relevant 3'-UTR mRNAs was identified by bioinformatics.

RESULTS

RNAi-mediated knockdown of hnRNPA2 reduced colony forming and proliferation, while hnRNPA2 overexpression increased proliferation of PCa cells. Nuclear hnRNPA2 is overexpressed in high-grade clinical PCa, and is also observed in the cytoplasm in some cases. Ectopic expression of a predominantly cytoplasmic variant hnRNPA2-ΔRGG also increased PCa cell proliferation, suggesting that cytoplasmic hnRNPA2 may also be functionally relevant in PCa. Consistent with its known cytoplasmic roles, hnRNPA2 was associated with 3'-UTR mRNAs of several cancer-relevant mRNAs including β-catenin (CTNNB1). Both wild-type hnRNPA2 and hnRNPA2-ΔRGG act on CTNNB1 3'-UTR mRNA, increasing endogenous CTNNB1 mRNA expression and β-catenin protein expression and nuclear localization.

CONCLUSION

Nuclear and cytoplasmic hnRNPA2 are present in PCa and appear to be functionally important. Cytoplasmic hnRNPA2 may affect the cancer cell phenotype through 3'-UTR mRNA-mediated regulation of β-catenin expression and other cancer-relevant genes.

A novel androgen-regulated isoform of the TSC2 tumour suppressor gene increases cell proliferation

TSC2 (Tuberous sclerosis complex 2) is an important tumour suppressor gene, mutations within which are linked to the development of tuberous sclerosis and implicated in multiple tumour types. TSC2 protein complexes with TSC1 and blocks the ability of the Rheb (Ras homolog enriched in brain) GTPase to activate mTOR (mammalian target of rapamycin), a crucial signal transducer which regulates protein synthesis and cell growth. Here, we report the characterisation of a novel isoform of TSC2 which is under direct control of the ligand-activated androgen receptor. TSC2 isoform A (TSC2A) is derived from an internal androgen-regulated alternative promoter and encodes a 508-amino acid cytoplasmic protein corresponding to the C-terminal region of full-length TSC2, lacking the interaction domain for TSC1 and containing an incomplete interaction domain required for Rheb inactivation. Expression of TSC2A is induced in response to androgens and full-length TSC2 is co-ordinately down-regulated, indicating an androgen-driven switch in TSC2 protein isoforms. In contrast to the well-characterised suppressive effect on cell proliferation of full-length TSC2 protein, both LNCaP and HEK293 cells over-expressing TSC2 isoform A proliferate more rapidly (measured by MTT assays) and have increased levels of cells in S-phase (measured by both Edu staining and FACS analysis). Our work indicates, for the first time, a novel role for this well-known tumour suppressor gene, which encodes an activator of cell proliferation in response to androgen stimulation.

A novel EGFR isoform confers increased invasiveness to cancer cells

As a validated therapeutic target in several human cancers, the EGF receptor (EGFR) provides a focus to gain deeper insights into cancer pathophysiology. In this study, we report the identification of a naturally occurring and widely expressed EGFR isoform termed EGFRvA, which substitutes a Ser/Thr-rich peptide for part of the carboxyl-terminal regulatory domain of the receptor. Intriguingly, EGFRvA expression relates more closely to histopathologic grade and poor prognosis in patients with glioma. Ectopic expression of EGFRvA in cancer cells conferred a higher invasive capacity than EGFR in vitro and in vivo. Mechanistically, EGFRvA stimulated expression of STAT3, which upregulated heparin-binding EGF (HB-EGF). Reciprocally, HB-EGF stimulated phosphorylation of EGFRvA at Y845 along with STAT3, generating a positive feedback loop that may reinforce invasive function. The significance of EGFRvA expression was reinforced by findings that it is attenuated by miR-542-5p, a microRNA that is a known tumor suppressor. Taken together, our findings define this newfound EGFR isoform as a key theranostic molecule.

Plasma proteomics, the Human Proteome Project, and cancer-associated alternative splice variant proteins

This article addresses three inter-related subjects: the development of the Human Plasma Proteome Peptide Atlas, the launch of the Human Proteome Project, and the emergence of alternative splice variant transcripts and proteins as important features of evolution and pathogenesis. The current Plasma Peptide Atlas provides evidence on which peptides have been detected for every protein confidently identified in plasma; there are links to their spectra and their estimated abundance, facilitating the planning of targeted proteomics for biomarker studies. The Human Proteome Project (HPP) combines a chromosome-centric C-HPP with a biology and disease-driven B/D-HPP, upon a foundation of mass spectrometry, antibody, and knowledgebase resource pillars. The HPP aims to identify the approximately 7000 "missing proteins" and to characterize all proteins and their many isoforms. Success will enable the larger research community to utilize newly-available peptides, spectra, informative MS transitions, and databases for targeted analyses of priority proteins for each organ and disease. Among the isoforms of proteins, splice variants have the special feature of greatly enlarging protein diversity without enlarging the genome; evidence is accumulating of striking differential expression of splice variants in cancers. In this era of RNA-sequencing and advanced mass spectrometry, it is no longer sufficient to speak simply of increased or decreased expression of genes or proteins without carefully examining the splice variants in the protein mixture produced from each multi-exon gene. This article is part of a Special Issue entitled: Biomarkers: A Proteomic Challenge.

Expression of Tra2 β in Cancer Cells as a Potential Contributory Factor to Neoplasia and Metastasis

The splicing regulator proteins SRSF1 (also known as ASF/SF2) and SRSF3 (also known as SRP20) belong to the SR family of proteins and can be upregulated in cancer. The SRSF1 gene itself is amplified in some cancer cells, and cancer-associated changes in the expression of MYC also increase SRSF1 gene expression. Increased concentrations of SRSF1 protein promote prooncogenic splicing patterns of a number of key regulators of cell growth. Here, we review the evidence that upregulation of the SR-related Tra2β protein might have a similar role in cancer cells. The TRA2B gene encoding Tra2β is amplified in particular tumours including those of the lung, ovary, cervix, stomach, head, and neck. Both TRA2B RNA and Tra2β protein levels are upregulated in breast, cervical, ovarian, and colon cancer, and Tra2β expression is associated with cancer cell survival. The TRA2B gene is a transcriptional target of the protooncogene ETS-1 which might cause higher levels of expression in some cancer cells which express this transcription factor. Known Tra2β splicing targets have important roles in cancer cells, where they affect metastasis, proliferation, and cell survival. Tra2β protein is also known to interact directly with the RBMY protein which is implicated in liver cancer.

New insights, recent advances, and current challenges in the biological treatment of multiple myeloma

INTRODUCTION

The availability of thalidomide, lenalidomide, and bortezomib has radically changed multiple myeloma (MM) treatment and significantly improved patients' outcome. Nevertheless, MM is still an incurable disease due to the development of resistance and relapse practically in all patients. Unraveling MM pathogenesis, identifying prognostically high-risk patient populations, and optimizing current treatment strategies are among the challenges we are facing to reach a cure for this disease.

AREAS COVERED

This article reviews recent advances of the genomic analysis of malignant plasma cells and summarizes new insights into the pathophysiologic role of the MM microenvironment and the clinical assessment of derived novel therapeutic strategies. Moreover, current efforts to improve risk stratification and drug development are discussed, and most recent results of Phase II and III clinical trials that aim to optimize existing treatment regimens and to assess the next-generation anti-MM strategies are discussed. A systematic search was conducted of the Pubmed Medline, Embase, and Cochrane Library databases for primary articles, as well as of conference abstracts (e.g., of the American Society of Hematology, the American Society of Clinical Oncology, the American Association of Cancer Research, the European Hematology Association, and the Multiple Myeloma Workshop 2013), practice guidelines, and registries of clinical trials.

EXPERT OPINION

Given continuing advances to overcome current treatment challenges in MM, we are confident that long-lasting responses can be expected in many of our patients within the next decade.

Isoform level expression profiles provide better cancer signatures than gene level expression profiles

BACKGROUND

The majority of mammalian genes generate multiple transcript variants and protein isoforms through alternative transcription and/or alternative splicing, and the dynamic changes at the transcript/isoform level between non-oncogenic and cancer cells remain largely unexplored. We hypothesized that isoform level expression profiles would be better than gene level expression profiles at discriminating between non-oncogenic and cancer cellsgene level.

METHODS

We analyzed 160 Affymetrix exon-array datasets, comprising cell lines of non-oncogenic or oncogenic tissue origins. We obtained the transcript-level and gene level expression estimates, and used unsupervised and supervised clustering algorithms to study the profile similarity between the samples at both gene and isoform levels.

RESULTS

Hierarchical clustering, based on isoform level expressions, effectively grouped the non-oncogenic and oncogenic cell lines with a virtually perfect homogeneity-grouping rate (97.5%), regardless of the tissue origin of the cell lines. However, gene levelthis rate was much lower, being 75% at best based on the gene level expressions. Statistical analyses of the difference between cancer and non-oncogenic samples identified the existence of numerous genes with differentially expressed isoforms, which otherwise were not significant at the gene level. We also found that canonical pathways of protein ubiquitination, purine metabolism, and breast-cancer regulation by stathmin1 were significantly enriched among genes thatshow differential expression at isoform level but not at gene level.

CONCLUSIONS

In summary, cancer cell lines, regardless of their tissue of origin, can be effectively discriminated from non-cancer cell lines at isoform level, but not at gene level. This study suggests the existence of an isoform signature, rather than a gene signature, which could be used to distinguish cancer cells from normal cells.

The proto-oncogene ERG is a target of microRNA miR-145 in prostate cancer

Prostate cancer is a leading cause of cancer mortality in men. One of the distinct characteristics of prostate cancer is over-expression of the ERG proto-oncogene. The TMPRSS2-ERG gene fusion, the most common gene fusion, is found in approximately 50% of prostate cancer cases. We show that certain microRNAs are extensively deregulated in prostate cancer cell lines and primary clinical cancer samples. MicroRNAs are capable of modulating post-transcriptional gene expression via inhibition of protein synthesis. Independent target prediction methods have indicated that the 3' untranslated region of the ERG mRNA is a potential target of miR-145. miR-145 is consistently down-regulated in prostate cancer. Here we show that the ERG 3' untranslated region is a regulative target of miR-145 in vitro. Ectopic expression of miR-145 led to a reduction in expression of the ERG protein. We analyzed 26 prostate cancer samples and corresponding normal tissue. ERG protein expression was found to be elevated in the tumor samples, together with increased expression of several ERG isoforms. We identified ERG proteins of 35 and 24 kDa, which may represent unknown ERG splice variants. Analyses of miR-145 and ERG mRNA expression revealed a general down-regulation of miR-145 irrespective of the presence or absence of translocations involving ERG. This observation indicates that down-regulation of miR-145 may contribute to the increased expression of most ERG splice variants sharing the miR-145 target sequence in their 3' untranslated region.