Germline pathogenic variants in neuroblastoma patients are enriched in BARD1 and predict worse survival.

BACKGROUND: Neuroblastoma is an embryonal cancer of the developing sympathetic nervous system. The genetic contribution of rare pathogenic or likely pathogenic germline variants in patients without a family history remains unclear. METHODS: Germline DNA sequencing was performed on 786 neuroblastoma patients. The frequency of rare cancer predisposition gene pathogenic or likely pathogenic variants in patients was compared with 2 cancer-free control cohorts. Matched tumor DNA sequencing was evaluated for second hits, and germline DNA array data from 5585 neuroblastoma patients and 23 505 cancer-free control children were analyzed to identify rare germline copy number variants. Patients with germline pathogenic or likely pathogenic variants were compared with those without to test for association with clinical characteristics, tumor features, and survival. RESULTS: We observed 116 pathogenic or likely pathogenic variants involving 13.9% (109 of 786) of neuroblastoma patients, representing a statistically significant excess burden compared with cancer-free participants (odds ratio [OR] = 1.60, 95% confidence interval [CI] = 1.27 to 2.00). BARD1 harbored the most statistically significant enrichment of pathogenic or likely pathogenic variants (OR = 32.30, 95% CI = 6.44 to 310.35). Rare germline copy number variants disrupting BARD1 were identified in patients but absent in cancer-free participants (OR = 29.47, 95% CI = 1.52 to 570.70). Patients harboring a germline pathogenic or likely pathogenic variant had a worse overall survival compared with those without (P = 8.6 x 10-3). CONCLUSIONS: BARD1 is an important neuroblastoma predisposition gene harboring both common and rare germline pathogenic or likely pathogenic variations. The presence of any germline pathogenic or likely pathogenic variant in a cancer predisposition gene was independently predictive of worse overall survival. As centers move toward paired tumor-normal sequencing at diagnosis, efforts should be made to centralize data and provide an infrastructure to support cooperative longitudinal prospective studies of germline pathogenic variation.

The Role of FAS Receptor Methylation in Osteosarcoma Metastasis.

Osteosarcoma is the most frequent primary malignant bone tumor with an annual incidence of about 400 cases in the United States. Osteosarcoma primarily metastasizes to the lungs, where FAS ligand (FASL) is constitutively expressed. The interaction of FASL and its cell surface receptor, FAS, triggers apoptosis in normal cells; however, this function is altered in cancer cells. DNA methylation has previously been explored as a mechanism for altering FAS expression, but no variability was identified in the CpG island (CGI) overlapping the promoter. Analysis of an expanded region, including CGI shores and shelves, revealed high variability in the methylation of certain CpG sites that correlated significantly with FAS mRNA expression in a negative manner. Bisulfite sequencing revealed additional CpG sites, which were highly methylated in the metastatic LM7 cell line but unmethylated in its parental non-metastatic SaOS-2 cell line. Treatment with the demethylating agent, 5-azacytidine, resulted in a loss of methylation in CpG sites located within the FAS promoter and restored FAS protein expression in LM7 cells, resulting in reduced migration. Orthotopic implantation of 5-azacytidine treated LM7 cells into severe combined immunodeficient mice led to decreased lung metastases. These results suggest that DNA methylation of CGI shore sites may regulate FAS expression and constitute a potential target for osteosarcoma therapy, utilizing demethylating agents currently approved for the treatment of other cancers.

Germline pathogenic variants in 786 neuroblastoma patients.

