Tian, Shulan

Background: Multiple myeloma (MM) is an incurable plasma cell malignancy with several recurrent primary cytogenetic abnormalities including t(11,14)(q13;q32), resulting in increased expression of CCND1. Although 40% of patients with t(11;14) MM respond favorably to the BCL-2 inhibitor venetoclax (ven), 60% fail to achieve an objective response. Previous studies have suggested upregulation of PMAIP1, encoding the pro-apoptotic NOXA protein, in association with ven sensitivity (Leblay, Blood, 2024). Here we sought to evaluate the role for NOXA in ven response in t(11;14) MM. Methods: We analyzed whole exome, whole genome and RNA sequencing data of CD138+ tumor cells from 1,359 patients with MM in the CoMMpass and Mayo Clinic cohorts. The histone mark H3K27ac was assessed in t(11;14)+ KMS12PE and U266B1 human myeloma cell lines using ChIP-seq. Gene expression and protein levels of PMAIP1 (NOXA) were also detected by qRT-PCR and Western blot, respectively. Live cell detection was evaluated by Annexin V and PI using flow cytometry and by YoYo3 dye using IncuCyte following treatment with ven. Results: The t(11;14) MM subtype was associated with a 1.6-5.5-fold increased expression of PMAIP1 compared to non-t(11,14) MM in the CoMMpass and Mayo Clinic cohorts (p<0.0001). These findings are consistent with a significant correlation between CCND1 and PMAIP1 expression in both cohorts (p<0.001). To evaluate whether cyclin D1 is directly associated with increasing expression of PMAIP1 in t(11;14) MM, we generated a cyclin D1 knock out (KO) version of the U266B1 line using CRISPR/cas9. Two independent U266B1 KO cell lines showed reduced expression of PMAIP1 and NOXA protein compared to the WT U266B1 consistent with a role for cyclin D1 in regulating PMAIP1 expression. In addition, the KO lines had a 4.8-fold reduced H3K27ac signal at the PMAIP1 locus in comparison to the WT control suggesting a role for cyclin D1 in promoting H3K27ac. We previously observed increased expression of CCND1 in association with the G allele compared to the A allele of rs9344 in both the CoMMpass and Mayo Clinic cohorts (Leukemia, in press). We next asked whether the G risk allele of rs9344 within CCND1 impacted gene expression and H3k27ac levels at the PMAIP1 locus. To study the allele specific effect of CCND1, we used the CRISPR/cas9 edited KMS12PE WT (rs9344 AA) and Knock In (KI) (rs9344 GG) cell lines. In the KI cells, the G allele was accompanied by an ~2-fold increased expression of PMAIP1 and NOXA protein detected by western blot (<0.05). Increased PMAIP1 expression in KI lines was also accompanied by a 2.2-fold increased H3K27ac further suggesting a role for the G allele of rs9344 in promoting H3K27ac and gene expression at the PMAIP1 locus. Finally, we evaluated the consequence of PMAIP1 expression levels in relation to ven response. Increased expression of PMAIP1 and elevated levels of NOXA protein in the KMS12PE KI lines was associated with increased sensitivity to ven resulting in ~50% increase in cell death of KI vs. WT lines (p<0.02). Conclusion: We demonstrate a role for cyclin D1 in regulating PMAIP1 expression through H3K27ac resulting in increased NOXA protein. Levels of NOXA are associated with the G allele of rs9344 and increased venetoclax response.

