Genetic insights into biological mechanisms governing human ovarian ageing.

Reproductive longevity is essential for fertility and influences healthy ageing in women1,2, but insights into its underlying biological mechanisms and treatments to preserve it are limited. Here we identify 290 genetic determinants of ovarian ageing, assessed using normal variation in age at natural menopause (ANM) in about 200,000 women of European ancestry. These common alleles were associated with clinical extremes of ANM; women in the top 1% of genetic susceptibility have an equivalent risk of premature ovarian insufficiency to those carrying monogenic FMR1 premutations3. The identified loci implicate a broad range of DNA damage response (DDR) processes and include loss-of-function variants in key DDR-associated genes. Integration with experimental models demonstrates that these DDR processes act across the life-course to shape the ovarian reserve and its rate of depletion. Furthermore, we demonstrate that experimental manipulation of DDR pathways highlighted by human genetics increases fertility and extends reproductive life in mice. Causal inference analyses using the identified genetic variants indicate that extending reproductive life in women improves bone health and reduces risk of type 2 diabetes, but increases the risk of hormone-sensitive cancers. These findings provide insight into the mechanisms that govern ovarian ageing, when they act, and how they might be targeted by therapeutic approaches to extend fertility and prevent disease.

The impact of coding germline variants on contralateral breast cancer risk and survival.

Evidence linking coding germline variants in breast cancer (BC)-susceptibility genes other than BRCA1, BRCA2, and CHEK2 with contralateral breast cancer (CBC) risk and breast cancer-specific survival (BCSS) is scarce. The aim of this study was to assess the association of protein-truncating variants (PTVs) and rare missense variants (MSVs) in nine known (ATM, BARD1, BRCA1, BRCA2, CHEK2, PALB2, RAD51C, RAD51D, and TP53) and 25 suspected BC-susceptibility genes with CBC risk and BCSS. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated with Cox regression models. Analyses included 34,401 women of European ancestry diagnosed with BC, including 676 CBCs and 3,449 BC deaths; the median follow-up was 10.9 years. Subtype analyses were based on estrogen receptor (ER) status of the first BC. Combined PTVs and pathogenic/likely pathogenic MSVs in BRCA1, BRCA2, and TP53 and PTVs in CHEK2 and PALB2 were associated with increased CBC risk [HRs (95% CIs): 2.88 (1.70-4.87), 2.31 (1.39-3.85), 8.29 (2.53-27.21), 2.25 (1.55-3.27), and 2.67 (1.33-5.35), respectively]. The strongest evidence of association with BCSS was for PTVs and pathogenic/likely pathogenic MSVs in BRCA2 (ER-positive BC) and TP53 and PTVs in CHEK2 [HRs (95% CIs): 1.53 (1.13-2.07), 2.08 (0.95-4.57), and 1.39 (1.13-1.72), respectively, after adjusting for tumor characteristics and treatment]. HRs were essentially unchanged when censoring for CBC, suggesting that these associations are not completely explained by increased CBC risk, tumor characteristics, or treatment. There was limited evidence of associations of PTVs and/or rare MSVs with CBC risk or BCSS for the 25 suspected BC genes. The CBC findings are relevant to treatment decisions, follow-up, and screening after BC diagnosis.

The sequence kernel association test for multicategorical outcomes.

Disease heterogeneity is ubiquitous in biomedical and clinical studies. In genetic studies, researchers are increasingly interested in understanding the distinct genetic underpinning of subtypes of diseases. However, existing set-based analysis methods for genome-wide association studies are either inadequate or inefficient to handle such multicategorical outcomes. In this paper, we proposed a novel set-based association analysis method, sequence kernel association test (SKAT)-MC, the sequence kernel association test for multicategorical outcomes (nominal or ordinal), which jointly evaluates the relationship between a set of variants (common and rare) and disease subtypes. Through comprehensive simulation studies, we showed that SKAT-MC effectively preserves the nominal type I error rate while substantially increases the statistical power compared to existing methods under various scenarios. We applied SKAT-MC to the Polish breast cancer study (PBCS), and identified gene FGFR2 was significantly associated with estrogen receptor (ER)+ and ER- breast cancer subtypes. We also investigated educational attainment using UK Biobank data ( N = 127 , 127 $N=127,127$ ) with SKAT-MC, and identified 21 significant genes in the genome. Consequently, SKAT-MC is a powerful and efficient analysis tool for genetic association studies with multicategorical outcomes. A freely distributed R package SKAT-MC can be accessed at https://github.com/Zhiwen-Owen-Jiang/SKATMC.