Importance: Neuroblastoma accounts for 12% of childhood cancer deaths. The genetic contribution of rare pathogenic germline variation in patients without a family history remains unclear. Objective: To define the prevalence, spectrum, and clinical significance of pathogenic germline variation in cancer predisposition genes (CPGs) in neuroblastoma patients. Design Setting and Participants: Germline DNA sequencing was performed on the peripheral blood from 786 neuroblastoma patients unselected for family history. Rare variants mapping to CPGs were evaluated for pathogenicity and the percentage of cases harboring pathogenic (P) or likely pathogenic (LP) variants was quantified. The frequency of CPG P-LP variants in neuroblastoma cases was compared to two distinct cancer-free control cohorts to assess enrichment. Matched tumor DNA sequencing was evaluated for "second hits" at CPGs and germline DNA array data from 5,585 neuroblastoma cases and 23,505 cancer-free control children was analyzed to identify rare germline copy number variants (CNVs) affecting genes with an excess burden of P-LP variants in neuroblastoma. Neuroblastoma patients with germline P-LP variants were compared to those without P-LP variants to test for association with clinical characteristics, tumor features, and patient survival. Main Outcomes and Measures: Rare variant prevalence, pathogenicity, enrichment, and association with clinical characteristics, tumor features, and patient survival. Results: We observed 116 P-LP variants in CPGs involving 13.9% (109/786) of patients, representing a significant excess burden of P-LP variants compared to controls (9.1%; P = 5.14 × 10-5, Odds Ratio: 1.60, 95% confidence interval: 1.27-2.00). BARD1 harbored the most significant burden of P-LP variants compared to controls (1.0% vs. 0.03%; P = 8.18 × 10-7; Odds Ratio: 32.30, 95% confidence interval: 6.44-310.35). Rare germline CNVs disrupting BARD1 were also identified in neuroblastoma patients (0.05%) but absent in controls (P = 7.08 × 10-3; Odds Ratio: 29.47, 95% confidence interval: 1.52 - 570.70). Overall, P-LP variants in DNA repair genes in this study were enriched in cases compared to controls (8.1% vs. 5.7%; P = 0.01; Odds Ratio: 1.45, 95% confidence interval: 1.08-1.92). Neuroblastoma patients harboring a germline P-LP variant had a worse overall survival when compared to patients without P-LP variants (P = 8.6 × 10-3), and this remained significant in a multivariate Cox proportional-hazards model (P = 0.01). Conclusions and Relevance: Neuroblastoma patients harboring germline P-LP variants in CPGs have worse overall survival and BARD1 is an important predisposition gene affected by both common and rare pathogenic variation. Germline sequencing should be performed for all neuroblastoma patients at diagnosis to inform genetic counseling and support future longitudinal and mechanistic studies. Patients with a germline P-LP variant should be closely monitored, regardless of risk group assignment.

Genetic subgroups inform on pathobiology in adult and pediatric Burkitt lymphoma.

Burkitt lymphoma (BL) accounts for most pediatric non-Hodgkin lymphomas, being less common but significantly more lethal when diagnosed in adults. Much of the knowledge of the genetics of BL thus far has originated from the study of pediatric BL (pBL), leaving its relationship to adult BL (aBL) and other adult lymphomas not fully explored. We sought to more thoroughly identify the somatic changes that underlie lymphomagenesis in aBL and any molecular features that associate with clinical disparities within and between pBL and aBL. Through comprehensive whole-genome sequencing of 230 BL and 295 diffuse large B-cell lymphoma (DLBCL) tumors, we identified additional significantly mutated genes, including more genetic features that associate with tumor Epstein-Barr virus status, and unraveled new distinct subgroupings within BL and DLBCL with 3 predominantly comprising BLs: DGG-BL (DDX3X, GNA13, and GNAI2), IC-BL (ID3 and CCND3), and Q53-BL (quiet TP53). Each BL subgroup is characterized by combinations of common driver and noncoding mutations caused by aberrant somatic hypermutation. The largest subgroups of BL cases, IC-BL and DGG-BL, are further characterized by distinct biological and gene expression differences. IC-BL and DGG-BL and their prototypical genetic features (ID3 and TP53) had significant associations with patient outcomes that were different among aBL and pBL cohorts. These findings highlight shared pathogenesis between aBL and pBL, and establish genetic subtypes within BL that serve to delineate tumors with distinct molecular features, providing a new framework for epidemiologic, diagnostic, and therapeutic strategies.

