Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Proteogenomic analysis of chemo-refractory high-grade serous ovarian cancer.

To improve the understanding of chemo-refractory high-grade serous ovarian cancers (HGSOCs), we characterized the proteogenomic landscape of 242 (refractory and sensitive) HGSOCs, representing one discovery and two validation cohorts across two biospecimen types (formalin-fixed paraffin-embedded and frozen). We identified a 64-protein signature that predicts with high specificity a subset of HGSOCs refractory to initial platinum-based therapy and is validated in two independent patient cohorts. We detected significant association between lack of Ch17 loss of heterozygosity (LOH) and chemo-refractoriness. Based on pathway protein expression, we identified 5 clusters of HGSOC, which validated across two independent patient cohorts and patient-derived xenograft (PDX) models. These clusters may represent different mechanisms of refractoriness and implicate putative therapeutic vulnerabilities.

A proteogenomic portrait of lung squamous cell carcinoma.

Lung squamous cell carcinoma (LSCC) remains a leading cause of cancer death with few therapeutic options. We characterized the proteogenomic landscape of LSCC, providing a deeper exposition of LSCC biology with potential therapeutic implications. We identify NSD3 as an alternative driver in FGFR1-amplified tumors and low-p63 tumors overexpressing the therapeutic target survivin. SOX2 is considered undruggable, but our analyses provide rationale for exploring chromatin modifiers such as LSD1 and EZH2 to target SOX2-overexpressing tumors. Our data support complex regulation of metabolic pathways by crosstalk between post-translational modifications including ubiquitylation. Numerous immune-related proteogenomic observations suggest directions for further investigation. Proteogenomic dissection of CDKN2A mutations argue for more nuanced assessment of RB1 protein expression and phosphorylation before declaring CDK4/6 inhibition unsuccessful. Finally, triangulation between LSCC, LUAD, and HNSCC identified both unique and common therapeutic vulnerabilities. These observations and proteogenomics data resources may guide research into the biology and treatment of LSCC.

Integrated Proteogenomic Characterization across Major Histological Types of Pediatric Brain Cancer.

We report a comprehensive proteogenomics analysis, including whole-genome sequencing, RNA sequencing, and proteomics and phosphoproteomics profiling, of 218 tumors across 7 histological types of childhood brain cancer: low-grade glioma (n = 93), ependymoma (32), high-grade glioma (25), medulloblastoma (22), ganglioglioma (18), craniopharyngioma (16), and atypical teratoid rhabdoid tumor (12). Proteomics data identify common biological themes that span histological boundaries, suggesting that treatments used for one histological type may be applied effectively to other tumors sharing similar proteomics features. Immune landscape characterization reveals diverse tumor microenvironments across and within diagnoses. Proteomics data further reveal functional effects of somatic mutations and copy number variations (CNVs) not evident in transcriptomics data. Kinase-substrate association and co-expression network analysis identify important biological mechanisms of tumorigenesis. This is the first large-scale proteogenomics analysis across traditional histological boundaries to uncover foundational pediatric brain tumor biology and inform rational treatment selection.

Proteogenomic Characterization Reveals Therapeutic Vulnerabilities in Lung Adenocarcinoma.

To explore the biology of lung adenocarcinoma (LUAD) and identify new therapeutic opportunities, we performed comprehensive proteogenomic characterization of 110 tumors and 101 matched normal adjacent tissues (NATs) incorporating genomics, epigenomics, deep-scale proteomics, phosphoproteomics, and acetylproteomics. Multi-omics clustering revealed four subgroups defined by key driver mutations, country, and gender. Proteomic and phosphoproteomic data illuminated biology downstream of copy number aberrations, somatic mutations, and fusions and identified therapeutic vulnerabilities associated with driver events involving KRAS, EGFR, and ALK. Immune subtyping revealed a complex landscape, reinforced the association of STK11 with immune-cold behavior, and underscored a potential immunosuppressive role of neutrophil degranulation. Smoking-associated LUADs showed correlation with other environmental exposure signatures and a field effect in NATs. Matched NATs allowed identification of differentially expressed proteins with potential diagnostic and therapeutic utility. This proteogenomics dataset represents a unique public resource for researchers and clinicians seeking to better understand and treat lung adenocarcinomas.

