Genome-wide association study of germline variants and breast cancer-specific mortality

BACKGROUND

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

METHODS

Meta-analyses included summary estimates based on Cox models of twelve datasets using ~10.4 million variants for 96,661 women with breast cancer and 7697 events (breast cancer-specific deaths). Oestrogen receptor (ER)-specific analyses were based on 64,171 ER-positive (4116) and 16,172 ER-negative (2125) patients. We evaluated the probability of a signal to be a true positive using the Bayesian false discovery probability (BFDP).

RESULTS

We did not find any variant associated with breast cancer-specific mortality at P < 5 × 10-8. For ER-positive disease, the most significantly associated variant was chr7:rs4717568 (BFDP = 7%, P = 1.28 × 10-7, hazard ratio [HR] = 0.88, 95% confidence interval [CI] = 0.84-0.92); the closest gene is AUTS2. For ER-negative disease, the most significant variant was chr7:rs67918676 (BFDP = 11%, P = 1.38 × 10-7, HR = 1.27, 95% CI = 1.16-1.39); located within a long intergenic non-coding RNA gene (AC004009.3), close to the HOXA gene cluster.

CONCLUSIONS

We uncovered germline variants on chromosome 7 at BFDP < 15% close to genes for which there is biological evidence related to breast cancer outcome. However, the paucity of variants associated with mortality at genome-wide significance underpins the challenge in providing genetic-based individualised prognostic information for breast cancer patients.

Antioxidant supplementation and breast cancer prognosis in postmenopausal women undergoing chemotherapy and radiation therapy

Background

There is a paucity of information on the prevalence of dietary supplement use in breast cancer survivors. Only a few studies have examined the impact of dietary supplements, particularly antioxidants, on breast cancer prognosis and the results are inconclusive.

Objective

We examined pre- and postdiagnosis use of supplements in postmenopausal breast cancer survivors in Germany and investigated associations between postdiagnosis use of antioxidants and other supplements, and prognosis (total and breast cancer mortality, and recurrence-free survival) both overall and in women who received chemotherapy and radiation therapy.

Design

Data from 2223 postmenopausal women diagnosed with nonmetastatic breast cancer from the population-based Mamma Carcinoma Risk Factor Investigation (MARIE) study were used. Women were interviewed at recruitment in 2002-2005 and again in 2009 and followed-up until 30 June 2015. Multivariate Cox regression analysis was used to estimate HRs and corresponding 95% CIs.

Results

Pre- and postdiagnosis supplement use was reported by 36% and 45% of the women, respectively. There were 240 deaths (134 from breast cancer) and 200 breast cancer recurrences after a median follow-up time of 6.0 y after the 2009 re-interview. After adjusting for relevant confounders, concurrent antioxidant use with chemotherapy or radiation therapy among 1940 women was associated with increased risk of total mortality (HR: 1.64; 95% CI: 1.01, 2.66) and worsened recurrence-free survival (HR: 1.84; 95% CI: 1.26, 2.68). Overall postdiagnosis supplement use was not associated with breast cancer prognosis.

Conclusions

Antioxidant use during chemotherapy or radiation therapy was associated with worsened breast cancer prognosis in postmenopausal women. There was no overall association between postdiagnosis supplement use and breast cancer prognosis. Results from our study align with the current recommendation to possibly avoid the use of antioxidants during chemotherapy or radiation therapy.

Circulating enterolactone concentrations and prognosis of postmenopausal breast cancer: assessment of mediation by inflammatory markers

Higher lignan exposure has been associated with lower all-cause mortality (ACM) and breast cancer-specific mortality (BCSM) for postmenopausal breast cancer patients. However, the biological mechanisms underpinning these associations are still unclear. We investigated for the first time whether and to what extent the association between enterolactone (ENL), the major lignan metabolite, and postmenopausal breast cancer prognosis is mediated by inflammatory biomarkers. Circulating concentrations of ENL and inflammatory markers were measured in a population-based prospective cohort of 1,743 breast cancer patients recruited between 2002 and 2005 and followed-up until 2009. Hazard ratios (HR) and 95% CIs were estimated using multivariable Cox regression. Mediation analysis was performed to estimate the percentage association between ENL (log2) and ACM, BCSM and distant disease-free survival (DDFS), which is mediated by C-reactive protein (CRP) (log2), as the strongest potential mediator, and also interleukin (IL)-10. Median serum/plasma ENL and CRP concentrations for all patients, including 180 deceased patients, were 23.2 and 17.5 nmol/L, and 3.2 and 6.5 mg/l, respectively. ENL concentrations were significantly inversely associated with ACM, BCSM and DDFS (per doubling of ENL concentrations: HRs 0.93 [0.87, 0.99], 0.91 [0.84, 0.99] and 0.92 [0.87, 0.99]), after adjusting for prognostic factors and BMI. Estimated 18, 14 and 12% of the effects of ENL on ACM, BCSM and DDFS, respectively, were mediated through CRP. No mediational effect of IL-10 was found. We provide first evidence that the proinflammatory marker CRP may partially mediate the association of ENL with postmenopausal breast cancer survival, which supports hormone-independent mechanisms.

A transcriptome-wide association study of 229,000 women identifies new candidate susceptibility genes for breast cancer

The breast cancer risk variants identified in genome-wide association studies explain only a small fraction of the familial relative risk, and the genes responsible for these associations remain largely unknown. To identify novel risk loci and likely causal genes, we performed a transcriptome-wide association study evaluating associations of genetically predicted gene expression with breast cancer risk in 122,977 cases and 105,974 controls of European ancestry. We used data from the Genotype-Tissue Expression Project to establish genetic models to predict gene expression in breast tissue and evaluated model performance using data from The Cancer Genome Atlas. Of the 8,597 genes evaluated, significant associations were identified for 48 at a Bonferroni-corrected threshold of P < 5.82 × 10-6, including 14 genes at loci not yet reported for breast cancer. We silenced 13 genes and showed an effect for 11 on cell proliferation and/or colony-forming efficiency. Our study provides new insights into breast cancer genetics and biology.

Publisher Correction: Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation

This corrects the article DOI: 10.1038/ncomms5999.

Prognostic impact of surgery for early-stage invasive breast cancer on breast cancer-specific survival, overall survival, and recurrence risk: a population-based analysis

PURPOSE

Recent cohort studies demonstrated better overall survival (OS) or breast cancer-specific survival (BCS) for breast-conserving therapy (BCT) followed by radiation (RT) compared to mastectomy alone (MT). This is the first observational study in which adjustments for a comprehensive set of prognostic factors, adjuvant therapies, mode of detection, and comorbidities were possible to investigate OS, BCS, as well as recurrence risk of patients undergoing BCT + RT, MT + RT, or MT.

METHODS

Women aged 50-74 years at diagnosis of early-stage invasive breast cancer (I-IIIa) between 2001 and 2005 at the German population-based case-control study (MARIE study) were recruited and followed prospectively as a case cohort until 2015. Kaplan-Meier estimates and stepwise adjusted multivariable Cox models were used to estimate hazard ratios (HRs) and 95% confidence intervals (95% CI).

RESULTS

The 2762 patients included were followed up for a median of 11.9 years (95% CI 11.8-12.0). 74.2% of patients underwent BCT + RT; 10.3% MT + RT and 15.6% MT alone. Compared to patients treated with MT alone, patients treated with BCT + RT showed non-statistically significant improved OS (HR 0.79, 95% CI 0.61-1.02), BCS (HR 0.79, 95% CI 0.55-1.12), and no difference in recurrence risks (HR 1.01, 95% CI 0.74-1.37). For patients treated with MT + RT, there were no differences in OS (HR 1.06, 95% CI 0.75-1.50), BCS (HR 1.17, 95% CI 0.75-1.82), or recurrence risk (HR 1.33, 95% CI 0.89-1.97).

