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.

Common variants at the CHEK2 gene locus and risk of epithelial ovarian cancer

Genome-wide association studies have identified 20 genomic regions associated with risk of epithelial ovarian cancer (EOC), but many additional risk variants may exist. Here, we evaluated associations between common genetic variants [single nucleotide polymorphisms (SNPs) and indels] in DNA repair genes and EOC risk. We genotyped 2896 common variants at 143 gene loci in DNA samples from 15 397 patients with invasive EOC and controls. We found evidence of associations with EOC risk for variants at FANCA, EXO1, E2F4, E2F2, CREB5 and CHEK2 genes (P ≤ 0.001). The strongest risk association was for CHEK2 SNP rs17507066 with serous EOC (P = 4.74 x 10(-7)). Additional genotyping and imputation of genotypes from the 1000 genomes project identified a slightly more significant association for CHEK2 SNP rs6005807 (r (2) with rs17507066 = 0.84, odds ratio (OR) 1.17, 95% CI 1.11-1.24, P = 1.1×10(-7)). We identified 293 variants in the region with likelihood ratios of less than 1:100 for representing the causal variant. Functional annotation identified 25 candidate SNPs that alter transcription factor binding sites within regulatory elements active in EOC precursor tissues. In The Cancer Genome Atlas dataset, CHEK2 gene expression was significantly higher in primary EOCs compared to normal fallopian tube tissues (P = 3.72×10(-8)). We also identified an association between genotypes of the candidate causal SNP rs12166475 (r (2) = 0.99 with rs6005807) and CHEK2 expression (P = 2.70×10(-8)). These data suggest that common variants at 22q12.1 are associated with risk of serous EOC and CHEK2 as a plausible target susceptibility gene.

Shared genetics underlying epidemiological association between endometriosis and ovarian cancer

Epidemiological studies have demonstrated associations between endometriosis and certain histotypes of ovarian cancer, including clear cell, low-grade serous and endometrioid carcinomas. We aimed to determine whether the observed associations might be due to shared genetic aetiology. To address this, we used two endometriosis datasets genotyped on common arrays with full-genome coverage (3194 cases and 7060 controls) and a large ovarian cancer dataset genotyped on the customized Illumina Infinium iSelect (iCOGS) arrays (10 065 cases and 21 663 controls). Previous work has suggested that a large number of genetic variants contribute to endometriosis and ovarian cancer (all histotypes combined) susceptibility. Here, using the iCOGS data, we confirmed polygenic architecture for most histotypes of ovarian cancer. This led us to evaluate if the polygenic effects are shared across diseases. We found evidence for shared genetic risks between endometriosis and all histotypes of ovarian cancer, except for the intestinal mucinous type. Clear cell carcinoma showed the strongest genetic correlation with endometriosis (0.51, 95% CI = 0.18-0.84). Endometrioid and low-grade serous carcinomas had similar correlation coefficients (0.48, 95% CI = 0.07-0.89 and 0.40, 95% CI = 0.05-0.75, respectively). High-grade serous carcinoma, which often arises from the fallopian tubes, showed a weaker genetic correlation with endometriosis (0.25, 95% CI = 0.11-0.39), despite the absence of a known epidemiological association. These results suggest that the epidemiological association between endometriosis and ovarian adenocarcinoma may be attributable to shared genetic susceptibility loci.

Network-Based Integration of GWAS and Gene Expression Identifies a HOX-Centric Network Associated with Serous Ovarian Cancer Risk

BACKGROUND

Genome-wide association studies (GWAS) have so far reported 12 loci associated with serous epithelial ovarian cancer (EOC) risk. We hypothesized that some of these loci function through nearby transcription factor (TF) genes and that putative target genes of these TFs as identified by coexpression may also be enriched for additional EOC risk associations.

METHODS

We selected TF genes within 1 Mb of the top signal at the 12 genome-wide significant risk loci. Mutual information, a form of correlation, was used to build networks of genes strongly coexpressed with each selected TF gene in the unified microarray dataset of 489 serous EOC tumors from The Cancer Genome Atlas. Genes represented in this dataset were subsequently ranked using a gene-level test based on results for germline SNPs from a serous EOC GWAS meta-analysis (2,196 cases/4,396 controls).

RESULTS

Gene set enrichment analysis identified six networks centered on TF genes (HOXB2, HOXB5, HOXB6, HOXB7 at 17q21.32 and HOXD1, HOXD3 at 2q31) that were significantly enriched for genes from the risk-associated end of the ranked list (P < 0.05 and FDR < 0.05). These results were replicated (P < 0.05) using an independent association study (7,035 cases/21,693 controls). Genes underlying enrichment in the six networks were pooled into a combined network.

CONCLUSION

We identified a HOX-centric network associated with serous EOC risk containing several genes with known or emerging roles in serous EOC development.

IMPACT

Network analysis integrating large, context-specific datasets has the potential to offer mechanistic insights into cancer susceptibility and prioritize genes for experimental characterization.

Cis-eQTL analysis and functional validation of candidate susceptibility genes for high-grade serous ovarian cancer

