Rare germline structural variants increase risk for pediatric solid tumors.

Pediatric solid tumors are rare malignancies that represent a leading cause of death by disease among children in developed countries. The early age-of-onset of these tumors suggests that germline genetic factors are involved, yet conventional germline testing for short coding variants in established predisposition genes only identifies pathogenic events in 10-15% of patients. Here, we examined the role of germline structural variants (SVs)-an underexplored form of germline variation-in pediatric extracranial solid tumors using germline genome sequencing of 1,766 affected children, their 943 unaffected relatives, and 6,665 adult controls. We discovered a sex-biased association between very large (>1 megabase) germline chromosomal abnormalities and a four-fold increased risk of solid tumors in male children. The overall impact of germline SVs was greatest in neuroblastoma, where we revealed burdens of ultra-rare SVs that cause loss-of-function of highly expressed, mutationally intolerant, neurodevelopmental genes, as well as noncoding SVs predicted to disrupt three-dimensional chromatin domains in neural crest-derived tissues. Collectively, our results implicate rare germline SVs as a predisposing factor to pediatric solid tumors that may guide future studies and clinical practice.

Somatic structural variants drive distinct modes of oncogenesis in melanoma.

The diversity of structural variants (SVs) in melanoma and how they impact oncogenesis are incompletely known. We performed harmonized analysis of SVs across melanoma histological and genomic subtypes, and we identified distinct global properties between subtypes. These included the frequency and size of SVs and SV classes, their relation to chromothripsis events, and the role of topologically associated domain (TAD) boundary altering SVs on cancer-related genes. Following our prior identification of double-stranded break repair deficiency in a subset of triple wild-type cutaneous melanoma, we identified MRE11 and NBN loss-of-function SVs in melanomas with this mutational signature. Experimental knockouts of MRE11 and NBN, followed by olaparib cell viability assays in melanoma cells, indicated that dysregulation of each of these genes may cause sensitivity to PARPi in cutaneous melanomas. Broadly, harmonized analysis of melanoma SVs revealed distinct global genomic properties and molecular drivers, which may have biological and therapeutic impact.

Overcoming Barriers to Tumor Genomic Profiling through Direct-to-Patient Outreach.

PURPOSE: To overcome barriers to genomic testing for patients with rare cancers, we initiated a program to offer free clinical tumor genomic testing worldwide to patients with select rare cancer subtypes. EXPERIMENTAL DESIGN: Patients were recruited through social media outreach and engagement with disease-specific advocacy groups, with a focus on patients with histiocytosis, germ cell tumors (GCT), and pediatric cancers. Tumors were analyzed using the MSK-IMPACT next-generation sequencing assay with the return of results to patients and their local physicians. Whole-exome recapture was performed for female patients with GCTs to define the genomic landscape of this rare cancer subtype. RESULTS: A total of 333 patients were enrolled, and tumor tissue was received for 288 (86.4%), with 250 (86.8%) having tumor DNA of sufficient quality for MSK-IMPACT testing. Eighteen patients with histiocytosis have received genomically guided therapy to date, of whom 17 (94%) have had clinical benefit with a mean treatment duration of 21.7 months (range, 6-40+). Whole-exome sequencing of ovarian GCTs identified a subset with haploid genotypes, a phenotype rarely observed in other cancer types. Actionable genomic alterations were rare in ovarian GCT (28%); however, 2 patients with ovarian GCTs with squamous transformation had high tumor mutational burden, one of whom had a complete response to pembrolizumab. CONCLUSIONS: Direct-to-patient outreach can facilitate the assembly of cohorts of rare cancers of sufficient size to define their genomic landscape. By profiling tumors in a clinical laboratory, results could be reported to patients and their local physicians to guide treatment. See related commentary by Desai and Subbiah, p. 2339.

Analysis of evolutionary dynamics and clonal architecture in prostate cancer.

The extent to which clinical and genomic characteristics associate with prostate cancer clonal architecture, tumor evolution, and therapeutic response remains unclear. Here, we reconstructed the clonal architecture and evolutionary trajectories of 845 prostate cancer tumors with harmonized clinical and molecular data. We observed that tumors from patients who self-reported as Black had more linear and monoclonal architectures, despite these men having higher rates of biochemical recurrence. This finding contrasts with prior observations relating polyclonal architecture to adverse clinical outcomes. Additionally, we utilized a novel approach to mutational signature analysis that leverages clonal architecture to uncover additional cases of homologous recombination and mismatch repair deficiency in primary and metastatic tumors and link the origin of mutational signatures to specific subclones. Broadly, prostate cancer clonal architecture analysis reveals novel biological insights that may be immediately clinically actionable and provide multiple opportunities for subsequent investigation. Statement of significance: Tumors from patients who self-reported as Black demonstrate linear and monoclonal evolutionary trajectories yet experience higher rates of biochemical recurrence. In addition, analysis of clonal and subclonal mutational signatures identifies additional tumors with potentially actionable alterations such as deficiencies in mismatch repair and homologous recombination.

