Physiological phenotypes have optimal values relevant to healthy aging: sweet spots deduced from the Canadian Longitudinal Study on Aging

Previous observations on a group of exceptionally healthy "Super-Seniors" showed a lower variance of multiple physiological measures relevant for health than did a less healthy group of the same age. The finding was interpreted as the healthier individuals having physiological measurement values closer to an optimal level, or "sweet spot." Here, we tested the generalizability of the sweet-spot hypothesis in a larger community sample, comparing differences in the variance between healthier and less healthy groups. We apply this method to the Canadian Longitudinal Study on Aging (CLSA) comprehensive cohort of 30,097 participants aged 45 to 85 years with deep phenotype data. Data from both sexes and four age ranges were analyzed. Five instruments were used to represent different aspects of health, physical, and cognitive functioning. We tested 231 phenotypic measures for lower variance in the most healthy vs. least healthy quartile of each sex and age group, as classified by the five instruments. Segmented regression was used to determine sex-specific optimal values. One hundred forty-two physiological measures (61%) showed lower variance in the healthiest than in the least healthy group, in at least one sex and age group. The difference in variance was most significant for hemoglobin A1c and was also significant for many body composition measurements, but not for bone mineral density. Ninety-four phenotypes showed a nonmonotonic relationship with health, consistent with the idea of a sweet spot; for these, we determined optimal values and 95% confidence intervals that were generally narrower than the ranges of current clinical reference intervals. These findings for sweet spot discovery validate the proposed approach for identifying traits important for healthy aging.

Concurrent RB1 loss and BRCA-deficiency predicts enhanced immunological response and long-term survival in tubo-ovarian high-grade serous carcinoma

Background

Somatic loss of the tumour suppressor RB1 is a common event in tubo-ovarian high-grade serous carcinoma (HGSC), which frequently co-occurs with alterations in homologous recombination DNA repair genes including BRCA1 and BRCA2 (BRCA). We examined whether tumour expression of RB1 was associated with survival across ovarian cancer histotypes (HGSC, endometrioid (ENOC), clear cell (CCOC), mucinous (MOC), low-grade serous carcinoma (LGSC)), and how co-occurrence of germline BRCA pathogenic variants and RB1 loss influences long-term survival in a large series of HGSC.

Patients and methods

RB1 protein expression patterns were classified by immunohistochemistry in epithelial ovarian carcinomas of 7436 patients from 20 studies participating in the Ovarian Tumor Tissue Analysis consortium and assessed for associations with overall survival (OS), accounting for patient age at diagnosis and FIGO stage. We examined RB1 expression and germline BRCA status in a subset of 1134 HGSC, and related genotype to survival, tumour infiltrating CD8+ lymphocyte counts and transcriptomic subtypes. Using CRISPR-Cas9, we deleted RB1 in HGSC cell lines with and without BRCA1 mutations to model co-loss with treatment response. We also performed genomic analyses on 126 primary HGSC to explore the molecular characteristics of concurrent homologous recombination deficiency and RB1 loss.

Results

RB1 protein loss was most frequent in HGSC (16.4%) and was highly correlated with RB1 mRNA expression. RB1 loss was associated with longer OS in HGSC (hazard ratio [HR] 0.74, 95% confidence interval [CI] 0.66-0.83, P = 6.8 ×10-7), but with poorer prognosis in ENOC (HR 2.17, 95% CI 1.17-4.03, P = 0.0140). Germline BRCA mutations and RB1 loss co-occurred in HGSC (P < 0.0001). Patients with both RB1 loss and germline BRCA mutations had a superior OS (HR 0.38, 95% CI 0.25-0.58, P = 5.2 ×10-6) compared to patients with either alteration alone, and their median OS was three times longer than non-carriers whose tumours retained RB1 expression (9.3 years vs. 3.1 years). Enhanced sensitivity to cisplatin (P < 0.01) and paclitaxel (P < 0.05) was seen in BRCA1 mutated cell lines with RB1 knockout. Among 126 patients with whole-genome and transcriptome sequence data, combined RB1 loss and genomic evidence of homologous recombination deficiency was correlated with transcriptional markers of enhanced interferon response, cell cycle deregulation, and reduced epithelial-mesenchymal transition in primary HGSC. CD8+ lymphocytes were most prevalent in BRCA-deficient HGSC with co-loss of RB1.

Conclusions

Co-occurrence of RB1 loss and BRCA mutation was associated with exceptionally long survival in patients with HGSC, potentially due to better treatment response and immune stimulation.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry

BACKGROUND

Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer.

METHODS

Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs.

RESULTS

Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94).

CONCLUSIONS

Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

Lifetime ovulatory years and risk of epithelial ovarian cancer: a multinational pooled analysis

BACKGROUND

The role of ovulation in epithelial ovarian cancer (EOC) is supported by the consistent protective effects of parity and oral contraceptive use. Whether these factors protect through anovulation alone remains unclear. We explored the association between lifetime ovulatory years (LOY) and EOC.

METHODS

LOY was calculated using 12 algorithms. Odds ratios (ORs) and 95% confidence intervals (CIs) estimated the association between LOY or LOY components and EOC among 26 204 control participants and 21 267 case patients from 25 studies. To assess whether LOY components act through ovulation suppression alone, we compared beta coefficients obtained from regression models with expected estimates assuming 1 year of ovulation suppression has the same effect regardless of source.

RESULTS

LOY was associated with increased EOC risk (OR per year increase = 1.014, 95% CI = 1.009 to 1.020 to OR per year increase = 1.044, 95% CI = 1.041 to 1.048). Individual LOY components, except age at menarche, also associated with EOC. The estimated model coefficient for oral contraceptive use and pregnancies were 4.45 times and 12- to 15-fold greater than expected, respectively. LOY was associated with high-grade serous, low-grade serous, endometrioid, and clear cell histotypes (ORs per year increase = 1.054, 1.040, 1.065, and 1.098, respectively) but not mucinous tumors. Estimated coefficients of LOY components were close to expected estimates for high-grade serous but larger than expected for low-grade serous, endometrioid, and clear cell histotypes.

CONCLUSIONS

LOY is positively associated with nonmucinous EOC. Differences between estimated and expected model coefficients for LOY components suggest factors beyond ovulation underlie the associations between LOY components and EOC in general and for non-HGSOC.

Increased FOXJ1 protein expression is associated with improved overall survival in high-grade serous ovarian carcinoma: an Ovarian Tumor Tissue Analysis Consortium Study

BACKGROUND

Recently, we showed a >60% difference in 5-year survival for patients with tubo-ovarian high-grade serous carcinoma (HGSC) when stratified by a 101-gene mRNA expression prognostic signature. Given the varied patient outcomes, this study aimed to translate prognostic mRNA markers into protein expression assays by immunohistochemistry and validate their survival association in HGSC.