CONCLUSIONS

Among patients with early-stage breast cancer, clinical outcomes more than 10 years after diagnosis did not differ between the primary treatment options BCT + RT, MT + RT versus MT alone after full adjustment.

Physical Activity and Postmenopausal Breast Cancer

Objectives: Epidemiological evidence suggests an inverse association between physical activity (PA) and postmenopausal breast cancer risk. Breast cancer is a heterogeneous disease, influenced by reproductive factors, lifestyle pattern, and predispositions. We investigated whether these risk factors modify the effect of PA on breast cancer risk.

Methods: We analyzed data from 2004 hormone-receptor-positive postmenopausal breast cancer cases and 6569 controls from the population-based MARIE study conducted 2002–2005 in Germany. Interaction was statistically tested using adjusted unconditional logistic regression models.

Results: The inverse association between leisure-time PA and risk of postmenopausal hormone-receptor-positive breast cancer was not heterogeneous by family history of breast cancer or by hormone therapy. PA showed a significant interaction with benign breast diseases (p = 0.023) and with breastfeeding (p = 0.045) but not with parity (p = 0.94), with clear risk reductions only for women who ever had breastfed or who ever had a benign breast disease (among ever breastfed: odds ratio = 0.63; 95% confidence interval = (0.52, 0.77), highest vs. lowest PA quartile). Interaction with BMI was weak (p = 0.053).

Conclusions: Breastfeeding and benign breast diseases modified the effect of PA on postmenopausal breast cancer risk. If other studies find similar modifications, increasing knowledge about these risk factors may contribute to a better understanding of the mode of action of PA on breast cancer risk. For women who are at higher risk for breast cancer due to family history or due to hormone therapy use, it is encouraging that they might lower their risk by being physically active.

Associations between adjuvant radiotherapy and different causes of death in a German breast cancer cohort

BACKGROUND

Studies of cohorts of breast cancer (BC) patients diagnosed before 1990 showed radiotherapy (RT) to be associated with increased cardiovascular (CVD) and lung cancer mortality many years after diagnosis. In the late 1990s, improvements in RT planning techniques reduced radiation doses to normal tissues. Recent studies did not consistently report higher RT-related mortality for CVD and second cancers. Aim of the study was to analyze specific causes of death after 3D-conformal RT in a recent BC cohort.

METHODS

Stage I-III BC patients diagnosed 2001-2005 and enrolled in the population based MARIEplus study were followed-up for 11.9 years (median). Associations between adjuvant RT and cause-specific mortality were analyzed by using competing risks models, yielding subdistribution hazard ratios (SHR) for RT directly related to cumulative incidences. Models were adjusted for differences in baseline characteristics applying inverse-probability-of-treatment-weighting (IPTW).

RESULTS

Of the 2951 patients, 2439 (83.0%) received RT. No significant association of RT with lung cancer mortality (SHR_IPTW 0.88, 0.35-2.12), other cancer mortality (SHR_IPTW 1.04, 95% CI 0.62-1.73) or cardiac mortality was observed (SHR_IPTW 1.57, 0.75-3.29). Mortality from lung and other diseases were significantly lower in irradiated women (SHR_IPTW 0.39, 95% CI 0.17-0.90 and SHR_IPTW 0.58, 95% CI 0.34-0.97, respectively).

CONCLUSION

In line with recent studies, 3D-conformal RT did not significantly increase mortality from non-BC causes in the German MARIEplus cohort. Since long-term data are still sparse and event rates low in BC-cohorts, who received modern RT, investigation of possible late RT effects on mortality beyond 14 years of follow-up is warranted.

Body mass index and breast cancer survival: a Mendelian randomization analysis

Background

There is increasing evidence that elevated body mass index (BMI) is associated with reduced survival for women with breast cancer. However, the underlying reasons remain unclear. We conducted a Mendelian randomization analysis to investigate a possible causal role of BMI in survival from breast cancer.

Methods

We used individual-level data from six large breast cancer case-cohorts including a total of 36 210 individuals (2475 events) of European ancestry. We created a BMI genetic risk score (GRS) based on genotypes at 94 known BMI-associated genetic variants. Association between the BMI genetic score and breast cancer survival was analysed by Cox regression for each study separately. Study-specific hazard ratios were pooled using fixed-effect meta-analysis.

Results

BMI genetic score was found to be associated with reduced breast cancer-specific survival for estrogen receptor (ER)-positive cases [hazard ratio (HR) = 1.11, per one-unit increment of GRS, 95% confidence interval (CI) 1.01-1.22, P = 0.03). We observed no association for ER-negative cases (HR = 1.00, per one-unit increment of GRS, 95% CI 0.89-1.13, P = 0.95).

Conclusions

Our findings suggest a causal effect of increased BMI on reduced breast cancer survival for ER-positive breast cancer. There is no evidence of a causal effect of higher BMI on survival for ER-negative breast cancer cases.

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.

PHIP - a novel candidate breast cancer susceptibility locus on 6q14.1

Most non-BRCA1/2 breast cancer families have no identified genetic cause. We used linkage and haplotype analyses in familial and sporadic breast cancer cases to identify a susceptibility locus on chromosome 6q. Two independent genome-wide linkage analysis studies suggested a 3 Mb locus on chromosome 6q and two unrelated Swedish families with a LOD >2 together seemed to share a haplotype in 6q14.1. We hypothesized that this region harbored a rare high-risk founder allele contributing to breast cancer in these two families. Sequencing of DNA and RNA from the two families did not detect any pathogenic mutations. Finally, 29 SNPs in the region were analyzed in 44,214 cases and 43,532 controls from BCAC, and the original haplotypes in the two families were suggested as low-risk alleles for European and Swedish women specifically. There was also some support for one additional independent moderate-risk allele in Swedish familial samples. The results were consistent with our previous findings in familial breast cancer and supported a breast cancer susceptibility locus at 6q14.1 around the PHIP gene.

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.

Anti-nuclear autoantibodies in the general German population: prevalence and lack of association with selected cardiovascular and metabolic disorders-findings of a multicenter population-based study

BACKGROUND

We determined the prevalence of anti-nuclear autoantibodies (ANAs) in the German adult population and examined the association between ANAs and cardiovascular and metabolic disorders.

METHODS

We used data and blood samples from the pretest phases of the German National Cohort, obtained from six of the 18 study centers (n = 1199). All centers applied standardized instruments including face-to-face interviews, anthropometric measurements and collection of blood samples. Self-reported histories of diabetes mellitus, heart attack and elevated blood cholesterol and/or lipids were recorded. Height, weight and blood pressure were measured. ANAs were detected using a semi-automated system (AKLIDES®; Medipan GmbH, Dahlewitz, Germany). A positive ANA was defined as a titer ≥ 1:80. ANA were classified as weakly (1:80 or 1:160), moderately (1:320 or 1:640) or strongly (≥1:1280) positive. Specific autoantibodies against nuclear antigens were detected with second-step assays according to the ANA staining pattern. Associations between the assessed disorders and ANA positivity and pattern were examined using sex and age-adjusted mixed-effects logistic regression models.

RESULTS

