moPepGen: Rapid and Comprehensive Proteoform Identification.

Gene expression is a multi-step transformation of biological information from its storage form (DNA) into functional forms (protein and some RNAs). Regulatory activities at each step of this transformation multiply a single gene into a myriad of proteoforms. Proteogenomics is the study of how genomic and transcriptomic variation creates this proteoform diversity, and is limited by the challenges of modeling the complexities of gene-expression. We therefore created moPepGen, a graph-based algorithm that comprehensively enumerates proteoforms in linear time. moPepGen works with multiple technologies, in multiple species and on all types of genetic and transcriptomic data. In human cancer proteomes, it detects and quantifies previously unobserved noncanonical peptides arising from germline and somatic genomic variants, noncoding open reading frames, RNA fusions and RNA circularization. By enabling efficient identification and quantitation of previously hidden proteins in both existing and new proteomic data, moPepGen facilitates all proteogenomics applications. It is available at: https://github.com/uclahs-cds/package-moPepGen.

The genomic and evolutionary landscapes of anaplastic thyroid carcinoma.

Anaplastic thyroid carcinoma is arguably the most lethal human malignancy. It often co-occurs with differentiated thyroid cancers, yet the molecular origins of its aggressivity are unknown. We sequenced tumor DNA from 329 regions of thyroid cancer, including 213 from patients with primary anaplastic thyroid carcinomas. We also whole genome sequenced 9 patients using multi-region sequencing of both differentiated and anaplastic thyroid cancer components. Using these data, we demonstrate thatanaplastic thyroid carcinomas have a higher burden of mutations than other thyroid cancers, with distinct mutational signatures and molecular subtypes. Further, different cancer driver genes are mutated in anaplastic and differentiated thyroid carcinomas, even those arising in a single patient. Finally, we unambiguously demonstrate that anaplastic thyroid carcinomas share a genomic origin with co-occurring differentiated carcinomas and emerge from a common malignant field through acquisition of characteristic clonal driver mutations.

A polygenic two-hit hypothesis for prostate cancer.

Prostate cancer is one of the most heritable cancers. Hundreds of germline polymorphisms have been linked to prostate cancer diagnosis and prognosis. Polygenic risk scores can predict genetic risk of a prostate cancer diagnosis. Although these scores inform the probability of developing a tumor, it remains unknown how germline risk influences the tumor molecular evolution. We cultivated a cohort of 1250 localized European-descent patients with germline and somatic DNA profiling. Men of European descent with higher genetic risk were diagnosed earlier and had less genomic instability and fewer driver genes mutated. Higher genetic risk was associated with better outcome. These data imply a polygenic "two-hit" model where germline risk reduces the number of somatic alterations required for tumorigenesis. These findings support further clinical studies of polygenic risk scores as inexpensive and minimally invasive adjuncts to standard risk stratification. Further studies are required to interrogate generalizability to more ancestrally and clinically diverse populations.

The proteomic landscape of glioblastoma recurrence reveals novel and targetable immunoregulatory drivers.

Glioblastoma (GBM) is characterized by extensive cellular and genetic heterogeneity. Its initial presentation as primary disease (pGBM) has been subject to exhaustive molecular and cellular profiling. By contrast, our understanding of how GBM evolves to evade the selective pressure of therapy is starkly limited. The proteomic landscape of recurrent GBM (rGBM), which is refractory to most treatments used for pGBM, are poorly known. We, therefore, quantified the transcriptome and proteome of 134 patient-derived pGBM and rGBM samples, including 40 matched pGBM-rGBM pairs. GBM subtypes transition from pGBM to rGBM towards a preferentially mesenchymal state at recurrence, consistent with the increasingly invasive nature of rGBM. We identified immune regulatory/suppressive genes as important drivers of rGBM and in particular 2-5-oligoadenylate synthase 2 (OAS2) as an essential gene in recurrent disease. Our data identify a new class of therapeutic targets that emerge from the adaptive response of pGBM to therapy, emerging specifically in recurrent disease and may provide new therapeutic opportunities absent at pGBM diagnosis.

Estimation of tumor cell total mRNA expression in 15 cancer types predicts disease progression.