A community challenge for a pancancer drug mechanism of action inference from perturbational profile data.

The Columbia Cancer Target Discovery and Development (CTD2) Center is developing PANACEA, a resource comprising dose-responses and RNA sequencing (RNA-seq) profiles of 25 cell lines perturbed with ∼400 clinical oncology drugs, to study a tumor-specific drug mechanism of action. Here, this resource serves as the basis for a DREAM Challenge assessing the accuracy and sensitivity of computational algorithms for de novo drug polypharmacology predictions. Dose-response and perturbational profiles for 32 kinase inhibitors are provided to 21 teams who are blind to the identity of the compounds. The teams are asked to predict high-affinity binding targets of each compound among ∼1,300 targets cataloged in DrugBank. The best performing methods leverage gene expression profile similarity analysis as well as deep-learning methodologies trained on individual datasets. This study lays the foundation for future integrative analyses of pharmacogenomic data, reconciliation of polypharmacology effects in different tumor contexts, and insights into network-based assessments of drug mechanisms of action.

An expanded universe of cancer targets.

The characterization of cancer genomes has provided insight into somatically altered genes across tumors, transformed our understanding of cancer biology, and enabled tailoring of therapeutic strategies. However, the function of most cancer alleles remains mysterious, and many cancer features transcend their genomes. Consequently, tumor genomic characterization does not influence therapy for most patients. Approaches to understand the function and circuitry of cancer genes provide complementary approaches to elucidate both oncogene and non-oncogene dependencies. Emerging work indicates that the diversity of therapeutic targets engendered by non-oncogene dependencies is much larger than the list of recurrently mutated genes. Here we describe a framework for this expanded list of cancer targets, providing novel opportunities for clinical translation.

Framework for quality assessment of whole genome cancer sequences.

Bringing together cancer genomes from different projects increases power and allows the investigation of pan-cancer, molecular mechanisms. However, working with whole genomes sequenced over several years in different sequencing centres requires a framework to compare the quality of these sequences. We used the Pan-Cancer Analysis of Whole Genomes cohort as a test case to construct such a framework. This cohort contains whole cancer genomes of 2832 donors from 18 sequencing centres. We developed a non-redundant set of five quality control (QC) measurements to establish a star rating system. These QC measures reflect known differences in sequencing protocol and provide a guide to downstream analyses and allow for exclusion of samples of poor quality. We have found that this is an effective framework of quality measures. The implementation of the framework is available at: https://dockstore.org/containers/quay.io/jwerner_dkfz/pancanqc:1.2.2 .

Somatic structural variation targets neurodevelopmental genes and identifies SHANK2 as a tumor suppressor in neuroblastoma.

Neuroblastoma is a malignancy of the developing sympathetic nervous system that accounts for 12% of childhood cancer deaths. Like many childhood cancers, neuroblastoma shows a relative paucity of somatic single-nucleotide variants (SNVs) and small insertions and deletions (indels) compared to adult cancers. Here, we assessed the contribution of somatic structural variation (SV) in neuroblastoma using a combination of whole-genome sequencing (WGS) of tumor-normal pairs (n = 135) and single-nucleotide polymorphism (SNP) genotyping of primary tumors (n = 914). Our study design allowed for orthogonal validation and replication across platforms. SV frequency, type, and localization varied significantly among high-risk tumors. MYCN nonamplified high-risk tumors harbored an increased SV burden overall, including a significant excess of tandem duplication events across the genome. Genes disrupted by SV breakpoints were enriched in neuronal lineages and associated with phenotypes such as autism spectrum disorder (ASD). The postsynaptic adapter protein-coding gene, SHANK2, located on Chromosome 11q13, was disrupted by SVs in 14% of MYCN nonamplified high-risk tumors based on WGS and 10% in the SNP array cohort. Expression of SHANK2 was low across human-derived neuroblastoma cell lines and high-risk neuroblastoma tumors. Forced expression of SHANK2 in neuroblastoma cells resulted in significant growth inhibition (P = 2.6 × 10-2 to 3.4 × 10-5) and accelerated neuronal differentiation following treatment with all-trans retinoic acid (P = 3.1 × 10-13 to 2.4 × 10-30). These data further define the complex landscape of somatic structural variation in neuroblastoma and suggest that events leading to deregulation of neurodevelopmental processes, such as inactivation of SHANK2, are key mediators of tumorigenesis in this childhood cancer.

