Pan-cancer proteogenomics connects oncogenic drivers to functional states.

Cancer driver events refer to key genetic aberrations that drive oncogenesis; however, their exact molecular mechanisms remain insufficiently understood. Here, our multi-omics pan-cancer analysis uncovers insights into the impacts of cancer drivers by identifying their significant cis-effects and distal trans-effects quantified at the RNA, protein, and phosphoprotein levels. Salient observations include the association of point mutations and copy-number alterations with the rewiring of protein interaction networks, and notably, most cancer genes converge toward similar molecular states denoted by sequence-based kinase activity profiles. A correlation between predicted neoantigen burden and measured T cell infiltration suggests potential vulnerabilities for immunotherapies. Patterns of cancer hallmarks vary by polygenic protein abundance ranging from uniform to heterogeneous. Overall, our work demonstrates the value of comprehensive proteogenomics in understanding the functional states of oncogenic drivers and their links to cancer development, surpassing the limitations of studying individual cancer types.

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

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

Oncogenic Signaling Pathways in The Cancer Genome Atlas.

Genetic alterations in signaling pathways that control cell-cycle progression, apoptosis, and cell growth are common hallmarks of cancer, but the extent, mechanisms, and co-occurrence of alterations in these pathways differ between individual tumors and tumor types. Using mutations, copy-number changes, mRNA expression, gene fusions and DNA methylation in 9,125 tumors profiled by The Cancer Genome Atlas (TCGA), we analyzed the mechanisms and patterns of somatic alterations in ten canonical pathways: cell cycle, Hippo, Myc, Notch, Nrf2, PI-3-Kinase/Akt, RTK-RAS, TGFß signaling, p53 and ß-catenin/Wnt. We charted the detailed landscape of pathway alterations in 33 cancer types, stratified into 64 subtypes, and identified patterns of co-occurrence and mutual exclusivity. Eighty-nine percent of tumors had at least one driver alteration in these pathways, and 57% percent of tumors had at least one alteration potentially targetable by currently available drugs. Thirty percent of tumors had multiple targetable alterations, indicating opportunities for combination therapy.

Comprehensive Characterization of Cancer Driver Genes and Mutations.

Identifying molecular cancer drivers is critical for precision oncology. Multiple advanced algorithms to identify drivers now exist, but systematic attempts to combine and optimize them on large datasets are few. We report a PanCancer and PanSoftware analysis spanning 9,423 tumor exomes (comprising all 33 of The Cancer Genome Atlas projects) and using 26 computational tools to catalog driver genes and mutations. We identify 299 driver genes with implications regarding their anatomical sites and cancer/cell types. Sequence- and structure-based analyses identified >3,400 putative missense driver mutations supported by multiple lines of evidence. Experimental validation confirmed 60%-85% of predicted mutations as likely drivers. We found that >300 MSI tumors are associated with high PD-1/PD-L1, and 57% of tumors analyzed harbor putative clinically actionable events. Our study represents the most comprehensive discovery of cancer genes and mutations to date and will serve as a blueprint for future biological and clinical endeavors.

Characteristics of The Cancer Genome Atlas cases relative to U.S. general population cancer cases.

BACKGROUND: Despite anecdotal reports of differences in clinical and demographic characteristics of The Cancer Genome Atlas (TCGA) relative to general population cancer cases, differences have not been systematically evaluated. METHODS: Data from 11,160 cases with 33 cancer types were ascertained from TCGA data portal. Corresponding data from the Surveillance, Epidemiology, and End Results (SEER) 18 and North American Association of Central Cancer Registries databases were obtained. Differences in characteristics were compared using Student's t, Chi-square, and Fisher's exact tests. Differences in mean survival months were assessed using restricted mean survival time analysis and generalised linear model. RESULTS: TCGA cases were 3.9 years (95% CI 1.7-6.2) younger on average than SEER cases, with a significantly younger mean age for 20/33 cancer types. Although most cancer types had a similar sex distribution, race and stage at diagnosis distributions were disproportional for 13/18 and 25/26 assessed cancer types, respectively. Using 12 months as an end point, the observed mean survival months were longer for 27 of 33 TCGA cancer types. CONCLUSIONS: Differences exist in the characteristics of TCGA vs. general population cancer cases. Our study highlights population subgroups where increased sample collection is warranted to increase the applicability of cancer genomic research results to all individuals.