Single-cell RNA sequencing studies have suggested that total mRNA content correlates with tumor phenotypes. Technical and analytical challenges, however, have so far impeded at-scale pan-cancer examination of total mRNA content. Here we present a method to quantify tumor-specific total mRNA expression (TmS) from bulk sequencing data, taking into account tumor transcript proportion, purity and ploidy, which are estimated through transcriptomic/genomic deconvolution. We estimate and validate TmS in 6,590 patient tumors across 15 cancer types, identifying significant inter-tumor variability. Across cancers, high TmS is associated with increased risk of disease progression and death. TmS is influenced by cancer-specific patterns of gene alteration and intra-tumor genetic heterogeneity as well as by pan-cancer trends in metabolic dysregulation. Taken together, our results indicate that measuring cell-type-specific total mRNA expression in tumor cells predicts tumor phenotypes and clinical outcomes.

Methylation Subtypes of Primary Prostate Cancer Predict Poor Prognosis.

BACKGROUND: Patients with prostate cancer experience heterogeneous outcomes after radical prostatectomy. Genomic studies including The Cancer Genome Atlas (TCGA) have reported molecular signatures of prostate cancer, but few studies have assessed the prognostic effects of DNA methylation profiles. METHODS: We conducted the largest methylome subtyping analysis for primary prostate tumors to date, using methylome data from three patient populations: TCGA, a prostate cancer cohort study conducted at the Fred Hutchinson Cancer Research Center (FH; Seattle, WA), and the Canadian International Cancer Genome Consortium (ICGC) cohort. Four subtypes were detected in the TCGA dataset, then independently assigned to FH and ICGC cohort data. The identified methylation subtypes were assessed for association with cancer prognosis in the above three patient populations. RESULTS: Using a set of hypermethylated CpG sites, four methylation subtypes were identified in TCGA. Compared with subtype 1, subtype 4 had an HR of 2.09 (P = 0.029) for biochemical recurrence (BCR) in TCGA patients. HRs of 2.76 (P = 0.002) for recurrence and 9.73 (P = 0.002) for metastatic-lethal (metastasis or prostate cancer-specific death) outcomes were observed in the FH cohort. A similar pattern of association was noted in the Canadian ICGC cohort, though HRs were not statistically significant. CONCLUSIONS: A hypermethylated subtype was associated with an increased hazard of recurrence and mortality in three studies with prostate tumor methylome data. Further molecular work is needed to understand the effect of methylation subtypes on cancer prognosis. IMPACT: This study identified a DNA methylation subtype that was associated with worse prostate cancer prognosis after radical prostatectomy.

The telomere length landscape of prostate cancer.

Replicative immortality is a hallmark of cancer, and can be achieved through telomere lengthening and maintenance. Although the role of telomere length in cancer has been well studied, its association to genomic features is less well known. Here, we report the telomere lengths of 392 localized prostate cancer tumours and characterize their relationship to genomic, transcriptomic and proteomic features. Shorter tumour telomere lengths are associated with elevated genomic instability, including single-nucleotide variants, indels and structural variants. Genes involved in cell proliferation and signaling are correlated with tumour telomere length at all levels of the central dogma. Telomere length is also associated with multiple clinical features of a tumour. Longer telomere lengths in non-tumour samples are associated with a lower rate of biochemical relapse. In summary, we describe the multi-level integration of telomere length, genomics, transcriptomics and proteomics in localized prostate cancer.

Somatic driver mutation prevalence in 1844 prostate cancers identifies ZNRF3 loss as a predictor of metastatic relapse.

Driver gene mutations that are more prevalent in metastatic, castration-resistant prostate cancer (mCRPC) than localized disease represent candidate prognostic biomarkers. We analyze 1,844 localized (1,289) or mCRPC (555) tumors and quantify the prevalence of 113 somatic driver single nucleotide variants (SNVs), copy number aberrations (CNAs), and structural variants (SVs) in each state. One-third are significantly more prevalent in mCRPC than expected while a quarter are less prevalent. Mutations in AR and its enhancer are more prevalent in mCRPC, as are those in TP53, MYC, ZNRF3 and PRKDC. ZNRF3 loss is associated with decreased ZNRF3 mRNA abundance, WNT, cell cycle & PRC1/2 activity, and genomic instability. ZNRF3 loss, RNA downregulation and hypermethylation are prognostic of metastasis and overall survival, independent of clinical and pathologic indices. These data demonstrate a strategy for identifying biomarkers of localized cancer aggression, with ZNRF3 loss as a predictor of metastasis in prostate cancer.

