Risk Prediction for Epithelial Ovarian Cancer in 11 United States-Based Case-Control Studies: Incorporation of Epidemiologic Risk Factors and 17 Confirmed Genetic Loci

Previously developed models for predicting absolute risk of invasive epithelial ovarian cancer have included a limited number of risk factors and have had low discriminatory power (area under the receiver operating characteristic curve (AUC) < 0.60). Because of this, we developed and internally validated a relative risk prediction model that incorporates 17 established epidemiologic risk factors and 17 genome-wide significant single nucleotide polymorphisms (SNPs) using data from 11 case-control studies in the United States (5,793 cases; 9,512 controls) from the Ovarian Cancer Association Consortium (data accrued from 1992 to 2010). We developed a hierarchical logistic regression model for predicting case-control status that included imputation of missing data. We randomly divided the data into an 80% training sample and used the remaining 20% for model evaluation. The AUC for the full model was 0.664. A reduced model without SNPs performed similarly (AUC = 0.649). Both models performed better than a baseline model that included age and study site only (AUC = 0.563). The best predictive power was obtained in the full model among women younger than 50 years of age (AUC = 0.714); however, the addition of SNPs increased the AUC the most for women older than 50 years of age (AUC = 0.638 vs. 0.616). Adapting this improved model to estimate absolute risk and evaluating it in prospective data sets is warranted.

REVEL: An Ensemble Method for Predicting the Pathogenicity of Rare Missense Variants

The vast majority of coding variants are rare, and assessment of the contribution of rare variants to complex traits is hampered by low statistical power and limited functional data. Improved methods for predicting the pathogenicity of rare coding variants are needed to facilitate the discovery of disease variants from exome sequencing studies. We developed REVEL (rare exome variant ensemble learner), an ensemble method for predicting the pathogenicity of missense variants on the basis of individual tools: MutPred, FATHMM, VEST, PolyPhen, SIFT, PROVEAN, MutationAssessor, MutationTaster, LRT, GERP, SiPhy, phyloP, and phastCons. REVEL was trained with recently discovered pathogenic and rare neutral missense variants, excluding those previously used to train its constituent tools. When applied to two independent test sets, REVEL had the best overall performance (p < 10^-12) as compared to any individual tool and seven ensemble methods: MetaSVM, MetaLR, KGGSeq, Condel, CADD, DANN, and Eigen. Importantly, REVEL also had the best performance for distinguishing pathogenic from rare neutral variants with allele frequencies <0.5%. The area under the receiver operating characteristic curve (AUC) for REVEL was 0.046-0.182 higher in an independent test set of 935 recent SwissVar disease variants and 123,935 putatively neutral exome sequencing variants and 0.027-0.143 higher in an independent test set of 1,953 pathogenic and 2,406 benign variants recently reported in ClinVar than the AUCs for other ensemble methods. We provide pre-computed REVEL scores for all possible human missense variants to facilitate the identification of pathogenic variants in the sea of rare variants discovered as sequencing studies expand in scale.

Association of vitamin D levels and risk of ovarian cancer: a Mendelian randomization study

BACKGROUND

In vitro and observational epidemiological studies suggest that vitamin D may play a role in cancer prevention. However, the relationship between vitamin D and ovarian cancer is uncertain, with observational studies generating conflicting findings. A potential limitation of observational studies is inadequate control of confounding. To overcome this problem, we used Mendelian randomization (MR) to evaluate the association between single nucleotide polymorphisms (SNPs) associated with circulating 25-hydroxyvitamin D [25(OH)D] concentration and risk of ovarian cancer.

METHODS

We employed SNPs with well-established associations with 25(OH)D concentration as instrumental variables for MR: rs7944926 (DHCR7), rs12794714 (CYP2R1) and rs2282679 (GC). We included 31 719 women of European ancestry (10 065 cases, 21 654 controls) from the Ovarian Cancer Association Consortium, who were genotyped using customized Illumina Infinium iSelect (iCOGS) arrays. A two-sample (summary data) MR approach was used and analyses were performed separately for all ovarian cancer (10 065 cases) and for high-grade serous ovarian cancer (4121 cases).

RESULTS

The odds ratio for epithelial ovarian cancer risk (10 065 cases) estimated by combining the individual SNP associations using inverse variance weighting was 1.27 (95% confidence interval: 1.06 to 1.51) per 20 nmol/L decrease in 25(OH)D concentration. The estimated odds ratio for high-grade serous epithelial ovarian cancer (4121 cases) was 1.54 (1.19, 2.01).

CONCLUSIONS

Genetically lowered 25-hydroxyvitamin D concentrations were associated with higher ovarian cancer susceptibility in Europeans. These findings suggest that increasing plasma vitamin D levels may reduce risk of ovarian cancer.

Functional mechanisms underlying pleiotropic risk alleles at the 19p13.1 breast-ovarian cancer susceptibility locus