Genome-wide association study of CSF levels of 59 alzheimer's disease candidate proteins: significant associations with proteins involved in amyloid processing and inflammation.

Cerebrospinal fluid (CSF) 42 amino acid species of amyloid beta (Aß42) and tau levels are strongly correlated with the presence of Alzheimer's disease (AD) neuropathology including amyloid plaques and neurodegeneration and have been successfully used as endophenotypes for genetic studies of AD. Additional CSF analytes may also serve as useful endophenotypes that capture other aspects of AD pathophysiology. Here we have conducted a genome-wide association study of CSF levels of 59 AD-related analytes. All analytes were measured using the Rules Based Medicine Human DiscoveryMAP Panel, which includes analytes relevant to several disease-related processes. Data from two independently collected and measured datasets, the Knight Alzheimer's Disease Research Center (ADRC) and Alzheimer's Disease Neuroimaging Initiative (ADNI), were analyzed separately, and combined results were obtained using meta-analysis. We identified genetic associations with CSF levels of 5 proteins (Angiotensin-converting enzyme (ACE), Chemokine (C-C motif) ligand 2 (CCL2), Chemokine (C-C motif) ligand 4 (CCL4), Interleukin 6 receptor (IL6R) and Matrix metalloproteinase-3 (MMP3)) with study-wide significant p-values (p<1.46×10-10) and significant, consistent evidence for association in both the Knight ADRC and the ADNI samples. These proteins are involved in amyloid processing and pro-inflammatory signaling. SNPs associated with ACE, IL6R and MMP3 protein levels are located within the coding regions of the corresponding structural gene. The SNPs associated with CSF levels of CCL4 and CCL2 are located in known chemokine binding proteins. The genetic associations reported here are novel and suggest mechanisms for genetic control of CSF and plasma levels of these disease-related proteins. Significant SNPs in ACE and MMP3 also showed association with AD risk. Our findings suggest that these proteins/pathways may be valuable therapeutic targets for AD. Robust associations in cognitively normal individuals suggest that these SNPs also influence regulation of these proteins more generally and may therefore be relevant to other diseases.

Population substructure in Cache County, Utah: the Cache County study.

BACKGROUND: Population stratification is a key concern for genetic association analyses. In addition, extreme homogeneity of ethnic origins of a population can make it difficult to interpret how genetic associations in that population may translate into other populations. Here we have evaluated the genetic substructure of samples from the Cache County study relative to the HapMap Reference populations and data from the Alzheimer's Disease Neuroimaging Initiative (ADNI). RESULTS: Our findings show that the Cache County study is similar in ethnic diversity to the self-reported "Whites" in the ADNI sample and less homogenous than the HapMap CEU population. CONCLUSIONS: We conclude that the Cache County study is genetically representative of the general European American population in the USA and is an appropriate population for conducting broadly applicable genetic studies.

Mitochondrial genomic variation associated with higher mitochondrial copy number: the Cache County Study on Memory Health and Aging.

BACKGROUND: The mitochondria are essential organelles and are the location of cellular respiration, which is responsible for the majority of ATP production. Each cell contains multiple mitochondria, and each mitochondrion contains multiple copies of its own circular genome. The ratio of mitochondrial genomes to nuclear genomes is referred to as mitochondrial copy number. Decreases in mitochondrial copy number are known to occur in many tissues as people age, and in certain diseases. The regulation of mitochondrial copy number by nuclear genes has been studied extensively. While mitochondrial variation has been associated with longevity and some of the diseases known to have reduced mitochondrial copy number, the role that the mitochondrial genome itself has in regulating mitochondrial copy number remains poorly understood. RESULTS: We analyzed the complete mitochondrial genomes from 1007 individuals randomly selected from the Cache County Study on Memory Health and Aging utilizing the inferred evolutionary history of the mitochondrial haplotypes present in our dataset to identify sequence variation and mitochondrial haplotypes associated with changes in mitochondrial copy number. Three variants belonging to mitochondrial haplogroups U5A1 and T2 were significantly associated with higher mitochondrial copy number in our dataset. CONCLUSIONS: We identified three variants associated with higher mitochondrial copy number and suggest several hypotheses for how these variants influence mitochondrial copy number by interacting with known regulators of mitochondrial copy number. Our results are the first to report sequence variation in the mitochondrial genome that causes changes in mitochondrial copy number. The identification of these variants that increase mtDNA copy number has important implications in understanding the pathological processes that underlie these phenotypes.

