Citations

Abstract Introduction: Axicabtagene ciloleucel (axi-cel) is an autologous anti-CD19 chimeric antigen receptor (CAR) T cell therapy. ZUMA-1 is a multicenter registrational trial of axi-cel in patients with refractory, aggressive B-cell non-Hodgkin lymphoma. In a pre-specified interim analysis, ZUMA-1 met its primary endpoint with a 76% objective response rate and 47% complete response. The incidence of Grade 3 and higher (Gr 3+) cytokine release syndrome (CRS) and neurologic events (NE) was 13% and 29%, respectively, in 93 patients with 1 month follow-up (Neelapu ASH 2016). We present novel immune signatures of Gr 3+ CRS and NE by analyzing the axi-cel-related biomarker profile in association with clinical outcomes. Methods: In this interim analysis of 62 patients, 44 serum analytes pre- and post-axi-cel treatment were measured via ELISA at multiple timepoints during the first month. The number of CAR+ cells in blood was determined by qPCR. Kinetics and association of these markers with CRS and NE were analyzed. The treatment-related profile was defined by analytes with at least double an increase over baseline in at least 50% of patients. These analytes were selected and evaluated for association with Gr 3+ CRS or NE. Immune signatures were determined based on p-values applied to peak or cumulative levels of analytes and adjusted for multiplicity (stepdown method). Results: Out of 44 analytes measured, 12 demonstrated an increase by double over baseline in at least 50% of patients. This axi-cel-related biomarker profile includes immune proliferative/modulating cytokines, pro-inflammatory cytokines, markers of myeloid activation, and chemokines. Markers with the highest median-fold increase were IFNγ (44x), IL-10 (31x), IL-6 (26x), IL-15 (20x), and GRZB (17x). These peaked within 7-14 days and generally returned to baseline within 1 month post-treatment. Within this panel of 12 analytes elevated post-axi-cel treatment, markers associated with Gr 3+ CRS (p<0.05) were IL-10, IL-6, IL-2Rα, IL-8, MCP-1, IFNγ, IL-1Rα, and IP-10. Markers associated with Gr 3+ NE (p<0.05) were IL-6, IL-15, IL-8, and MCP-1. Peak levels of CAR T cells in blood, occurring within 14 days post-infusion, were associated with Gr 3+ NE (p=0.02) but not Gr 3+ CRS (p=0.6). A majority of patients with Gr 3+ NE had only Gr 0-2 CRS (16/22) while nearly 50% of patients with Gr 3+ CRS had Gr 0-2 NE (5/11). Conclusions: Biomarker and clinical evaluation suggests mechanistic differences underlying CRS and NE following axi-cel treatment. While Gr 3+ NE was associated with CAR T cell expansion, Gr 3+ CRS was instead associated with a broader array of serum markers, many linked to systemic myeloid activation. Further understanding of these immune signatures may provide insight into managing AEs associated with anti-CD19 CAR T cell therapy. Citation Format: Frederick L. Locke, John Rossi, Xiaodong Xue, Sattva S. Neelapu, Daniel H. Ryan, Armin Ghobadi, Lazaros J. Lekakis, David Miklos, Caron A. Jacobson, Ira Braunschweig, Olalekan Oluwole, Tanya Siddiqi, Yi Lin, John Timmerman, Patrick M. Reagan, Marika Sherman, Janice Nagatani, Xiao Zhang, Lynn Navale, William Y. Go, Jeff Wiezorek, Adrian Bot. Immune signatures of cytokine release syndrome and neurologic events in a multicenter registrational trial (ZUMA-1) in subjects with refractory diffuse large B cell lymphoma treated with axicabtagene ciloleucel (KTE-C19) [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 CT020. doi:10.1158/1538-7445.AM2017-CT020

Abstract Lung cancer is the most lethal malignancy around the world. In Hong Kong, chronic tobacco use is prevalent in male but <20% of female patients smoke. Epidermal growth factor receptor (EGFR) mutations affect 30% and 67% of male and female lung adenocarcinomas, respectively. Treatment responses to tyrosine kinase inhibitors (TKI) are not uniform and all patients eventually develop drug resistance. Mutation screening of immediate pre-treatment tumors is the most ideal for assessing drug targets but in practice, excision specimens of primary or metastatic sites are more efficient in terms of tumor quantity, timeliness and patient acceptance. However, information on the applicability of this approach and correlation with the clinical outcome is limited. For this purpose, we compared the whole exome mutation profiles of 39 excised EGFR mutant cancers treated with first line TKI with the clinical response type, including 15 stage I/II and 24 stage III/IV tumors. They included 16 non-responders (NR) defined clinically as those with stable or progressive disease, and 23 responders (R) defined by partial or complete tumor shrinkage. Considering only the non-synonymous SNV and INDEL mutations of coding regions of known actionable targets and genes listed in both the CGC and COSMIC cancer databases, the tumors harbored 54 mutated genes including 26 recurrently mutated and 28 genes involving only 1 tumor. Excluding EGFR, TP53 was the most commonly mutated gene occurring in 25/39 (64.1%) cases. Most known resistant genes involved in the EGFR and bypass signaling network showed mutations only in the NR group including recurrent mutations of PTEN (3/16), PIK3CA (2/16) and NF1 (2/16), and single case mutations of AKT1, ALK, RAF1 and KDR. Two EGFR network candidates showed mutations in both the NR and R groups, including HGF (2/16 NR, 1/23 R) and ROS1 (1 case in either group). Four post-treatment (post-TKI) tumors of acquired resistance were also analyzed, all of which showed EGFR T790M while no pre-treatment NR or R tumor harbored this mutation. Notably, mutations in the β-catenin pathway were prominent in the NR and post-TKI tumors, including APC (2/16 NR), CTNNB1 (1/16 NR, 1/4 post-TKI) and c-MYC (1/16 NR, 1/4 post-TKI) while they were not detected in the R tumors. Also, nonsense mutations of ARID1A was observed in 2/16 NR but none of the R tumors. In summary, our findings revealed candidate TKI resistant mutations involving the EGFR and bypass signaling networks in pre-treatment excision specimens particularly in non-responding patients. The EGFR T790M was not detectable in pre-treatment samples but was prevalent in post-TKI treated cancers. While this study is limited by its small cohort size, the findings indicate deep sequencing analysis of pre-treatment excision specimens of EGFR-mutant lung adenocarcinomas is warranted for detection of resistant mutations and predicting treatment response. Note: This abstract was not presented at the meeting. Citation Format: Xu-yuan Gao, Hang Xu, James CM Ho, Oscar SH Chan, Feng Xu, Junwen Wang, Victor HF Lee, Vicky PC Tin, Zhijie Xiao, Siqi Wang, Judy WP Yam, Maria P. Wong. Whole-exome sequencing reveals tyrosine kinase-resistant mutations in pretreatment EGFR-mutant lung adenocarcinomas [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 LB-274. doi:10.1158/1538-7445.AM2017-LB-274

Abstract Objective: To elucidate the genetic architecture of gene expression in pancreatic tissues. Design: We performed expression quantitative trait locus (eQTL) and allele specific expression (ASE) analyses using RNA-sequence data and 1000 Genomes (1000G) imputed GWAS genotypes from 95 fresh frozen histologically normal pancreatic tissue samples. Data from 115 pancreatic tumor-derived tissue samples from The Cancer Genome Atlas (TCGA) was included for comparison. Results: We identified 38,615 cis-eQTLs (corresponding to 484 Genes) in histologically normal tissues and 39,713 cis-eQTL (corresponding to 237 Genes) in tumor tissues (FDR<0.1), with the strongest effects seen near transcriptional start sites (TSS). Approximately 23% and 42% of genes with significant cis-eQTLs (eGenes) appeared to be specific for tumor and normal derived tissues, respectively. Significant enrichment of cis-eQTL variants was noted in noncoding regulatory regions marked by modified histones, DNAse hypersensitivity and bound transcription factors, in particular for pancreatic tissues (1.53-3.12 fold, P≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 in the ABO gene (rs687289) was associated with ABO expression in histologically normal (P=5.8x10-8) and tumor-derived (P=8.3x10-5) pancreatic tissues. The high linkage disequilibrium (LD) between this variant and the O blood group generating deletion variant in exon 6 of ABO suggested that nonsense-mediated decay (NMD) of the “O” mRNA could explain the eQTL. However, knock-down of crucial NMD regulators did not influence decay of the ABO “O” mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL. Conclusions: We have identified cis-eQTLs representing potential functional regulatory variants in the pancreas and generated a rich dataset for further studies on gene expression and regulation in pancreatic tissues. Citation Format: Laufey T. Amundadottir, Soren Lykke Andresen, Wenming Xiao, Jason Hoskins, Ashley Jermusyk, Lauren Rost, Irene Collins, Jinping Jia, Michael Mobaraki, Bin Zhu, Robert Kurtz, Hemang Parikh, Lei Song, Meredith Yeager, Torben Jensen, William Bamlet, Nilanjan Chatterjee, Brian Wolpin, Jill Smith, Sara Olson, Gloria Petersen, Jianxin Shi, Mingfeng Zhang. Analysis of cis-eQTLs in normal and tumor-derived pancreatic tissues reveals functional insights, including for the 9q34.1 ABO pancreatic cancer risk locus [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 1442. doi:10.1158/1538-7445.AM2017-1442

Tumor-initiating cells (TIC) are dynamic cancer cell subsets that display enhanced tumor functions and resilience to treatment but the mechanism of TIC induction or maintenance in lung cancer is not fully understood. In this study, we show the calcium pathway transcription factor NFATc2 is a novel regulator of lung TIC phenotypes, including tumorspheres, cell motility, tumorigenesis, as well as in vitro and in vivo responses to chemotherapy and targeted therapy. In human lung cancers, high NFATc2 expression predicted poor tumor differentiation, adverse recurrence-free and cancer-specific overall survivals. Mechanistic investigations identified NFATc2 response elements in the 3’ enhancer region of SOX2, and NFATc2/SOX2 coupling upregulates ALDH1A1 by binding to its 5’ enhancer. Through this axis, oxidative stress induced by cancer drug treatment is attenuated, leading to increased resistance in a mutation-independent manner. Targeting this axis provides a novel approach for the long-term treatment of lung cancer through TIC elimination.

AbstractPituitary stalk interruption syndrome (PSIS) is a rare type of hypopituitarism manifesting various degrees of pituitary hormone deficiency. Although mutations have been identified in some familial cases, the underpinning mechanisms of sporadic patients with PSIS who are in a vast majority remain elusive, necessitating a comprehensive study using systemic approaches. We postulate that other genetic mechanisms may be responsible for the sporadic PSIS. To test this hypothesis, we conducted a study in 24 patients with PSIS of Han Chinese with no family history using whole‐exome sequencing (WES) and bioinformatic analysis. We identified a group of heterozygous mutations in 92% (22 of 24) of the patients, and these genes are mostly associated with Notch, Shh, Wnt signalling pathways. Importantly, 83% (20 of 24) of the patients had more than one mutation in those pathways suggesting synergy of compound mutations underpin the pathogenesis of sporadic PSIS.

Objective To elucidate the genetic architecture of gene expression in pancreatic tissues. Design We performed expression quantitative trait locus (eQTL) analysis in histologically normal pancreatic tissue samples (n=95) using RNA sequencing and the corresponding 1000 genomes imputed germline genotypes. Data from pancreatic tumour-derived tissue samples (n=115) from The Cancer Genome Atlas were included for comparison. Results We identified 38 615 cis -eQTLs (in 484 genes) in histologically normal tissues and 39 713 cis -eQTL (in 237 genes) in tumour-derived tissues (false discovery rate <0.1), with the strongest effects seen near transcriptional start sites. Approximately 23% and 42% of genes with significant cis -eQTLs appeared to be specific for tumour-derived and normal-derived tissues, respectively. Significant enrichment of cis -eQTL variants was noted in non-coding regulatory regions, in particular for pancreatic tissues (1.53-fold to 3.12-fold, p≤0.0001), indicating tissue-specific functional relevance. A common pancreatic cancer risk locus on 9q34.2 (rs687289) was associated with ABO expression in histologically normal (p=5.8×10 −8 ) and tumour-derived (p=8.3×10 −5 ) tissues. The high linkage disequilibrium between this variant and the O blood group generating deletion variant in ABO (exon 6) suggested that nonsense-mediated decay (NMD) of the ‘O’ mRNA might explain this finding. However, knockdown of crucial NMD regulators did not influence decay of the ABO ‘O’ mRNA, indicating that a gene regulatory element influenced by pancreatic cancer risk alleles may underlie the eQTL. Conclusions We have identified cis -eQTLs representing potential functional regulatory variants in the pancreas and generated a rich data set for further studies on gene expression and its regulation in pancreatic tissues.

Abstract FAM46C is one of the most recurrently mutated genes in multiple myeloma; however its role in disease pathogenesis has not been determined. Here we demonstrate that wild-type (WT) FAM46C overexpression induces substantial cytotoxicity in multiple myeloma cells. In contrast, FAM46C mutations found in multiple myeloma patients abrogate this cytotoxicity, indicating a survival advantage conferred by the FAM46C mutant phenotype. WT FAM46C overexpression downregulated IRF4, CEBPB, and MYC and upregulated immunoglobulin (Ig) light chain and HSPA5/BIP. Furthermore, pathway analysis suggests that enforced FAM46C expression activated the unfolded protein response pathway and induced mitochondrial dysfunction. CRISPR-mediated depletion of endogenous FAM46C enhanced multiple myeloma cell growth, decreased Ig light chain and HSPA5/BIP expression, activated ERK and antiapoptotic signaling, and conferred relative resistance to dexamethasone and lenalidomide treatments. Genes altered in FAM46C-depleted cells were enriched for signaling pathways regulating estrogen, glucocorticoid, B-cell receptor signaling, and ATM signaling. Together these results implicate FAM46C in myeloma cell growth and survival and identify FAM46C mutation as a contributor to myeloma pathogenesis and disease progression via perturbation in plasma cell differentiation and endoplasmic reticulum homeostasis. Cancer Res; 77(16); 4317–27. ©2017 AACR.

Abstract Purpose: We sought to investigate whether B cell receptor immunoglobulin (BcR IG) stereotypy is associated with particular clinicobiological features among chronic lymphocytic leukemia (CLL) patients expressing mutated BcR IG (M-CLL) encoded by the IGHV4-34 gene, and also ascertain whether these associations could refine prognostication. Experimental Design: In a series of 19,907 CLL cases with available immunogenetic information, we identified 339 IGHV4-34–expressing cases assigned to one of the four largest stereotyped M-CLL subsets, namely subsets #4, #16, #29 and #201, and investigated in detail their clinicobiological characteristics and disease outcomes. Results: We identified shared and subset-specific patterns of somatic hypermutation (SHM) among patients assigned to these subsets. The greatest similarity was observed between subsets #4 and #16, both including IgG-switched cases (IgG-CLL). In contrast, the least similarity was detected between subsets #16 and #201, the latter concerning IgM/D-expressing CLL. Significant differences between subsets also involved disease stage at diagnosis and the presence of specific genomic aberrations. IgG subsets #4 and #16 emerged as particularly indolent with a significantly (P &lt; 0.05) longer time-to-first-treatment (TTFT; median TTFT: not yet reached) compared with the IgM/D subsets #29 and #201 (median TTFT: 11 and 12 years, respectively). Conclusions: Our findings support the notion that BcR IG stereotypy further refines prognostication in CLL, superseding the immunogenetic distinction based solely on SHM load. In addition, the observed distinct genetic aberration landscapes and clinical heterogeneity suggest that not all M-CLL cases are equal, prompting further research into the underlying biological background with the ultimate aim of tailored patient management. Clin Cancer Res; 23(17); 5292–301. ©2017 AACR.

AbstractMotivationNext generation sequencing technologies have enabled the study of the human microbiome through direct sequencing of microbial DNA, resulting in an enormous amount of microbiome sequencing data. One unique characteristic of microbiome data is the phylogenetic tree that relates all the bacterial species. Closely related bacterial species have a tendency to exhibit a similar relationship with the environment or disease. Thus, incorporating the phylogenetic tree information can potentially improve the detection power for microbiome-wide association studies, where hundreds or thousands of tests are conducted simultaneously to identify bacterial species associated with a phenotype of interest. Despite much progress in multiple testing procedures such as false discovery rate (FDR) control, methods that take into account the phylogenetic tree are largely limited.ResultsWe propose a new FDR control procedure that incorporates the prior structure information and apply it to microbiome data. The proposed procedure is based on a hierarchical model, where a structure-based prior distribution is designed to utilize the phylogenetic tree. By borrowing information from neighboring bacterial species, we are able to improve the statistical power of detecting associated bacterial species while controlling the FDR at desired levels. When the phylogenetic tree is mis-specified or non-informative, our procedure achieves a similar power as traditional procedures that do not take into account the tree structure. We demonstrate the performance of our method through extensive simulations and real microbiome datasets. We identified far more alcohol-drinking associated bacterial species than traditional methods.Availability and implementationR package StructFDR is available from CRAN.Supplementary informationSupplementary data are available at Bioinformatics online.