A locus at 19p13 is associated with breast cancer (BC) and ovarian cancer (OC) risk. Here we analyse 438 SNPs in this region in 46,451 BC and 15,438 OC cases, 15,252 BRCA1 mutation carriers and 73,444 controls and identify 13 candidate causal SNPs associated with serous OC (P=9.2 × 10(-20)), ER-negative BC (P=1.1 × 10(-13)), BRCA1-associated BC (P=7.7 × 10(-16)) and triple negative BC (P-diff=2 × 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P=2 × 10(-3)) and ABHD8 (P<2 × 10(-3)). Chromosome conformation capture identifies interactions between four candidate SNPs and ABHD8, and luciferase assays indicate six risk alleles increased transactivation of the ADHD8 promoter. Targeted deletion of a region containing risk SNP rs56069439 in a putative enhancer induces ANKLE1 downregulation; and mRNA stability assays indicate functional effects for an ANKLE1 3'-UTR SNP. Altogether, these data suggest that multiple SNPs at 19p13 regulate ABHD8 and perhaps ANKLE1 expression, and indicate common mechanisms underlying breast and ovarian cancer risk.

Breast Cancer Risk and Mammographic Density Assessed with Semiautomated and Fully Automated Methods and BI-RADS

Purpose To compare three metrics of breast density on full-field digital mammographic (FFDM) images as predictors of future breast cancer risk. Materials and Methods This institutional review board-approved study included 125 women with invasive breast cancer and 274 age- and race-matched control subjects who underwent screening FFDM during 2004-2013 and provided informed consent. The percentage of density and dense area were assessed semiautomatically with software (Cumulus 4.0; University of Toronto, Toronto, Canada), and volumetric percentage of density and dense volume were assessed automatically with software (Volpara; Volpara Solutions, Wellington, New Zealand). Clinical Breast Imaging Reporting and Data System (BI-RADS) classifications of breast density were extracted from mammography reports. Odds ratios and 95% confidence intervals (CIs) were estimated by using conditional logistic regression stratified according to age and race and adjusted for body mass index, parity, and menopausal status, and the area under the receiver operating characteristic curve (AUC) was computed. Results The adjusted odds ratios and 95% CIs for each standard deviation increment of the percentage of density, dense area, volumetric percentage of density, and dense volume were 1.61 (95% CI: 1.19, 2.19), 1.49 (95% CI: 1.15, 1.92), 1.54 (95% CI: 1.12, 2.10), and 1.41 (95% CI: 1.11, 1.80), respectively. Odds ratios for women with extremely dense breasts compared with those with scattered areas of fibroglandular density were 2.06 (95% CI: 0.85, 4.97) and 2.05 (95% CI: 0.90, 4.64) for BI-RADS and Volpara density classifications, respectively. Clinical BI-RADS was more accurate (AUC, 0.68; 95% CI: 0.63, 0.74) than Volpara (AUC, 0.64; 95% CI: 0.58, 0.70) and continuous measures of percentage of density (AUC, 0.66; 95% CI: 0.60, 0.72), dense area (AUC, 0.66; 95% CI: 0.60, 0.72), volumetric percentage of density (AUC, 0.64; 95% CI: 0.58, 0.70), and density volume (AUC, 0.65; 95% CI: 0.59, 0.71), although the AUC differences were not statistically significant. Conclusion Mammographic density on FFDM images was positively associated with breast cancer risk by using the computer assisted methods and BI-RADS. BI-RADS classification was as accurate as computer-assisted methods for discrimination of patients from control subjects. ^© RSNA, 2016.

Genome-Wide Meta-Analyses of Breast, Ovarian, and Prostate Cancer Association Studies Identify Multiple New Susceptibility Loci Shared by at Least Two Cancer Types

Breast, ovarian, and prostate cancers are hormone-related and may have a shared genetic basis, but this has not been investigated systematically by genome-wide association (GWA) studies. Meta-analyses combining the largest GWA meta-analysis data sets for these cancers totaling 112,349 cases and 116,421 controls of European ancestry, all together and in pairs, identified at P < 10(-8) seven new cross-cancer loci: three associated with susceptibility to all three cancers (rs17041869/2q13/BCL2L11; rs7937840/11q12/INCENP; rs1469713/19p13/GATAD2A), two breast and ovarian cancer risk loci (rs200182588/9q31/SMC2; rs8037137/15q26/RCCD1), and two breast and prostate cancer risk loci (rs5013329/1p34/NSUN4; rs9375701/6q23/L3MBTL3). Index variants in five additional regions previously associated with only one cancer also showed clear association with a second cancer type. Cell-type-specific expression quantitative trait locus and enhancer-gene interaction annotations suggested target genes with potential cross-cancer roles at the new loci. Pathway analysis revealed significant enrichment of death receptor signaling genes near loci with P < 10(-5) in the three-cancer meta-analysis.

SIGNIFICANCE