Analysis of Ugandan cervical carcinomas identifies human papillomavirus clade-specific epigenome and transcriptome landscapes.

Cervical cancer is the most common cancer affecting sub-Saharan African women and is prevalent among HIV-positive (HIV+) individuals. No comprehensive profiling of cancer genomes, transcriptomes or epigenomes has been performed in this population thus far. We characterized 118 tumors from Ugandan patients, of whom 72 were HIV+, and performed extended mutation analysis on an additional 89 tumors. We detected human papillomavirus (HPV)-clade-specific differences in tumor DNA methylation, promoter- and enhancer-associated histone marks, gene expression and pathway dysregulation. Changes in histone modification at HPV integration events were correlated with upregulation of nearby genes and endogenous retroviruses.

Identification and Analyses of Extra-Cranial and Cranial Rhabdoid Tumor Molecular Subgroups Reveal Tumors with Cytotoxic T Cell Infiltration.

Extra-cranial malignant rhabdoid tumors (MRTs) and cranial atypical teratoid RTs (ATRTs) are heterogeneous pediatric cancers driven primarily by SMARCB1 loss. To understand the genome-wide molecular relationships between MRTs and ATRTs, we analyze multi-omics data from 140 MRTs and 161 ATRTs. We detect similarities between the MYC subgroup of ATRTs (ATRT-MYC) and extra-cranial MRTs, including global DNA hypomethylation and overexpression of HOX genes and genes involved in mesenchymal development, distinguishing them from other ATRT subgroups that express neural-like features. We identify five DNA methylation subgroups associated with anatomical sites and SMARCB1 mutation patterns. Groups 1, 3, and 4 exhibit cytotoxic T cell infiltration and expression of immune checkpoint regulators, consistent with a potential role for immunotherapy in rhabdoid tumor patients.

Before and After: Comparison of Legacy and Harmonized TCGA Genomic Data Commons' Data.

We present a systematic analysis of the effects of synchronizing a large-scale, deeply characterized, multi-omic dataset to the current human reference genome, using updated software, pipelines, and annotations. For each of 5 molecular data platforms in The Cancer Genome Atlas (TCGA)-mRNA and miRNA expression, single nucleotide variants, DNA methylation and copy number alterations-comprehensive sample, gene, and probe-level studies were performed, towards quantifying the degree of similarity between the 'legacy' GRCh37 (hg19) TCGA data and its GRCh38 (hg38) version as 'harmonized' by the Genomic Data Commons. We offer gene lists to elucidate differences that remained after controlling for confounders, and strategies to mitigate their impact on biological interpretation. Our results demonstrate that the hg19 and hg38 TCGA datasets are very highly concordant, promote informed use of either legacy or harmonized omics data, and provide a rubric that encourages similar comparisons as new data emerge and reference data evolve.

Genetic mechanisms of primary chemotherapy resistance in pediatric acute myeloid leukemia.

