Genomic landscape of Down syndrome-associated acute lymphoblastic leukemia

Trisomy 21, the genetic cause of Down syndrome (DS), is the most common congenital chromosomal anomaly. It is associated with a 20-fold increased risk of acute lymphoblastic leukemia (ALL) during childhood and results in distinctive leukemia biology. To comprehensively define the genomic landscape of DS-ALL, we performed whole-genome sequencing and whole-transcriptome sequencing (RNA-Seq) on 295 cases. Our integrated genomic analyses identified 15 molecular subtypes of DS-ALL, with marked enrichment of CRLF2-r, IGH::IGF2BP1, and C/EBP altered (C/EBPalt) subtypes compared with 2257 non-DS-ALL cases. We observed abnormal activation of the CEBPD, CEBPA, and CEBPE genes in 10.5% of DS-ALL cases via a variety of genomic mechanisms, including chromosomal rearrangements and noncoding mutations leading to enhancer hijacking. A total of 42.3% of C/EBP-activated DS-ALL also have concomitant FLT3 point mutations or insertions/deletions, compared with 4.1% in other subtypes. CEBPD overexpression enhanced the differentiation of mouse hematopoietic progenitor cells into pro-B cells in vitro, particularly in a DS genetic background. Notably, recombination-activating gene-mediated somatic genomic abnormalities were common in DS-ALL, accounting for a median of 27.5% of structural alterations, compared with 7.7% in non-DS-ALL. Unsupervised hierarchical clustering analyses of CRLF2-rearranged DS-ALL identified substantial heterogeneity within this group, with the BCR::ABL1-like subset linked to an inferior event-free survival, even after adjusting for known clinical risk factors. These results provide important insights into the biology of DS-ALL and point to opportunities for targeted therapy and treatment individualization.

Racial and ethnic disparities in childhood and young adult acute lymphocytic leukaemia: secondary analyses of eight Children's Oncology Group cohort trials

BACKGROUND

Previous studies have identified racial and ethnic disparities in childhood acute lymphocytic leukaemia survival. We aimed to establish whether disparities persist in contemporaneous cohorts and, if present, are attributable to differences in leukaemia biology or insurance status.

METHODS

Patients with newly diagnosed acute lymphocytic leukaemia in inpatient and outpatient centres in the USA, Canada, Australia, and New Zealand, aged 0-30 years, who had race or ethnicity data available, enrolled on eight completed Children's Oncology Group trials (NCT00103285, NCT00075725, NCT00408005, NCT01190930, NCT02883049, NCT02112916, NCT02828358, and NCT00557193) were included in this secondary analysis. Race and ethnicity were categorised as non-Hispanic White, Hispanic, non-Hispanic Black, non-Hispanic Asian, and non-Hispanic other. Event-free survival and overall survival were compared across race and ethnicity groups. The relative contribution of clinical and biological disease prognosticators and insurance status was examined through multivariable regression models, both among the entire cohort and among those with B-cell lineage versus T-cell lineage disease.

FINDINGS

Between Jan 1, 2004, and Dec 31, 2019, 24 979 eligible children, adolescents, and young adults with acute lymphocytic leukaemia were enrolled, of which 21 152 had race or ethnicity data available. 11 849 (56·0%) were male and 9303 (44·0%) were female. Non-Hispanic White patients comprised the largest racial or ethnic group (13 872 [65·6%]), followed by Hispanic patients (4354 [20·6%]), non-Hispanic Black patients (1517 [7·2%]), non-Hispanic Asian (n=1071 [5·1%]), and non-Hispanic other (n=338 [1·6%]). 5-year event-free survival was 87·4% (95% CI 86·7-88·0%) among non-Hispanic White patients compared with 82·8% (81·4-84·1%; hazard ratio [HR] 1·37, 95% CI 1·26-1·49; p<0·0001) among Hispanic patients and 81·8% (79·3-84·0; HR 1·45, 1·28-1·65; p<0·0001) among non-Hispanic Black patients. Non-hispanic Asian patients had a 5-year event-free survival of 88·1% (95% CI 85·5-90·3%) and non-Hispanic other patients had a survival of 82·8% (76·4-87·6%). Inferior event-free survival among Hispanic patients was substantially attenuated by disease prognosticators and insurance status (HR decreased from 1·37 [1·26-1·49; p<0·0001] to 1·11 [1·00-1·22; p=0·045]). The increased risk among non-Hispanic Black patients was minimally attenuated (HR 1·45 [1·28-1·65; p<0·0001] to 1·32 [1·14-1·52; p<0·0001]). 5-year overall survival was 93·6% (91·5-95·1%) in non-Hispanic Asian patients, 93·3% (92·8-93·7%) in non-Hispanic White patients, 89·9% (88·7-90·9%) in Hispanic, 89·7% (87·6-91·4%) in non-Hispanic Black patients, 88·9% (83·2-92·7%) in non-Hispanic other patients. Disparities in overall survival were wider than event-free survival (eg, among non-Hispanic other patients, the HR for event-free survival was 1·43 [1·10-1·85] compared with 1·74 [1·27-2·40] for overall survival). Disparities were restricted to patients with B-cell acute lymphocytic leukaemia, no differences in event-free survival or overall survival were seen in the T-cell acute lymphocytic leukaemia group.

INTERPRETATION

Substantial disparities in outcome for B-cell acute lymphocytic leukaemia persist by race and ethnicity, but are not observed in T-cell acute lymphocytic leukaemia. Future studies of relapsed patients, access to and quality of care, and other potential aspects of structural racism are warranted to inform interventions aimed at dismantling racial and ethnic disparities.

FUNDING

National Cancer Institute and St Baldrick's Foundation.

Multiomics in primary and metastatic breast tumors from the AURORA US network finds microenvironment and epigenetic drivers of metastasis

The AURORA US Metastasis Project was established with the goal to identify molecular features associated with metastasis. We assayed 55 females with metastatic breast cancer (51 primary cancers and 102 metastases) by RNA sequencing, tumor/germline DNA exome and low-pass whole-genome sequencing and global DNA methylation microarrays. Expression subtype changes were observed in ~30% of samples and were coincident with DNA clonality shifts, especially involving HER2. Downregulation of estrogen receptor (ER)-mediated cell-cell adhesion genes through DNA methylation mechanisms was observed in metastases. Microenvironment differences varied according to tumor subtype; the ER+/luminal subtype had lower fibroblast and endothelial content, while triple-negative breast cancer/basal metastases showed a decrease in B and T cells. In 17% of metastases, DNA hypermethylation and/or focal deletions were identified near HLA-A and were associated with reduced expression and lower immune cell infiltrates, especially in brain and liver metastases. These findings could have implications for treating individuals with metastatic breast cancer with immune- and HER2-targeting therapies.

The genomic landscape of pediatric acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is the most common childhood cancer. Here, using whole-genome, exome and transcriptome sequencing of 2,754 childhood patients with ALL, we find that, despite a generally low mutation burden, ALL cases harbor a median of four putative somatic driver alterations per sample, with 376 putative driver genes identified varying in prevalence across ALL subtypes. Most samples harbor at least one rare gene alteration, including 70 putative cancer driver genes associated with ubiquitination, SUMOylation, noncoding transcripts and other functions. In hyperdiploid B-ALL, chromosomal gains are acquired early and synchronously before ultraviolet-induced mutation. By contrast, ultraviolet-induced mutations precede chromosomal gains in B-ALL cases with intrachromosomal amplification of chromosome 21. We also demonstrate the prognostic significance of genetic alterations within subtypes. Intriguingly, DUX4- and KMT2A-rearranged subtypes separate into CEBPA/FLT3- or NFATC4-expressing subgroups with potential clinical implications. Together, these results deepen understanding of the ALL genomic landscape and associated outcomes.

Noncoding genetic variation in GATA3 increases acute lymphoblastic leukemia risk through local and global changes in chromatin conformation