Introduction: Cyclin D1, encoded by the CCND1 gene, is primarily known as a major regulator of cell cycle progression. Emerging studies have demonstrated a role for cyclin D1 in regulating gene transcription. In mantle cell lymphoma with t(11;14), cyclin D1 strongly binds to the promoters of highly transcribed genes and its overexpression led to genome-wide transcriptional downregulation. However, the transcriptional role for cyclin D1 in t(11;14) MM remains undefined. This study aims to elucidate the global chromatin and transcriptional changes induced by cyclin D1 in t(11;14) MM. Methods: Knockout (KO) of CCND1 in U266B1, a MM cell line carrying t(11;14), was performed by the Synthego (Redwood City, CA) using CRISPR/Cas9, which generated a frameshift insertion of base pair (bp) “A” after 69,641,478 bp within exon 1. Two independent KO clones were used for downstream analysis. ChIP-seq and ATAC-seq were performed for KO (n=2) and wild-type clones (WT, n=2). ChIP-seq used antibodies against H3K4me3 (promoter) and H3K4me1 and H3K27ac (enhancer). Libraries were sequenced to 2x51 bp on the Illumina NextSeq 2000 platform. Total RNA libraries (in triplicate) were sequenced to 2x150 bp. ChIP-seq and ATAC-seq reads were mapped to hg38 with BWA, and peaks were identified by MACS2 at FDR 0.01. RNA-seq reads were mapped to hg38 with STAR and gene raw counts were estimated using featureCounts. Differential peaks and differentially expressed genes between KO and WT were identified using edgeR at FDR 0.05 and fold-change 2. Only peaks present in at least 2 samples and genes having reads-per-million value 0.1 in at least 1 sample were included, excluding those from chrY. Pathway analysis was performed with the Enrichr package. Results: Differential expression analysis identified 2-fold more genes that were down-regulated than up-regulated in KO compared to WT (467 vs. 233). Pathway analysis revealed an enrichment of TNF-alpha signaling via NF-kappaB, inflammatory response, and apoptosis in association with down-regulated genes, while up-regulated genes are enriched in myogenesis, KRAS signaling, and epithelial mesenchymal transition pathway. In parallel, differential chromatin analysis identified 2.5% (1,161) of the peaks from H3K27ac and 0.5-1% (359-610) of the peaks for the other 3 marks that showed significant changes of signal levels between KO and WT. For the promoter mark H3K4me3, there are largely a comparable number of peaks showing increased (191 peaks) and decreased (168 peaks) signals, respectively, in KO compared to WT. However, for ATAC-seq and H3K4me1, about 70% of the differential peaks showed an increased signal in KO, in contrast with H3K27ac for which 60% of differential peaks showed decreased signal in KO. The preferential loss of H3K27ac occupancy is consistent with the observation of a higher proportion of genes being down-regulated in KO. These results indicated that knockout of CCND1 led to a global increase of chromatin accessibility, but a reduction of H3K27ac and gene expression. Further annotation revealed that, for H3K27ac and H3K4me1/ATAC-seq, 92% and 97% of the differential peaks were located >2kb away from the transcription start sites (TSSs) of expressed genes. Even for H3K4me3, 88% of the differential peaks are >2kb away from TSSs. The results suggest that the changes of chromatin upon CCND1 knockout predominantly occurred at distal regulatory regions. For up- and down-regulated genes, we estimated ChIP-seq and ATAC-seq signals in the TSS 2kb regions, revealing that up-regulated genes were overall associated with increased signals. Further analysis found that, for ATAC-seq and H3K4me3, 44% and 59% of the differential peaks located in the promoter regions were associated with differentially expressed genes, which dropped to ~20% for the two enhancer marks. These data indicate that changes in chromatin state impact the expression of the associated genes. Conclusion: Knockout of CCND1 impacts local chromatin state, particularly at enhancer regions, and the transcriptional program. In promoter-proximal regulatory regions, the changes of chromatin accessibility and H3K4me3 are strongly correlated with the expression changes of nearby genes. Further studies will identify the genes targeted by the differential enhancers and their possible roles in MM pathogenesis mediated by cyclin D1 following the t(11;14) event.

