Kubo, Michiaki

Abstract Purpose: Germline pathogenic variants in CHEK2 confer moderately elevated breast cancer risk (odds ratio, OR ∼ 2.5), qualifying carriers for enhanced breast cancer screening. Besides pathogenic variants, dozens of missense CHEK2 variants of uncertain significance (VUS) have been identified, hampering the clinical utility of germline genetic testing (GGT). Experimental Design: We collected 460 CHEK2 missense VUS identified by the ENIGMA consortium in 15 countries. Their functional characterization was performed using CHEK2-complementation assays quantifying KAP1 phosphorylation and CHK2 autophosphorylation in human RPE1–CHEK2-knockout cells. Concordant results in both functional assays were used to categorize CHEK2 VUS from 12 ENIGMA case–control datasets, including 73,048 female patients with breast cancer and 88,658 ethnicity-matched controls. Results: A total of 430/460 VUS were successfully analyzed, of which 340 (79.1%) were concordant in both functional assays and categorized as functionally impaired (N = 102), functionally intermediate (N = 12), or functionally wild-type (WT)–like (N = 226). We then examined their association with breast cancer risk in the case–control analysis. The OR and 95% CI (confidence intervals) for carriers of functionally impaired, intermediate, and WT-like variants were 2.83 (95% CI, 2.35–3.41), 1.57 (95% CI, 1.41–1.75), and 1.19 (95% CI, 1.08–1.31), respectively. The meta-analysis of population-specific datasets showed similar results. Conclusions: We determined the functional consequences for the majority of CHEK2 missense VUS found in patients with breast cancer (3,660/4,436; 82.5%). Carriers of functionally impaired missense variants accounted for 0.5% of patients with breast cancer and were associated with a moderate risk similar to that of truncating CHEK2 variants. In contrast, 2.2% of all patients with breast cancer carried functionally wild-type/intermediate missense variants with no clinically relevant breast cancer risk in heterozygous carriers.

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.

Anastrozole is a widely prescribed aromatase inhibitor for the therapy of estrogen receptor positive (ER+) breast cancer. We performed a genome‐wide association study (GWAS) for plasma anastrozole concentrations in 687 postmenopausal women with ER+ breast cancer. The top single‐nucleotide polymorphism (SNP) signal mapped across SLC38A7 (rs11648166, P = 2.3E‐08), which we showed to encode an anastrozole influx transporter. The second most significant signal (rs28845026, P = 5.4E‐08) mapped near ALPPL2 and displayed epistasis with the SLC38A7 signal. Both of these SNPs were cis expression quantitative trait loci (eQTL)s for these genes, and patients homozygous for variant genotypes for both SNPs had the highest drug concentrations, the highest SLC38A7 expression, and the lowest ALPPL2 expression. In summary, our GWAS identified a novel gene encoding an anastrozole transporter, SLC38A7, as well as epistatic interaction between SNPs in that gene and SNPs near ALPPL2 that influenced both the expression of the transporter and anastrozole plasma concentrations.

Background: Coronary heart disease (CHD) is a leading cause of death globally. Although therapy with statins decreases circulating levels of low-density lipoprotein cholesterol and the incidence of CHD, additional events occur despite statin therapy in some individuals. The genetic determinants of this residual cardiovascular risk remain unknown. Methods: We performed a 2-stage genome-wide association study of CHD events during statin therapy. We first identified 3099 cases who experienced CHD events (defined as acute myocardial infarction or the need for coronary revascularization) during statin therapy and 7681 controls without CHD events during comparable intensity and duration of statin therapy from 4 sites in the Electronic Medical Records and Genomics Network. We then sought replication of candidate variants in another 160 cases and 1112 controls from a fifth Electronic Medical Records and Genomics site, which joined the network after the initial genome-wide association study. Finally, we performed a phenome-wide association study for other traits linked to the most significant locus. Results: The meta-analysis identified 7 single nucleotide polymorphisms at a genome-wide level of significance within the LPA/PLG locus associated with CHD events on statin treatment. The most significant association was for an intronic single nucleotide polymorphism within LPA/PLG (rs10455872; minor allele frequency, 0.069; odds ratio, 1.58; 95% confidence interval, 1.35–1.86; P =2.6×10 − 10 ). In the replication cohort, rs10455872 was also associated with CHD events (odds ratio, 1.71; 95% confidence interval, 1.14–2.57; P =0.009). The association of this single nucleotide polymorphism with CHD events was independent of statin-induced change in low-density lipoprotein cholesterol (odds ratio, 1.62; 95% confidence interval, 1.17–2.24; P =0.004) and persisted in individuals with low-density lipoprotein cholesterol ≤70 mg/dL (odds ratio, 2.43; 95% confidence interval, 1.18–4.75; P =0.015). A phenome-wide association study supported the effect of this region on coronary heart disease and did not identify noncardiovascular phenotypes. Conclusions: Genetic variations at the LPA locus are associated with CHD events during statin therapy independently of the extent of low-density lipoprotein cholesterol lowering. This finding provides support for exploring strategies targeting circulating concentrations of lipoprotein(a) to reduce CHD events in patients receiving statins.

