A comprehensive evaluation of interaction between genetic variants and use of menopausal hormone therapy on mammographic density.

INTRODUCTION: Mammographic density is an established breast cancer risk factor with a strong genetic component and can be increased in women using menopausal hormone therapy (MHT). Here, we aimed to identify genetic variants that may modify the association between MHT use and mammographic density. METHODS: The study comprised 6,298 postmenopausal women from the Mayo Mammography Health Study and nine studies included in the Breast Cancer Association Consortium. We selected for evaluation 1327 single nucleotide polymorphisms (SNPs) showing the lowest P-values for interaction (P int) in a meta-analysis of genome-wide gene-environment interaction studies with MHT use on risk of breast cancer, 2541 SNPs in candidate genes (AKR1C4, CYP1A1-CYP1A2, CYP1B1, ESR2, PPARG, PRL, SULT1A1-SULT1A2 and TNF) and ten SNPs (AREG-rs10034692, PRDM6-rs186749, ESR1-rs12665607, ZNF365-rs10995190, 8p11.23-rs7816345, LSP1-rs3817198, IGF1-rs703556, 12q24-rs1265507, TMEM184B-rs7289126, and SGSM3-rs17001868) associated with mammographic density in genome-wide studies. We used multiple linear regression models adjusted for potential confounders to evaluate interactions between SNPs and current use of MHT on mammographic density. RESULTS: No significant interactions were identified after adjustment for multiple testing. The strongest SNP-MHT interaction (unadjusted P int <0.0004) was observed with rs9358531 6.5kb 5' of PRL. Furthermore, three SNPs in PLCG2 that had previously been shown to modify the association of MHT use with breast cancer risk were found to modify also the association of MHT use with mammographic density (unadjusted P int <0.002), but solely among cases (unadjusted P int SNP×MHT×case-status <0.02). CONCLUSIONS: The study identified potential interactions on mammographic density between current use of MHT and SNPs near PRL and in PLCG2, which require confirmation. Given the moderate size of the interactions observed, larger studies are needed to identify genetic modifiers of the association of MHT use with mammographic density.This is the final version of the article. It first appeared from BioMed Central via http://dx.doi.org/10.1186/s13058-015-0625-

Evidence that breast cancer risk at the 2q35 locus is mediated through IGFBP5 regulation.

GWAS have identified a breast cancer susceptibility locus on 2q35. Here we report the fine mapping of this locus using data from 101,943 subjects from 50 case-control studies. We genotype 276 SNPs using the 'iCOGS' genotyping array and impute genotypes for a further 1,284 using 1000 Genomes Project data. All but two, strongly correlated SNPs (rs4442975 G/T and rs6721996 G/A) are excluded as candidate causal variants at odds against >100:1. The best functional candidate, rs4442975, is associated with oestrogen receptor positive (ER+) disease with an odds ratio (OR) in Europeans of 0.85 (95% confidence interval=0.84-0.87; P=1.7 × 10(-43)) per t-allele. This SNP flanks a transcriptional enhancer that physically interacts with the promoter of IGFBP5 (encoding insulin-like growth factor-binding protein 5) and displays allele-specific gene expression, FOXA1 binding and chromatin looping. Evidence suggests that the g-allele confers increased breast cancer susceptibility through relative downregulation of IGFBP5, a gene with known roles in breast cell biology

Breast cancer risk variants at 6q25 display different phenotype associations and regulate ESR1, RMND1 and CCDC170.

We analyzed 3,872 common genetic variants across the ESR1 locus (encoding estrogen receptor α) in 118,816 subjects from three international consortia. We found evidence for at least five independent causal variants, each associated with different phenotype sets, including estrogen receptor (ER(+) or ER(-)) and human ERBB2 (HER2(+) or HER2(-)) tumor subtypes, mammographic density and tumor grade. The best candidate causal variants for ER(-) tumors lie in four separate enhancer elements, and their risk alleles reduce expression of ESR1, RMND1 and CCDC170, whereas the risk alleles of the strongest candidates for the remaining independent causal variant disrupt a silencer element and putatively increase ESR1 and RMND1 expression.This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/ng.352

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

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