Abstract Glioblastoma (GBM) remains one of the most therapy-resistant malignancies with frequent local failures despite aggressive surgery, chemotherapy, and ionizing radiation (IR). Small molecule inhibitors of DNA-dependent protein kinase (DNA-PKi) are potent radiosensitizers currently in clinical trials. Determining which patients may benefit from radiosensitization with DNA-PKi is critical to avoid unnecessary increased risk of normal tissue toxicity. In this study, we used GBM patient-derived xenografts (PDX) in orthotopic murine models to study the relationship between molecular features, pharmacokinetics, and the radiosensitizing potential of the DNA-PKi peposertib. We show that peposertib radiosensitizes established and PDX GBM lines in vitro at 300 nmol/L and above, with a significant increase in radiosensitization by maintaining post-IR exposure for >12 hours. Radiosensitization by peposertib is mediated by catalytic inhibition of DNA-PK, and knockdown of DNA-PK by short hairpin RNA (shRNA) largely abolished the radiosensitizing effect. Peposertib decreased auto-phosphorylation of DNA-PKcs after IR in a dose-dependent manner with a delay in resolution of γH2AX foci at 24 hours. The addition of peposertib to IR significantly increased survival in GBM120 orthotopic xenografts, but not in GBM10. There was no difference in plasma or average tumor concentrations of peposertib in the two cohorts. Although the mechanism underpinning this discordant effect in vitro versus in vivo is not clear, there was an association for greater sensitization in TP53 mutant lines. Transfection of a dominant-negative TP53 mutant in baseline TP53 wild-type GBM lines significantly delayed growth and decreased nonhomologous end joining efficiency (but not homologous recombination), after peposertib exposure. See related commentary by Buchsbaum, p. 840

Abstract Multiple myeloma (MM) is the most common blood cancer in African American (AA) individuals. AAs have an increased risk of developing Monoclonal Gammopathy of Undetermined Significance (MGUS) and MM compared to non-Hispanic White (NHW) individuals. While several studies have shown a pivotal role for the immune tumor microenvironment (iTME) in MM pathogenesis, the composition of the iTME in AAs with MM is under investigated. Viably frozen whole bone marrow (BM) aspirates from 128 individuals (53 AA and 75 NHW) with either MM (n=65), MGUS (n=28), smoldering MM (SMM, n=13), or healthy controls (n=22) were evaluated by single cell RNA-sequencing (scRNA-seq) or mass cytometry (CyTOF). We characterized 34 (14 AA and 20 NHW) by scRNA-seq and 94 (39 AA and 55 NHW) by CyTOF. CyTOF samples were stained with Maxpar Direct Immune Profiling Assay Kit (Fluidigm), acquired on a Helios mass cytometer and analyzed in Cytobank. Data were analyzed by logistic equivalent regression model with p<0.05 for significance. ScRNA-seq samples were processed as in Yao, 2022. Data were normalized using scvi-tools and clustered by Uniform manifold approximation and projection (UMAP). Differentially expressed markers between clusters were identified using Model-based Analysis of Single-cell Transcriptomics (MAST) in Seurat R package at FDR < 0.05 and log2 fold change > 0.55. ScRNA-seq of the CD45+ non-tumor lymphocytes revealed 17 immune populations (HSC, erythroblast, MDSC, GMP, CD1c+ DC, CD4+ and CD8+ T cell, CD14+ M1 macrophage, CD16+monocyte, monocyte, pDC, naïve B, pre-B, pro-B, memory B, NK/MAIT). CD3E/D+cells were clustered into naïve/central memory (CCR7+, LEF1+, SELL+), memory (GZMK+), effector (GZMB+, GZMH+, FGFBP2+), g/d (TRDC+), MAIT (SLC4A10+, KLRB1+) and IFN CD8+ T cells (IFIT1+, IFIT2+,IFIT3+). Disease progression was associated with increases in median CD4+ and CD8+ GZMK+ T cells, memory B cells and a reduction in all monocyte clusters. Within the diseased cohort in comparison to NHW patients, AA patients had increased CD3+ T cells (CD4 and CD8+) (AA 57.5% vs. NHW 44.3%) and a reduction in overall monocytes (AA 10.5% vs. 26.5%). These results were validated by CyTOF, which confirmed an increase in CD3+ T cells (AA 53.7% vs. NHW 35.7%, p<0.001) including increases in CD8+ (AA 7.0% vs. NHW 3.2%, p<0.05) and CD4+ terminal effector T cells (AA 5.7% vs. NHW 2.7%, p<0.05) and a reduction in monocytes (AA 18.3% vs. NHW 26.6%) among AAs. After adjusting for age and sex, AAs retained an increase in T cells (OR 1.59 (95%CI 1.2-2.01, p=0.00023)) and a reduction in monocytes (OR 0.88 (95%CI 0.78-1, p=0.045)). Despite the increase in CD8+ T cells in AA patients, we identified reduced GZMK+ T cells and reductions in exhaustion markers (TIGIT, LAG3, TOX) in CD8+ T cells in comparison to NHWs. Nevertheless, the distribution of the immune populations was similar between healthy AA and NHW individuals. Our multiomics approach to characterize the iTME identified significant differences among AA compared to NHW patients with MM. Citation Format: Suganti Shivaram, Hongwei Tang, Huihuang Yan, Shulan Tian, Guyu Qin, Emilie Anderson, Neeraj Sharma, Abiola Bolarinwa, Cinthya Zepeda-Mendoza, Stacey Lehman, Felicia Gomez, Amit Mitra, Surendra Dasari, Taxiarchis Kourelis, Shaji Kumar, Linda Baughn. Characterization of the Immune Tumor Microenvironment in a Diverse Cohort of Patients with Multiple Myeloma [abstract]. In: Proceedings of the Blood Cancer Discovery Symposium; 2024 Mar 4-6; Boston, MA. Philadelphia (PA): AACR; Blood Cancer Discov 2024;5(2_Suppl):Abstract nr P06.