Thirty-three percent (95% confidence interval; 31-36%) of the 1196 participants (measurements could not be obtained from three samples) were ANA positive (titer ≥ 1:80). The proportions of weakly, moderately and strongly positive ANA were 29%, 3.3% and 1.3%, respectively. ANA positivity was more common among women than men across all titers (χ², p = 0.03). ANA positivity, even when stratified according to height of titer or immunofluorescent pattern, was not associated with diabetes, elevated blood cholesterol and/or lipids, obesity or hypertension. Second-step autoantibody assays were positive in 41 of the 83 samples (49%) tested, with anti-DFS70 (n = 13) and anti-dsDNA (n = 7) being most frequent. These subgroups were too small to test for associations with the disorders assessed.

CONCLUSIONS

The prevalence of ANA positivity in the German general population was similar to values reported from other countries. Contrary to other studies, there was no association with selected self-reported and objectively measured cardiovascular and metabolic variables.

Genomic analyses identify hundreds of variants associated with age at menarche and support a role for puberty timing in cancer risk

The timing of puberty is a highly polygenic childhood trait that is epidemiologically associated with various adult diseases. Using 1000 Genomes Project-imputed genotype data in up to ∼370,000 women, we identify 389 independent signals (P < 5 × 10-8) for age at menarche, a milestone in female pubertal development. In Icelandic data, these signals explain ∼7.4% of the population variance in age at menarche, corresponding to ∼25% of the estimated heritability. We implicate ∼250 genes via coding variation or associated expression, demonstrating significant enrichment in neural tissues. Rare variants near the imprinted genes MKRN3 and DLK1 were identified, exhibiting large effects when paternally inherited. Mendelian randomization analyses suggest causal inverse associations, independent of body mass index (BMI), between puberty timing and risks for breast and endometrial cancers in women and prostate cancer in men. In aggregate, our findings highlight the complexity of the genetic regulation of puberty timing and support causal links with cancer susceptibility.

Determination of nasal and oropharyngeal microbiomes in a multicenter population-based study - findings from Pretest 1 of the German National Cohort

We examined acceptability, preference and feasibility of collecting nasal and oropharyngeal swabs, followed by microbiome analysis, in a population-based study with 524 participants. Anterior nasal and oropharyngeal swabs were collected by certified personnel. In addition, participants self-collected nasal swabs at home four weeks later. Four swab types were compared regarding (1) participants' satisfaction and acceptance and (2) detection of microbial community structures based on deep sequencing of the 16 S rRNA gene V1-V2 variable regions. All swabbing methods were highly accepted. Microbial community structure analysis revealed 846 phylotypes, 46 of which were unique to oropharynx and 164 unique to nares. The calcium alginate tipped swab was found unsuitable for microbiome determinations. Among the remaining three swab types, there were no differences in oropharyngeal microbiomes detected and only marginal differences in nasal microbiomes. Microbial community structures did not differ between staff-collected and self-collected nasal swabs. These results suggest (1) that nasal and oropharyngeal swabbing are highly feasible methods for human population-based studies that include the characterization of microbial community structures in these important ecological niches, and (2) that self-collection of nasal swabs at home can be used to reduce cost and resources needed, particularly when serial measurements are to be taken.

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.

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

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

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

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

SIGNIFICANCE

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

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

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

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

PURPOSE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

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

OBJECTIVE

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Identification of four novel susceptibility loci for oestrogen receptor negative breast cancer

Common variants in 94 loci have been associated with breast cancer including 15 loci with genome-wide significant associations (P<5 × 10(-8)) with oestrogen receptor (ER)-negative breast cancer and BRCA1-associated breast cancer risk. In this study, to identify new ER-negative susceptibility loci, we performed a meta-analysis of 11 genome-wide association studies (GWAS) consisting of 4,939 ER-negative cases and 14,352 controls, combined with 7,333 ER-negative cases and 42,468 controls and 15,252 BRCA1 mutation carriers genotyped on the iCOGS array. We identify four previously unidentified loci including two loci at 13q22 near KLF5, a 2p23.2 locus near WDR43 and a 2q33 locus near PPIL3 that display genome-wide significant associations with ER-negative breast cancer. In addition, 19 known breast cancer risk loci have genome-wide significant associations and 40 had moderate associations (P<0.05) with ER-negative disease. Using functional and eQTL studies we implicate TRMT61B and WDR43 at 2p23.2 and PPIL3 at 2q33 in ER-negative breast cancer aetiology. All ER-negative loci combined account for ∼11% of familial relative risk for ER-negative disease and may contribute to improved ER-negative and BRCA1 breast cancer risk prediction.

A Comprehensive Multistate Model Analyzing Associations of Various Risk Factors With the Course of Breast Cancer in a Population-Based Cohort of Breast Cancer Cases

We employed a semi-Markov multistate model for the simultaneous analysis of various endpoints describing the course of breast cancer. Results were compared with those from standard analyses using a Cox proportional hazards model. We included 3,012 patients with invasive breast cancer newly diagnosed between 2001 and 2005 who were recruited in Germany for a population-based study, the Mamma Carcinoma Risk Factor Investigation (MARIE Study), and prospectively followed up until the end of 2009. Locoregional recurrence and distant metastasis were included as intermediate states, and deaths from breast cancer, secondary cancer, and other causes were included as competing absorbing states. Tumor characteristics were significantly associated with all breast cancer-related endpoints. Nodal involvement was significantly related to local recurrence but more strongly related to distant metastases. Smoking was significantly associated with mortality from second cancers and other causes, whereas menopausal hormone use was significantly associated with reduced distant metastasis and death from causes other than cancer. The presence of cardiovascular disease at diagnosis was solely associated with mortality from other causes. Compared with separate Cox models, multistate models allow for dissection of prognostic factors and intermediate events in the analysis of cause-specific mortality and can yield new insights into disease progression and associated pathways.

BRCA2 Polymorphic Stop Codon K3326X and the Risk of Breast, Prostate, and Ovarian Cancers

BACKGROUND

The K3326X variant in BRCA2 (BRCA2*c.9976A>T; p.Lys3326*; rs11571833) has been found to be associated with small increased risks of breast cancer. However, it is not clear to what extent linkage disequilibrium with fully pathogenic mutations might account for this association. There is scant information about the effect of K3326X in other hormone-related cancers.

METHODS

Using weighted logistic regression, we analyzed data from the large iCOGS study including 76 637 cancer case patients and 83 796 control patients to estimate odds ratios (ORw) and 95% confidence intervals (CIs) for K3326X variant carriers in relation to breast, ovarian, and prostate cancer risks, with weights defined as probability of not having a pathogenic BRCA2 variant. Using Cox proportional hazards modeling, we also examined the associations of K3326X with breast and ovarian cancer risks among 7183 BRCA1 variant carriers. All statistical tests were two-sided.

RESULTS

The K3326X variant was associated with breast (ORw = 1.28, 95% CI = 1.17 to 1.40, P = 5.9x10(-) (6)) and invasive ovarian cancer (ORw = 1.26, 95% CI = 1.10 to 1.43, P = 3.8x10(-3)). These associations were stronger for serous ovarian cancer and for estrogen receptor-negative breast cancer (ORw = 1.46, 95% CI = 1.2 to 1.70, P = 3.4x10(-5) and ORw = 1.50, 95% CI = 1.28 to 1.76, P = 4.1x10(-5), respectively). For BRCA1 mutation carriers, there was a statistically significant inverse association of the K3326X variant with risk of ovarian cancer (HR = 0.43, 95% CI = 0.22 to 0.84, P = .013) but no association with breast cancer. No association with prostate cancer was observed.

CONCLUSIONS

Our study provides evidence that the K3326X variant is associated with risk of developing breast and ovarian cancers independent of other pathogenic variants in BRCA2. Further studies are needed to determine the biological mechanism of action responsible for these associations.