Background The majority of hypertensive individuals require combination antihypertensive therapy to achieve adequate blood pressure ( BP ) control. This study aimed to identify genetic variants associated with uncontrolled BP on combination therapy with a thiazide diuretic and a β‐blocker. Methods and Results A genome‐wide association study of uncontrolled BP on combination therapy was conducted among 314 white participants of the PEAR (Pharmacogenomic Evaluation of Antihypertensive Responses) trial. Multivariable logistic regression analysis was used. Genetic variants meeting a suggestive level of significance ( P <1.0E‐05) were tested for replication in an external cohort, INVEST (International Verapamil‐ SR Trandolapril study). We also examined genome‐wide variant associations with systolic and diastolic BP response on combination therapy and tested for replication. We discovered a single nucleotide polymorphism, the rs261316 major allele, at chromosome 15 in the gene ALDH1A2 associated with an increased odds of having uncontrolled BP on combination therapy (odds ratio: 2.56, 95% confidence interval, 1.69–3.88, P =8.64E‐06). This single nucleotide polymorphism replicated (odds ratio: 1.86, 95% confidence interval, 1.35–2.57, P =0.001) and approached genome‐wide significance in the meta‐analysis between discovery and replication cohorts (odds ratio: 2.16, 95% confidence interval, 1.63–2.86, P =8.60E‐08). Other genes in the region surrounding rs261316 ( ALDH1A2) include AQP9 and LIPC . Conclusions A single nucleotide polymorphism in the gene ALDH1A2 may be associated with uncontrolled BP following treatment with a thiazide diuretic/β‐blocker combination. Clinical Trial Registration URL: https://www.clinicaltrials.gov . Unique identifier: NCT 00246519.

Abstract BACKGROUND: Based on prospective clinical trials, there is no evidence for differences in efficacy between the 3 aromatase inhibitors (AIs) anastrozole, exemestane, and letrozole. The purpose of this study was to identify germline genetic variants associated with response to AIs and to help identify novel mechanisms associated with drug disease efficacy. METHODS: A genome-wide association study (GWAS) was performed for 624 patients (Steroids 2015;99:32-38) to identify SNPs associated with estrogen level change in women with estrogen receptor (ER) positive breast cancer treated with anastrozole. Replication of associated SNPs was performed in a GWAS from the MA.27 trial that compared adjuvant anastrozole and exemestane treatment of post-menopausal women with ER+ breast cancer. Functional studies were subsequently performed to determine SNP effects and underlying mechanisms. RESULTS: Our initial GWAS identified SNPs within CSMD1 that were associated with changes in estrogen levels during anastrozole therapy. An additional SNP in CSMD1 was also associated with breast cancer events in CCTG MA.27. Functionally, we showed that CSMD1 regulates CYP19 expression in a SNP-, and in an anastrozole- dependent fashion. These phenomena were not observed for either letrozole or exemestane. In MA.27, an anastrozole- specific effect was also seen with the minor allele having a protective effect on time to distant metastasis (HR=0.49, p=0.00259), but this was not the case for exemestane (HR=0.71, p=0.111). Our in vitro functional studies indicated that overexpression of CSMD1 sensitized anastrozole or letrozole resistant cells to anastrozole but not to the other two AIs. The SNP in CSMD1 that was associated with increased CSMD1 and CYP19 expression levels increased anastrozole sensitivity, but not letrozole or exemestane in lymphoblastoid cell lines (LCLs) homozygous for either WT or variant CSMD1 SNP genotypes. Based on these observations, we explored whether anastrozole has additional mechanisms beyond its function as a CYP19 inhibitor. Utilizing an estrogen response element (ERE) luciferase reporter assay in a CYP19 CRISPR knockout breast cancer T47D cell line and a surface plasmon resonance (SPR) assay, we found that anastrozole can also function as an ERα agonist, and can bind to, and result in, proteasome dependent ERα degradation, especially in the presence of E2. Treatment of these CYP19 CRISPR knockout cells with anastrozole in the presence of increasing concentrations of E2 results in greater sensitivity compared with anastrozole alone, while the addition of E2, as expected, does not improve letrozole or exemestane sensitivity. These same observations were also seen in letrozole and anastrazole resistant cells. CONCLUSIONS: Our findings suggest that anastrozole might be more effective than letrozole or exemestane in patients with the CSMD1 SNP. Furthermore, anastrozole can function as an ERα agonist, binding to ERα and resulting in its degradation, especially in the presence of E2. These findings should help to make it possible to develop precision endocrine therapies for women who are candidates for AIs. Citation Format: Cairns J, Ingle J, Dudenkov T, Kalari K, Buzdar A, Kubo M, Robson M, Ellis M, Goss P, Shepherd L, Goetz M, Weinshilboum R, Wang L. CSMD1 SNPs selectively affect anastrozole response in postmenopausal breast cancer patients [abstract]. In: Proceedings of the 2016 San Antonio Breast Cancer Symposium; 2016 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2017;77(4 Suppl):Abstract nr PD1-04.