This corrects the article DOI: 10.1038/ncomms5999

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 x 10(-20)), ER-negative BC (P = 1.1 x 10(-13)), BRCA1-associated BC (P = 7.7 x 10(-16)) and triple negative BC (P-diff = 2 x 10(-5)). Genotype-gene expression associations are identified for candidate target genes ANKLE1 (P = 2 x 10(-3)) and ABHD8 (PPeer reviewe

Is RNASEL:p.Glu265* a modifier of early-onset breast cancer risk for carriers of high-risk mutations?

Abstract Background Breast cancer risk for BRCA1 and BRCA2 pathogenic mutation carriers is modified by risk factors that cluster in families, including genetic modifiers of risk. We considered genetic modifiers of risk for carriers of high-risk mutations in other breast cancer susceptibility genes. Methods In a family known to carry the high-risk mutation PALB2:c.3113G>A (p.Trp1038*), whole-exome sequencing was performed on germline DNA from four affected women, three of whom were mutation carriers. Results RNASEL:p.Glu265* was identified in one of the PALB2 carriers who had two primary invasive breast cancer diagnoses before 50 years. Gene-panel testing of BRCA1, BRCA2, PALB2 and RNASEL in the Australian Breast Cancer Family Registry identified five carriers of RNASEL:p.Glu265* in 591 early onset breast cancer cases. Three of the five women (60%) carrying RNASEL:p.Glu265* also carried a pathogenic mutation in a breast cancer susceptibility gene compared with 30 carriers of pathogenic mutations in the 586 non-carriers of RNASEL:p.Glu265* (5%) (p < 0.002). Taqman genotyping demonstrated that the allele frequency of RNASEL:p.Glu265* was similar in affected and unaffected Australian women, consistent with other populations. Conclusion Our study suggests that RNASEL:p.Glu265* may be a genetic modifier of risk for early-onset breast cancer predisposition in carriers of high-risk mutations. Much larger case-case and case-control studies are warranted to test the association observed in this report

A comprehensive evaluation of interaction between genetic variants and use of menopausal hormone therapy on mammographic density

Abstract Introduction Mammographic density is an established breast cancer risk factor with a strong genetic component and can be increased in women using menopausal hormone therapy (MHT). Here, we aimed to identify genetic variants that may modify the association between MHT use and mammographic density. Methods The study comprised 6,298 postmenopausal women from the Mayo Mammography Health Study and nine studies included in the Breast Cancer Association Consortium. We selected for evaluation 1327 single nucleotide polymorphisms (SNPs) showing the lowest P-values for interaction (P int) in a meta-analysis of genome-wide gene-environment interaction studies with MHT use on risk of breast cancer, 2541 SNPs in candidate genes (AKR1C4, CYP1A1-CYP1A2, CYP1B1, ESR2, PPARG, PRL, SULT1A1-SULT1A2 and TNF) and ten SNPs (AREG-rs10034692, PRDM6-rs186749, ESR1-rs12665607, ZNF365-rs10995190, 8p11.23-rs7816345, LSP1-rs3817198, IGF1-rs703556, 12q24-rs1265507, TMEM184B-rs7289126, and SGSM3-rs17001868) associated with mammographic density in genome-wide studies. We used multiple linear regression models adjusted for potential confounders to evaluate interactions between SNPs and current use of MHT on mammographic density. Results No significant interactions were identified after adjustment for multiple testing. The strongest SNP-MHT interaction (unadjusted P int <0.0004) was observed with rs9358531 6.5kb 5′ of PRL. Furthermore, three SNPs in PLCG2 that had previously been shown to modify the association of MHT use with breast cancer risk were found to modify also the association of MHT use with mammographic density (unadjusted P int <0.002), but solely among cases (unadjusted P int SNP×MHT×case-status <0.02). Conclusions The study identified potential interactions on mammographic density between current use of MHT and SNPs near PRL and in PLCG2, which require confirmation. Given the moderate size of the interactions observed, larger studies are needed to identify genetic modifiers of the association of MHT use with mammographic density

Genome-wide association and transcriptome studies identify target genes and risk loci for breast cancer

Abstract Genome-wide association studies (GWAS) have identified more than 170 breast cancer susceptibility loci. Here we hypothesize that some risk-associated variants might act in non-breast tissues, specifically adipose tissue and immune cells from blood and spleen. Using expression quantitative trait loci (eQTL) reported in these tissues, we identify 26 previously unreported, likely target genes of overall breast cancer risk variants, and 17 for estrogen receptor (ER)-negative breast cancer, several with a known immune function. We determine the directional effect of gene expression on disease risk measured based on single and multiple eQTL. In addition, using a gene-based test of association that considers eQTL from multiple tissues, we identify seven (and four) regions with variants associated with overall (and ER-negative) breast cancer risk, which were not reported in previous GWAS. Further investigation of the function of the implicated genes in breast and immune cells may provide insights into the etiology of breast cancer