METHODS

Two prognostic genes, FOXJ1 and GMNN, were selected based on high-quality antibodies, correlation with protein expression and variation in immunohistochemical scores in a preliminary cohort (n = 134 and n = 80, respectively). Six thousand four hundred and thirty-four (FOXJ1) and 5470 (GMNN) formalin-fixed, paraffin-embedded ovarian neoplasms (4634 and 4185 HGSC, respectively) represented on tissue microarrays from the Ovarian Tumor Tissue Analysis consortium underwent immunohistochemical staining and scoring, then univariate and multivariate survival analysis.

RESULTS

Consistent with mRNA, FOXJ1 protein expression exhibited a linear, increasing association with improved overall survival in HGSC patients. Women with >50% expression had the most favourable outcomes (HR = 0.78, 95% CI 0.67-0.91, p < 0.0001). GMNN protein expression was not significantly associated with overall HSGC patient survival. However, HGSCs with >35% GMNN expression showed a trend for better outcomes, though this was not significant.

CONCLUSION

We provide foundational evidence for the prognostic value of FOXJ1 in HGSC, validating the prior mRNA-based prognostic association by immunohistochemistry.

Profiling the immune landscape in mucinous ovarian carcinoma

OBJECTIVE

Mucinous ovarian carcinoma (MOC) is a rare histotype of ovarian cancer, with low response rates to standard chemotherapy, and very poor survival for patients diagnosed at advanced stage. There is a limited understanding of the MOC immune landscape, and consequently whether immune checkpoint inhibitors could be considered for a subset of patients.

METHODS

We performed multicolor immunohistochemistry (IHC) and immunofluorescence (IF) on tissue microarrays in a cohort of 126 MOC patients. Cell densities were calculated in the epithelial and stromal components for tumor-associated macrophages (CD68+/PD-L1+, CD68+/PD-L1-), T cells (CD3+/CD8-, CD3+/CD8+), putative T-regulatory cells (Tregs, FOXP3+), B cells (CD20+/CD79A+), plasma cells (CD20-/CD79a+), and PD-L1+ and PD-1+ cells, and compared these values with clinical factors. Univariate and multivariable Cox Proportional Hazards assessed overall survival. Unsupervised k-means clustering identified patient subsets with common patterns of immune cell infiltration.

RESULTS

Mean densities of PD1+ cells, PD-L1- macrophages, CD4+ and CD8+ T cells, and FOXP3+ Tregs were higher in the stroma compared to the epithelium. Tumors from advanced (Stage III/IV) MOC had greater epithelial infiltration of PD-L1- macrophages, and fewer PD-L1+ macrophages compared with Stage I/II cancers (p = 0.004 and p = 0.014 respectively). Patients with high epithelial density of FOXP3+ cells, CD8+/FOXP3+ cells, or PD-L1- macrophages, had poorer survival, and high epithelial CD79a + plasma cells conferred better survival, all upon univariate analysis only. Clustering showed that most MOC (86%) had an immune depleted (cold) phenotype, with only a small proportion (11/76,14%) considered immune inflamed (hot) based on T cell and PD-L1 infiltrates.

CONCLUSION

In summary, MOCs are mostly immunogenically 'cold', suggesting they may have limited response to current immunotherapies.

Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes

INTRODUCTION

Increasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero. For a given m, log-F-penalized logistic regression may be easily implemented using data augmentation and standard software.

METHOD

We propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F-penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation.

RESULTS

We evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F-penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls.

DISCUSSION/CONCLUSION

We have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach.

B-Cell NHL Subtype Risk Associated with Autoimmune Conditions and PRS

BACKGROUND

A previous International Lymphoma Epidemiology (InterLymph) Consortium evaluation of joint associations between five immune gene variants and autoimmune conditions reported interactions between B-cell response-mediated autoimmune conditions and the rs1800629 genotype on risk of B-cell non-Hodgkin lymphoma (NHL) subtypes. Here, we extend that evaluation using NHL subtype-specific polygenic risk scores (PRS) constructed from loci identified in genome-wide association studies of three common B-cell NHL subtypes.

METHODS

In a pooled analysis of NHL cases and controls of Caucasian descent from 14 participating InterLymph studies, we evaluated joint associations between B-cell-mediated autoimmune conditions and tertile (T) of PRS for risk of diffuse large B-cell lymphoma (DLBCL; n = 1,914), follicular lymphoma (n = 1,733), and marginal zone lymphoma (MZL; n = 407), using unconditional logistic regression.

RESULTS

We demonstrated a positive association of DLBCL PRS with DLBCL risk [T2 vs. T1: OR = 1.24; 95% confidence interval (CI), 1.08-1.43; T3 vs. T1: OR = 1.81; 95% CI, 1.59-2.07; P-trend (Ptrend) < 0.0001]. DLBCL risk also increased with increasing PRS tertile among those with an autoimmune condition, being highest for those with a B-cell-mediated autoimmune condition and a T3 PRS [OR = 6.46 vs. no autoimmune condition and a T1 PRS, Ptrend < 0.0001, P-interaction (Pinteraction) = 0.49]. Follicular lymphoma and MZL risk demonstrated no evidence of joint associations or significant Pinteraction.

CONCLUSIONS

Our results suggest that PRS constructed from currently known subtype-specific loci may not necessarily capture biological pathways shared with autoimmune conditions.

IMPACT

Targeted genetic (PRS) screening among population subsets with autoimmune conditions may offer opportunities for identifying those at highest risk for (and early detection from) DLBCL.

