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

BRCA2 polymorphic stop codon K3326X and the risk of breast, prostate, and ovarian cancers

Background: The K3326X variant in BRCA2 (BRCA2*c.9976A&gt;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

Reproductive profiles and risk of breast cancer subtypes : a multi-center case-only study

Background: Previous studies have shown that reproductive factors are differentially associated with breast cancer (BC) risk by subtypes. The aim of this study was to investigate associations between reproductive factors and BC subtypes, and whether these vary by age at diagnosis. Methods: We used pooled data on tumor markers (estrogen and progesterone receptor, human epidermal growth factor receptor-2 (HER2)) and reproductive risk factors (parity, age at first full-time pregnancy (FFTP) and age at menarche) from 28,095 patients with invasive BC from 34 studies participating in the Breast Cancer Association Consortium (BCAC). In a case-only analysis, we used logistic regression to assess associations between reproductive factors and BC subtype compared to luminal A tumors as a reference. The interaction between age and parity in BC subtype risk was also tested, across all ages and, because age was modeled non-linearly, specifically at ages 35, 55 and 75 years. Results: Parous women were more likely to be diagnosed with triple negative BC (TNBC) than with luminal A BC, irrespective of age (OR for parity = 1.38, 95% CI 1.16-1.65, p = 0.0004; p for interaction with age = 0.076). Parous women were also more likely to be diagnosed with luminal and non-luminal HER2-like BCs and this effect was slightly more pronounced at an early age (p for interaction with age = 0.037 and 0. 030, respectively). For instance, women diagnosed at age 35 were 1.48 (CI 1.01-2.16) more likely to have luminal HER2-like BC than luminal A BC, while this association was not significant at age 75 (OR = 0.72, CI 0.45-1.14). While age at menarche was not significantly associated with BC subtype, increasing age at FFTP was non-linearly associated with TNBC relative to luminal A BC. An age at FFTP of 25 versus 20 years lowered the risk for TNBC (OR = 0.78, CI 0.70-0.88, p <0.0001), but this effect was not apparent at a later FFTP. Conclusions: Our main findings suggest that parity is associated with TNBC across all ages at BC diagnosis, whereas the association with luminal HER2-like BC was present only for early onset BC.Peer reviewe

A likelihood ratio approach for utilizing case-control data in the clinical classification of rare sequence variants: Application to BRCA1 and BRCA2

A large number of variants identified through clinical genetic testing in disease susceptibility genes are of uncertain significance (VUS). Following the recommendations of the American College of Medical Genetics and Genomics (ACMG) and Association for Molecular Pathology (AMP), the frequency in case-control datasets (PS4 criterion) can inform their interpretation. We present a novel case-control likelihood ratio-based method that incorporates gene-specific age-related penetrance. We demonstrate the utility of this method in the analysis of simulated and real datasets. In the analysis of simulated data, the likelihood ratio method was more powerful compared to other methods. Likelihood ratios were calculated for a case-control dataset of BRCA1 and BRCA2 variants from the Breast Cancer Association Consortium (BCAC) and compared with logistic regression results. A larger number of variants reached evidence in favor of pathogenicity, and a substantial number of variants had evidence against pathogenicity findings that would not have been reached using other case-control analysis methods. Our novel method provides greater power to classify rare variants compared with classical case-control methods. As an initiative from the ENIGMA Analytical Working Group, we provide user-friendly scripts and preformatted Excel calculators for implementation of the method for rare variants in BRCA1, BRCA2, and other high-risk genes with known penetrance

Association of breast cancer risk with genetic variants showing differential allelic expression: Identification of a novel breast cancer susceptibility locus at 4q21

There are significant inter-individual differences in the levels of gene expression. Through modulation of gene expression, $\textit{cis}$-acting variants represent an important source of phenotypic variation. Consequently, $\textit{cis}$-regulatory SNPs associated with differential allelic expression are functional candidates for further investigation as disease-causing variants. To investigate whether common variants associated with differential allelic expression were involved in breast cancer susceptibility, a list of genes was established on the basis of their involvement in cancer related pathways and/or mechanisms. Thereafter, using data from a genome-wide map of allelic expression associated SNPs, 313 genetic variants were selected and their association with breast cancer risk was then evaluated in 46,451 breast cancer cases and 42,599 controls of European ancestry ascertained from 41 studies participating in the Breast Cancer Association Consortium. The associations were evaluated with overall breast cancer risk and with estrogen receptor negative and positive disease. One novel breast cancer susceptibility locus on 4q21 (rs11099601) was identified (OR = 1.05, $\textit{P}$ = 5.6x10$^{-6}$). rs11099601 lies in a 135 kb linkage disequilibrium block containing several genes, including, $\textit{HELQ}$, encoding the protein HEL308 a DNA dependant ATPase and DNA Helicase involved in DNA repair, $\textit{MRPS18C}$ encoding the Mitochondrial Ribosomal Protein S18C and $\textit{FAM175A (ABRAXAS)}$, encoding a $\textit{BRCA1}$ BRCT domain-interacting protein involved in DNA damage response and double-strand break (DSB) repair. Expression QTL analysis in breast cancer tissue showed rs11099601 to be associated with $\textit{HELQ }$ ($\textit{P}$ = 8.28x10$^{-14}$), $\textit{MRPS18C}$ ($\textit{P}$ = 1.94x10$^{-27}$) and $\textit{FAM175A }$ ($\textit{P}$ = 3.83x10$^{-3}$), explaining about 20%, 14% and 1%, respectively of the variance inexpression of these genes in breast carcinomas.Information regarding funding can be found in the published article or the publisher's website. Funders include Cancer Research UK and the National Institute for Health Research

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.Peer reviewe

Evaluation of European-based polygenic risk score for breast cancer in Ashkenazi Jewish women in Israel

To date, most BC GWASs have been performed Background Polygenic risk score (PRS), calculated in individuals of European (EUR) ancestry, and based on genome-wide association studies (GWASs), the generalisation of EUR-based PRS to other can improve breast cancer (BC) risk assessment. populations is a major challenge. In this study, we examined the performance of EUR-based BC PRS models in Ashkenazi Jewish (AJ) women. Methods We generated PRSs based on data on EUR women from the Breast Cancer Association Consortium (BCAC). We tested the performance of the PRSs in a cohort of 2161 AJ women from Israel (1437 cases and 724 controls) from BCAC (BCAC cohort from Israel (BCAC-IL)). In addition, we tested the performance of these EUR-based BC PRSs, as well as the established 313-SNP EUR BC PRS, in an independent cohort of 181 AJ women from Hadassah Medical Center (HMC) in Israel. Results In the BCAC-IL cohort, the highest OR per 1 SD was 1.56 (±0.09). The OR for AJ women at the top 10% of the PRS distribution compared with the middle quintile was 2.10 (±0.24). In the HMC cohort, the OR per 1 SD of the EUR-based PRS that performed best in the BCAC-IL cohort was 1.58±0.27. The OR per 1 SD of the commonly used 313-SNP BC PRS was 1.64 (±0.28). Conclusions Extant EUR GWAS data can be used for generating PRSs that identify AJ women with markedly elevated risk of BC and therefore hold promise for improving BC risk assessment in AJ women

An intergenic risk locus containing an enhancer deletion in 2q35 modulates breast cancer risk by deregulating IGFBP5 expression.

Breast cancer is the most diagnosed malignancy and the second leading cause of cancer mortality in females. Previous association studies have identified variants on 2q35 associated with the risk of breast cancer. To identify functional susceptibility loci for breast cancer, we interrogated the 2q35 gene desert for chromatin architecture and functional variation correlated with gene expression. We report a novel intergenic breast cancer risk locus containing an enhancer copy number variation (enCNV; deletion) located approximately 400Kb upstream to IGFBP5, which overlaps an intergenic ERα-bound enhancer that loops to the IGFBP5 promoter. The enCNV is correlated with modified ERα binding and monoallelic-repression of IGFBP5 following estrogen treatment. We investigated the association of enCNV genotype with breast cancer in 1,182 cases and 1,362 controls, and replicate our findings in an independent set of 62,533 cases and 60,966 controls from 41 case control studies and 11 GWAS. We report a dose-dependent inverse association of 2q35 enCNV genotype (percopy OR=0.68 95%CI 0.55-0.83, P=0.0002; replication OR=0.77 95%CI 0.73-0.82, P=2.1x10(-19)) and identify 13 additional linked variants (r(2)>0.8) in the 20Kb linkage block containing the enCNV (P=3.2x10(-15) - 5.6x10(-17)). These associations were independent of previously reported 2q35 variants, rs13387042/rs4442975 and rs16857609, and were stronger for ER-positive than ER-negative disease. Together, these results suggest that 2q35 breast cancer risk loci may be mediating their effect through IGFBP5