Glutathione supports lipid abundance in vivo.

Cells rely on antioxidants to survive. The most abundant antioxidant is glutathione (GSH). The synthesis of GSH is non-redundantly controlled by the glutamate-cysteine ligase catalytic subunit (GCLC). GSH imbalance is implicated in many diseases, but the requirement for GSH in adult tissues is unclear. To interrogate this, we developed a series of in vivo models to induce Gclc deletion in adult animals. We find that GSH is essential to lipid abundance in vivo. GSH levels are reported to be highest in liver tissue, which is also a hub for lipid production. While the loss of GSH did not cause liver failure, it decreased lipogenic enzyme expression, circulating triglyceride levels, and fat stores. Mechanistically, we found that GSH promotes lipid abundance by repressing NRF2, a transcription factor induced by oxidative stress. These studies identify GSH as a fulcrum in the liver's balance of redox buffering and triglyceride production.

A patient-driven clinicogenomic partnership for metastatic prostate cancer.

Molecular profiling studies have enabled discoveries for metastatic prostate cancer (MPC) but have predominantly occurred in academic medical institutions and involved non-representative patient populations. We established the Metastatic Prostate Cancer Project (MPCproject, mpcproject.org), a patient-partnered initiative to involve patients with MPC living anywhere in the US and Canada in molecular research. Here, we present results from our partnership with the first 706 MPCproject participants. While 41% of patient partners live in rural, physician-shortage, or medically underserved areas, the MPCproject has not yet achieved racial diversity, a disparity that demands new initiatives detailed herein. Among molecular data from 333 patient partners (572 samples), exome sequencing of 63 tumor and 19 cell-free DNA (cfDNA) samples recapitulated known findings in MPC, while inexpensive ultra-low-coverage sequencing of 318 cfDNA samples revealed clinically relevant AR amplifications. This study illustrates the power of a growing, longitudinal partnership with patients to generate a more representative understanding of MPC.

Clonal populations of a human TNBC model display significant functional heterogeneity and divergent growth dynamics in distinct contexts.

Intratumoral heterogeneity has been described for various tumor types and models of human cancer, and can have profound effects on tumor progression and drug resistance. This study describes an in-depth analysis of molecular and functional heterogeneity among subclonal populations (SCPs) derived from a single triple-negative breast cancer cell line, including copy number analysis, whole-exome and RNA sequencing, proteome analysis, and barcode analysis of clonal dynamics, as well as functional assays. The SCPs were found to have multiple unique genetic alterations and displayed significant variation in anchorage independent growth and tumor forming ability. Analyses of clonal dynamics in SCP mixtures using DNA barcode technology revealed selection for distinct clonal populations in different in vitro and in vivo environmental contexts, demonstrating that in vitro propagation of cancer cell lines using different culture conditions can contribute to the establishment of unique strains. These analyses also revealed strong enrichment of a single SCP during the development of xenograft tumors in immune-compromised mice. This SCP displayed attenuated interferon signaling in vivo and reduced sensitivity to the antiproliferative effects of type I interferons. Reduction in interferon signaling was found to provide a selective advantage within the xenograft microenvironment specifically. In concordance with the previously described role of interferon signaling as tumor suppressor, these findings suggest that similar selective pressures may be operative in human cancer and patient-derived xenograft models.

Subtype heterogeneity and epigenetic convergence in neuroendocrine prostate cancer.

Neuroendocrine carcinomas (NEC) are tumors expressing markers of neuronal differentiation that can arise at different anatomic sites but have strong histological and clinical similarities. Here we report the chromatin landscapes of a range of human NECs and show convergence to the activation of a common epigenetic program. With a particular focus on treatment emergent neuroendocrine prostate cancer (NEPC), we analyze cell lines, patient-derived xenograft (PDX) models and human clinical samples to show the existence of two distinct NEPC subtypes based on the expression of the neuronal transcription factors ASCL1 and NEUROD1. While in cell lines and PDX models these subtypes are mutually exclusive, single-cell analysis of human clinical samples exhibits a more complex tumor structure with subtypes coexisting as separate sub-populations within the same tumor. These tumor sub-populations differ genetically and epigenetically contributing to intra- and inter-tumoral heterogeneity in human metastases. Overall, our results provide a deeper understanding of the shared clinicopathological characteristics shown by NECs. Furthermore, the intratumoral heterogeneity of human NEPCs suggests the requirement of simultaneous targeting of coexisting tumor populations as a therapeutic strategy.