A Pan-Cancer Patient-Derived Xenograft Histology Image Repository with Genomic and Pathologic Annotations Enables Deep Learning Analysis.

Patient-derived xenografts (PDX) model human intra- and intertumoral heterogeneity in the context of the intact tissue of immunocompromised mice. Histologic imaging via hematoxylin and eosin (H&E) staining is routinely performed on PDX samples, which could be harnessed for computational analysis. Prior studies of large clinical H&E image repositories have shown that deep learning analysis can identify intercellular and morphologic signals correlated with disease phenotype and therapeutic response. In this study, we developed an extensive, pan-cancer repository of >1,000 PDX and paired parental tumor H&E images. These images, curated from the PDX Development and Trial Centers Research Network Consortium, had a range of associated genomic and transcriptomic data, clinical metadata, pathologic assessments of cell composition, and, in several cases, detailed pathologic annotations of neoplastic, stromal, and necrotic regions. The amenability of these images to deep learning was highlighted through three applications: (i) development of a classifier for neoplastic, stromal, and necrotic regions; (ii) development of a predictor of xenograft-transplant lymphoproliferative disorder; and (iii) application of a published predictor of microsatellite instability. Together, this PDX Development and Trial Centers Research Network image repository provides a valuable resource for controlled digital pathology analysis, both for the evaluation of technical issues and for the development of computational image-based methods that make clinical predictions based on PDX treatment studies. Significance: A pan-cancer repository of >1,000 patient-derived xenograft hematoxylin and eosin-stained images will facilitate cancer biology investigations through histopathologic analysis and contributes important model system data that expand existing human histology repositories.

Comparative Genomic Analysis of Pancreatic Acinar Cell Carcinoma (PACC) and Pancreatic Ductal Adenocarcinoma (PDAC) Unveils New Actionable Genomic Aberrations in PACC.

PURPOSE: Pure pancreatic acinar cell carcinomas (PACC) are rare malignancies with no established treatment. PACC demonstrates significant genetic intertumoral heterogeneity with multiple pathways involved, suggesting using targeted cancer therapeutics to treat this disease. We aggregated one of the largest datasets of pure PACC to examine the genomic variability and explore patient-specific therapeutic targets. EXPERIMENTAL DESIGN: PACC specimens (n = 51) underwent next-generation sequencing of DNA (n = 29) or whole exome (n = 22) and RNA (whole transcriptome, n = 29) at a commercial laboratory. We performed comparative analyses of a genomic cohort of pancreatic ductal adenocarcinomas (PDAC; n = 4,205). In parallel, we conducted a retrospective review of patients with PACC treated at Huntsman Cancer Institute (HCI). RESULTS: The real-world dataset included samples from 51 patients with PACC. We found key molecular differences between pure PACC and PDAC, highlighting the unique characteristics of pure PACC. Major differences in PACC include lower MAPK signaling and less stromal cell abundance compared with PDAC. Pure PACC showed genomic loss-of-heterozygosity to largely coincide with mutations in BRCA1, BRCA2, and PALB2. Of the 7 patients treated at HCI, one had a tumor that harbored a BRAF-V600E mutation. Leveraging precision oncology, this patient is being treated with encorafenib plus binimetinib, achieving an exceptionally durable and ongoing complete response of more than 3 years. CONCLUSIONS: There are major differences between PACC and PDAC, including downregulation of the MAPK signaling pathway, and less stromal cell abundance. In addition, genomic characterization of pure PACC revealed frequent targetable alterations, which can guide patient treatment.

