Alternative splicing and biological heterogeneity in prostate cancer

Personal genomes, quantitative dynamic omics and personalized medicine

The rapid technological developments following the Human Genome Project have made possible the availability of personalized genomes. As the focus now shifts from characterizing genomes to making personalized disease associations, in combination with the availability of other omics technologies, the next big push will be not only to obtain a personalized genome, but to quantitatively follow other omics. This will include transcriptomes, proteomes, metabolomes, antibodyomes, and new emerging technologies, enabling the profiling of thousands of molecular components in individuals. Furthermore, omics profiling performed longitudinally can probe the temporal patterns associated with both molecular changes and associated physiological health and disease states. Such data necessitates the development of computational methodology to not only handle and descriptively assess such data, but also construct quantitative biological models. Here we describe the availability of personal genomes and developing omics technologies that can be brought together for personalized implementations and how these novel integrated approaches may effectively provide a precise personalized medicine that focuses on not only characterization and treatment but ultimately the prevention of disease.

Promise of personalized omics to precision medicine

The rapid development of high-throughput technologies and computational frameworks enables the examination of biological systems in unprecedented detail. The ability to study biological phenomena at omics levels in turn is expected to lead to significant advances in personalized and precision medicine. Patients can be treated according to their own molecular characteristics. Individual omes as well as the integrated profiles of multiple omes, such as the genome, the epigenome, the transcriptome, the proteome, the metabolome, the antibodyome, and other omics information are expected to be valuable for health monitoring, preventative measures, and precision medicine. Moreover, omics technologies have the potential to transform medicine from traditional symptom-oriented diagnosis and treatment of diseases toward disease prevention and early diagnostics. We discuss here the advances and challenges in systems biology-powered personalized medicine at its current stage, as well as a prospective view of future personalized health care at the end of this review.

Beyond comparing means: the usefulness of analyzing interindividual variation in gene expression for identifying genes associated with cancer development

Identifying genes associated with cancer development is typically accomplished by comparing mean expression values in normal and tumor tissues, which identifies differentially expressed (DE) genes. Interindividual variation (IV) in gene expression is indirectly included in DE gene identification because given the same absolute differences in means, genes with lower variance tend to have lower p-values. We explored the direct use of IV in gene expression to identify candidate genes associated with cancer development. We focused on prostate (PCa) and lung (LC) cancers and compared IV in the expression level of genes shown to be cancer related with that in all other genes in the human genome. Compared with all those other genes, cancer-related genes tended to have greater IV in normal tissues and a greater increase in IV during the transition from normal to tumorous tissue. Genes without significantly different mean expression values between tumor and normal tissues but with greater IV in tumor than in normal tissue (note: the DE-based approach completely ignores those genes) had stronger associations with clinically important features like Gleason score in PCa or tumor histology in LC than all other genes were. Our results suggest that analyzing IV in gene expression level is useful in identifying novel candidate genes associated with cancer development.

Androgen receptor (AR) aberrations in castration-resistant prostate cancer

Genetic aberrations affecting the androgen receptor (AR) are rare in untreated prostate cancers (PCs) but have been found in castration-resistant prostate cancers (CRPCs). Further, successful treatment with novel endocrine therapies indicates that CRPCs remain androgen-sensitive. Known AR aberrations include amplification of the AR gene leading to the overexpression of the receptor, point mutations of AR resulting in promiscuous ligand usage, and constitutively active AR splice variants. Gain, or amplification, of the AR gene is one of the most frequent genetic alterations observed in CRPCs. Up to 80% of CRPCs have been reported to carry an elevated AR gene copy number, and about 30% have a high-level amplification of the gene. AR mutations are also commonly observed and have been found in approximately 10-30% of the CRPC treated with antiandrogens; however, the frequency and significance of AR splice variants is still unclear. Because AR aberrations are found almost exclusively in CRPC, these alterations must have been selected for during therapy. Interestingly, these aberrations lead to activation of the receptor, despite treatment-induced emergence of therapy-resistant tumor clones. Therefore, future novel treatment strategies should focus on suppressing AR activity in CRPC.

Transcriptome-wide detection of differentially expressed coding and non-coding transcripts and their clinical significance in prostate cancer

