A meta-analysis of genome-wide association studies of breast cancer identifies two novel susceptibility loci at 6q14 and 20q11

Pathway-based discovery of genetic interactions in breast cancer

Breast cancer is the second largest cause of cancer death among U.S. women and the leading cause of cancer death among women worldwide. Genome-wide association studies (GWAS) have identified several genetic variants associated with susceptibility to breast cancer, but these still explain less than half of the estimated genetic contribution to the disease. Combinations of variants (i.e. genetic interactions) may play an important role in breast cancer susceptibility. However, due to a lack of statistical power, the current tests for genetic interactions from GWAS data mainly leverage prior knowledge to focus on small sets of genes or SNPs that are known to have an association with breast cancer. Thus, many genetic interactions, particularly among novel variants, remain understudied. Reverse-genetic interaction screens in model organisms have shown that genetic interactions frequently cluster into highly structured motifs, where members of the same pathway share similar patterns of genetic interactions. Based on this key observation, we recently developed a method called BridGE to search for such structured motifs in genetic networks derived from GWAS studies and identify pathway-level genetic interactions in human populations. We applied BridGE to six independent breast cancer cohorts and identified significant pathway-level interactions in five cohorts. Joint analysis across all five cohorts revealed a high confidence consensus set of genetic interactions with support in multiple cohorts. The discovered interactions implicated the glutathione conjugation, vitamin D receptor, purine metabolism, mitotic prometaphase, and steroid hormone biosynthesis pathways as major modifiers of breast cancer risk. Notably, while many of the pathways identified by BridGE show clear relevance to breast cancer, variants in these pathways had not been previously discovered by traditional single variant association tests, or single pathway enrichment analysis that does not consider SNP-SNP interactions.

Mechanisms of Breast Cancer in Shift Workers: DNA Methylation in Five Core Circadian Genes in Nurses Working Night Shifts

Shift work has been suggested to be associated with breast cancer risk, and circadian disruption in shift workers is hypothesized as one of the mechanisms of increased cancer risk. There is, however, insufficient molecular evidence supporting this hypothesis. Using the quantitative methodology of pyrosequencing, epigenetic changes in 5-methyl cytosine (5mC) in five circadian genes CLOCK, BMAL1, CRY1, PER1 and PER2 in female nurses working night shift work (278 breast cancer cases, 280 controls) were analyzed. In breast cancer cases, a medium exposure to night work was associated with increased methylation levels of the CLOCK (p=0.050), BMAL1 (p=0.001) and CRY1 (p=0.040) genes, compared with controls. Within the cases, analysis of the effects of shift work on the methylation patterns showed that methylation of CRY1 was lower in those who had worked night shift and had a high exposure (p=0.006) compared with cases that had worked only days. For cases with a medium exposure to night work, an increase in BMAL1 (p=0.003) and PER1 (p=0.035) methylation was observed compared with day working (unexposed) cases. The methylation levels of the five core circadian genes were also analyzed in relation to the estrogen and progesterone receptors status of the tumors in the cases, and no correlations were observed. Furthermore, nineteen polymorphisms in the five circadian genes were assessed for their effects on the methylation levels of the respective genes, but no associations were found. In summary, our data suggest that epigenetic regulation of CLOCK, BMAL1, CRY1 and PER1 may contribute to breast cancer in shift workers.

Identification of candidate protective variants for common diseases and evaluation of their protective potential

BACKGROUND

Human polymorphisms with derived alleles that are protective against disease may provide powerful translational opportunities. Here we report a method to identify such candidate polymorphisms and apply it to common non-synonymous SNPs (nsSNPs) associated with common diseases. Our study also sought to establish which of the identified protective nsSNPs show evidence of positive selection, taking this as indirect evidence that the protective variant has a beneficial effect on phenotype. Further, we performed an analysis to quantify the predicted effect of each protective variant on protein function/structure.

RESULTS

An initial analysis of eight SNPs previously identified as associated with age-related macular degeneration (AMD), revealed that two of them have a derived allele that is protective against developing the disease. One is in the complement component 2 gene (C2; E318D) and the other is in the complement factor B gene (CFB; R32Q). Then, combining genomewide ancestral allele information with known common disease-associated nsSNPs from the GWAS catalog, we found 32 additional SNPs which have a derived allele that is disease protective. Out of the total 34 identified candidate protective variants (CPVs), we found that 30 show stronger evidence of positive selection than the protective variant in lipoprotein lipase (LPL; S447X), which has already been translated into gene therapy. Furthermore, 11 of these CPVs have a higher probability of affecting protein structure than the lipoprotein lipase protective variant (LPL; S447X).

CONCLUSIONS

We identify 34 CPVs from the human genome. Diseases they confer protection against include, but are not limited to, type 2 diabetes, inflammatory bowel disease, age-related macular degeneration, multiple sclerosis and rheumatoid arthritis. We propose that those 30 CPVs with evidence of stronger positive selection than the LPL protective variant, may be considered as priority candidates for therapeutic approaches. The next step towards translation will require testing the hypotheses generated by our analyses, specifically whether the CPV arose from a gain-of-function or a loss-of-function mutation.

Genome-Wide Association Studies of Cancer in Diverse Populations

Genome-wide association studies (GWAS) of cancer have identified more than 700 risk loci, of which approximately 80% were first discovered in European ancestry populations, approximately 15% in East Asians, 3% in multiethnic scans, and less than 1% in African and Latin American populations. These percentages closely mirror the distribution of samples included in the discovery phase of cancer GWAS to date (84% European, 11% East Asian, 4% African, and 1% Latin American ancestry). GWAS in non-European ancestry populations have provided insight into ancestry-specific variation in cancer and have pointed to regions of susceptibility that are of particular importance in certain populations. Uncovering and characterizing cancer risk loci in diverse populations is critical for understanding underlying biological mechanisms and developing future genetic risk prediction models in non-European ancestry populations. New GWAS and continued collaborations will be required to eliminate population inequalities in the number of studies, sample sizes, and variant content on GWAS arrays, and to better align genetic research in cancer to the global distribution of race/ethnicity Cancer Epidemiol Biomarkers Prev; 27(4); 405-17. ©2018 AACRSee all articles in this CEBP Focus section, "Genome-Wide Association Studies in Cancer."

Genetic Variants in Epigenetic Pathways and Risks of Multiple Cancers in the GAME-ON Consortium

Background: Epigenetic disturbances are crucial in cancer initiation, potentially with pleiotropic effects, and may be influenced by the genetic background.Methods: In a subsets (ASSET) meta-analytic approach, we investigated associations of genetic variants related to epigenetic mechanisms with risks of breast, lung, colorectal, ovarian and prostate carcinomas using 51,724 cases and 52,001 controls. False discovery rate-corrected P values (q values < 0.05) were considered statistically significant.Results: Among 162,887 imputed or genotyped variants in 555 candidate genes, SNPs in eight genes were associated with risk of more than one cancer type. For example, variants in BABAM1 were confirmed as a susceptibility locus for squamous cell lung, overall breast, estrogen receptor (ER)-negative breast, and overall prostate, and overall serous ovarian cancer; the most significant variant was rs4808076 [OR = 1.14; 95% confidence interval (CI) = 1.10-1.19; q = 6.87 × 10-5]. DPF1 rs12611084 was inversely associated with ER-negative breast, endometrioid ovarian, and overall and aggressive prostate cancer risk (OR = 0.93; 95% CI = 0.91-0.96; q = 0.005). Variants in L3MBTL3 were associated with colorectal, overall breast, ER-negative breast, clear cell ovarian, and overall and aggressive prostate cancer risk (e.g., rs9388766: OR = 1.06; 95% CI = 1.03-1.08; q = 0.02). Variants in TET2 were significantly associated with overall breast, overall prostate, overall ovarian, and endometrioid ovarian cancer risk, with rs62331150 showing bidirectional effects. Analyses of subpathways did not reveal gene subsets that contributed disproportionately to susceptibility.Conclusions: Functional and correlative studies are now needed to elucidate the potential links between germline genotype, epigenetic function, and cancer etiology.Impact: This approach provides novel insight into possible pleiotropic effects of genes involved in epigenetic processes. Cancer Epidemiol Biomarkers Prev; 26(6); 816-25. ©2017 AACR.

