Back to Authors

Tian, Shulan

https://orcid.org/0000-0002-3348-7439

Mayo Clinic

Mayo Clinic

Citations 52

Abstract 6466: Serum immunoglobulins are an independent prognostic marker of time to first therapy in newly diagnosed chronic lymphocytic leukemia (CLL) and monoclonal B-cell lymphocytosis (MBL)
Allmer, Cristine, Tian, Shulan , Campos, Cecília Bonolo de , Bruins, Laura A., O’Brien, Daniel , Kabat, Brian, Yan, Huihuang , Call, Timothy G. , Leis, José F., Kenderian, Saad J., Vachon, Celine M. , Parikh, Sameer A. , Kleinstern, Geffen , Boddicker, Nicholas , Rabe, Kari G. , Norman, Aaron D. , Li, Xing , Ding, Wei , Vachon, Celine M. , Cerhan, James R. , Kay, Neil E. , Slager, Susan L. , Braggio, Esteban

Cancer Research | September 4, 2020 | Volume: 80 | Issue: 16_Supplement | Pages: 6466-6466

10.1158/1538-7445.am2020-6466 Get It

Abstract Immune dysfunction is a hallmark of CLL/MBL; ~25% individuals with newly diagnosed CLL have hypogammaglobulinemia at diagnosis. In this study of newly diagnosed CLL/MBL, we sought to a) investigate the association between serum immunoglobulins and specific mutations identified by next generation sequencing, and b) determine if serum immunoglobulin levels predict TFT, independent of the CLL-International Prognostic Index (CLL-IPI) and tumor mutational load (TML). We used the Mayo Clinic CLL Resource to identify newly diagnosed CLL/MBL individuals consented within 2 years of diagnosis. Baseline clinical characteristics and the CLL-IPI score was ascertained on all individuals. DNA was extracted from peripheral blood CD5+/CD19+ clonal B-cells; and the entire coding regions of 60 somatically recurring mutated CLL driver genes were enriched using a customized Agilent SureSelect gene panel and sequenced in Illumina Hiseq4000. Serum immunoglobulins were quantitated by radial immunodiffusion using Immunoplates. The relationship between serum immunoglobulins was assessed: a) with individual gene mutations (mutated vs. not); and b) using the TML score (using a value of 0, 1, and 2 or greater [2+] mutated genes). The impact of serum immunoglobulin levels on TFT was analyzed for all individuals. Of 356 individuals identified, 80 (22%) had MBL, median age was 63 years (28-87), 261 (73%) were male, and the CLL-IPI distribution was: low risk (n=149, 42%), intermediate risk (n=112, 32%), and high/very high risk (n=94, 26%). We observed 43%, 39%, and 18% individuals with 0, 1, or 2+ mutated genes, respectively. The most commonly mutated genes were NOTCH1 (13%), SF3B1 (13%), TP53 (8%), ATM (8%), BIRC3 (7%), and NFKBIE (7%). Using standard cutoffs, 135 (38%) individuals had low serum IgG, 73 (21%) had low serum IgA, and 154 (44%) had low serum IgM at diagnosis. There were no differences in individual gene mutation frequencies between individuals with low and normal serum IgG. In contrast, individuals with low serum IgA were more likely to have mutations in NOTCH1 (25% vs 10%) and ATM (14% vs 7%) compared to those with normal IgA, and individuals with low serum IgM were more likely to have mutations in NOTCH1 (19% vs 9%) and BIRC3 (11% vs 4%) compared to those with normal IgM (all p-values <0.05). There were no significant differences in individuals with low serum IgG based on TML scores of 0, 1, and 2+ (31%, 35%, 40%, respectively, p=0.28). However, individuals with low serum IgA and low serum IgM were more likely to have a higher TML (8%, 20%, 36% for IgA, p<0.0001; and 35%, 45%, 54% for IgM, p<0.0001, for TML scores of 0, 1, and 2+, respectively). After a median follow-up of 6.9 years, 157 individuals required CLL therapy. After adjusting for CLL-IPI and TML, low serum IgG, low serum IgA, and low serum IgM (hazard ratio [HR]: 1.8, 2.3, and 2.4, respectively, all p-values <0.05) were independently associated with a shorter TFT. Our study shows that individuals with CLL/MBL who harbor mutations in NOTCH1, ATM and BIRC3 are more likely to have low serum IgA and IgM at diagnosis. Serum IgA and IgM (but not serum IgG) levels decrease with an increase in the TML score. Finally, low serum immunoglobulins can predict TFT, independent of the CLL-IPI and TML, which has important implications for risk stratification in newly diagnosed CLL/MBL. Citation Format: Sameer A. Parikh, Cristine Allmer, Geffen Kleinstern, Nicholas J. Boddicker, Shulan Tian, Cecilia B. De Campos, Laura A. Bruins, Daniel R. o"brien, Brian F. Kabat, Kari G. Rabe, Aaron D. Norman, Huihuang Yan, Xing Li, Timothy G. Call, Wei Ding, Jose F. Leis, Saad J. Kenderian, Celine M. Vachon, James R. Cerhan, Neil E. Kay, Susan L. Slager, Esteban Braggio. Serum immunoglobulins are an independent prognostic marker of time to first therapy in newly diagnosed chronic lymphocytic leukemia (CLL) and monoclonal B-cell lymphocytosis (MBL) [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 6466.
Abstract 1256: Widespread alterations of chromatin accessibility in chronic lymphocytic leukemia (CLL)
Yan, Huihuang , Tian, Shulan , Kleinstern, Geffen , Li, Xing , Wang, Zhiquan , Boddicker, Nicholas , Slager, Susan L. , Braggio, Esteban