Prostate cancer is a clinically and biologically heterogeneous disease. Deregulation of splice variants has been shown to contribute significantly to this complexity. High-throughput technologies such as oligonucleotide microarrays allow for the detection of transcripts that play a role in disease progression in a transcriptome-wide level. In this study, we use a publicly available dataset of normal adjacent, primary tumor, and metastatic prostate cancer samples (GSE21034) to detect differentially expressed coding and non-coding transcripts between these disease states. To achieve this, we focus on transcript-specific probe selection regions, that is, those probe sets that correspond unambiguously to a single transcript. Based on this, we are able to pinpoint at the transcript-specific level transcripts that are differentially expressed throughout prostate cancer progression. We confirm previously reported cases and find novel transcripts for which no prior implication in prostate cancer progression has been made. Furthermore, we show that transcript-specific differential expression has unique prognostic potential and provides a clinically significant source of biomarker signatures for prostate cancer risk stratification. The results presented here serve as a catalog of differentially expressed transcript-specific markers throughout prostate cancer progression that can be used as basis for further development and translation into the clinic.

Alternative transcription and alternative splicing in cancer

In recent years, the notion of "one gene makes one protein that functions in one signaling pathway" in mammalian cells has been shown to be overly simplistic. Recent genome-wide studies suggest that at least half of the human genes, including many therapeutic target genes, produce multiple protein isoforms through alternative splicing and alternative usage of transcription initiation and/or termination. For example, alternative splicing of the vascular endothelial growth factor gene (VEGFA) produces multiple protein isoforms, which display either pro-angiogenic or anti-angiogenic activities. Similarly, for the majority of human genes, the inclusion or exclusion of exonic sequences enhances the generation of transcript variants and/or protein isoforms that can vary in structure and functional properties. Many of the isoforms produced in this manner are tightly regulated during normal development but are misregulated in cancer cells. Altered expression of transcript variants and protein isoforms for numerous genes is linked with disease and its prognosis, and cancer cells manipulate regulatory mechanisms to express specific isoforms that confer drug resistance and survival advantages. Emerging insights indicate that modulating the expression of transcript and protein isoforms of a gene may hold the key to impeding tumor growth and act as a model for efficient targeting of disease-associated genes at the isoform level. This review highlights the role and regulation of alternative transcription and splicing mechanisms in generating the transcriptome, and the misuse and diagnostic/prognostic potential of alternative transcription and splicing in cancer.

RNA-seq analysis of prostate cancer in the Chinese population identifies recurrent gene fusions, cancer-associated long noncoding RNAs and aberrant alternative splicings

There are remarkable disparities among patients of different races with prostate cancer; however, the mechanism underlying this difference remains unclear. Here, we present a comprehensive landscape of the transcriptome profiles of 14 primary prostate cancers and their paired normal counterparts from the Chinese population using RNA-seq, revealing tremendous diversity across prostate cancer transcriptomes with respect to gene fusions, long noncoding RNAs (long ncRNA), alternative splicing and somatic mutations. Three of the 14 tumors (21.4%) harbored a TMPRSS2-ERG fusion, and the low prevalence of this fusion in Chinese patients was further confirmed in an additional tumor set (10/54=18.5%). Notably, two novel gene fusions, CTAGE5-KHDRBS3 (20/54=37%) and USP9Y-TTTY15 (19/54=35.2%), occurred frequently in our patient cohort. Further systematic transcriptional profiling identified numerous long ncRNAs that were differentially expressed in the tumors. An analysis of the correlation between expression of long ncRNA and genes suggested that long ncRNAs may have functions beyond transcriptional regulation. This study yielded new insights into the pathogenesis of prostate cancer in the Chinese population.

Current and emerging approaches to define intestinal epithelium-specific transcriptional networks

Upon developmental or environmental cues, the composition of transcription factors in a transcriptional regulatory network is deeply implicated in controlling the signature of the gene expression and thereby specifies the cell or tissue type. Novel methods including ChIP-chip and ChIP-Seq have been applied to analyze known transcription factors and their interacting regulatory DNA elements in the intestine. The intestine is an example of a dynamic tissue where stem cells in the crypt proliferate and undergo a differentiation process toward the villus. During this differentiation process, specific regulatory networks of transcription factors are activated to target specific genes, which determine the intestinal cell fate. The expanding genomewide mapping of transcription factor binding sites and construction of transcriptional regulatory networks provide new insight into how intestinal differentiation occurs. This review summarizes the current overview of the transcriptional regulatory networks driving epithelial differentiation in adult intestine. The novel technologies that have been implied to study these networks are presented and their prospects for implications in future research are also addressed.