Genome-wide association studies have reported 11 regions conferring risk of high-grade serous epithelial ovarian cancer (HGSOC). Expression quantitative trait locus (eQTL) analyses can identify candidate susceptibility genes at risk loci. Here we evaluate cis-eQTL associations at 47 regions associated with HGSOC risk (P≤10(-5)). For three cis-eQTL associations (P<1.4 × 10(-3), FDR<0.05) at 1p36 (CDC42), 1p34 (CDCA8) and 2q31 (HOXD9), we evaluate the functional role of each candidate by perturbing expression of each gene in HGSOC precursor cells. Overexpression of HOXD9 increases anchorage-independent growth, shortens population-doubling time and reduces contact inhibition. Chromosome conformation capture identifies an interaction between rs2857532 and the HOXD9 promoter, suggesting this SNP is a leading causal variant. Transcriptomic profiling after HOXD9 overexpression reveals enrichment of HGSOC risk variants within HOXD9 target genes (P=6 × 10(-10) for risk variants (P<10(-4)) within 10 kb of a HOXD9 target gene in ovarian cells), suggesting a broader role for this network in genetic susceptibility to HGSOC.

Common Genetic Variation in Circadian Rhythm Genes and Risk of Epithelial Ovarian Cancer (EOC)

Disruption in circadian gene expression, whether due to genetic variation or environmental factors (e.g., light at night, shiftwork), is associated with increased incidence of breast, prostate, gastrointestinal and hematologic cancers and gliomas. Circadian genes are highly expressed in the ovaries where they regulate ovulation; circadian disruption is associated with several ovarian cancer risk factors (e.g., endometriosis). However, no studies have examined variation in germline circadian genes as predictors of ovarian cancer risk and invasiveness. The goal of the current study was to examine single nucleotide polymorphisms (SNPs) in circadian genes BMAL1, CRY2, CSNK1E, NPAS2, PER3, REV1 and TIMELESS and downstream transcription factors KLF10 and SENP3 as predictors of risk of epithelial ovarian cancer (EOC) and histopathologic subtypes. The study included a test set of 3,761 EOC cases and 2,722 controls and a validation set of 44,308 samples including 18,174 (10,316 serous) cases and 26,134 controls from 43 studies participating in the Ovarian Cancer Association Consortium (OCAC). Analysis of genotype data from 36 genotyped SNPs and 4600 imputed SNPs indicated that the most significant association was rs117104877 in BMAL1 (OR = 0.79, 95% CI = 0.68-0.90, p = 5.59 × 10-4]. Functional analysis revealed a significant down regulation of BMAL1 expression following cMYC overexpression and increasing transformation in ovarian surface epithelial (OSE) cells as well as alternative splicing of BMAL1 exons in ovarian and granulosa cells. These results suggest that variation in circadian genes, and specifically BMAL1, may be associated with risk of ovarian cancer, likely through disruption of hormonal pathways.

Abstract 2792: Polymorphisms in vitamin D-related genes and risk of breast cancer

Vitamin D is hypothesized to decrease breast cancer risk by controlling cell proliferation, differentiation and apoptosis, and estrogen synthesis and signalling in breast tissue. Transport of circulating vitamin D to sites of action is facilitated by vitamin D binding protein, encoded by the GC gene, whereupon it binds to the vitamin D receptor, encoded by the VDR gene. Given that polymorphisms in the GC and VDR gene have not been consistently associated with breast cancer risk, our objective was to assess the associations of vitamin D-related genes with breast cancer risk, accounting for ethnicity and menopausal status, and to determine these associations according to tumour receptor status defined by expression of estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor-2 (HER2). A case-control study of 1,037 breast cancer cases and 1,050 controls age-frequency matched in Vancouver, British Columbia and Kingston, Ontario was used. The distribution and associations of 21 single nucleotide polymorphisms (SNPs) of the GC and VDR gene with breast cancer risk and breast tumour subtypes were investigated among premenopausal and postmenopausal European and Asian women in stratified analyses. Following false discovery rate adjustment, no SNP was associated with breast cancer risk among any subgroup. However, among European women, two SNPs in the VDR gene, rs1544410 (OR = 0.38, 95% CI: 0.21 - 0.70) and rs7967152 (OR = 2.80, 95% CI: 1.62 - 4.85) were associated with risk of ER-/PR-/HER2+ breast tumour subtype, but not with other breast tumour subtypes. Breast tumour subtypes are proposed to be etiologically distinct disease states, with different risk factors and prognostic outcomes. Interactions between variations in vitamin D-related genes and breast tumour subtypes are biologically plausible given our current understanding of anticarcinogenic properties of vitamin D elicited through the regulation of cell growth and reproductive hormone synthesis. These findings provide additional insight into the etiologic differences among breast tumour subtypes.

Citation Format: Joy Shi, Anne Grundy, Harriet Richardson, Angela Brooks-Wilson, John Spinelli, Kristan J. Aronson. Polymorphisms in vitamin D-related genes and risk of breast cancer. [abstract]. In: Proceedings of the 106th Annual Meeting of the American Association for Cancer Research; 2015 Apr 18-22; Philadelphia, PA. Philadelphia (PA): AACR; Cancer Res 2015;75(15 Suppl):Abstract nr 2792. doi:10.1158/1538-7445.AM2015-2792