Evaluating genetic variants associated with breast cancer risk in high and moderate-penetrance genes in Asians

Over the past 20 years, high-penetrance pathogenic mutations in genes BRCA1, BRCA2, TP53, PTEN, STK11 and CDH1 and moderate-penetrance mutations in genes CHEK2, ATM, BRIP1, PALB2, RAD51C, RAD50 and NBN have been identified for breast cancer. In this study, we investigated whether there are additional variants in these 13 genes associated with breast cancer among women of Asian ancestry. We analyzed up to 654 single nucleotide polymorphisms (SNPs) from 6269 cases and 6624 controls of Asian descent included in the Breast Cancer Association Consortium (BCAC), and up to 236 SNPs from 5794 cases and 5529 controls included in the Shanghai Breast Cancer Genetics Study (SBCGS). We found three missense variants with minor allele frequency (MAF) <0.05: rs80358978 (Gly2508Ser), rs80359065 (Lys2729Asn) and rs11571653 (Met784Val) in the BRCA2 gene, showing statistically significant associations with breast cancer risk, with P-values of 1.2 × 10-4, 1.0 × 10-3 and 5.0 × 10-3, respectively. In addition, we found four low-frequency variants (rs8176085, rs799923, rs8176173 and rs8176258) in the BRCA1 gene, one common variant in the CHEK2 gene (rs9620817), and one common variant in the PALB2 gene (rs13330119) associated with breast cancer risk at P < 0.01. Our study identified several new risk variants in BRCA1, BRCA2, CHEK2, and PALB2 genes in relation to breast cancer risk in Asian women. These results provide further insights that, in addition to the high/moderate penetrance mutations, other low-penetrance variants in these genes may also contribute to breast cancer risk.

Characterizing Genetic Susceptibility to Breast Cancer in Women of African Ancestry

Background: Genome-wide association studies have identified approximately 100 common genetic variants associated with breast cancer risk, the majority of which were discovered in women of European ancestry. Because of different patterns of linkage disequilibrium, many of these genetic markers may not represent signals in populations of African ancestry.Methods: We tested 74 breast cancer risk variants and conducted fine-mapping of these susceptibility regions in 6,522 breast cancer cases and 7,643 controls of African ancestry from three genetic consortia (AABC, AMBER, and ROOT).Results: Fifty-four of the 74 variants (73%) were found to have ORs that were directionally consistent with those previously reported, of which 12 were nominally statistically significant (P < 0.05). Through fine-mapping, in six regions (3p24, 12p11, 14q13, 16q12/FTO, 16q23, 19p13), we observed seven markers that better represent the underlying risk variant for overall breast cancer or breast cancer subtypes, whereas in another two regions (11q13, 16q12/TOX3), we identified suggestive evidence of signals that are independent of the reported index variant. Overlapping chromatin features and regulatory elements suggest that many of the risk alleles lie in regions with biological functionality.Conclusions: Through fine-mapping of known susceptibility regions, we have revealed alleles that better characterize breast cancer risk in women of African ancestry.Impact: The risk alleles identified represent genetic markers for modeling and stratifying breast cancer risk in women of African ancestry. Cancer Epidemiol Biomarkers Prev; 26(7); 1016-26. ©2017 AACR.

Genomic Disparities in Breast Cancer Among Latinas

BACKGROUND

Breast cancer is the most common cancer diagnosed among Latinas in the United States and the leading cause of cancer-related death among this population. Latinas tend to be diagnosed at a later stage and have worse prognostic features than their non-Hispanic white counterparts. Genetic and genomic factors may contribute to observed breast cancer health disparities in Latinas.

METHODS

We provide a landscape of our current understanding and the existing gaps that need to be filled across the cancer prevention and control continuum.

RESULTS

We summarize available data on mutations in high and moderate penetrance genes for inherited risk of breast cancer and the associated literature on disparities in awareness of and uptake of genetic counseling and testing in Latina populations. We also discuss common genetic polymorphisms and risk of breast cancer in Latinas. In the treatment setting, we examine tumor genomics and pharmacogenomics in Latina patients with breast cancer.

CONCLUSIONS

As the US population continues to diversify, extending genetic and genomic research into this underserved and understudied population is critical. By understanding the risk of breast cancer among ethnically diverse populations, we will be better positioned to make treatment advancements for earlier stages of cancer, identify more effective and ideally less toxic treatment regimens, and increase rates of survival.

A somatic-mutational process recurrently duplicates germline susceptibility loci and tissue-specific super-enhancers in breast cancers

Somatic rearrangements contribute to the mutagenized landscape of cancer genomes. Here, we systematically interrogated rearrangements in 560 breast cancers by using a piecewise constant fitting approach. We identified 33 hotspots of large (>100 kb) tandem duplications, a mutational signature associated with homologous-recombination-repair deficiency. Notably, these tandem-duplication hotspots were enriched in breast cancer germline susceptibility loci (odds ratio (OR) = 4.28) and breast-specific 'super-enhancer' regulatory elements (OR = 3.54). These hotspots may be sites of selective susceptibility to double-strand-break damage due to high transcriptional activity or, through incrementally increasing copy number, may be sites of secondary selective pressure. The transcriptomic consequences ranged from strong individual oncogene effects to weak but quantifiable multigene expression effects. We thus present a somatic-rearrangement mutational process affecting coding sequences and noncoding regulatory elements and contributing a continuum of driver consequences, from modest to strong effects, thereby supporting a polygenic model of cancer development.

Cis-eQTL-based trans-ethnic meta-analysis reveals novel genes associated with breast cancer risk

Breast cancer is the most common solid organ malignancy and the most frequent cause of cancer death among women worldwide. Previous research has yielded insights into its genetic etiology, but there remains a gap in the understanding of genetic factors that contribute to risk, and particularly in the biological mechanisms by which genetic variation modulates risk. The National Cancer Institute's "Up for a Challenge" (U4C) competition provided an opportunity to further elucidate the genetic basis of the disease. Our group leveraged the seven datasets made available by the U4C organizers and data from the publicly available UK Biobank cohort to examine associations between imputed gene expression and breast cancer risk. In particular, we used reference datasets describing the breast tissue and whole blood transcriptomes to impute expression levels in breast cancer cases and controls. In trans-ethnic meta-analyses of U4C and UK Biobank data, we found significant associations between breast cancer risk and the expression of RCCD1 (joint p-value: 3.6x10-06) and DHODH (p-value: 7.1x10-06) in breast tissue, as well as a suggestive association for ANKLE1 (p-value: 9.3x10-05). Expression of RCCD1 in whole blood was also suggestively associated with disease risk (p-value: 1.2x10-05), as were expression of ACAP1 (p-value: 1.9x10-05) and LRRC25 (p-value: 5.2x10-05). While genome-wide association studies (GWAS) have implicated RCCD1 and ANKLE1 in breast cancer risk, they have not identified the remaining three genes. Among the genetic variants that contributed to the predicted expression of the five genes, we found 23 nominally (p-value < 0.05) associated with breast cancer risk, among which 15 are not in high linkage disequilibrium with risk variants previously identified by GWAS. In summary, we used a transcriptome-based approach to investigate the genetic underpinnings of breast carcinogenesis. This approach provided an avenue for deciphering the functional relevance of genes and genetic variants involved in breast cancer.