Engineering of targeted nanoparticles for cancer therapy using internalizing aptamers isolated by cell-uptake selection

One of the major challenges in the development of targeted nanoparticles (NPs) for cancer therapy is to discover targeting ligands that allow for differential binding and uptake by the target cancer cells. Using prostate cancer (PCa) as a model disease, we developed a cell-uptake selection strategy to isolate PCa-specific internalizing 2'-O-methyl RNA aptamers (Apts) for NP incorporation. Twelve cycles of selection and counter-selection were done to obtain a panel of internalizing Apts, which can distinguish PCa cells from nonprostate and normal prostate cells. After Apt characterization, size minimization, and conjugation of the Apts with fluorescently labeled polymeric NPs, the NP-Apt conjugates exhibit PCa specificity and enhancement in cellular uptake when compared to nontargeted NPs lacking the internalizing Apts. Furthermore, when docetaxel, a chemotherapeutic agent used for the treatment of PCa, was encapsulated within the NP-Apt, a significant improvement in cytotoxicity was achieved in targeted PCa cells. Rather than isolating high-affinity Apts as reported in previous selection processes, our selection strategy was designed to enrich cancer cell-specific internalizing Apts. A similar cell-uptake selection strategy may be used to develop specific internalizing ligands for a myriad of other diseases and can potentially facilitate delivering various molecules, including drugs and siRNAs, into target cells.

Identification of novel androgen-regulated pathways and mRNA isoforms through genome-wide exon-specific profiling of the LNCaP transcriptome

Androgens drive the onset and progression of prostate cancer (PCa) by modulating androgen receptor (AR) transcriptional activity. Although several microarray-based studies have identified androgen-regulated genes, here we identify in-parallel global androgen-dependent changes in both gene and alternative mRNA isoform expression by exon-level analyses of the LNCaP transcriptome. While genome-wide gene expression changes correlated well with previously-published studies, we additionally uncovered a subset of 226 novel androgen-regulated genes. Gene expression pathway analysis of this subset revealed gene clusters associated with, and including the tyrosine kinase LYN, as well as components of the mTOR (mammalian target of rapamycin) pathway, which is commonly dysregulated in cancer. We also identified 1279 putative androgen-regulated alternative events, of which 325 (∼25%) mapped to known alternative splicing events or alternative first/last exons. We selected 30 androgen-dependent alternative events for RT-PCR validation, including mRNAs derived from genes encoding tumour suppressors and cell cycle regulators. Of seven positively-validating events (∼23%), five events involved transcripts derived from alternative promoters of known AR gene targets. In particular, we found a novel androgen-dependent mRNA isoform derived from an alternative internal promoter within the TSC2 tumour suppressor gene, which is predicted to encode a protein lacking an interaction domain required for mTOR inhibition. We confirmed that expression of this alternative TSC2 mRNA isoform was directly regulated by androgens, and chromatin immunoprecipitation indicated recruitment of AR to the alternative promoter region at early timepoints following androgen stimulation, which correlated with expression of alternative transcripts. Together, our data suggest that alternative mRNA isoform expression might mediate the cellular response to androgens, and may have roles in clinical PCa.

Discovery of circulating microRNAs associated with human prostate cancer using a mouse model of disease

Circulating microRNAs (miRNAs) are emerging as useful non-invasive markers of disease. The objective of this study was to use a mouse model of prostate cancer as a tool to discover serum miRNAs that could be assessed in a clinical setting. Global miRNA profiling identified 46 miRNAs at significantly altered levels (p ≤ 0.05) in the serum of TRansgenic Adenocarcinoma of Mouse Prostate (TRAMP) mice with advanced prostate cancer compared to healthy controls. A subset of these miRNAs with known human homologues were validated in an independent cohort of mice and then measured in serum from men with metastatic castration-resistant prostate cancer (mCRPC; n = 25) or healthy men (n = 25). Four miRNAs altered in mice, mmu-miR-141, mmu-miR-298, mmu-miR-346 and mmu-miR-375, were also found to be at differential levels in the serum of men with mCRPC. Three of these (hsa-miR-141, hsa-miR-298 and hsa-miR-375) were upregulated in prostate tumors compared with normal prostate tissue, suggesting that they are released into the blood as disease progresses. Moreover, the intra-tumoral expression of hsa-miR-141 and hsa-miR-375 were predictors of biochemical relapse after surgery. This study is the first to demonstrate that specific serum miRNAs are common between human prostate cancer and a mouse model of the disease, highlighting the potential of such models for the discovery of novel biomarkers.