We demonstrate that combining large-scale GWA meta-analysis findings across cancer types can identify completely new risk loci common to breast, ovarian, and prostate cancers. We show that the identification of such cross-cancer risk loci has the potential to shed new light on the shared biology underlying these hormone-related cancers. Cancer Discov; 6(9); 1052-67. ©2016 AACR.This article is highlighted in the In This Issue feature, p. 932.

Cross-Cancer Genome-Wide Analysis of Lung, Ovary, Breast, Prostate, and Colorectal Cancer Reveals Novel Pleiotropic Associations

Identifying genetic variants with pleiotropic associations can uncover common pathways influencing multiple cancers. We took a two-stage approach to conduct genome-wide association studies for lung, ovary, breast, prostate, and colorectal cancer from the GAME-ON/GECCO Network (61,851 cases, 61,820 controls) to identify pleiotropic loci. Findings were replicated in independent association studies (55,789 cases, 330,490 controls). We identified a novel pleiotropic association at 1q22 involving breast and lung squamous cell carcinoma, with eQTL analysis showing an association with ADAM15/THBS3 gene expression in lung. We also identified a known breast cancer locus CASP8/ALS2CR12 associated with prostate cancer, a known cancer locus at CDKN2B-AS1 with different variants associated with lung adenocarcinoma and prostate cancer, and confirmed the associations of a breast BRCA2 locus with lung and serous ovarian cancer. This is the largest study to date examining pleiotropy across multiple cancer-associated loci, identifying common mechanisms of cancer development and progression. Cancer Res; 76(17); 5103-14. ©2016 AACR.

Post hoc Analysis for Detecting Individual Rare Variant Risk Associations Using Probit Regression Bayesian Variable Selection Methods in Case-Control Sequencing Studies

Rare variants (RVs) have been shown to be significant contributors to complex disease risk. By definition, these variants have very low minor allele frequencies and traditional single-marker methods for statistical analysis are underpowered for typical sequencing study sample sizes. Multimarker burden-type approaches attempt to identify aggregation of RVs across case-control status by analyzing relatively small partitions of the genome, such as genes. However, it is generally the case that the aggregative measure would be a mixture of causal and neutral variants, and these omnibus tests do not directly provide any indication of which RVs may be driving a given association. Recently, Bayesian variable selection approaches have been proposed to identify RV associations from a large set of RVs under consideration. Although these approaches have been shown to be powerful at detecting associations at the RV level, there are often computational limitations on the total quantity of RVs under consideration and compromises are necessary for large-scale application. Here, we propose a computationally efficient alternative formulation of this method using a probit regression approach specifically capable of simultaneously analyzing hundreds to thousands of RVs. We evaluate our approach to detect causal variation on simulated data and examine sensitivity and specificity in instances of high RV dimensionality as well as apply it to pathway-level RV analysis results from a prostate cancer (PC) risk case-control sequencing study. Finally, we discuss potential extensions and future directions of this work.

Age- and Tumor Subtype-Specific Breast Cancer Risk Estimates for CHEK2*1100delC Carriers

PURPOSE

CHEK2*1100delC is a well-established breast cancer risk variant that is most prevalent in European populations; however, there are limited data on risk of breast cancer by age and tumor subtype, which limits its usefulness in breast cancer risk prediction. We aimed to generate tumor subtype- and age-specific risk estimates by using data from the Breast Cancer Association Consortium, including 44,777 patients with breast cancer and 42,997 controls from 33 studies genotyped for CHEK2*1100delC.

PATIENTS AND METHODS

CHEK2*1100delC genotyping was mostly done by a custom Taqman assay. Breast cancer odds ratios (ORs) for CHEK2*1100delC carriers versus noncarriers were estimated by using logistic regression and adjusted for study (categorical) and age. Main analyses included patients with invasive breast cancer from population- and hospital-based studies.

RESULTS

Proportions of heterozygous CHEK2*1100delC carriers in controls, in patients with breast cancer from population- and hospital-based studies, and in patients with breast cancer from familial- and clinical genetics center-based studies were 0.5%, 1.3%, and 3.0%, respectively. The estimated OR for invasive breast cancer was 2.26 (95%CI, 1.90 to 2.69; P = 2.3 × 10(-20)). The OR was higher for estrogen receptor (ER)-positive disease (2.55 [95%CI, 2.10 to 3.10; P = 4.9 × 10(-21)]) than it was for ER-negative disease (1.32 [95%CI, 0.93 to 1.88; P = .12]; P interaction = 9.9 × 10(-4)). The OR significantly declined with attained age for breast cancer overall (P = .001) and for ER-positive tumors (P = .001). Estimated cumulative risks for development of ER-positive and ER-negative tumors by age 80 in CHEK2*1100delC carriers were 20% and 3%, respectively, compared with 9% and 2%, respectively, in the general population of the United Kingdom.

CONCLUSION

These CHEK2*1100delC breast cancer risk estimates provide a basis for incorporating CHEK2*1100delC into breast cancer risk prediction models and into guidelines for intensified screening and follow-up.