Acute myeloid leukemias (AML) are characterized by mutations of tumor suppressor and oncogenes, involving distinct genes in adults and children. While certain mutations have been associated with the increased risk of AML relapse, the genomic landscape of primary chemotherapy-resistant AML is not well defined. As part of the TARGET initiative, we performed whole-genome DNA and transcriptome RNA and miRNA sequencing analysis of pediatric AML with failure of induction chemotherapy. We identified at least three genetic groups of patients with induction failure, including those with NUP98 rearrangements, somatic mutations of WT1 in the absence of apparent NUP98 mutations, and additional recurrent variants including those in KMT2C and MLLT10. Comparison of specimens before and after chemotherapy revealed distinct and invariant gene expression programs. While exhibiting overt therapy resistance, these leukemias nonetheless showed diverse forms of clonal evolution upon chemotherapy exposure. This included selection for mutant alleles of FRMD8, DHX32, PIK3R1, SHANK3, MKLN1, as well as persistence of WT1 and TP53 mutant clones, and elimination of FLT3, PTPN11, and NRAS mutant clones. These findings delineate genetic mechanisms of primary chemotherapy resistance in pediatric AML, which should inform improved approaches for its diagnosis and therapy.

Publisher Correction: The molecular landscape of pediatric acute myeloid leukemia reveals recurrent structural alterations and age-specific mutational interactions.

In the version of this article originally published, the color key in Fig. 1a was wrong. In the Cytogenetics key, the box over t(8;21) originally was green. It should have been red, matching the color of the sections of the pie graphs below the key that were labeled with 15% and 19%.

Genome-wide discovery of somatic coding and noncoding mutations in pediatric endemic and sporadic Burkitt lymphoma.

Although generally curable with intensive chemotherapy in resource-rich settings, Burkitt lymphoma (BL) remains a deadly disease in older patients and in sub-Saharan Africa. Epstein-Barr virus (EBV) positivity is a feature in more than 90% of cases in malaria-endemic regions, and up to 30% elsewhere. However, the molecular features of BL have not been comprehensively evaluated when taking into account tumor EBV status or geographic origin. Through an integrative analysis of whole-genome and transcriptome data, we show a striking genome-wide increase in aberrant somatic hypermutation in EBV-positive tumors, supporting a link between EBV and activation-induced cytidine deaminase (AICDA) activity. In addition to identifying novel candidate BL genes such as SIN3A, USP7, and CHD8, we demonstrate that EBV-positive tumors had significantly fewer driver mutations, especially among genes with roles in apoptosis. We also found immunoglobulin variable region genes that were disproportionally used to encode clonal B-cell receptors (BCRs) in the tumors. These include IGHV4-34, known to produce autoreactive antibodies, and IGKV3-20, a feature described in other B-cell malignancies but not yet in BL. Our results suggest that tumor EBV status defines a specific BL phenotype irrespective of geographic origin, with particular molecular properties and distinct pathogenic mechanisms. The novel mutation patterns identified here imply rational use of DNA-damaging chemotherapy in some patients with BL and targeted agents such as the CDK4/6 inhibitor palbociclib in others, whereas the importance of BCR signaling in BL strengthens the potential benefit of inhibitors for PI3K, Syk, and Src family kinases among these patients.

A unique subset of low-risk Wilms tumors is characterized by loss of function of TRIM28 (KAP1), a gene critical in early renal development: A Children's Oncology Group study.

This study explores the genomic alterations that contribute to the formation of a unique subset of low-risk, epithelial differentiated, favorable histology Wilms tumors (WT), tumors that have been characterized by their expression of post-induction renal developmental genes (Subset 1 WT). We demonstrate copy neutral loss of heterozygosity involving 19q13.32-q13.43, unaccompanied by evidence for imprinting by DNA methylation. We further identified loss-of-function somatic mutations in TRIM28 (also known as KAP1), located at 19q13, in 8/9 Subset 1 tumors analyzed. An additional germline TRIM28 mutation was identified in one patient. Retrospective evaluation of previously analyzed WT outside of Subset 1 identified an additional tumor with anaplasia and both TRIM28 and TP53 mutations. A major function of TRIM28 is the repression of endogenous retroviruses early in development. We depleted TRIM28 in HEK293 cells, which resulted in increased expression of endogenous retroviruses, a finding also demonstrated in TRIM28-mutant WT. TRIM28 has been shown by others to be active during early renal development, and to interact with WTX, another gene recurrently mutated in WT. Our findings suggest that inactivation of TRIM28 early in renal development contributes to the formation of this unique subset of FHWTs, although the precise manner in which TRIM28 impacts both normal renal development and oncogenesis remains elusive.