Splicing programs and cancer

Numerous studies report splicing alterations in a multitude of cancers by using gene-by-gene analysis. However, understanding of the role of alternative splicing in cancer is now reaching a new level, thanks to the use of novel technologies allowing the analysis of splicing at a large-scale level. Genome-wide analyses of alternative splicing indicate that splicing alterations can affect the products of gene networks involved in key cellular programs. In addition, many splicing variants identified as being misregulated in cancer are expressed in normal tissues. These observations suggest that splicing programs contribute to specific cellular programs that are altered during cancer initiation and progression. Supporting this model, recent studies have identified splicing factors controlling cancer-associated splicing programs. The characterization of splicing programs and their regulation by splicing factors will allow a better understanding of the genetic mechanisms involved in cancer initiation and progression and the development of new therapeutic targets.

Alternative transcript initiation and splicing as a response to DNA damage

Humans are exposed to the DNA damaging agent, ionizing radiation (IR), from background radiation, medical treatments, occupational and accidental exposures. IR causes changes in transcription, but little is known about alternative transcription in response to IR on a genome-wide basis. These investigations examine the response to IR at the exon level in human cells, using exon arrays to comprehensively characterize radiation-induced transcriptional expression products. Previously uncharacterized alternative transcripts that preferentially occur following IR exposure have been discovered. A large number of genes showed alternative transcription initiation as a response to IR. Dose-response and time course kinetics have also been characterized. Interestingly, most genes showing alternative transcript induction maintained these isoforms over the dose range and times tested. Finally, clusters of co-ordinately up- and down-regulated radiation response genes were identified at specific chromosomal loci. These data provide the first genome-wide view of the transcriptional response to ionizing radiation at the exon level. This study provides novel insights into alternative transcripts as a mechanism for response to DNA damage and cell stress responses in general.

A-to-I RNA editing: the "ADAR" side of human cancer

Carcinogenesis is a complex, multi-stage process depending on both endogenous and exogenous factors. In the past years, DNA mutations provided important clues to the comprehension of the molecular pathways involved in numerous cancers. Recently, post-transcriptional modification events, such as RNA editing, are emerging as new players in several human diseases, including tumours. A-to-I RNA editing changes the nucleotide sequence of target RNAs, introducing A-to-I/G "mutations". Since ADAR enzymes catalyse this nucleotide conversion, their expression/activity is essential and finely regulated in normal cells. This review summarizes the available knowledge on A-to-I RNA editing in the cancer field, giving a new view on how ADARs may play a role in carcinogenesis.

Sensitive gene fusion detection using ambiguously mapping RNA-Seq read pairs

MOTIVATION

Paired-end whole transcriptome sequencing provides evidence for fusion transcripts. However, due to the repetitiveness of the transcriptome, many reads have multiple high-quality mappings. Previous methods to find gene fusions either ignored these reads or required additional longer single reads. This can obscure up to 30% of fusions and unnecessarily discards much of the data.

RESULTS

We present a method for using paired-end reads to find fusion transcripts without requiring unique mappings or additional single read sequencing. Using simulated data and data from tumors and cell lines, we show that our method can find fusions with ambiguously mapping read pairs without generating numerous spurious fusions from the many mapping locations.

AVAILABILITY

A C++ and Python implementation of the method demonstrated in this article is available at http://exon.ucsd.edu/ShortFuse.

CONTACT

mckinsel@ucsd.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Multiple sclerosis etiology: beyond genes and environment

