Eckel‐Passow, Jeanette E.

Abstract Background The four most commonly-mutated genes in clear cell renal cell carcinoma (ccRCC) tumors are BAP1 , PBRM1 , SETD2 and VHL . And, there are currently 14 known RCC germline variants that have been reproducibly shown to be associated with RCC risk. However, the association of germline genetics with tumor genetics and clinical aggressiveness are unknown. Methods We analyzed 420 ccRCC patients from The Cancer Genome Atlas (TCGA). Molecular subtype was determined based on acquired mutations in BAP1 , PBRM1 , SETD2 and VHL . Aggressive subtype was defined clinically using Mayo SSIGN score and molecularly using the ccA/ccB gene expression subtype. Publically-available Hi-C data were used to link germline risk variants with candidate target genes. Results The 8q24 variant rs35252396 was significantly associated with VHL mutation status (OR=1.6, p=0.0037) and SSIGN score (OR=1.9, p=0.00094), after adjusting for multiple comparisons. We observed that, while some germline variants have interactions with nearby genes, some variants demonstrate long-range interactions with target genes. Conclusions These data further demonstrate the link between rs35252396, HIF pathway and ccRCC clinical aggressiveness, providing a more comprehensive picture of how germline genetics and tumor genetics interact with respect to tumor development and progression.