Proteogenomic Characterization of Endometrial Carcinoma.

We undertook a comprehensive proteogenomic characterization of 95 prospectively collected endometrial carcinomas, comprising 83 endometrioid and 12 serous tumors. This analysis revealed possible new consequences of perturbations to the p53 and Wnt/ß-catenin pathways, identified a potential role for circRNAs in the epithelial-mesenchymal transition, and provided new information about proteomic markers of clinical and genomic tumor subgroups, including relationships to known druggable pathways. An extensive genome-wide acetylation survey yielded insights into regulatory mechanisms linking Wnt signaling and histone acetylation. We also characterized aspects of the tumor immune landscape, including immunogenic alterations, neoantigens, common cancer/testis antigens, and the immune microenvironment, all of which can inform immunotherapy decisions. Collectively, our multi-omic analyses provide a valuable resource for researchers and clinicians, identify new molecular associations of potential mechanistic significance in the development of endometrial cancers, and suggest novel approaches for identifying potential therapeutic targets.

Integrated Proteogenomic Characterization of Clear Cell Renal Cell Carcinoma.

To elucidate the deregulated functional modules that drive clear cell renal cell carcinoma (ccRCC), we performed comprehensive genomic, epigenomic, transcriptomic, proteomic, and phosphoproteomic characterization of treatment-naive ccRCC and paired normal adjacent tissue samples. Genomic analyses identified a distinct molecular subgroup associated with genomic instability. Integration of proteogenomic measurements uniquely identified protein dysregulation of cellular mechanisms impacted by genomic alterations, including oxidative phosphorylation-related metabolism, protein translation processes, and phospho-signaling modules. To assess the degree of immune infiltration in individual tumors, we identified microenvironment cell signatures that delineated four immune-based ccRCC subtypes characterized by distinct cellular pathways. This study reports a large-scale proteogenomic analysis of ccRCC to discern the functional impact of genomic alterations and provides evidence for rational treatment selection stemming from ccRCC pathobiology.

Mutations in cancer-relevant genes are ubiquitous in histologically normal endometrial tissue.

OBJECTIVE: Endometrial cancer (EndoCA) is the most common gynecologic cancer and incidence and mortality rate continue to increase. Despite well-characterized knowledge of EndoCA-defining mutations, no effective diagnostic or screening tests exist. To lay the foundation for testing development, our study focused on defining the prevalence of somatic mutations present in non-cancerous uterine tissue. METHODS: We obtained ≥8 uterine samplings, including separate endometrial and myometrial layers, from each of 22 women undergoing hysterectomy for non-cancer conditions. We ultra-deep sequenced (>2000× coverage) samples using a 125 cancer-relevant gene panel. RESULTS: All women harbored complex mutation patterns. In total, 308 somatic mutations were identified with mutant allele frequencies ranging up to 96.0%. These encompassed 56 unique mutations from 24 genes. The majority of samples possessed predicted functional cancer mutations but curiously no growth advantage over non-functional mutations was detected. Functional mutations were enriched with increasing patient age (p = 0.045) and BMI (p = 0.0007) and in endometrial versus myometrial layers (68% vs 39%, p = 0.0002). Finally, while the somatic mutation landscape shared similar mutation prevalence in key TCGA-defined EndoCA genes, notably PIK3CA, significant differences were identified, including NOTCH1 (77% vs 10%), PTEN (9% vs 61%), TP53 (0% vs 37%) and CTNNB1 (0% vs 26%). CONCLUSIONS: An important caveat for future liquid biopsy/DNA-based cancer diagnostics is the repertoire of shared and distinct mutation profiles between histologically unremarkable and EndoCA tissues. The lack of selection pressure between functional and non-functional mutations in histologically unremarkable uterine tissue may offer a glimpse into an unrecognized EndoCA protective mechanism.

ProTrack: An Interactive Multi-Omics Data Browser for Proteogenomic Studies.