Multiple sclerosis (MS) is a disorder of the CNS with inflammatory and neurodegenerative components. The etiology is unknown, but there is evidence for a role of both genetic and environmental factors. Among the heritable factors, MHC class II genes are strongly involved, as well as genes coding for others molecules of immunological relevance, genes controlling neurobiological pathways and genes of unknown function. Among nonheritable factors, many infectious agents (mainly viruses) and environmental factors (e.g., smoke, sun exposition and diet) seem to be of etiologic importance. Here, we report and discuss recent findings in MS on largely unexplored fields: the alternative splicing of mRNAs and regulatory noncoding RNAs, the major sources of transcriptome diversity; and epigenetic changes with special attention paid to DNA methylation and histone acetylation, the main regulators of gene expression.

Alterations in LMTK2, MSMB and HNF1B gene expression are associated with the development of prostate cancer

Background

Genome wide association studies (GWAS) have identified several genetic variants that are associated with prostate cancer. Most of these variants, like other GWAS association signals, are located in non-coding regions of potential candidate genes, and thus could act at the level of the mRNA transcript.

Methods

We measured the expression and isoform usage of seven prostate cancer candidate genes in benign and malignant prostate by real-time PCR, and correlated these factors with cancer status and genotype at the GWAS risk variants.

Results

We determined that levels of LMTK2 transcripts in prostate adenocarcinomas were only 32% of those in benign tissues (p = 3.2 × 10^-7), and that an independent effect of genotype at variant rs6465657 on LMTK2 expression in benign (n = 39) and malignant tissues (n = 21) was also evident (P = 0.002). We also identified that whilst HNF1B(C) and MSMB2 comprised the predominant isoforms in benign tissues (90% and 98% of total HNF1B or MSMB expression), HNF1B(B) and MSMB1 were predominant in malignant tissue (95% and 96% of total HNF1B or MSMB expression; P = 1.7 × 10^-7 and 4 × 10^-4 respectively), indicating major shifts in isoform usage.

Conclusions

Our results indicate that the amount or nature of mRNA transcripts expressed from the LMTK2, HNF1B and MSMB candidate genes is altered in prostate cancer, and provides further evidence for a role for these genes in this disorder. The alterations in isoform usage we detect highlights the potential importance of alternative mRNA processing and moderation of mRNA stability as potentially important disease mechanisms.

Spatiotemporal controlled delivery of nanoparticles to injured vasculature

There are a number of challenges associated with designing nanoparticles for medical applications. We define two challenges here: (i) conventional targeting against up-regulated cell surface antigens is limited by heterogeneity in expression, and (ii) previous studies suggest that the optimal size of nanoparticles designed for systemic delivery is approximately 50-150 nm, yet this size range confers a high surface area-to-volume ratio, which results in fast diffusive drug release. Here, we achieve spatial control by biopanning a phage library to discover materials that target abundant vascular antigens exposed in disease. Next, we achieve temporal control by designing 60-nm hybrid nanoparticles with a lipid shell interface surrounding a polymer core, which is loaded with slow-eluting conjugates of paclitaxel for controlled ester hydrolysis and drug release over approximately 12 days. The nanoparticles inhibited human aortic smooth muscle cell proliferation in vitro and showed greater in vivo vascular retention during percutaneous angioplasty over nontargeted controls. This nanoparticle technology may potentially be used toward the treatment of injured vasculature, a clinical problem of primary importance.

Transcription variants of the prostate-specific PrLZ gene and their interaction with 14-3-3 proteins

We have reported isolation and characterization of the prostate-specific and androgen-regulated PrLZ gene abnormally expressed in prostate cancer. PrLZ is a potential biomarker for prostate cancer and a candidate oncogene promoting cell proliferation and survival in prostate cancer cells. A full delineation of the PrLZ gene and its gene products may provide clues to the mechanisms regulating its expression and function. In this report, we identified three additional exons in the PrLZ gene and recognized five transcript variants from alternative splicing that could be detected by RT-PCR and Western blotting. Structural comparison demonstrated that the PrLZ proteins are highly conserved among species. PrLZ contains multiple potential sites for interaction with other proteins. We used mammalian two-hybrid assays to demonstrate that PrLZ isoforms interact with 14-3-3 proteins, and multiple sites in the PrLZ may be involved in the interaction. Alternative splicing may contribute to abnormally enhanced PrLZ levels in prostate cancer, and interaction with 14-3-3 proteins may be a mechanism by which PrLZ promotes cell proliferation and survival during prostate cancer development and progression. This information is a valuable addition to the investigation of the oncogenic properties of the PrLZ gene.