Polygenic risk scores for prediction of breast cancer risk in women of African ancestry: a cross-ancestry approach.

Polygenic risk scores (PRSs) are useful for predicting breast cancer risk, but the prediction accuracy of existing PRSs in women of African ancestry (AA) remains relatively low. We aim to develop optimal PRSs for the prediction of overall and estrogen receptor (ER) subtype-specific breast cancer risk in AA women. The AA dataset comprised 9235 cases and 10 184 controls from four genome-wide association study (GWAS) consortia and a GWAS study in Ghana. We randomly divided samples into training and validation sets. We built PRSs using individual-level AA data by a forward stepwise logistic regression and then developed joint PRSs that combined (1) the PRSs built in the AA training dataset and (2) a 313-variant PRS previously developed in women of European ancestry. PRSs were evaluated in the AA validation set. For overall breast cancer, the odds ratio per standard deviation of the joint PRS in the validation set was 1.34 [95% confidence interval (CI): 1.27-1.42] with the area under receiver operating characteristic curve (AUC) of 0.581. Compared with women with average risk (40th-60th PRS percentile), women in the top decile of the PRS had a 1.98-fold increased risk (95% CI: 1.63-2.39). For PRSs of ER-positive and ER-negative breast cancer, the AUCs were 0.608 and 0.576, respectively. Compared with existing methods, the proposed joint PRSs can improve prediction of breast cancer risk in AA women.

Functional annotation of the 2q35 breast cancer risk locus implicates a structural variant in influencing activity of a long-range enhancer element.

A combination of genetic and functional approaches has identified three independent breast cancer risk loci at 2q35. A recent fine-scale mapping analysis to refine these associations resulted in 1 (signal 1), 5 (signal 2), and 42 (signal 3) credible causal variants at these loci. We used publicly available in silico DNase I and ChIP-seq data with in vitro reporter gene and CRISPR assays to annotate signals 2 and 3. We identified putative regulatory elements that enhanced cell-type-specific transcription from the IGFBP5 promoter at both signals (30- to 40-fold increased expression by the putative regulatory element at signal 2, 2- to 3-fold by the putative regulatory element at signal 3). We further identified one of the five credible causal variants at signal 2, a 1.4 kb deletion (esv3594306), as the likely causal variant; the deletion allele of this variant was associated with an average additional increase in IGFBP5 expression of 1.3-fold (MCF-7) and 2.2-fold (T-47D). We propose a model in which the deletion allele of esv3594306 juxtaposes two transcription factor binding regions (annotated by estrogen receptor alpha ChIP-seq peaks) to generate a single extended regulatory element. This regulatory element increases cell-type-specific expression of the tumor suppressor gene IGFBP5 and, thereby, reduces risk of estrogen receptor-positive breast cancer (odds ratio = 0.77, 95% CI 0.74-0.81, p = 3.1 × 10-31).

Mosquito control exposures and breast cancer risk: analysis of 1071 cases and 2096 controls from the Ghana Breast Health Study.