The Clinical Proteomic Tumor Analysis Consortium (CPTAC) initiative has generated extensive multi-omics data resources of deep proteogenomic profiles for multiple cancer types. To enable the broader community of biological and medical researchers to intuitively query, explore, and download data and analysis results from various CPTAC projects, a prototype user-friendly web application called "ProTrack" is built with the CPTAC clear cell renal cell carcinoma (ccRCC) data set (http://ccrcc.cptac-data-view.org). Here the salient features of this application which provides a dynamic, comprehensive, and granular visualization of the rich proteogenomic data is described.

Nonsyndromic craniosynostosis: novel coding variants.

BACKGROUND: Craniosynostosis (CS), the premature fusion of one or more neurocranial sutures, is associated with approximately 200 syndromes; however, about 65-85% of patients present with no additional major birth defects. METHODS: We conducted targeted next-generation sequencing of 60 known syndromic and other candidate genes in patients with sagittal nonsyndromic CS (sNCS, n = 40) and coronal nonsyndromic CS (cNCS, n = 19). RESULTS: We identified 18 previously published and 5 novel pathogenic variants, including three de novo variants. Novel variants included a paternally inherited c.2209C>G:p.(Leu737Val) variant in BBS9 of a patient with cNCS. Common variants in BBS9, a gene required for ciliogenesis during cranial suture development, have been associated with sNCS risk in a previous genome-wide association study. We also identified c.313G>T:p.(Glu105*) variant in EFNB1 and c.435G>C:p.(Lys145Asn) variant in TWIST1, both in patients with cNCS. Mutations in EFNB1 and TWIST1 have been linked to craniofrontonasal and Saethre-Chotzen syndrome, respectively; both present with coronal CS. CONCLUSIONS: We provide additional evidence that variants in genes implicated in syndromic CS play a role in isolated CS, supporting their inclusion in genetic panels for screening patients with NCS. We also identified a novel BBS9 variant that further shows the potential involvement of BBS9 in the pathogenesis of CS.

A new molecular signature method for prediction of driver cancer pathways from transcriptional data.

Assigning cancer patients to the most effective treatments requires an understanding of the molecular basis of their disease. While DNA-based molecular profiling approaches have flourished over the past several years to transform our understanding of driver pathways across a broad range of tumors, a systematic characterization of key driver pathways based on RNA data has not been undertaken. Here we introduce a new approach for predicting the status of driver cancer pathways based on signature functions derived from RNA sequencing data. To identify the driver cancer pathways of interest, we mined DNA variant data from TCGA and nominated driver alterations in seven major cancer pathways in breast, ovarian and colon cancer tumors. The activation status of these driver pathways were then characterized using RNA sequencing data by constructing classification signature functions in training datasets and then testing the accuracy of the signatures in test datasets. The signature functions differentiate well tumors with nominated pathway activation from tumors with no signs of activation: average AUC equals to 0.83. Our results confirm that driver genomic alterations are distinctively displayed at the transcriptional level and that the transcriptional signatures can generally provide an alternative to DNA sequencing methods in detecting specific driver pathways.

Genomic Analysis of Uterine Lavage Fluid Detects Early Endometrial Cancers and Reveals a Prevalent Landscape of Driver Mutations in Women without Histopathologic Evidence of Cancer: A Prospective Cross-Sectional Study.

BACKGROUND: Endometrial cancer is the most common gynecologic malignancy, and its incidence and associated mortality are increasing. Despite the immediate need to detect these cancers at an earlier stage, there is no effective screening methodology or protocol for endometrial cancer. The comprehensive, genomics-based analysis of endometrial cancer by The Cancer Genome Atlas (TCGA) revealed many of the molecular defects that define this cancer. Based on these cancer genome results, and in a prospective study, we hypothesized that the use of ultra-deep, targeted gene sequencing could detect somatic mutations in uterine lavage fluid obtained from women undergoing hysteroscopy as a means of molecular screening and diagnosis. METHODS AND FINDINGS: Uterine lavage and paired blood samples were collected and analyzed from 107 consecutive patients who were undergoing hysteroscopy and curettage for diagnostic evaluation from this single-institution study. The lavage fluid was separated into cellular and acellular fractions by centrifugation. Cellular and cell-free DNA (cfDNA) were isolated from each lavage. Two targeted next-generation sequencing (NGS) gene panels, one composed of 56 genes and the other of 12 genes, were used for ultra-deep sequencing. To rule out potential NGS-based errors, orthogonal mutation validation was performed using digital PCR and Sanger sequencing. Seven patients were diagnosed with endometrial cancer based on classic histopathologic analysis. Six of these patients had stage IA cancer, and one of these cancers was only detectable as a microscopic focus within a polyp. All seven patients were found to have significant cancer-associated gene mutations in both cell pellet and cfDNA fractions. In the four patients in whom adequate tumor sample was available, all tumor mutations above a specific allele fraction were present in the uterine lavage DNA samples. Mutations originally only detected in lavage fluid fractions were later confirmed to be present in tumor but at allele fractions significantly less than 1%. Of the remaining 95 patients diagnosed with benign or non-cancer pathology, 44 had no significant cancer mutations detected. Intriguingly, 51 patients without histopathologic evidence of cancer had relatively high allele fraction (1.0%-30.4%), cancer-associated mutations. Participants with detected driver and potential driver mutations were significantly older (mean age mutated = 57.96, 95% confidence interval [CI]: 3.30-∞, mean age no mutations = 50.35; p-value = 0.002; Benjamini-Hochberg [BH] adjusted p-value = 0.015) and more likely to be post-menopausal (p-value = 0.004; BH-adjusted p-value = 0.015) than those without these mutations. No associations were detected between mutation status and race/ethnicity, body mass index, diabetes, parity, and smoking status. Long-term follow-up was not presently available in this prospective study for those women without histopathologic evidence of cancer. CONCLUSIONS: Using ultra-deep NGS, we identified somatic mutations in DNA extracted both from cell pellets and a never previously reported cfDNA fraction from the uterine lavage. Using our targeted sequencing approach, endometrial driver mutations were identified in all seven women who received a cancer diagnosis based on classic histopathology of tissue curettage obtained at the time of hysteroscopy. In addition, relatively high allele fraction driver mutations were identified in the lavage fluid of approximately half of the women without a cancer diagnosis. Increasing age and post-menopausal status were associated with the presence of these cancer-associated mutations, suggesting the prevalent existence of a premalignant landscape in women without clinical evidence of cancer. Given that a uterine lavage can be easily and quickly performed even outside of the operating room and in a physician's office-based setting, our findings suggest the future possibility of this approach for screening women for the earliest stages of endometrial cancer. However, our findings suggest that further insight into development of cancer or its interruption are needed before translation to the clinic.

Computational approaches to identify functional genetic variants in cancer genomes.

The International Cancer Genome Consortium (ICGC) aims to catalog genomic abnormalities in tumors from 50 different cancer types. Genome sequencing reveals hundreds to thousands of somatic mutations in each tumor but only a minority of these drive tumor progression. We present the result of discussions within the ICGC on how to address the challenge of identifying mutations that contribute to oncogenesis, tumor maintenance or response to therapy, and recommend computational techniques to annotate somatic variants and predict their impact on cancer phenotype.

Necdin, a p53 target gene, regulates the quiescence and response to genotoxic stress of hematopoietic stem/progenitor cells.

We recently defined a critical role for p53 in regulating the quiescence of adult hematopoietic stem cells (HSCs) and identified necdin as a candidate p53 target gene. Necdin is a growth-suppressing protein and the gene encoding it is one of several that are deleted in patients with Prader-Willi syndrome. To define the intrinsic role of necdin in adult hematopoiesis, in the present study, we transplanted necdin-null fetal liver cells into lethally irradiated recipients. We show that necdin-null adult HSCs are less quiescent and more proliferative than normal HSCs, demonstrating the similar role of necdin and p53 in promoting HSC quiescence during steady-state conditions. However, wild-type recipients repopulated with necdin-null hematopoietic stem/progenitor cells show enhanced sensitivity to irradiation and chemotherapy, with increased p53-dependent apoptosis, myelosuppression, and mortality. Necdin controls the HSC response to genotoxic stress via both cell-cycle-dependent and cell-cycle-independent mechanisms, with the latter occurring in a Gas2L3-dependent manner. We conclude that necdin functions as a molecular switch in adult hematopoiesis, acting in a p53-like manner to promote HSC quiescence in the steady state, but suppressing p53-dependent apoptosis in response to genotoxic stress.

The molecular diversity of Luminal A breast tumors.

Breast cancer is a collection of diseases with distinct molecular traits, prognosis, and therapeutic options. Luminal A breast cancer is the most heterogeneous, both molecularly and clinically. Using genomic data from over 1,000 Luminal A tumors from multiple studies, we analyzed the copy number and mutational landscape of this tumor subtype. This integrated analysis revealed four major subtypes defined by distinct copy-number and mutation profiles. We identified an atypical Luminal A subtype characterized by high genomic instability, TP53 mutations, and increased Aurora kinase signaling; these genomic alterations lead to a worse clinical prognosis. Aberrations of chromosomes 1, 8, and 16, together with PIK3CA, GATA3, AKT1, and MAP3K1 mutations drive the other subtypes. Finally, an unbiased pathway analysis revealed multiple rare, but mutually exclusive, alterations linked to loss of activity of co-repressor complexes N-Cor and SMRT. These rare alterations were the most prevalent in Luminal A tumors and may predict resistance to endocrine therapy. Our work provides for a further molecular stratification of Luminal A breast tumors, with potential direct clinical implications.

Predicting the functional impact of protein mutations: application to cancer genomics.

As large-scale re-sequencing of genomes reveals many protein mutations, especially in human cancer tissues, prediction of their likely functional impact becomes important practical goal. Here, we introduce a new functional impact score (FIS) for amino acid residue changes using evolutionary conservation patterns. The information in these patterns is derived from aligned families and sub-families of sequence homologs within and between species using combinatorial entropy formalism. The score performs well on a large set of human protein mutations in separating disease-associated variants (∼19 200), assumed to be strongly functional, from common polymorphisms (∼35 600), assumed to be weakly functional (area under the receiver operating characteristic curve of ∼0.86). In cancer, using recurrence, multiplicity and annotation for ∼10 000 mutations in the COSMIC database, the method does well in assigning higher scores to more likely functional mutations ('drivers'). To guide experimental prioritization, we report a list of about 1000 top human cancer genes frequently mutated in one or more cancer types ranked by likely functional impact; and, an additional 1000 candidate cancer genes with rare but likely functional mutations. In addition, we estimate that at least 5% of cancer-relevant mutations involve switch of function, rather than simply loss or gain of function.

The most likely but largely ignored triggering factor for breast (or all) cancer invasion.

Breast cancer development and progression are believed to be a sequential process, from normal to hyperplastic, to in situ, and to invasive and metastatic stages. Given that over 90% of cancer deaths are caused by invasive and metastatic lesions, countless factors and multiple theories have been proposed as the triggering factor for the cascade of actions of cancer invasion. However, those factors and theories are largely based on the studies of cell lines or animal models. In addition, corresponding interventions based on these factors and theories have failed to reduce the incidence rate of invasive and metastatic lesions, suggesting that previous efforts may have failed to arm at the right target. Considering these facts and observations, we are proposing "A focal aberrant degeneration in the myoepithelial cell layer (MECL) as the most likely triggering factor for breast cancer invasion". Our hypothesis is based on our recent studies of breast and multiple other cancers. Our commentary provides the rationale, morphologic, immunohistochemical, and molecular data to support our hypotheses. As all epithelium-derived cancers share a very similar architecture, our hypothesis is likely to be applicable to invasion of all cancer types. We believe that human tissue-derived data may provide a more realistic roadmap to guide the clinic practice.

Proteogenomic data and resources for pan-cancer analysis.

The National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium (CPTAC) investigates tumors from a proteogenomic perspective, creating rich multi-omics datasets connecting genomic aberrations to cancer phenotypes. To facilitate pan-cancer investigations, we have generated harmonized genomic, transcriptomic, proteomic, and clinical data for >1000 tumors in 10 cohorts to create a cohesive and powerful dataset for scientific discovery. We outline efforts by the CPTAC pan-cancer working group in data harmonization, data dissemination, and computational resources for aiding biological discoveries. We also discuss challenges for multi-omics data integration and analysis, specifically the unique challenges of working with both nucleotide sequencing and mass spectrometry proteomics data.