Common Genetic Variation In Cellular Transport Genes and Epithelial Ovarian Cancer (EOC) Risk

BACKGROUND

Defective cellular transport processes can lead to aberrant accumulation of trace elements, iron, small molecules and hormones in the cell, which in turn may promote the formation of reactive oxygen species, promoting DNA damage and aberrant expression of key regulatory cancer genes. As DNA damage and uncontrolled proliferation are hallmarks of cancer, including epithelial ovarian cancer (EOC), we hypothesized that inherited variation in the cellular transport genes contributes to EOC risk.

METHODS

In total, DNA samples were obtained from 14,525 case subjects with invasive EOC and from 23,447 controls from 43 sites in the Ovarian Cancer Association Consortium (OCAC). Two hundred seventy nine SNPs, representing 131 genes, were genotyped using an Illumina Infinium iSelect BeadChip as part of the Collaborative Oncological Gene-environment Study (COGS). SNP analyses were conducted using unconditional logistic regression under a log-additive model, and the FDR q<0.2 was applied to adjust for multiple comparisons.

RESULTS

The most significant evidence of an association for all invasive cancers combined and for the serous subtype was observed for SNP rs17216603 in the iron transporter gene HEPH (invasive: OR = 0.85, P = 0.00026; serous: OR = 0.81, P = 0.00020); this SNP was also associated with the borderline/low malignant potential (LMP) tumors (P = 0.021). Other genes significantly associated with EOC histological subtypes (p<0.05) included the UGT1A (endometrioid), SLC25A45 (mucinous), SLC39A11 (low malignant potential), and SERPINA7 (clear cell carcinoma). In addition, 1785 SNPs in six genes (HEPH, MGST1, SERPINA, SLC25A45, SLC39A11 and UGT1A) were imputed from the 1000 Genomes Project and examined for association with INV EOC in white-European subjects. The most significant imputed SNP was rs117729793 in SLC39A11 (per allele, OR = 2.55, 95% CI = 1.5-4.35, p = 5.66x10-4).

CONCLUSION

These results, generated on a large cohort of women, revealed associations between inherited cellular transport gene variants and risk of EOC histologic subtypes.

Etiologic heterogeneity among non-Hodgkin lymphoma subtypes: the InterLymph Non-Hodgkin Lymphoma Subtypes Project

BACKGROUND

Non-Hodgkin lymphoma (NHL) comprises biologically and clinically heterogeneous subtypes. Previously, study size has limited the ability to compare and contrast the risk factor profiles among these heterogeneous subtypes.

METHODS

We pooled individual-level data from 17 471 NHL cases and 23 096 controls in 20 case-control studies from the International Lymphoma Epidemiology Consortium (InterLymph). We estimated the associations, measured as odds ratios, between each of 11 NHL subtypes and self-reported medical history, family history of hematologic malignancy, lifestyle factors, and occupation. We then assessed the heterogeneity of associations by evaluating the variability (Q value) of the estimated odds ratios for a given exposure among subtypes. Finally, we organized the subtypes into a hierarchical tree to identify groups that had similar risk factor profiles. Statistical significance of tree partitions was estimated by permutation-based P values (P NODE).

RESULTS

Risks differed statistically significantly among NHL subtypes for medical history factors (autoimmune diseases, hepatitis C virus seropositivity, eczema, and blood transfusion), family history of leukemia and multiple myeloma, alcohol consumption, cigarette smoking, and certain occupations, whereas generally homogeneous risks among subtypes were observed for family history of NHL, recreational sun exposure, hay fever, allergy, and socioeconomic status. Overall, the greatest difference in risk factors occurred between T-cell and B-cell lymphomas (P NODE < 1.0×10(-4)), with increased risks generally restricted to T-cell lymphomas for eczema, T-cell-activating autoimmune diseases, family history of multiple myeloma, and occupation as a painter. We further observed substantial heterogeneity among B-cell lymphomas (P NODE < 1.0×10(-4)). Increased risks for B-cell-activating autoimmune disease and hepatitis C virus seropositivity and decreased risks for alcohol consumption and occupation as a teacher generally were restricted to marginal zone lymphoma, Burkitt/Burkitt-like lymphoma/leukemia, diffuse large B-cell lymphoma, and/or lymphoplasmacytic lymphoma/Waldenström macroglobulinemia.

CONCLUSIONS

Using a novel approach to investigate etiologic heterogeneity among NHL subtypes, we identified risk factors that were common among subtypes as well as risk factors that appeared to be distinct among individual or a few subtypes, suggesting both subtype-specific and shared underlying mechanisms. Further research is needed to test putative mechanisms, investigate other risk factors (eg, other infections, environmental exposures, and diet), and evaluate potential joint effects with genetic susceptibility.