Epidemiologic data on insecticide exposures and breast cancer risk are inconclusive and mostly from high-income countries. Using data from 1071 invasive pathologically confirmed breast cancer cases and 2096 controls from the Ghana Breast Health Study conducted from 2013 to 2015, we investigated associations with mosquito control products to reduce the spread of mosquito-borne diseases, such as malaria. These mosquito control products were insecticide-treated nets, mosquito coils, repellent room sprays, and skin creams for personal protection against mosquitos. Multivariable and polytomous logistic regression models were used to estimate odds ratios (ORadj) and 95% confidence intervals (CI) with breast cancer risk-adjusted for potential confounders and known risk factors. Among controls, the reported use of mosquito control products were mosquito coils (65%), followed by insecticide-treated nets (56%), repellent room sprays (53%), and repellent skin creams (15%). Compared to a referent group of participants unexposed to mosquito control products, there was no significant association between breast cancer risk and mosquito coils. There was an association in breast cancer risk with reported use of insecticide-treated nets; however, that association was weak and not statistically significant. Participants who reported using repellent sprays were at elevated risks compared to women who did not use any mosquito control products, even after adjustment for all other mosquito control products (OR = 1.42, 95% CI=1.15-1.75). We had limited power to detect an association with repellent skin creams. Although only a few participants reported using repellent room sprays weekly/daily or < month-monthly, no trends were evident with increased frequency of use of repellent sprays, and there was no statistical evidence of heterogeneity by estrogen receptor (ER) status (p-het > 0.25). Our analysis was limited when determining if an association existed with repellent skin creams; therefore, we cannot conclude an association. We found limited evidence of risk associations with widely used mosquito coils and insecticide-treated nets, which are reassuring given their importance for malaria prevention. Our findings regarding specific breast cancer risk associations, specifically those observed between repellent sprays, require further study.

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.

Moving Toward Findable, Accessible, Interoperable, Reusable Practices in Epidemiologic Research.

Data sharing is essential for reproducibility of epidemiologic research, replication of findings, pooled analyses in consortia efforts, and maximizing study value to address multiple research questions. However, barriers related to confidentiality, costs, and incentives often limit the extent and speed of data sharing. Epidemiological practices that follow Findable, Accessible, Interoperable, Reusable (FAIR) principles can address these barriers by making data resources findable with the necessary metadata, accessible to authorized users, and interoperable with other data, to optimize the reuse of resources with appropriate credit to its creators. We provide an overview of these principles and describe approaches for implementation in epidemiology. Increasing degrees of FAIRness can be achieved by moving data and code from on-site locations to remote, accessible ("Cloud") data servers, using machine-readable and nonproprietary files, and developing open-source code. Adoption of these practices will improve daily work and collaborative analyses and facilitate compliance with data sharing policies from funders and scientific journals. Achieving a high degree of FAIRness will require funding, training, organizational support, recognition, and incentives for sharing research resources, both data and code. However, these costs are outweighed by the benefits of making research more reproducible, impactful, and equitable by facilitating the reuse of precious research resources by the scientific community.

Measured body size and serum estrogen metabolism in postmenopausal women: the Ghana Breast Health Study.

BACKGROUND: Several anthropometric measures have been associated with hormone-related cancers, and it has been shown that estrogen metabolism in postmenopausal women plays an important role in these relationships. However, little is known about circulating estrogen levels in African women, and the relevance to breast cancer or breast cancer risk factors. To shed further light on the relationship of anthropometric factors and estrogen levels in African women, we examined whether measured body mass index (BMI), waist-to-hip ratio (WHR), height, and self-reported body size were associated with serum estrogens/estrogen metabolites in a cross-sectional analysis among postmenopausal population-based controls of the Ghana Breast Health Study. METHODS: Fifteen estrogens/estrogen metabolites were quantified using liquid chromatography-tandem mass spectrometry in serum samples collected from postmenopausal female controls enrolled in the Ghana Breast Health Study, a population-based case-control study conducted in Accra and Kumasi. Geometric means (GMs) of estrogens/estrogen metabolites were estimated using linear regression, adjusting for potential confounders. RESULTS: Measured BMI (≥ 30 vs. 18.5-24.9 kg/m2) was positively associated with parent estrogens (multivariable adjusted GM for unconjugated estrone: 78.90 (66.57-93.53) vs. 50.89 (43.47-59.59), p-value < 0.0001; and unconjugated estradiol: 27.83 (21.47-36.07) vs. 13.26 (10.37-16.95), p-value < 0.0001). Independent of unconjugated estradiol, measured BMI was associated with lower levels of 2-pathway metabolites and higher levels of 16-ketoestradriol. Similar patterns of association were found with WHR; however, the associations were not entirely independent of BMI. Height was not associated with postmenopausal estrogens/estrogen metabolite levels in African women. CONCLUSIONS: We observed strong associations between measured BMI and parent estrogens and estrogen metabolite patterns that largely mirrored relations that have previously been associated with higher breast cancer risk in postmenopausal White women. The consistency of the BMI-estrogen metabolism associations in our study with those previously noted among White women suggests that estrogens likely explain part of the BMI-postmenopausal breast cancer risk in both groups. These findings merit evaluation in Black women, including prospective studies.

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.