A functionally significant SNP in TP53 and breast cancer risk in African-American women

A coding region polymorphism exists in the TP53 gene (Pro47Ser; rs1800371) in individuals of African descent, which reduces p53 tumor suppressor function in a mouse model. It has been unclear whether this functionally significant polymorphism alters cancer risk in humans. This analysis included 6907 women with breast cancer and 7644 controls from the AMBER, ROOT, and AABC consortia. We used multivariable logistic regression to estimate associations between the TP53 Pro47Ser allele and overall breast cancer risk. Because polymorphisms in TP53 tend to be associated with cancer risk in pre-menopausal women, we also limited our analyses to this population in the AMBER and ROOT consortia, where menopausal status was known, and conducted a fixed effects meta-analysis. In an analysis of all women in the AMBER, ROOT, and AABC consortia, we found no evidence for association of the Pro47Ser variant with breast cancer risk. However, when we restricted our analysis to only pre-menopausal women from the AMBER and ROOT consortia, there was a per allele odds ratio of 1.72 (95% confidence interval 1.08-2.76; p-value = 0.023). Although the Pro47Ser variant was not associated with overall breast cancer risk, it may increase risk among pre-menopausal women of African ancestry. Following up on more studies in human populations may better elucidate the role of this variant in breast cancer etiology. However, because of the low frequency of the polymorphism in women of African ancestry, its impact at a population level may be minimal.

Using Breast Cancer Risk Associated Polymorphisms to Identify Women for Breast Cancer Chemoprevention

BACKGROUND

Breast cancer can be prevented with selective estrogen receptor modifiers (SERMs) and aromatase inhibitors (AIs). The US Preventive Services Task Force recommends that women with a 5-year breast cancer risk ≥3% consider chemoprevention for breast cancer. More than 70 single nucleotide polymorphisms (SNPs) have been associated with breast cancer. We sought to determine how to best integrate risk information from SNPs with other risk factors to risk stratify women for chemoprevention.

METHODS

We used the risk distribution among women ages 35-69 estimated by the Breast Cancer Surveillance Consortium (BCSC) risk model. We modeled the effect of adding 70 SNPs to the BCSC model and examined how this would affect how many women are reclassified above and below the threshold for chemoprevention.

RESULTS

We found that most of the benefit of SNP testing a population is achieved by testing a modest fraction of the population. For example, if women with a 5-year BCSC risk of >2.0% are tested (~21% of all women), ~75% of the benefit of testing all women (shifting women above or below 3% 5-year risk) would be derived. If women with a 5-year risk of >1.5% are tested (~36% of all women), ~90% of the benefit of testing all women would be derived.

CONCLUSION

SNP testing is effective for reclassification of women for chemoprevention, but is unlikely to reclassify women with <1.5% 5-year risk. These results can be used to implement an efficient two-step testing approach to identify high risk women who may benefit from chemoprevention.

Previous GWAS hits in relation to young-onset breast cancer

PURPOSE

Genome-wide association studies (GWAS) have identified dozens of single-nucleotide polymorphisms (SNPs) associated with breast cancer. Few studies focused on young-onset breast cancer, which exhibits etiologic and tumor-type differences from older-onset disease. Possible confounding by prenatal effects of the maternal genome has also not been considered.

METHODS

Using a family-based design for breast cancer before age 50, we assessed the relationship between breast cancer and 77 GWAS-identified breast cancer risk SNPs. We estimated relative risks (RR) for inherited and maternally mediated genetic effects. We also used published RR estimates to calculate genetic risk scores and model joint effects.

RESULTS

Seventeen of the candidate SNPs were nominally associated with young-onset breast cancer in our 1296 non-Hispanic white affected families (uncorrected p value <0.05). Top-ranked SNPs included rs3803662-A (TOX3, RR = 1.39; p = 7.0 × 10-6), rs12662670-G (ESR1, RR = 1.56; p = 5.7 × 10-4), rs2981579-A (FGFR2, RR = 1.24; p = 0.002), and rs999737-G (RAD51B, RR = 1.37; p = 0.003). No maternally mediated effects were found. A risk score based on all 77 SNPs indicated that their overall relationship to young-onset breast cancer risk was more than additive (additive-fit p = 2.2 × 10-7) and consistent with a multiplicative joint effect (multiplicative-fit p = 0.27). With the multiplicative formulation, the case sister's genetic risk score exceeded that of her unaffected sister in 59% of families.

CONCLUSIONS

The results of this family-based study indicate that no effects of previously identified risk SNPs were explained by prenatal effects of maternal variants. Many of the known breast cancer risk variants were associated with young-onset breast cancer, with evidence that TOX3, ESR1, FGFR2, and RAD51B are important for young-onset disease.

Telomere structure and maintenance gene variants and risk of five cancer types

Telomeres cap chromosome ends, protecting them from degradation, double-strand breaks, and end-to-end fusions. Telomeres are maintained by telomerase, a reverse transcriptase encoded by TERT, and an RNA template encoded by TERC. Loci in the TERT and adjoining CLPTM1L region are associated with risk of multiple cancers. We therefore investigated associations between variants in 22 telomere structure and maintenance gene regions and colorectal, breast, prostate, ovarian, and lung cancer risk. We performed subset-based meta-analyses of 204,993 directly-measured and imputed SNPs among 61,851 cancer cases and 74,457 controls of European descent. Independent associations for SNP minor alleles were identified using sequential conditional analysis (with gene-level p value cutoffs ≤3.08 × 10-5 ). Of the thirteen independent SNPs observed to be associated with cancer risk, novel findings were observed for seven loci. Across the DCLRE1B region, rs974494 and rs12144215 were inversely associated with prostate and lung cancers, and colorectal, breast, and prostate cancers, respectively. Across the TERC region, rs75316749 was positively associated with colorectal, breast, ovarian, and lung cancers. Across the DCLRE1B region, rs974404 and rs12144215 were inversely associated with prostate and lung cancers, and colorectal, breast, and prostate cancers, respectively. Near POT1, rs116895242 was inversely associated with colorectal, ovarian, and lung cancers, and RTEL1 rs34978822 was inversely associated with prostate and lung cancers. The complex association patterns in telomere-related genes across cancer types may provide insight into mechanisms through which telomere dysfunction in different tissues influences cancer risk.

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, cis-acting variants represent an important source of phenotypic variation. Consequently, 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, P = 5.6x10-6). rs11099601 lies in a 135 kb linkage disequilibrium block containing several genes, including, HELQ, encoding the protein HEL308 a DNA dependant ATPase and DNA Helicase involved in DNA repair, MRPS18C encoding the Mitochondrial Ribosomal Protein S18C and FAM175A (ABRAXAS), encoding a 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 HELQ (P = 8.28x10-14), MRPS18C (P = 1.94x10-27) and FAM175A (P = 3.83x10-3), explaining about 20%, 14% and 1%, respectively of the variance inexpression of these genes in breast carcinomas.

Association between breast cancer genetic susceptibility variants and terminal duct lobular unit involution of the breast