Cancer Research | September 4, 2020 | Volume: 80 | Issue: 16_Supplement | Pages: 1256-1256

10.1158/1538-7445.am2020-1256 Get It

Abstract Chronic lymphocytic leukemia (CLL) is the most common adult leukemia in western countries that has a strong genetic basis. Using data from both chromatin immunoprecipitation and sequencing (ChIP-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq), recent studies have uncovered alterations of chromatin landscape in CLL. However, the spectrum of epigenetic changes and the mechanisms of casualtyare only partially understood. In this study, we used ATAC-seq to profile chromatin accessibility across the two major CLL subtypes, IGHV-mutated (n=19) and IGHV-unmutated (n=22). DNA library of accessible chromatin regions was generated from purified CD5+/CD19+ leukemiacells. All libraries were assessed to have high data quality based on the ENCODE quality metrics. Additionally, publicly available ATAC-seq data (GSE118189 and GSE111015) from total (n=5), memory (n=6) and naïve B cells (n=6) were analyzed in parallel and used as normal reference. ATAC-seq peaks from CLL and normal B cells were merged to generate a consensus list of 107,928 high-confidence peaks that did not overlap ENCODE blacklisted regions, were present in at least two CLL or two normal B-cell samples, and had ≥1 counts-per-million in at least 10% of the samples. Raw read counts were Trimmed Mean of M-values normalized and log2 transformed. K-means clustering of differential peaks revealed an extensive reprogramming of chromatin accessibility landscape in CLL, with over 11,000 lost and nearly 9,000 gained open chromatin regions relative to normal B-cells. In particular, we identified two clusters comprising 1,401 regions that were largely open either in IGHV-mutated(571 peaks) or in IGHV-unmutated CLL (830 peaks). Open chromatin regions gained in IGHV-mutatedwere preferentially linked to genes in the B cell receptor signaling pathway and those in IGHV-mutatedto genes enriched in the interferon type I signaling pathway. These discriminatory regions are highly enriched with DNA binding motifs for transcription factors (TFs) with key roles in B cell development, such as PU.1, AP-1, c-Jun, BATF, and TCF family. In addition, 94% of the 1,401 peaks were located in distal regulatory regions, being ≥2.5 kb away from transcription start sites of genes; two-third of these distal regulatory regions overlapped H3K27ac peaks from CLL patients or immune cell types in the Roadmap reference epigenomes,indicating that enhancers are the prime sites of epigenetic differentiation between the two subtypes. Together, our analyses provide new insights into the global alterations of regulatory regions in CLL. Given the noticeable heterogeneity of chromatin landscape, further integration with histone modifications, TF binding, chromatin interaction, gene expression and genetic variation data are needed to understand the roles of epigenetic dysregulation and its interaction with genetic lesions in CLL pathogenesis. Citation Format: Huihuang Yan, Zhiquan Wang, Shulan Tian, Geffen Kleinstern, Nicholas Boddicker, Xing Li, Susan Slager, Esteban Braggio. Widespread alterations of chromatin accessibility in chronic lymphocytic leukemia (CLL) [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 1256.
Abstract 1113: Evaluation of novel therapeutics using models from the Mayo Clinic GBM patient-derived xenograft (PDX) collection
Tian, Shulan , Hu, Zeng, Carlson, Brett L. , Decker, Paul A. , Eckel‐Passow, Jeanette E. , Burgenske, Danielle M. , Mladek, Ann C. , Vaubel, Rachael A. , Schroeder, Mark A. , Decker, Paul A. , Sarkaria, Jann N.

Cancer Research | September 4, 2020 | Volume: 80 | Issue: 16_Supplement | Pages: 1113-1113

10.1158/1538-7445.am2020-1113 Get It

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.
Chronic lymphocytic leukemia (CLL) risk is mediated by multiple enhancer variants within CLL risk loci
Yan, Huihuang , Tian, Shulan , Kleinstern, Geffen , Wang, Zhiquan , Lee, Jeong‐Heon , Cerhan, James R. , Braggio, Esteban , Boddicker, Nicholas J, Cerhan, James R , Lee, Jeong-Heon , Braggio, Esteban , Boddicker, Nicholas , Kay, Neil E. , Slager, Susan L. , Wang, Zhiquan

Human Molecular Genetics | August 3, 2020 | Volume: 29 | Issue: 16 | Pages: 2761-2774

10.1093/hmg/ddaa165 Get It

Abstract Chronic lymphocytic leukemia (CLL) is the most common adult leukemia in Western countries. It has a strong genetic basis, showing a ~ 8-fold increased risk of CLL in first-degree relatives. Genome-wide association studies (GWAS) have identified 41 risk variants across 41 loci. However, for a majority of the loci, the functional variants and the mechanisms underlying their causal roles remain undefined. Here, we examined the genetic and epigenetic features associated with 12 index variants, along with any correlated (r2 ≥ 0.5) variants, at the CLL risk loci located outside of gene promoters. Based on publicly available ChIP-seq and chromatin accessibility data as well as our own ChIP-seq data from CLL patients, we identified six candidate functional variants at six loci and at least two candidate functional variants at each of the remaining six loci. The functional variants are predominantly located within enhancers or super-enhancers, including bi-directionally transcribed enhancers, which are often restricted to immune cell types. Furthermore, we found that, at 78% of the functional variants, the alternative alleles altered the transcription factor binding motifs or histone modifications, indicating the involvement of these variants in the change of local chromatin state. Finally, the enhancers carrying functional variants physically interacted with genes enriched in the type I interferon signaling pathway, apoptosis, or TP53 network that are known to play key roles in CLL. These results support the regulatory roles for inherited noncoding variants in the pathogenesis of CLL.
Tumor mutational load predicts time to first treatment in chronic lymphocytic leukemia ( CLL ) and monoclonal B‐cell lymphocytosis beyond the CLL internat
Kleinstern, Geffen , O’Brien, Daniel R. , Tian, Shulan , Kabat, Brian, Rabe, Kari G. , Yan, Huihuang , Vachon, Celine M. , Call, Timothy G. , Bruins, Laura A., Campos, Cecília Bonolo de , Leis, José F., Shanafelt, Tait D. , Ding, Wei , Braggio, Esteban , Kabat, Brian F, Chaffee, Kari G. , Norman, Aaron D. , Boddicker, Nicholas J, Parikh, Sameer A , Bruins, Laura, Bonolo de Campos, Cecilia, Ding, Wei , Cerhan, James R , Shanafelt, Tait D , Braggio, Esteban , Li, Xing , Vachon, Celine M. , Boddicker, Nicholas , Parikh, Sameer A. , Cerhan, James R. , Kay, Neil E. , Slager, Susan L.

American Journal of Hematology | April 12, 2020 | Volume: 95 | Issue: 8 | Pages: 906-917

10.1002/ajh.25831 Get It

Abstract Next‐generation sequencing identified about 60 genes recurrently mutated in chronic lymphocytic leukemia (CLL). We examined the additive prognostic value of the total number of recurrently mutated CLL genes (i.e., tumor mutational load [TML]) or the individually mutated genes beyond the CLL international prognostic index (CLL‐IPI) in newly diagnosed CLL and high‐count monoclonal B‐cell lymphocytosis (HC MBL). We sequenced 59 genes among 557 individuals (112 HC MBL/445 CLL) in a multi‐stage design, to estimate hazard ratios (HR) and 95% confidence intervals (CI) for time‐to‐first treatment (TTT), adjusted for CLL‐IPI and sex. TML was associated with shorter TTT in the discovery and validation cohorts, with a combined estimate of continuous HR = 1.27 (CI:1.17‐1.39, P = 2.6 × 10 −8 ; c ‐statistic = 0.76). When stratified by CLL‐IPI, the association of TML with TTT was stronger and validated within low/intermediate risk (combined HR = 1.54, CI:1.37‐1.72, P = 7.0 × 10 −14 ). Overall, 80% of low/intermediate CLL‐IPI cases with two or more mutated genes progressed to require therapy within 5 years, compared to 24% among those without mutations. TML was also associated with shorter TTT in the HC MBL cohort (HR = 1.53, CI:1.12‐2.07, P = .007; c ‐statistic = 0.71). TML is a strong prognostic factor for TTT independent of CLL‐IPI, especially among low/intermediate CLL‐IPI risk, and a better predictor than any single gene. Mutational screening at early stages may improve risk stratification and better predict TTT.
HLA class-I and class-II restricted neoantigen loads predict overall survival in breast cancer
Ren, Yingxue , Cherukuri, Yesesri , Wickland, Daniel P. , Sarangi, Vivekananda , Tian, Shulan , Carter, Jodi M. , Mansfield, Aaron S. , Block, Matthew S. , Knutson, Keith L. , Lin, Yi , Asmann, Yan W. , Mansfield, Aaron S , Lin, Yi , Carter, Jodi , Sherman, Mark E. , Knutson, Keith L.

OncoImmunology | April 1, 2020 | Volume: 9 | Issue: 1 | Pages: No page available

10.1080/2162402x.2020.1744947 Get It

Tumor Mutational Load and Germline Polygenic Risk Score Predicts Time-to-First Treatment in Chronic Lymphocytic Leukemia (CLL) and High-Count Monoclonal B Cell Lymphocytosis (MBL)
Kleinstern, Geffen , O’Brien, Daniel R. , Kabat, Brian, Yan, Huihuang , Vachon, Celine M. , Tian, Shulan , Call, Timothy G. , Leis, José F., Rabe, Kari G. , Norman, Aaron D. , Vachon, Celine M. , Parikh, Sameer A. , Ding, Wei , Cerhan, James R. , Kay, Neil E. , Braggio, Esteban , Slager, Susan L.

Blood | February 14, 2020 | Volume: 134 | Issue: Supplement_1 | Pages: 852-852

10.1182/blood-2019-124346 Get It

Background The CLL international prognostic index (CLL-IPI) is a validated staging system for newly diagnosed CLL patients that combines age, clinical stage, serum Beta-2 microglobulin, TP53, and IGHV mutation status into a single score and stratifies CLL patients into four prognostic risk groups with a c-statistic=0.72. Next-generation sequencing identified ~60 putative driver genes recurrently mutated in CLL, and genome-wide association studies (GWAS) identified 41 single nucleotide polymorphisms (SNPs) associated with CLL risk. We previously showed that a polygenic risk score (PRS) of the weighted average of the number of risk alleles among these 41 SNPs is associated with risk of both CLL and MBL. We and others have also shown that the tumor mutational load (TML), i.e., the total number of somatically mutated driver genes, is associated with time to first treatment (TTT). Herein, we examine whether TML and PRS add prognostic value beyond CLL-IPI for disease progression in a cohort of CLL and high-count MBLs. Methods We used the Mayo Clinic CLL resource to identify newly diagnosed (<2 years) previously untreated CLL and high-count MBL patients based on the 2008 International Workshop CLL (IWCLL) criteria . We extracted DNA from CD5+/CD19+ clonal B-cells and sequenced the coding regions of 59 recurrently mutated CLL driver genes using a custom SureSelect panel and Illumina HiSeq 4000. The median coverage depth per sample was 2737X. Somatic mutations were called using MuTect2 in tumor-only mode. Germline variants were eliminated based on minor allele frequencies >0.01%, and by those identified in 1000 Genomes Project, ExAC and/or ESP6500 databases, unless present in known mutation hotspots or COSMIC. After filtering high/moderate impact mutations, we calculated the TML. Peripheral blood DNA was genotyped on Illumina Omni Express array; after extensive quality control metrics, we calculated the germline PRS as the weighted average of the number of risk alleles among the CLL susceptibility SNPs. Cox regression was used to estimate hazard ratios (HR) and 95% confidence intervals (CI) to test for associations with TTT. We analyzed the low/intermediate risk groups and the high/very high risk groups in stratified analyses and tested for interactions between CLL-IPI and TML or PRS. Survival curves are displayed using the Kaplan-Meier method using P-values from the log-rank test; the PRS is displayed by tertile categories. Results Of 166 patients (120 CLL, 46 MBL), 67% were male and the median age at diagnosis was 62 years (range: 28-85). The median follow-up was 5.3 years, and 73 were treated for progressive disease as defined by the 2018 IWCLL criteria. Low risk CLL-IPI comprised 51%, intermediate risk 27%, high risk 18%, and very-high risk 4%. The most commonly somatically mutated genes were TP53 (8%), CHD2 (8%), BIRC3 (8%), SF3B1 (7%), and NOTCH1 (7%). We found that 42%, 32%, and 26% of patients had 0, 1, or 2+ genes with somatic mutations, respectively. The median PRS was 8.3 (range: 6.1-11.3) which is consistent with our prior PRS study on CLL/MBL and higher than 7.53 found in 7983 controls from our prior work. When we modeled TML, PRS, sex, CLL/MBL status and CLL-IPI together with the joint effects of CLL-IPI with TML and PRS, both interactions were significant in the model (P=0.008 for TML*CLL-IPI and P=0.014 for PRS*CLL-IPI), with c-statistic=0.81 (CI:0.77-0.85). When we stratified by CLL-IPI, TML was associated with shorter TTT in the low/intermediate risk (N=129, continuous HR=2.36 per TML category, CI: 1.61-3.47, P=1.3x10-5, Figure 1.A), but there was no evidence of an association in the high/very-high risk (N=37, continuous HR=1.12 per TML category, CI: 0.70-1. 78, P=0.64, Figure 1.B), after adjusting for PRS, sex, and CLL/MBL status. However, the PRS was associated with a shorter TTT in the high/very-high risk (continuous HR=1.77 per SD, CI: 1.00-3.12, P=0.048, Figure 1.D), with no evidence of an association in the low/intermediate risk (continuous HR=0.95 per SD, CI: 0.95-1.29, P=0.74, Figure 1.C), after adjusting for TML, sex, and CLL/MBL status. Conclusions In our study, TML was prognostic among the low/intermediate CLL-IPI risk groups while the PRS was prognostic among the high/very-high CLL-IPI risk groups. Sequencing CLL driver genes and genotyping relevant germline SNPs at time of diagnosis may provide incremental value over the CLL-IPI in predicting TTT among patients with CLL/MBL. Disclosures Parikh: MorphoSys: Research Funding; Genentech: Honoraria; AbbVie: Honoraria, Research Funding; AstraZeneca: Honoraria, Research Funding; Pharmacyclics: Honoraria, Research Funding; Janssen: Research Funding; Ascentage Pharma: Research Funding; Acerta Pharma: Research Funding. Ding:Merck: Research Funding; DTRM Biopharma: Research Funding. Cerhan:Celgene: Research Funding; Janssen: Membership on an entity's Board of Directors or advisory committees; NanoString: Research Funding. Kay:MorphoSys: Other: Data Safety Monitoring Board; Infinity Pharmaceuticals: Other: DSMB; Celgene: Other: Data Safety Monitoring Board; Agios: Other: DSMB.
Integrative Analysis of FISH, Transcriptomics and Mutational Status Predicts Responsiveness to Novel Agents in Multiple Myeloma
Campos, Cecília Bonolo de , Meurice, Nathalie , Bruins, Laura A., Polito, Alysia, Ahmann, Gregory , Tian, Shulan , Stewart, A. Keith, Stein, Caleb K. , Petit, Joachim L. , Bergsagel, P. Leif , Braggio, Esteban , Stewart, A. Keith

Blood | February 14, 2020 | Volume: 134 | Issue: Supplement_1 | Pages: 574-574

10.1182/blood-2019-130977 Get It

Introduction Despite continuous improvement of clinical outcome in multiple myeloma (MM), disease relapse remains a major challenge, leading to progressively shorter remissions and fewer treatment options. Strategies attempting to counteract this challenge include recent efforts resulting in an increase in the availability of novel promising anti-MM agents and targeting specific genetic profiles of the disease. In this context, we aim to develop predictive models of sensitivity and resistance to novel compounds by connecting an ex vivo high-throughput drug screen with genetic, transcriptomics, FISH, and clinical features. Methods Twenty compounds (afatinib, afuresertib, belinostat, buparlisib, cobimetinib, CPI-0610, crenolanib, dinaciclib, dovitinib, JQ1, LGH447, osimertinib, OTX015, panobinostat, romidepsin, selinexor, sunitinib, trametinib, venetoclax, and vorinostat) were selected based on overall promising anti-MM activity from an ex vivo high throughput drug screen with a panel of 79 single agents incubated for 24 hours. The area under the curve (AUC) was used to rank order the ex vivo responses for each compound and the lowest and highest quartile samples were identified for further analysis. Clinical data and FISH data, including t(11;14), t(4;14), t(14;16), del(17p), +1q, monosomy 13, and MYC rearrangement, were collected. Targeted DNA sequencing was performed using a 2.3 Mb custom capture panel covering 139 MM-relevant genes. mRNA-sequencing was performed and differential gene expression analysis in the highest and lowest quartile identified subsets of markers positively and negatively associated with the AUC response for a given compound. An additional unbiased selection of markers using lasso techniques was performed, resulting in predictive generalized linear models (GLM) for each agent. Responses from the remaining intermediate samples were estimated with the predictive models, with overall predictive ability assessed by correlating predicted AUCs with their actual counterparts. Results Our integrative analysis was performed on 50 primary patient samples (36% untreated and 64% relapsed MM). Venetoclax, dinaciclib, romidepsin, panobinostat, osimertinib, belinostat and selinexor were the most active compounds in the cohort. Interestingly, LGH447, dovitinib, selinexor, JQ1, OTX-015, cobimetinib, and trametinib showed increased activity in relapsed MM when compared to untreated samples (Wilcoxon Test; p<0.05). We generated GLMs using an average of 92 markers (range 64-107) per compound, combining mRNA-sequencing expression with FISH and mutation data. The analysis proposed in the present study was validated through the unbiased selection of BCL2 among the subset of markers included in the GLM predicting sensitivity to venetoclax, a first-in-class orally bioavailable selective BCL2 inhibitor. Expression level of critical NF-kB and cell cycle genes, such as BIRC3, CKS1B, PAX5, NFKB2, and CCND2, were included in 60% of our predictive models. Mutations of DNA repair genes (ATM,TP53) were included in the GLMs of three epigenetic therapies, one histone deacetylase inhibitor and two BET inhibitors, associated to ex vivo resistance to the drugs. The presence of monosomy 13 was also a marker for ex vivo resistance for five epigenetic therapies, four HDAC inhibitors and one BET inhibitor. The three BET inhibitors, JQ1, CPI-0610, and OTX015, were among the compounds most accurately predicted by our integrative approach, with Spearman correlation values between 0.773-0.858. Overall, our models accurately predicted the ex vivo response for 16 (80%) of the compounds (r>0.7). Five (25%) of these compounds displayed a remarkably accurate prediction model in both training (highest and lowest quartiles) and validation (intermediate quartiles) samples (r>0.8). Conclusions The GLM data integration approach enabled the establishment of effective predictive models, identifying FISH, transcriptomics, and mutations of putative driver genes important in anti-MM agent responsiveness. In addition, the resulting dataset is promising for future research focusing on the discovery of novel mechanisms of action and establishing markers of sensitivity and resistance to novel compounds. We are currently increasing our dataset and seek to create an omnibus approach that predicts responses to multiple anti-MM agents simultaneously. Disclosures Bergsagel: Celgene: Consultancy; Ionis Pharmaceuticals: Consultancy; Janssen Pharmaceuticals: Consultancy. Stewart:Amgen: Consultancy, Research Funding; Bristol Myers-Squibb: Consultancy; Celgene: Consultancy, Research Funding; Ionis: Consultancy; Janssen: Consultancy, Research Funding; Oncopeptides: Consultancy; Ono: Consultancy; Roche: Consultancy; Seattle Genetics: Consultancy; Takeda: Consultancy.
Genomic and Phenotypic Characterization of a Broad Panel of Patient-Derived Xenografts Reflects the Diversity of Glioblastoma
Vaubel, Rachael A. , Tian, Shulan , Remonde, Dioval , Mladek, Ann C. , Kitange, Gaspar J. , Caron, Alissa, Kollmeyer, Thomas M. , Grove, Rebecca, Peng, Sen, Carlson, Brett L. , Sarkar, Gobinda , Evers, Lisa, Decker, Paul A. , Yan, Huihuang , Berens, Michael E. , Wang, Qianghu , Marin, Bianca M. , Klee, Eric W. , Eckel‐Passow, Jeanette E. , Verhaak, Roel G.W. , Sulman, Erik P. , Meyer, Fredrick B., Tran, Nhan L. , Jenkins, Robert B. , Schroeder, Mark A. , Sarkar, Gobinda , Berens, Michael E , Wang, Qianghu , Marin, Bianca M , Califano, Andrea , Sulman, Erik P , O'Neill, Brian P, Parney, Ian F , Kollmeyer, Thomas , Vaubel, Rachael A , Mladek, Ann C , Tran, Nhan L , Grove, Rebecca, Daniel, J. , Decker, Paul A. , Dhruv, Harshil , Lachance, Daniel Honore , Burns, Terry C. , O’Neill, Brian Patrick , Giannini, Caterina , Parney, Ian F. , Sarkaria, Jann N.

Clinical Cancer Research | December 18, 2019 | Volume: 26 | Issue: 5 | Pages: 1094-1104

10.1158/1078-0432.ccr-19-0909 Get It

Abstract Purpose: Glioblastoma is the most frequent and lethal primary brain tumor. Development of novel therapies relies on the availability of relevant preclinical models. We have established a panel of 96 glioblastoma patient-derived xenografts (PDX) and undertaken its genomic and phenotypic characterization. Experimental Design: PDXs were established from glioblastoma, IDH-wildtype (n = 93), glioblastoma, IDH-mutant (n = 2), diffuse midline glioma, H3 K27M-mutant (n = 1), and both primary (n = 60) and recurrent (n = 34) tumors. Tumor growth rates, histopathology, and treatment response were characterized. Integrated molecular profiling was performed by whole-exome sequencing (WES, n = 83), RNA-sequencing (n = 68), and genome-wide methylation profiling (n = 76). WES data from 24 patient tumors was compared with derivative models. Results: PDXs recapitulate many key phenotypic and molecular features of patient tumors. Orthotopic PDXs show characteristic tumor morphology and invasion patterns, but largely lack microvascular proliferation and necrosis. PDXs capture common and rare molecular drivers, including alterations of TERT, EGFR, PTEN, TP53, BRAF, and IDH1, most at frequencies comparable with human glioblastoma. However, PDGFRA amplification was absent. RNA-sequencing and genome-wide methylation profiling demonstrated broad representation of glioblastoma molecular subtypes. MGMT promoter methylation correlated with increased survival in response to temozolomide. WES of 24 matched patient tumors showed preservation of most genetic driver alterations, including EGFR amplification. However, in four patient–PDX pairs, driver alterations were gained or lost on engraftment, consistent with clonal selection. Conclusions: Our PDX panel captures the molecular heterogeneity of glioblastoma and recapitulates many salient genetic and phenotypic features. All models and genomic data are openly available to investigators.
Association between a Polygenic Risk Score for Multiple Myeloma Risk and Overall Survival
Clay‐Gilmour, Alyssa , Asmann, Yan W. , Brown, Elizabeth E. , Hofmann, Jonathan N, Spinelli, John J. , Robinson, Dennis P. , O’Brien, Daniel R. , Tian, Shulan , Berndt, Sonja I. , Machiela, Mitchell J. , Wu, Xifeng , Waller, Rosalie, Kumar, Shaji , Camp, Nicola J. , Vachon, Celine M. , Machiela, Mitchell J , Hildebrandt, Michelle A.T. , Hofmann, Jonathan N, Giles, Graham G. , Bhatti, Parveen , Cozen, Wendy , Rajkumar, S. Vincent , Chanock, Stephen Jacob , Norman, Aaron D. , Sinnwell, Jason P. , Milne, Roger L. , Slager, Susan L. , Ziv, Elad , Vachon, Celine M.

Blood | November 14, 2019 | Volume: 134 | Issue: Supplement_1 | Pages: 4366-4366

10.1182/blood-2019-126088 Get It

Background Genome-wide association studies (GWAS) conducted in populations of European ancestry (EA) have identified and confirmed 23 germline susceptibility loci for multiple myeloma (MM). The effect sizes of single nucleotide polymorphisms (SNPs) at these loci are small, therefore combining them into a single summary measure, known as a polygenic risk score (PRS), may provide a more meaningful risk factor. We have previously shown a PRS comprised of the 23 SNPs for MM contributes to increased risk of MM, with a 2.7-fold increase for highest vs. lowest PRS quintiles. Whether the MM-PRS is also associated with overall survival (OS) in MM cases has not been evaluated. We examined the association between MM-PRS and OS in two EA studies. Methods The first study consisted of 2,179 EA MM cases from ten studies included in the Multiple Myeloma Working Group within the International Lymphoma Consortium (InterLymph). Cases were diagnosed between 1970 and 2015 and genotyped using multiple platforms (Oncoarray, Affymetrix, Human660W-quad Beadchip, and Illumina arrays); 885 cases also had stage [based on International Staging System (ISS)] available. Each of the GWAS was subjected to rigorous standard quality control independently (prior to imputation via the Michigan imputation server based on the Haplotype Reference Consortium (HRC). The second study consisted of 515 newly diagnosed EA MM cases from CoMMpass (Relating Clinical Outcomes in Multiple Myeloma to Personal Assessment of Genetic Profile), diagnosed from 2011-2013, who had whole genome sequencing (WGS) performed on germline DNA. The WGS data was used to call common germline genetic variants through the Mayo Clinic bioinformatics pipeline. Briefly, genetic variants were detected with GenomeGPS, aligned to the hg19 reference genome, called using the GATK (V3.6) Haplotype Caller, and merged for multiple-sample joint calling. To reduce the false positive variants, variant quality score recalibration (VQSR) was applied for both SNPs and INDELs. After quality control, 458 EA samples remained. Follow-up was available for both studies and consisted of time from MM diagnosis date until death or date of last known follow-up. The PRS was constructed from the 23 MM SNPs using the published per allele odds ratio associated with MM risk. The published log odds ratios for each SNP were multiplied by the number of risk alleles (0, 1, 2) for the corresponding SNP, and summed, resulting in a unique score per person. Kaplan-Meier curves and Cox proportional hazard models were used to assess the association between PRS with MM OS considering two models: 1) adjusted for age, sex, study and 2) additional adjustment by stage (ISS). Hazard ratios (OR) and 95% confidence intervals (CI) were estimated. The PRS was evaluated both as a continuous variable, per standard deviation (SD), and as a categorical variable (quintiles). Results MM cases (N=2,179) in the InterLymph study were 59% male and 41% female and the median age was 61.0 years (26-90 years). Median follow-up time was 57.2 months (1.0-509.0 months) with 868 reported deaths. MM cases with stage information available consisted of 20% stage I (n=178), 53% stage II (n=466), and 27% stage III (n=241). No association was observed between PRS and OS in MM patients regardless of adjustment for stage (continuous PRS (HR: 1.03, 95% CI: 0.83-1.28, P=0.80) or by quintile PRS (p>0.05)) (Table). The CoMMpass EA MM cases (n=458) had similar distributions for sex (61% male and 39% females) but were slightly older 65 years (27-93 years) and had shorter follow-up time (median=39.75 months (0.13-77.2)) with 117 deaths. Stage was available for 96% of CoMMpass cases including 36% stage I (n=159), 33% stage II (n=146), and 31% stage III (n=134). We also observed no association of PRS and OS in the CoMMpass study (HR=1.02, 95% CI: 0.72 -1.46, P= 0.89), adjusted for age, sex, and stage (Table). Discussion A PRS score for MM risk is not associated with OS for MM cases in two EA populations. Given that prior studies have shown association of genetic variation with MM survival, efforts to identify additional loci associated with OS or MM specific survival are warranted. Future studies should also consider germline variants impact on molecular subtypes, specific therapies, and outcomes. Disclosures Kumar: Celgene: Consultancy, Research Funding; Janssen: Consultancy, Research Funding; Takeda: Research Funding.