AbstractGene‐set analyses have been widely used in gene expression studies, and some of the developed methods have been extended to genome wide association studies (GWAS). Yet, complications due to linkage disequilibrium (LD) among single nucleotide polymorphisms (SNPs), and variable numbers of SNPs per gene and genes per gene‐set, have plagued current approaches, often leading to ad hoc “fixes.” To overcome some of the current limitations, we developed a general approach to scan GWAS SNP data for both gene‐level and gene‐set analyses, building on score statistics for generalized linear models, and taking advantage of the directed acyclic graph structure of the gene ontology when creating gene‐sets. However, other types of gene‐set structures can be used, such as the popular Kyoto Encyclopedia of Genes and Genomes (KEGG). Our approach combines SNPs into genes, and genes into gene‐sets, but assures that positive and negative effects of genes on a trait do not cancel. To control for multiple testing of many gene‐sets, we use an efficient computational strategy that accounts for LD and provides accurate step‐down adjustedP‐values for each gene‐set. Application of our methods to two different GWAS provide guidance on the potential strengths and weaknesses of our proposed gene‐set analyses.
SULT1A1
,
CYP2C19
and Disease-Free Survival in early Breast Cancer Patients Receiving Tamoxifen
Evaluation of CYP2D6 and Efficacy of Tamoxifen and Raloxifene in Women Treated for Breast Cancer Chemoprevention: Results from the NSABP P1 and P2 Clinical Trials
Abstract Background: Mammographic density, the variation in fat, epithelial and stromal tissues seen on screening mammography, is a strong risk factor for breast cancer and can be modified by hormonal agents. Changes in density from tamoxifen or postmenopausal hormone (PMH) use are associated with risk, suggesting that density may be a surrogate marker of therapeutic efficacy. Aromatase inhibitors (AIs) are given as adjuvant therapy in hormone receptor positive postmenopausal breast cancer and are known to decrease levels of estrone and estradiol in both serum and breast tissue. Our goal here was to examine the influence of AIs on mammographic density in women with early breast cancer. Methods: We conducted a case-control study of postmenopausal breast cancer patients initiating adjuvant AI therapy (anastrozole or exemestane) on protocols NCIC CTG MA27, NCCTG N063I and MC (Mayo Clinic) 0532. Eligibility included; an intact contralateral breast with no prior surgery; a screening mammogram within twelve months before AI initiation and at 9-15 months on therapy; no prior endocrine therapy and informed consent. Controls were sampled from the Mayo Mammography Health Study, a cohort of 19,924 receiving screening mammography at the Mayo Clinic, and matched to cases on age, prior PMH use, baseline body mass index (BMI) and interval between mammograms. Pre-treatment and on-study mammograms for cases (corresponding mammograms for controls) were digitized. Change in percent density was estimated on the craniocaudal view of the non-cancerous breast using two methods: a subjective assessment of change by an expert radiologist (within 5%; 5-10% increase, 10-25% increase, 25%+ increase, 5-10% decrease, 10-25% decrease and 25%+ decrease) and a quantitative assessment of absolute change using a computer-assisted thresholding program (Cumulus). Analyses compared magnitude of change in density by both the subjective and quantitative methods between cases and matched controls. Results: 574 pairs were eligible for analyses (MA27-505 cases; N063I-12 cases; MC0532-57 cases). Characteristics of the two groups are shown in the table below. Using either density estimation method, there was a greater decrease in density among women on AI therapy vs. matched controls. In 33% (95% CI: 29-37%) of pairs, there was at least a one greater category decrease for the case relative to her control by subjective estimation. In 14% (95% CI: 11-18%) of the pairs, there was at least a 5% greater decrease for the case relative to her control by quantitative estimation. Data will be available according to AI class (non-steroidal versus steroidal) in November. Conclusions: In the largest report to date to examine the influence of AI therapy on mammographic density, we provide evidence that AI is associated with decreases in density in a small subgroup of women. We are currently examining factors that influence these AI-associated decreases in density and whether these differences are unique to one class of AI. (Supported in part by NIH grants P50CA116201, U01GM61388, U10CA77202, U10CA25224) Citation Information: Cancer Res 2010;70(24 Suppl):Abstract nr P2-09-03.
A Genome-Wide Association Study in Patients Experiencing Musculoskeletal Adverse Events on Aromatase Inhibitors as Adjuvant Therapy in Early Breast Cancer Entered on NCIC CTG Trial MA.27. A Pharmacogenetics Research Network-RIKEN Collaboration.
Abstract Background: The third-generation aromatase inhibitors (AIs) are established adjuvant therapy for postmenopausal women with early stage breast cancer. Musculoskeletal adverse events (MS AE) have become increasingly recognized as an issue for many women receiving AIs resulting in decreased quality of life and, in some instances, discontinuation of the AI. MA.27 is an ongoing randomized clinical trial comparing anastrozole with exemestane that accrued 7576 patients (pts), including 6827 from North America, about 79% of whom provided DNA and consent for genetic testing. A review of withdrawals related to adverse events revealed that the most common cause was MS AE. The hypothesis for our study was that a genome-wide association (GWA) case control study would identify single nucleotide polymorphisms (SNPs) associated with MS AE.Methods: Cases were defined as pts with grade 3 or 4 MS AE (NCI CTCAE v3.0) or who went off treatment for any grade MS AE within the first 2 years and were matched to 2 controls based on treatment arm, length of follow-up, prior chemotherapy, age and prior celecoxib (at least 3 months prior to MS AE). The study was restricted to the 94% of pts self-identified as Caucasian. Genotyping was performed with the Illumina Human610-Quad. Genotype quality control prior to analysis required call rates ≥98%. Eigenstrat analyses were performed to control for population stratification.Results: The GWA study included 293 cases and 586 controls with 6 pts (0.7%) (4 cases, 2 controls) excluded for call rates <98%. 551,358 SNPs were used in the analyses after removing 29,439 for a minor allele frequency <1% and 82 with Hardy Weinberg Equilibrium p-values <10-6. Eigenvectors did not impact the results when used as covariates. Four SNPs with the lowest p-values (3.3-8.2 x 10-6, Armitage test), 3 in high LD (r2>.8) on chromosome (ch) 14 (rs7158782, rs7159713, rs2369049) and 1 on the X ch (rs6637820) were identified. The gene closest (4000-7000 bp) to the 3 SNPs on ch 14 was T-Cell Leukemia 1A (TCL1A), encoding a protein known to augment AKT kinase activity. Expression array data from lymphoblastoid cell lines from 300 subjects of 3 ethnicities revealed that all 3 SNPs were associated with decreased TCL1A expression after adjusting for ethnicity using 2 probe sets (p=0.006 to 0.03). Electromobility shift assay (EMSA) showed that 2 of these SNPs (rs7158782, rs7159713) displayed a shift and the former SNP showed substantially less binding for variant than wild type. The gene closest to the X ch SNP was Immunoglobulin Superfamily 1 (IGSF1), but EMSA did not show a shift for this SNP.Conclusions: MS AEs represent a major impediment to optimal use of AIs in women with breast cancer and this GWA study identified SNPs on ch 14 that provide a focus for further research to identify pts at risk for, and means to ameliorate, this adverse event.(Supported in part by NIH grants U01GM61388, U01GM63173, P50CA116201 and U10CA77202) Citation Information: Cancer Res 2009;69(24 Suppl):Abstract nr 15.
Pharmacogenetic (CYP2D6) and gene expression profiles (HOXB13/IL17BR and molecular grade index) for prediction of adjuvant endocrine therapy benefit in the ABCSG 8 trial.
Abstract Abstract #57 Background: ABCSG 8 is the only prospective trial examining the value of switching patients (pts) to an aromatase inhibitor after tamoxifen (Tam) in which randomization occurred prior to therapy. This trial demonstrated that 2 years (yrs) of Tam followed by 3 yrs of anastrozole (TAM to A) was superior to 5 yrs of Tam (Jakesz SABC 2005). Tam is metabolically activated by CYP2D6 to endoxifen and pts with impaired CYP2D6 metabolism have a higher risk of recurrence (Goetz JCO 2005). The HOXB13/IL17BR (H/I) index (Ma Cancer Cell 2004) is associated with Tam resistance while a 5-gene molecular grade index (MGI) distinguishes outcome in grade II tumors. H/I and MGI are complementary prognostic factors that define low (low MGI), intermediate (low H/I, high MGI) and high risk groups (high H/I and MGI) (Ma Clin Cancer Res 2008). We sought to validate these biomarkers in ABCSG Trial 8, and hypothesized that they may predict the benefit of sequencing TAM to A versus Tam alone.
 Methods: From 1996-2004, ABCSG 8 randomized 3,901 pts with ER positive breast cancer. The primary endpoint was event-free survival. Assuming 40% of pts had reduced CYP2D6 metabolism, we designed a nested case control study in which 110 Tam treated cases matched to 2 controls would yield 80% power to detect a relative risk of 2.3 for pts with reduced CYP2D6 metabolism. Additionally, we sought to perform an exploratory analysis of the benefit of switching from Tam to A in pts with impaired vs non-impaired CYP2D6 metabolism. DNA and RNA derived from paraffin embedded sections were used to genotype for CYP2D6 and to obtain H/I and MGI RT-PCR profiles.
 Findings: Thus far, CYP2D6 genotyping is complete in 345 pts corresponding to 53 cases/106 controls (Tam) and 62 cases/124 controls (Tam to A). DNA and RNA yields were excellent with a CYP2D6 genotype call rate of >95%. CYP2D6 allele frequencies were as expected: *3 (.01), *4, (.23), *6 (.01), *10 (.02), *17 (.00), and *41 (.10). 44% of pts had reduced CYP2D6 metabolism (at least one CYP2D6 null allele or intermediate metabolizer) and 8% were CYP2D6 poor metabolizers. For H/I and MGI, 340 pts were profiled with 64%, 20% and 16% in the low, intermediate, and high risk groups, respectively. Tissue sectioning/genotyping for additional cases/controls is ongoing.
 Discussion: This is the first retrospective analysis of a large Tam versus Tam to A trial to evaluate the association between 1) CYP2D6, HOXB13/IL17BR, and MGI with TAM treatment outcome and 2) the predictive benefit of these biomarkers in pts randomized to switch from Tam to AI versus Tam alone. Final results will be presented at the meeting. Citation Information: Cancer Res 2009;69(2 Suppl):Abstract nr 57.
