Source of Patients
The source of the patients for this genome-wide case control study was MA.27, which was conducted as a multi-cooperative group effort under the auspices of the NCI Breast Cancer Intergroup of North America. The NCIC Clinical Trials Group (CTG) serves as the coordinating group, with participation by the NCI-sponsored North Central Cancer Treatment Group, Eastern Cooperative Oncology Group (ECOG), Southwest Oncology Group, and Cancer and Leukemia Group B (CALGB). MA.27 involved postmenopausal women with histologically confirmed and completely resected invasive breast cancer with surgical margins clear of invasive carcinoma and DCIS in the following TMN categories (AJCC Version 6): pT1, pT2, pT3; pNx, pN0, pN1, pN2, pN3 (only when the sole basis is presence of 10 or more involved axillary nodes); MO. The primary tumor must have been ER and/or progesterone receptor positive. Patients were stratified by lymph node status at diagnosis, prior adjuvant chemotherapy, and trastuzumab use and were randomized to 5 years of adjuvant therapy ... (Show More)

Source of Patients
The source of the patients for this genome-wide case control study was MA.27, which was conducted as a multi-cooperative group effort under the auspices of the NCI Breast Cancer Intergroup of North America. The NCIC Clinical Trials Group (CTG) serves as the coordinating group, with participation by the NCI-sponsored North Central Cancer Treatment Group, Eastern Cooperative Oncology Group (ECOG), Southwest Oncology Group, and Cancer and Leukemia Group B (CALGB). MA.27 involved postmenopausal women with histologically confirmed and completely resected invasive breast cancer with surgical margins clear of invasive carcinoma and DCIS in the following TMN categories (AJCC Version 6): pT1, pT2, pT3; pNx, pN0, pN1, pN2, pN3 (only when the sole basis is presence of 10 or more involved axillary nodes); MO. The primary tumor must have been ER and/or progesterone receptor positive. Patients were stratified by lymph node status at diagnosis, prior adjuvant chemotherapy, and trastuzumab use and were randomized to 5 years of adjuvant therapy with anastrozole or exemestane. The trial was activated on May 26, 2003, and reached its accrual objectives on July 31, 2008, after the randomization of 6827 North American patients, with the majority (79%) providing DNA and consent for genetic testing. Non-North American patients were also entered by the International Breast Cancer Study Group but they did not contribute DNA. From 2003 to December 21, 2004, patients also underwent a second randomization to celecoxib 400 mg twice daily or placebo but, after the entry of 1,622 patients, this treatment was discontinued because of reports of increased cardiovascular risk associated with celecoxib.

Case Definition for Musculoskeletal Adverse Events
Cases were defined to have had at least one of the following six MSAE): 1. pain - joint, 2. pain - muscle, 3. pain - bone, 4. arthritis, 5. joint - function, 6. musculoskeletal - other. Cases were required to have at least grade 3 toxicity, according to the NCI Common Terminology Criteria for Adverse Events v3.0, or go off treatment for any grade of the six MSAEs; and were required to have had an MSAE within the first two years (i.e., a MSAE occurring after 2 years was not be considered a case). Subjects were excluded from the case group who met the case definition while on celecoxib or in the three months after stopping celecoxib but subjects were allowed to be cases if they developed their first MSAE more than three months after stopping celecoxib.

Study Design
The design was a nested matched case-control study with two controls selected for each case. Random matching from the potential control group (i.e., risk-set sampling such that controls had at least the same length of followup as the timing of the case event) was done by matching the case on the following factors: 1. treatment arm (exemestane vs. anastrozole, blinded), 2. prior adjuvant chemotherapy (yes vs. no), 3. age at start of treatment (+/- 5 years), 4. Celecoxib allocation (yes vs. no). Attempts were made to exactly match each case with two controls according to the above criteria. If, however, controls could not be exactly matched, "close" matching was used to define a distance between each case and all potential controls in an optimal matching algorithm. Race is a potential confounding factor, and 94% of the MA.27 population is Caucasian, 3% are Black, 1% Asian, 1% American Indian, and 1% are of unknown race. Because there are so few non-Caucasian participants, was restricted to only Caucasians. Following identification of cases and controls, the list was sent to the CALGB Pathology Coordinating Office (PCO) at Ohio State University Medical Center for plating of the samples. Other patient factors that might influence the association of SNP genotypes with a MSAE, either as confounding factors or as effect modifiers, were included as investigational covariates. These include 1. body mass index (BMI), 2. bisphosphonate use, 3. fractures in past 10 years, 4. baseline ECOG performance status, 5. prior hormone replacement therapy, 6. prior adjuvant radiotherapy, and 7. prior taxane treatment. An Honest Broker Process was utilized to ensure patient confidentiality.

Genotype Assays and Quality Control
A GWAS was performed with the Illumina Human610-Quad at the RIKEN Center for Genomic Medicine. The DNA for all of the samples were adjusted to a concentration of 50 to 100 ng/µL in a final volume of 50 µL. The 96-well plates were designed to have a unique distribution pattern with x duplicates per plate to evaluate genotype concordance. Samples with <98% call rate were removed from analyses. SNPs with ambiguous call will be excluded using the following criteria: 1. low signal intensity of each allele, 2. low cluster separation score, 3. SNP call rate <0.98, 4. SNPs with Hardy Weinberg Equilibrium p-value <0.001 were excluded. X-chromosome SNPs were used to validate female sex and quality control measures were examined across the different submitting centers. For all pairs of subjects, the percentage of matching genotypes were computed to screen for unexpected duplicate samples. All SNPs found to be associated with MSAE were subjected to visual inspection of the cluster plots to verify that the genotype clustering was adequate.