Etheridge, Amy S.

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Abstract Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Abstract Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Enhancers are critical for regulating tissue-specific gene expression, and genetic variants within enhancer regions have been suggested to contribute to various cancer-related processes, including therapeutic resistance. However, the precise mechanisms remain elusive. Using a well-defined drug-gene pair, we identified an enhancer region for dihydropyrimidine dehydrogenase (DPD, DPYD gene) expression that is relevant to the metabolism of the anti-cancer drug 5-fluorouracil (5-FU). Using reporter systems, CRISPR genome-edited cell models, and human liver specimens, we demonstrated in vitro and vivo that genotype status for the common germline variant (rs4294451; 27% global minor allele frequency) located within this novel enhancer controls DPYD transcription and alters resistance to 5-FU. The variant genotype increases recruitment of the transcription factor CEBPB to the enhancer and alters the level of direct interactions between the enhancer and DPYD promoter. Our data provide insight into the regulatory mechanisms controlling sensitivity and resistance to 5-FU.

Abstract Genetic variants in the DPYD gene that are deleterious to dihydropyrimidine dehydrogenase (DPD) enzyme function have demonstrated clinical validity as predictive biomarkers for severe (i.e., CTCAE grade ≥3) and potentially lethal toxicity to fluorouracil (5-FU) and capecitabine chemotherapies. Variants within the coding region only explain 5-10% of 5-FU-related toxicities; and GWAS approaches have been unable to identify validated predictive biomarkers to date. Our previous studies demonstrated that DPD expression was controlled by multiple epigenetic processes; considerable variation was also noted in liver DPYD/DPD levels. Using CRISPR-interference (CRISPRi) and CRISPR-activation (CRISPRa), we identified functional enhancer regions within liver tissues and cells. Epigenetic status at these regions was shown to correlate with DPYD expression in cell models and human liver tissues. Chromosome conformation capture (3C) confirmed direct promoter-enhancer interactions. A multi-omics analysis of >1000 human liver specimens identified a proximal eQTL for DPYD expression within an enhancer, with a global minor allele frequency of 27%. This variant was subsequently shown to directly alter DPYD expression in knock-in cell line models. The variant disrupts transcription factor binding at the DPYD enhancer, leading to reduced enhancer-promoter interaction, reduced DPD expression, and increased cellular sensitivity to 5-FU. Collectively, our data suggest that the variant could represent a novel predictive biomarker for severe 5-FU toxicity. The evaluation of this marker in the clinical setting is currently underway. Citation Format: Ting Zhang, Alisa Ambrodji, Huixing Huang, Kelly Bouchonville, Amy Etheridge, Remington Schmidt, Jose Cardiel Nunez, Zoey Temesgen, Federico Innocenti, Robert Diasio, Carlo Largiadèr, Steven M. Offer. DPYD enhancer variant confers increased risk of severe toxicity in 5-FU chemotherapy [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 5420.

Hypertension is a common bevacizumab‐induced toxicity. No markers are available to predict patients at risk of developing hypertension. We hypothesized that genetic risk of essential hypertension, as measured by a blood pressure polygenic risk score (PRS), would be associated with risk of severe bevacizumab‐induced hypertension. PRSs were calculated for 1,027 bevacizumab‐treated patients of European descent with cancer from four clinical trials (Alliance for Clinical Trials in Oncology (Alliance) / Cancer and Leukemia Group B (CALGB) 80303, 40503, 90401, 40502) using summary systolic blood pressure (SBP) and diastolic blood pressure (DBP) genome‐wide association results obtained from 757,601 individuals of European descent. The association between PRS and grade 3 bevacizumab‐induced hypertension (Common Toxicity Criteria for Adverse Events version 3) in each trial was performed by multivariable logistic regression. Fixed‐effect meta‐analyses odds ratios (ORs) per standard deviation (SD) of the association of PRS (quantitative) and hypertension across trials were estimated by inverse‐variance weighting. PRSs were additionally stratified into quintiles, with the bottom quintile as the referent group. The OR of the association between hypertension and each quintile vs. the referent group was determined by logistic regression. The most significant PRS (quantitative)‐hypertension association included up to 67 single‐nucleotide variants (SNPs) associated with SBP (P = 0.0077, OR per SD = 1.31, 95% confidence interval (CI), 1.07–1.60), and up to 53 SNPs associated with DBP (P = 0.0209, OR per SD = 1.27, 95% CI, 1.04–1.56). Patients in the top quintile had a higher risk of developing bevacizumab‐induced hypertension compared with patients in the bottom quintile using SNPs associated with SBP (P = 4.75 × 10−4, OR = 3.72, 95% CI, 1.84–8.16) and DBP (P = 0.076, OR = 1.83, 95% CI, 0.95–3.64). Genetic variants associated with essential hypertension, mainly SBP, increase the risk of severe bevacizumab‐induced hypertension.

3030 Background: Diarrhea, HFS, and hypertension are common toxicities of sorafenib. No markers are validated to predict patients at risk of these toxicities. This study aimed to identify genetic predictors of sorafenib-induced toxicities. Methods: A two-step, discovery-validation approach was used. The discovery set included 140 renal cell carcinoma patients from the TARGET study treated with sorafenib (400 mg twice daily) and genotyped for 1040 single-nucleotide polymorphisms (SNPs) in 56 genes. The three most statistically significant SNPs associated with grade ≥2 composite toxicity (either hypertension, diarrhea, HFS, or other skin toxicities, CTCAE v.3.0) were tested for association with grade 3 composite toxicity (either hypertension, diarrhea, or HFS, CTCAE v.4.0) in a validation set of 240 hepatocellular carcinoma patients from Alliance/CALGB 80802 treated with sorafenib (400 mg twice daily) alone or with doxorubicin. Associations between SNPs and composite toxicity was performed by logistic regression, with adjusting covariates (age, gender, race, and treatment arm, the latter two covariates for the validation set only). A meta-analysis odds ratio (OR) of each SNP-grade 3 toxicity association between the discovery and validation sets was obtained by inverse variance to point toward effects specific to a type of toxicity. Results: In the discovery set, the top three SNPs associated with grade ≥2 composite toxicity were rs12366035 (C>T, minor allele frequency, MAF 0.34) in VEGFB (p 0.0007), rs4035887 (G>A, MAF 0.49) in EPAS1 (p 0.0021), and rs4864950 (T>A, MAF 0.23) in KDR (p 0.0058). These SNPs were genotyped in the validation set and only rs4864950 in KDR was replicated. No grade 4 toxicities were reported. Similar to the discovery set (OR 2.41, 95% CI 1.29-4.51), the A allele of rs4864950 increased the risk of grade 3 composite toxicity (p 0.032, OR 2.12, 95% CI 1.70-4.27) in the validation set. Grade 3 toxicity prevalence in the discovery and validation sets were 3.6% and 7.4% diarrhea, 8.6% and 12.3% HFS, 3.6% and 8.8% hypertension, respectively. The meta-analysis of the two datasets showed that the A allele of rs4864950 increased the risk of grade 3 diarrhea (p 0.045, OR 3.09, 95% CI 1.03-9.29), grade 3 HFS (p 0.012, OR 2.57, 95% CI 1.24-5.37), but not grade 3 hypertension (p 0.207, OR 0.51, 95% CI 0.18-1.45). Conclusions: We provide the first evidence of clinical validity of a marker of sorafenib-induced diarrhea and HFS. Sorafenib inhibits VEGFR2 (coded by KDR), leading to epithelial hypoxia and causing diarrhea and HFS. Variant rs4864950 might affect the function VEGFR2, which, during VEGFR2 inhibition, increases the risk of diarrhea and HFS. This SNP is common and can be genotyped in patients before receiving sorafenib for a better risk assessment. Support: U10CA180821, U10CA180882, U24CA196171; https://acknowledgments.alliancefound.org ClinicalTrials.gov Id: NCT01015833.

Abstract Purpose: Irinotecan/5-fluorouracil (5-FU; FOLFIRI) or oxaliplatin/5-FU (FOLFOX), combined with bevacizumab or cetuximab, are approved, first-line treatments for metastatic colorectal cancer (mCRC). We aimed at identifying germline variants associated with survival in patients with mCRC treated with these regimens in Cancer and Leukemia Group B/SWOG 80405. Experimental Design: Patients with mCRC receiving either FOLFOX or FOLFIRI were randomized to either cetuximab or bevacizumab. DNA from peripheral blood was genotyped for approximately 700,000 SNPs. The association between SNPs and overall survival (OS) was tested in 613 patients of genetically estimated European ancestry using Cox proportional hazards models. Results: The four most significant SNPs associated with OS were three haplotypic SNPs between microsomal glutathione S-transferase 1 (MGST1) and LIM domain only 3 (LMO3, representative HR, 1.56; P = 1.30 × 10−6), and rs11644916 in AXIN1 (HR, 1.39, P = 4.26 × 10−6). AXIN1 is a well-established tumor suppressor gene in colorectal cancer, and rs11644916 (G>A) conferred shorter OS. Median OS for patients with the AA, AG, or GG genotypes was 18.4, 25.6, or 36.4 months, respectively. In 90 patients with stage IV colorectal cancer from The Cancer Genome Atlas (TCGA), rs11649255 in AXIN1 [in almost complete linkage disequilibrium (LD) with rs11644916], was associated with shorter OS (HR, 2.24, P = 0.0096). Using rs11648673 in AXIN1 (in very high LD with rs11644916 and with functional evidence), luciferase activity in three colorectal cancer cell lines was reduced. Conclusions: This is the first large genome-wide association study ever conducted in patients with mCRC treated with first-line standard treatment in a randomized phase III trial. A common SNP in AXIN1 conferred worse OS and the effect was replicated in TCGA. Further studies in colorectal cancer experimental models are required.