Genome-wide interaction analysis of menopausal hormone therapy use and breast cancer risk among 62,370 women

Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values < 5 × 10-8 as genome-wide significant, and p-values < 1 × 10-5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants. None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values < 1 × 105. The strongest evidence was found for rs4674019 (p-value = 2.27 × 10-7), which showed genome-wide significant interaction (p-value = 3.8 × 10-8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen-progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT-breast cancer risk association.

Functional MRI evaluation of cognitive effects of carotid stenosis revascularization

INTRODUCTION

Severe internal carotid stenosis, if left untreated, can pose serious risks for ischemic stroke and cognitive impairments. The effects of revascularization on any aspects of cognition, however, are not well understood, as conflicting results are reported, which have mainly been centered on paper-based cognitive analyses. Here, we summarized and evaluated the publications to date of functional MRI (fMRI) studies that examined the mechanisms of functional brain activation and connectivity as a way to reflect cognitive effects of revascularization on patients with carotid stenosis.

METHODS

A PubMed and Google Scholar (covering the relevant literature until November 1, 2021) search yielded eight original studies of the research line, including seven resting-state and one task-based fMRI reports.

RESULTS

Findings demonstrated treatment-related alterations in fMRI signal intensity and symmetry level, regional fMRI activation pattern, and functional brain network connectivity. The functional brain changes were associated largely with improvement in cognitive function assessed using standard cognitive test scores.

CONCLUSIONS

These findings support the contribution of fMRI to the understanding of brain functional activation and connectivity changes revealing cognitive effects of revascularization in the management of severe carotid stenosis. The review also highlighted the importance of reproducibility through enhancing experimental designs and cognitive task applications with future research for potential clinical translation.

Common variants in breast cancer risk loci predispose to distinct tumor subtypes

BACKGROUND

Genome-wide association studies (GWAS) have identified multiple common breast cancer susceptibility variants. Many of these variants have differential associations by estrogen receptor (ER) status, but how these variants relate with other tumor features and intrinsic molecular subtypes is unclear.

METHODS

Among 106,571 invasive breast cancer cases and 95,762 controls of European ancestry with data on 173 breast cancer variants identified in previous GWAS, we used novel two-stage polytomous logistic regression models to evaluate variants in relation to multiple tumor features (ER, progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and grade) adjusting for each other, and to intrinsic-like subtypes.

RESULTS

Eighty-five of 173 variants were associated with at least one tumor feature (false discovery rate < 5%), most commonly ER and grade, followed by PR and HER2. Models for intrinsic-like subtypes found nearly all of these variants (83 of 85) associated at p < 0.05 with risk for at least one luminal-like subtype, and approximately half (41 of 85) of the variants were associated with risk of at least one non-luminal subtype, including 32 variants associated with triple-negative (TN) disease. Ten variants were associated with risk of all subtypes in different magnitude. Five variants were associated with risk of luminal A-like and TN subtypes in opposite directions.

CONCLUSION

This report demonstrates a high level of complexity in the etiology heterogeneity of breast cancer susceptibility variants and can inform investigations of subtype-specific risk prediction.

Clinical and pathological associations of PTEN expression in ovarian cancer: a multicentre study from the Ovarian Tumour Tissue Analysis Consortium

BACKGROUND

PTEN loss is a putative driver in histotypes of ovarian cancer (high-grade serous (HGSOC), endometrioid (ENOC), clear cell (CCOC), mucinous (MOC), low-grade serous (LGSOC)). We aimed to characterise PTEN expression as a biomarker in epithelial ovarian cancer in a large population-based study.

METHODS

Tumours from 5400 patients from a multicentre observational, prospective cohort study of the Ovarian Tumour Tissue Analysis Consortium were used to evaluate associations between immunohistochemical PTEN patterns and overall survival time, age, stage, grade, residual tumour, CD8+ tumour-infiltrating lymphocytes (TIL) counts, expression of oestrogen receptor (ER), progesterone receptor (PR) and androgen receptor (AR) by means of Cox proportional hazard models and generalised Cochran-Mantel-Haenszel tests.

RESULTS

Downregulation of cytoplasmic PTEN expression was most frequent in ENOC (most frequently in younger patients; p value = 0.0001) and CCOC and was associated with longer overall survival in HGSOC (hazard ratio: 0.78, 95% CI: 0.65-0.94, p value = 0.022). PTEN expression was associated with ER, PR and AR expression (p values: 0.0008, 0.062 and 0.0002, respectively) in HGSOC and with lower CD8 counts in CCOC (p value < 0.0001). Heterogeneous expression of PTEN was more prevalent in advanced HGSOC (p value = 0.019) and associated with higher CD8 counts (p value = 0.0016).

CONCLUSIONS

PTEN loss is a frequent driver in ovarian carcinoma associating distinctly with expression of hormonal receptors and CD8+ TIL counts in HGSOC and CCOC histotypes.

SimRVSequences: an R package to simulate genetic sequence data for pedigrees

SUMMARY

We present the R package SimRVSequences to simulate sequence data for pedigrees. SimRVSequences allows for simulations of large numbers of single-nucleotide variants (SNVs) and scales well with increasing numbers of pedigrees. Users provide a sample of pedigrees and SNV data from a sample of unrelated individuals.

AVAILABILITY AND IMPLEMENTATION

SimRVSequences is publicly-available on CRAN https://cran.r-project.org/web/packages/SimRVSequences/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Dietary patterns in the healthy oldest old in the healthy aging study and the Canadian longitudinal study of aging: a cohort study

Background

Very few people live to eighty-five years and older (the ‘oldest old’), and even fewer live to this age without developing chronic diseases. It is important to understand the relationship, if any, of modifiable factors such as diet on healthy aging. However, there are few studies of diet among healthy oldest old, especially in North American populations. We aimed to characterize dietary patterns among ‘super-seniors’ (SS) within the Canadian Healthy Aging Study.

Methods

122 SS aged 85 years or older and free of cancer, cardiovascular or pulmonary disease, dementia and diabetes were recruited. Comparisons were made to 12,626 participants aged 65–86 in the Canadian Longitudinal Study of Aging who completed the same 36-item food frequency questionnaire that queried consumption over the prior 12 months of nutrients and foods thought to be important for aging. Dietary patterns were identified with principal component analysis. The odds of being a SS were determined for quartiles of each dietary pattern with logistic regression.

Results

Two dietary patterns were identified; a western diet characterized by french fries, red meat, processed meat and a nutrient-rich diet which included fruits, vegetables, whole grains, nuts and seeds among other healthy food choices. Higher scores for both dietary patterns were associated with increased odds of being a SS, however, only the western dietary pattern remained associated with adjustment for covariates (Quartile 4: OR = 3.21, 95% CI 1.91–5.51).

Conclusions

Our finding adds to the limited evidence on dietary intake among the healthiest oldest old but it is unclear whether assocations reflect generational differences between groups or possible contributions to longevity.

Fine-mapping of 150 breast cancer risk regions identifies 191 likely target genes

Genome-wide association studies have identified breast cancer risk variants in over 150 genomic regions, but the mechanisms underlying risk remain largely unknown. These regions were explored by combining association analysis with in silico genomic feature annotations. We defined 205 independent risk-associated signals with the set of credible causal variants in each one. In parallel, we used a Bayesian approach (PAINTOR) that combines genetic association, linkage disequilibrium and enriched genomic features to determine variants with high posterior probabilities of being causal. Potentially causal variants were significantly over-represented in active gene regulatory regions and transcription factor binding sites. We applied our INQUSIT pipeline for prioritizing genes as targets of those potentially causal variants, using gene expression (expression quantitative trait loci), chromatin interaction and functional annotations. Known cancer drivers, transcription factors and genes in the developmental, apoptosis, immune system and DNA integrity checkpoint gene ontology pathways were over-represented among the highest-confidence target genes.

A comprehensive gene-environment interaction analysis in Ovarian Cancer using genome-wide significant common variants

As a follow-up to genome-wide association analysis of common variants associated with ovarian carcinoma (cancer), our study considers seven well-known ovarian cancer risk factors and their interactions with 28 genome-wide significant common genetic variants. The interaction analyses were based on data from 9971 ovarian cancer cases and 15,566 controls from 17 case-control studies. Likelihood ratio and Wald tests for multiplicative interaction and for relative excess risk due to additive interaction were used. The top multiplicative interaction was noted between oral contraceptive pill (OCP) use (ever vs. never) and rs13255292 (p value = 3.48 × 10-4 ). Among women with the TT genotype for this variant, the odds ratio for OCP use was 0.53 (95% CI = 0.46-0.60) compared to 0.71 (95%CI = 0.66-0.77) for women with the CC genotype. When stratified by duration of OCP use, women with 1-5 years of OCP use exhibited differential protective benefit across genotypes. However, no interaction on either the multiplicative or additive scale was found to be statistically significant after multiple testing correction. The results suggest that OCP use may offer increased benefit for women who are carriers of the T allele in rs13255292. On the other hand, for women carrying the C allele in this variant, longer (5+ years) use of OCP may reduce the impact of carrying the risk allele of this SNP. Replication of this finding is needed. The study presents a comprehensive analytic framework for conducting gene-environment analysis in ovarian cancer.

Polygenic Risk Scores for Prediction of Breast Cancer and Breast Cancer Subtypes

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.

rs495139 in the TYMS-ENOSF1 Region and Risk of Ovarian Carcinoma of Mucinous Histology

Thymidylate synthase (TYMS) is a crucial enzyme for DNA synthesis. TYMS expression is regulated by its antisense mRNA, ENOSF1. Disrupted regulation may promote uncontrolled DNA synthesis and tumor growth. We sought to replicate our previously reported association between rs495139 in the TYMS-ENOSF1 3' gene region and increased risk of mucinous ovarian carcinoma (MOC) in an independent sample. Genotypes from 24,351 controls to 15,000 women with invasive OC, including 665 MOC, were available. We estimated per-allele odds ratios (OR) and 95% confidence intervals (CI) using unconditional logistic regression, and meta-analysis when combining these data with our previous report. The association between rs495139 and MOC was not significant in the independent sample (OR = 1.09; 95% CI = 0.97⁻1.22; p = 0.15; N = 665 cases). Meta-analysis suggested a weak association (OR = 1.13; 95% CI = 1.03⁻1.24; p = 0.01; N = 1019 cases). No significant association with risk of other OC histologic types was observed (p = 0.05 for tumor heterogeneity). In expression quantitative trait locus (eQTL) analysis, the rs495139 allele was positively associated with ENOSF1 mRNA expression in normal tissues of the gastrointestinal system, particularly esophageal mucosa (r = 0.51, p = 1.7 × 10-28), and nonsignificantly in five MOC tumors. The association results, along with inconclusive tumor eQTL findings, suggest that a true effect of rs495139 might be small.

HLA Class I and II Diversity Contributes to the Etiologic Heterogeneity of Non-Hodgkin Lymphoma Subtypes

A growing number of loci within the human leukocyte antigen (HLA) region have been implicated in non-Hodgkin lymphoma (NHL) etiology. Here, we test a complementary hypothesis of "heterozygote advantage" regarding the role of HLA and NHL, whereby HLA diversity is beneficial and homozygous HLA loci are associated with increased disease risk. HLA alleles at class I and II loci were imputed from genome-wide association studies (GWAS) using SNP2HLA for 3,617 diffuse large B-cell lymphomas (DLBCL), 2,686 follicular lymphomas (FL), 2,878 chronic lymphocytic leukemia/small lymphocytic lymphomas (CLL/SLL), 741 marginal zone lymphomas (MZL), and 8,753 controls of European descent. Both DLBCL and MZL risk were elevated with homozygosity at class I HLA-B and -C loci (OR DLBCL = 1.31, 95% CI = 1.06-1.60; OR MZL = 1.45, 95% CI = 1.12-1.89) and class II HLA-DRB1 locus (OR DLBCL = 2.10, 95% CI = 1.24-3.55; OR MZL = 2.10, 95% CI = 0.99-4.45). Increased FL risk was observed with the overall increase in number of homozygous HLA class II loci (P trend < 0.0001, FDR = 0.0005). These results support a role for HLA zygosity in NHL etiology and suggests that distinct immune pathways may underly the etiology of the different NHL subtypes.Significance: HLA gene diversity reduces risk for non-Hodgkin lymphoma. Cancer Res; 78(14); 4086-96. ©2018 AACR.

Variants in genes encoding small GTPases and association with epithelial ovarian cancer susceptibility

Epithelial ovarian cancer (EOC) is the fifth leading cause of cancer mortality in American women. Normal ovarian physiology is intricately connected to small GTP binding proteins of the Ras superfamily (Ras, Rho, Rab, Arf, and Ran) which govern processes such as signal transduction, cell proliferation, cell motility, and vesicle transport. We hypothesized that common germline variation in genes encoding small GTPases is associated with EOC risk. We investigated 322 variants in 88 small GTPase genes in germline DNA of 18,736 EOC patients and 26,138 controls of European ancestry using a custom genotype array and logistic regression fitting log-additive models. Functional annotation was used to identify biofeatures and expression quantitative trait loci that intersect with risk variants. One variant, ARHGEF10L (Rho guanine nucleotide exchange factor 10 like) rs2256787, was associated with increased endometrioid EOC risk (OR = 1.33, p = 4.46 x 10-6). Other variants of interest included another in ARHGEF10L, rs10788679, which was associated with invasive serous EOC risk (OR = 1.07, p = 0.00026) and two variants in AKAP6 (A-kinase anchoring protein 6) which were associated with risk of invasive EOC (rs1955513, OR = 0.90, p = 0.00033; rs927062, OR = 0.94, p = 0.00059). Functional annotation revealed that the two ARHGEF10L variants were located in super-enhancer regions and that AKAP6 rs927062 was associated with expression of GTPase gene ARHGAP5 (Rho GTPase activating protein 5). Inherited variants in ARHGEF10L and AKAP6, with potential transcriptional regulatory function and association with EOC risk, warrant investigation in independent EOC study populations.

Adult height is associated with increased risk of ovarian cancer: a Mendelian randomisation study

BACKGROUND

Observational studies suggest greater height is associated with increased ovarian cancer risk, but cannot exclude bias and/or confounding as explanations for this. Mendelian randomisation (MR) can provide evidence which may be less prone to bias.

METHODS

We pooled data from 39 Ovarian Cancer Association Consortium studies (16,395 cases; 23,003 controls). We applied two-stage predictor-substitution MR, using a weighted genetic risk score combining 609 single-nucleotide polymorphisms. Study-specific odds ratios (OR) and 95% confidence intervals (CI) for the association between genetically predicted height and risk were pooled using random-effects meta-analysis.

RESULTS

Greater genetically predicted height was associated with increased ovarian cancer risk overall (pooled-OR (pOR) = 1.06; 95% CI: 1.01-1.11 per 5 cm increase in height), and separately for invasive (pOR = 1.06; 95% CI: 1.01-1.11) and borderline (pOR = 1.15; 95% CI: 1.02-1.29) tumours.

CONCLUSIONS

Women with a genetic propensity to being taller have increased risk of ovarian cancer. This suggests genes influencing height are involved in pathways promoting ovarian carcinogenesis.

Racial/ethnic differences in the epidemiology of ovarian cancer: a pooled analysis of 12 case-control studies

Background

Ovarian cancer incidence differs substantially by race/ethnicity, but the reasons for this are not well understood. Data were pooled from the African American Cancer Epidemiology Study (AACES) and 11 case-control studies in the Ovarian Cancer Association Consortium (OCAC) to examine racial/ethnic differences in epidemiological characteristics with suspected involvement in epithelial ovarian cancer (EOC) aetiology.

Methods

We used multivariable logistic regression to estimate associations for 17 reproductive, hormonal and lifestyle characteristics and EOC risk by race/ethnicity among 10 924 women with invasive EOC (8918 Non-Hispanic Whites, 433 Hispanics, 911 Blacks, 662 Asian/Pacific Islanders) and 16 150 controls (13 619 Non-Hispanic Whites, 533 Hispanics, 1233 Blacks, 765 Asian/Pacific Islanders). Likelihood ratio tests were used to evaluate heterogeneity in the risk factor associations by race/ethnicity.

Results

We observed statistically significant racial/ethnic heterogeneity for hysterectomy and EOC risk (P = 0.008), where the largest odds ratio (OR) was observed in Black women [OR = 1.64, 95% confidence interval (CI) = 1.34-2.02] compared with other racial/ethnic groups. Although not statistically significant, the associations for parity, first-degree family history of ovarian or breast cancer, and endometriosis varied by race/ethnicity. Asian/Pacific Islanders had the greatest magnitude of association for parity (≥3 births: OR = 0.38, 95% CI = 0.28-0.54), and Black women had the largest ORs for family history (OR = 1.77, 95% CI = 1.42-2.21) and endometriosis (OR = 2.42, 95% CI = 1.65-3.55).

Conclusions

Although racial/ethnic heterogeneity was observed for hysterectomy, our findings support the validity of EOC risk factors across all racial/ethnic groups, and further suggest that any racial/ethnic population with a higher prevalence of a modifiable risk factor should be targeted to disseminate information about prevention.

Identification of ten variants associated with risk of estrogen-receptor-negative breast cancer

Most common breast cancer susceptibility variants have been identified through genome-wide association studies (GWAS) of predominantly estrogen receptor (ER)-positive disease. We conducted a GWAS using 21,468 ER-negative cases and 100,594 controls combined with 18,908 BRCA1 mutation carriers (9,414 with breast cancer), all of European origin. We identified independent associations at P < 5 × 10-8 with ten variants at nine new loci. At P < 0.05, we replicated associations with 10 of 11 variants previously reported in ER-negative disease or BRCA1 mutation carrier GWAS and observed consistent associations with ER-negative disease for 105 susceptibility variants identified by other studies. These 125 variants explain approximately 16% of the familial risk of this breast cancer subtype. There was high genetic correlation (0.72) between risk of ER-negative breast cancer and breast cancer risk for BRCA1 mutation carriers. These findings may lead to improved risk prediction and inform further fine-mapping and functional work to better understand the biological basis of ER-negative breast cancer.

Effect of age on chronic inflammation and responsiveness to bacterial and viral challenges

To identify reliable biomarkers of age-related changes in chronic inflammation and responsiveness to bacterial and viral challenges, we evaluated endogenous and ex vivo stimulated levels of 18 inflammatory markers, using whole blood collected in EDTA and sodium heparin tubes from 41 healthy volunteers, i.e., 11 men + 10 women aged 20–35 and 10 men + 10 women aged 50–77. These studies revealed significant differences in the levels of inflammatory markers when blood was collected in EDTA versus sodium heparin and age related differences in these biomarkers were confirmed with blood collected in EDTA from 120 healthy volunteers in 3 age categories, ie, 20 men + 20 women, aged 20–35, 36–49 and 50–77. Studies with unstimulated blood samples, to measure levels of chronic inflammation, revealed a significant increase with age in IL-12p70, CRP and PGE₂, consistent with the concept of “inflammaging”, and a decrease in G-CSF in both men and women. Interestingly, in response to E. coli stimulation, PGE₂ levels were markedly reduced in the 50–77 year old cohort while they were increased following Herpes Simplex virus-1 (HSV-1) stimulation, along with IL-8. In addition, unlike E. coli, HSV-1 potently stimulated IFNα production, but levels were dramatically reduced in the older cohort, consistent with a reduced ability to generate an anti-viral response. We also found platelets and CD8⁺ T cells were reduced with age while CD4⁺ T cells were significantly increased, resulting in a substantially higher CD4/CD8 ratio in the older cohort. Surprisingly, however, we found that the older cohort exhibited more T cell proliferation and IFNγ production in response to anti-CD3+anti-CD28 stimulation. Importantly, there was considerable person-to-person variation in these inflammatory markers in all age groups, making possible comparisons between a person’s “inflammage” and chronological age. These assays should help to identify individuals at high risk of autoimmune disorders and cancer.

Association analysis identifies 65 new breast cancer risk loci

Breast cancer risk is influenced by rare coding variants in susceptibility genes, such as BRCA1, and many common, mostly non-coding variants. However, much of the genetic contribution to breast cancer risk remains unknown. Here we report the results of a genome-wide association study of breast cancer in 122,977 cases and 105,974 controls of European ancestry and 14,068 cases and 13,104 controls of East Asian ancestry. We identified 65 new loci that are associated with overall breast cancer risk at P < 5 × 10-8. The majority of credible risk single-nucleotide polymorphisms in these loci fall in distal regulatory elements, and by integrating in silico data to predict target genes in breast cells at each locus, we demonstrate a strong overlap between candidate target genes and somatic driver genes in breast tumours. We also find that heritability of breast cancer due to all single-nucleotide polymorphisms in regulatory features was 2-5-fold enriched relative to the genome-wide average, with strong enrichment for particular transcription factor binding sites. These results provide further insight into genetic susceptibility to breast cancer and will improve the use of genetic risk scores for individualized screening and prevention.

Identification of 12 new susceptibility loci for different histotypes of epithelial ovarian cancer

To identify common alleles associated with different histotypes of epithelial ovarian cancer (EOC), we pooled data from multiple genome-wide genotyping projects totaling 25,509 EOC cases and 40,941 controls. We identified nine new susceptibility loci for different EOC histotypes: six for serous EOC histotypes (3q28, 4q32.3, 8q21.11, 10q24.33, 18q11.2 and 22q12.1), two for mucinous EOC (3q22.3 and 9q31.1) and one for endometrioid EOC (5q12.3). We then performed meta-analysis on the results for high-grade serous ovarian cancer with the results from analysis of 31,448 BRCA1 and BRCA2 mutation carriers, including 3,887 mutation carriers with EOC. This identified three additional susceptibility loci at 2q13, 8q24.1 and 12q24.31. Integrated analyses of genes and regulatory biofeatures at each locus predicted candidate susceptibility genes, including OBFC1, a new candidate susceptibility gene for low-grade and borderline serous EOC.

Nonrandom occurrence of lymphoid cancer types in 140 families

We studied 140 families with two or more lymphoid cancers, including non-Hodgkin lymphoma (NHL), Hodgkin lymphoma (HL), chronic lymphocytic leukemia (CLL), and multiple myeloma (MM), for deviation from the population age of onset and lymphoid cancer co-occurrence patterns. Median familial NHL, HL, CLL and MM ages of onset are substantially earlier than comparable population data. NHL, HL and CLL (but not MM) also show earlier age of onset in later generations, known as anticipation. The co-occurrence of lymphoid cancers is significantly different from that expected based on population frequencies (p < .0001), and the pattern differs more in families with more affected members (p < .0001), suggesting specific lymphoid cancer combinations have a shared genetic basis. These families provide evidence for inherited factors that increase the risk of multiple lymphoid cancers. This study was approved by the BC Cancer Agency - University of British Columbia Clinical Research Ethics Board.

A targeted genetic association study of epithelial ovarian cancer susceptibility

BACKGROUND

Genome-wide association studies have identified several common susceptibility alleles for epithelial ovarian cancer (EOC). To further understand EOC susceptibility, we examined previously ungenotyped candidate variants, including uncommon variants and those residing within known susceptibility loci.

RESULTS

At nine of eleven previously published EOC susceptibility regions (2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13), novel variants were identified that were more strongly associated with risk than previously reported variants. Beyond known susceptibility regions, no variants were found to be associated with EOC risk at genome-wide statistical significance (p <5x10(-8)), nor were any significant after Bonferroni correction for 17,000 variants (p< 3x10-6).

METHODS

A customized genotyping array was used to assess over 17,000 variants in coding, non-coding, regulatory, and known susceptibility regions in 4,973 EOC cases and 5,640 controls from 13 independent studies. Susceptibility for EOC overall and for select histotypes was evaluated using logistic regression adjusted for age, study site, and population substructure.

CONCLUSION

Given the novel variants identified within the 2q31, 3q25, 5p15, 8q21, 8q24, 10p12, 17q12, 17q21.31, and 19p13 regions, larger follow-up genotyping studies, using imputation where necessary, are needed for fine-mapping and confirmation of low frequency variants that fall below statistical significance.

Combined oral contraceptive use before the first birth and epithelial ovarian cancer risk

Background:

Combined oral contraceptive (COC) use reduces epithelial ovarian cancer (EOC) risk. However, little is known about risk with COC use before the first full-term pregnancy (FFTP).

Methods:

This Canadian population-based case–control study (2001–2012) included 854 invasive cases/2139 controls aged ⩾40 years who were parous and had information on COC use. We estimated odds ratios (aORs) and 95% confidence intervals (CI) adjusted for study site, age, parity, breastfeeding, age at FFTP, familial breast/ovarian cancer, tubal ligation, and body mass.

Results:

Among parous women, per year of COC use exclusively before the FFTP was associated with a 9% risk reduction (95% CI=0.86–0.96). Results were similar for high-grade serous and endometrioid/clear cell EOC. In contrast, per year of use exclusively after the FFTP was not associated with risk (aOR=0.98, 95% CI=0.95–1.02).

Conclusions:

Combined oral contraceptive use before the FFTP may provide a risk reduction that remains for many years, informing possible prevention strategies.

PALB2, CHEK2 and ATM rare variants and cancer risk: data from COGS

BACKGROUND

The rarity of mutations in PALB2, CHEK2 and ATM make it difficult to estimate precisely associated cancer risks. Population-based family studies have provided evidence that at least some of these mutations are associated with breast cancer risk as high as those associated with rare BRCA2 mutations. We aimed to estimate the relative risks associated with specific rare variants in PALB2, CHEK2 and ATM via a multicentre case-control study.

METHODS

We genotyped 10 rare mutations using the custom iCOGS array: PALB2 c.1592delT, c.2816T>G and c.3113G>A, CHEK2 c.349A>G, c.538C>T, c.715G>A, c.1036C>T, c.1312G>T, and c.1343T>G and ATM c.7271T>G. We assessed associations with breast cancer risk (42 671 cases and 42 164 controls), as well as prostate (22 301 cases and 22 320 controls) and ovarian (14 542 cases and 23 491 controls) cancer risk, for each variant.

RESULTS

For European women, strong evidence of association with breast cancer risk was observed for PALB2 c.1592delT OR 3.44 (95% CI 1.39 to 8.52, p=7.1×10-5), PALB2 c.3113G>A OR 4.21 (95% CI 1.84 to 9.60, p=6.9×10-8) and ATM c.7271T>G OR 11.0 (95% CI 1.42 to 85.7, p=0.0012). We also found evidence of association with breast cancer risk for three variants in CHEK2, c.349A>G OR 2.26 (95% CI 1.29 to 3.95), c.1036C>T OR 5.06 (95% CI 1.09 to 23.5) and c.538C>T OR 1.33 (95% CI 1.05 to 1.67) (p≤0.017). Evidence for prostate cancer risk was observed for CHEK2 c.1343T>G OR 3.03 (95% CI 1.53 to 6.03, p=0.0006) for African men and CHEK2 c.1312G>T OR 2.21 (95% CI 1.06 to 4.63, p=0.030) for European men. No evidence of association with ovarian cancer was found for any of these variants.

CONCLUSIONS

This report adds to accumulating evidence that at least some variants in these genes are associated with an increased risk of breast cancer that is clinically important.

Inherited variants affecting RNA editing may contribute to ovarian cancer susceptibility: results from a large-scale collaboration

RNA editing in mammals is a form of post-transcriptional modification in which adenosine is converted to inosine by the adenosine deaminases acting on RNA (ADAR) family of enzymes. Based on evidence of altered ADAR expression in epithelial ovarian cancers (EOC), we hypothesized that single nucleotide polymorphisms (SNPs) in ADAR genes modify EOC susceptibility, potentially by altering ovarian tissue gene expression. Using directly genotyped and imputed data from 10,891 invasive EOC cases and 21,693 controls, we evaluated the associations of 5,303 SNPs in ADAD1, ADAR, ADAR2, ADAR3, and SND1. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI), with adjustment for European ancestry. We conducted gene-level analyses using the Admixture Maximum Likelihood (AML) test and the Sequence-Kernel Association test for common and rare variants (SKAT-CR). Association analysis revealed top risk-associated SNP rs77027562 (OR (95% CI)= 1.39 (1.17-1.64), P=1.0x10-4) in ADAR3 and rs185455523 in SND1 (OR (95% CI)= 0.68 (0.56-0.83), P=2.0x10-4). When restricting to serous histology (n=6,500), the magnitude of association strengthened for rs185455523 (OR=0.60, P=1.0x10-4). Gene-level analyses revealed that variation in ADAR was associated (P<0.05) with EOC susceptibility, with PAML=0.022 and PSKAT-CR=0.020. Expression quantitative trait locus analysis in EOC tissue revealed significant associations (P<0.05) with ADAR expression for several SNPs in ADAR, including rs1127313 (G/A), a SNP in the 3' untranslated region. In summary, germline variation involving RNA editing genes may influence EOC susceptibility, warranting further investigation of inherited and acquired alterations affecting RNA editing.