Abstract Genetic risks in breast cancer remain only partly understood. Here, we report the results of a genome-wide association study of germline DNA from 4,658 women, including 252 women experiencing a breast cancer recurrence, who were entered on the MA.27 adjuvant trial comparing the aromatase inhibitors (AI) anastrozole and exemestane. Single-nucleotide polymorphisms (SNP) of top significance were identified in the gene encoding MIR2052HG, a long noncoding RNA of unknown function. Heterozygous or homozygous individuals for variant alleles exhibited a ∼40% or ∼63% decrease, respectively, in the hazard of breast cancer recurrence relative to homozygous wild-type individuals. Functional genomic studies in lymphoblastoid cell lines and ERα-positive breast cancer cell lines showed that expression from MIR2052HG and the ESR1 gene encoding estrogen receptor-α (ERα) was induced by estrogen and AI in a SNP-dependent manner. Variant SNP genotypes exhibited increased ERα binding to estrogen response elements, relative to wild-type genotypes, a pattern that was reversed by AI treatment. Further, variant SNPs were associated with lower expression of MIR2052HG and ERα. RNAi-mediated silencing of MIR2052HG in breast cancer cell lines decreased ERα expression, cell proliferation, and anchorage-independent colony formation. Mechanistic investigations revealed that MIR2052HG sustained ERα levels both by promoting AKT/FOXO3-mediated ESR1 transcription and by limiting ubiquitin-mediated, proteasome-dependent degradation of ERα. Taken together, our results define MIR2052HS as a functionally polymorphic gene that affects risks of breast cancer recurrence in women treated with AI. More broadly, our results offer a pharmacogenomic basis to understand differences in the response of breast cancer patients to AI therapy. Cancer Res; 76(23); 7012–23. ©2016 AACR.

Background Intracranial aneurysms ( IAs ), abdominal aortic aneurysms ( AAAs ), and thoracic aortic aneurysms ( TAAs ) all have a familial predisposition. Given that aneurysm types are known to co‐occur, we hypothesized that there may be shared genetic risk factors for IAs , AAAs, and TAAs . Methods and Results We performed a mega‐analysis of 1000 Genomes Project‐imputed genome‐wide association study ( GWAS ) data of 4 previously published aneurysm cohorts: 2 IA cohorts (in total 1516 cases, 4305 controls), 1 AAA cohort (818 cases, 3004 controls), and 1 TAA cohort (760 cases, 2212 controls), and observed associations of 4 known IA , AAA , and/or TAA risk loci (9p21, 18q11, 15q21, and 2q33) with consistent effect directions in all 4 cohorts. We calculated polygenic scores based on IA ‐, AAA ‐, and TAA ‐associated SNP s and tested these scores for association to case‐control status in the other aneurysm cohorts; this revealed no shared polygenic effects. Similarly, linkage disequilibrium–score regression analyses did not show significant correlations between any pair of aneurysm subtypes. Last, we evaluated the evidence for 14 previously published aneurysm risk single‐nucleotide polymorphisms through collaboration in extended aneurysm cohorts, with a total of 6548 cases and 16 843 controls ( IA ) and 4391 cases and 37 904 controls ( AAA ), and found nominally significant associations for IA risk locus 18q11 near RBBP 8 to AAA (odds ratio [ OR ]=1.11; P =4.1×10 −5 ) and for TAA risk locus 15q21 near FBN 1 to AAA ( OR =1.07; P =1.1×10 −3 ). Conclusions Although there was no evidence for polygenic overlap between IAs , AAAs , and TAAs , we found nominally significant effects of two established risk loci for IAs and TAAs in AAAs . These two loci will require further replication.