Terminal duct lobular units (TDLUs) are the predominant source of future breast cancers, and lack of TDLU involution (higher TDLU counts, higher acini count per TDLU and the product of the two) is a breast cancer risk factor. Numerous breast cancer susceptibility single nucleotide polymorphisms (SNPs) have been identified, but whether they are associated with TDLU involution is unknown. In a pooled analysis of 872 women from two studies, we investigated 62 established breast cancer SNPs and relationships with TDLU involution. Poisson regression models with robust variance were used to calculate adjusted per-allele relative risks (with the non-breast cancer risk allele as the referent) and 95% confidence intervals between TDLU measures and each SNP. All statistical tests were two-sided; P < 0.05 was considered statistically significant. Overall, 36 SNPs (58.1%) were related to higher TDLU counts although this was not statistically significant (p = 0.25). Six of the 62 SNPs (9.7%) were nominally associated with at least one TDLU measure: rs616488 (PEX14), rs11242675 (FOXQ1) and rs6001930 (MKL1) were associated with higher TDLU count (p = 0.047, 0.045 and 0.031, respectively); rs1353747 (PDE4D) and rs6472903 (8q21.11) were associated with higher acini count per TDLU (p = 0.007 and 0.027, respectively); and rs1353747 (PDE4D) and rs204247 (RANBP9) were associated with the product of TDLU and acini counts (p = 0.024 and 0.017, respectively). Our findings suggest breast cancer SNPs may not strongly influence TDLU involution. Agnostic genome-wide association studies of TDLU involution may provide new insights on its biologic underpinnings and breast cancer susceptibility.

Deep targeted sequencing of 12 breast cancer susceptibility regions in 4611 women across four different ethnicities

Background

Although genome-wide association studies (GWASs) have identified thousands of disease susceptibility regions, the underlying causal mechanism in these regions is not fully known. It is likely that the GWAS signal originates from one or many as yet unidentified causal variants.

Methods

Using next-generation sequencing, we characterized 12 breast cancer susceptibility regions identified by GWASs in 2288 breast cancer cases and 2323 controls across four populations of African American, European, Japanese, and Hispanic ancestry.

Results

After genotype calling and quality control, we identified 137,530 single-nucleotide variants (SNVs); of those, 87.2 % had a minor allele frequency (MAF) <0.005. For SNVs with MAF >0.005, we calculated the smallest number of SNVs needed to obtain a posterior probability set (PPS) such that there is 90 % probability that the causal SNV is included. We found that the PPS for two regions, 2q35 and 11q13, contained less than 5 % of the original SNVs, dramatically decreasing the number of potentially causal SNVs. However, we did not find strong evidence supporting a causal role for any individual SNV. In addition, there were no significant gene-based rare SNV associations after correcting for multiple testing.

Conclusions

This study illustrates some of the challenges faced in fine-mapping studies in the post-GWAS era, most importantly the large sample sizes needed to identify rare-variant associations or to distinguish the effects of strongly correlated common SNVs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13058-016-0772-7) contains supplementary material, which is available to authorized users.

The genetics of breast cancer risk in the post-genome era: thoughts on study design to move past BRCA and towards clinical relevance

More than 12 % of women will be diagnosed with breast cancer in their lifetime. Although there have been tremendous advances in elucidating genetic risk factors underlying both familial and sporadic breast cancer, much of the genetic contribution to breast cancer etiology remains unknown. The discovery of BRCA1 and BRCA2 over 20 years ago remains the seminal event in the field and has paved the way for the discovery of other high-penetrance susceptibility genes by linkage analysis. The advent of genome-wide association studies made possible the next wave of discoveries, in which over 80 low-penetrance and moderate-penetrance variants were identified. Although these studies were highly successful at discovering variants associated with both familial and sporadic breast cancer, the variants identified to date explain only 50 % of the heritability of breast cancer. In this review, we look back at the investigative strategies that have led to our current understanding of breast cancer genetics, consider the challenges of performing association studies in heterogeneous complex diseases such as breast cancer, and look ahead toward the types of study designs that may lead to the identification of the genetic variation accounting for the remaining missing heritability.

Modifiers of breast and ovarian cancer risks for BRCA1 and BRCA2 mutation carriers

Pathogenic mutations in BRCA1 and BRCA2 are associated with high risks of breast and ovarian cancer. However, penetrance estimates for mutation carriers have been found to vary substantially between studies, and the observed differences in risk are consistent with the hypothesis that genetic and environmental factors modify cancer risks for women with these mutations. Direct evidence that this is the case has emerged in the past decade, through large-scale international collaborative efforts. Here, we describe the methodological challenges in the identification and characterisation of these risk-modifying factors, review the latest evidence on genetic and lifestyle/hormonal risk factors that modify breast and ovarian cancer risks for women with BRCA1 and BRCA2 mutations and outline the implications of these findings for cancer risk prediction. We also review the unresolved issues in this area of research and identify strategies of clinical implementation so that women with BRCA1 and BRCA2 mutations are no longer counselled on the basis of 'average' risk estimates.

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.

DcR3 gene polymorphisms are associated with sporadic breast infiltrating ductal carcinoma in Northeast Chinese women

Decoy Receptor 3 (DcR3), also called TNFRSF6β, is a member of the tumor necrosis factor receptor superfamily and is a soluble receptor for FasL. DcR3 is overexpressed in cancers and contributes to tumorigenesis through immune suppression and promotion of angiogenesis. We found that DcR3 is overexpressed in breast infiltrating ductal carcinoma (IDC) cells as compared with normal controls. We also conducted a case-control study analyzing associations of DcR3 polymorphisms with breast IDC risk. Subjects included 531 females with breast IDC and 592 age-matched healthy controls. Four DcR3 single nucleotide polymorphism loci with minor frequencies of more than 5% (rs3208008, rs41309931, rs2297441 and rs1291207) were genotyped using polymerase chain reaction restriction fragment length polymorphism and sequencing. Our results revealed significant differences in rs41309931genotypes and alleles (P < 0.01). Based on Haploview software analysis, the haplotype block A_rs3208008 G_rs41309931 G_rs2297441 A_rs1291207 exhibited the highest frequency, but, haplotype blocks A_rs3208008 T_rs41309931 G_rs2297441 A_rs1291207 and C_rs3208008 G_rs41309931 G_rs2297441 A_rs1291207 were associated with breast IDC risk. This study also detected associations between DcR3 gene polymorphisms and the clinicopathological features of breast IDC, including lymph node metastasis and C-erbB2, P53, estrogen receptor and progesterone receptor status. These data indicate that DcR3 gene polymorphisms are associated with sporadic breast IDC risk in Northeast Chinese females.

An investigation of the association of genetic susceptibility risk with somatic mutation burden in breast cancer

BACKGROUND

Genome-wide association studies have reported nearly 100 common germline susceptibility loci associated with the risk for breast cancer. Tumour sequencing studies have characterised somatic mutation profiles in breast cancer patients. The relationship between breast cancer susceptibility loci and somatic mutation patterns in breast cancer remains largely unexplored.

METHODS

We used single-nucleotide polymorphism (SNP) genotyping array data and tumour exome sequencing data available from 638 breast cancer patients of European ancestry from The Cancer Genome Atlas (TCGA) project. We analysed both genotype data and, when necessary, imputed genotypes for 90 known breast cancer susceptibility loci. We performed linear regression models to investigate possible associations between germline risk variants with total somatic mutation count (TSMC), as well as specific mutation types. We examined individual SNP genotypes, as well as a multi-SNP polygenic risk score (PRS). Models were statistically adjusted for age at diagnosis, stage, oestrogen-receptor (ER) and progesterone-receptor (PR) status of breast cancer. We also performed stratified analyses by ER and PR status.

RESULTS

We observed a significant inverse association (P=8.75 × 10(-6); FDR=0.001) between the risk allele in rs2588809 of the gene RAD51B and TSMC across all breast cancer patients, for both ER(+) and ER(-) tumours. This association was also evident for different types of mutations. The PRS analysis for all patients, with or without rs2588809, showed a significant inverse association (P=0.01 and 0.04, respectively) with TSMC. This inverse association was significant in ER(+) patients with the ER(+)-specific PRS (P=0.02), but not among ER(-) patients for the ER(-)-specific PRS (P=0.39).