Genetic variation in the immunosuppression pathway genes and breast cancer susceptibility: a pooled analysis of 42,510 cases and 40,577 controls from the Breast Cancer Association Consortium

Immunosuppression plays a pivotal role in assisting tumors to evade immune destruction and promoting tumor development. We hypothesized that genetic variation in the immunosuppression pathway genes may be implicated in breast cancer tumorigenesis. We included 42,510 female breast cancer cases and 40,577 controls of European ancestry from 37 studies in the Breast Cancer Association Consortium (2015) with available genotype data for 3595 single nucleotide polymorphisms (SNPs) in 133 candidate genes. Associations between genotyped SNPs and overall breast cancer risk, and secondarily according to estrogen receptor (ER) status, were assessed using multiple logistic regression models. Gene-level associations were assessed based on principal component analysis. Gene expression analyses were conducted using RNA sequencing level 3 data from The Cancer Genome Atlas for 989 breast tumor samples and 113 matched normal tissue samples. SNP rs1905339 (A>G) in the STAT3 region was associated with an increased breast cancer risk (per allele odds ratio 1.05, 95 % confidence interval 1.03-1.08; p value = 1.4 × 10(-6)). The association did not differ significantly by ER status. On the gene level, in addition to TGFBR2 and CCND1, IL5 and GM-CSF showed the strongest associations with overall breast cancer risk (p value = 1.0 × 10(-3) and 7.0 × 10(-3), respectively). Furthermore, STAT3 and IL5 but not GM-CSF were differentially expressed between breast tumor tissue and normal tissue (p value = 2.5 × 10(-3), 4.5 × 10(-4) and 0.63, respectively). Our data provide evidence that the immunosuppression pathway genes STAT3, IL5, and GM-CSF may be novel susceptibility loci for breast cancer in women of European ancestry.

Relationship between menopausal hormone therapy and mortality after breast cancer The MARIEplus study, a prospective case cohort

Cohort studies of breast cancer (BC) patients, but not of disease-free women at inclusion, have found menopausal hormone therapy (MHT) to be associated with decreased BC specific mortality (BCM). Here, the German population-based MARIEplus BC cohort was analyzed to further elucidate associations of prediagnostic MHT with BCM (and modification by tumor characteristics), recurrence, and secondarily with other cause and overall mortality. Enrolled 2002-2005, incident invasive BC cases (N = 3,321) were followed up for a median of 6.1 years. Cox proportional hazards models adjusted for tumor characteristics, mammography and lifestyle were applied. Compared with never users of MHT, current users at date of diagnosis had significantly lower BCM (Hazard ratio (HR) 0.72, 95% CI 0.53-0.97) and risk of recurrence (HR 0.61, 95% CI 0.46-0.82). The MHT related reduced BCM was confined to patients with low grade tumors (HR 0.44, 95% CI 0.28-0.70; phet  = 0.01) and not modified by estrogen receptor or nodal status. BCM decreased with MHT duration in current and increased in past users (phet  = 0.015). Mortality due to causes other than BC and overall mortality were also reduced in current MHT users (HR 0.51, 95% CI 0.32-0.81, HR 0.66, 95% CI 0.52-0.86, respectively). Favorable tumor characteristics and mammographic surveillance could not fully explain associations of current MHT use with BCM and recurrence risk. Thus, the study contributes to the evidence that prediagnostic MHT does not have a negative impact on prognosis after BC. The restriction of a reduced BCM to low grade tumors should be confirmed in independent studies.

A polymorphism in the base excision repair gene PARP2 is associated with differential prognosis by chemotherapy among postmenopausal breast cancer patients

BACKGROUND

Personalized therapy considering clinical and genetic patient characteristics will further improve breast cancer survival. Two widely used treatments, chemotherapy and radiotherapy, can induce oxidative DNA damage and, if not repaired, cell death. Since base excision repair (BER) activity is specific for oxidative DNA damage, we hypothesized that germline genetic variation in this pathway will affect breast cancer-specific survival depending on treatment.

METHODS

We assessed in 1,408 postmenopausal breast cancer patients from the German MARIE study whether cancer specific survival after adjuvant chemotherapy, anthracycline chemotherapy, and radiotherapy is modulated by 127 Single Nucleotide Polymorphisms (SNPs) in 21 BER genes. For SNPs with interaction terms showing p<0.1 (likelihood ratio test) using multivariable Cox proportional hazard analyses, replication in 6,392 patients from nine studies of the Breast Cancer Association Consortium (BCAC) was performed.

RESULTS

rs878156 in PARP2 showed a differential effect by chemotherapy (p=0.093) and was replicated in BCAC studies (p=0.009; combined analysis p=0.002). Compared to non-carriers, carriers of the variant G allele (minor allele frequency=0.07) showed better survival after chemotherapy (combined allelic hazard ratio (HR)=0.75, 95% 0.53-1.07) and poorer survival when not treated with chemotherapy (HR=1.42, 95% 1.08-1.85). A similar effect modification by rs878156 was observed for anthracycline-based chemotherapy in both MARIE and BCAC, with improved survival in carriers (combined allelic HR=0.73, 95% CI 0.40-1.32). None of the SNPs showed significant differential effects by radiotherapy.

CONCLUSIONS

Our data suggest for the first time that a SNP in PARP2, rs878156, may together with other genetic variants modulate cancer specific survival in breast cancer patients depending on chemotherapy. These germline SNPs could contribute towards the design of predictive tests for breast cancer patients.

Height and Breast Cancer Risk: Evidence From Prospective Studies and Mendelian Randomization

BACKGROUND

Epidemiological studies have linked adult height with breast cancer risk in women. However, the magnitude of the association, particularly by subtypes of breast cancer, has not been established. Furthermore, the mechanisms of the association remain unclear.

METHODS

We performed a meta-analysis to investigate associations between height and breast cancer risk using data from 159 prospective cohorts totaling 5216302 women, including 113178 events. In a consortium with individual-level data from 46325 case patients and 42482 control patients, we conducted a Mendelian randomization analysis using a genetic score that comprised 168 height-associated variants as an instrument. This association was further evaluated in a second consortium using summary statistics data from 16003 case patients and 41335 control patients.

RESULTS

The pooled relative risk of breast cancer was 1.17 (95% confidence interval [CI] = 1.15 to 1.19) per 10cm increase in height in the meta-analysis of prospective studies. In Mendelian randomization analysis, the odds ratio of breast cancer per 10cm increase in genetically predicted height was 1.22 (95% CI = 1.13 to 1.32) in the first consortium and 1.21 (95% CI = 1.05 to 1.39) in the second consortium. The association was found in both premenopausal and postmenopausal women but restricted to hormone receptor-positive breast cancer. Analyses of height-associated variants identified eight new loci associated with breast cancer risk after adjusting for multiple comparisons, including three loci at 1q21.2, DNAJC27, and CCDC91 at genome-wide significance level P < 5×10(-8).

CONCLUSIONS

Our study provides strong evidence that adult height is a risk factor for breast cancer in women and certain genetic factors and biological pathways affecting adult height have an important role in the etiology of breast cancer.

Large-scale genomic analyses link reproductive aging to hypothalamic signaling, breast cancer susceptibility and BRCA1-mediated DNA repair

Menopause timing has a substantial impact on infertility and risk of disease, including breast cancer, but the underlying mechanisms are poorly understood. We report a dual strategy in ∼70,000 women to identify common and low-frequency protein-coding variation associated with age at natural menopause (ANM). We identified 44 regions with common variants, including two regions harboring additional rare missense alleles of large effect. We found enrichment of signals in or near genes involved in delayed puberty, highlighting the first molecular links between the onset and end of reproductive lifespan. Pathway analyses identified major association with DNA damage response (DDR) genes, including the first common coding variant in BRCA1 associated with any complex trait. Mendelian randomization analyses supported a causal effect of later ANM on breast cancer risk (∼6% increase in risk per year; P = 3 × 10(-14)), likely mediated by prolonged sex hormone exposure rather than DDR mechanisms.