Tea, coffee, and caffeinated beverage consumption and risk of epithelial ovarian cancers

BACKGROUND

The risk for epithelial ovarian cancer associated with the consumption of caffeinated beverages (tea, coffee, and soft drinks) and green tea is inconclusive. However, few studies have investigated the type of caffeinated beverage or the type of tea.

OBJECTIVE

We assessed consumption of tea (black/caffeinated tea and green tea separately), coffee, and caffeinated soft drinks, as well as level of consumption, and the risk for epithelial ovarian cancer and its histotypes.

STUDY DESIGN

This study was conducted within a population-based case-control study in Alberta and British Columbia, Canada from 2001 to 2012. After restricting to cases of epithelial invasive cancers and controls aged 40-79 years who completed an interview that included coffee, soft drink, and tea consumption (ascertained starting in 2005 in British Columbia and 2008 in Alberta), there were a total of 524 cases and 1587 controls. Those that did not meet the threshold for beverage consumption (at least once per month for 6 months or more) were classified as non-drinkers. Adult lifetime cumulative consumption (cup-years=cups/day*years) was calculated. Unconditional logistic regression was used to estimate adjusted odds ratios (aOR) and 95% confidence intervals (CI) to describe the association between the relevant drink consumption and risk.

RESULTS

No excess risk was seen for coffee or caffeinated soft drinks. Similarly, any tea consumption was not associated with risk, but when stratified by the type of tea, there was an increase in risk in black tea only drinkers (aOR=1.56; 95% CI:1.07-2.28 for >40 cup-years), but no excess risk for the exclusive green tea drinkers. Similar findings were observed for post-menopausal women. The association for black tea only consumption was mainly seen in the endometrioid histotype (aOR=3.19; 95% CI: 1.32-7.69).