Medical history, lifestyle, family history, and occupational risk factors for lymphoplasmacytic lymphoma/Waldenström's macroglobulinemia: the InterLymph Non-Hodgkin Lymphoma Subtypes Project

BACKGROUND

Lymphoplasmacytic lymphoma/Waldenström's macroglobulinemia (LPL/WM), a rare non-Hodgkin lymphoma subtype, shows strong familial aggregation and a positive association with chronic immune stimulation, but evidence regarding other risk factors is very limited.

METHODS

The International Lymphoma Epidemiology Consortium (InterLymph) pooled data from 11 predominantly population-based case-control studies from North America, Europe, and Australia to examine medical history, lifestyle, family history, and occupational risk factors for LPL/WM. Age-, sex-, race/ethnicity-, and study-adjusted odds ratios (ORs) and 95% confidence intervals (CIs) were estimated using logistic regression for a total of 374 LPL/WM cases and 23 096 controls.

RESULTS

In multivariate analysis including all putative risk factors, LPL/WM risk was associated with history of Sjögren's syndrome (OR = 14.0, 95% CI = 3.60 to 54.6), systemic lupus erythematosus (OR = 8.23, 95% CI = 2.69 to 25.2), hay fever (OR = 0.73, 95% CI = 0.54 to 0.99), positive hepatitis C serology (OR = 2.51, 95% CI = 1.03 to 6.17), hematologic malignancy in a first-degree relative (OR = 1.64, 95% CI = 1.02 to 2.64), adult weight (OR = 0.61, 95% CI = 0.44 to 0.85 for highest vs. lowest quartile), duration of cigarette smoking (OR = 1.46, 95% CI = 1.04 to 2.05 for ≥ 40 years vs. nonsmokers), and occupation as a medical doctor (OR = 5.54, 95% CI = 2.19 to 14.0). There was no association with other medical conditions, lifestyle factors, or occupations.

CONCLUSIONS

This pooled analysis confirmed associations with immune conditions and family history of hematologic malignancy, and identified new associations with hay fever, weight, smoking, and occupation, and no association with other lifestyle factors. These findings offer clues to LPL/WM biology and prevention.

Associations of non-Hodgkin Lymphoma (NHL) risk with autoimmune conditions according to putative NHL loci

Autoimmune conditions and immune system-related genetic variations are associated with risk of non-Hodgkin lymphoma (NHL). In a pooled analysis of 8,692 NHL cases and 9,260 controls from 14 studies (1988-2007) within the International Lymphoma Epidemiology Consortium, we evaluated the interaction between immune system genetic variants and autoimmune conditions in NHL risk. We evaluated the immunity-related single nucleotide polymorphisms rs1800629 (tumor necrosis factor gene (TNF) G308A), rs1800890 (interleukin-10 gene (IL10) T3575A), rs6457327 (human leukocyte antigen gene (HLA) class I), rs10484561 (HLA class II), and rs2647012 (HLA class II)) and categorized autoimmune conditions as primarily mediated by B-cell or T-cell responses. We constructed unconditional logistic regression models to measure associations between autoimmune conditions and NHL with stratification by genotype. Autoimmune conditions mediated by B-cell responses were associated with increased NHL risk, specifically diffuse large B-cell lymphoma (odds ratio (OR) = 3.11, 95% confidence interval (CI): 2.25, 4.30) and marginal zone lymphoma (OR = 5.80, 95% CI: 3.82, 8.80); those mediated by T-cell responses were associated with peripheral T-cell lymphoma (OR = 2.14, 95% CI: 1.35, 3.38). In the presence of the rs1800629 AG/AA genotype, B-cell-mediated autoimmune conditions increased NHL risk (OR = 3.27, 95% CI: 2.07, 5.16; P-interaction = 0.03) in comparison with the GG genotype (OR = 1.82, 95% CI: 1.31, 2.53). This interaction was consistent across major B-cell NHL subtypes, including marginal zone lymphoma (P-interaction = 0.02) and follicular lymphoma (P-interaction = 0.04).

Evaluating the ovarian cancer gonadotropin hypothesis: a candidate gene study

OBJECTIVE

Ovarian cancer is a hormone-related disease with a strong genetic basis. However, none of its high-penetrance susceptibility genes and GWAS-identified variants to date are known to be involved in hormonal pathways. Given the hypothesized etiologic role of gonadotropins, an assessment of how variability in genes involved in the gonadotropin signaling pathway impacts disease risk is warranted.

METHODS

Genetic data from 41 ovarian cancer study sites were pooled and unconditional logistic regression was used to evaluate whether any of the 2185 SNPs from 11 gonadotropin signaling pathway genes was associated with ovarian cancer risk. A burden test using the admixture likelihood (AML) method was also used to evaluate gene-level associations.

RESULTS

We did not find any genome-wide significant associations between individual SNPs and ovarian cancer risk. However, there was some suggestion of gene-level associations for four gonadotropin signaling pathway genes: INHBB (p=0.045, mucinous), LHCGR (p=0.046, high-grade serous), GNRH (p=0.041, high-grade serous), and FSHB (p=0.036, overall invasive). There was also suggestive evidence for INHA (p=0.060, overall invasive).