Introduction: Despite the improved clinical outcomes with novel therapies in patients with chronic lymphocytic leukemia (CLL), Richter transformation (RT) of CLL remains a major clinical challenge. Our group showed a selective response (~40%) with the PD-1 blocking antibody pembrolizumab in patients with RT, particularly after prior exposure to ibrutinib (Ding et al, Blood, 2017). It has been established that PD-1 inhibits cytotoxic T-cell functions via the interaction with its ligands, PD-L1 and PD-L2. However, the role of PD-1 in tumor B-cells is understudied. Here we show the expression and functional implications of the PD-1 signaling axis in B-cell pathobiology. Methods: CLL-involved lymph node (LN) (n=26) and RT (with diffuse large B-cell lymphoma [DLBCL] histology) -involved LN (n=20) samples were tested for PD-1 expression by immunohistochemistry (IHC) (mouse clone NAT105, Abcam). Then, PD-1 expression was examined in 10 lymphoma cell lines and 1 CLL-like cell line (MEC-1) by both flow cytometry (FC) and Western blot (WB) analysis. PD-1 overexpression using pLEX-lentiviral packaging system (Thermo Scientific) was evaluated for its impact on in vitro cell cycle regulation and in vivo tumor growth. CRISPR/CAS9 mediated PD-1 knockdown evaluated on cell growth in MEC-1 and Mino cells. Role of PD-1 was further examined on the regulation of cell signaling pathways using human phospho kinase array kit (R&D) and WB analysis. Additionally, gene expression signatures were evaluated by RNA-seq using pLEX-lentiviral-PD-1 transfected DLBCL cell lines OCI-LY7 and OCI-LY19. Results: PD-1 expression was significantly increased in RT-LN compared to CLL-LN (mean ± SEM in RT vs. CLL, 30.6% ± 4.7% vs. 11.5% ± 2.8%, p<0.001) as determined by IHC. PD-1 expression was highest in patients with RT where the last prior CLL therapy was ibrutinib. To test the role of PD-1 in tumor B-cells, its expression was assessed in lymphoma and CLL cell lines. The expression of PD-1 was found to be variable by both FC and WB analysis, but at very low-levels in DLBCL cell lines OCI-LY7 and OCI-LY19. Constitutive lentiviral (pLEX-PD-1)-mediated overexpression of PD-1 in OCI-LY7 and OCI-LY19 cells led to increased cell growth (1.4 and 1.9 fold, respectively, compared to original lines after 4 days of culturing of equal number of cells), which was further confirmed by cell cycle analysis which showed an increase in S phase by 20.4% (p<0.05) and 24.5% (p<0.05) in OCI-LY7 and OCI-LY19 cells, respectively. When the luciferase+ PD-1 overexpressing OCI-LY7 cells (1X10 6) were injected intravenously into NSG mice, increased tumor growth (p<0.01) was observed by bioluminescence imaging after 22 days of follow-up. PD-1 knockdown in MEC-1 and Mino cells (high basal PD-1 expression) led to significant inhibition of cell growth (85% and 79.5%, respectively, p<0.05). Using phospho-kinase array, an overall decrease of phosphorylation was detected on multiple sites of p53 (S392, S15, S46) and Chk-2 (T68) in PD-1 overexpressing OCI-LY7 and OCI-LY19 cell lines. This finding was further confirmed by WB analyses. In addition, decreased total p53 and increased phosphorylation on both SHP-1 and SHP-2 were found in PD-1 overexpressing OCI-LY7 and OCI-LY19 cell lines. SHP-1 (TPI-1) and SHP-1+SHP-2 (TPI-1+TNO155) inhibitor treatment promoted OCI-LY19 cell death and led to recovery of phosphorylation on ATM, Chk-2 and p53 in these cells. RNA sequencing data in PD-1 overexpressing OCI-LY7 and OCI-LY19 cells revealed increased expression of genes that promote cell cycle progression and DNA replication. Conclusion: Our data showed robust PD-1 expression in patients with Richter transformation from CLL to DLBCL. PD-1 overexpression in DLBCL lymphoma cell lines enhanced cell proliferation in vitro and in vivo. Further investigation identified that PD-1 modified the phosphorylation/function of SHP phosphatases and thereby regulated p53 pathways ( Figure 1). These findings indicate that RT with DLBCL cells have the ability to “hijack” the PD-1 pathway which results in downregulation of p53 mediated DNA repair. This novel finding provides for new strategies to further evaluate the interaction of PD-1 signals with DNA repair pathways and suggest potential novel targets for RT.