CONCLUSIONS

We observed an inverse association between common germline risk variants and TSMC, which, if confirmed, could provide new insights into how germline variation informs our understanding of somatic mutation patterns in breast cancer.

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.

Noncoding somatic and inherited single-nucleotide variants converge to promote ESR1 expression in breast cancer

Sustained expression of the oestrogen receptor alpha (ESR1) drives two-thirds of breast cancer and defines the ESR1-positive subtype. ESR1 engages enhancers upon oestrogen stimulation to establish an oncogenic expression program¹. Somatic copy number alterations involving the ESR1 gene occur in approximately 1% of ESR1-positive breast cancers^2–5, implying that other mechanisms underlie the persistent expression of ESR1. We report the significant enrichment of somatic mutations within the set of regulatory elements (SRE) regulating ESR1 in 7% of ESR1-positive breast cancers. These mutations regulate ESR1 expression by modulating transcription factor binding to the DNA. The SRE includes a recurrently mutated enhancer whose activity is also affected by a functional inherited single nucleotide variant (SNV) rs9383590 that accounts for several breast cancer risk-loci. Our work highlights the importance of considering the combinatorial activity of regulatory elements as a single unit to delineate the impact of noncoding genetic alterations on single genes in cancer.

Breast cancer risk prediction using a clinical risk model and polygenic risk score

Breast cancer risk assessment can inform the use of screening and prevention modalities. We investigated the performance of the Breast Cancer Surveillance Consortium (BCSC) risk model in combination with a polygenic risk score (PRS) comprised of 83 single nucleotide polymorphisms identified from genome-wide association studies. We conducted a nested case-control study of 486 cases and 495 matched controls within a screening cohort. The PRS was calculated using a Bayesian approach. The contributions of the PRS and variables in the BCSC model to breast cancer risk were tested using conditional logistic regression. Discriminatory accuracy of the models was compared using the area under the receiver operating characteristic curve (AUROC). Increasing quartiles of the PRS were positively associated with breast cancer risk, with OR 2.54 (95 % CI 1.69-3.82) for breast cancer in the highest versus lowest quartile. In a multivariable model, the PRS, family history, and breast density remained strong risk factors. The AUROC of the PRS was 0.60 (95 % CI 0.57-0.64), and an Asian-specific PRS had AUROC 0.64 (95 % CI 0.53-0.74). A combined model including the BCSC risk factors and PRS had better discrimination than the BCSC model (AUROC 0.65 versus 0.62, p = 0.01). The BCSC-PRS model classified 18 % of cases as high-risk (5-year risk ≥3 %), compared with 7 % using the BCSC model. The PRS improved discrimination of the BCSC risk model and classified more cases as high-risk. Further consideration of the PRS's role in decision-making around screening and prevention strategies is merited.

Fine-Mapping of the 1p11.2 Breast Cancer Susceptibility Locus

The Cancer Genetic Markers of Susceptibility genome-wide association study (GWAS) originally identified a single nucleotide polymorphism (SNP) rs11249433 at 1p11.2 associated with breast cancer risk. To fine-map this locus, we genotyped 92 SNPs in a 900kb region (120,505,799-121,481,132) flanking rs11249433 in 45,276 breast cancer cases and 48,998 controls of European, Asian and African ancestry from 50 studies in the Breast Cancer Association Consortium. Genotyping was done using iCOGS, a custom-built array. Due to the complicated nature of the region on chr1p11.2: 120,300,000-120,505,798, that lies near the centromere and contains seven duplicated genomic segments, we restricted analyses to 429 SNPs excluding the duplicated regions (42 genotyped and 387 imputed). Per-allelic associations with breast cancer risk were estimated using logistic regression models adjusting for study and ancestry-specific principal components. The strongest association observed was with the original identified index SNP rs11249433 (minor allele frequency (MAF) 0.402; per-allele odds ratio (OR) = 1.10, 95% confidence interval (CI) 1.08-1.13, P = 1.49 x 10-21). The association for rs11249433 was limited to ER-positive breast cancers (test for heterogeneity P≤8.41 x 10-5). Additional analyses by other tumor characteristics showed stronger associations with moderately/well differentiated tumors and tumors of lobular histology. Although no significant eQTL associations were observed, in silico analyses showed that rs11249433 was located in a region that is likely a weak enhancer/promoter. Fine-mapping analysis of the 1p11.2 breast cancer susceptibility locus confirms this region to be limited to risk to cancers that are ER-positive.

The PLCO Biorepository: Creating, Maintaining, and Administering a Unique Biospecimen Resource

Inclusion of biospecimens in population-based studies is an integral part of understanding disease etiology, identifying biomarkers and developing prevention and treatment strategies. The Prostate, Lung, Colorectal, and Ovarian (PLCO) cancer screening trial collected, processed and stored biospecimens from participants to create a biorepository of specimens which serves as a useful resource for a broad research community. PLCO collected blood samples from consented screening arm participants at six screening rounds and a buccal sample from consented control arm participants. In addition, formalin-fixed paraffin embedded tumor tissue specimens were collected for participants in both arms for selected cancer sites. Collection of biospecimens at multiple timepoints (i.e. serial samples) and prior to cancer diagnosis, paired with rich epidemiologic and screening data, makes the PLCO collection of biospecimens a uniquely valuable resource. As such, access to the PLCO biorepository is granted to investigators by a rigorous scientific review process and guided by a steering committee which is responsible for developing and implementing the biospecimen use policies. Here, we describe the procedures for biospecimen collection, processing, storage, requisition, and distribution, as well as data management employed in PLCO. We also provide examples of how the biospecimens have been used to advance cancer research and describe relevant lessons learned to help inform cohorts wishing to add or modify biospecimen collection.

Genome-wide association study in East Asians identifies two novel breast cancer susceptibility loci

Breast cancer is one of the most common malignancies among women worldwide. Genetic factors have been shown to play an important role in breast cancer aetiology. We conducted a two-stage genome-wide association study (GWAS) including 14 224 cases and 14 829 controls of East Asian women to search for novel genetic susceptibility loci for breast cancer. Single nucleotide polymorphisms (SNPs) in two loci were found to be associated with breast cancer risk at the genome-wide significance level. The first locus, represented by rs12118297 at 1p22.3 (near the LMO4 gene), was associated with breast cancer risk with odds ratio (OR) and (95% confidence interval (CI)) of 0.91 (0.88-0.94) and a P-value of 4.48 × 10^{-  8} This association was replicated in another study, DRIVE GAME-ON Consortium, including 16 003 cases and 41 335 controls of European ancestry (OR = 0.95, 95% CI = 0.91-0.99, P-value = 0.019). The second locus, rs16992204 at 21q22.12 (near the LINC00160 gene), was associated with breast cancer risk with OR (95% CI) of 1.13 (1.07-1.18) and a P-value of 4.63 × 10 ^{-  8} The risk allele frequency for this SNP is zero in European-ancestry populations in 1000 Genomes Project and thus its association with breast cancer risk cannot be assessed in DRIVE GAME-ON Consortium. Functional annotation using the ENCODE data indicates that rs12118297 might be located in a repressed element and locus 21q22.12 may affect breast cancer risk through regulating LINC00160 expressions and interaction with oestrogen receptor signalling. Our findings provide additional insights into the genetics of breast cancer.

A Review of Whole-Exome Sequencing Efforts Toward Hereditary Breast Cancer Susceptibility Gene Discovery