Inherited noncoding genetic variants confer significant disease susceptibility to childhood acute lymphoblastic leukemia (ALL) but the molecular processes linking germline polymorphisms with somatic lesions in this cancer are poorly understood. Through targeted sequencing in 5,008 patients, we identified a key regulatory germline variant in GATA3 associated with Philadelphia chromosome-like ALL (Ph-like ALL). Using CRISPR-Cas9 editing and samples from patients with Ph-like ALL, we showed that this variant activated a strong enhancer that upregulated GATA3 transcription. This, in turn, reshaped global chromatin accessibility and three-dimensional genome organization, including regions proximal to the ALL oncogene CRLF2. Finally, we showed that GATA3 directly regulated CRLF2 and potentiated the JAK-STAT oncogenic effects during leukemogenesis. Taken together, we provide evidence for a distinct mechanism by which a germline noncoding variant contributes to oncogene activation, epigenetic regulation and three-dimensional genome reprogramming.

Abstract 98: Custom Gene Fusion Assays for the Rapid Diagnosis of Pediatric Cancers in Low-Resourced Settings

Purpose: Risk stratification and molecular targeting have been key to increasing cure rates for pediatric cancers in high-income countries. Precise diagnosis and successful treatment of pediatric cancer in low-resourced settings is often hindered by insufficient pathology infrastructure, including lack of laboratory platforms for molecular analysis. Given the high frequency of gene fusions in pediatric cancers, identifying such fusions would greatly aid cost-effective pediatric cancer diagnosis, risk stratification, and precision medicine in low-resourced settings.

Methods: To allow for implementation of gene fusion detection at Global HOPE sites in Sub-Saharan Africa (SSA), methodologies were reviewed to consider minimal technical expertise required, the ability to utilize samples with sub-optimal RNA quality, and rapid turn-around-time. Literature review, clinical laboratory results, public databases, and large-scale genomic studies were used to obtain exact breakpoint sequence information for gene fusions associated with pediatric and adolescent cancers.

Results: Two custom pediatric gene fusions panels were designed using the NanoString Elements technology. The hematologic malignancy panel was designed to detect 439 breakpoints for 223 non-IGH/TCR fusions reported in ALL, AML, lymphomas, and histiocytosis. The solid tumor panel was designed to detect 204 breakpoints for 436 fusions associated with pediatric sarcomas, brain tumors, and renal malignancies. Each panel was tested using 96 samples with known fusion status at Texas Children's Hospital to determine specificity, sensitivity, precision, and ease of workflow.

Conclusions: The design, testing, and implementation of a rapid assay to detect gene fusions with diagnostic, prognostic, and therapeutic impact would be transformational in the care of pediatric cancer patients in low-resourced settings. The custom designed panels will allow for large-scale fusion detection in 2-3 days with only 15 minutes of technician time after RNA isolation. Additional steps are needed to identify and address any challenges upon initiation in SSA to fully realize the potential of such technology.

Citation Format: Julie Gastier-Foster, Fredrick Lutwama, Joseph Lubega, Kevin Fisher, Dolores Lopez-Terrada, Angshumoy Roy, Nmazuo Ozuah, Jeremy Slone, Peter Wasswa, Carl Allen, David Poplack. Custom Gene Fusion Assays for the Rapid Diagnosis of Pediatric Cancers in Low-Resourced Settings [abstract]. In: Proceedings of the 9th Annual Symposium on Global Cancer Research; Global Cancer Research and Control: Looking Back and Charting a Path Forward; 2021 Mar 10-11. Philadelphia (PA): AACR; Cancer Epidemiol Biomarkers Prev 2021;30(7 Suppl):Abstract nr 98.

Abstract 2118: Non-coding germline GATA3 variants alter chromatin topology and contribute to pathogenesis of acute lymphoblastic leukemia

Inherited non-coding genetic variants confer significant disease susceptibility in many cancers. However, the molecular processes of by which germline variants contribute to somatic lesions are poorly understood. We performed targeted sequencing in 5,008 patients and identified a key regulatory germline variant in GATA3 strongly associated with Philadelphia chromosome-like acute lymphoblastic leukemia (Ph-like ALL). By creating an isogenic cellular model with CRISPR-Cas9 system, we showed that this variant activated a strong enhancer that significantly upregulated GATA3 transcription, which in turn reshaped the global chromatin accessibility and 3D genome organization. Remarkably, this genotype switch induced a chromatin loop between the CRLF2 oncogene and a distal enhancer, similar to the somatically acquired super-enhancer hijacking event in patients. GATA3 genotype-related alterations in transcriptional control and 3D chromatin organization were further validated in Ph-like ALL patients. Finally, we showed that GATA3 directly regulates CRLF2 and potentiates the oncogenic effects of JAK-STAT signaling in leukemogenesis. Altogether, our results provide evidence for a novel mechanism by which a germline non-coding variant contributes to oncogene activation epigenetic regulation and 3D genome reprogramming.