CONCLUSION

Black tea consumption may be associated with an increased risk epithelial ovarian carcinoma. The excess risk is seen only in the endometrioid histotype but not in serous or clear cell. Further studies are required to confirm these findings and identify the constituents in black tea that may increase the risk.

Assessment of variation in immunosuppressive pathway genes reveals TGFBR2 to be associated with risk of clear cell ovarian cancer

BACKGROUND

Regulatory T (Treg) cells, a subset of CD4+ T lymphocytes, are mediators of immunosuppression in cancer, and, thus, variants in genes encoding Treg cell immune molecules could be associated with ovarian cancer.

METHODS

In a population of 15,596 epithelial ovarian cancer (EOC) cases and 23,236 controls, we measured genetic associations of 1,351 SNPs in Treg cell pathway genes with odds of ovarian cancer and tested pathway and gene-level associations, overall and by histotype, for the 25 genes, using the admixture likelihood (AML) method. The most significant single SNP associations were tested for correlation with expression levels in 44 ovarian cancer patients.

RESULTS

The most significant global associations for all genes in the pathway were seen in endometrioid ( p = 0.082) and clear cell ( p = 0.083), with the most significant gene level association seen with TGFBR2 ( p = 0.001) and clear cell EOC. Gene associations with histotypes at p < 0.05 included: IL12 ( p = 0.005 and p = 0.008, serous and high-grade serous, respectively), IL8RA ( p = 0.035, endometrioid and mucinous), LGALS1 ( p = 0.03, mucinous), STAT5B ( p = 0.022, clear cell), TGFBR1 ( p = 0.021 endometrioid) and TGFBR2 ( p = 0.017 and p = 0.025, endometrioid and mucinous, respectively).