Evaluation of a questionnaire to assess selected infectious diseases and their risk factors : findings of a multicenter study

Background/objectives

The risk to die from an infectious disease in Germany has been continuously decreasing over the last century. Since infections are, however, not only causes of death but risk factors for diseases like cardiovascular diseases, it is essential to monitor and analyze their prevalence and frequency, especially in consideration of the increased life expectancy. To gain more knowledge about infectious diseases as risk factors and their implications on the condition and change of the immune status, the German National Cohort (GNC), a population-based prospective cohort study, will recruit 200,000 subjects between 2014 and 2017. In Pretest 1, a feasibility study for the GNC, we evaluated a self-administered and self-report questionnaire on infectious diseases and on the use of health care facilities (hereinafter called “ID Screen”) for feasibility and validity.

Methods

From August–November 2011, 435 participants between the ages of 20–69 completed the ID Screen. All subjects had been recruited via a random sample from the local residents’ registration offices by 4 of the 18 participating study centers. The questionnaire encompasses 77 variables in six sections assessing items such as 12-month prevalence of infections, cumulative prevalence of infectious diseases, visit of health care facilities and vaccination. The feasibility was amongst others evaluated by assessing the completeness and comprehensiveness of the questionnaire. To assess the questionnaires ability to measure “immune status” and “susceptibility to infections”, multivariate analysis was used.

Results

The overall practicability was good and most items were well understood, demonstrated by < 2/33 missing questions per questionnaire and only three variables: vaccination for influenza and pneumococci and infection with chickenpox had a frequency > 5 % of missing values. However, direct comparison of the items 12-month prevalence and lifetime prevalence of nephritis/pyelitis showed poor agreement and thereby poor understanding by 80 % of the participants, illustrating the necessity for a clear, lay person appropriate description of rare diseases to increase comprehensibility. The questionnaire will be used to support the assessment of immune dysfunction and frequency of infection. An analysis of these constructs in an exploratory factor analysis revealed limited applicability due to low interitem correlation (Cronbach’s α < 0.5). This is corroborated by the extraction of more than one factor with a Kaiser–Meyer–Olkin measure of 0.6 instead of a unidimensional latent construct for “immune status”.

Conclusion

All in all, the ID Screen is a good and reliable tool to measure infectious diseases as risk factors and outcome in general, but requires a better translation of infection specific terms into lay person terms. For the assessment of the overall immune status, the tool has strong limitations. Vaccinations status should also rather be assessed based on vaccination certificates than on participants’ recall.

Electronic supplementary material

The online version of this article (doi: 10.1007/s00103-014-2052-y) contains supplementary material, which is available to authorized users.

Feasibility and acceptance of cervicovaginal self-sampling within the German National Cohort (Pretest 2)

Background and objectives

Within the German National Cohort (GNC) 100,000 adult women in Germany will be comprehensively interviewed and examined. While women’s health is addressed in the basic interview, direct detection of cervicovaginal microbial colonisation or infection is not part of the examination protocol. In a pilot project the feasibility of female study participants of the GNC collecting a cervicovaginal lavage at home without having to involve a gynecologist or other medical personnel was thus investigated. The ability of the procedure to detect vaginal microbes and conditions including human papillomavirus (HPV), Chlamydia trachomatis and bacterial vaginosis (BV) were also explored.

Methods

This cross-sectional study was conducted in two study centers (Hamburg and Hanover) of the GNC during Pretest 2 in 2012 as an add-on module to the main program of the National Cohort. Participants were randomly selected through the population registration office. After providing written informed consent at the study center, participants self-collected a cervicovaginal lavage (Delphi Screener™) at home following written instructions. Participants mailed samples and acceptability questionnaires to the laboratory and the study center, respectively. Acceptability of self-sampling was categorized as consent, partial consent and rejection. The samples were analyzed by multiplex HPV genotyping for the presence of 27 mucosal HPV subtypes. To detect other pathogens “Sexually Transmitted Infection Profiling” (STIP) was used, a novel multiplex polymerase chain reaction (PCR) for various vaginally occurring pathogens/conditions coupled with subsequent bead-based Luminex^® hybridization. Human beta-globin and DNA polymerase alpha (PolA) sequences were used as positive controls for the detection of human DNA during HPV detection and STIP, respectively.

Results

The participation based on the proportion of all women in Pretest 2 who could take part in the add-on Pretest 2 was 67.3 % (109 out of 162). The age of participants ranged from 20 to 69 years. The self-reported median duration of the collection of the lavage was 5 min. Analysis of the questionnaires (n = 108) revealed that the self-sampling of a cervicovaginal lavage was acceptable to 98 % of women (106 out of 108), and considered to be easy by 89 % (96 out of 108) as well as user-friendly by 96 % of the women (104 out of 108). Human beta-globin and PolA as markers for human DNA and sample quality were detected in all samples analyzed while HPV as a marker for pathogen detectability was identified in 18 out of 109 samples. Of the 107 samples tested with STIP as a second marker for pathogen detectability, 5 samples were excluded from statistical analyses on bacterial colonization because of signs in the laboratory results of the use of antibiotics. For the computation of the possible occurrence of bacterial vaginosis and candidiasis 7 and 8 samples, respectively, were excluded because of low signal intensities resulting in an evaluation of 95 or 94 samples, respectively. Ureaplasma parvum was detected in 22 out of 102 samples, BV in 14 out of 95 samples and candidiasis in 13 out of 94 samples. Chlamydia trachomatis was not detected in any sample.

Conclusion

The feasibility study on cervicovaginal self-sampling indicates that this form of biosampling was very well accepted within the framework of the GNC and feasible in terms of pathogen detection. Its further application in the GNC would allow investigation of transience and persistence, or long-term effects of vaginal (co)infections and colonization.

Dietary inflammation potential and postmenopausal breast cancer risk in a German case-control study

Unhealthy dietary habits can increase the risk for serious medical conditions, such as cancer, yet the association between diet and breast cancer remains unclear. We investigated whether individual diets based on their inflammatory potential are associated with postmenopausal breast cancer risk by employing an energy-adjusted dietary inflammation index. In a German population-based case-control study, 2887 postmenopausal breast cancer patients (aged 50-74 years, first diagnosed between 2002 and 2005) and 5512 healthy age-matched controls provided information on dietary habits for the year prior to diagnosis (cases) or recruitment (controls) using a 176-items food frequency questionnaire. Associations between the energy-adjusted dietary inflammation index (E-DII) score (both as continuous variable and in quintiles) and risk for breast cancer were assessed using conditional logistic regression adjusted for potential confounders. No significant associations between the E-DII score and postmenopausal breast cancer risk were observed (adjusted OR Q5 vs Q1: 1.01, 95% CI: 0.86-1.17). Associations did not differ by estrogen receptor/progesterone receptor status (ER + PR+: adjusted OR Q5 vs Q1: 1.06, 95% CI: 0.88-1.27; ER + or PR+: OR Q5 vs Q1: 1,07, 95% CI: 0.79-1.45; ER-PR-: OR Q5 vs Q1: 0.87 95% CI: 0.63-1.20). Our results regarding E-DII are consistent with previous studies reporting a lack of association between C-reactive protein, a marker of systemic inflammation, and postmenopausal breast cancer risk. The findings may reflect a real absence of association between dietary inflammatory potential and postmenopausal cancer risk or an underestimation of association due to recall bias. Further investigation is warranted in cohort studies.