Citation Format: Hongbo Yang, Hui Zhang, Yu Luan, Tingting Liu, Kathryn Roberts, Mao-xiang Qian, Bo Zhang, Wenjian Yang, Virginia Perez-Andreu, Jie Xu, Sriranga lyyanki, Da Kuang, Shalini Reshmi, Julie Gastier-Foster, Colton Smith, Ching-Hon Pui, William Evans, Stephen Hunger, Stephen Hunger, Leonidas Platanias, Mary Relling, Charles Mullighan, Mignon Loh, Feng Yue, Jun Yang. Non-coding germline GATA3 variants alter chromatin topology and contribute to pathogenesis of acute lymphoblastic leukemia [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2118.

Abstract PO-055: Molecular characterization of papillary thyroid cancer in relation to ionizing radiation dose following the Chernobyl accident

The 1986 Chernobyl nuclear power plant accident increased papillary thyroid cancer (PTC) incidence in surrounding regions, particularly for 131I-exposed children. To investigate the contribution of environmental radiation to PTC characteristics and improve understanding of radiation-induced carcinogenesis, we analyzed genomic, transcriptomic, and epigenomic characteristics of 440 pathologically-confirmed fresh-frozen PTCs from Ukraine (359 with estimated childhood or in utero 131I exposure and 81 from unexposed children born after March 1987) and matched normal tissue (non-tumor thyroid tissue and/or blood). Mean age at PTC was 28.0 years (range: 10.0-45.6). Among 131I-exposed individuals, mean radiation dose was 250 mGy (range: 11.0-8,800). In multivariable models adjusted for age at PTC and sex, we observed radiation dose-dependent enrichment of fusion drivers (P=6.6 × 10−8), nearly all occurring in the MAPK pathway, as well as increases in small deletions (P=8.0 × 10−9) and simple/balanced structural variants (P=1.2 × 10−14). Further analyses demonstrated even stronger associations for those small deletions and simple/balanced structural variants that were clonal and bore hallmarks of non-homologous end-joining repair (deletions: P=4.9 × 10−31; simple/balanced structural variants: P=5.5 × 10−19). In contrast, radiation dose was not associated with subclonal small deletions (P=0.82) or subclonal simple/balanced structural variants (P=0.91). Additionally, radiation dose was not associated with TINS (locally templated insertions), which are characteristic of alt-end-joining repair (P=0.69). The effects of radiation on genomic alterations with more pronounced for those younger at exposure. Analyses generally were consistent with a linear radiation dose-response for all molecular characteristics except clonal small deletions. Analyses of transcriptomic and epigenomic features demonstrated strong associations with the PTC driver gene but not radiation dose. Our results point to DNA double-strand breaks as early carcinogenic events that subsequently enable PTC growth following environmental radiation exposure.

Citation Format: Lindsay M. Morton, Danielle Karyadi, Chip Stewart, Tetiana Bogdanova, Eric Dawson, Mia Steinberg, Jieqiong Dai, Stephen Hartley, Sara Schonfeld, Joshua Sampson, Yosi Maruvka, Vidushi Kapoor, Dale Ramsden, Juan Carvajal-Garcia, Chuck Perou, Joel Parker, Marko Krznaric, Meredith Yeager, Joseph Boland, Amy Hutchinson, Belynda Hicks, Casey Dagnall, Julie Gastier-Foster, Jay Bowen, Olivia Lee, Mitchell Machiela, Elizabeth Cahoon, Alina Brenner, Kiyohiko Mabuchi, Vladimir Drozdovitch, Sergii Masiuk, Mykola Chepurny, Liudmyla Yu Zurnadzhy, Maureen Hatch, Amy Berrington de Gonzalez, Gerry Thomas, Mykola Tronko, Gad Getz, Stephen Chanock. Molecular characterization of papillary thyroid cancer in relation to ionizing radiation dose following the Chernobyl accident [abstract]. In: Proceedings of the AACR Virtual Special Conference on Radiation Science and Medicine; 2021 Mar 2-3. Philadelphia (PA): AACR; Clin Cancer Res 2021;27(8_Suppl):Abstract nr PO-055.

Molecular Features of Cancers Exhibiting Exceptional Responses to Treatment

A small fraction of cancer patients with advanced disease survive significantly longer than patients with clinically comparable tumors. Molecular mechanisms for exceptional responses to therapy have been identified by genomic analysis of tumor biopsies from individual patients. Here, we analyzed tumor biopsies from an unbiased cohort of 111 exceptional responder patients using multiple platforms to profile genetic and epigenetic aberrations as well as the tumor microenvironment. Integrative analysis uncovered plausible mechanisms for the therapeutic response in nearly a quarter of the patients. The mechanisms were assigned to four broad categories-DNA damage response, intracellular signaling, immune engagement, and genetic alterations characteristic of favorable prognosis-with many tumors falling into multiple categories. These analyses revealed synthetic lethal relationships that may be exploited therapeutically and rare genetic lesions that favor therapeutic success, while also providing a wealth of testable hypotheses regarding oncogenic mechanisms that may influence the response to cancer therapy.

The Exceptional Responders Initiative: Feasibility of a National Cancer Institute Pilot Study

BACKGROUND

Tumor molecular profiling from patients experiencing exceptional responses to systemic therapy may provide insights into cancer biology and improve treatment tailoring. This pilot study evaluates the feasibility of identifying exceptional responders retrospectively, obtaining pre-exceptional response treatment tumor tissues, and analyzing them with state-of-the-art molecular analysis tools to identify potential molecular explanations for responses.

METHODS

Exceptional response was defined as partial (PR) or complete (CR) response to a systemic treatment with population PR or CR rate less than 10% or an unusually long response (eg, duration >3 times published median). Cases proposed by patients' clinicians were reviewed by clinical and translational experts. Tumor and normal tissue (if possible) were profiled with whole exome sequencing and, if possible, targeted deep sequencing, RNA sequencing, methylation arrays, and immunohistochemistry. Potential germline mutations were tracked for relevance to disease.

RESULTS

Cases reflected a variety of tumors and standard and investigational treatments. Of 520 cases, 476 (91.5%) were accepted for further review, and 222 of 476 (46.6%) proposed cases met requirements as exceptional responders. Clinical data were obtained from 168 of 222 cases (75.7%). Tumor was provided from 130 of 168 cases (77.4%). Of 117 of the 130 (90.0%) cases with sufficient nucleic acids, 109 (93.2%) were successfully analyzed; 6 patients had potentially actionable germline mutations.

CONCLUSION

Exceptional responses occur with standard and investigational treatment. Retrospective identification of exceptional responders, accessioning, and sequencing of pretreatment archived tissue is feasible. Data from molecular analyses of tumors, particularly when combining results from patients who received similar treatments, may elucidate molecular bases for exceptional responses.

Transcriptome analysis of desmoplastic small round cell tumors identifies actionable therapeutic targets: a report from the Children's Oncology Group

To further understand the molecular pathogenesis of desmoplastic small round cell tumor (DSRCT), a fatal malignancy occurring primarily in adolescent/young adult males, we used next-generation RNA sequencing to investigate the gene expression profiles intrinsic to this disease. RNA from DSRCT specimens obtained from the Children's Oncology Group was sequenced using the Illumina HiSeq 2000 system and subjected to bioinformatic analyses. Validation and functional studies included WT1 ChIP-seq, EWS-WT1 knockdown using JN-DSRCT-1 cells and immunohistochemistry. A panel of immune signature genes was also evaluated to identify possible immune therapeutic targets. Twelve of 14 tumor samples demonstrated presence of the diagnostic EWSR1-WT1 translocation and these 12 samples were used for the remainder of the analysis. RNA sequencing confirmed the lack of full-length WT1 in all fusion positive samples as well as the JN-DSRCT-1 cell line. ChIP-seq for WT1 showed significant overlap with genes found to be highly expressed, including IGF2 and FGFR4, which were both highly expressed and targets of the EWS-WT1 fusion protein. In addition, we identified CD200 and CD276 as potentially targetable immune checkpoints whose expression is independent of the EWS-WT1 fusion gene in cultured DSCRT cells. In conclusion, we identified IGF2, FGFR4, CD200, and CD276 as potential therapeutic targets with clinical relevance for patients with DSRCT.

Strategies to Improve Pediatric Cancer Pathology and Tissue Handling for Clinical Care and Research in Low-Resource Settings

PURPOSE

Evidence-based pediatric cancer control is dependent on accurate, timely cancer diagnosis, as well as maintenance of tissue integrity and annotation for clinical diagnosis and clinical/translational research. Scarcity of cancer pathology infrastructure and the delinkage of clinical pathology and research workflows in low- and middle-income countries thwarts real-time, high-resolution diagnosis, risk stratification, and maintenance of tissue integrity for future research. This renders the application of evidence-based treatments and conduct of translational research and clinical trials, particularly across multiple centers in low- and middle-income countries, impractical. We present the innovative approach and impact of a Texas Children’s Global HOPE model designed to overcome these systemic obstacles in a sub-Saharan Africa setting.

METHODS

We designed a pathology algorithm at Mulago Hospital Complex in Uganda with the following goals: integrate clinical and research workflows, ensure safe and optimal tissue biopsies, and maintain tissue integrity for clinical and research assays. We introduced new assays built on locally available technology platforms for anatomic pathology and immunophenotyping.

RESULTS

The novel workflow and assays enabled consistent real-time hematoxylin and eosin microscopy, immunohistochemistry, and multiparametric flow cytometry. The proportion of cancer diagnoses based on tissue pathology increased from 67.5% to 100%. Streamlined pathology processes improved the timeliness of commencement of treatment from 4 to 6 weeks to 24 hours in the case of acute leukemias. The refined diagnosis process has enabled implementation of risk-adapted treatment regimens. Overall, these strategies contributed to improved overall survival from 30% to 85% at 1 month, and have revealed hitherto unrecognized clinical patterns of pediatric cancers in sub-Saharan Africa.

CONCLUSION

Accurate real-time pathology for pediatric cancers improves survival. Improvements should aim at enhancing the capabilities of existing in-country platforms and skilled personnel. Future efforts will focus on scaling up these practices to other sites and integrating molecular diagnostics that reduce overall cost and improve diagnostic accuracy or risk stratification.

Variant Interpretation for Dilated Cardiomyopathy: Refinement of the American College of Medical Genetics and Genomics/ClinGen Guidelines for the DCM Precision Medicine Study

BACKGROUND

The hypothesis of the Dilated Cardiomyopathy Precision Medicine Study is that most dilated cardiomyopathy has a genetic basis. The study returns results to probands and, when indicated, to relatives. While both the American College of Medical Genetics and Genomics/Association for Molecular Pathology and ClinGen's MYH7-cardiomyopathy specifications provide relevant guidance for variant interpretation, further gene- and disease-specific considerations were required for dilated cardiomyopathy. To this end, we tailored the ClinGen MYH7-cardiomyopathy variant interpretation framework; the specifications implemented for the study are presented here.

METHODS

Modifications were created and approved by an external Variant Adjudication Oversight Committee. After a pilot using 81 probands, further adjustments were made, resulting in 27 criteria (9 modifications of the ClinGen MYH7 framework and reintroduction of 2 American College of Medical Genetics and Genomics/Association of Molecular Pathology criteria that were deemed not applicable by the ClinGen MYH7 working group).

RESULTS

These criteria were applied to 2059 variants in a test set of 97 probands. Variants were classified as benign (n=1702), likely benign (n=33), uncertain significance (n=71), likely pathogenic (likely pathogenic; n=12), and pathogenic (P; n=3). Only 2/15 likely pathogenic/P variants were identified in Non-Hispanic African ancestry probands.

CONCLUSIONS

We tailored the ClinGen MYH7 criteria for our study. Our preliminary data show that 15/97 (15.5%) probands have likely pathogenic/P variants, most of which were identified in probands of Non-Hispanic European ancestry. We anticipate continued evolution of our approach, one that will be informed by new insights on variant interpretation and a greater understanding of the genetic architecture of dilated cardiomyopathy.

CLINICAL TRIAL REGISTRATION

URL: https://www.clinicaltrials.gov; Unique identifier: NCT03037632.

Abstract GS3-08: Multiplatform analysis of matched primary and metastatic breast tumors from the AURORA US Network

Background: It has become increasingly clear that effective treatment of metastatic breast cancer (MBC) requires an in-depth understanding of the molecular differences between primary tumors and metastases. The AURORA US Network was established to collect primary breast cancer-metastasis pairs for multi-platform genomic profiling in order to identify the molecular drivers of metastatic disease. AURORA US has both a retrospective and prospective phase. This is the first report of the retrospective phase.

Methods: Archived tissue samples from the primary tumor and at least one distant metastasis were retrospectively collected from 83 MBC patients. Following internal quality assessment, samples from 55 pts, including 105 distinct metastatic lesions, were subject to DNA low pass whole genome and exome sequencing, DNA methylation arrays, and RNA sequencing. Early analyses of these multi-platform data include: DNA methylation, tumor gene expression and microenvironmental signatures, somatic and germline variants, DNA copy number changes, and structural variants between breast primaries and matched metastases.

Results: Median age at diagnosis was 49 years (25-76); 32 (58%) were Stage I or II at presentation, 27 (49%) had a family history of breast cancer, and 20 (36%) had a second breast primary. Median disease-free interval before developing MBC was 2 years (range 0-36, 5 patients presented with Stage IV). Median overall survival from initial presentation was 4 years (range 0-37). Median survival after developing MBC was 1 year (range 0-13), with a median of three treatments. Primary tissue samples were banked from 1977-2017 and metastases were banked from 1999-2017. Clinical phenotypes of the primaries included 27 HR+ (49%), 15 triple negative (TNBC, 27%), and 11 HER2+ (20%, 12 missing HER2 status). Intrinsic subtype distribution of the primaries included 17 Basal-like (31%), 17 Luminal A (31%), 7 Normal-like (13%), 5 HER2-enriched (9%), and 1 Luminal B, with 8 pending. All metastases from the Basal-like cases remained Basal-like, while metastases from luminal primaries tended to gain HER2-Enriched subtype features (5/18, p = 0.01). Overall, we identified significant metastasis-enriched alterations in metabolism pathways, an increase in proliferation, and the loss of differentiation signatures and immune infiltrates with progression; the latter being the most pronounced in brain metastases. The most frequent somatic mutations in this cohort were in TP53, NCOR1, and RUNX1. Interestingly, ERBB2, EGFR, and ATM were also mutated in ≥10% of the tumors sequenced. In almost all cases, CpG island hypermethylation was clonally present in the primary tumor and persisted stably in the majority of metastatic lesions. Promoter CpG island hypermethylation was also identified in some metastatic lesions at JAM3, an important cellular adhesion molecule,and this was accompanied by reduced mRNA expression.

CONCLUSIONS: Collection of banked primary and metastatic tissue pairs identified a young MBC cohort with a high frequency of breast cancer family history and second breast primaries. Molecular characterization of luminal tumor pairs highlighted acquisition of aggressive traits including increased proliferation and loss of differentiation in the metastases. In contrast, basal-like pairs remained relatively unchanged, except for the loss of immune activation. Ongoing analyses to be presented include clonal heterogeneity and phylogeny, novel metastasis signature discovery, gene fusion, and endogenous retrovirus detection.

Citation Format: Tari A King, Minetta C Liu, Marni B McClure, Toshinori Hinoue, Benjamin J Kelly, Chad J Creighton, Jay Bowen, Kristen Leraas, Robyn T Burns, Sara Coppens, Salma Rezk, Amy L Garrett, Justin M Balko, Joel S Parker, Ben H Park, Ian Krop, Carey Anders, Katherine A Hoadley, Julie Gastier-Foster, Mothaffar F Rimawi, Rita Nanda, Nancy U Lin, Claudine Isaacs, P. Kelly Marcom, Anna Maria Storniolo, Fergus J Couch, Elaine R Mardis, Adrian V Lee, Uma Chandran, Peter W Laird, Susan G Hilsenbeck, Larry Norton, Andrea L Richardson, W. Fraser Symmans, Lisa A Carey, Antonio C Wolff, Nancy E Davidson, Charles M Perou, the AURORA US Network. Multiplatform analysis of matched primary and metastatic breast tumors from the AURORA US Network [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr GS3-08.

Augmentation of Therapy for Combined Loss of Heterozygosity 1p and 16q in Favorable Histology Wilms Tumor: A Children's Oncology Group AREN0532 and AREN0533 Study Report

PURPOSE

In National Wilms Tumor Study 5 (NWTS-5), tumor-specific combined loss of heterozygosity of chromosomes 1p and 16q (LOH1p/16q) was associated with adverse outcomes in patients with favorable histology Wilms tumor. The AREN0533/AREN0532 studies assessed whether augmenting therapy improved event-free survival (EFS) for these patients. Patients with stage I/II disease received regimen DD4A (vincristine, dactinomycin and doxorubicin) but no radiation therapy. Patients with stage III/IV disease received regimen M (vincristine, dactinomycin, and doxorubicin alternating with cyclophosphamide and etoposide) and radiation therapy.

METHODS

Patients were enrolled through the AREN03B2 Biology study between October 2006 and October 2013; all underwent central review of pathology, surgical reports, and imaging. Tumors were evaluated for LOH1p/16q by microsatellite testing. EFS and overall survival were compared using the log-rank test between NWTS-5 and current studies.

RESULTS

LOH1p/16q was detected in 49 of 1,147 evaluable patients with stage I/II disease (4.27%) enrolled in AREN03B2; 32 enrolled in AREN0532. LOH1p/16q was detected in 82 of 1,364 evaluable patients with stage III/IV disease (6.01%) in AREN03B2; 51 enrolled in AREN0533. Median follow-up for 83 eligible patients enrolled in AREN0532/0533 was 5.73 years (range, 2.84 to 9.63 years). The 4-year EFS for patients with stage I/II and stage III/IV disease with LOH1p/16 was 87.3% (95% CI, 75.1% to 99.5%) and 90.2% (95% CI, 81.8% to 98.6%), respectively. These results are improved compared with the NWTS-5 updated 4-year EFS of 68.8% for patients with stage I/II disease (P = .042), and 61.3% for patients with stage III/IV disease (P = .001), with trends toward improved 4-year overall survival. The most common grade 3 or higher nonhematologic toxicities with regimen M were febrile neutropenia (39.2%) and infections (21.6%).

CONCLUSION

Augmentation of therapy improved EFS for patients with favorable histology Wilms tumor and LOH1p/16q compared with the historical NWTS-5 comparison group, with an expected toxicity profile.

Maintaining Outstanding Outcomes Using Response- and Biology-Based Therapy for Intermediate-Risk Neuroblastoma: A Report From the Children's Oncology Group Study ANBL0531

PURPOSE

The primary objective of the Children's Oncology Group study ANBL0531 (ClinicalTrials.gov identifier: NCT00499616) was to reduce therapy for subsets of patients with intermediate-risk neuroblastoma using a biology- and response-based algorithm to assign treatment duration while maintaining a 3-year overall survival (OS) of 95% or more for the entire cohort.

PATIENTS AND METHODS

Children younger than age 12 years with intermediate-risk stage 2A/2B or stage 3 tumors with favorable histology; infants younger than age 365 days with stage 3, 4 or 4S disease; and toddlers from 365 to younger than 547 days with favorable histology, hyperdiploid stage 4, or unfavorable histology stage 3 tumors were eligible. Patients with MYCN-amplified tumors were excluded. Patients were assigned to initially receive two (group 2), four (group 3), or eight (group 4) cycles of chemotherapy with or without surgery on the basis of prognostic markers, including allelic status of chromosomes 1p and 11q; ultimate duration of therapy was determined by overall response.

RESULTS

Between 2007 and 2011, 404 evaluable patients were enrolled. Compared with legacy Children's Oncology Group studies, subsets of patients had a reduction in treatment. The 3-year event-free survival and OS rates were 83.2% (95% CI, 79.4% to 87.0%) and 94.9% (95% CI, 92.7% to 97.2%), respectively. Infants with stage 4 tumors with favorable biology (n = 61) had superior 3-year event-free survival compared with patients with one or more unfavorable biologic features (n = 47; 86.9% [95% CI, 78.3% to 95.4%] v 66.8% [95% CI, 53.1% to 80.6%]; P = .02), with a trend toward OS advantage (95.0% [95% CI, 89.5% to 100%] v 86.7% [95% CI, 76.6% to 96.7%], respectively; P = .08). OS for patients with localized disease was 100%.

CONCLUSION

Excellent survival was achieved with this treatment algorithm, with reduction of therapy for subsets of patients. More-effective treatment strategies still are needed for infants with unfavorable biology stage 4 disease.

Masked hypodiploidy: Hypodiploid acute lymphoblastic leukemia (ALL) mimicking hyperdiploid ALL in children: A report from the Children's Oncology Group

Hyperdiploidy with greater than 50 chromosomes is usually associated with favorable prognosis in pediatric acute lymphoblastic leukemia (ALL), whereas hypodiploidy with ≤43 chromosomes is associated with extremely poor prognosis. Sometimes, hypodiploidy is "masked" and patients do not have a karyotypically visible clone with ≤43 chromosomes. Instead, their abnormal karyotypes contain 50-78 or more chromosomes from doubling of previously hypodiploid cells. When the hypodiploid and doubled hyperdiploid clones are both present, patients can be identified by traditional test methods [karyotype, DNA Index (DI), fluorescence in situ hybridization (FISH)], but the incidence of masked hypodiploid cases in which only the doubled clone is visible is unknown. We analyzed 7013 patients with B-ALL enrolled in COG AALL03B1 (2003-2011) for whom chromosome studies were available. Of 115 patients with hypodiploidy (25-39 chromosomes), karyotypes of 40 showed only the hypodiploid clone, 47 showed mosaicism with both hypodiploid and hyperdiploid (doubled) karyotypes, and 28 with masked hypodiploidy showed only a hyperdiploid (doubled) clone. Unique karyotypic signatures were identified, and widespread loss of heterozygosity (LOH) was seen in the microsatellite panel for all patients with masked hypodiploidy. An increased awareness of the unusual karyotypic profile associated with a doubled hypodiploid clone and coordinated use of DI, FISH, and LOH studies when indicated can identify patients with masked hypodiploidy and allow appropriate treatment selection.

No evidence that G6PD deficiency affects the efficacy or safety of daunorubicin in acute lymphoblastic leukemia induction therapy

BACKGROUND/OBJECTIVES

Anthracyclines are used in induction therapy of pediatric acute lymphoblastic leukemia (ALL) and are known to generate oxidative stress; whether this translates into enhanced antileukemic activity or hemolytic effects in patients with glucose-6-phosphate dehydrogenase (G6PD) deficiency is unknown.

DESIGN/METHODS

Among 726 pediatric patients with newly diagnosed ALL treated at St. Jude Children's Research Hospital, 22 had deficient G6PD activity. We compared the prevalence of positive minimal residual disease (MRD) ≥1% at Day 15/Day 19 of induction or ≥0.01% at Day 42/Day 46 (end of induction) and the number of red blood cell (RBC) transfusions after daunorubicin in induction between patients with or without G6PD deficiency, adjusting for ALL risk group, treatment protocol, age, and gender.

RESULTS

There was no difference in Day 15/19 (P = 1) or end of induction MRD (P = 0.76) nor in the number of RBC transfusions (P = 0.73); the lack of association with MRD was confirmed in a dataset of 1192 newly diagnosed male patients enrolled in a Children's Oncology Group trial (P = 0.78).

CONCLUSION

We found no evidence that G6PD deficiency affects daunorubicin activity during induction treatment for ALL.

PAX5-driven subtypes of B-progenitor acute lymphoblastic leukemia

Recent genomic studies have identified chromosomal rearrangements defining new subtypes of B-progenitor acute lymphoblastic leukemia (B-ALL), however many cases lack a known initiating genetic alteration. Using integrated genomic analysis of 1,988 childhood and adult cases, we describe a revised taxonomy of B-ALL incorporating 23 subtypes defined by chromosomal rearrangements, sequence mutations or heterogeneous genomic alterations, many of which show marked variation in prevalence according to age. Two subtypes have frequent alterations of the B lymphoid transcription-factor gene PAX5. One, PAX5alt (7.4%), has diverse PAX5 alterations (rearrangements, intragenic amplifications or mutations); a second subtype is defined by PAX5 p.Pro80Arg and biallelic PAX5 alterations. We show that p.Pro80Arg impairs B lymphoid development and promotes the development of B-ALL with biallelic Pax5 alteration in vivo. These results demonstrate the utility of transcriptome sequencing to classify B-ALL and reinforce the central role of PAX5 as a checkpoint in B lymphoid maturation and leukemogenesis.

Abstract 4871: Whole-exome sequencing identifies a high frequency of germline deleterious variants in cancer predisposition genes in individuals with osteosarcoma

In children and adolescents, osteosarcoma (OS) is the most common malignant bone tumor. OS occurs in certain cancer predisposition syndromes at a higher than expected frequency. However, outside of these rare syndromes, OS is significantly more common and its genetic etiology remains poorly understood. We conducted an evaluation of rare exonic variants in 545 unselected OS cases compared with 1061 cancer-free controls using whole-exome sequencing of blood or buccal cell DNA to estimate the prevalence and burden of rare deleterious germline variants.

We assessed potentially pathogenic rare variants in 126 established cancer predisposing genes (CPG) and known somatically mutated genes. We also surveyed the exome for genes with a higher burden of rare variants in 342 EUR cases with 994 EUR controls sequenced at the same time. Rare genetic variants, defined by a MAF &lt;0.01 in 1000G, ESP or ExAC, that passed quality control filters were classified based on well-annotated databases (IARC-TP53, Telomerase Database, ARUP, BIC, Insight and COSMIC-hotspots), variant databases (e.g. Clinvar, HGMD and LOVD), and in silico predictions (Sift, Polyphen, Mutation assessor, MutationTaster, FATHMM and LTR). Rare variants matching known pathogenic variants, or inexact matches at the same residue, in selected well-annotated databases were labeled as putative pathogenic variants. High impact variants (frameshift, stop gain/loss, or known splice sites) not matching a previously identified pathogenic variant in the selected databases were categorized as probably pathogenic. Predictions based solely on in silico predictions were categorized as possibly pathogenic, uncertain or likely benign.

A total of 14.5% of cases had a predicted pathogenic or high impact variant in an autosomal dominant CPG. A significantly higher pathogenic rare variant burden was present in EUR cases compared with EUR controls (Pburden=3.7x10-05). TP53 had the highest prevalence of pathogenic mutations (5% of EUR cases; Pburden=1.05x10-09). CDKN2A, MEN1, MLH1, MUTYH, PALB2 and VHL had a significantly higher rare variant burden in the EUR cases. 7.5% of cases had a predicted pathogenic or probably pathogenic variant in an autosomal recessive CPG (Pburden=5.4x10-03). Additionally, two males had a pathogenic variant in the X-linked genes, DKC1 and GPC3. In total, 21% of cases had a predicted pathogenic or high impact variant in a CPG (Pburden=2.4x10-07). Exome-wide analysis identified two novel genes with significantly higher rare variant burdens in the cases compared to controls.

In conclusion, several CPGs and two novel genes, not previously associated with OS, had an enrichment of rare variants in OS and warrant further follow-up. Our results indicate that a clinically significant fraction of OS cases may harbor one or more mutations worthy of consideration for further investigation and genetic counseling.

Citation Format: Roelof Koster, Bin Zhu, Meredith Yeager, Michael Dean, Matthew Gianferante, Lei Song, Joshua Sampson, NCI DCEG Cancer Genomics Research Laboratory, Julie Gastier-Foster, Richard Gorlick, Silvia Regina Caminada de Toledo, Antonio Petrilli, Ana Patiño-Garcia, Fernando Lecanda, Massimo Serra, Claudia Hattinger, Piero Picci, Katia Scotlandi, Adrienne Flanagan, Roberto Tirabosco, Maria Amary, Nilgün Kurucu, Inci Ergurhan Ilhan, Neriman Sari, Mandy Ballinger, David Thomas, Donald Barkauskas, Belynda Hicks, Margaret Tucker, Neil Caporaso, Robert Hoover, Stephen Chanock, Sharon Savage, Lisa Mirabello. Whole-exome sequencing identifies a high frequency of germline deleterious variants in cancer predisposition genes in individuals with osteosarcoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4871. doi:10.1158/1538-7445.AM2017-4871

Truncating Erythropoietin Receptor Rearrangements in Acute Lymphoblastic Leukemia

Chromosomal rearrangements are a hallmark of acute lymphoblastic leukemia (ALL) and are important ALL initiating events. We describe four different rearrangements of the erythropoietin receptor gene EPOR in Philadelphia chromosome-like (Ph-like) ALL. All of these rearrangements result in truncation of the cytoplasmic tail of EPOR at residues similar to those mutated in primary familial congenital polycythemia, with preservation of the proximal tyrosine essential for receptor activation and loss of distal regulatory residues. This resulted in deregulated EPOR expression, hypersensitivity to erythropoietin stimulation, and heightened JAK-STAT activation. Expression of truncated EPOR in mouse B cell progenitors induced ALL in vivo. Human leukemic cells with EPOR rearrangements were sensitive to JAK-STAT inhibition, suggesting a therapeutic option in high-risk ALL.

Analysis of Heritability and Shared Heritability Based on Genome-Wide Association Studies for Thirteen Cancer Types

BACKGROUND

Studies of related individuals have consistently demonstrated notable familial aggregation of cancer. We aim to estimate the heritability and genetic correlation attributable to the additive effects of common single-nucleotide polymorphisms (SNPs) for cancer at 13 anatomical sites.

METHODS

Between 2007 and 2014, the US National Cancer Institute has generated data from genome-wide association studies (GWAS) for 49 492 cancer case patients and 34 131 control patients. We apply novel mixed model methodology (GCTA) to this GWAS data to estimate the heritability of individual cancers, as well as the proportion of heritability attributable to cigarette smoking in smoking-related cancers, and the genetic correlation between pairs of cancers.

RESULTS

GWAS heritability was statistically significant at nearly all sites, with the estimates of array-based heritability, hl (2), on the liability threshold (LT) scale ranging from 0.05 to 0.38. Estimating the combined heritability of multiple smoking characteristics, we calculate that at least 24% (95% confidence interval [CI] = 14% to 37%) and 7% (95% CI = 4% to 11%) of the heritability for lung and bladder cancer, respectively, can be attributed to genetic determinants of smoking. Most pairs of cancers studied did not show evidence of strong genetic correlation. We found only four pairs of cancers with marginally statistically significant correlations, specifically kidney and testes (ρ = 0.73, SE = 0.28), diffuse large B-cell lymphoma (DLBCL) and pediatric osteosarcoma (ρ = 0.53, SE = 0.21), DLBCL and chronic lymphocytic leukemia (CLL) (ρ = 0.51, SE =0.18), and bladder and lung (ρ = 0.35, SE = 0.14). Correlation analysis also indicates that the genetic architecture of lung cancer differs between a smoking population of European ancestry and a nonsmoking Asian population, allowing for the possibility that the genetic etiology for the same disease can vary by population and environmental exposures.

CONCLUSION

Our results provide important insights into the genetic architecture of cancers and suggest new avenues for investigation.

Clinical Application of Prognostic Gene Expression Signature in Fusion Gene-Negative Rhabdomyosarcoma: A Report from the Children's Oncology Group

PURPOSE

Pediatric rhabdomyosarcoma (RMS) has two common histologic subtypes: embryonal (ERMS) and alveolar (ARMS). PAX-FOXO1 fusion gene status is a more reliable prognostic marker than alveolar histology, whereas fusion gene-negative (FN) ARMS patients are clinically similar to ERMS patients. A five-gene expression signature (MG5) previously identified two diverse risk groups within the fusion gene-negative RMS (FN-RMS) patients, but this has not been independently validated. The goal of this study was to test whether expression of the MG5 metagene, measured using a technical platform that can be applied to routine pathology material, would correlate with outcome in a new cohort of patients with FN-RMS.

EXPERIMENTAL DESIGN

Cases were taken from the Children's Oncology Group (COG) D9803 study of children with intermediate-risk RMS, and gene expression profiling for the MG5 genes was performed using the nCounter assay. The MG5 score was correlated with clinical and pathologic characteristics as well as overall and event-free survival.

RESULTS

MG5 standardized score showed no significant association with any of the available clinicopathologic variables. The MG5 signature score showed a significant correlation with overall (N = 57; HR, 7.3; 95% CI, 1.9-27.0; P = 0.003) and failure-free survival (N = 57; HR, 6.1; 95% CI, 1.9-19.7; P = 0.002).

CONCLUSIONS

This represents the first, validated molecular prognostic signature for children with FN-RMS who otherwise have intermediate-risk disease. The capacity to measure the expression of a small number of genes in routine pathology material and apply a simple mathematical formula to calculate the MG5 metagene score provides a clear path toward better risk stratification in future prospective clinical trials.

Abstract 5574: High prevalence of germline TP53 mutations in young osteosarcoma cases

Osteosarcoma, the most common primary bone malignancy, has a bimodal age distribution, with a primary peak in adolescence and a smaller peak in the elderly. The etiologic contribution of germline genetic variation to osteosarcoma is not well-understood. It occurs at higher than expected frequencies in individuals with the Li-Fraumeni syndrome (LFS) cancer predisposition syndrome. Approximately 70% of classic LFS families have germline TP53 mutations. Two previous studies reported that 3% of young osteosarcoma cases (&lt;20 years old, N = 235) and 7% of all aged cases (N = 95) had germline TP53 mutations.

We determined the prevalence of germline TP53 mutations in 765 unselected osteosarcoma cases. DNA was extracted from blood and TP53 sequenced using custom Ampliseq panels. Variants were validated with Sanger sequencing. The IARC germline TP53 database was used to identify TP53 mutations reported in families with LFS (LFS-associated mutations). Variants were considered “likely LFS-associated mutations” if absent from publically available databases (ESP and 1,000 Genomes Project) and predicted non-functional or deleterious using in silico algorithms. Variants were considered “rare exonic variants” if their minor allele frequency (MAF) was &lt;2% in public databases and they had uncertain clinical significance.

There were 32 LFS-associated or rare TP53 variants in 62 osteosarcoma cases. The frequency of cases with an LFS or likely LFS-associated mutation and/or rare exonic variant was 8.1%. Notably, all 32 TP53 variants were present in cases &lt;30 years of age (“young cases”, N = 505), 9.5% of young cases, compared with none in older cases (N = 51; P&lt;0.001). TP53 variants did not have an even distribution within the first three age decades in young cases. LFS or likely LFS-associated mutations, which confer significant cancer risk, had the highest frequency in patients aged 0-9 years (4.8% of all cases, and 6.1% of European cases), and rare exonic variants were most frequent in patients aged 10-19 years (6.1% of all cases, and 5.6% of European cases).

A logistic regression case-case analysis identified a novel significant association between a rare TP53 variant, rs1800372 (p.R213R), and metastasis at diagnosis in cases of European ancestry (odds ratio [OR] 4.27, 95% CI 1.2-15.5, P = 0.026). We additionally confirmed that a common exonic variant, rs1042522 (p.P72R), was significantly associated with osteosarcoma risk (OR 1.22, 95% CI 1.1-1.4, P = 0.0098) and poorer survival (HR 1.35, 95% CI 1.00-1.83, P = 0.048).

Our data suggest that genetic susceptibility to young onset osteosarcoma is distinct from adult onset osteosarcoma. The high TP53 mutation prevalence we identified in osteosarcoma cases aged &lt;20 years of 9.7% is significantly greater than the previously reported prevalence of 3% in unselected cases (P = 0.002). Based on these findings, genetic counseling and TP53-mutation testing of all young patients with osteosarcoma should be considered.

Citation Format: Lisa J. Mirabello, Meredith Yeager, Phuong L. Mai, Julie Gastier-Foster, Richard Gorlick, Chand Khanna, Ana Patiño-Garcia, Luis Sierrasesúmaga, Fernando Lecanda, Irene L. Andrulis, Jay S. Wunder, Nalan Gokgoz, Donald A. Barkauskas, Xijun Zhang, Aurelie Vogt, Kristine Jones, Joseph F. Boland, Stephen J. Chanock, Sharon A. Savage. High prevalence of germline TP53 mutations in young osteosarcoma cases. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 5574. doi:10.1158/1538-7445.AM2015-5574

Abstract 484: Fidelity of subclonal representation in human neuroblastoma-derived cell line and patient-derived xenograft models: A report from the NCI-TARGET project

Patient-derived xenografts and cell lines have been the underpinning of functional characterization and drug discovery efforts in cancer. The use of these models is often under the assumption that these systems are renewable, faithful representations of original progenitor tumor cell populations. To test this assumption, we performed whole exome (92X median coverage) and genome sequencing (36X median coverage) analysis of cell line and patient-derived mouse xenografts (PDXs) originating from 7 neuroblastoma patients. Data are from1 primary tumor, 3 PDXs, and 15 neuroblastoma cell lines cultured from tumor, bone marrow, or blood. The cell lines consisted of 4 pre-/post-therapy pairs and 3 pairs established and maintained in either hyperoxia (room air i.e. “standard” cell culture) or physiologic (bone marrow hypoxia = 5%) oxygen. 7 lymphoblastoid or fibroblast cell lines were used as matched normals to identify somatic mutations. Subclonal population structures were inferred from somatic mutation calls calibrated for copy number state and tumor purity.

In all cell lines and xenografts, we observed 1-2 additional subclonal populations, primarily supported by deep coverage from exome sequencing. In nearly every case, we observed shifts in the proportional representation of genetic subclones and many subclones showed additional mutations not evident in the progenitor tissue or cancer line derived in parallel. Comparison of three cell line pairs established in bone marrow level hypoxia versus room air found only ∼40% of coding mutations in each line were shared (average 82 mutations per line), suggesting significant genetic impact of growing tumor cells in the two different culture conditions. Matched PDXs from these cases had only ∼17% of coding mutations shared across all three models. The greatest genetic similarity was seen between paired cell lines established from tissue obtained pre-/post-therapy from the same patient (36 coding mutations shared, 14 private to diagnosis and 13 private to progression). However, a second pre/post-therapy cell line pair did not share any coding mutations, although they did have 585 non-coding mutations in common (of 4,033 and 2,480 in each line), assuring that the relapse was derived from a diagnostic tumor clone. These results highlight a need for comprehensive subclonal analysis of human cancer laboratory models to better inform design and interpretation of biological and preclinical therapeutic studies.

Citation Format: Maya Schonbach, Arnavaz Danesh, Jeff Bruce, Tito Woodburn, Tanja Davidsen, Leandro Hermida, Patee Gesuwan, Jaime Guidry Auvil, Oliver Hampton, David Wheeler, Julie Gastier-Foster, Malcolm Smith, Daniela Gerhard, John M. Maris, Patrick Reynolds, Trevor J. Pugh. Fidelity of subclonal representation in human neuroblastoma-derived cell line and patient-derived xenograft models: A report from the NCI-TARGET project. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 484. doi:10.1158/1538-7445.AM2015-484

Comprehensive, Integrative Genomic Analysis of Diffuse Lower-Grade Gliomas

BACKGROUND

Diffuse low-grade and intermediate-grade gliomas (which together make up the lower-grade gliomas, World Health Organization grades II and III) have highly variable clinical behavior that is not adequately predicted on the basis of histologic class. Some are indolent; others quickly progress to glioblastoma. The uncertainty is compounded by interobserver variability in histologic diagnosis. Mutations in IDH, TP53, and ATRX and codeletion of chromosome arms 1p and 19q (1p/19q codeletion) have been implicated as clinically relevant markers of lower-grade gliomas.

METHODS

We performed genomewide analyses of 293 lower-grade gliomas from adults, incorporating exome sequence, DNA copy number, DNA methylation, messenger RNA expression, microRNA expression, and targeted protein expression. These data were integrated and tested for correlation with clinical outcomes.

RESULTS

Unsupervised clustering of mutations and data from RNA, DNA-copy-number, and DNA-methylation platforms uncovered concordant classification of three robust, nonoverlapping, prognostically significant subtypes of lower-grade glioma that were captured more accurately by IDH, 1p/19q, and TP53 status than by histologic class. Patients who had lower-grade gliomas with an IDH mutation and 1p/19q codeletion had the most favorable clinical outcomes. Their gliomas harbored mutations in CIC, FUBP1, NOTCH1, and the TERT promoter. Nearly all lower-grade gliomas with IDH mutations and no 1p/19q codeletion had mutations in TP53 (94%) and ATRX inactivation (86%). The large majority of lower-grade gliomas without an IDH mutation had genomic aberrations and clinical behavior strikingly similar to those found in primary glioblastoma.

CONCLUSIONS

The integration of genomewide data from multiple platforms delineated three molecular classes of lower-grade gliomas that were more concordant with IDH, 1p/19q, and TP53 status than with histologic class. Lower-grade gliomas with an IDH mutation either had 1p/19q codeletion or carried a TP53 mutation. Most lower-grade gliomas without an IDH mutation were molecularly and clinically similar to glioblastoma. (Funded by the National Institutes of Health.).

Risk factors and timing of relapse after allogeneic transplantation in pediatric ALL: for whom and when should interventions be tested?

We previously showed that minimal residual disease (MRD) detection pre-hematopoietic cell transplant (HCT) and acute GvHD (aGvHD) independently predicted risk of relapse in pediatric ALL. In this study we further define risk by assessing timing of relapse and the effects of leukemia risk category and post-HCT MRD. By multivariate analysis, pre-HCT MRD <0.1% and aGvHD by day +55 were associated with decreased relapse and improved event-free survival (EFS). Intermediate leukemia risk status predicted decreased relapse, and improved EFS and overall survival (OS). Patients with pre-HCT MRD ⩾0.1% who did not develop aGvHD compared with those with MRD <0.1% who did develop aGvHD had much worse survival (2 years EFS 18% vs 71%; P=0.001, 2 years OS 46 vs 74%; P=0.04). Patients with pre-HCT MRD <0.1% who did not experience aGvHD had higher rates of relapse than those who did develop aGvHD (40% vs 13%; P= 0.008). Post-HCT MRD led to a substantial increase in relapse risk (HR=4.5, P<0.01). Patients at high risk of relapse can be defined after transplant using leukemia risk category, presence of MRD pre or post HCT, and occurrence of aGvHD. An optimal window to initiate intervention to prevent relapse occurs between day +55 and +200 after HCT.

A slowly progressive form of limb-girdle muscular dystrophy type 2C associated with founder mutation in the SGCG gene in Puerto Rican Hispanics

Limb-girdle muscular dystrophy type 2C (LGMD2C) is considered one of the severe forms of childhood-onset muscular dystrophy. The geographical distribution of founder mutations in the SGCG gene has a prominent effect on the prevalence of LGMD2C in certain populations. The aim of this study was to confirm the hypothesis that the c.787G>A (p.E263K) mutation in the SGCG gene is a founder mutation among Puerto Rican Hispanics and to characterize the associated clinical and immunohistochemical phenotype. Genotyping of six polymorphic microsatellite markers internal to (D13S232) and flanking (D13S175, D13S292, D13S787, D13S1243, D13S283) the SGCG gene was performed on four unrelated Puerto Rican patients with LGMD2C. Preserved ambulation to the second decade of life was observed in at least two subjects. Immunostaining of skeletal muscle demonstrated absence of γ-sarcoglycan in all affected subjects. Two markers, D13S232 and D13S292, were highly informative and confirmed that all four families share the haplotype of the mutant allele. Our findings confirm that the E263K missense mutation in the SGCG gene is a founder mutation in Puerto Rican Hispanics. A slowly progressive disease course with prolonged preservation of ambulation can be seen in association with this mutation, providing evidence for phenotypic variability.

Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia

BACKGROUND

Many mutations that contribute to the pathogenesis of acute myeloid leukemia (AML) are undefined. The relationships between patterns of mutations and epigenetic phenotypes are not yet clear.

METHODS

We analyzed the genomes of 200 clinically annotated adult cases of de novo AML, using either whole-genome sequencing (50 cases) or whole-exome sequencing (150 cases), along with RNA and microRNA sequencing and DNA-methylation analysis.

RESULTS

AML genomes have fewer mutations than most other adult cancers, with an average of only 13 mutations found in genes. Of these, an average of 5 are in genes that are recurrently mutated in AML. A total of 23 genes were significantly mutated, and another 237 were mutated in two or more samples. Nearly all samples had at least 1 nonsynonymous mutation in one of nine categories of genes that are almost certainly relevant for pathogenesis, including transcription-factor fusions (18% of cases), the gene encoding nucleophosmin (NPM1) (27%), tumor-suppressor genes (16%), DNA-methylation-related genes (44%), signaling genes (59%), chromatin-modifying genes (30%), myeloid transcription-factor genes (22%), cohesin-complex genes (13%), and spliceosome-complex genes (14%). Patterns of cooperation and mutual exclusivity suggested strong biologic relationships among several of the genes and categories.

CONCLUSIONS

We identified at least one potential driver mutation in nearly all AML samples and found that a complex interplay of genetic events contributes to AML pathogenesis in individual patients. The databases from this study are widely available to serve as a foundation for further investigations of AML pathogenesis, classification, and risk stratification. (Funded by the National Institutes of Health.).

Assessment of 2q23.1 microdeletion syndrome implicates MBD5 as a single causal locus of intellectual disability, epilepsy, and autism spectrum disorder

Persons with neurodevelopmental disorders or autism spectrum disorder (ASD) often harbor chromosomal microdeletions, yet the individual genetic contributors within these regions have not been systematically evaluated. We established a consortium of clinical diagnostic and research laboratories to accumulate a large cohort with genetic alterations of chromosomal region 2q23.1 and acquired 65 subjects with microdeletion or translocation. We sequenced translocation breakpoints; aligned microdeletions to determine the critical region; assessed effects on mRNA expression; and examined medical records, photos, and clinical evaluations. We identified a single gene, methyl-CpG-binding domain 5 (MBD5), as the only locus that defined the critical region. Partial or complete deletion of MBD5 was associated with haploinsufficiency of mRNA expression, intellectual disability, epilepsy, and autistic features. Fourteen alterations, including partial deletions of noncoding regions not typically captured or considered pathogenic by current diagnostic screening, disrupted MBD5 alone. Expression profiles and clinical characteristics were largely indistinguishable between MBD5-specific alteration and deletion of the entire 2q23.1 interval. No copy-number alterations of MBD5 were observed in 7878 controls, suggesting MBD5 alterations are highly penetrant. We surveyed MBD5 coding variations among 747 ASD subjects compared to 2043 non-ASD subjects analyzed by whole-exome sequencing and detected an association with a highly conserved methyl-CpG-binding domain missense variant, p.79Gly>Glu (c.236G>A) (p = 0.012). These results suggest that genetic alterations of MBD5 cause features of 2q23.1 microdeletion syndrome and that this epigenetic regulator significantly contributes to ASD risk, warranting further consideration in research and clinical diagnostic screening and highlighting the importance of chromatin remodeling in the etiology of these complex disorders.

Abstract 4756: Exome sequencing of 81 neuroblastomas identifies a wide diversity of somatic mutation

Neuroblastoma (NBL) is a cancer arising from the peripheral nervous system and is the most common solid tumor outside the brain in infants and children. Despite improved treatment over the past 20 years, only 45% of those with high-risk disease survive emphasizing the need for new molecular targets for therapy. To identify genetic mutations that encode such targets, we present an analysis of 162 whole exome sequences from 81 tumor/normal pairs from high-risk NBL patients (N=27 with MYCN amplification). From each sample, genomic DNA fragments encoding 185,961 exons from 18,380 genes were isolated using in-solution hybrid capture and sequenced using the Illumina sequencing platform to a median of 195X coverage. A wide diversity of somatic variation was observed with non-silent mutation counts ranging from 5 to 149 per tumor (median=24). Novel hyper-mutation phenotypes (&gt; 125 candidate coding somatic mutations) were observed in four tumors including one case with both a deletion and a truncating mutation in MLH1. Several genes were mutated across the set at a statistically significant level including ALK (6 tumors, q=0.006) and PTPN11 (3 tumors, q=0.012). There are 27 other genes with candidate mutation frequencies in the 4-10% range, but no single gene identified to date with a mutation frequency &gt;10%. Importantly, multiple paired blood samples (acquired at diagnosis) have sufficient numbers of circulating tumor DNA to potentially interfere with our SNV detection algorithm; refinement of the analytic pipeline using SNP-derived “contamination” metrics is ongoing. We anticipate validating candidate somatic candidates using an orthogonal genotyping method and further defining their true frequency in an expanded sample set that includes non-high-risk cases as well. Finally, we plan a concerted effort to catalogue potentially detrimental germline variations in order to further extend our understanding of the heritable component of NBL.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4756. doi:10.1158/1538-7445.AM2011-4756

Abstract 5249: Identification of genome-wide alternative splicing events in high-risk MYCN amplified neuroblastomas

Introduction: High-risk neuroblastoma (H-NB) is an aggressive tumor of the sympathetic nervous system for which there are few predictors of therapeutic response. We sought to determine the association of global alternative splicing (AS) gene expression patterns with NB aggressiveness and MYCN gene status.

Methods: As part of the NCI-sponsored Therapeutically Applicable Research to Generate Effective Treatments (TARGET) initiative, we profiled 190 primary human NBs using Affymetrix Human Exon arrays. Tumors were obtained from children who participated in Children's Oncology Group trials with informed consent. The high-risk cohort included tumors with (H-NBA n=48) and without MYCN amplification (H-NBN n=111), a key molecular determinant of clinical aggressiveness. Twenty-four samples from low-risk NB lacking MYCN amplification (L-NB) were also profiled as a comparison group. Robust multichip analysis (RMA) was used for summarization of probe selection regions (PSR) and transcript regions. Finding isoforms using RMA (FIRMA) was used to generate PSR isoform scores. After an unsupervised pre-filtering step based on coefficient of variation, ANOVA analyses were conducted comparing H-NBA or H-NBN against the L-NB group for differential transcript and isoform expression. Last, gene set enrichment analyses (GSEA) comparing high-risk to low-risk NBs were conducted using unfiltered transcript data and the Molecular Signature Database's (MsigDB) canonical gene sets.

Results: 1914 out of the 22,011 core transcripts passed our strict prefiltering criteria that required a high coefficient of variation. There was a larger number of genes with significant differential expression between H-NBA vs. L-NB than H-NBN vs. L-NB (23.5% versus 4.6%, P&lt;0.001). There was also significantly larger number of PSRs exhibiting differential FIRMA scores between H-NBA vs. L-NB compared to H-NBN vs. L-NB (14.6% versus 2.3% of PSRs; P&lt;0.001). NTRK1, a gene with known alternative splicing pattern among neuroblastomas, was identified in our data to have lower expression in both H-NBA and H-NBN vs. the L-NB subgroup but showed significant evidence of alternative splicing only in the H-NBA group. GSEA analyses also identified more differentially expressed pathways in H-NBA vs. L-NB than in H-NBN versus L-NB neuroblastomas (66 versus 31 pathways; P&lt;0.001). Pathways related to cell cycle and axon guidance were among the top differentially expressed pathways in the H-NBA vs. L-NB comparison.

Conclusion: Significant differences in gene expression and pathway expression are identified among H-NBA compared to H-NBN tumors. In addition, we show for the first time a high degree of alternative splicing events in the H-NBA tumors. These findings could aid in understanding the regulatory networks in the two subgroups of high-risk neuroblastomas and suggest possible therapeutic targets for further validation.

Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 101st Annual Meeting of the American Association for Cancer Research; 2010 Apr 17-21; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2010;70(8 Suppl):Abstract nr 5249.