Abstract Background High-grade gliomas (HGG) are aggressive brain tumors associated with short median patient survival and limited response to therapies, driving the need to develop tools to improve patient outcomes. Patient-derived xenograft (PDX) models, such as mouse PDX, have emerged as potential Avatar platforms for personalized oncology approaches, but the difficulty for some human grafts to grow successfully and the long time required for mice to develop tumors preclude their use for HGG. Methods We used a rapid and efficient ex-ovo chicken embryo chorioallantoic membrane (CAM) culture system to evaluate the efficacy of oncologic drug options for HGG patients. Results Implantation of fresh glioma tissue fragments from 59 of 60 patients, that include difficult-to-grow IDH-mutated samples, successfully established CAM tumor xenografts within 7 days, with a tumor take rate of 98.3%. These xenografts faithfully recapitulate the histological and molecular characteristics of the primary tumor, and the ability of individual fragments to form tumors was predictive of poor patient prognosis. Treatment of drug-sensitive or drug-resistant xenografts indicates that the CAM-glioma assay enables testing tumor sensitivity to temozolomide and carboplatin at doses consistent with those administered to patients. In a proof-of-concept study involving 14 HGG patients, we observed a correlation of 100% between the CAM xenograft response to temozolomide or carboplatin and the clinical response of patients. Conclusion The CAM-glioma model is a fast and reliable assay that has the potential to serve as a complementary model to drug discovery and a real-time Avatar platform to predict the best treatment for HGG patients.