Inherited genetic risk factors contribute toward breast cancer (BC) onset. BC risk variants can be divided into three categories of penetrance (high, moderate, and low) that reflect the probability of developing the disease. Traditional BC susceptibility gene discovery approaches that searched for high- and moderate-risk variants in familial BC cases have had limited success; to date, these risk variants explain only ∼30% of familial BC cases. Next-generation sequencing technologies can be used to search for novel high and moderate BC risk variants, and this manuscript reviews 12 familial BC whole-exome sequencing efforts. Study design, filtering strategies, and segregation and validation analyses are discussed. Overall, only a modest number of novel BC risk genes were identified, and 90% and 97% of the exome-sequenced families and cases, respectively, had no BC risk variants reported. It is important to learn from these studies and consider alternate strategies in order to make further advances. The discovery of new BC susceptibility genes is critical for improved risk assessment and to provide insight toward disease mechanisms for the development of more effective therapies.

An exome-wide analysis of low frequency and rare variants in relation to risk of breast cancer in African American Women: the AMBER Consortium

A large percentage of breast cancer heritability remains unaccounted for, and most of the known susceptibility loci have been established in European and Asian populations. Rare variants may contribute to the unexplained heritability of this disease, including in women of African ancestry (AA). We conducted an exome-wide analysis of rare variants in relation to risk of overall and subtype-specific breast cancer in the African American Breast Cancer Epidemiology and Risk (AMBER) Consortium, which includes data from four large studies of AA women. Genotyping on the Illumina Human Exome Beadchip yielded data for 170 812 SNPs and 8287 subjects: 3629 cases (1093 estrogen receptor negative (ER-), 1968 ER+, 568 ER unknown) and 4658 controls, the largest exome chip study to date for AA breast cancer. Pooled gene-based association analyses were performed using the unified optimal sequence kernel association test (SKAT-O) for variants with minor allele frequency (MAF) ≤ 5%. In addition, each variant with MAF >0.5% was tested for association using logistic regression. There were no significant associations with overall breast cancer. However, a novel gene, FBXL22 (P = 8.2×10(-6)), and a gene previously identified in GWAS of European ancestry populations, PDE4D (P = 1.2×10(-6)), were significantly associated with ER- breast cancer after correction for multiple testing. Cases with the associated rare variants were also negative for progesterone and human epidermal growth factor receptors-thus, triple-negative cancer. Replication is required to confirm these gene-level associations, which are based on very small counts at extremely rare SNPs.

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.

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.

Genes associated with histopathologic features of triple negative breast tumors predict molecular subtypes

Distinct subtypes of triple negative (TN) breast cancer have been identified by tumor expression profiling. However, little is known about the relationship between histopathologic features of TN tumors, which reflect aspects of both tumor behavior and tumor microenvironment, and molecular TN subtypes. The histopathologic features of TN tumors were assessed by central review and 593 TN tumors were subjected to whole genome expression profiling using the Illumina Whole Genome DASL array. TN molecular subtypes were defined based on gene expression data associated with histopathologic features of TN tumors. Gene expression analysis yielded signatures for four TN subtypes (basal-like, androgen receptor positive, immune, and stromal) consistent with previous studies. Expression analysis also identified genes significantly associated with the 12 histological features of TN tumors. Development of signatures using these markers of histopathological features resulted in six distinct TN subtype signatures, including an additional basal-like and stromal signature. The additional basal-like subtype was distinguished by elevated expression of cell motility and glucose metabolism genes and reduced expression of immune signaling genes, whereas the additional stromal subtype was distinguished by elevated expression of immunomodulatory pathway genes. Histopathologic features that reflect heterogeneity in tumor architecture, cell structure, and tumor microenvironment are related to TN subtype. Accounting for histopathologic features in the development of gene expression signatures, six major subtypes of TN breast cancer were identified.

Analysis of functional germline variants in APOBEC3 and driver genes on breast cancer risk in Moroccan study population

BACKGROUND

Breast cancer (BC) is the most prevalent cancer in women and a major public health problem in Morocco. Several Moroccan studies have focused on studying this disease, but more are needed, especially at the genetic and molecular levels. Therefore, we investigated the potential association of several functional germline variants in the genes commonly mutated in sporadic breast cancer.

METHODS

In this case-control study, we examined 36 single nucleotide polymorphisms (SNPs) in 13 genes (APOBEC3A, APOBEC3B, ARID1B, ATR, MAP3K1, MLL2, MLL3, NCOR1, RUNX1, SF3B1, SMAD4, TBX3, TTN), which were located in the core promoter, 5'-and 3'UTR or which were nonsynonymous SNPs to assess their potential association with inherited predisposition to breast cancer development. Additionally, we identified a ~29.5-kb deletion polymorphism between APOBEC3A and APOBEC3B and explored possible associations with BC. A total of 226 Moroccan breast cancer cases and 200 matched healthy controls were included in this study.

RESULTS

The analysis showed that12 SNPs in 8 driver genes, 4 SNPs in APOBEC3B gene and 1 SNP in APOBEC3A gene were associated with BC risk and/or clinical outcome at P ≤ 0.05 level. RUNX1_rs8130963 (odds ratio (OR) = 2.25; 95 % CI 1.42-3.56; P = 0.0005; dominant model), TBX3_rs8853 (OR = 2.04; 95 % CI 1.38-3.01; P = 0.0003; dominant model), TBX3_rs1061651 (OR= 2.14; 95 % CI1.43-3.18; P = 0.0002; dominant model), TTN_rs12465459 (OR = 2.02; 95 % confidence interval 1.33-3.07; P = 0.0009; dominant model), were the most significantly associated SNPs with BC risk. A strong association with clinical outcome were detected for the genes SMAD4 _rs3819122 with tumor size (OR = 0.45; 95 % CI 0.25-0.82; P = 0.009) and TTN_rs2244492 with estrogen receptor (OR = 0.45; 95 % CI 0.25-0.82; P = 0.009).

CONCLUSION

Our results suggest that genetic variations in driver and APOBEC3 genes were associated with the risk of BC and may have impact on clinical outcome. However, the reported association between the deletion polymorphism and BC risk was not confirmed in the Moroccan population. These preliminary findings require replication in larger studies.

Genetic predisposition to ductal carcinoma in situ of the breast

BACKGROUND

Ductal carcinoma in situ (DCIS) is a non-invasive form of breast cancer. It is often associated with invasive ductal carcinoma (IDC), and is considered to be a non-obligate precursor of IDC. It is not clear to what extent these two forms of cancer share low-risk susceptibility loci, or whether there are differences in the strength of association for shared loci.

METHODS

To identify genetic polymorphisms that predispose to DCIS, we pooled data from 38 studies comprising 5,067 cases of DCIS, 24,584 cases of IDC and 37,467 controls, all genotyped using the iCOGS chip.

RESULTS

Most (67 %) of the 76 known breast cancer predisposition loci showed an association with DCIS in the same direction as previously reported for invasive breast cancer. Case-only analysis showed no evidence for differences between associations for IDC and DCIS after considering multiple testing. Analysis by estrogen receptor (ER) status confirmed that loci associated with ER positive IDC were also associated with ER positive DCIS. Analysis of DCIS by grade suggested that two independent SNPs at 11q13.3 near CCND1 were specific to low/intermediate grade DCIS (rs75915166, rs554219). These associations with grade remained after adjusting for ER status and were also found in IDC. We found no novel DCIS-specific loci at a genome wide significance level of P < 5.0x10(-8).

CONCLUSION

In conclusion, this study provides the strongest evidence to date of a shared genetic susceptibility for IDC and DCIS. Studies with larger numbers of DCIS are needed to determine if IDC or DCIS specific loci exist.

GENESIS: a French national resource to study the missing heritability of breast cancer

Background