Clinically Relevant Cytotoxic Immune Cell Signatures and Clonal Expansion of T-Cell Receptors in High-Risk MYCN-Not-Amplified Human Neuroblastoma.

Purpose: High-risk neuroblastoma is an aggressive disease. DNA sequencing studies have revealed a paucity of actionable genomic alterations and a low mutation burden, posing challenges to develop effective novel therapies. We used RNA sequencing (RNA-seq) to investigate the biology of this disease, including a focus on tumor-infiltrating lymphocytes (TIL).Experimental Design: We performed deep RNA-seq on pretreatment diagnostic tumors from 129 high-risk and 21 low- or intermediate-risk patients with neuroblastomas. We used single-sample gene set enrichment analysis to detect gene expression signatures of TILs in tumors and examined their association with clinical and molecular parameters, including patient outcome. The expression profiles of 190 additional pretreatment diagnostic neuroblastomas, a neuroblastoma tissue microarray, and T-cell receptor (TCR) sequencing were used to validate our findings.Results: We found that MYCN-not-amplified (MYCN-NA) tumors had significantly higher cytotoxic TIL signatures compared with MYCN-amplified (MYCN-A) tumors. A reported MYCN activation signature was significantly associated with poor outcome for high-risk patients with MYCN-NA tumors; however, a subgroup of these patients who had elevated activated natural killer (NK) cells, CD8+ T cells, and cytolytic signatures showed improved outcome and expansion of infiltrating TCR clones. Furthermore, we observed upregulation of immune exhaustion marker genes, indicating an immune-suppressive microenvironment in these neuroblastomas.Conclusions: This study provides evidence that RNA signatures of cytotoxic TIL are associated with the presence of activated NK/T cells and improved outcomes in high-risk neuroblastoma patients harboring MYCN-NA tumors. Our findings suggest that these high-risk patients with MYCN-NA neuroblastoma may benefit from additional immunotherapies incorporated into the current therapeutic strategies. Clin Cancer Res; 24(22); 5673-84. ©2018 AACR.

Erratum: The molecular landscape of pediatric acute myeloid leukemia reveals recurrent structural alterations and age-specific mutational interactions.

This corrects the article DOI: 10.1038/nm.4439.

Pan-cancer genome and transcriptome analyses of 1,699 paediatric leukaemias and solid tumours.

Analysis of molecular aberrations across multiple cancer types, known as pan-cancer analysis, identifies commonalities and differences in key biological processes that are dysregulated in cancer cells from diverse lineages. Pan-cancer analyses have been performed for adult but not paediatric cancers, which commonly occur in developing mesodermic rather than adult epithelial tissues. Here we present a pan-cancer study of somatic alterations, including single nucleotide variants, small insertions or deletions, structural variations, copy number alterations, gene fusions and internal tandem duplications in 1,699 paediatric leukaemias and solid tumours across six histotypes, with whole-genome, whole-exome and transcriptome sequencing data processed under a uniform analytical framework. We report 142 driver genes in paediatric cancers, of which only 45% match those found in adult pan-cancer studies; copy number alterations and structural variants constituted the majority (62%) of events. Eleven genome-wide mutational signatures were identified, including one attributed to ultraviolet-light exposure in eight aneuploid leukaemias. Transcription of the mutant allele was detectable for 34% of protein-coding mutations, and 20% exhibited allele-specific expression. These data provide a comprehensive genomic architecture for paediatric cancers and emphasize the need for paediatric cancer-specific development of precision therapies.