CONCLUSIONS

Common inherited gene variation in Treg cell pathways shows some evidence of germline genetic contribution to odds of EOC that varies by histologic subtype and may be associated with mRNA expression of immune-complex receptor in EOC patients.

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.

Exome genotyping arrays to identify rare and low frequency variants associated with epithelial ovarian cancer risk

Rare and low frequency variants are not well covered in most germline genotyping arrays and are understudied in relation to epithelial ovarian cancer (EOC) risk. To address this gap, we used genotyping arrays targeting rarer protein-coding variation in 8,165 EOC cases and 11,619 controls from the international Ovarian Cancer Association Consortium (OCAC). Pooled association analyses were conducted at the variant and gene level for 98,543 variants directly genotyped through two exome genotyping projects. Only common variants that represent or are in strong linkage disequilibrium (LD) with previously-identified signals at established loci reached traditional thresholds for exome-wide significance (P < 5.0 × 10 -  7). One of the most significant signals (Pall histologies = 1.01 × 10 -  13;Pserous = 3.54 × 10 -  14) occurred at 3q25.31 for rs62273959, a missense variant mapping to the LEKR1 gene that is in LD (r2 = 0.90) with a previously identified 'best hit' (rs7651446) mapping to an intron of TIPARP. Suggestive associations (5.0 × 10 -  5 > P≥5.0 ×10 -  7) were detected for rare and low-frequency variants at 16 novel loci. Four rare missense variants were identified (ACTBL2 rs73757391 (5q11.2), BTD rs200337373 (3p25.1), KRT13 rs150321809 (17q21.2) and MC2R rs104894658 (18p11.21)), but only MC2R rs104894668 had a large effect size (OR = 9.66). Genes most strongly associated with EOC risk included ACTBL2 (PAML = 3.23 × 10 -  5; PSKAT-o = 9.23 × 10 -  4) and KRT13 (PAML = 1.67 × 10 -  4; PSKAT-o = 1.07 × 10 -  5), reaffirming variant-level analysis. In summary, this large study identified several rare and low-frequency variants and genes that may contribute to EOC susceptibility, albeit with possible small effects. Future studies that integrate epidemiology, sequencing, and functional assays are needed to further unravel the unexplained heritability and biology of this disease.

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.