Less than 20 % of familial breast cancer patients who undergo genetic testing for BRCA1 and BRCA2 carry a pathogenic mutation in one of these two genes. The GENESIS (GENE SISter) study was designed to identify new breast cancer susceptibility genes in women attending cancer genetics clinics and with no BRCA1/2 mutation.

Methods

The study involved the French national network of family cancer clinics. It was based on enrichment in genetic factors of the recruited population through case selection relying on familial criteria, but also on the consideration of environmental factors and endophenotypes like mammary density or tumor characteristics to assess potential genetic heterogeneity. One of the initial aims of GENESIS was to recruit affected sibpairs. Siblings were eligible when index cases and at least one affected sister were diagnosed with infiltrating mammary or ductal adenocarcinoma, with no BRCA1/2 mutation. In addition, unrelated controls and unaffected sisters were recruited. The enrolment of patients, their relatives and their controls, the collection of the clinical, epidemiological, familial and biological data were centralized by a coordinating center.

Results

Inclusion of participants started in February 2007 and ended in December 2013. A total of 1721 index cases, 826 affected sisters, 599 unaffected sisters and 1419 controls were included. 98 % of participants completed the epidemiological questionnaire, 97 % provided a blood sample, and 76 % were able to provide mammograms. Index cases were on average 59 years old at inclusion, were born in 1950, and were 49.7 years of age at breast cancer diagnosis. The mean age at diagnosis of affected sisters was slightly higher (51.4 years). The representativeness of the control group was verified.

Conclusions

The size of the study, the availability of biological specimens and the clinical data collection together with the detailed and complete epidemiological questionnaire make this a unique national resource for investigation of the missing heritability of breast cancer, by taking into account environmental and life style factors and stratifying data on endophenotypes to decrease genetic heterogeneity.

A Cross-Cancer Genetic Association Analysis of the DNA Repair and DNA Damage Signaling Pathways for Lung, Ovary, Prostate, Breast, and Colorectal Cancer

BACKGROUND

DNA damage is an established mediator of carcinogenesis, although genome-wide association studies (GWAS) have identified few significant loci. This cross-cancer site, pooled analysis was performed to increase the power to detect common variants of DNA repair genes associated with cancer susceptibility.

METHODS

We conducted a cross-cancer analysis of 60,297 single nucleotide polymorphisms, at 229 DNA repair gene regions, using data from the NCI Genetic Associations and Mechanisms in Oncology (GAME-ON) Network. Our analysis included data from 32 GWAS and 48,734 controls and 51,537 cases across five cancer sites (breast, colon, lung, ovary, and prostate). Because of the unavailability of individual data, data were analyzed at the aggregate level. Meta-analysis was performed using the Association analysis for SubSETs (ASSET) software. To test for genetic associations that might escape individual variant testing due to small effect sizes, pathway analysis of eight DNA repair pathways was performed using hierarchical modeling.

RESULTS

We identified three susceptibility DNA repair genes, RAD51B (P < 5.09 × 10(-6)), MSH5 (P < 5.09 × 10(-6)), and BRCA2 (P = 5.70 × 10(-6)). Hierarchical modeling identified several pleiotropic associations with cancer risk in the base excision repair, nucleotide excision repair, mismatch repair, and homologous recombination pathways.

CONCLUSIONS

Only three susceptibility loci were identified, which had all been previously reported. In contrast, hierarchical modeling identified several pleiotropic cancer risk associations in key DNA repair pathways.

IMPACT

Results suggest that many common variants in DNA repair genes are likely associated with cancer susceptibility through small effect sizes that do not meet stringent significance testing criteria.

Identification of shared and unique susceptibility pathways among cancers of the lung, breast, and prostate from genome-wide association studies and tissue-specific protein interactions

Results from genome-wide association studies (GWAS) have indicated that strong single-gene effects are the exception, not the rule, for most diseases. We assessed the joint effects of germline genetic variations through a pathway-based approach that considers the tissue-specific contexts of GWAS findings. From GWAS meta-analyses of lung cancer (12 160 cases/16 838 controls), breast cancer (15 748 cases/18 084 controls) and prostate cancer (14 160 cases/12 724 controls) in individuals of European ancestry, we determined the tissue-specific interaction networks of proteins expressed from genes that are likely to be affected by disease-associated variants. Reactome pathways exhibiting enrichment of proteins from each network were compared across the cancers. Our results show that pathways associated with all three cancers tend to be broad cellular processes required for growth and survival. Significant examples include the nerve growth factor (P = 7.86 × 10(-33)), epidermal growth factor (P = 1.18 × 10(-31)) and fibroblast growth factor (P = 2.47 × 10(-31)) signaling pathways. However, within these shared pathways, the genes that influence risk largely differ by cancer. Pathways found to be unique for a single cancer focus on more specific cellular functions, such as interleukin signaling in lung cancer (P = 1.69 × 10(-15)), apoptosis initiation by Bad in breast cancer (P = 3.14 × 10(-9)) and cellular responses to hypoxia in prostate cancer (P = 2.14 × 10(-9)). We present the largest comparative cross-cancer pathway analysis of GWAS to date. Our approach can also be applied to the study of inherited mechanisms underlying risk across multiple diseases in general.

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.

Cross Cancer Genomic Investigation of Inflammation Pathway for Five Common Cancers: Lung, Ovary, Prostate, Breast, and Colorectal Cancer

BACKGROUND

Inflammation has been hypothesized to increase the risk of cancer development as an initiator or promoter, yet no large-scale study of inherited variation across cancer sites has been conducted.

METHODS

We conducted a cross-cancer genomic analysis for the inflammation pathway based on 48 genome-wide association studies within the National Cancer Institute GAME-ON Network across five common cancer sites, with a total of 64 591 cancer patients and 74 467 control patients. Subset-based meta-analysis was used to account for possible disease heterogeneity, and hierarchical modeling was employed to estimate the effect of the subcomponents within the inflammation pathway. The network was visualized by enrichment map. All statistical tests were two-sided.

RESULTS

We identified three pleiotropic loci within the inflammation pathway, including one novel locus in Ch12q24 encoding SH2B3 (rs3184504), which reached GWAS significance with a P value of 1.78 x 10(-8), and it showed an association with lung cancer (P = 2.01 x 10(-6)), colorectal cancer (GECCO P = 6.72x10(-6); CORECT P = 3.32x10(-5)), and breast cancer (P = .009). We also identified five key subpathway components with genetic variants that are relevant for the risk of these five cancer sites: inflammatory response for colorectal cancer (P = .006), inflammation related cell cycle gene for lung cancer (P = 1.35x10(-6)), and activation of immune response for ovarian cancer (P = .009). In addition, sequence variations in immune system development played a role in breast cancer etiology (P = .001) and innate immune response was involved in the risk of both colorectal (P = .022) and ovarian cancer (P = .003).

CONCLUSIONS

Genetic variations in inflammation and its related subpathway components are keys to the development of lung, colorectal, ovary, and breast cancer, including SH2B3, which is associated with lung, colorectal, and breast cancer.

Modeling Linkage Disequilibrium Increases Accuracy of Polygenic Risk Scores

Polygenic risk scores have shown great promise in predicting complex disease risk and will become more accurate as training sample sizes increase. The standard approach for calculating risk scores involves linkage disequilibrium (LD)-based marker pruning and applying a p value threshold to association statistics, but this discards information and can reduce predictive accuracy. We introduce LDpred, a method that infers the posterior mean effect size of each marker by using a prior on effect sizes and LD information from an external reference panel. Theory and simulations show that LDpred outperforms the approach of pruning followed by thresholding, particularly at large sample sizes. Accordingly, predicted R(2) increased from 20.1% to 25.3% in a large schizophrenia dataset and from 9.8% to 12.0% in a large multiple sclerosis dataset. A similar relative improvement in accuracy was observed for three additional large disease datasets and for non-European schizophrenia samples. The advantage of LDpred over existing methods will grow as sample sizes increase.

