Novel genetic markers of breast cancer survival identified by a genome-wide association study

The relationship between eight GWAS-identified single-nucleotide polymorphisms and primary breast cancer outcomes

BACKGROUND

Several single-nucleotide polymorphisms (SNPs) associated with breast cancer risk have been identified through genome-wide association studies (GWAS). We investigated whether eight risk SNPs identified in GWAS were associated with breast cancer disease-free survival (DFS) and overall survival (OS) rates.

PATIENTS AND METHODS

A cohort of 739 white women with early-stage breast cancer was genotyped for eight GWAS-identified SNPs (rs2981582, rs1219648 [FGFR2], rs3803662, rs12443621, rs8051542 [TOX3], rs999737 [RAD51L1], rs6504950 [17q23], and rs4973768 [3p24]). Relationships between SNPs and breast cancer outcomes were evaluated using Cox proportional hazard regression models. The cumulative effects of SNPs on breast cancer outcomes were assessed by computing the number of at-risk genotypes.

RESULTS

At a median follow-up of 121 months (range: 188-231 months) for survivors, 237 deaths (32%) and 186 breast cancer events (25%) were identified among the 739 patients. After adjusting for age, clinical stage, and treatment, rs12443621 (16q12; p = .03) and rs6504950 (17q23; p = .008) were prognostic for OS but not DFS. A higher risk for death was also found in the multivariable analysis of patients harboring three or four at-risk genotypes of the GWAS SNPs compared to patients carrying two or less at-risk genotypes (hazard ratio: 1.60, 95% confidence interval: 1.23-2.24; p = .0008).

CONCLUSION

The study results suggest that previously identified breast cancer risk susceptibility loci, rs12443621 (16q12) and rs6504950 (17q23), may influence breast cancer prognosis or comorbid conditions associated with overall survival. The precise molecular mechanisms through which these risk SNPs, as well as others that were not included in the analysis, influence clinical outcomes remain to be determined.

Identification of inherited genetic variations influencing prognosis in early-onset breast cancer

Genome-Wide Association Studies (GWAS) have begun to investigate associations between inherited genetic variations and breast cancer prognosis. Here, we report our findings from a GWAS conducted in 536 patients with early-onset breast cancer aged 40 or less at diagnosis and with a mean follow-up period of 4.1 years (SD = 1.96). Patients were selected from the Prospective Study of Outcomes in Sporadic versus Hereditary breast cancer. A Bonferroni correction for multiple testing determined that a P value of 1.0 × 10(-7) was a statistically significant association signal. Following quality control, we identified 487,496 single nucleotide polymorphisms (SNP) for association tests in stage 1. In stage 2, 35 SNPs with the most significant associations were genotyped in 1,516 independent cases from the same early-onset cohort. In stage 2, 11 SNPs remained associated in the same direction (P ≤ 0.05). Fixed effects meta-analysis models identified one SNP associated at close to genome wide level of significance 556 kb upstream of the ARRDC3 locus [HR = 1.61; 95% confidence interval (CI), 1.33-1.96; P = 9.5 × 10(-7)]. Four further associations at or close to the PBX1, RORα, NTN1, and SYT6 loci also came close to genome-wide significance levels (P = 10(-6)). In the first ever GWAS for the identification of SNPs associated with prognosis in patients with early-onset breast cancer, we report a SNP upstream of the ARRDC3 locus as potentially associated with prognosis (median follow-up time for genotypes: CC = 4 years, CT = 3 years, and TT = 2.7 years; Wilcoxon rank-sum test CC vs. CT, P = 4 × 10(-4) and CT vs. TT, P = 0.76). Four further loci may also be associated with prognosis.

Heterozygous genetic variations of FOXP3 in Xp11.23 elevate breast cancer risk in Chinese population via skewed X-chromosome inactivation

FOXP3 (forkhead box P3: also known as IPEX, XPID) is not only a hallmark of immunosuppressive regulatory T cells (Tregs), but also an X-linked breast cancer suppressor gene expressed in tumor cells. A two-stage investigation was conducted in individuals from northern, southern and eastern China. Individuals carrying a FOXP3 rs2294021CT genotype showed about 1.5-fold increased risk of breast cancer compared with TT carriers. In a related biochemical assay, the rs2294021C allele was found to significantly enhance transcription activity, leading to higher mRNA levels of FOXP3 compared with T allele. Moreover, the number of Tregs and its corresponding interleukin-10 (IL-10) secretion were elevated whereas the proliferation of antitumor T cells was decreased in the C-allele carriers. The breast cancer oncogenes Her-2/ErbB2 and Skp2 were also found to be significantly inhibited in C-allele carriers. Moreover, skewed X-chromosome inactivation (SXCI) analysis showed that rs2294021CT carriers with SXCI showed higher risk than the homozygous carriers and CT carriers without SXCI, suggesting a possible interaction between the rs2294021CT genotype and SXCI. Taken together, these findings indicate that the rs2294021CT genotype may increase an individual's susceptibility to breast cancer by breaking the balance between Treg-mediated immune tolerance and FOXP3-controlled tumor-suppressive effect.

Germline DNA copy number aberrations identified as potential prognostic factors for breast cancer recurrence