Fine-mapping identifies two additional breast cancer susceptibility loci at 9q31.2

We recently identified a novel susceptibility variant, rs865686, for estrogen-receptor positive breast cancer at 9q31.2. Here, we report a fine-mapping analysis of the 9q31.2 susceptibility locus using 43 160 cases and 42 600 controls of European ancestry ascertained from 52 studies and a further 5795 cases and 6624 controls of Asian ancestry from nine studies. Single nucleotide polymorphism (SNP) rs676256 was most strongly associated with risk in Europeans (odds ratios [OR] = 0.90 [0.88-0.92]; P-value = 1.58 × 10(-25)). This SNP is one of a cluster of highly correlated variants, including rs865686, that spans ∼14.5 kb. We identified two additional independent association signals demarcated by SNPs rs10816625 (OR = 1.12 [1.08-1.17]; P-value = 7.89 × 10(-09)) and rs13294895 (OR = 1.09 [1.06-1.12]; P-value = 2.97 × 10(-11)). SNP rs10816625, but not rs13294895, was also associated with risk of breast cancer in Asian individuals (OR = 1.12 [1.06-1.18]; P-value = 2.77 × 10(-05)). Functional genomic annotation using data derived from breast cancer cell-line models indicates that these SNPs localise to putative enhancer elements that bind known drivers of hormone-dependent breast cancer, including ER-α, FOXA1 and GATA-3. In vitro analyses indicate that rs10816625 and rs13294895 have allele-specific effects on enhancer activity and suggest chromatin interactions with the KLF4 gene locus. These results demonstrate the power of dense genotyping in large studies to identify independent susceptibility variants. Analysis of associations using subjects with different ancestry, combined with bioinformatic and genomic characterisation, can provide strong evidence for the likely causative alleles and their functional basis.

Identification of novel genetic markers of breast cancer survival

BACKGROUND

Survival after a diagnosis of breast cancer varies considerably between patients, and some of this variation may be because of germline genetic variation. We aimed to identify genetic markers associated with breast cancer-specific survival.

METHODS

We conducted a large meta-analysis of studies in populations of European ancestry, including 37954 patients with 2900 deaths from breast cancer. Each study had been genotyped for between 200000 and 900000 single nucleotide polymorphisms (SNPs) across the genome; genotypes for nine million common variants were imputed using a common reference panel from the 1000 Genomes Project. We also carried out subtype-specific analyses based on 6881 estrogen receptor (ER)-negative patients (920 events) and 23059 ER-positive patients (1333 events). All statistical tests were two-sided.

RESULTS

We identified one new locus (rs2059614 at 11q24.2) associated with survival in ER-negative breast cancer cases (hazard ratio [HR] = 1.95, 95% confidence interval [CI] = 1.55 to 2.47, P = 1.91 x 10(-8)). Genotyping a subset of 2113 case patients, of which 300 were ER negative, provided supporting evidence for the quality of the imputation. The association in this set of case patients was stronger for the observed genotypes than for the imputed genotypes. A second locus (rs148760487 at 2q24.2) was associated at genome-wide statistical significance in initial analyses; the association was similar in ER-positive and ER-negative case patients. Here the results of genotyping suggested that the finding was less robust.

CONCLUSIONS

This is currently the largest study investigating genetic variation associated with breast cancer survival. Our results have potential clinical implications, as they confirm that germline genotype can provide prognostic information in addition to standard tumor prognostic factors.

Prediction of breast cancer risk based on profiling with common genetic variants

BACKGROUND

Data for multiple common susceptibility alleles for breast cancer may be combined to identify women at different levels of breast cancer risk. Such stratification could guide preventive and screening strategies. However, empirical evidence for genetic risk stratification is lacking.

METHODS

We investigated the value of using 77 breast cancer-associated single nucleotide polymorphisms (SNPs) for risk stratification, in a study of 33 673 breast cancer cases and 33 381 control women of European origin. We tested all possible pair-wise multiplicative interactions and constructed a 77-SNP polygenic risk score (PRS) for breast cancer overall and by estrogen receptor (ER) status. Absolute risks of breast cancer by PRS were derived from relative risk estimates and UK incidence and mortality rates.

RESULTS

There was no strong evidence for departure from a multiplicative model for any SNP pair. Women in the highest 1% of the PRS had a three-fold increased risk of developing breast cancer compared with women in the middle quintile (odds ratio [OR] = 3.36, 95% confidence interval [CI] = 2.95 to 3.83). The ORs for ER-positive and ER-negative disease were 3.73 (95% CI = 3.24 to 4.30) and 2.80 (95% CI = 2.26 to 3.46), respectively. Lifetime risk of breast cancer for women in the lowest and highest quintiles of the PRS were 5.2% and 16.6% for a woman without family history, and 8.6% and 24.4% for a woman with a first-degree family history of breast cancer.

CONCLUSIONS

The PRS stratifies breast cancer risk in women both with and without a family history of breast cancer. The observed level of risk discrimination could inform targeted screening and prevention strategies. Further discrimination may be achievable through combining the PRS with lifestyle/environmental factors, although these were not considered in this report.

Common germline polymorphisms associated with breast cancer-specific survival

INTRODUCTION

Previous studies have identified common germline variants nominally associated with breast cancer survival. These associations have not been widely replicated in further studies. The purpose of this study was to evaluate the association of previously reported SNPs with breast cancer-specific survival using data from a pooled analysis of eight breast cancer survival genome-wide association studies (GWAS) from the Breast Cancer Association Consortium.

METHODS

A literature review was conducted of all previously published associations between common germline variants and three survival outcomes: breast cancer-specific survival, overall survival and disease-free survival. All associations that reached the nominal significance level of P value <0.05 were included. Single nucleotide polymorphisms that had been previously reported as nominally associated with at least one survival outcome were evaluated in the pooled analysis of over 37,000 breast cancer cases for association with breast cancer-specific survival. Previous associations were evaluated using a one-sided test based on the reported direction of effect.

RESULTS

Fifty-six variants from 45 previous publications were evaluated in the meta-analysis. Fifty-four of these were evaluated in the full set of 37,954 breast cancer cases with 2,900 events and the two additional variants were evaluated in a reduced sample size of 30,000 samples in order to ensure independence from the previously published studies. Five variants reached nominal significance (P <0.05) in the pooled GWAS data compared to 2.8 expected under the null hypothesis. Seven additional variants were associated (P <0.05) with ER-positive disease.

CONCLUSIONS

Although no variants reached genome-wide significance (P <5 x 10(-8)), these results suggest that there is some evidence of association between candidate common germline variants and breast cancer prognosis. Larger studies from multinational collaborations are necessary to increase the power to detect associations, between common variants and prognosis, at more stringent significance levels.

Genome-wide association analysis of more than 120,000 individuals identifies 15 new susceptibility loci for breast cancer

Genome-wide association studies (GWAS) and large-scale replication studies have identified common variants in 79 loci associated with breast cancer, explaining ∼14% of the familial risk of the disease. To identify new susceptibility loci, we performed a meta-analysis of 11 GWAS, comprising 15,748 breast cancer cases and 18,084 controls together with 46,785 cases and 42,892 controls from 41 studies genotyped on a 211,155-marker custom array (iCOGS). Analyses were restricted to women of European ancestry. We generated genotypes for more than 11 million SNPs by imputation using the 1000 Genomes Project reference panel, and we identified 15 new loci associated with breast cancer at P < 5 × 10(-8). Combining association analysis with ChIP-seq chromatin binding data in mammary cell lines and ChIA-PET chromatin interaction data from ENCODE, we identified likely target genes in two regions: SETBP1 at 18q12.3 and RNF115 and PDZK1 at 1q21.1. One association appears to be driven by an amino acid substitution encoded in EXO1.