Genome-Wide Associations and Functional Genomic Studies of Musculoskeletal Adverse Events in Women Receiving Aromatase Inhibitors
Purpose We performed a case-control genome-wide association study (GWAS) to identify single nucleotide polymorphisms (SNPs) associated with musculoskeletal adverse events (MS-AEs) in women treated with aromatase inhibitors (AIs) for early breast cancer. Patients and Methods A nested case-control design was used to select patients enrolled onto the MA.27 phase III trial comparing anastrozole with exemestane. Cases were matched to two controls and were defined as patients with grade 3 or 4 MS-AEs (according to the National Cancer Institute's Common Terminology Criteria for Adverse Events v3.0) or those who discontinued treatment for any grade of MS-AE within the first 2 years. Genotyping was performed with the Illumina Human610-Quad BeadChip. Results The GWAS included 293 cases and 585 controls. A total of 551,358 SNPs were analyzed, followed by imputation and fine mapping of a region of interest on chromosome 14. Four SNPs on chromosome 14 had the lowest P values (2.23E-06 to 6.67E-07). T-cell leukemia 1A (TCL1A) was the gene closest (926-7000 bp) to the four SNPs. Functional genomic studies revealed that one of these SNPs (rs11849538) created an estrogen response element and that TCL1A expression was estrogen dependent, was associated with the variant SNP genotypes in estradiol-treated lymphoblastoid cells transfected with estrogen receptor alpha and was directly related to interleukin 17 receptor A (IL17RA) expression. Conclusion This GWAS identified SNPs associated with MS-AEs in women treated with AIs and with a gene (TCL1A) which, in turn, was related to a cytokine (IL17). These findings provide a focus for further research to identify patients at risk for MS-AEs and to explore the mechanisms for these adverse events.
Strategic Plan for Lung Vascular Research: An NHLBI-ORDR Workshop Report
Abstract The Division of Lung Diseases of the National Heart, Lung, and Blood Institute, with the Office of Rare Diseases Research, held a workshop to identify priority areas and strategic goals to enhance and accelerate research that will result in improved understanding of the lung vasculature, translational research needs, and ultimately the care of patients with pulmonary vascular diseases. Multidisciplinary experts with diverse experience in laboratory, translational, and clinical studies identified seven priority areas and discussed limitations in our current knowledge, technologies, and approaches. The focus for future research efforts include the following: (1) better characterizing vascular genotype–phenotype relationships and incorporating systems biology approaches when appropriate; (2) advancing our understanding of pulmonary vascular metabolic regulatory signaling in health and disease; (3) expanding our knowledge of the biologic relationships between the lung circulation and circulating elements, systemic vascular function, and right heart function and disease; (4) improving translational research for identifying disease-modifying therapies for the pulmonary hypertensive diseases; (5) establishing an appropriate and effective platform for advancing translational findings into clinical studies testing; and (6) developing the specific technologies and tools that will be enabling for these goals, such as question-guided imaging techniques and lung vascular investigator training programs. Recommendations from this workshop will be used within the Lung Vascular Biology and Disease Extramural Research Program for planning and strategic implementation purposes.
Coprescription of Tamoxifen and Medications That Inhibit CYP2D6
Evidence has emerged that the clinical benefit of tamoxifen is related to the functional status of the hepatic metabolizing enzyme cytochrome P450 2D6 (CYP2D6). CYP2D6 is the key enzyme responsible for the generation of the potent tamoxifen metabolite, endoxifen. Multiple studies have examined the relationship of CYP2D6 status to breast cancer outcomes in tamoxifen-treated women; the majority of studies demonstrated that women with impaired CYP2D6 metabolism have lower endoxifen concentrations and a greater risk of breast cancer recurrence. As a result, practitioners must be aware that some of the most commonly prescribed medications coadministered with tamoxifen interfere with CYP2D6 function, thereby reducing endoxifen concentrations and potentially increasing the risk of breast cancer recurrence. After reviewing the published data regarding tamoxifen metabolism and the evidence relating CYP2D6 status to breast cancer outcomes in tamoxifen-treated patients, we are providing a guide for the use of medications that inhibit CYP2D6 in patients administered tamoxifen.
Variation in Anastrozole Metabolism and Pharmacodynamics in Women with Early Breast Cancer