Simulation of defects, flexibility and rupture in biopolymer networks.

Networks of biopolymers occur often in nature, and are vulnerable to damage over time. In this work, a coarse grained model of collagen IV molecules is applied in a 2D hexagonal network and the mechanisms by which these networks can rupture are explored. The networks are stretched linearly in order to study their structural limits and mechanism of rupture over timescale of up to 100 µs. Metrics are developed to track the damage networks suffer over time, and qualitatively analyse ruptures that occur. Further simulations repeatedly stretch the networks sinusoidally to mimic the in vivo strains. Defects of increasing levels of complexity are introduced into an ordered network, and their effect on the rupturing behaviour of the biopolymer networks studied. The effect of introducing holes of varying size in the network, as well as strips of finite width to mimic surgical damage are studied. These demonstrate the importance of the flexibility of the networks to preventing damage.

Immunologically "cold" triple negative breast cancers engraft at a higher rate in patient derived xenografts.

TNBC is a heterogeneous subtype of breast cancer, and only a subset of TNBC can be established as PDXs. Here, we show that there is an engraftment bias toward TNBC with low levels of immune cell infiltration. Additionally, TNBC that failed to engraft show gene expression consistent with a cancer-promoting immunological state, leading us to hypothesize that the immunological state of the tumor and possibly the state of the immune system of the host may be essential for engraftment.

A human breast cancer-derived xenograft and organoid platform for drug discovery and precision oncology.

Models that recapitulate the complexity of human tumors are urgently needed to develop more effective cancer therapies. We report a bank of human patient-derived xenografts (PDXs) and matched organoid cultures from tumors that represent the greatest unmet need: endocrine-resistant, treatment-refractory and metastatic breast cancers. We leverage matched PDXs and PDX-derived organoids (PDxO) for drug screening that is feasible and cost-effective with in vivo validation. Moreover, we demonstrate the feasibility of using these models for precision oncology in real time with clinical care in a case of triple-negative breast cancer (TNBC) with early metastatic recurrence. Our results uncovered a Food and Drug Administration (FDA)-approved drug with high efficacy against the models. Treatment with this therapy resulted in a complete response for the individual and a progression-free survival (PFS) period more than three times longer than their previous therapies. This work provides valuable methods and resources for functional precision medicine and drug development for human breast cancer.

Self assembly of model polymers into biological random networks.

The properties of biological networks, such as those found in the ocular lens capsule, are difficult to study without simplified models. Model polymers are developed, inspired by "worm-like" curve models, that are shown to spontaneously self assemble to form networks similar to those observed experimentally in biological systems. These highly simplified coarse-grained models allow the self assembly process to be studied on near-realistic time-scales. Metrics are developed (using a polygon-based framework) which are useful for describing simulated networks and can also be applied to images of real networks. These metrics are used to show the range of control that the computational polymer model has over the networks, including the polygon structure and short range order. The structure of the simulated networks are compared to previous simulation work and microscope images of real networks. The network structure is shown to be a function of the interaction strengths, cooling rates and external pressure. In addition, "pre-tangled" network structures are introduced and shown to significantly influence the subsequent network structure. The network structures obtained fit into a region of the network landscape effectively inaccessible to random (entropically-driven) networks but which are occupied by experimentally-derived configurations.

Moving pan-cancer studies from basic research toward the clinic.

Although all cancers share common hallmarks, we have long realized that there is no silver-bullet treatment for the disease. Many clinical oncologists specialize in a single cancer type, based predominantly on the tissue of origin. With advances brought by genetics and cancer genomic research, we now know that cancers are profoundly different, both in origins and in genetic alterations. At the same time, commonalities such as key driver mutations, altered pathways, mutational, immune and microbial signatures and other areas (many revealed by pan-cancer studies) point to the intriguing possibility of targeting common traits across diverse cancer types with the same therapeutic strategies. Studies designed to delineate differences and similarities across cancer types are thus critical in discerning the basic dynamics of oncogenesis, as well as informing diagnoses, prognoses and therapies. We anticipate growing emphases on the development and application of therapies targeting underlying commonalities of different cancer types, while tailoring to the unique tissue environment and intrinsic molecular fingerprints of each cancer type and subtype. Here we summarize the facets of pan-cancer research and how they are pushing progress toward personalized medicine.