Investigation of gene-environment interactions between 47 newly identified breast cancer susceptibility loci and environmental risk factors

A large genotyping project within the Breast Cancer Association Consortium (BCAC) recently identified 41 associations between single nucleotide polymorphisms (SNPs) and overall breast cancer (BC) risk. We investigated whether the effects of these 41 SNPs, as well as six SNPs associated with estrogen receptor (ER) negative BC risk are modified by 13 environmental risk factors for BC. Data from 22 studies participating in BCAC were pooled, comprising up to 26,633 cases and 30,119 controls. Interactions between SNPs and environmental factors were evaluated using an empirical Bayes-type shrinkage estimator. Six SNPs showed interactions with associated p-values (pint ) <1.1 × 10(-3) . None of the observed interactions was significant after accounting for multiple testing. The Bayesian False Discovery Probability was used to rank the findings, which indicated three interactions as being noteworthy at 1% prior probability of interaction. SNP rs6828523 was associated with increased ER-negative BC risk in women ≥170 cm (OR = 1.22, p = 0.017), but inversely associated with ER-negative BC risk in women <160 cm (OR = 0.83, p = 0.039, pint = 1.9 × 10(-4) ). The inverse association between rs4808801 and overall BC risk was stronger for women who had had four or more pregnancies (OR = 0.85, p = 2.0 × 10(-4) ), and absent in women who had had just one (OR = 0.96, p = 0.19, pint = 6.1 × 10(-4) ). SNP rs11242675 was inversely associated with overall BC risk in never/former smokers (OR = 0.93, p = 2.8 × 10(-5) ), but no association was observed in current smokers (OR = 1.07, p = 0.14, pint = 3.4 × 10(-4) ). In conclusion, recently identified BC susceptibility loci are not strongly modified by established risk factors and the observed potential interactions require confirmation in independent studies.

Inherited variants in the inner centromere protein (INCENP) gene of the chromosomal passenger complex contribute to the susceptibility of ER-negative breast cancer

The chromosomal passenger complex (CPC) plays a pivotal role in the regulation of cell division. Therefore, inherited CPC variability could influence tumor development. The present candidate gene approach investigates the relationship between single nucleotide polymorphisms (SNPs) in genes encoding key CPC components and breast cancer risk. Fifteen SNPs in four CPC genes (INCENP, AURKB, BIRC5 and CDCA8) were genotyped in 88 911 European women from 39 case-control studies of the Breast Cancer Association Consortium. Possible associations were investigated in fixed-effects meta-analyses. The synonymous SNP rs1675126 in exon 7 of INCENP was associated with overall breast cancer risk [per A allele odds ratio (OR) 0.95, 95% confidence interval (CI) 0.92-0.98, P = 0.007] and particularly with estrogen receptor (ER)-negative breast tumors (per A allele OR 0.89, 95% CI 0.83-0.95, P = 0.0005). SNPs not directly genotyped were imputed based on 1000 Genomes. The SNPs rs1047739 in the 3' untranslated region and rs144045115 downstream of INCENP showed the strongest association signals for overall (per T allele OR 1.03, 95% CI 1.00-1.06, P = 0.0009) and ER-negative breast cancer risk (per A allele OR 1.06, 95% CI 1.02-1.10, P = 0.0002). Two genotyped SNPs in BIRC5 were associated with familial breast cancer risk (top SNP rs2071214: per G allele OR 1.12, 95% CI 1.04-1.21, P = 0.002). The data suggest that INCENP in the CPC pathway contributes to ER-negative breast cancer susceptibility in the European population. In spite of a modest contribution of CPC-inherited variants to the total burden of sporadic and familial breast cancer, their potential as novel targets for breast cancer treatment should be further investigated.

Fine-scale mapping of the 5q11.2 breast cancer locus reveals at least three independent risk variants regulating MAP3K1

Genome-wide association studies (GWASs) have revealed SNP rs889312 on 5q11.2 to be associated with breast cancer risk in women of European ancestry. In an attempt to identify the biologically relevant variants, we analyzed 909 genetic variants across 5q11.2 in 103,991 breast cancer individuals and control individuals from 52 studies in the Breast Cancer Association Consortium. Multiple logistic regression analyses identified three independent risk signals: the strongest associations were with 15 correlated variants (iCHAV1), where the minor allele of the best candidate, rs62355902, associated with significantly increased risks of both estrogen-receptor-positive (ER(+): odds ratio [OR] = 1.24, 95% confidence interval [CI] = 1.21-1.27, ptrend = 5.7 × 10(-44)) and estrogen-receptor-negative (ER(-): OR = 1.10, 95% CI = 1.05-1.15, ptrend = 3.0 × 10(-4)) tumors. After adjustment for rs62355902, we found evidence of association of a further 173 variants (iCHAV2) containing three subsets with a range of effects (the strongest was rs113317823 [pcond = 1.61 × 10(-5)]) and five variants composing iCHAV3 (lead rs11949391; ER(+): OR = 0.90, 95% CI = 0.87-0.93, pcond = 1.4 × 10(-4)). Twenty-six percent of the prioritized candidate variants coincided with four putative regulatory elements that interact with the MAP3K1 promoter through chromatin looping and affect MAP3K1 promoter activity. Functional analysis indicated that the cancer risk alleles of four candidates (rs74345699 and rs62355900 [iCHAV1], rs16886397 [iCHAV2a], and rs17432750 [iCHAV3]) increased MAP3K1 transcriptional activity. Chromatin immunoprecipitation analysis revealed diminished GATA3 binding to the minor (cancer-protective) allele of rs17432750, indicating a mechanism for its action. We propose that the cancer risk alleles act to increase MAP3K1 expression in vivo and might promote breast cancer cell survival.

Identification and characterization of novel associations in the CASP8/ALS2CR12 region on chromosome 2 with breast cancer risk

Previous studies have suggested that polymorphisms in CASP8 on chromosome 2 are associated with breast cancer risk. To clarify the role of CASP8 in breast cancer susceptibility, we carried out dense genotyping of this region in the Breast Cancer Association Consortium (BCAC). Single-nucleotide polymorphisms (SNPs) spanning a 1 Mb region around CASP8 were genotyped in 46 450 breast cancer cases and 42 600 controls of European origin from 41 studies participating in the BCAC as part of a custom genotyping array experiment (iCOGS). Missing genotypes and SNPs were imputed and, after quality exclusions, 501 typed and 1232 imputed SNPs were included in logistic regression models adjusting for study and ancestry principal components. The SNPs retained in the final model were investigated further in data from nine genome-wide association studies (GWAS) comprising in total 10 052 case and 12 575 control subjects. The most significant association signal observed in European subjects was for the imputed intronic SNP rs1830298 in ALS2CR12 (telomeric to CASP8), with per allele odds ratio and 95% confidence interval [OR (95% confidence interval, CI)] for the minor allele of 1.05 (1.03-1.07), P = 1 × 10(-5). Three additional independent signals from intronic SNPs were identified, in CASP8 (rs36043647), ALS2CR11 (rs59278883) and CFLAR (rs7558475). The association with rs1830298 was replicated in the imputed results from the combined GWAS (P = 3 × 10(-6)), yielding a combined OR (95% CI) of 1.06 (1.04-1.08), P = 1 × 10(-9). Analyses of gene expression associations in peripheral blood and normal breast tissue indicate that CASP8 might be the target gene, suggesting a mechanism involving apoptosis.