Interaction between common breast cancer susceptibility variants, genetic ancestry, and nongenetic risk factors in Hispanic women

BACKGROUND

Most genetic variants associated with breast cancer risk have been discovered in women of European ancestry, and only a few genome-wide association studies (GWAS) have been conducted in minority groups. This research disparity persists in post-GWAS gene-environment interaction analyses. We tested the interaction between hormonal and lifestyle risk factors for breast cancer, and ten GWAS-identified SNPs among 2,107 Hispanic women with breast cancer and 2,587 unaffected controls, to gain insight into a previously reported gene by ancestry interaction in this population.

METHODS

We estimated genetic ancestry with a set of 104 ancestry-informative markers selected to discriminate between Indigenous American and European ancestry. We used logistic regression models to evaluate main effects and interactions.

RESULTS

We found that the rs13387042-2q35(G/A) SNP was associated with breast cancer risk only among postmenopausal women who never used hormone therapy [per A allele OR: 0.94 (95% confidence intervals, 0.74-1.20), 1.20 (0.94-1.53), and 1.49 (1.28-1.75) for current, former, and never hormone therapy users, respectively, Pinteraction 0.002] and premenopausal women who breastfed >12 months [OR: 1.01 (0.72-1.42), 1.19 (0.98-1.45), and 1.69 (1.26-2.26) for never, <12 months, and >12 months breastfeeding, respectively, Pinteraction 0.014].

CONCLUSIONS

The correlation between genetic ancestry, hormone replacement therapy use, and breastfeeding behavior partially explained a previously reported interaction between a breast cancer risk variant and genetic ancestry in Hispanic women.

IMPACT

These results highlight the importance of understanding the interplay between genetic ancestry, genetics, and nongenetic risk factors and their contribution to breast cancer risk.

Concurrent DNA Copy-Number Alterations and Mutations in Genes Related to Maintenance of Genome Stability in Uninvolved Mammary Glandular Tissue from Breast Cancer Patients

Somatic mosaicism for DNA copy-number alterations (SMC-CNAs) is defined as gain or loss of chromosomal segments in somatic cells within a single organism. As cells harboring SMC-CNAs can undergo clonal expansion, it has been proposed that SMC-CNAs may contribute to the predisposition of these cells to genetic disease including cancer. Herein, the gross genomic alterations (>500 kbp) were characterized in uninvolved mammary glandular tissue from 59 breast cancer patients and matched samples of primary tumors and lymph node metastases. Array-based comparative genomic hybridization showed 10% (6/59) of patients harbored one to 359 large SMC-CNAs (mean: 1,328 kbp; median: 961 kbp) in a substantial portion of glandular tissue cells, distal from the primary tumor site. SMC-CNAs were partially recurrent in tumors, albeit with considerable contribution of stochastic SMC-CNAs indicating genomic destabilization. Targeted resequencing of 301 known predisposition and somatic driver loci revealed mutations and rare variants in genes related to maintenance of genomic integrity: BRCA1 (p.Gln1756Profs*74, p.Arg504Cys), BRCA2 (p.Asn3124Ile), NCOR1 (p.Pro1570Glnfs*45), PALB2 (p.Ser500Pro), and TP53 (p.Arg306*). Co-occurrence of gross SMC-CNAs along with point mutations or rare variants in genes responsible for safeguarding genomic integrity highlights the temporal and spatial neoplastic potential of uninvolved glandular tissue in breast cancer patients.

Association of breast cancer risk loci with breast cancer survival

The survival of breast cancer patients is largely influenced by tumor characteristics, such as TNM stage, tumor grade and hormone receptor status. However, there is growing evidence that inherited genetic variation might affect the disease prognosis and response to treatment. Several lines of evidence suggest that alleles influencing breast cancer risk might also be associated with breast cancer survival. We examined the associations between 35 breast cancer susceptibility loci and the disease over-all survival (OS) in 10,255 breast cancer patients from the National Cancer Institute Breast and Prostate Cancer Cohort Consortium (BPC3) of which 1,379 died, including 754 of breast cancer. We also conducted a meta-analysis of almost 35,000 patients and 5,000 deaths, combining results from BPC3 and the Breast Cancer Association Consortium (BCAC) and performed in silico analyses of SNPs with significant associations. In BPC3, the C allele of LSP1-rs3817198 was significantly associated with improved OS (HRper-allele =0.70; 95% CI: 0.58-0.85; ptrend  = 2.84 × 10(-4) ; HRheterozygotes  = 0.71; 95% CI: 0.55-0.92; HRhomozygotes  = 0.48; 95% CI: 0.31-0.76; p2DF  = 1.45 × 10(-3) ). In silico, the C allele of LSP1-rs3817198 was predicted to increase expression of the tumor suppressor cyclin-dependent kinase inhibitor 1C (CDKN1C). In the meta-analysis, TNRC9-rs3803662 was significantly associated with increased death hazard (HRMETA =1.09; 95% CI: 1.04-1.15; ptrend  = 6.6 × 10(-4) ; HRheterozygotes  = 0.96 95% CI: 0.90-1.03; HRhomozygotes  = 1.21; 95% CI: 1.09-1.35; p2DF =1.25 × 10(-4) ). In conclusion, we show that there is little overlap between the breast cancer risk single nucleotide polymorphisms (SNPs) identified so far and the SNPs associated with breast cancer prognosis, with the possible exceptions of LSP1-rs3817198 and TNRC9-rs3803662.

Genetic risk variants associated with in situ breast cancer

INTRODUCTION

Breast cancer in situ (BCIS) diagnoses, a precursor lesion for invasive breast cancer, comprise about 20 % of all breast cancers (BC) in countries with screening programs. Family history of BC is considered one of the strongest risk factors for BCIS.

METHODS

To evaluate the association of BC susceptibility loci with BCIS risk, we genotyped 39 single nucleotide polymorphisms (SNPs), associated with risk of invasive BC, in 1317 BCIS cases, 10,645 invasive BC cases, and 14,006 healthy controls in the National Cancer Institute's Breast and Prostate Cancer Cohort Consortium (BPC3). Using unconditional logistic regression models adjusted for age and study, we estimated the association of SNPs with BCIS using two different comparison groups: healthy controls and invasive BC subjects to investigate whether BCIS and BC share a common genetic profile.

RESULTS

We found that five SNPs (CDKN2BAS-rs1011970, FGFR2-rs3750817, FGFR2-rs2981582, TNRC9-rs3803662, 5p12-rs10941679) were significantly associated with BCIS risk (P value adjusted for multiple comparisons <0.0016). Comparing invasive BC and BCIS, the largest difference was for CDKN2BAS-rs1011970, which showed a positive association with BCIS (OR = 1.24, 95 % CI: 1.11-1.38, P = 1.27 x 10(-4)) and no association with invasive BC (OR = 1.03, 95 % CI: 0.99-1.07, P = 0.06), with a P value for case-case comparison of 0.006. Subgroup analyses investigating associations with ductal carcinoma in situ (DCIS) found similar associations, albeit less significant (OR = 1.25, 95 % CI: 1.09-1.42, P = 1.07 x 10(-3)). Additional risk analyses showed significant associations with invasive disease at the 0.05 level for 28 of the alleles and the OR estimates were consistent with those reported by other studies.

CONCLUSIONS

Our study adds to the knowledge that several of the known BC susceptibility loci are risk factors for both BCIS and invasive BC, with the possible exception of rs1011970, a putatively functional SNP situated in the CDKN2BAS gene that may be a specific BCIS susceptibility locus.

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.

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.