Abstract The Mayo Clinic collection of patient-derived xenograft (PDXs) has been developed from the direct implantation of glioblastoma (GBM) samples into immunodeficient mice. Over 105 models have been extensively characterized with regards to histology, invasion, flank and orthotopic growth rates, and molecular profiling. Evaluation of orthotopic tumors demonstrated infiltration of tumor cells into the brain parenchymal in 90% of models with contralateral hemisphere involvement in 64%. Extensive molecular profiling confirmed the PDX collection captures the genetic heterogeneity of GBM, with the majority of typical molecular alterations seen at a frequency similar to the TCGA including gain of chromosome 7 and loss of chromosome 10. TERT promoter mutations (86%) and homozygous deletion of CDKN2A (70%) were the most frequent. Additional alterations in the p53 pathway included TP53 mutations (36%) and MDM2/4 amplification (10%/2%). IDH mutations were found in three models while MGMT promoter methylation was observed in 45% of PDXs. The available molecular sequencing results are viewable on cBioPortal. These genetically diverse models are extremely useful for biomarker discovery and efficacy evaluations. Using an innovative 1x1x1 PDX preclinical trial design, we are establishing a workflow to evaluate control and drug treatment in an individual animal across PDX lines in the flank. As our first use of this strategy, a small molecule MDM2 inhibitor was evaluated in 17 PDX lines with variable MDM2 and TP53 statuses. Tumors harboring MDM2 amplification without TP53 mutation (n=5) showed profound response, with 4 PDXs showing a complete response that lasted >150 days. Response in TP53 wild-type tumors lacking MDM2 amplification was variable, with a 1.9 to 4.8 fold longer time to endpoint (mean 3; median 2.7). A more traditionally powered intracranial efficacy study performed with one MDM2 amplified (GBM108) and two non-amplified (GBM14, GBM10) PDX lines confirmed these findings of reduced and more variable response to therapy in MDM2 non-amplified TP53 wild-type tumors. Future studies will aim to examine novel biomarkers indicative of therapeutic response to better understand molecular features associated with MDM2 inhibitor efficacy. With over 20 years of experience, our group has extensive expertise in development of clinically relevant study designs focused on examining efficacy, tolerability, and pharmacokinetic/pharmacodynamic assessments. Our diverse collection of clinically relevant GBM PDX models provides a highly sought-after resource which has been utilized by the pharmaceutical community with 58 working relationships to date (26 contracts, 32 collaborations). Through these collaborations, we aim to expand the use of our highly characterized GBM PDX panel for biomarker discovery using similar strategies. Citation Format: Danielle M. Burgenske, Ann C. Mladek, Shulan Tian, Mark A. Schroeder, Zeng Hu, Brett L. Carlson, Paul A. Decker, Jeanette E. Eckel-Passow, Rachael Vaubel, Jann N. Sarkaria. Evaluation of novel therapeutics using the Mayo Clinic GBM Patient-Derived Xenograft (PDX) collection [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 4001.

Abstract Background RBBP4 activates transcription by histone acetylation, but the partner histone acetyltransferases are unknown. Thus, we investigated the hypothesis that RBBP4 interacts with p300 in a complex in glioblastoma (GBM). Methods shRNA silencing of RBBP4 or p300 and RNAseq was used to identify genes co-regulated by RBBP4 and p300 in GBM43 patient-derived xenograft (PDX). RBBP4/p300 complex was demonstrated using proximity ligation assay (PLA) and ChIPseq delineated histone H3 acetylation and RBBP4/p300 complex binding in promoters/enhancers. Temozolomide (TMZ)-induced DNA double strand breaks (DSBs) were evaluated by γ-H2AX and proliferation by CyQuant and live cell monitoring assays. In vivo efficacy was based on survival of mice with orthotopic tumors. Results shRBBP4 and shp300 downregulated 4768 genes among which 1485 (31%) were commonly downregulated by both shRNAs, while upregulated genes were 2484, including 863 (35%) common genes. The pro-survival genes were the top-ranked among the downregulated genes, including C-MYC. RBBP4/p300 complex was demonstrated in the nucleus, and shRBBP4 or shp300 significantly sensitized GBM cells to TMZ compared to the control shNT in vitro (P < .05). Moreover, TMZ significantly prolonged the survival of mice bearing GBM22-shRBBP4 orthotopic tumors compared with control shNT tumors (median shNT survival 52 days vs. median shRBBP4 319 days; P = .001). CREB-binding protein (CBP)/p300 inhibitor CPI-1612 suppressed H3K27Ac and RBBP4/p300 complex target proteins, including C-MYC, and synergistically sensitized TMZ in vitro. Pharmacodynamic evaluation confirmed brain penetration by CPI-1612 supporting further investigation to evaluate efficacy to sensitize TMZ. Conclusions RBBP4/p300 complex is present in GBM cells and is a potential therapeutic target.