Refined histopathological predictors of BRCA1 and BRCA2 mutation status: a large-scale analysis of breast cancer characteristics from the BCAC, CIMBA, and ENIGMA consortia

INTRODUCTION

The distribution of histopathological features of invasive breast tumors in BRCA1 or BRCA2 germline mutation carriers differs from that of individuals with no known mutation. Histopathological features thus have utility for mutation prediction, including statistical modeling to assess pathogenicity of BRCA1 or BRCA2 variants of uncertain clinical significance. We analyzed large pathology datasets accrued by the Consortium of Investigators of Modifiers of BRCA1/2 (CIMBA) and the Breast Cancer Association Consortium (BCAC) to reassess histopathological predictors of BRCA1 and BRCA2 mutation status, and provide robust likelihood ratio (LR) estimates for statistical modeling.

METHODS

Selection criteria for study/center inclusion were estrogen receptor (ER) status or grade data available for invasive breast cancer diagnosed younger than 70 years. The dataset included 4,477 BRCA1 mutation carriers, 2,565 BRCA2 mutation carriers, and 47,565 BCAC breast cancer cases. Country-stratified estimates of the likelihood of mutation status by histopathological markers were derived using a Mantel-Haenszel approach.

RESULTS

ER-positive phenotype negatively predicted BRCA1 mutation status, irrespective of grade (LRs from 0.08 to 0.90). ER-negative grade 3 histopathology was more predictive of positive BRCA1 mutation status in women 50 years or older (LR = 4.13 (3.70 to 4.62)) versus younger than 50 years (LR = 3.16 (2.96 to 3.37)). For BRCA2, ER-positive grade 3 phenotype modestly predicted positive mutation status irrespective of age (LR = 1.7-fold), whereas ER-negative grade 3 features modestly predicted positive mutation status at 50 years or older (LR = 1.54 (1.27 to 1.88)). Triple-negative tumor status was highly predictive of BRCA1 mutation status for women younger than 50 years (LR = 3.73 (3.43 to 4.05)) and 50 years or older (LR = 4.41 (3.86 to 5.04)), and modestly predictive of positive BRCA2 mutation status in women 50 years or older (LR = 1.79 (1.42 to 2.24)).

CONCLUSIONS

These results refine likelihood-ratio estimates for predicting BRCA1 and BRCA2 mutation status by using commonly measured histopathological features. Age at diagnosis is an important variable for most analyses, and grade is more informative than ER status for BRCA2 mutation carrier prediction. The estimates will improve BRCA1 and BRCA2 variant classification and inform patient mutation testing and clinical management.

Common non-synonymous SNPs associated with breast cancer susceptibility: findings from the Breast Cancer Association Consortium

Candidate variant association studies have been largely unsuccessful in identifying common breast cancer susceptibility variants, although most studies have been underpowered to detect associations of a realistic magnitude. We assessed 41 common non-synonymous single-nucleotide polymorphisms (nsSNPs) for which evidence of association with breast cancer risk had been previously reported. Case-control data were combined from 38 studies of white European women (46 450 cases and 42 600 controls) and analyzed using unconditional logistic regression. Strong evidence of association was observed for three nsSNPs: ATXN7-K264R at 3p21 [rs1053338, per allele OR = 1.07, 95% confidence interval (CI) = 1.04-1.10, P = 2.9 × 10(-6)], AKAP9-M463I at 7q21 (rs6964587, OR = 1.05, 95% CI = 1.03-1.07, P = 1.7 × 10(-6)) and NEK10-L513S at 3p24 (rs10510592, OR = 1.10, 95% CI = 1.07-1.12, P = 5.1 × 10(-17)). The first two associations reached genome-wide statistical significance in a combined analysis of available data, including independent data from nine genome-wide association studies (GWASs): for ATXN7-K264R, OR = 1.07 (95% CI = 1.05-1.10, P = 1.0 × 10(-8)); for AKAP9-M463I, OR = 1.05 (95% CI = 1.04-1.07, P = 2.0 × 10(-10)). Further analysis of other common variants in these two regions suggested that intronic SNPs nearby are more strongly associated with disease risk. We have thus identified a novel susceptibility locus at 3p21, and confirmed previous suggestive evidence that rs6964587 at 7q21 is associated with risk. The third locus, rs10510592, is located in an established breast cancer susceptibility region; the association was substantially attenuated after adjustment for the known GWAS hit. Thus, each of the associated nsSNPs is likely to be a marker for another, non-coding, variant causally related to breast cancer risk. Further fine-mapping and functional studies are required to identify the underlying risk-modifying variants and the genes through which they act.

Genetic variation in mitotic regulatory pathway genes is associated with breast tumor grade

Mitotic index is an important component of histologic grade and has an etiologic role in breast tumorigenesis. Several small candidate gene studies have reported associations between variation in mitotic genes and breast cancer risk. We measured associations between 2156 single nucleotide polymorphisms (SNPs) from 194 mitotic genes and breast cancer risk, overall and by histologic grade, in the Breast Cancer Association Consortium (BCAC) iCOGS study (n = 39 067 cases; n = 42 106 controls). SNPs in TACC2 [rs17550038: odds ratio (OR) = 1.24, 95% confidence interval (CI) 1.16-1.33, P = 4.2 × 10(-10)) and EIF3H (rs799890: OR = 1.07, 95% CI 1.04-1.11, P = 8.7 × 10(-6)) were significantly associated with risk of low-grade breast cancer. The TACC2 signal was retained (rs17550038: OR = 1.15, 95% CI 1.07-1.23, P = 7.9 × 10(-5)) after adjustment for breast cancer risk SNPs in the nearby FGFR2 gene, suggesting that TACC2 is a novel, independent genome-wide significant genetic risk locus for low-grade breast cancer. While no SNPs were individually associated with high-grade disease, a pathway-level gene set analysis showed that variation across the 194 mitotic genes was associated with high-grade breast cancer risk (P = 2.1 × 10(-3)). These observations will provide insight into the contribution of mitotic defects to histological grade and the etiology of breast cancer.

MicroRNA related polymorphisms and breast cancer risk

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88-0.96), rs1052532 (OR 0.97; 95% CI: 0.95-0.99), rs10719 (OR 0.97; 95% CI: 0.94-0.99), rs4687554 (OR 0.97; 95% CI: 0.95-0.99, and rs3134615 (OR 1.03; 95% CI: 1.01-1.05) located in the 3' UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.

Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche

Age at menarche is a marker of timing of puberty in females. It varies widely between individuals, is a heritable trait and is associated with risks for obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality. Studies of rare human disorders of puberty and animal models point to a complex hypothalamic-pituitary-hormonal regulation, but the mechanisms that determine pubertal timing and underlie its links to disease risk remain unclear. Here, using genome-wide and custom-genotyping arrays in up to 182,416 women of European descent from 57 studies, we found robust evidence (P < 5 × 10(-8)) for 123 signals at 106 genomic loci associated with age at menarche. Many loci were associated with other pubertal traits in both sexes, and there was substantial overlap with genes implicated in body mass index and various diseases, including rare disorders of puberty. Menarche signals were enriched in imprinted regions, with three loci (DLK1-WDR25, MKRN3-MAGEL2 and KCNK9) demonstrating parent-of-origin-specific associations concordant with known parental expression patterns. Pathway analyses implicated nuclear hormone receptors, particularly retinoic acid and γ-aminobutyric acid-B2 receptor signalling, among novel mechanisms that regulate pubertal timing in humans. Our findings suggest a genetic architecture involving at least hundreds of common variants in the coordinated timing of the pubertal transition.