Loss of CDCP1 triggers FAK activation in detached prostate cancer cells.

A major metastasis suppressing mechanism is the rapid apoptotic death of cancer cells upon detachment from extracellular matrix, a process called anoikis. Focal adhesion kinase (PTK2/FAK) is a key enzyme involved in evasion of anoikis. We show that loss of the Cub-domain containing protein-1 (CDCP1), paradoxically stimulates FAK activation in the detached state of prostate cancer cells. In CDCP1low DU145 and PC3 prostate cancer cells, detachment-activation of FAK occurs through local production of PI(4,5)P2. PI(4,5)P2 is generated by the PIP5K1c-201 splicing isoform of PIP5K1c, which contains a unique SRC phosphorylation site. In the detached state, reduced expression of CDCP1 and an alternative CDCP1-independent SRC activation mechanism triggers PIP5K1c-pY644 phosphorylation by SRC. This causes a switch of Talin binding from ß1-integrin to PIP5K1c-pY644 and leads to activation of PIP5K1c-FAK. Reduced CDCP1 expression also inactivates CDK5, a negative regulator of PIP5K1c. Furthermore, immersion of prostate cancer cells in 10% human plasma or fetal bovine serum is required for activation of PIP5K1c-FAK. The PIP5K1c induced detachment-activation of FAK in preclinical models sensitizes CDCP1low prostate cancer cells to FAK inhibitors. In patients, CDCP1High versus CDCP1low circulating tumor cells differ in expression of AR-v7, ONECUT2 and HOXB13 oncogenes and TMPRSS2 and display intra-patient heterogeneity of FAK-pY397 expression. Taken together, CDCP1low and CDCP1high detached prostate cancer cells activate distinct cytoplasmic kinase complexes and targetable transcription factors, which has important therapeutic implications.

Decreased ATM Protein Expression Is Substantiated with PTEN Loss in Defining Aggressive Phenotype of Prostate Cancer Associated with Lethal Disease.