Abstract Background: Genome-wide association studies (GWAS) revealed that 25 regions in 24 genes are associated with adult diffuse glioma development. These regions were identified by performing GWAS of glioma overall and by pathology (GBM and nonGBM). The 2016 WHO Classification of Tumors of the Central Nervous System utilizes two somatic alterations to molecularly-classify adult diffuse glioma: IDH mutation and 1p/19q codeletion. We hypothesized that germline variants may increase susceptibility to, or interact with, these somatic alterations to accelerate the development of specific molecular subtypes of glioma. We further hypothesize that germline variants associated with IDH-mutated glioma might be associated with other IDH-mutated tumors, namely, cholangiocarcinoma, acute myeloid leukemia (AML) and melanoma. Methods: We performed a GWAS by glioma molecular subtype - as defined by presence or absence of IDH somatic mutation and 1p/19q codeletion. A total of 1320 glioma cases and 1889 controls were used in the discovery set, and 799 glioma cases and 808 controls in the validation set. A meta-analysis was performed with a genome-wide p-value threshold of 5 × 10−8. GTEx data were used to perform an expression quantitative trait loci (eQTL) analysis. For germline variants that were significantly associated with IDH-mutated glioma, we evaluated pleiotropy with cholangiocarcinoma, AML and melanoma using TCGA and Mayo Biobank controls. Results: Variants in or near D2HGDH on chromosome 2 were genome-wide significant in IDH-mutated glioma (meta p-value = 2.82 × 10−10). TCGA reported that the D2HGDH region was commonly deleted in IDH-mutated gliomas that do not have 1p/19q codeletion. In TCGA data for IDH-mutated, non-codeleted glioma, we observed that the D2HGDH variant was inversely associated with tumor deletions of D2HGDH (odds ratio=0.57, p-value=0.015). The eQTL analyses demonstrated significant associations between D2HGDH germline variant and expression of D2HGDH (p=2.2 × 10−11). Further stratifying IDH-mutated glioma by 1p/19q codeletion status, one variant near FAM20C on chromosome 7 was genome-wide significant in gliomas that have IDH mutation and 1p/19q codeletion (meta p-value=9.56 × 10−9). Analyses are currently underway to evaluate pleiotropy of these IDH-mutated glioma germline variants with other IDH-mutated tumors including cholangiocarcinoma, AML and melanoma. Variants in or near GMEB2 on chromosome 20 were genome-wide significant in IDH wild-type glioma (meta p-value=2.60 × 10−10). The most significant variant in the GMEB2 region remained significant after adjustment for the known RTEL1 glioma risk variant nearby on chromosome 20 (p=0.029). Conclusions: We identified and validated novel germline variants in two genes that are associated with etiology of IDH-mutated and one gene that is associated with IDH wild-type adult diffuse glioma. Citation Format: Jeanette E. Eckel-Passow, Kristen L. Drucker, Thomas M. Kollmeyer, Matthew L. Kosel, Paul A. Decker, Annette M. Molinaro, Terri Rice, Corrine E. Praska, Lauren E. Clark, Alissa A. Caron, Alexej Abyzov, Anthony Batzler, Jun S. Song, Melike Pekmezci, Helen M. Hansen, Lucie S. McCoy, Paige M. Bracci, Joseph Wiemels, John K. Wiencke, Stephen Francis, Terence C. Burns, Caterina Giannini, Daniel H. Lachance, Margaret Wrensch, Robert B. Jenkins. Adult diffuse glioma GWAS by molecular subtype identifies variants in D2HGDH, FAM20C and GMEB2 [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1193.

Abstract The Mayo Clinic collection of patient-derived xenograft (PDX) models has been developed from the direct implantation of fresh glioblastoma samples into immunodeficient mice. Over 95 models have been extensively characterized with regards to histology, invasion, flank and orthotopic growth rates, and molecular profiling. Evaluation of orthotopic tumors demonstrated clear infiltration of tumor cells into the brain parenchymal in 90% of PDXs with contralateral hemisphere involvement in 64%. Extensive molecular profiling confirmed that our collection of PDXs captures the genetic heterogeneity of GBM, with the majority of known alterations represented at a frequency similar to the TCGA including gain of chromosome 7 and loss of chromosome 10. TERT promoter mutations were the most frequent overall (86%) followed by homozygous deletion of CDKN2A (70%). Multiple alterations were identified in the p53 pathway, including TP53 mutations (36%) and MDM2/4 (10%/2%) amplification. PI3-kinase pathway alterations mainly involved PTEN (48%) while MAP kinase pathway alterations were present in a subset of PDXs (NF1 alterations 17%, BRAF 4%). IDH mutations were found in only two models. Gene expression analysis revealed that our PDX collection includes all GBM subtypes with classical being the most represented (46% vs mesenchymal 31% and proneural 23%). MGMT promoter methylation was also observed in 45% of PDXs. The high fidelity of PDX models to their human counterpart makes these models ideal for testing standard of care and novel targeted therapies. To assess if PDXs are predictive of patient treatment response, 37 orthotopic PDX lines were treated with radiation (RT), temozolomide (TMZ), or concurrent RT/TMZ. Across all PDXs, RT and TMZ demonstrated median survival extension of 90 and 46 days relative to vehicle. Concurrent RT/TMZ further extended survival by 146 days. MGMT promoter methylation was strongly associated with response to TMZ-based strategies but not RT alone.This observation suggests that our PDX collection possesses clinically relevant predictive value. To date, we have characterized treatment responses to over 135 regimens including both GBM standards of care and targeted strategies in our PDX models. Recent lab efforts strive to define the effect of blood brain barrier disruption on drug delivery. With over 20 years of experience, our group has extensive expertise in the development of clinically relevant study designs focused on examining efficacy, tolerability, and pharmacokinetic/pharmacodynamic assessments. Since 2016, we have shared over 5200 individual samples (tumor tissue, cells, DNA/RNA, protein extracts, tissue microarray slides) in both internal and external collaborators. Our models have been utilized in over 41 NIH funded grants, including 33 funded RO1s, as well as by the pharmaceutical community with 50 working relationships to date (20 contracts, 30 collaborations). Finally, our PDX collection has been used in over 120 peer-reviewed published manuscripts. Citation Format: Danielle Burgenske, Ann C. Mladek, Rachael A. Vaubel, Shulan Tian, Mark A. Schroeder, Zeng Hu, Brett L. Carlson, Paul A. Decker, Jeanette E. Eckel-Passow, Jann N. Sarkaria. Evaluation of novel therapeutics using models from the Mayo Clinic GBM patient-derived xenograft (PDX) collection [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 1113.

Abstract Background: The four most commonly-mutated genes in clear cell renal cell carcinoma (ccRCC) tumors are BAP1, PBRM1, SETD2 and VHL. And, there are currently 14 known RCC germline variants that have been reproducibly shown to be associated with RCC risk. However, the association of germline genetics with tumor genetics and clinical aggressiveness are unknown.Methods: We analyzed 420 ccRCC patients from The Cancer Genome Atlas (TCGA). Molecular subtype was determined based on acquired mutations in BAP1, PBRM1, SETD2 and VHL. Aggressive subtype was defined clinically using Mayo SSIGN score and molecularly using the ccA/ccB gene expression subtype. Publically-available Hi-C data were used to link germline risk variants with candidate target genes. Results: The 8q24 variant rs35252396 was significantly associated with VHL mutation status (OR=1.6, p=0.0037) and SSIGN score (OR=1.9, p=0.00094), after adjusting for multiple comparisons. We observed that, while some germline variants have interactions with nearby genes, some variants demonstrate long-range interactions with target genes.Conclusions: These data further demonstrate the link between rs35252396, HIF pathway and ccRCC clinical aggressiveness, providing a more comprehensive picture of how germline genetics and tumor genetics interact with respect to tumor development and progression.

Abstract Background Twenty-five germline variants are associated with adult diffuse glioma, and some of these variants have been shown to be associated with particular subtypes of glioma. We hypothesized that additional germline variants could be identified if a genome-wide association study (GWAS) were performed by molecular subtype. Methods A total of 1320 glioma cases and 1889 controls were used in the discovery set and 799 glioma cases and 808 controls in the validation set. Glioma cases were classified into molecular subtypes based on combinations of isocitrate dehydrogenase (IDH) mutation, telomerase reverse transcriptase (TERT) promoter mutation, and 1p/19q codeletion. Logistic regression was applied to the discovery and validation sets to test for associations of variants with each of the subtypes. A meta-analysis was subsequently performed using a genome-wide P-value threshold of 5 × 10−8. Results Nine variants in or near D-2-hydroxyglutarate dehydrogenase (D2HGDH) on chromosome 2 were genome-wide significant in IDH-mutated glioma (most significant was rs5839764, meta P = 2.82 × 10−10). Further stratifying by 1p/19q codeletion status, one variant in D2HGDH was genome-wide significant in IDH-mutated non-codeleted glioma (rs1106639, meta P = 4.96 × 10−8). Further stratifying by TERT mutation, one variant near FAM20C (family with sequence similarity 20, member C) on chromosome 7 was genome-wide significant in gliomas that have IDH mutation, TERT mutation, and 1p/19q codeletion (rs111976262, meta P = 9.56 × 10−9). Thirty-six variants in or near GMEB2 on chromosome 20 near regulator of telomere elongation helicase 1 (RTEL1) were genome-wide significant in IDH wild-type glioma (most significant was rs4809313, meta P = 2.60 × 10−10). Conclusions Performing a GWAS by molecular subtype identified 2 new regions and a candidate independent region near RTEL1, which were associated with specific glioma molecular subtypes.

Abstract Background: B-chronic lymphocytic leukemia (B-CLL) is the most common leukemia in the Western world but its pathogenesis remains largely unknown. Our strategy was to employ gene expression profiling (GEP) to discover genetic differences between CLL B cells and blood B cells from healthy adults. Previously, we and others identified the transcription factor (TF) lymphocyte enhancer factor-1 (LEF-1) as one of several genes significantly over expressed in CLL B cells suggesting a role for the oncogenic Wnt signaling pathway. LEF-1 acts as a central mediator of Wnt signaling and is crucial for the proliferation and survival of pro-B cells during development. Recently, deregulated LEF-1 activation has been directly linked with leukemogenesis. Methodology: The goals of our study were to: validate LEF-1 over expression in CLL B cells in a new cohort of B-CLL patients and normal controls; study CLL B cell LEF-1 protein expression; and investigate CLL B cell expression levels of other genes linked to Wnt signaling. To accomplish these goals, Affymetrix U133 GEP was performed on B cells from 41 B-CLL patients who were enrolled on a trial of combined pentostatin, cyclophosphamide, and rituximab and 11 healthy adults over age 60. Only the perfect match data were utilized for all analyses; the non-background corrected perfect match data were normalized using an intensity-dependent procedure and analyses were done using the base-2 logarithm transformed normalized values. Results: GEP analysis confirmed CLL B cells expressed the LEF-1 gene at ~28-fold higher levels than normal B cells (p<.0001). These results were validated by quantitative PCR and all CLL B cells studied expressed LEF-1 mRNA. By contrast, LEF-1 mRNA expression was undetectable in blood B cells from healthy adults both before and after in vitro mitogenic stimulation suggesting aberrant LEF-1 expression in B-CLL. Moreover, LEF-1 protein was readily detected by flow cytometry in CLL B cell samples and all clonal B cells expressed a uniform level of LEF-1. Full length LEF-1 is more oncogenic and predominates in certain cancers, while a dominant negative short isoform is found at greater levels in normal lymphocytes. Of interest, western blot analysis revealed that CLL B cells predominantly expressed full length oncogenic LEF-1. Ongoing studies are focused on silencing LEF-1 expression to determine LEF-1 target genes. In this regard, IGFBP4 (IGF binding protein 4) mRNA expression levels are 19-fold higher in CLL vs control B cells (p<0.0001) and the IGFBP4 protein was recently shown to antagonize Wnt signaling. A scan of the IGFBP4 promoter element identified 3 possible LEF/TCF consensus sites. Thus, there is a possible negative feedback loop in CLL, with activation of the Wnt pathway leading to expression of a negative Wnt regulator. Finally, we queried the GEP data for evidence of differential expression of 63 other Wnt-related genes. Surprisingly, in addition to LEF1, only five genes met our selection threshold of ≥1.5 fold change in expression between CLL and control B cells and p value <.001. These genes were CSNK1D (casein kinase, delta 1; 1.9-fold higher in control vs CLL; p<0.0001); CCND2 (cyclin D2; 2.4-fold higher in CLL vs control; p<0.0001); JUN (3.2-fold higher in control vs CLL); WNT3 (6.4-fold higher in CLL vs control; p<0.0001); and TCF4 (transcription factor 4, alias ITF-2; 4-fold higher in CLL vs control; p<0.0001). To our knowledge, we are the first to report that the TCF4 (ITF-2) gene is over expressed in B-CLL. This TF is of great interest because it is a known downstream target of the Wnt/TCF pathway, is activated in human cancers with b-catenin defects, and it promotes neoplastic transformation. Using Spearman correlation analysis to explore the relationship between LEF1, TCF4, and WNT3 expression levels, we observed a marginal correlation between WNT3 expression and LEF1 (ρ=0.31; p value=0.05) and TCF4 (ρ=0.29; p value=0.07). Summary: These studies add new support for a critical role of the Wnt pathway in CLL. The WNT3 results are consistent with the literature; however, it is notable that we failed to corroborate other reports demonstrating that other Wnts and Fzd genes are differentially expressed in CLL B cells vs normal B cells. It is striking that CLL B cells over express two TFs in this pathway that have been linked with neoplastic transformation. Ongoing studies are aimed at further elucidating the roles of these genes in B-CLL.

Abstract Lisavanbulin (LIS; BAL101553) is the prodrug of BAL27862, a microtubule-binding, tumor checkpoint controller and potential radiosensitizer. These studies evaluated optimal integration of LIS with standard of care radiation therapy (RT) and/or temozolomide (TMZ) using GBM PDX models. Distribution across the blood brain barrier was evaluated after a single 30 mg/kg oral LIS dose, and concentrations of the active metabolite BAL27862 were measured by liquid chromatography-tandem mass spectrometry. Similar BAL27862 concentrations were detected in the brain (B) and plasma (P) at both two (B:P ratio 1.29) and six hours (B:P ratio 1.64) post-dose. An in vivo screen of LIS monotherapy across 14 orthotopic GBM PDX models showed significant survival benefit (p<0.01) in seven models (median survival extension 24-87%). Extending from these results, LIS was evaluated in several of the sensitive models in combination with RT +/- TMZ. Two MGMT unmethylated PDXs, GBM6 and GBM150, were treated with vehicle or two weeks of RT +/- LIS. LIS dosing during the RT dosing period did not significantly improve median survival in either line (GBM6 survival with RT 54 days vs RT/LIS 58 days, p=0.16; GBM150 RT 86 days vs RT/LIS 101 days, p=0.21). However, prolonged LIS dosing from the start of RT until mice reached a moribund state demonstrated added benefit (GBM6 median 90 days vs RT 69 days, p=0.0001; GBM150 median 143 days vs RT 73 days, p=0.06). In GBM6, prolonged LIS dosing also significantly extended survival when combined with 2 weeks of RT/TMZ (median 101 days vs 66 days, p<0.0001), while LIS alone or RT/TMZ resulted in similar median survivals (63 days vs 66 days, respectively; p=0.68). This same RT/TMZ/LIS benefit was not seen in the MGMT methylated GBM12. Subsequent experiments were performed to evaluate integration of prolonged LIS dosing with concurrent RT/TMZ followed by 3 cycles of adjuvant TMZ (‘Stupp’ regimen). In MGMT methylated GBM39, LIS alone did not significantly extend survival, but LIS addition to the Stupp regimen doubled median survival (Stupp 249 days vs Stupp/LIS 502 days, p=0.0001). GBM150 demonstrated equal benefit from LIS alone or Stupp regimen (median 118 days vs 123 days, p=0.49). Stupp/LIS showed no additional survival benefit (median 98 days, p=0.97). In a second MGMT unmethylated, TMZ-resistant GBM26 PDX, LIS alone or combined with the Stupp regimen provided significant survival benefit: median survival 53 days for vehicle vs. 80 days for LIS (p=0.0001), 114 days for RT only (p<0.0001), 147 days for RT/LIS (p=0.30 relative to RT), 121 days for ‘Stupp’ regimen alone (p=0.57 relative to RT), and 172 days for Stupp/LIS (p=0.04 relative to Stupp). A follow-up GBM39 study revealed a significant increase in the mitotic marker phospho-histone H3 with LIS treatment relative to vehicle-treated controls (p=0.01) while Ki67 levels were similar (p=0.15). This suggests that LIS induces a mitotic arrest associated with microtubule deregulation. Collectively, these data provide a strong rationale to evaluate lisavanbulin (BAL101553) with RT +/- TMZ in GBM and provided the basis for an ongoing Phase I clinical trial. Citation Format: Danielle M Burgenske, Ann C Mladek, Jenny L Pokorny, Heidi A Lane, Felix Bachmann, Rachael A Vaubel, Mark A Schroeder, Katrina K Bakken, Lihong He, Zeng Hu, Brett L Carlson, Surabhi Talele, Gautham Gampa, Matthew L Kosel, Paul A Decker, Jeanette E Eckel-Passow, William F Elmquist, Jann Sarkaria. Modeling the clinical paradigm of lisavanbulin (BAL101553) deployment in patient-derived xenografts (PDX) of glioblastoma (GBM) [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference on Molecular Targets and Cancer Therapeutics; 2019 Oct 26-30; Boston, MA. Philadelphia (PA): AACR; Mol Cancer Ther 2019;18(12 Suppl):Abstract nr C096. doi:10.1158/1535-7163.TARG-19-C096

4567 Background: Trials combining anti-PDL1 and anti-VEGF therapies for clear cell renal cell carcinoma (ccRCC) are underway; however the relationship between the expression of PDL1 and genes in the VEGF pathway is poorly understood. Using the Affymetrix platform, we observed an inverse association between the expression of PDL1 and key genes in the VEGF pathway. Herein, we validate this inverse association using an independent set of 100 primary ccRCC tumors. Methods: From our registry database we sampled 100 ccRCC tumors with varying PDL1 protein expression (0%-100%) determined by immunohistochemistry (IHC). We extracted RNA from FFPE slides and performed RT-PCR to quantify gene expression of PDL1, VEGF, VEGFR1, and VEGFR2. All genes were normalized to the POLR2a gene. We evaluated the association of PDL1 protein expression and VEGF gene expression using Spearman rank correlation. In addition, we employed a linear mixed effects model to compare the fold-change in expression of VEGF genes between PDL1 low (0-5%, n=68) and PDL1 high (>5%, n=32) tumors. Results: As expected, PDL1 protein expression positively correlates with PDL1 geneexpression (corr=0.42, p<0.001). Validating our array data, PDL1 protein expression inversely correlates with expression of key VEGF genes: VEGF (corr=-0.23, p=0.01), VEGFR1 (corr=-0.34, p<0.001), and VEGFR2 (corr=-0.23, p=0.01). In our dichotomized analysis, we noted significantly higher expression of VEGF genes in the PDL1 low compared to the PDL1 high group: VEGF (fold change=1.82, p<0.001), VEGFR1 (FC=2.63, p<0.001), and VEGFR2 (FC=2.13, p=0.001). Conclusions: We independently validate an inverse association between the expression of PDL1 and key genes in the VEGF pathway in primary ccRCC. If validated further in larger studies, the existence of an immune evasive and an angiogenic phenotype within ccRCC could inform current clinical trials targeting these two pathways. Ultimately, whether VEGF signaling affects PDL1 expression, or whether angiogenic or immune evasive phenotypes predict response to anti-PDL1, anti-VEGF, or a combination of the two therapies remains unclear.

e15539 Background: Despite the development of prognostic algorithms based on clinico-pathologic features, the ability to identify aggressive forms of clear cell renal cell carcinoma (ccRCC) remains suboptimal. Topoisomerase IIA (TOP2a) is a biomarker of DNA replication and a target for antineoplastic agents. Herein, we evaluate the association of TOP2a expression in ccRCC tumors with pathologic features of aggressiveness and risk of cancer-specific death. Methods: We identified 1,380 patients who underwent nephrectomy to treat clinically localized ccRCC between 1/16/1990 and 4/14/2009. TOP2a expression was assessed using IHC and scored as number of positive cells per mm2. We evaluated TOP2a expression using a continuous variable and tertile categories. For associations with pathologic features we employed Kruskal-Wallis tests and for associations with cancer-specific survival we generated Cox proportional hazard regression models. Results: HigherTOP2a expression is associated with later stage, higher grade and higher MayoSSIGN score (all p < 0.001). The risk of death from RCC increases with increasing TOP2a expression (p trend < 0.0001), and this association remained strong after after multivariate adjustment for well-known predictors of RCC aggressiveness. Compared to patients in the lowest tertile, those patients with tumors in the highest tertile of TOP2a expression were at increased risk of RCC death (HR=2.62 95% CI 1.95-3.54; p<0.0001). Interestingly, among those patients with low risk disease (SSIGN score 0-3; ~95% 10 year survival), those with high TOP2a were at increased risk of RCC death (HR=3.48 95% CI 1.56-7.76; p =0.002). Conclusions: Higher TOP2a expression is associated with more aggressive pathologic features and increased risk of cancer-specific death among patients undergoing surgery for localized ccRCC. If confirmed, these data support further inquiry for TOP2a as a prognostic and predictive biomarker for ccRCC patients.

2072 Background: Human glioma stem cells and xenograft lines are common translational models in neuro-oncology but it has not been established if they are genetically and phenotypically comparable. This study aimed to determine if human glioma xenografts and stem cell-derived tumors had similar genotypes. Methods: Matched glioma stem cell cultures and subcutaneous xenograft lines were generated from four human glioblastoma specimens (BT114, BT116, BT120, BT132). Comparison was made between subcutaneous stem cell-derived tumors (flank) and xenografts established in nude mice. Copy number variation (CNV) and gene expression microarray studies were performed. Results: Various differences in copy number and gene expression were seen. Observed CNVs included regions within EGFR, myc, and p16 (INK). For example, EGFR copy number was two fold higher in xenografts vs. stem cell-derived tumor in one line (BT114). This difference was corroborated by western blot. Other differences included a heat shock protein homolog (DNAJA4), tetraspanin 13, and a p53 family target gene (ISG20L1). Two lines (BT114, BT116) had a greater than two fold increase in DNAJA4 expression in xenografts vs. stem cell-derived tumors (p = 0.04, 0.01). Two cell lines (BT116, BT120) had a two to eight fold increase in tetraspanin 13 expression in xenografts (p = 0.02, 0.05). However, neither copy number nor gene expression variations were consistent across all cell lines. Conclusions: Xenografts and glioma stem cell-derived tumors established from the same patient specimens have distinct genotypes. Further work is needed to establish if these differences are random or represent characteristic changes selected by different in vitro or in vivo pressures. However, these variations raise questions regarding which model is ideal for studying glioma biology, and which ones best replicate glioma characteristics in human patients.