Biologically informed deep neural network for prostate cancer discovery.

The determination of molecular features that mediate clinically aggressive phenotypes in prostate cancer remains a major biological and clinical challenge1,2. Recent advances in interpretability of machine learning models as applied to biomedical problems may enable discovery and prediction in clinical cancer genomics3-5. Here we developed P-NET-a biologically informed deep learning model-to stratify patients with prostate cancer by treatment-resistance state and evaluate molecular drivers of treatment resistance for therapeutic targeting through complete model interpretability. We demonstrate that P-NET can predict cancer state using molecular data with a performance that is superior to other modelling approaches. Moreover, the biological interpretability within P-NET revealed established and novel molecularly altered candidates, such as MDM4 and FGFR1, which were implicated in predicting advanced disease and validated in vitro. Broadly, biologically informed fully interpretable neural networks enable preclinical discovery and clinical prediction in prostate cancer and may have general applicability across cancer types.

Molecular features of exceptional response to neoadjuvant anti-androgen therapy in high-risk localized prostate cancer.

High-risk localized prostate cancer (HRLPC) is associated with a substantial risk of recurrence and disease mortality. Recent clinical trials have shown that intensifying anti-androgen therapies administered before prostatectomy can induce pathologic complete responses or minimal residual disease, called exceptional response, although the molecular determinants of these clinical outcomes are largely unknown. Here, we perform whole-exome and transcriptome sequencing on pre-treatment multi-regional tumor biopsies from exceptional responders (ERs) and non-responders (NRs, pathologic T3 or lymph node-positive disease) to intensive neoadjuvant anti-androgen therapies. Clonal SPOP mutation and SPOPL copy-number loss are exclusively observed in ERs, while clonal TP53 mutation and PTEN copy-number loss are exclusively observed in NRs. Transcriptional programs involving androgen signaling and TGF-ß signaling are enriched in ERs and NRs, respectively. These findings may guide prospective validation studies of these molecular features in large HRLPC clinical cohorts treated with neoadjuvant anti-androgens to improve patient stratification.

Gene Fusions Create Partner and Collateral Dependencies Essential to Cancer Cell Survival.

Gene fusions frequently result from rearrangements in cancer genomes. In many instances, gene fusions play an important role in oncogenesis; in other instances, they are thought to be passenger events. Although regulatory element rearrangements and copy number alterations resulting from these structural variants are known to lead to transcriptional dysregulation across cancers, the extent to which these events result in functional dependencies with an impact on cancer cell survival is variable. Here we used CRISPR-Cas9 dependency screens to evaluate the fitness impact of 3,277 fusions across 645 cell lines from the Cancer Dependency Map. We found that 35% of cell lines harbored either a fusion partner dependency or a collateral dependency on a gene within the same topologically associating domain as a fusion partner. Fusion-associated dependencies revealed numerous novel oncogenic drivers and clinically translatable alterations. Broadly, fusions can result in partner and collateral dependencies that have biological and clinical relevance across cancer types. SIGNIFICANCE: This study provides insights into how fusions contribute to fitness in different cancer contexts beyond partner-gene activation events, identifying partner and collateral dependencies that may have direct implications for clinical care.

Phase II Multicenter Study of Enzalutamide in Metastatic Castration-Resistant Prostate Cancer to Identify Mechanisms Driving Resistance.

PURPOSE: Enzalutamide is a second-generation androgen receptor (AR) inhibitor that has improved overall survival (OS) in metastatic castration-resistant prostate cancer (CRPC). However, nearly all patients develop resistance. We designed a phase II multicenter study of enzalutamide in metastatic CRPC incorporating tissue and blood biomarkers to dissect mechanisms driving resistance. PATIENTS AND METHODS: Eligible patients with metastatic CRPC underwent a baseline metastasis biopsy and then initiated enzalutamide 160 mg daily. A repeat metastasis biopsy was obtained at radiographic progression from the same site when possible. Blood for circulating tumor cell (CTC) analysis was collected at baseline and progression. The primary objective was to analyze mechanisms of resistance in serial biopsies. Whole-exome sequencing was performed on tissue biopsies. CTC samples underwent RNA sequencing. RESULTS: A total of 65 patients initiated treatment, of whom 22 (33.8%) had received prior abiraterone. Baseline biopsies were enriched for alterations in AR (mutations, amplifications) and tumor suppression genes (PTEN, RB1, and TP53), which were observed in 73.1% and 92.3% of baseline biopsies, respectively. Progression biopsies revealed increased AR amplifications (64.7% at progression vs. 53.9% at baseline) and BRCA2 alterations (64.7% at progression vs. 38.5% at baseline). Genomic analysis of baseline and progression CTC samples demonstrated increased AR splice variants, AR-regulated genes, and neuroendocrine markers at progression. CONCLUSIONS: Our results demonstrate that a large proportion of enzalutamide-treated patients have baseline and progression alterations in the AR pathway and tumor suppressor genes. We demonstrate an increased number of BRCA2 alterations post-enzalutamide, highlighting the importance of serial tumor sampling in CRPC.