BACKGROUND: Ataxia Telangiectasia Mutated (ATM) serine/threonine protein kinase is a known tumor suppressor, involved in DNA damage repair. It has prognostic and predictive therapeutic implications and is associated with aggressive prostate cancer (PCa). OBJECTIVE: To investigate the prognostic value of ATM protein expression in PCa patients and assessed the combined value of ATM, ERG, and PTEN status. DESIGN SETTING AND PARTICIPANTS: This study consisted of 303 patients with incidental, locally advanced, and castrate-resistant PCa by transurethral resection of the prostate (TURP). OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: TURP samples from 303 PCa patients were assessed by immunohistochemistry (IHC for ATM, ERG, and PTEN. Individual and combined marker status were correlated with International Society of Urological Pathology Gleason grade group, overall survival (OS), and PCa-specific mortality (PCSM). RESULTS AND LIMITATIONS: Decreased ATM expression (negative/weak intensity) occurred in 164/303 (54.1%) patients, and was associated with shorter OS and higher PCSM (p = 0.015 and p = 0.001, respectively). Negative/weak ATM expression was significantly associated with PCSM with a hazard ratio of 2.09 (95% confidence interval 1.34-3.27, p = 0.001). Assessment of Combined ATM/PTEN expression showed improved prognostic power to predict OS and PCSM, independent of Gleason grade groups. CONCLUSIONS: Decreased ATM protein expression is associated with poor outcomes in advanced PCa patients. Patients with combined low ATM/PTEN negative expression are at the highest risk for reduced OS and PCSM. Assessing the combined status of ATM/PTEN by IHC in PCa patients may aid in risk stratification relative to OS and PCSM. Moreover, since ATM plays an integral role in DNA damage response pathways, future studies will enhance our understanding of how outcomes of patients with altered ATM and PTEN expression can be improved further with poly-ADP ribose polymerase inhibitors (PARPi), combinations of PARPi and androgen receptor-targeted therapies, as well as platinum-based chemotherapies. PATIENT SUMMARY: Lower ATM intensity is associated with increased cancer-specific mortality in prostate cancer patients. Patients with lower ATM and PTEN negative expression showed decreased overall survival and increased cancer mortality compared with controls.

Copy Number Profiles of Prostate Cancer in Men of Middle Eastern Ancestry.

Our knowledge of prostate cancer (PCa) genomics mainly reflects European (EUR) and Asian (ASN) populations. Our understanding of the influence of Middle Eastern (ME) and African (AFR) ancestry on the mutational profiles of prostate cancer is limited. To characterize genomic differences between ME, EUR, ASN, and AFR ancestry, fluorescent in situ hybridization (FISH) studies for NKX3-1 deletion and MYC amplification were carried out on 42 tumors arising in individuals of ME ancestry. These were supplemented by analysis of genome-wide copy number profiles of 401 tumors of all ancestries. FISH results of NKX3-1 and MYC were assessed in the ME cohort and compared to other ancestries. Gene level copy number aberrations (CNAs) for each sample were statistically compared between ancestry groups. NKX3-1 deletions by FISH were observed in 17/42 (17.5%) prostate tumors arising in men of ME ancestry, while MYC amplifications were only observed in 1/42 (2.3%). Using CNAs called from arrays, the incidence of NKX3-1 deletions was significantly lower in ME vs. other ancestries (20% vs. 52%; p = 2.3 × 10-3). Across the genome, tumors arising in men of ME ancestry had fewer CNAs than those in men of other ancestries (p = 0.014). Additionally, the somatic amplification of 21 specific genes was more frequent in tumors arising in men of ME vs. EUR ancestry (two-sided proportion test; Q < 0.05). Those included amplifications in the glutathione S-transferase family on chromosome 1 (GSTM1, GSTM2, GSTM5) and the IQ motif-containing family on chromosome 3 (IQCF1, IQCF2, IQCF13, IQCF4, IQCF5, IQCF6). Larger studies investigating ME populations are warranted to confirm these observations.

Immune-focused multi-omics analysis of prostate cancer: leukocyte Ig-Like receptors are associated with disease progression.

Prostate cancer (PCa) immunotherapy has shown limited efficacy so far, even in advanced-stage cancers. The success rate of PCa immunotherapy might be improved by approaches more adapted to the immunobiology of the disease. The objective of this study was to perform a multi-omics analysis to identify immune genes associated with PCa progression to better characterize PCa immunobiology and propose new immunotherapeutic targets. mRNA, miRNA, methylation, copy number aberration, and single nucleotide variant datasets from The Cancer Genome Atlas PRAD cohort were analyzed after filtering for genes associated with immunity. Sparse partial least squares-discriminant analyses were performed to identify features associated with biochemical recurrence (BCR) in each type of omics data. Selected features predicted BCR with a balanced error rate (BER) of 0.20 to 0.51 in single-omics and of 0.05 in multi-omics analyses. Amongst features associated with BCR were genes from the Immunoglobulin Ig-like Receptor (LILR) family which are immune checkpoints with immunotherapeutic potential. Using Multivariate INTegrative (MINT) analysis, the association of five LILR genes with BCR was quantified in a combination of three RNA-seq datasets and confirmed with Kaplan-Meier analysis in both these and in an independent RNA-seq dataset. Finally, immunohistochemistry showed that a high number of LILRB1 positive cells within the tumors predicted long-term adverse outcomes. Thus, tumors characterized by abnormal expression of LILR genes have an elevated risk of recurring after definitive local therapy. The immunotherapeutic potential of these regulators to stimulate the immune response against PCa should be evaluated in pre-clinical models.

Transcriptomic Heterogeneity of Gleason Grade Group 5 Prostate Cancer.

Gleason grade group (GG) 5 prostate cancer has been associated with an aggressive natural history, and retrospective data support a role for treatment intensification. However, clinical outcomes remain heterogeneous in this cohort, and intensified treatments carry an increased risk of adverse events. We sought to explore the transcriptomic heterogeneity of GG 5 tumors by querying transcriptomic data from the tumors of 2138 patients with GG 5 disease who underwent prostatectomy. Four distinct consensus clusters were identified with respect to differential transcriptional activation of hallmark pathways, with distinct molecular subtyping profiles and different average genomic risks (AGRs). One cluster, accounting for 325 tumors (15.2% of the population), was enriched for genes related to the cell cycle/proliferation, metabolic pathways, androgen response pathways, and DNA repair, and had a higher AGR than the other clusters (p < 0.001). This clustering, with an identification of a high genomic risk cluster, was subsequently validated in a separate cohort of 1921 patients as well as a third cohort of 201 patients. The latter cohort had outcomes available, and it was found that patients in the high genomic risk cluster had significantly worse distant metastasis-free survival than the other clusters. Tumors in this high genomic risk cluster of GG 5 disease may be particularly likely to benefit from treatment intensification. PATIENT SUMMARY: In this report, we examined differences in gene expression in tumors from men with Gleason grade group 5 prostate cancer. We identified significant diversity, with one specific subgroup of tumors associated with expression profiles that suggest a worse prognosis.

Noncoding mutations target cis-regulatory elements of the FOXA1 plexus in prostate cancer.

Prostate cancer is the second most commonly diagnosed malignancy among men worldwide. Recurrently mutated in primary and metastatic prostate tumors, FOXA1 encodes a pioneer transcription factor involved in disease onset and progression through both androgen receptor-dependent and androgen receptor-independent mechanisms. Despite its oncogenic properties however, the regulation of FOXA1 expression remains unknown. Here, we identify a set of six cis-regulatory elements in the FOXA1 regulatory plexus harboring somatic single-nucleotide variants in primary prostate tumors. We find that deletion and repression of these cis-regulatory elements significantly decreases FOXA1 expression and prostate cancer cell growth. Six of the ten single-nucleotide variants mapping to FOXA1 regulatory plexus significantly alter the transactivation potential of cis-regulatory elements by modulating the binding of transcription factors. Collectively, our results identify cis-regulatory elements within the FOXA1 plexus mutated in primary prostate tumors as potential targets for therapeutic intervention.

The genomic landscape of metastatic castration-resistant prostate cancers reveals multiple distinct genotypes with potential clinical impact.

Metastatic castration-resistant prostate cancer (mCRPC) has a highly complex genomic landscape. With the recent development of novel treatments, accurate stratification strategies are needed. Here we present the whole-genome sequencing (WGS) analysis of fresh-frozen metastatic biopsies from 197 mCRPC patients. Using unsupervised clustering based on genomic features, we define eight distinct genomic clusters. We observe potentially clinically relevant genotypes, including microsatellite instability (MSI), homologous recombination deficiency (HRD) enriched with genomic deletions and BRCA2 aberrations, a tandem duplication genotype associated with CDK12-/- and a chromothripsis-enriched subgroup. Our data suggests that stratification on WGS characteristics may improve identification of MSI, CDK12-/- and HRD patients. From WGS and ChIP-seq data, we show the potential relevance of recurrent alterations in non-coding regions identified with WGS and highlight the central role of AR signaling in tumor progression. These data underline the potential value of using WGS to accurately stratify mCRPC patients into clinically actionable subgroups.

Genome-wide germline correlates of the epigenetic landscape of prostate cancer.

Oncogenesis is driven by germline, environmental and stochastic factors. It is unknown how these interact to produce the molecular phenotypes of tumors. We therefore quantified the influence of germline polymorphisms on the somatic epigenome of 589 localized prostate tumors. Predisposition risk loci influence a tumor's epigenome, uncovering a mechanism for cancer susceptibility. We identified and validated 1,178 loci associated with altered methylation in tumoral but not nonmalignant tissue. These tumor methylation quantitative trait loci influence chromatin structure, as well as RNA and protein abundance. One prominent tumor methylation quantitative trait locus is associated with AKT1 expression and is predictive of relapse after definitive local therapy in both discovery and validation cohorts. These data reveal intricate crosstalk between the germ line and the epigenome of primary tumors, which may help identify germline biomarkers of aggressive disease to aid patient triage and optimize the use of more invasive or expensive diagnostic assays.

The Proteogenomic Landscape of Curable Prostate Cancer.

DNA sequencing has identified recurrent mutations that drive the aggressiveness of prostate cancers. Surprisingly, the influence of genomic, epigenomic, and transcriptomic dysregulation on the tumor proteome remains poorly understood. We profiled the genomes, epigenomes, transcriptomes, and proteomes of 76 localized, intermediate-risk prostate cancers. We discovered that the genomic subtypes of prostate cancer converge on five proteomic subtypes, with distinct clinical trajectories. ETS fusions, the most common alteration in prostate tumors, affect different genes and pathways in the proteome and transcriptome. Globally, mRNA abundance changes explain only ∼10% of protein abundance variability. As a result, prognostic biomarkers combining genomic or epigenomic features with proteomic ones significantly outperform biomarkers comprised of a single data type.

Widespread and Functional RNA Circularization in Localized Prostate Cancer.

The cancer transcriptome is remarkably complex, including low-abundance transcripts, many not polyadenylated. To fully characterize the transcriptome of localized prostate cancer, we performed ultra-deep total RNA-seq on 144 tumors with rich clinical annotation. This revealed a linear transcriptomic subtype associated with the aggressive intraductal carcinoma sub-histology and a fusion profile that differentiates localized from metastatic disease. Analysis of back-splicing events showed widespread RNA circularization, with the average tumor expressing 7,232 circular RNAs (circRNAs). The degree of circRNA production was correlated to disease progression in multiple patient cohorts. Loss-of-function screening identified 11.3% of highly abundant circRNAs as essential for cell proliferation; for ∼90% of these, their parental linear transcripts were not essential. Individual circRNAs can have distinct functions, with circCSNK1G3 promoting cell growth by interacting with miR-181. These data advocate for adoption of ultra-deep RNA-seq without poly-A selection to interrogate both linear and circular transcriptomes.

Molecular Hallmarks of Multiparametric Magnetic Resonance Imaging Visibility in Prostate Cancer.

Multiparametric magnetic resonance imaging (mpMRI) has transformed the management of localized prostate cancer by improving identification of clinically significant disease at diagnosis. Approximately 20% of primary prostate tumors are invisible to mpMRI, and we hypothesize that this invisibility reflects fundamental molecular properties of the tumor. We therefore profiled the genomes and transcriptomes of 40 International Society of Urological Pathology grade 2 tumors: 20 mpMRI-invisible (Prostate Imaging-Reporting and Data System [PI-RADS] v2 <3) and 20 mpMRI-visible (PI-RADS v2 5) tumors. mpMRI-visible tumors were enriched in hallmarks of nimbosus, an aggressive pathological, molecular, and microenvironmental phenomenon in prostate cancer. These hallmarks included genomes with increased mutation density, a higher prevalence of intraductal carcinoma/cribriform architecture pathology, and altered abundance of 102 transcripts, including overexpression of noncoding RNAs such as SCHLAP1. Multiple small nucleolar RNAs (snoRNAs) were identified, and a snoRNA signature synergized with nimbosus hallmarks to discriminate visible from invisible tumors. These data suggest a confluence of aggressive molecular and microenvironmental phenomena underlie mpMRI visibility of localized prostate cancer. PATIENT SUMMARY: We examined the correlation between tumor biology and magnetic resonance imaging (MRI) visibility in a group of patients with low- intermediate-risk prostate cancer. We observed that MRI findings are associated with biological features of aggressive prostate cancer.

Molecular landmarks of tumor hypoxia across cancer types.

Many primary-tumor subregions have low levels of molecular oxygen, termed hypoxia. Hypoxic tumors are at elevated risk for local failure and distant metastasis, but the molecular hallmarks of tumor hypoxia remain poorly defined. To fill this gap, we quantified hypoxia in 8,006 tumors across 19 tumor types. In ten tumor types, hypoxia was associated with elevated genomic instability. In all 19 tumor types, hypoxic tumors exhibited characteristic driver-mutation signatures. We observed widespread hypoxia-associated dysregulation of microRNAs (miRNAs) across cancers and functionally validated miR-133a-3p as a hypoxia-modulated miRNA. In localized prostate cancer, hypoxia was associated with elevated rates of chromothripsis, allelic loss of PTEN and shorter telomeres. These associations are particularly enriched in polyclonal tumors, representing a constellation of features resembling tumor nimbosus, an aggressive cellular phenotype. Overall, this work establishes that tumor hypoxia may drive aggressive molecular features across cancers and shape the clinical trajectory of individual tumors.