The oral microbiome and breast cancer and nonmalignant breast disease, and its relationship with the fecal microbiome in the Ghana Breast Health Study.

The oral microbiome, like the fecal microbiome, may be related to breast cancer risk. Therefore, we investigated whether the oral microbiome was associated with breast cancer and nonmalignant breast disease, and its relationship with the fecal microbiome in a case-control study in Ghana. A total of 881 women were included (369 breast cancers, 93 nonmalignant cases and 419 population-based controls). The V4 region of the 16S rRNA gene was sequenced from oral and fecal samples. Alpha-diversity (observed amplicon sequence variants [ASVs], Shannon index and Faith's Phylogenetic Diversity) and beta-diversity (Bray-Curtis, Jaccard and weighted and unweighted UniFrac) metrics were computed. MiRKAT and logistic regression models were used to investigate the case-control associations. Oral sample alpha-diversity was inversely associated with breast cancer and nonmalignant breast disease with odds ratios (95% CIs) per every 10 observed ASVs of 0.86 (0.83-0.89) and 0.79 (0.73-0.85), respectively, compared to controls. Beta-diversity was also associated with breast cancer and nonmalignant breast disease compared to controls (P ≤ .001). The relative abundances of Porphyromonas and Fusobacterium were lower for breast cancer cases compared to controls. Alpha-diversity and presence/relative abundance of specific genera from the oral and fecal microbiome were strongly correlated among breast cancer cases, but weakly correlated among controls. Particularly, the relative abundance of oral Porphyromonas was strongly, inversely correlated with fecal Bacteroides among breast cancer cases (r = -.37, P ≤ .001). Many oral microbial metrics were strongly associated with breast cancer and nonmalignant breast disease, and strongly correlated with fecal microbiome among breast cancer cases, but not controls.

A mixed-model approach for powerful testing of genetic associations with cancer risk incorporating tumor characteristics.

Cancers are routinely classified into subtypes according to various features, including histopathological characteristics and molecular markers. Previous genome-wide association studies have reported heterogeneous associations between loci and cancer subtypes. However, it is not evident what is the optimal modeling strategy for handling correlated tumor features, missing data, and increased degrees-of-freedom in the underlying tests of associations. We propose to test for genetic associations using a mixed-effect two-stage polytomous model score test (MTOP). In the first stage, a standard polytomous model is used to specify all possible subtypes defined by the cross-classification of the tumor characteristics. In the second stage, the subtype-specific case-control odds ratios are specified using a more parsimonious model based on the case-control odds ratio for a baseline subtype, and the case-case parameters associated with tumor markers. Further, to reduce the degrees-of-freedom, we specify case-case parameters for additional exploratory markers using a random-effect model. We use the Expectation-Maximization algorithm to account for missing data on tumor markers. Through simulations across a range of realistic scenarios and data from the Polish Breast Cancer Study (PBCS), we show MTOP outperforms alternative methods for identifying heterogeneous associations between risk loci and tumor subtypes. The proposed methods have been implemented in a user-friendly and high-speed R statistical package called TOP (https://github.com/andrewhaoyu/TOP).

Physical activity, sedentary time and breast cancer risk: a Mendelian randomisation study.

OBJECTIVES: Physical inactivity and sedentary behaviour are associated with higher breast cancer risk in observational studies, but ascribing causality is difficult. Mendelian randomisation (MR) assesses causality by simulating randomised trial groups using genotype. We assessed whether lifelong physical activity or sedentary time, assessed using genotype, may be causally associated with breast cancer risk overall, pre/post-menopause, and by case-groups defined by tumour characteristics. METHODS: We performed two-sample inverse-variance-weighted MR using individual-level Breast Cancer Association Consortium case-control data from 130 957 European-ancestry women (69 838 invasive cases), and published UK Biobank data (n=91 105-377 234). Genetic instruments were single nucleotide polymorphisms (SNPs) associated in UK Biobank with wrist-worn accelerometer-measured overall physical activity (nsnps=5) or sedentary time (nsnps=6), or accelerometer-measured (nsnps=1) or self-reported (nsnps=5) vigorous physical activity. RESULTS: Greater genetically-predicted overall activity was associated with lower breast cancer overall risk (OR=0.59; 95% confidence interval (CI) 0.42 to 0.83 per-standard deviation (SD;~8 milligravities acceleration)) and for most case-groups. Genetically-predicted vigorous activity was associated with lower risk of pre/perimenopausal breast cancer (OR=0.62; 95% CI 0.45 to 0.87,≥3 vs. 0 self-reported days/week), with consistent estimates for most case-groups. Greater genetically-predicted sedentary time was associated with higher hormone-receptor-negative tumour risk (OR=1.77; 95% CI 1.07 to 2.92 per-SD (~7% time spent sedentary)), with elevated estimates for most case-groups. Results were robust to sensitivity analyses examining pleiotropy (including weighted-median-MR, MR-Egger). CONCLUSION: Our study provides strong evidence that greater overall physical activity, greater vigorous activity, and lower sedentary time are likely to reduce breast cancer risk. More widespread adoption of active lifestyles may reduce the burden from the most common cancer in women.

Association of germline genetic variants with breast cancer-specific survival in patient subgroups defined by clinic-pathological variables related to tumor biology and type of systemic treatment.

BACKGROUND: Given the high heterogeneity among breast tumors, associations between common germline genetic variants and survival that may exist within specific subgroups could go undetected in an unstratified set of breast cancer patients. METHODS: We performed genome-wide association analyses within 15 subgroups of breast cancer patients based on prognostic factors, including hormone receptors, tumor grade, age, and type of systemic treatment. Analyses were based on 91,686 female patients of European ancestry from the Breast Cancer Association Consortium, including 7531 breast cancer-specific deaths over a median follow-up of 8.1 years. Cox regression was used to assess associations of common germline variants with 15-year and 5-year breast cancer-specific survival. We assessed the probability of these associations being true positives via the Bayesian false discovery probability (BFDP < 0.15). RESULTS: Evidence of associations with breast cancer-specific survival was observed in three patient subgroups, with variant rs5934618 in patients with grade 3 tumors (15-year-hazard ratio (HR) [95% confidence interval (CI)] 1.32 [1.20, 1.45], P = 1.4E-08, BFDP = 0.01, per G allele); variant rs4679741 in patients with ER-positive tumors treated with endocrine therapy (15-year-HR [95% CI] 1.18 [1.11, 1.26], P = 1.6E-07, BFDP = 0.09, per G allele); variants rs1106333 (15-year-HR [95% CI] 1.68 [1.39,2.03], P = 5.6E-08, BFDP = 0.12, per A allele) and rs78754389 (5-year-HR [95% CI] 1.79 [1.46,2.20], P = 1.7E-08, BFDP = 0.07, per A allele), in patients with ER-negative tumors treated with chemotherapy. CONCLUSIONS: We found evidence of four loci associated with breast cancer-specific survival within three patient subgroups. There was limited evidence for the existence of associations in other patient subgroups. However, the power for many subgroups is limited due to the low number of events. Even so, our results suggest that the impact of common germline genetic variants on breast cancer-specific survival might be limited.

CYP3A7*1C allele: linking premenopausal oestrone and progesterone levels with risk of hormone receptor-positive breast cancers.

BACKGROUND: Epidemiological studies provide strong evidence for a role of endogenous sex hormones in the aetiology of breast cancer. The aim of this analysis was to identify genetic variants that are associated with urinary sex-hormone levels and breast cancer risk. METHODS: We carried out a genome-wide association study of urinary oestrone-3-glucuronide and pregnanediol-3-glucuronide levels in 560 premenopausal women, with additional analysis of progesterone levels in 298 premenopausal women. To test for the association with breast cancer risk, we carried out follow-up genotyping in 90,916 cases and 89,893 controls from the Breast Cancer Association Consortium. All women were of European ancestry. RESULTS: For pregnanediol-3-glucuronide, there were no genome-wide significant associations; for oestrone-3-glucuronide, we identified a single peak mapping to the CYP3A locus, annotated by rs45446698. The minor rs45446698-C allele was associated with lower oestrone-3-glucuronide (-49.2%, 95% CI -56.1% to -41.1%, P = 3.1 × 10-18); in follow-up analyses, rs45446698-C was also associated with lower progesterone (-26.7%, 95% CI -39.4% to -11.6%, P = 0.001) and reduced risk of oestrogen and progesterone receptor-positive breast cancer (OR = 0.86, 95% CI 0.82-0.91, P = 6.9 × 10-8). CONCLUSIONS: The CYP3A7*1C allele is associated with reduced risk of hormone receptor-positive breast cancer possibly mediated via an effect on the metabolism of endogenous sex hormones in premenopausal women.

Mendelian randomisation study of smoking exposure in relation to breast cancer risk.

BACKGROUND: Despite a modest association between tobacco smoking and breast cancer risk reported by recent epidemiological studies, it is still equivocal whether smoking is causally related to breast cancer risk. METHODS: We applied Mendelian randomisation (MR) to evaluate a potential causal effect of cigarette smoking on breast cancer risk. Both individual-level data as well as summary statistics for 164 single-nucleotide polymorphisms (SNPs) reported in genome-wide association studies of lifetime smoking index (LSI) or cigarette per day (CPD) were used to obtain MR effect estimates. Data from 108,420 invasive breast cancer cases and 87,681 controls were used for the LSI analysis and for the CPD analysis conducted among ever-smokers from 26,147 cancer cases and 26,072 controls. Sensitivity analyses were conducted to address pleiotropy. RESULTS: Genetically predicted LSI was associated with increased breast cancer risk (OR 1.18 per SD, 95% CI: 1.07-1.30, P = 0.11 × 10-2), but there was no evidence of association for genetically predicted CPD (OR 1.02, 95% CI: 0.78-1.19, P = 0.85). The sensitivity analyses yielded similar results and showed no strong evidence of pleiotropic effect. CONCLUSION: Our MR study provides supportive evidence for a potential causal association with breast cancer risk for lifetime smoking exposure but not cigarettes per day among smokers.

Discovery of structural deletions in breast cancer predisposition genes using whole genome sequencing data from > 2000 women of African-ancestry.

Single germline nucleotide pathogenic variants have been identified in 12 breast cancer predisposition genes, but structural deletions in these genes remain poorly characterized. We conducted in-depth whole genome sequencing (WGS) in genomic DNA samples obtained from 1340 invasive breast cancer cases and 675 controls of African ancestry. We identified 25 deletions in the intragenic regions of ten established breast cancer predisposition genes based on a consensus call from six state-of-the-art SV callers. Overall, no significant case-control difference was found in the frequency of these deletions. However, 1.0% of cases and 0.3% of controls carried any of the eight putative protein-truncating rare deletions located in BRCA1, BRCA2, CDH1, TP53, NF1, RAD51D, RAD51C and CHEK2, resulting in an odds ratio (OR) of 3.29 (95% CI 0.74-30.16). We also identified a low-frequency deletion in NF1 associated with breast cancer risk (OR 1.93, 95% CI 1.14-3.42). In addition, we detected 56 deletions, including six putative protein-truncating deletions, in suspected breast predisposition genes. This is the first large study to systematically search for structural deletions in breast cancer predisposition genes. Many of the deletions, particularly those resulting in protein truncations, are likely to be pathogenic. Results from this study, if confirmed in future large-scale studies, could have significant implications for genetic testing for this common cancer.

Family history of prostate cancer and the incidence of ERG- and phosphatase and tensin homolog-defined prostate cancer.

Family history is among the strongest known risk factors for prostate cancer (PCa). Emerging data suggest molecular subtypes of PCa, including two somatic genetic aberrations: fusions of androgen-regulated promoters with ERG and, separately, phosphatase and tensin homolog (PTEN) loss. We examined associations between family history and incidence of these subtypes in 44,126 men from the prospective Health Professionals Follow-up Study. ERG and PTEN status were assessed by immunohistochemistry. Multivariable competing risks models were used to estimate hazard ratios (HR) and 95% confidence intervals (CI) for associations between self-reported family history of PCa and molecular subtypes of disease. Thirteen percent of men had a positive family history of PCa at baseline. During a median follow-up of 18.5 years, 5,511 PCa cases were diagnosed. Among them, 888 were assayed for ERG status (47% ERG-positive) and 715 were assayed for PTEN loss (14% PTEN null). Family history was more strongly associated with risk of ERG-negative (HR: 2.15; 95% CI: 1.71-2.70) than ERG-positive (HR: 1.49; 95% CI: 1.13-1.95) disease (pheterogeneity : 0.04). The strongest difference was among men with an affected father (HRERG-negative : 2.09; 95% CI: 1.64-2.66; HRERG-positive : 1.30; 95% CI: 0.96-1.76; pheterogeneity : 0.01). Family history of PCa was positively associated with both PTEN null (HR: 2.10; 95% CI: 1.26-3.49) and PTEN intact (HR: 1.72; 95% CI: 1.39-2.13) PCa (pheterogeneity : 0.47). Our results indicate that PCa family history may be positively associated with PCa in all ERG and PTEN subtypes, suggesting a role of genetic susceptibility in their development. It is possible that ERG-negative disease could be especially associated with positive family history.

Molecular mechanisms linking high body mass index to breast cancer etiology in post-menopausal breast tumor and tumor-adjacent tissues.

PURPOSE: In post-menopausal women, high body mass index (BMI) is an established breast cancer risk factor and is associated with worse breast cancer prognosis. We assessed the associations between BMI and gene expression of both breast tumor and adjacent tissue in estrogen receptor-positive (ER+) and estrogen receptor-negative (ER-) diseases to help elucidate the mechanisms linking obesity with breast cancer biology in 519 post-menopausal women from the Nurses' Health Study (NHS) and NHSII. METHODS: Differential gene expression was analyzed separately in ER+ and ER- disease both comparing overweight (BMI ≥ 25 to < 30) or obese (BMI ≥ 30) women to women with normal BMI (BMI < 25), and per 5 kg/m2 increase in BMI. Analyses controlled for age and year of diagnosis, physical activity, alcohol consumption, and hormone therapy use. Gene set enrichment analyses were performed and validated among a subset of post-menopausal cases in The Cancer Genome Atlas (for tumor) and Polish Breast Cancer Study (for tumor-adjacent). RESULTS: No gene was differentially expressed by BMI (FDR < 0.05). BMI was significantly associated with increased cellular proliferation pathways, particularly in ER+ tumors, and increased inflammation pathways in ER- tumor and ER- tumor-adjacent tissues (FDR < 0.05). High BMI was associated with upregulation of genes involved in epithelial-mesenchymal transition in ER+ tumor-adjacent tissues. CONCLUSIONS: This study provides insights into molecular mechanisms of BMI influencing post-menopausal breast cancer biology. Tumor and tumor-adjacent tissues provide independent information about potential mechanisms.

Detectible mosaic truncating PPM1D mutations, age and breast cancer risk.

Mosaic protein truncating variants (PTVs) in the phosphatase, Mg2+/Mn2+dependent 1D (PPM1D) gene in blood-derived DNA have been associated with increased risk of breast cancer. We analyzed PPM1D PTVs in blood from 3817 breast cancer cases and 3058 controls by deep sequencing of a previously defined region in exon 6 of PPM1D. We identified 50 of 6875 (0.73%) participants having a mosaic PPM1D PTV. We observed a higher frequency of mosaic PPM1D PTVs with increasing age (Ptrend = 2.9 × 10-6). We did not observe an overall association between PPM1D PTVs and increased breast cancer risk (OR = 1.51, 95% CI = 0.84-2.71). Evidence for an association was observed in a subset of cases with DNA collected 1-year or more before breast cancer diagnosis (OR = 3.44, 95% CI = 1.62-7.30, P-value = 0.001); however, no significant association was observed for the larger series of cases with DNA collected post diagnosis (OR = 1.01, 95% CI = 0.51-2.01, P-value = 0.98). Our study indicates that the PPM1D PTVs are present at higher rates than previously reported and the frequency of PPM1D PTVs increases with age. We observed limited evidence for an association between mosaic PPM1D PTVs and breast cancer risk, suggesting mosaic PPM1D PTVs in the blood likely do not influence risk of breast cancer.