Polymorphisms of Insulin-Like Growth Factor 1 Pathway Genes and Breast Cancer Risk

Genetic variants of insulin-like growth factor 1 (IGF1) pathway genes have been shown to be associated with breast density and IGF1 levels and, therefore, may also influence breast cancer risk via pro-survival signaling cascades. The aim of this study was to investigate associations between IGF1 pathway single nucleotide polymorphisms (SNPs) and breast cancer risk among European and East Asian women, and potential interactions with menopausal status and breast tumor subtype. Stratified analyses of 1,037 cases and 1,050 controls from a population-based case–control study were conducted to assess associations with breast cancer for 22 SNPs across 5 IGF1 pathway genes in European and East Asian women. Odds ratios were calculated using logistic regression in additive genetic models. Polytomous logistic regression was used to assess heterogeneity by breast tumor subtype. Two SNPs of the IGF1 gene (rs1019731 and rs12821878) were associated with breast cancer risk among European women. Four highly linked IGF1 SNPs (rs2288378, rs17727841, rs7136446, and rs7956547) were modified by menopausal status among East Asian women only and associated with postmenopausal breast cancers. The association between rs2288378 and breast cancer risk was also modified by breast tumor subtype among East Asian women. Several IGF1 polymorphisms were found to be associated with breast cancer risk and some of these associations were modified by menopausal status or breast tumor subtype. Such interactions should be considered when assessing the role of these variants in breast cancer etiology.

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.

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.