CONCLUSIONS

Ovarian cancer studies have limited sample numbers, thus fewer genome-wide susceptibility alleles, with only modest associations, have been identified relative to breast and prostate cancers. We have evaluated the majority of ovarian cancer studies with biological samples, to our knowledge, leaving no opportunity for replication. Using both our understanding of biology and powerful gene-level tests, we have identified four putative ovarian cancer loci near INHBB, LHCGR, GNRH, and FSHB that warrant a second look if larger sample sizes and denser genotype chips become available.

Identification of six new susceptibility loci for invasive epithelial ovarian cancer

Genome-wide association studies (GWAS) have identified 12 epithelial ovarian cancer (EOC) susceptibility alleles. The pattern of association at these loci is consistent in BRCA1 and BRCA2 mutation carriers who are at high risk of EOC. After imputation to 1000 Genomes Project data, we assessed associations of 11 million genetic variants with EOC risk from 15,437 cases unselected for family history and 30,845 controls and from 15,252 BRCA1 mutation carriers and 8,211 BRCA2 mutation carriers (3,096 with ovarian cancer), and we combined the results in a meta-analysis. This new study design yielded increased statistical power, leading to the discovery of six new EOC susceptibility loci. Variants at 1p36 (nearest gene, WNT4), 4q26 (SYNPO2), 9q34.2 (ABO) and 17q11.2 (ATAD5) were associated with EOC risk, and at 1p34.3 (RSPO1) and 6p22.1 (GPX6) variants were specifically associated with the serous EOC subtype, all with P < 5 × 10(-8). Incorporating these variants into risk assessment tools will improve clinical risk predictions for BRCA1 and BRCA2 mutation carriers.

Consortium analysis of gene and gene-folate interactions in purine and pyrimidine metabolism pathways with ovarian carcinoma risk

SCOPE

We reevaluated previously reported associations between variants in pathways of one-carbon (1-C) (folate) transfer genes and ovarian carcinoma (OC) risk, and in related pathways of purine and pyrimidine metabolism, and assessed interactions with folate intake.

METHODS AND RESULTS

Odds ratios (OR) for 446 genetic variants were estimated among 13,410 OC cases and 22,635 controls, and among 2281 cases and 3444 controls with folate information. Following multiple testing correction, the most significant main effect associations were for dihydropyrimidine dehydrogenase (DPYD) variants rs11587873 (OR = 0.92; p = 6 × 10⁻⁵) and rs828054 (OR = 1.06; p = 1 × 10⁻⁴). Thirteen variants in the pyrimidine metabolism genes, DPYD, DPYS, PPAT, and TYMS, also interacted significantly with folate in a multivariant analysis (corrected p = 9.9 × 10⁻⁶) but collectively explained only 0.2% of OC risk. Although no other associations were significant after multiple testing correction, variants in SHMT1 in 1-C transfer, previously reported with OC, suggested lower risk at higher folate (p(interaction) = 0.03-0.006).

CONCLUSION

Variation in pyrimidine metabolism genes, particularly DPYD, which was previously reported to be associated with OC, may influence risk; however, stratification by folate intake is unlikely to modify disease risk appreciably in these women. SHMT1 SNP-by-folate interactions are plausible but require further validation. Polymorphisms in selected genes in purine metabolism were not associated with OC.

Genome-wide association study of subtype-specific epithelial ovarian cancer risk alleles using pooled DNA

Epithelial ovarian cancer (EOC) is a heterogeneous cancer with both genetic and environmental risk factors. Variants influencing the risk of developing the less-common EOC subtypes have not been fully investigated. We performed a genome-wide association study (GWAS) of EOC according to subtype by pooling genomic DNA from 545 cases and 398 controls of European descent, and testing for allelic associations. We evaluated for replication 188 variants from the GWAS [56 variants for mucinous, 55 for endometrioid and clear cell, 53 for low-malignant potential (LMP) serous, and 24 for invasive serous EOC], selected using pre-defined criteria. Genotypes from 13,188 cases and 23,164 controls of European descent were used to perform unconditional logistic regression under the log-additive genetic model; odds ratios (OR) and 95 % confidence intervals are reported. Nine variants tagging six loci were associated with subtype-specific EOC risk at P < 0.05, and had an OR that agreed in direction of effect with the GWAS results. Several of these variants are in or near genes with a biological rationale for conferring EOC risk, including ZFP36L1 and RAD51B for mucinous EOC (rs17106154, OR = 1.17, P = 0.029, n = 1,483 cases), GRB10 for endometrioid and clear cell EOC (rs2190503, P = 0.014, n = 2,903 cases), and C22orf26/BPIL2 for LMP serous EOC (rs9609538, OR = 0.86, P = 0.0043, n = 892 cases). In analyses that included the 75 GWAS samples, the association between rs9609538 (OR = 0.84, P = 0.0007) and LMP serous EOC risk remained statistically significant at P < 0.0012 adjusted for multiple testing. Replication in additional samples will be important to verify these results for the less-common EOC subtypes.