Transcriptional mediators of treatment resistance in lethal prostate cancer.

Metastatic castration-resistant prostate cancer is typically lethal, exhibiting intrinsic or acquired resistance to second-generation androgen-targeting therapies and minimal response to immune checkpoint inhibitors1. Cellular programs driving resistance in both cancer and immune cells remain poorly understood. We present single-cell transcriptomes from 14 patients with advanced prostate cancer, spanning all common metastatic sites. Irrespective of treatment exposure, adenocarcinoma cells pervasively coexpressed multiple androgen receptor isoforms, including truncated isoforms hypothesized to mediate resistance to androgen-targeting therapies2,3. Resistance to enzalutamide was associated with cancer cell-intrinsic epithelial-mesenchymal transition and transforming growth factor-ß signaling. Small cell carcinoma cells exhibited divergent expression programs driven by transcriptional regulators promoting lineage plasticity and HOXB5, HOXB6 and NR1D2 (refs. 4-6). Additionally, a subset of patients had high expression of dysfunction markers on cytotoxic CD8+ T cells undergoing clonal expansion following enzalutamide treatment. Collectively, the transcriptional characterization of cancer and immune cells from human metastatic castration-resistant prostate cancer provides a basis for the development of therapeutic approaches complementing androgen signaling inhibition.

CoMut: visualizing integrated molecular information with comutation plots.

MOTIVATION: Large-scale sequencing studies have created a need to succinctly visualize genomic characteristics of patient cohorts linked to widely variable phenotypic information. This is often done by visualizing the co-occurrence of variants with comutation plots. Current tools lack the ability to create highly customizable and publication quality comutation plots from arbitrary user data. RESULTS: We developed CoMut, a stand-alone, object-oriented Python package that creates comutation plots from arbitrary input data, including categorical data, continuous data, bar graphs, side bar graphs and data that describes relationships between samples. AVAILABILITY AND IMPLEMENTATION: The CoMut package is open source and is available at https://github.com/vanallenlab/comut under the MIT License, along with documentation and examples. A no installation, easy-to-use implementation is available on Google Colab (see GitHub).

Seasonal changes in diet and chemical defense in the Climbing Mantella frog (Mantella laevigata).

Poison frogs acquire chemical defenses from the environment for protection against potential predators. These defensive chemicals are lipophilic alkaloids that are sequestered by poison frogs from dietary arthropods and stored in skin glands. Despite decades of research focusing on identifying poison frog alkaloids, we know relatively little about how environmental variation and subsequent arthropod availability impacts alkaloid loads in poison frogs. We investigated how seasonal environmental variation influences poison frog chemical profiles through changes in the diet of the Climbing Mantella (Mantella laevigata). We collected M. laevigata females on the Nosy Mangabe island reserve in Madagascar during the wet and dry seasons and tested the hypothesis that seasonal differences in rainfall is associated with changes in diet composition and skin alkaloid profiles of M. laevigata. The arthropod diet of each frog was characterized into five groups (i.e. ants, termites, mites, insect larvae, or 'other') using visual identification and cytochrome oxidase 1 DNA barcoding. We found that frog diet differed between the wet and dry seasons, where frogs had a more diverse diet in the wet season and consumed a higher percentage of ants in the dry season. To determine if seasonality was associated with variation in frog defensive chemical composition, we used gas chromatography / mass spectrometry to quantify alkaloids from individual skin samples. Although the assortment of identified alkaloids was similar across seasons, we detected significant differences in the abundance of certain alkaloids, which we hypothesize reflects seasonal variation in the diet of M. laevigata. We suggest that these variations could originate from seasonal changes in either arthropod leaf litter composition or changes in frog behavioral patterns. Although additional studies are needed to understand the consequences of long-term environmental shifts, this work suggests that alkaloid profiles are relatively robust against short-term environmental perturbations.