Genetically Predicted Body Mass Index and Breast Cancer Risk: Mendelian Randomization Analyses of Data from 145,000 Women of European Descent

BACKGROUND

Observational epidemiological studies have shown that high body mass index (BMI) is associated with a reduced risk of breast cancer in premenopausal women but an increased risk in postmenopausal women. It is unclear whether this association is mediated through shared genetic or environmental factors.

METHODS

We applied Mendelian randomization to evaluate the association between BMI and risk of breast cancer occurrence using data from two large breast cancer consortia. We created a weighted BMI genetic score comprising 84 BMI-associated genetic variants to predicted BMI. We evaluated genetically predicted BMI in association with breast cancer risk using individual-level data from the Breast Cancer Association Consortium (BCAC) (cases  =  46,325, controls  =  42,482). We further evaluated the association between genetically predicted BMI and breast cancer risk using summary statistics from 16,003 cases and 41,335 controls from the Discovery, Biology, and Risk of Inherited Variants in Breast Cancer (DRIVE) Project. Because most studies measured BMI after cancer diagnosis, we could not conduct a parallel analysis to adequately evaluate the association of measured BMI with breast cancer risk prospectively.

RESULTS

In the BCAC data, genetically predicted BMI was found to be inversely associated with breast cancer risk (odds ratio [OR]  =  0.65 per 5 kg/m2 increase, 95% confidence interval [CI]: 0.56-0.75, p = 3.32 × 10-10). The associations were similar for both premenopausal (OR   =   0.44, 95% CI:0.31-0.62, p  =  9.91 × 10-8) and postmenopausal breast cancer (OR  =  0.57, 95% CI: 0.46-0.71, p  =  1.88 × 10-8). This association was replicated in the data from the DRIVE consortium (OR  =  0.72, 95% CI: 0.60-0.84, p   =   1.64 × 10-7). Single marker analyses identified 17 of the 84 BMI-associated single nucleotide polymorphisms (SNPs) in association with breast cancer risk at p < 0.05; for 16 of them, the allele associated with elevated BMI was associated with reduced breast cancer risk.

CONCLUSIONS

BMI predicted by genome-wide association studies (GWAS)-identified variants is inversely associated with the risk of both pre- and postmenopausal breast cancer. The reduced risk of postmenopausal breast cancer associated with genetically predicted BMI observed in this study differs from the positive association reported from studies using measured adult BMI. Understanding the reasons for this discrepancy may reveal insights into the complex relationship of genetic determinants of body weight in the etiology of breast cancer.

Assessing the genetic architecture of epithelial ovarian cancer histological subtypes

Epithelial ovarian cancer (EOC) is one of the deadliest common cancers. The five most common types of disease are high-grade and low-grade serous, endometrioid, mucinous and clear cell carcinoma. Each of these subtypes present distinct molecular pathogeneses and sensitivities to treatments. Recent studies show that certain genetic variants confer susceptibility to all subtypes while other variants are subtype-specific. Here, we perform an extensive analysis of the genetic architecture of EOC subtypes. To this end, we used data of 10,014 invasive EOC patients and 21,233 controls from the Ovarian Cancer Association Consortium genotyped in the iCOGS array (211,155 SNPs). We estimate the array heritability (attributable to variants tagged on arrays) of each subtype and their genetic correlations. We also look for genetic overlaps with factors such as obesity, smoking behaviors, diabetes, age at menarche and height. We estimated the array heritabilities of high-grade serous disease ([Formula: see text] = 8.8 ± 1.1 %), endometrioid ([Formula: see text] = 3.2 ± 1.6 %), clear cell ([Formula: see text] = 6.7 ± 3.3 %) and all EOC ([Formula: see text] = 5.6 ± 0.6 %). Known associated loci contributed approximately 40 % of the total array heritability for each subtype. The contribution of each chromosome to the total heritability was not proportional to chromosome size. Through bivariate and cross-trait LD score regression, we found evidence of shared genetic backgrounds between the three high-grade subtypes: serous, endometrioid and undifferentiated. Finally, we found significant genetic correlations of all EOC with diabetes and obesity using a polygenic prediction approach.

Parity and Oral Contraceptive Use in Relation to Ovarian Cancer Risk in Older Women

BACKGROUND

Several studies have suggested that the ovarian cancer risk reductions associated with parity and oral contraceptive use are weaker in postmenopausal than premenopausal women, yet little is known about the persistence of these reductions as women age. This question gains importance with the increasing numbers of older women in the population.

METHODS

We addressed the question using data from three large U.S. cohort studies involving 310,290 white women aged 50+ years at recruitment, of whom 1,815 developed subsequent incident invasive epithelial ovarian cancer. We used Cox regression, stratified by cohort, to examine age-related trends in the HRs per full-term pregnancy and per year of oral contraceptive use.

RESULTS