Breast cancer recurrence (BCR) is a common treatment outcome despite curative-intent primary treatment of non-metastatic breast cancer. Currently used prognostic and predictive factors utilize tumor-based markers, and are not optimal determinants of risk of BCR. Germline-based copy number aberrations (CNAs) have not been evaluated as determinants of predisposition to experience BCR. In this study, we accessed germline DNA from 369 female breast cancer subjects who received curative-intent primary treatment following diagnosis. Of these, 155 experienced BCR and 214 did not, after a median duration of follow up after breast cancer diagnosis of 6.35 years (range = 0.60–21.78) and 8.60 years (range = 3.08–13.57), respectively. Whole genome CNA genotyping was performed on the Affymetrix SNP array 6.0 platform. CNAs were identified using the SNP-Fast Adaptive States Segmentation Technique 2 algorithm implemented in Nexus Copy Number 6.0. Six samples were removed due to poor quality scores, leaving 363 samples for further analysis. We identified 18,561 CNAs with ≥1 kb as a predefined cut-off for observed aberrations. Univariate survival analyses (log-rank tests) identified seven CNAs (two copy number gains and five copy neutral-loss of heterozygosities, CN-LOHs) showing significant differences (P<2.01×10⁻⁵) in recurrence-free survival (RFS) probabilities with and without CNAs.We also observed three additional but distinct CN-LOHs showing significant differences in RFS probabilities (P<2.86×10⁻⁵) when analyses were restricted to stratified cases (luminal A, n = 208) only. After adjusting for tumor stage and grade in multivariate analyses (Cox proportional hazards models), all the CNAs remained strongly associated with the phenotype of BCR. Of these, we confirmed three CNAs at 17q11.2, 11q13.1 and 6q24.1 in representative samples using independent genotyping platforms. Our results suggest further investigations on the potential use of germline DNA variations as prognostic markers in cancer-associated phenotypes.

Genetic links between post-reproductive lifespan and family size in Framingham

BACKGROUND AND OBJECTIVES

Is there a trade-off between children ever born (CEB) and post-reproductive lifespan in humans? Here, we report a comprehensive analysis of reproductive trade-offs in the Framingham Heart Study (FHS) dataset using phenotypic and genotypic correlations and a genome-wide association study (GWAS) to look for single-nucleotide polymorphisms (SNPs) that are related to the association between CEB and lifespan.

METHODOLOGY

We calculated the phenotypic and genetic correlations of lifespan with CEB for men and women in the Framingham dataset, and then performed a GWAS to search for SNPs that might affect the relationship between post-reproductive lifespan and CEB.

RESULTS

We found significant negative phenotypic correlations between CEB and lifespan in both women (rP = -0.133, P < 0.001) and men (rP = -0. 079, P = 0.036). The genetic correlation was large, highly significant and strongly negative in women (rG = -0.877, P = 0.009) in a model without covariates, but not in men (P = 0.777). The GWAS identified five SNPs associated with the relationship between CEB and post-reproductive lifespan in women; some are near genes that have been linked to cancer. None were identified in men.

CONCLUSIONS AND IMPLICATIONS

We identified several SNPs for which the relationship between CEB and post-reproductive lifespan differs by genotype in women in the FHS who were born between 1889 and 1958. That result was not robust to changes in the sample. Further studies on larger samples are needed to validate the antagonistic pleiotropy of these genes.

A common polymorphism near the ESR1 gene is associated with risk of breast cancer: evidence from a case-control study and a meta-analysis

Background

Genome-wide association studies have reported that a polymorphism near the estrogen receptor gene (ESR1) (rs2046210) is associated with a risk of breast cancer, with the A allele conferring an increased risk. However, considering the controversial results from more recent replicated studies, we conducted a case-control study in an independent Chinese Han population and a meta-analysis to clarify the association of this polymorphism with breast cancer risk.

Method and Findings

A hospital-based case-control study including 461 cases and 537 controls from a Chinese Han population was conducted initially, and this study showed that the rs2046210 A allele was significantly associated with breast cancer risk, with an OR of 1.32 (95% CI  = 1.10–1.59). Subsequently, a meta-analysis integrating the current study and previous publications with a total of 53,379 cases and 55,493 controls was performed to further confirm our findings. Similarly, a significant association between this polymorphism and breast cancer risk was also observed in the overall population especially among Asians, with ORs for per A allele of 1.14 (95% CI  = 1.10–1.18) in the overall population and 1.27 (95% CI  = 1.23–1.31) in the Asian population.

Conclusion

Our results provide strong evidence to support that the common polymorphism near the ESR1 gene, rs2046210, is associated with an increased risk of breast cancer in Asian and European populations but not in Africans, although the biological mechanisms need to be further investigated.

The SLC4A7 variant rs4973768 is associated with breast cancer risk: evidence from a case-control study and a meta-analysis

Recent genome-wide association study has identified a genetic variant rs4973768, located in 3'-UTR of solute carrier family 4, sodium bicarbonate cotransporter, member 7 (SLC4A7), was associated with increased risk of breast cancer (BC). However, several following replication studies cannot yield consistent results. We thus conducted a hospital-based case-control study including 485 patients and 514 controls, combined a meta-analysis including 108,632 cases and 135,818 controls to explore the relationship between this variant and BC risk. Our case-control study showed that rs4973768 was significantly associated with increased BC risk with the odds ratio (OR) of 1.29 (95 % confidence interval [CI]: 1.04-1.60) under the allelic model. In addition, the meta-analysis also indicated that the variant slightly increased the risk of BC with the pooled OR of the per-allele effect being 1.08 (95 % CI: 1.04-1.11) although with significant heterogeneity between studies. Stratified analyses showed that ethnicity, sample size, and study design may explain part of the heterogeneity. Moreover, the bioinformatics analysis suggested that this variant may influence the transcriptional capacity of SLC4A7. In summary, our results showed that the SLC4A7 variant, rs4973768, is associated with risk of BC although the underlying biologic mechanism warrants further studies.