Comprehensive characterization of 536 patient-derived xenograft models prioritizes candidatesfor targeted treatment.

Development of candidate cancer treatments is a resource-intensive process, with the research community continuing to investigate options beyond static genomic characterization. Toward this goal, we have established the genomic landscapes of 536 patient-derived xenograft (PDX) models across 25 cancer types, together with mutation, copy number, fusion, transcriptomic profiles, and NCI-MATCH arms. Compared with human tumors, PDXs typically have higher purity and fit to investigate dynamic driver events and molecular properties via multiple time points from same case PDXs. Here, we report on dynamic genomic landscapes and pharmacogenomic associations, including associations between activating oncogenic events and drugs, correlations between whole-genome duplications and subclone events, and the potential PDX models for NCI-MATCH trials. Lastly, we provide a web portal having comprehensive pan-cancer PDX genomic profiles and source code to facilitate identification of more druggable events and further insights into PDXs' recapitulation of human tumors.

Conservation of copy number profiles during engraftment and passaging of patient-derived cancer xenografts.

Patient-derived xenografts (PDXs) are resected human tumors engrafted into mice for preclinical studies and therapeutic testing. It has been proposed that the mouse host affects tumor evolution during PDX engraftment and propagation, affecting the accuracy of PDX modeling of human cancer. Here, we exhaustively analyze copy number alterations (CNAs) in 1,451 PDX and matched patient tumor (PT) samples from 509 PDX models. CNA inferences based on DNA sequencing and microarray data displayed substantially higher resolution and dynamic range than gene expression-based inferences, and they also showed strong CNA conservation from PTs through late-passage PDXs. CNA recurrence analysis of 130 colorectal and breast PT/PDX-early/PDX-late trios confirmed high-resolution CNA retention. We observed no significant enrichment of cancer-related genes in PDX-specific CNAs across models. Moreover, CNA differences between patient and PDX tumors were comparable to variations in multiregion samples within patients. Our study demonstrates the lack of systematic copy number evolution driven by the PDX mouse host.

Simplified computational model for generating biological networks.

A method to generate and simulate biological networks is discussed. An expanded Wooten-Winer-Weaire bond switching methods is proposed which allows for a distribution of node degrees in the network while conserving the mean average node degree. The networks are characterised in terms of their polygon structure and assortativities (a measure of local ordering). A wide range of experimental images are analysed and the underlying networks quantified in an analogous manner. Limitations in obtaining the network structure are discussed. A "network landscape" of the experimentally observed and simulated networks is constructed from the underlying metrics. The enhanced bond switching algorithm is able to generate networks spanning the full range of experimental observations.

Discovery of driver non-coding splice-site-creating mutations in cancer.

Non-coding mutations can create splice sites, however the true extent of how such somatic non-coding mutations affect RNA splicing are largely unexplored. Here we use the MiSplice pipeline to analyze 783 cancer cases with WGS data and 9494 cases with WES data, discovering 562 non-coding mutations that lead to splicing alterations. Notably, most of these mutations create new exons. Introns associated with new exon creation are significantly larger than the genome-wide average intron size. We find that some mutation-induced splicing alterations are located in genes important in tumorigenesis (ATRX, BCOR, CDKN2B, MAP3K1, MAP3K4, MDM2, SMAD4, STK11, TP53 etc.), often leading to truncated proteins and affecting gene expression. The pattern emerging from these exon-creating mutations suggests that splice sites created by non-coding mutations interact with pre-existing potential splice sites that originally lacked a suitable splicing pair to induce new exon formation. Our study suggests the importance of investigating biological and clinical consequences of noncoding splice-inducing mutations that were previously neglected by conventional annotation pipelines. MiSplice will be useful for automatically annotating the splicing impact of coding and non-coding mutations in future large-scale analyses.