The parity-associated risk reductions waned with age (Ptrend < 0.001 in HR with increasing age), particularly among women aged 75 years or more, for whom we observed no association with parity. However, we observed no such attenuation in the oral contraceptive-associated risk reductions (P = 0.79 for trend in HR with increasing age).

CONCLUSION

These findings suggest that prior oral contraceptive use is important for ovarian cancer risk assessment among women of all ages, while the benefits of parity wane as women age.

IMPACT

This information, if duplicated in other studies, will be useful to preventive counseling and risk prediction, particularly for women at increased ovarian cancer risk due to a personal history of breast cancer or a family history of ovarian cancer. Cancer Epidemiol Biomarkers Prev; 25(7); 1059-63. ©2016 AACR.

Identification of independent association signals and putative functional variants for breast cancer risk through fine-scale mapping of the 12p11 locus

BACKGROUND

Multiple recent genome-wide association studies (GWAS) have identified a single nucleotide polymorphism (SNP), rs10771399, at 12p11 that is associated with breast cancer risk.

METHOD

We performed a fine-scale mapping study of a 700 kb region including 441 genotyped and more than 1300 imputed genetic variants in 48,155 cases and 43,612 controls of European descent, 6269 cases and 6624 controls of East Asian descent and 1116 cases and 932 controls of African descent in the Breast Cancer Association Consortium (BCAC; http://bcac.ccge.medschl.cam.ac.uk/ ), and in 15,252 BRCA1 mutation carriers in the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA). Stepwise regression analyses were performed to identify independent association signals. Data from the Encyclopedia of DNA Elements project (ENCODE) and the Cancer Genome Atlas (TCGA) were used for functional annotation.

RESULTS

Analysis of data from European descendants found evidence for four independent association signals at 12p11, represented by rs7297051 (odds ratio (OR) = 1.09, 95 % confidence interval (CI) = 1.06-1.12; P = 3 × 10(-9)), rs805510 (OR = 1.08, 95 % CI = 1.04-1.12, P = 2 × 10(-5)), and rs1871152 (OR = 1.04, 95 % CI = 1.02-1.06; P = 2 × 10(-4)) identified in the general populations, and rs113824616 (P = 7 × 10(-5)) identified in the meta-analysis of BCAC ER-negative cases and BRCA1 mutation carriers. SNPs rs7297051, rs805510 and rs113824616 were also associated with breast cancer risk at P < 0.05 in East Asians, but none of the associations were statistically significant in African descendants. Multiple candidate functional variants are located in putative enhancer sequences. Chromatin interaction data suggested that PTHLH was the likely target gene of these enhancers. Of the six variants with the strongest evidence of potential functionality, rs11049453 was statistically significantly associated with the expression of PTHLH and its nearby gene CCDC91 at P < 0.05.

CONCLUSION

This study identified four independent association signals at 12p11 and revealed potentially functional variants, providing additional insights into the underlying biological mechanism(s) for the association observed between variants at 12p11 and breast cancer risk.

Adult body mass index and risk of ovarian cancer by subtype: a Mendelian randomization study

BACKGROUND

Observational studies have reported a positive association between body mass index (BMI) and ovarian cancer risk. However, questions remain as to whether this represents a causal effect, or holds for all histological subtypes. The lack of association observed for serous cancers may, for instance, be due to disease-associated weight loss. Mendelian randomization (MR) uses genetic markers as proxies for risk factors to overcome limitations of observational studies. We used MR to elucidate the relationship between BMI and ovarian cancer, hypothesizing that genetically predicted BMI would be associated with increased risk of non-high grade serous ovarian cancers (non-HGSC) but not HGSC.

METHODS

We pooled data from 39 studies (14 047 cases, 23 003 controls) in the Ovarian Cancer Association Consortium. We constructed a weighted genetic risk score (GRS, partial F-statistic = 172), summing alleles at 87 single nucleotide polymorphisms previously associated with BMI, weighting by their published strength of association with BMI. Applying two-stage predictor-substitution MR, we used logistic regression to estimate study-specific odds ratios (OR) and 95% confidence intervals (CI) for the association between genetically predicted BMI and risk, and pooled these using random-effects meta-analysis.

RESULTS

Higher genetically predicted BMI was associated with increased risk of non-HGSC (pooled OR = 1.29, 95% CI 1.03-1.61 per 5 units BMI) but not HGSC (pooled OR = 1.06, 95% CI 0.88-1.27). Secondary analyses stratified by behaviour/subtype suggested that, consistent with observational data, the association was strongest for low-grade/borderline serous cancers (OR = 1.93, 95% CI 1.33-2.81).

CONCLUSIONS

Our data suggest that higher BMI increases risk of non-HGSC, but not the more common and aggressive HGSC subtype, confirming the observational evidence.

Case-control study of mammographic density and breast cancer risk using processed digital mammograms

Background

Full-field digital mammography (FFDM) has largely replaced film-screen mammography in the US. Breast density assessed from film mammograms is strongly associated with breast cancer risk, but data are limited for processed FFDM images used for clinical care.

Methods

We conducted a case-control study nested among non-Hispanic white female participants of the Research Program in Genes, Environment and Health of Kaiser Permanente Northern California who were aged 40 to 74 years and had screening mammograms acquired on Hologic FFDM machines. Cases (n = 297) were women with a first invasive breast cancer diagnosed after a screening FFDM. For each case, up to five controls (n = 1149) were selected, matched on age and year of FFDM and image batch number, and who were still under follow-up and without a history of breast cancer at the age of diagnosis of the matched case. Percent density (PD) and dense area (DA) were assessed by a radiological technologist using Cumulus. Conditional logistic regression was used to estimate odds ratios (ORs) for breast cancer associated with PD and DA, modeled continuously in standard deviation (SD) increments and categorically in quintiles, after adjusting for body mass index, parity, first-degree family history of breast cancer, breast area, and menopausal hormone use.

Results

Median intra-reader reproducibility was high with a Pearson’s r of 0.956 (range 0.902 to 0.983) for replicate PD measurements across 23 image batches. The overall mean was 20.02 (SD, 14.61) for PD and 27.63 cm² (18.22 cm²) for DA. The adjusted ORs for breast cancer associated with each SD increment were 1.70 (95 % confidence interval, 1.41–2.04) for PD, and 1.54 (1.34–1.77) for DA. The adjusted ORs for each quintile were: 1.00 (ref.), 1.49 (0.91–2.45), 2.57 (1.54–4.30), 3.22 (1.91–5.43), 4.88 (2.78–8.55) for PD, and 1.00 (ref.), 1.43 (0.85–2.40), 2.53 (1.53–4.19), 2.85 (1.73–4.69), 3.48 (2.14–5.65) for DA.

Conclusions

PD and DA measured using Cumulus on processed FFDM images are positively associated with breast cancer risk, with similar magnitudes of association as previously reported for film-screen mammograms. Processed digital mammograms acquired for routine clinical care in a general practice setting are suitable for breast density and cancer research.

Association of genetic susceptibility variants for type 2 diabetes with breast cancer risk in women of European ancestry

PURPOSE

Type 2 diabetes (T2D) has been reported to be associated with an elevated risk of breast cancer. It is unclear, however, whether this association is due to shared genetic factors.

METHODS

We constructed a genetic risk score (GRS) using risk variants from 33 known independent T2D susceptibility loci and evaluated its relation to breast cancer risk using the data from two consortia, including 62,328 breast cancer patients and 83,817 controls of European ancestry. Unconditional logistic regression models were used to derive adjusted odds ratios (ORs) and 95 % confidence intervals (CIs) to measure the association of breast cancer risk with T2D GRS or T2D-associated genetic risk variants. Meta-analyses were conducted to obtain summary ORs across all studies.

RESULTS

The T2D GRS was not found to be associated with breast cancer risk, overall, by menopausal status, or for estrogen receptor positive or negative breast cancer. Three T2D associated risk variants were individually associated with breast cancer risk after adjustment for multiple comparisons using the Bonferroni method (at p < 0.001), rs9939609 (FTO) (OR 0.94, 95 % CI = 0.92-0.95, p = 4.13E-13), rs7903146 (TCF7L2) (OR 1.04, 95 % CI = 1.02-1.06, p = 1.26E-05), and rs8042680 (PRC1) (OR 0.97, 95 % CI = 0.95-0.99, p = 8.05E-04).

CONCLUSIONS

We have shown that several genetic risk variants were associated with the risk of both T2D and breast cancer. However, overall genetic susceptibility to T2D may not be related to breast cancer risk.

No clinical utility of KRAS variant rs61764370 for ovarian or breast cancer

OBJECTIVE

Clinical genetic testing is commercially available for rs61764370, an inherited variant residing in a KRAS 3' UTR microRNA binding site, based on suggested associations with increased ovarian and breast cancer risk as well as with survival time. However, prior studies, emphasizing particular subgroups, were relatively small. Therefore, we comprehensively evaluated ovarian and breast cancer risks as well as clinical outcome associated with rs61764370.

METHODS

Centralized genotyping and analysis were performed for 140,012 women enrolled in the Ovarian Cancer Association Consortium (15,357 ovarian cancer patients; 30,816 controls), the Breast Cancer Association Consortium (33,530 breast cancer patients; 37,640 controls), and the Consortium of Modifiers of BRCA1 and BRCA2 (14,765 BRCA1 and 7904 BRCA2 mutation carriers).

RESULTS

We found no association with risk of ovarian cancer (OR=0.99, 95% CI 0.94-1.04, p=0.74) or breast cancer (OR=0.98, 95% CI 0.94-1.01, p=0.19) and results were consistent among mutation carriers (BRCA1, ovarian cancer HR=1.09, 95% CI 0.97-1.23, p=0.14, breast cancer HR=1.04, 95% CI 0.97-1.12, p=0.27; BRCA2, ovarian cancer HR=0.89, 95% CI 0.71-1.13, p=0.34, breast cancer HR=1.06, 95% CI 0.94-1.19, p=0.35). Null results were also obtained for associations with overall survival following ovarian cancer (HR=0.94, 95% CI 0.83-1.07, p=0.38), breast cancer (HR=0.96, 95% CI 0.87-1.06, p=0.38), and all other previously-reported associations.

CONCLUSIONS

rs61764370 is not associated with risk of ovarian or breast cancer nor with clinical outcome for patients with these cancers. Therefore, genotyping this variant has no clinical utility related to the prediction or management of these cancers.

The LEGACY Girls Study: Growth and Development in the Context of Breast Cancer Family History

BACKGROUND

Although the timing of pubertal milestones has been associated with breast cancer risk, few studies of girls' development include girls at increased breast cancer risk due to their family history.

METHODS

The Lessons in Epidemiology and Genetics of Adult Cancer from Youth (LEGACY) Girls Study was initiated in 2011 in the USA and Canada to assess the relation between early life exposures and intermediate markers of breast cancer risk (e.g., pubertal development, breast tissue characteristics) and to investigate psychosocial well being and health behaviors in the context of family history. We describe the methods used to establish and follow a cohort of 1,040 girls ages 6-13 years at baseline, half with a breast cancer family history, and the collection of questionnaire data (family history, early life exposures, growth and development, psychosocial and behavioral), anthropometry, biospecimens, and breast tissue characteristics using optical spectroscopy.

RESULTS

During this initial 5-year phase of the study, follow-up visits are conducted every 6 months for repeated data and biospecimen collection. Participation in baseline components was high (98% for urine, 97.5% for blood or saliva, and 98% for anthropometry). At enrollment, 77% of girls were premenarcheal and 49% were at breast Tanner stage T1.

CONCLUSIONS

This study design allows thorough examination of events affecting girls' growth and development and how they differ across the spectrum of breast cancer risk. A better understanding of early life breast cancer risk factors will be essential to enhance prevention across the lifespan for those with and without a family history of the disease.

Assessment of Multifactor Gene-Environment Interactions and Ovarian Cancer Risk: Candidate Genes, Obesity, and Hormone-Related Risk Factors

BACKGROUND

Many epithelial ovarian cancer (EOC) risk factors relate to hormone exposure and elevated estrogen levels are associated with obesity in postmenopausal women. Therefore, we hypothesized that gene-environment interactions related to hormone-related risk factors could differ between obese and non-obese women.

METHODS

We considered interactions between 11,441 SNPs within 80 candidate genes related to hormone biosynthesis and metabolism and insulin-like growth factors with six hormone-related factors (oral contraceptive use, parity, endometriosis, tubal ligation, hormone replacement therapy, and estrogen use) and assessed whether these interactions differed between obese and non-obese women. Interactions were assessed using logistic regression models and data from 14 case-control studies (6,247 cases; 10,379 controls). Histotype-specific analyses were also completed.

RESULTS

SNPs in the following candidate genes showed notable interaction: IGF1R (rs41497346, estrogen plus progesterone hormone therapy, histology = all, P = 4.9 × 10(-6)) and ESR1 (rs12661437, endometriosis, histology = all, P = 1.5 × 10(-5)). The most notable obesity-gene-hormone risk factor interaction was within INSR (rs113759408, parity, histology = endometrioid, P = 8.8 × 10(-6)).

CONCLUSIONS

We have demonstrated the feasibility of assessing multifactor interactions in large genetic epidemiology studies. Follow-up studies are necessary to assess the robustness of our findings for ESR1, CYP11A1, IGF1R, CYP11B1, INSR, and IGFBP2 Future work is needed to develop powerful statistical methods able to detect these complex interactions.

IMPACT

Assessment of multifactor interaction is feasible, and, here, suggests that the relationship between genetic variants within candidate genes and hormone-related risk factors may vary EOC susceptibility. Cancer Epidemiol Biomarkers Prev; 25(5); 780-90. ©2016 AACR.

The association between socioeconomic status and tumour stage at diagnosis of ovarian cancer: A pooled analysis of 18 case-control studies

PURPOSE

Socioeconomic status (SES) is a known predictor of survival for several cancers and it has been suggested that SES differences affecting tumour stage at diagnosis may be the most important explanatory factor for this. However, only a limited number of studies have investigated SES differences in tumour stage at diagnosis of ovarian cancer. In a pooled analysis, we investigated whether SES as represented by level of education is predictive for advanced tumour stage at diagnosis of ovarian cancer, overall and by histotype. The effect of cigarette smoking and body mass index (BMI) on the association was also evaluated.

METHODS

From 18 case-control studies, we obtained information on 10,601 women diagnosed with epithelial ovarian cancer. Study specific odds ratios (ORs) with corresponding 95% confidence intervals (CI) were obtained from logistic regression models and combined into a pooled odds ratio (pOR) using a random effects model.

RESULTS

Overall, women who completed ≤high school had an increased risk of advanced tumour stage at diagnosis compared with women who completed >high school (pOR 1.15; 95% CI 1.03-1.28). The risk estimates for the different histotypes of ovarian cancer resembled that observed for ovarian cancers combined but did not reach statistical significance. Our results were unchanged when we included BMI and cigarette smoking.

CONCLUSION

Lower level of education was associated with an increased risk of advanced tumour stage at diagnosis of ovarian cancer. The observed socioeconomic difference in stage at diagnosis of ovarian cancer calls for further studies on how to reduce this diagnostic delay.

PPM1D Mosaic Truncating Variants in Ovarian Cancer Cases May Be Treatment-Related Somatic Mutations

Mosaic truncating mutations in the protein phosphatase, Mg(2+)/Mn(2+)-dependent, 1D (PPM1D) gene have recently been reported with a statistically significantly greater frequency in lymphocyte DNA from ovarian cancer case patients compared with unaffected control patients. Using massively parallel sequencing (MPS) we identified truncating PPM1D mutations in 12 of 3236 epithelial ovarian cancer (EOC) case patients (0.37%) but in only one of 3431 unaffected control patients (0.03%) (P = .001). All statistical tests were two-sided. A combination of Sanger sequencing, pyrosequencing, and MPS data suggested that 12 of the 13 mutations were mosaic. All mutations were identified in post-chemotherapy treatment blood samples from case patients (n = 1827) (average 1234 days post-treatment in carriers) rather than from cases collected pretreatment (less than 14 days after diagnosis, n = 1384) (P = .002). These data suggest that PPM1D variants in EOC cases are primarily somatic mosaic mutations caused by treatment and are not associated with germline predisposition to EOC.

Identification of Susceptibility Loci for Cutaneous Squamous Cell Carcinoma

We report a genome-wide association study of cutaneous squamous cell carcinoma conducted among non-Hispanic white members of the Kaiser Permanente Northern California health care system. The study includes a genome-wide screen of 61,457 members (6,891 cases and 54,566 controls) genotyped on the Affymetrix Axiom European array and a replication phase involving an independent set of 6,410 additional members (810 cases and 5,600 controls). Combined analysis of screening and replication phases identified 10 loci containing single-nucleotide polymorphisms (SNPs) with P-values < 5 × 10(-8). Six loci contain genes in the pigmentation pathway; SNPs at these loci appear to modulate squamous cell carcinoma risk independently of the pigmentation phenotypes. Another locus contains HLA class II genes studied in relation to elevated squamous cell carcinoma risk following immunosuppression. SNPs at the remaining three loci include an intronic SNP in FOXP1 at locus 3p13, an intergenic SNP at 3q28 near TP63, and an intergenic SNP at 9p22 near BNC2. These findings provide insights into the genetic factors accounting for inherited squamous cell carcinoma susceptibility.

Epithelial-Mesenchymal Transition (EMT) Gene Variants and Epithelial Ovarian Cancer (EOC) Risk

Epithelial-mesenchymal transition (EMT) is a process whereby epithelial cells assume mesenchymal characteristics to facilitate cancer metastasis. However, EMT also contributes to the initiation and development of primary tumors. Prior studies that explored the hypothesis that EMT gene variants contribute to epithelial ovarian carcinoma (EOC) risk have been based on small sample sizes and none have sought replication in an independent population. We screened 15,816 single-nucleotide polymorphisms (SNPs) in 296 genes in a discovery phase using data from a genome-wide association study of EOC among women of European ancestry (1,947 cases and 2,009 controls) and identified 793 variants in 278 EMT-related genes that were nominally (P < 0.05) associated with invasive EOC. These SNPs were then genotyped in a larger study of 14,525 invasive-cancer patients and 23,447 controls. A P-value <0.05 and a false discovery rate (FDR) <0.2 were considered statistically significant. In the larger dataset, GPC6/GPC5 rs17702471 was associated with the endometrioid subtype among Caucasians (odds ratio (OR) = 1.16, 95% CI = 1.07-1.25, P = 0.0003, FDR = 0.19), whereas F8 rs7053448 (OR = 1.69, 95% CI = 1.27-2.24, P = 0.0003, FDR = 0.12), F8 rs7058826 (OR = 1.69, 95% CI = 1.27-2.24, P = 0.0003, FDR = 0.12), and CAPN13 rs1983383 (OR = 0.79, 95% CI = 0.69-0.90, P = 0.0005, FDR = 0.12) were associated with combined invasive EOC among Asians. In silico functional analyses revealed that GPC6/GPC5 rs17702471 coincided with DNA regulatory elements. These results suggest that EMT gene variants do not appear to play a significant role in the susceptibility to EOC.