Gender and telomere length: systematic review and meta-analysis

BACKGROUND

It is widely believed that females have longer telomeres than males, although results from studies have been contradictory.

METHODS

We carried out a systematic review and meta-analyses to test the hypothesis that in humans, females have longer telomeres than males and that this association becomes stronger with increasing age. Searches were conducted in EMBASE and MEDLINE (by November 2009) and additional datasets were obtained from study investigators. Eligible observational studies measured telomeres for both females and males of any age, had a minimum sample size of 100 and included participants not part of a diseased group. We calculated summary estimates using random-effects meta-analyses. Heterogeneity between studies was investigated using sub-group analysis and meta-regression.

RESULTS

Meta-analyses from 36 cohorts (36,230 participants) showed that on average females had longer telomeres than males (standardised difference in telomere length between females and males 0.090, 95% CI 0.015, 0.166; age-adjusted). There was little evidence that these associations varied by age group (p=1.00) or cell type (p=0.29). However, the size of this difference did vary by measurement methods, with only Southern blot but neither real-time PCR nor Flow-FISH showing a significant difference. This difference was not associated with random measurement error.

CONCLUSIONS

Telomere length is longer in females than males, although this difference was not universally found in studies that did not use Southern blot methods. Further research on explanations for the methodological differences is required.

Risk of ovarian cancer and the NF-κB pathway: genetic association with IL1A and TNFSF10

A missense single-nucleotide polymorphism (SNP) in the immune modulatory gene IL1A has been associated with ovarian cancer risk (rs17561). Although the exact mechanism through which this SNP alters risk of ovarian cancer is not clearly understood, rs17561 has also been associated with risk of endometriosis, an epidemiologic risk factor for ovarian cancer. Interleukin-1α (IL1A) is both regulated by and able to activate NF-κB, a transcription factor family that induces transcription of many proinflammatory genes and may be an important mediator in carcinogenesis. We therefore tagged SNPs in more than 200 genes in the NF-κB pathway for a total of 2,282 SNPs (including rs17561) for genotype analysis of 15,604 cases of ovarian cancer in patients of European descent, including 6,179 of high-grade serous (HGS), 2,100 endometrioid, 1,591 mucinous, 1,034 clear cell, and 1,016 low-grade serous, including 23,235 control cases spanning 40 studies in the Ovarian Cancer Association Consortium. In this large population, we confirmed the association between rs17561 and clear cell ovarian cancer [OR, 0.84; 95% confidence interval (CI), 0.76-0.93; P = 0.00075], which remained intact even after excluding participants in the prior study (OR, 0.85; 95% CI, 0.75-0.95; P = 0.006). Considering a multiple-testing-corrected significance threshold of P < 2.5 × 10(-5), only one other variant, the TNFSF10 SNP rs6785617, was associated significantly with a risk of ovarian cancer (low malignant potential tumors OR, 0.85; 95% CI, 0.79-0.91; P = 0.00002). Our results extend the evidence that borderline tumors may have a distinct genetic etiology. Further investigation of how these SNPs might modify ovarian cancer associations with other inflammation-related risk factors is warranted.

Inherited common variants in mitochondrial DNA and invasive serous epithelial ovarian cancer risk

Background

Mitochondria are the site of oxidative phosphorylation, a process which generates reactive oxygen species (ROS). Elevated ROS levels can lead to oxidative stress, a cellular state implicated in carcinogenesis. It is hypothesized that alternations in mitochondrial (MT) DNA, including heritable MT single nucleotide polymorphisms (MT-SNPs), have the potential to change the capacity of MT function, leading to increased oxidative stress and cancer risk. We investigated if common MT-SNPs and/or haplogroups and are associated with invasive serous ovarian cancer (OvCa) risk.

Methods

A panel of 64 MT-SNPs designed to tag all common variation in the European MT genome (minor allele frequency (MAF) >1%, r^2 >0.8) was genotyped in study participants of European descent using the Sequenom MassARRAY iPlex Gold® system (Sequenom Inc, CA, USA). Invasive serous OvCa cases (n = 405) and frequency age-matched controls (n = 445) were drawn from a population-based case-control study of OvCa in western Canada. Binary logistic regression was used to estimate the odds ratio (OR) and 95% confidence intervals (C.I.) for carriage of the minor versus major allele by case-control status. MitoTool was used to test the relationship between European haplogroup status and case-control status using Fisher’s exact test.

Results

The most significant disease-SNP association was for rs2857285, a synonymous MT-SNP in ND4 (OR = 4.84, 95% CI: 1.03–22.68, P = 0.045). After adjustment for multiple testing using a Bonferroni correction of the Type 1 error this MT-SNP was not significant. No other MT-SNP had a P-value < 0.05. European haplogroup status was not associated with case status. Most MT-SNPs (73%) genotyped had a MAF <5%.

Conclusion

Common European MT-SNPs (MAF > 5%) and haplogroups were not associated with invasive serous OvCa risk in this study; however, most European MT-SNPs have a low MAF (<5%), which we were underpowered to adequately assess. Larger studies are needed to clarify the role of low MAF MT-SNPs (MAF < 5%) in invasive serous OvCa risk.