DNA repair variants and breast cancer risk

A functional DNA repair system has been identified as important in the prevention of tumour development. Previous studies have hypothesized that common polymorphisms in DNA repair genes could play a role in breast cancer risk and also identified the potential for interactions between these polymorphisms and established breast cancer risk factors such as physical activity. Associations with breast cancer risk for 99 single nucleotide polymorphisms (SNPs) from genes in ten DNA repair pathways were examined in a case-control study including both Europeans (644 cases, 809 controls) and East Asians (299 cases, 160 controls). Odds ratios in both additive and dominant genetic models were calculated separately for participants of European and East Asian ancestry using multivariate logistic regression. The impact of multiple comparisons was assessed by correcting for the false discovery rate within each DNA repair pathway. Interactions between several breast cancer risk factors and DNA repair SNPs were also evaluated. One SNP (rs3213282) in the gene XRCC1 was associated with an increased risk of breast cancer in the dominant model of inheritance following adjustment for the false discovery rate (P < 0.05), although no associations were observed for other DNA repair SNPs. Interactions of six SNPs in multiple DNA repair pathways with physical activity were evident prior to correction for FDR, following which there was support for only one of the interaction terms (P < 0.05). No consistent associations between variants in DNA repair genes and breast cancer risk or their modification by breast cancer risk factors were observed.

Adult lifetime alcohol consumption and invasive epithelial ovarian cancer risk in a population-based case-control study

OBJECTIVE

Meta-analyses report a null association between recent alcohol consumption and ovarian cancer risk. However, because few studies investigated different types of alcohol over adult ages, we investigated adult lifetime and type (beer, wine, spirits) of consumption and risk.

METHODS

Consumption after age 20years was ascertained in 1144 invasive epithelial ovarian cancer cases and 2513 controls in a population-based case-control study (Alberta and British Columbia, Canada, 2001-2012). Non-drinkers consumed any types of alcohol <12 times per year on average. Logistic regression was use to estimate adjusted odds ratios [aOR] and 95% confidence intervals [CIs].

RESULTS

Wine consumption was associated with a risk reduction (aOR=0.67, 95% CI: 0.50-0.88) relative to non-drinkers, but not beer (aOR=1.06, 95% CI: 0.71-1.58) or spirits (aOR=0.98, 95% CI: 0.69-1.39). The reduced risk was stronger for exclusive red wine drinkers (aOR=0.44, 95% CI: 0.19-0.92) than white wine drinkers (aOR=0.79, 95% CI: 0.46-1.34), although most women drank both types of wine. Risk decreased with increasing cumulative consumption of any wine (P-trend<0.05) and was evident for the serous histotype. Wine consumption initiated prior to age 50 was associated with a risk reduction (e.g., at 40-49years, aOR=0.58, 95% CI: 0.42-0.78), but not drinking initiated after 50years of age. For any type, level, or age at initiation of alcohol consumption, we found no increased risks.

CONCLUSIONS

For the moderate consumption in this study, higher levels of wine consumption were generally associated with risk reductions; reductions may be stronger for red wine. Our results suggest that alcohol consumption that is guideline concordant will not increase epithelial ovarian cancer risk.

Genetic variation in vitamin D-related genes and risk of breast cancer among women of European and East Asian descent

Studies of vitamin D-related genetic variants and breast cancer have been inconsistent. This study aimed to investigate associations of vitamin D-related polymorphisms and breast cancer risk among European and East Asian women and potential interactions with menopausal status and breast tumour subtypes. Data from a case-control study of breast cancer (1037 cases and 1050 controls) were used to assess relationships between 21 polymorphisms in two vitamin D-related genes (GC and VDR) and breast cancer risk. Odds ratios were calculated in stratified analyses of European and East Asian women, using logistic regression in an additive genetic model. An interaction term was used to explore modification by menopausal status. Polytomous regression was used to assess heterogeneity by breast tumour subtype. False discovery rate adjustments were conducted to account for multiple testing. No association was observed between GC or VDR polymorphisms and breast cancer risk. Modification of these relationships by menopausal status was observed for select polymorphisms in both Europeans (VDR rs4328262 and rs11168292) and East Asians (GC rs7041 and VDR rs11168287). Heterogeneity by tumour subtype was seen for three VDR polymorphisms (rs1544410, rs7967152 and rs2239186) among Europeans, in which associations with ER-/PR-/HER2+ tumours, but not with other subtypes, were observed. In conclusion, associations between vitamin D-related genetic variants and breast cancer were not observed overall, although the relationships between vitamin D pathway polymorphisms and breast cancer may be modified by menopausal status and breast tumour subtype.

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.