Citations

Abstract Background:. Poor patient outcomes in triple negative breast cancer (TNBC) largely stem from a lack of understanding of therapeutic vulnerabilities and an insufficient armamentarium of effective drugs. The standard of care for most early stage TNBC patients continues to be (neo)adjuvant chemotherapy and radiation. Identifying additional therapeutic options for this subset of patients represents a major unmet need. We and others demonstrated that estrogen receptor beta (ERβ) is expressed in about 20% of TN tumors. Prior research has shown that ligand-mediated activation of ERβ decreases proliferation, invasion, and migration in TNBC cell lines. In vivo, ERβ suppresses the growth of cell line xenograft models and prevents the development of metastatic lesions in a ligand dependent manner. Mechanistically, we found that ERβ repurposes EZH2 to suppress oncogenic NFκB signaling and that pharmacologic inhibition of EZH2 diminishes ERβ function and its anti-cancer effects. This research has led to an ongoing clinical trial (NCT03941730) assessing the efficacy of estradiol in ERβ+ TNBC patients with chemorefractory disease. As with all therapies, patients will exhibit de novo and acquired resistance and therefore we sought to understand the mechanisms of resistance to ERβ targeted therapies. Methods:. Using multiple models of ERβ positive TNBC, we developed ERβ resistant cell lines through chronic exposure to estradiol and the ERβ specific agonist LY500307 over a period of 8 months. We employed RNA sequencing to characterize the transcriptomic changes which occurred following ERβ resistance. We profiled XIST expression across multiple publicly available datasets, including TCGA, Metabric, BEAUTY, and GTEx. XIST expression was also modulated using CRISPR/Cas9 to assess subsequent effects on TNBC cell biology, ERβ function, and response to ERβ targeted therapies. Results:. Our resistant cell line models of ERβ positive TNBC maintained ERβ expression, but were no longer growth inhibited by ERβ agonists. RNAseq revealed substantial differences comparing the transcriptome of resistant versus sensitive cell lines, of which the most increased transcript in the resistant setting was the lncRNA XIST. XIST is best known for its role in X-chromosome inactivation through recruitment and association with the PRC2 complex, and thus EZH2. However, little is known about its functions in breast cancer. We therefore assessed the expression levels of XIST in breast tumors and cell lines. XIST expression was highly variable in breast cancer, was found in a proportion of all breast cancer sub-types and did not correlate with ERβ status. However, XIST expression was up-regulated in ERβ expressing cell lines following long term estrogen treatment. No effects on XIST expression were identified in multiple ERα positive models following estrogen, SERM/SERD treatment, or estrogen deprivation. Furthermore, upregulation of XIST was not associated with resistance to Paclitaxel or Doxorubicin, suggesting XIST is not a part of a broad resistance phenotype. Strikingly, CRISPR mediated knockout of XIST in ERβ resistant cells completely re-sensitized cells to ERβ targeted therapies suggesting that XIST expression is a critical component of ERβ resistance. Conclusions:. ERβ represents a relevant therapeutic target that is being tested in clinical trials. Using multiple in vitro models, we provide evidence that XIST expression is sufficient to induce resistance to. ERβ targeted therapies in TNBC and may therefore represent a relevant biomarker for patient stratification. Further strategies to suppress XIST expression may elicit anti-cancer effects on their own and may resensitize a sub-set of ERβ resistant tumors to ERβ agonists. Citation Format: Michael J Emch, Kirsten GM Aspros, Elizabeth S Bruinsma, Krishna R Kalari, Calley J Jones, Brandon W Simone, Matthew P Goetz, John R Hawse. The lncRNA XIST mediates sensitivity to ERβ targeted therapies in triple negative breast cancer [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P4-02-09.

Abstract Background: FGFR dysregulation is observed in multiple cancers and targeting FGFR is an emerging therapeutic strategy with FDA approved treatments in bladder and cholangiocarcinoma. Here we examined the prevalence of FGFR mutations, fusions, and high-level amplifications in breast cancer, stratified by receptor subtype and local/metastatic status, in both Foundation Medicine (FM) and institutional Mayo Clinic (MC) cohorts. Methods: For the FM cohort, comprehensive genomic profiling (CGP) examining at least 324 genes for all classes of alterations, including FGFR1-4 was carried out for 32,048 breast cancers during the course of routine clinical care in a Clinical Laboratory Improvement Amendments (CLIA)-certified lab (Foundation Medicine Inc., Cambridge, MA, USA). Tumor mutational burden (TMB) was determined on up to 1.1 Mb, microsatellite instability high (MSI-High) was determined on up to 114 loci and predicted ancestry from >10,000 SNPs. Estrogen receptor (ER) and HER2 status were available for a subset of FM samples. Additionally, 131 patients with metastatic breast cancer from a subset of patients at three Mayo Clinic sites (MC cohort) with clinical characteristics and cancer-panel DNA sequencing data from a CLIA-certified lab (Tempus, Chicago, IL) were included. Results: In the FM cohort, the prevalence of FGFR1-4 high-level amplification (CN≥10) was 10.1%, while mutations (1.5%) and fusions (0.72%) were rare. Most amplifications occurred in FGFR1 (9.2%); most fusions and mutations occurred in FGFR2 (0.46%, 0.77%). FGFR alteration prevalence was highest in ER+/HER2- subtype (14.4%) and lowest in HER2+ disease (7.7%). FGFR alterations were more common in IDC (11.7%) than ILC (7.7%), p<3E-08. FGFR alterations were more prevalent in the metastatic setting relative to breast-biopsied disease (13.6% v 10.1%; OR = 1.4; p=2E-17), especially in the HER2+ (OR =1.9, p=0.004) and ER-/HER2- (OR = 1.9, p = 0.05) disease; no enrichment was seen in the ER+/HER2- metastases (OR =1.0, p = 1). FGFR amplifications were observed at a higher prevalence in patients with predicted East Asian ancestry, relative to patients with European ancestry (12.1% v 10.0%; p = 0.03). Overall, the most common activating mutations in FGFR were FGFR2 N549K (n=85), FGFR1 N546K (n=78), FGFR4 V510M (n=28), FGFR2 K659E (n=28), FGFR4 V510L (n=20), and FGFR2 Y375C (n=15). The most common recurrent fusions were FGFR3:TACC3 (n=36), FGFR2:TACC2 (n=17), FGFR1:TACC1 (n=9), FGFR1:BAG4 (n=6), and FGFR2:ATE1 (n=5). In patients with FGFR amplifications, the most frequently co-occurring alterations were ZNF703 (78.4%), TP53 (51.5%), CCND1 (36.1%), FGF3/4/19 (32.9 - 34.4%), PIK3CA (30.7%), MYC (29.6%), ESR1 (17.2%), EMSY (16.3%), and PTEN (10.6%). Significant co-occurrence was observed for a number of genes including FGF3/4/19, CDK4, and CDK8 (all OR>2, p<1E-07); mutual exclusivity was observed with PIK3R1, BRCA1, and BRCA2 (all OR <0.5, p<4E-13), among other genes. In the 131 metastatic tumors from MC, the prevalence of FGFR1-4 high-level amplifications was 19.8% [FGFR1 (12.4%), FGFR2 (7.4%), and FGFR3 (0.8%)]. The prevalence of high-level FGFR amplifications did not differ by clinical subtypes: HR-/HER2- (7/31), HR+/HER2- (15/79), and HER2+ (2/11), p=0.68. Conclusions: High-level FGFR amplifications are observed in >11% of breast cancers, especially the ER+/HER2- subtype, while mutations/fusions are rare. These data support the ongoing studies evaluating targeted therapies for FGFR amplified ER + breast cancer. Correlations with clinical information (MC cohort) and associations with actionable alterations are ongoing and may inform potential combination strategies. Citation Format: Karthik V Giridhar, Ethan S Sokol, Peter T Vedell, Jason P Sinnwell, Aakash Desai, Tufia C Haddad, Ciara C O’Sullivan, Roberto A Leon-Ferre, Siddhartha Yadav, Kostandinos Sideras, Brenda Ernst, Minetta C Liu, Abe Eyman Casey, Xiaojia Tang, Zoe Fleischmann, Karthikeyan Murugesan, Krishna R Kalari, Matthew P Goetz. The frequency and somatic mutation landscape of Fibroblast growth factor receptor (FGFR) alterations in breast cancer [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P3-08-02.

Abstract Background: The importance of the antitumor immune response in TNBC is well established. TNBC with higher TILs are less likely to recur and more responsive to systemic therapy. Likewise, PD-L1+ TNBC are more likely to benefit from chemoimmunotherapy. However, TNBC is highly heterogeneous. Of the TNBC molecular subtypes, LAR TNBC is less sensitive to systemic therapy, has lower TILs and lower rates of PD-L1 positivity. The role of other immune related proteins in LAR TNBC is not well established. Here, we evaluated differentially expressed (DE) immune related proteins in the stromal and intratumoral compartments of LAR vs non-LAR TNBC tumors. Methods: We used the Nanostring GeoMX DSP platform to quantitate 58 proteins within spatially distinct intraepithelial, cytokeratin (CK)-positive tumor segments and adjacent CK-negative/nuclei-positive stromal segments in 248 TNBC tumors included in a tissue microarray generated from a cohort of pts with centrally confirmed TNBC who underwent breast surgery without prior neoadjuvant therapy. A subset (n=111) underwent bulk tumor RNA sequencing and were classified as LAR or non-LAR TNBC. DE proteins were identified using a negative binomial generalized linear model (SNR>2, p<0.05). A targeted set of DE proteins was dichotomized at the 80th percentile. Results: Of 111 TNBC tumors, 17 (15%) were LAR and 94 (85%) non-LAR. Compared to non-LAR TNBC, pts with LAR TNBC were older (age ≥50: 82% vs 52%, p<0.01), with tumors that were more often of apocrine histology (35% vs 0%, p <0.01), grade 1-2 (24% vs 1%, p<0.01), and had lower Ki67 (Ki67 ≤15: 24% vs 11%, p=0.06). Most tumors were T1-2 (94% vs 93%, p=0.82) and N0 (53% vs 62%, p=0.09), respectively. As expected, expression of most immune-related proteins was higher in the stromal vs the intratumoral compartment for both LAR and non-LAR TNBC. When focusing on the stromal compartment, expression of multiple immune related proteins was significantly lower in LAR compared to non-LAR TNBC, including the pan-leukocyte marker CD45 (log-2 fold change [log2FC]: 0.552, p=0.05), the macrophage marker CD14 (log2FC: 0.834, p=0.06), CD44 (lof2FC: 0.637, p=0.07), and the immune checkpoint proteins IDO1 (log2FC: 0.914, p=0.04), VISTA (log2FC: 0.471, p=0.07), ICOS (log2FC: 0.444, p=0.08), and STING (log2FC: 0.544, p=0.09). Proteins with expression levels too low for comparisons included PD-L1, LAG3, FOXP3 and BCL-2. When focusing on the intratumoral compartment, expression of most immune-related proteins was very low in both LAR and non-LAR TNBC. Like in the stromal compartment, CD45 expression was lower in LAR TNBC (log2FC: 0.78, p=0.02). Expression of the immune checkpoint B7-H3 was lower in LAR TNBC (log2FC: 0.737, p=0.02), while expression of the T cell marker CD127 was higher (log2FC: -0.528, p=0.34). With regards to relevant non-immune markers, expression of Ki67 was lower in LAR TNBC (log2FC: 0.5498, p=0.05), consistent with the clinical assay. Conclusion: In this ultra high-plex spatial analysis, we provide first insights into the differential expression at the protein level of several targetable immune checkpoint molecules in LAR vs non-LAR TNBC. The lower expression of several immune related proteins in LAR TNBC is consistent with the hypothesis that LAR TNBC exhibits a “cold” immune microenvironment compared to other TNBC subtypes, potentially rendering itself less susceptible to immunotherapy-based strategies. These data support the need to consider TNBC molecular subtypes in future evaluations of immune-based therapeutic approaches. Funding: This work was supported by NIH grant P50CA116201 to RLF, JMC, KRK, FJC, DZ, JNI, and MPG; BCRF grant 19-161 to EAT and NCATS grant CTSA KL2 TR002379 to RLF. The contents are solely the responsibility of the authors and do not necessarily represent the official views of the NIH Citation Format: Roberto A Leon-Ferre, Jodi M. Carter, David M. Zahrieh, David W. Hillman, Saranya Chumsri, Yaohua Ma, Jennifer M. Kachergus, Xue Wang, Judy C. Boughey, Minetta C. Liu, James N. Ingle, Krishna R. Kalari, Jose C. Villasboas Bisneto, Fergus J. Couch, E. Aubrey Thompson, Matthew P. Goetz. Digital spatial profiling of immune-related proteins in luminal androgen receptor (LAR) vs non-LAR triple-negative breast cancer (TNBC) [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr P1-04-01.

Abstract Background: Although major advances have been made in the treatment of HER2+ metastatic breast cancer (MBC), the goal of care remains largely palliative, therefore better treatments are needed. Given encouraging preclinical and clinical data, the combination of cyclin dependent 4/6 kinase inhibitors and HER2-directed therapy is further being evaluated in this trial. Trial Design: In this phase I/II multicenter trial, we will determine the maximum tolerated dose (MTD) of abemaciclib (Abem) combined with T-DM1. Three Abem dose levels will be examined, 50 mg, 100 mg and 150 mg. The phase 2 portion of this trial will examine whether PFS is increased with addition of Abem to T-DM1 in two pt cohorts - those with ER+ HER2+ MBC and those with ER-HER2+ MBC. For phase 2, a pre-registration biopsy is required to confirm ER and HER2-status and to determine levels of tumor infiltrating lymphocytes; vimentin (an epithelial-mesenchymal transition marker); and CD8 and FOXP3 expression. Blood samples will be collected pretreatment , at 6 weeks, and at progression for all pts. Eligibility Criteria: Phase I&II: All pts must have HER2+ MBC per ASCO-CAP guidelines and prior treatment with a taxane, trastuzumab and pertuzumab. For the phase Iportion, pts can have measurable or non-measurable disease with no restriction on the number of prior lines of therapy. In the phase II portion, pts must have measurable disease with ≤1-2 prior lines of chemotherapy alone, ≤1 HER2-directed therapy alone, and/or chemotherapy with HER2-directed therapies. There is no limit on prior endocrine therapy. Specific Aims: The primary objective of the phase II trial is to assess whether addition of Abem to T-DM1 increases PFS in one or both patient cohorts. Secondary objectives include an assessment of toxicity, objective response rates and overall survival. Correlative studies will assess the association between baseline TIL levels, vimentin expression, and CD8/FOXP3 expression with PFS. Changes in peripheral blood mononuclear cells, CTCs, ctDNA and serum thymidine kinase 1 during the course of treatment will be examined to determine if there is a link with PFS outcomes overall and separately for each cohort. Pharmacogenomic studies will determine if pts with the FCGR3A-158 polymorphism derive less benefit from T-DM1 and have inferior PFS outcomes compared with pts who do not have this polymorphism. Statistical Methods:Phase I: Standard 3+3 design, with dose limiting toxicities as per protocol. Phase II: For each pt cohort, a stratified randomization scheme will be used to assign pts to treatment with liver mets as a stratification factor. For each pt cohort, a stratified log rank test will be used to assess whether PFS is increased with the addition of Abem to T-DM1. A non-binding futility analysis will be applied in each cohort after 58 events in the ER+ HER2+ MBC study cohort and 48 events in the ER- HER2+ MBC study cohort. Present Accrual: 0; target accrual: minimal 120 pts., maximal 140 pts CohortOne-sidedalphaPowerAccrual Period (accrual rate)Follow-up after close of enrollmentPFS with T-DM1PFS with abema and T-DM1Number of eligible patientsER+/HER2+0.100.912 months (5-6 pts per month)12 months12 weeks24 weeks64 (32 per arm)ER-/HER2+0.100.8512 months (3-4 pts per month)12 months6 weeks12 weeks50 (25 pts per arm) Citation Format: Ciara C O Sullivan, Jun He, Vera J Suman, Krishna R Kalari, Roberto A Leon-Ferre, Jose C Villasboas-Bisneto, Pratima Chalasani, Demet Gokalp Yasar, Daniel M Anderson, Philip J Stella, Anthony J Jaslowski, Susan H Tannenbaum, Angela Saverimuthu, Donald Northfelt, Alvaro Moreno-Aspitia, Jodi M Carter, Minetta C Liu, Liewei Wang, Zhenkun Lou, Matthew P Goetz. A phase I/II trial of abemaciclib and T-DM1 in women and men with HER2-positive advanced or metastatic breast cancer that has progressed on treatment with a taxane, trastuzumab and pertuzumab (THP) (ACCRU-BR-1801) [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr OT2-19-05.

Background: Age is the most common risk factor for Alzheimer’s disease (AD), a neurodegenerative disorder characterized by the hallmarks of toxic amyloid-β (Aβ) plaques and hyperphosphorylated tau tangles. Moreover, sub-physiological brain insulin levels have emerged as a pathological manifestation of AD. Objective: Identify age-related changes in the plasma disposition and blood-brain barrier (BBB) trafficking of Aβ peptides and insulin in mice. Methods: Upon systemic injection of 125I-Aβ40, 125I-Aβ42, or 125I-insulin, the plasma pharmacokinetics and brain influx were assessed in wild-type (WT) or AD transgenic (APP/PS1) mice at various ages. Additionally, publicly available single-cell RNA-Seq data [GSE129788] was employed to investigate pathways regulating BBB transport in WT mice at different ages. Results: The brain influx of 125I-Aβ40, estimated as the permeability-surface area product, decreased with age, accompanied by an increase in plasma AUC. In contrast, the brain influx of 125I-Aβ42 increased with age, accompanied by a decrease in plasma AUC. The age-dependent changes observed in WT mice were accelerated in APP/PS1 mice. As seen with 125I-Aβ40, the brain influx of 125I-insulin decreased with age in WT mice, accompanied by an increase in plasma AUC. This finding was further supported by dynamic single-photon emission computed tomography (SPECT/CT) imaging studies. RAGE and PI3K/AKT signaling pathways at the BBB, which are implicated in Aβ and insulin transcytosis, respectively, were upregulated with age in WT mice, indicating BBB insulin resistance. Conclusion: Aging differentially affects the plasma pharmacokinetics and brain influx of Aβ isoforms and insulin in a manner that could potentially augment AD risk.

Objective To test the ability of machine learning (ML) approaches with clinical and genomic biomarkers to predict methotrexate treatment response in patients with early rheumatoid arthritis (RA). Methods Demographic, clinical, and genomic data from 643 patients of European ancestry with early RA (mean age 54 years; 70% female) subdivided into a training (n = 336) and validation cohort (n = 307) were used. The genomic data comprised 160 single‐nucleotide polymorphisms (SNPs) previously associated with RA or methotrexate metabolism. Response to methotrexate monotherapy was defined as good or moderate by the European Alliance of Associations for Rheumatology (EULAR) response criteria at the 3‐month follow‐up. Supervised ML methods were trained with 5 repeats and 10‐fold cross‐validation using the training cohort. Prediction performance was validated in the independent validation cohort. Results Supervised ML methods combining age, sex, smoking, rheumatoid factor, baseline Disease Activity Score in 28 joints (DAS28) scores and 160 SNPs predicted EULAR response at 3 months with the area under the receiver operating curve of 0.84 ( P = 0.05) in the training cohort and achieved a prediction accuracy of 76% ( P = 0.05) in the validation cohort (sensitivity 72%, specificity 77%). Intergenic SNPs rs12446816, rs13385025, rs113798271, and ATIC (rs2372536) had variable importance above 60.0 and along with baseline DAS28 scores were among the top predictors of methotrexate response. Conclusion Pharmacogenomic biomarkers combined with baseline DAS28 scores can be useful in predicting response to methotrexate in patients with early RA. Applying ML to predict treatment response holds promise for guiding effective RA treatment choices, including timely escalation of RA therapies.

Abstract Our previous matched case–control study of postmenopausal women with resected early-stage breast cancer revealed that only anastrozole, but not exemestane or letrozole, showed a significant association between the 6-month estrogen concentrations and risk of breast cancer. Anastrozole, but not exemestane or letrozole, is a ligand for estrogen receptor α. The mechanisms of endocrine resistance are heterogenous and with the new mechanism of anastrozole, we have found that treatment of anastrozole maintains fatty acid synthase (FASN) protein level by limiting the ubiquitin-mediated FASN degradation, leading to increased breast cancer cell growth. Mechanistically, anastrozole decreases the guided entry of tail-anchored proteins factor 4 (GET4) expression, resulting in decreased BCL2-associated athanogene cochaperone 6 (BAG6) complex activity, which in turn, prevents RNF126-mediated degradation of FASN. Increased FASN protein level can induce a negative feedback loop mediated by the MAPK pathway. High levels of FASN are associated with poor outcome only in patients with anastrozole-treated breast cancer, but not in patients treated with exemestane or letrozole. Repressing FASN causes regression of breast cancer cell growth. The anastrozole-FASN signaling pathway is eminently targetable in endocrine-resistant breast cancer.

Abstract The role of maternal microbiome transmitted at birth in cancer control is poorly understood. We have developed the first germfree B6 mouse model of breast cancer to investigate the role of the maternal microbiome in controlling oncogenic/metastatic frequencies of pro-oncogenic mammary cells. In this model, a DNA barcoded, primitive normal mouse mammary epithelial cell encoding MMTV-PyMT oncogene, was transplanted in large numbers into conventional or germfree B6 mice. Next-gen sequencing analysis of the DNA barcodes in tissues enabled us to clonally track millions of cells and measure the frequency and growth dynamics of clones at the primary site and their systemic distribution in circulation and all vital organs, generating an unprecedented high-definition map of cancer progression. Our results show that in conventional B6 mice with maternal microbial transmission at-birth, a small fraction (~0.01%) of transplanted cells transform de novo and produce slow growing, late-onset benign tumors (mouse median survival of > 1 year). In contrast, in germfree B6 mice, a >10-fold higher frequency of cells transform de novo and generate early-onset, highly aggressive metastatic clones, and are frequently associated with features leading to early euthanasia or endpoint (i.e. sudden death, intracardiac metastasis, paralysis, swollen abdomen, early multiorgan aggressive metastasis) (median survival of ~4 months; p<0.005 vs conventional mice). However, postnatal intervention with healthy fecal implants reconstituted gut microbiota in germfree B6 mice and significantly restored protective effects against cancer progression and sudden death, as well as improved quality/duration of life (median survival of >1 year; p<0.005 vs germfree mice; p=ns vs conventional mice). We found maternal microbial transmission fully protected against systemic metastasis. GC-MS analysis identified significantly distinct primary tumor metabolome associated with conventional (non-metastatic) mice vs germfree (metastatic) mice. Metagenomic analysis of fecal samples from a subset of tumor bearing conventional (non-metastatic) mice vs fecal implanted (late-onset, less aggressive metastatic) germfree mice identified significant association between metastasis and altered microbial communities. These findings using novel germfree mouse model of metastasis and cellular DNA barcoding technology provide first direct evidence that diverse microbial communities transmitted at birth play a critical role in protecting the host from clonal de novo mammary cell transformation and progression to systemic disease and thereby improving quality and duration of life. Citation Format: Nagarajan Kannan, Syed Mohammed Musheer Aalam, Xiaojia Tang, Krishna Kalari, Purna Kashyap, Jun Chen, Stephen Johnson, Anguraj Sadanandam, Mark Sherman. Maternal microbiome protects host against clonal de novo transformation, early onset systemic metastasis, and sudden death [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr LB226.

Thiopurine drugs such as 6‐mercaptopurine (6‐MP) and 6‐thioguanine (6‐TG) are effective immunosuppressants and antineoplastic agents. We used a genome‐wide approach to identify genes that might contribute to individual variation in thiopurine response. Basal gene expression profiles, genome‐wide SNP data, and 6‐MP and 6‐TG Cell growth inhibition (GI50) values were generated for 200 Coriell Institute “Human Variation Panel” lymphoblastoid cell lines. Gene expression and genome‐wide SNPs were correlated with GI50 values to identify candidate genes that might contribute to thiopurine response. We identified 104 expression probes and 161 SNPs, as well as 315 gene probes and 579 SNPs that correlated significantly with GI50 values for 6‐MP and 6‐TG, respectively (P<10 −5 and 10 −6 , respectively). Three genes, CD83 , TXNDC13 and SNAP23 , contained SNPs that contributed significantly to both 6‐MP and 6‐TG cytotoxicity and were also significantly associated with expression levels (P<10 −4 ), suggesting a “three‐way model” involving genotype, gene expression and cytotoxicity. GBP5 , a gene with an expression level that correlated negatively with 6‐MP and 6‐TG GI50 values, sensitized leukemia cells when overexpressed. Functional studies of the remainder of the candidate genes are being performed. These results may enhance our ability to predict thiopurine response phenotypes. NIH GM28157.

Objectives Betaine‐homocysteine methyltransferase (BHMT) catalyzes the remethylation of homocysteine. BHMT in the human liver is responsible for up to 25% of homocysteine metabolism. Knowledge of the relationship between genetic polymorphisms and BHMT phenotypes could enhance our understanding of the role of this reaction in drug metabolism. Methods Levels of BHMT enzyme activity and immunoreactive protein were measured in 268 human liver biopsies. Twelve tag BHMT SNPs from our own gene resequencing and HapMap data were used to genotype DNA from these samples. Genotype‐phenotype association analysis was then performed. BHMT enzyme activity averaged 0.077 ± 0.028 units/μg, with a 12 fold variation. Results BHMT protein levels averaged 6.9 ± 3.3 units/μg, with a 50‐fold variation. Three SNPs (rs41272270, rs16876512, and rs6875201), located 28kb upstream, in the 5′‐UTR and in intron 1 of BHMT , respectively, were significantly correlated with both phenotypes (p<10 −8 for enzyme activity and p<10 −5 for protein level). rs41272270 was associated with altered protein binding during electrophoresis mobility shift assay. Conclusions Three tag SNPs were significantly associated with human hepatic BHMT enzyme activity and immunoreactive protein. These observations may help us understand the regulation of BHMT in the human liver and its role in drug metabolism. Supported by NIH U01 GM61388, R01 GM28157, R01 CA132780.