Shift work, circadian gene variants and risk of breast cancer

Circadian (clock) genes have been linked with several functions relevant to cancer, and epidemiologic research has suggested relationships with breast cancer risk for variants in NPAS2, CLOCK, CRY2 and TIMELESS. Increased breast cancer risk has also been observed among shift workers, suggesting potential interactions in relationships of circadian genes with breast cancer. Relationships with breast cancer of 100 SNPs in 14 clock-related genes, as well as potential interactions with shift work history, were investigated in a case-control study (1042 cases, 1051 controls). Odds ratios in an additive genetic model for European-ancestry participants (645 cases, 806 controls) were calculated, using a two-step correction for multiple testing: within each gene through permutation testing (10,000 permutations), and correcting for the false discovery rate across genes. Interactions of genotypes with ethnicity and shift work (<2 years vs ≥2 years) were evaluated individually. Following permutation analysis, two SNPs (rs3816360 in ARNTL and rs11113179 in CRY1) displayed significant associations with breast cancer and one SNP (rs3027188 in PER1) was marginally significant; however, none were significant following adjustment for the false discovery rate. No significant interaction with shift work history was detected. If shift work causes circadian disruption, this was not reflected in associations between clock gene variants and breast cancer risk in this study. Larger studies are needed to assess interactions with longer durations (>30 years) of shift work that have been associated with breast cancer.

Shared genetic factors for age at natural menopause in Iranian and European women

STUDY QUESTION

Do differences in heritable genetic factors explain some of the difference in age at natural menopause (ANM) among populations?

SUMMARY ANSWER

One single nucleotide polymorphism (SNP)-ANM association (rs16991615) detected in European women was replicated in Iranian women.

WHAT IS KNOWN ALREADY

Genetics plays an important role in ANM, and well-powered genome-wide association studies (GWAS) of ANM performed in European women have discovered many statistically significant SNP-ANM associations. Average ANM varies by ethnicity, and population-specific differences in ANM-associated alleles may in part explain these differences.

STUDY DESIGN, SIZE, DURATION

After quality control procedures, 97 SNPs were analyzed in genotype data of 828 Iranian women who experienced natural menopause. SNP genotyping data were used to perform linear regression analyses with ANM as a quantitative trait. Study participants were drawn from the population-based Tehran Lipid and Glucose Study based in Tehran, Iran. This study was performed between February 2009 and March 2012.

PARTICIPANTS/MATERIALS, SETTING AND METHODS

Based on an ANM-GWAS literature review, eight SNPs at four loci previously associated with ANM in European women were tested for replication in Iranian women. Linear regression analyses were performed including (n = 828) and excluding (n = 783) women who experience premature ovarian failure (ANM before 40 years of age). In addition, to search for novel population-specific ANM risk alleles, a pool-based GWAS was performed using this collection of Iranian women. Two DNA pools were constructed and compared: an 'early' ANM pool (lower 20(th) percentile of menopause ages, 40-45 years, n = 165) and a 'late' ANM pool (upper 20(th) percentile of menopause ages, 54-65 years, n = 187). Each DNA pool was assayed on four Illumina Human1M-Duo arrays, and allele-based tests of association were used to rank SNPs. One hundred and two highly ranked SNPs were chosen for individual genotyping by Sequenom MassARRAY and association analysis in the Iranian women.

MAIN RESULTS AND THE ROLE OF CHANCE

One SNP-ANM association previously detected in European women was replicated in Iranian women (rs16991615; β = 1.07, standard error (SE): 0.49, P = 0.02). SNPs at the previously reported 19q13.42 and 6p24.2 loci also approached statistical significance and had consistent SNP effects (magnitude and direction) in Iranian women (rs1172822; β = -0.39, SE: 0.22, P = 0.08; and rs2153157, β = 0.41, SE: 0.21, P = 0.05). We found little evidence for novel SNP-ANM associations in Iranian women; no SNP selected based on the pool-based GWAS achieved genome-wide significance.

LIMITATIONS, REASONS FOR CAUTION

Due to small sample size this study was powered to reliably detect only moderate-to-large SNP effect sizes. This limited our ability to replicate many of the previously reported SNP-ANM risk alleles and to discover novel SNP-ANM associations' specific to the Iranian population. In performing our pool-based GWAS, a reduction in power was introduced relative to a conventional GWAS.

WIDER IMPLICATIONS OF THE FINDINGS

Our results imply that European and Iranian women share ANM-associated genetic variants. Our study was underpowered but for all SNPs tested the direction of the effect was consistent with data from the European study. Therefore, we anticipate that many (if not all) of the ANM-associated SNPs discovered in European women will replicate in Iranian women upon genotyping a sufficient number of women. Our data do not support the hypothesis that population-specific SNP-ANM associations explain population-specific differences in the mean ANM.

GWAS meta-analysis and replication identifies three new susceptibility loci for ovarian cancer

Genome-wide association studies (GWAS) have identified four susceptibility loci for epithelial ovarian cancer (EOC), with another two suggestive loci reaching near genome-wide significance. We pooled data from a GWAS conducted in North America with another GWAS from the UK. We selected the top 24,551 SNPs for inclusion on the iCOGS custom genotyping array. We performed follow-up genotyping in 18,174 individuals with EOC (cases) and 26,134 controls from 43 studies from the Ovarian Cancer Association Consortium. We validated the two loci at 3q25 and 17q21 that were previously found to have associations close to genome-wide significance and identified three loci newly associated with risk: two loci associated with all EOC subtypes at 8q21 (rs11782652, P = 5.5 × 10(-9)) and 10p12 (rs1243180, P = 1.8 × 10(-8)) and another locus specific to the serous subtype at 17q12 (rs757210, P = 8.1 × 10(-10)). An integrated molecular analysis of genes and regulatory regions at these loci provided evidence for functional mechanisms underlying susceptibility and implicated CHMP4C in the pathogenesis of ovarian cancer.

Multiple independent variants at the TERT locus are associated with telomere length and risks of breast and ovarian cancer

TERT-locus SNPs and leukocyte telomere measures are reportedly associated with risks of multiple cancers. Using the Illumina custom genotyping array iCOGs, we analyzed ∼480 SNPs at the TERT locus in breast (n = 103,991), ovarian (n = 39,774) and BRCA1 mutation carrier (n = 11,705) cancer cases and controls. Leukocyte telomere measurements were also available for 53,724 participants. Most associations cluster into three independent peaks. The minor allele at the peak 1 SNP rs2736108 associates with longer telomeres (P = 5.8 × 10(-7)), lower risks for estrogen receptor (ER)-negative (P = 1.0 × 10(-8)) and BRCA1 mutation carrier (P = 1.1 × 10(-5)) breast cancers and altered promoter assay signal. The minor allele at the peak 2 SNP rs7705526 associates with longer telomeres (P = 2.3 × 10(-14)), higher risk of low-malignant-potential ovarian cancer (P = 1.3 × 10(-15)) and greater promoter activity. The minor alleles at the peak 3 SNPs rs10069690 and rs2242652 increase ER-negative (P = 1.2 × 10(-12)) and BRCA1 mutation carrier (P = 1.6 × 10(-14)) breast and invasive ovarian (P = 1.3 × 10(-11)) cancer risks but not via altered telomere length. The cancer risk alleles of rs2242652 and rs10069690, respectively, increase silencing and generate a truncated TERT splice variant.

Genetic polymorphisms at TIMP3 are associated with survival of adenocarcinoma of the gastroesophageal junction

The poor survival of adenocarcinomas of the gastroesophageal junction (GEJ) makes them clinically important. Discovery of host genetic factors that affect outcome may guide more individualized treatment. This study tests whether constitutional genetic variants in matrix metalloproteinases (MMP) and tissue inhibitors of metalloproteinases (TIMP) genes are associated with outcome of GEJ adenocarcinoma. Single nucleotide polymorphisms (SNPs) at four TIMP (TIMP1-4) and three MMP genes (MMP2, MMP7 and MMP9) were genotyped in DNA samples from a prospective cohort of patients with primary adenocarcinoma of the GEJ admitted to the British Columbia Cancer Agency. Cox proportional hazards regression, with adjustment for patient, disease and treatment variables, was used to estimate the association of SNPs with survival. Genotypes for 85 samples and 48 SNPs were analyzed. Four SNPs across TIMP3, (rs130274, rs715572, rs1962223 and rs5754312) were associated with survival. Interaction analyses revealed that the survival associations with rs715572 and rs5754312 are specific and significant for 5FU+cisplatin treated patients. Sanger sequencing of the TIMP3 coding and promoter regions revealed an additional SNP, rs9862, also associated with survival. TIMP3 genetic variants are associated with survival and may be potentially useful in optimizing treatment strategies for individual patients.

Identification and molecular characterization of a new ovarian cancer susceptibility locus at 17q21.31

Epithelial ovarian cancer (EOC) has a heritable component that remains to be fully characterized. Most identified common susceptibility variants lie in non-protein-coding sequences. We hypothesized that variants in the 3' untranslated region at putative microRNA (miRNA)-binding sites represent functional targets that influence EOC susceptibility. Here, we evaluate the association between 767 miRNA-related single-nucleotide polymorphisms (miRSNPs) and EOC risk in 18,174 EOC cases and 26,134 controls from 43 studies genotyped through the Collaborative Oncological Gene-environment Study. We identify several miRSNPs associated with invasive serous EOC risk (odds ratio=1.12, P=10(-8)) mapping to an inversion polymorphism at 17q21.31. Additional genotyping of non-miRSNPs at 17q21.31 reveals stronger signals outside the inversion (P=10(-10)). Variation at 17q21.31 is associated with neurological diseases, and our collaboration is the first to report an association with EOC susceptibility. An integrated molecular analysis in this region provides evidence for ARHGAP27 and PLEKHM1 as candidate EOC susceptibility genes.