Carcinogenesis Advance Access originally published online on May 19, 2006
Carcinogenesis 2006 27(11):2209-2216; doi:10.1093/carcin/bgl077
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Polymorphisms and haplotypes of the NBS1 gene are associated with risk of sporadic breast cancer in non-Hispanic white women
55 years
Department of Epidemiology, The University of Texas M. D. Anderson Cancer Center Houston, TX 77030, USA
1 Department of Gastrointestinal Medical Oncology, The University of Texas M. D. Anderson Cancer Center Houston, TX 77030, USA
2 Department of Clinical Cancer Prevention, The University of Texas M. D. Anderson Cancer Center Houston, TX 77030, USA
3 Department of Radiation Oncology, The University of Texas M. D. Anderson Cancer Center Houston, TX 77030, USA
4 Department of Pathology, The University of Texas M. D. Anderson Cancer Center Houston, TX 77030, USA
5 Department of Surgical Oncology, The University of Texas M. D. Anderson Cancer Center Houston, TX 77030, USA
*To whom requests for reprints should be addressed at: Department of Epidemiology, Unit 1365, The University of Texas M. D. Anderson Cancer Center, 1515 Holcombe Boulevard, Houston, TX 77030, USA. Tel: +1 713 792 3020; Fax: +1 713 563 0999; Email: lwang{at}mdanderson.org
| Abstract |
|---|
|
|
|---|
DNA double-strand breaks (DSBs) may lead to genomic instability and cancer if unrepaired. Nijmegen breakage syndrome 1 (NBS1) protein is one of the key proteins that participates in recognition and repair of DSBs in humans. We hypothesized that polymorphisms of NBS1 are associated with breast cancer risk. We selected three NBS1 haplotype-tagging polymorphisms (i.e. 924T>C, 8360G>C and 30537G>C) to represent all common (
5%) haplotypes reported in the National Institute of Environmental Health Sciences database and to reconstruct haplotypes. In a hospital-based casecontrol study of 421 non-Hispanic white patients with sporadic breast cancer (
55 years) and 423 cancer-free controls who were frequency-matched with the cases by age (±5 years and
55), we tested our hypothesis and found that compared with 924TT homozygotes the variant homozygote 924CC carriers had a 4.55-fold increased risk of breast cancer [95% confidence interval (CI) = 1.5113.7] and that compared with the 8360GG genotype the variant genotypes were also associated with a significantly increased risk [adjusted odds ratio (OR) = 1.33, 95% CI = 1.001.78 for 8360CG; adjusted OR = 1.83, 95% CI = 1.142.94 for 8360CC]. However, these effects were not observed for the 30537G>C polymorphism. Furthermore, the derived haplotypes were associated with risk in a doseresponse manner as the number of variant (risk) alleles (i.e. 8360C, 924C or 30537C) increased (adjusted OR = 1.07, 95% CI = 0.781.46 for 12 variant alleles; adjusted OR = 2.47, 95% CI = 1.484.14 for 36 variant alleles; Ptrend = 0.006). These findings suggest that NBS1 polymorphisms and haplotypes may contribute to the etiology of sporadic breast cancer in young non-Hispanic white women. Large studies are warranted to confirm these findings.
Abbreviations: BPDE, benzo[a]pyrene diol epoxide; CI, confidence interval; DSB, double-strand break; htSNP, haplotype-tagging SNP; LD, linkage disequilibrium; NBS1, Nijmegen breakage syndrome 1; OR, odds ratio; PCR, polymerase chain reaction; SNP, single nucleotide polymorphism
| Introduction |
|---|
|
|
|---|
Breast cancer is the leading incident cancer among women in the United States, and the incidence rate has increased since the 1980s, with an estimated 212 920 new cases and 40 970 deaths in 2006 (1). Epidemiological studies have established many etiologic factors, including some environmental factors, such as ionizing radiation, high-fat diets, alcohol consumption and use of hormone replacement or oral contraceptives (2). These environmental factors can potentially cause DNA damage and thus lead to a higher risk of breast cancer. However, only a small fraction of exposed individuals develops breast cancer, suggesting that genetic susceptibility plays a role in an individual's risk of developing breast cancer. Therefore, it is important to identify the at-risk populations so that they may be targeted for prevention and early detection.
One of the genetic factors for breast cancer development is a family history of the disease. Germline mutations in high-penetrance breast cancer susceptibility genes, such as BRCA1, BRCA2, p53, ATM and PTEN, have been found to account for the majority of familial hereditary breast cancers (2,3). These genes are all more or less involved in DNA repair; however, because of their low allele frequencies, they only account for 510% of all breast cancer cases (2,3). For most sporadic breast cancer cases, genetic susceptibility is still poorly understood. Low-penetrance but highly prevalent susceptibility genes, acting together with exogenous and endogenous risk factors, are likely to account for many sporadic breast cancers (4).
It has been suggested that DNA repair capacity phenotype may be markers for breast cancer susceptibility (5,6). We have published two pilot studies, in which we found that increased breast cancer risks in young women (
55 years) were significantly associated with a reduced capacity for repairing DNA damage induced by the tobacco carcinogen benzo[a]pyrene diol epoxide (BPDE) (7) and a higher level of BPDE-induced chromatid breaks in peripheral blood lymphocytes (8). Ionizing radiation, one of the established etiologic contributors to breast cancer, can cause DNA damage that is characteristic of DNA double-strand breaks (DSBs). If not repaired, DSBs may cause chromosomal breaks and genomic instability, thus increasing one's probability of developing cancer.
In humans, there are two pathways that repair DSBs, the homologous recombination (HR) and non-homologous end-joining (NHEJ), in which the initial step is to recognize DSBs by a protein complex (MRE11/RAD50/NBS1 complex, MRN). Nijmegen breakage syndrome 1 (NBS1) is the key regulator of this protein complex (911). If DSBs persist, NBS1 promotes the enzymatic activity of RAD50 and MRE11 and directs the MRN complexes to the site of DNA damage. NBS1 also interacts with
H2ax to form nuclear foci and activates the ATM kinase that may activate cell-cycle checkpoint response to initiate the repair of DNA damage and maintain genomic stability (1215). A marked impaired DSB repair has been observed in cells from NBS patients (16). An NBS1 null mutant murine model showed that the mutants had increased chromosome damage, radiomimetic sensitivity and premature cell death (17). Human cell lines with NBS1 knockdown by small interfering RNA showed increased ionizing radiation-induced mutagenesis and telomere instability (18).
The NBS1 gene has been mapped to chromosome 8q21 (1921), and it spans over 50 kb, contains 16 exons and encodes the 754-amino acid protein (Nibrin) (GenBank accession no. AF069291
[GenBank]
). Mutations in NBS1 cause NBS, a chromosome instability syndrome characterized by microcephaly, growth retardation, immunodeficiency, hypersensitivity to ionizing radiation and a predisposition to malignancy (22,23). The primary cause of deaths among NBS patients is cancer; 40% of NBS patients develop cancer before the age of 21 (23). A homozygous mutation, a 5 bp deletion in exon 6 (657del5), has been reported to be associated with an elevated risk of breast cancer in Polish and Russian populations (24,25). However, the 657del5 mutation appears to be of Slavic origin with a low frequency of
0.5% in Eastern European populations and even rarer in other ethnic groups (26). Furthermore, several mutations and polymorphisms in NBS1 were identified in breast cancer and ovarian cancer cell lines (27). On the basis of these studies, it is likely that common polymorphisms of NBS1 may contribute to the risk of sporadic breast cancer.
Several studies have investigated the Glu185Gln polymorphism in exon 5 of the NBS1 gene (2830). One study investigated four single nucleotide polymorphisms (SNPs) (i.e. Leu34Leu in exon 2, Asp399Asp in exon 10, Pro672Pro in exon 13 and Glu185Gln) simultaneously (31) in breast cancer, and the results suggested that these polymorphisms were not associated with risk of breast cancer. However, numerous polymorphic loci that exist in the exons of NBS1 are in complete linkage disequilibrium (LD), although not necessarily functional, and a few SNPs may not be representative of genetic variation in the gene. A more holistic approach using multiple, unlinked genetic variants is needed for such association studies. The use of haplotype-tagging SNPs (htSNPs) in association studies for identification of common disease-associated genes is not only hypothesis-free in selecting the SNPs and greatly reduces genotyping efforts but also appears to have a greater power over single-allele studies (3234). Therefore, we hypothesized that polymorphisms and haplotypes of the NBS1 gene are associated with risk of sporadic breast cancer, particularly in young women, and tested this hypothesis in a casecontrol study that was built on our previous studies of young women (
55 years old) (7,8) with three representative htSNPs based on prior pair-wise LD analysis.
| Materials and methods |
|---|
|
|
|---|
Study population
From August 1998 to August 2005, histopathologically confirmed, untreated primary female breast cancer patients who registered at the Nellie B. Connally Breast Center at The University of Texas M. D. Anderson Cancer Center were recruited. To facilitate the study of genetic susceptibility in a homogeneous population, the case subjects were limited to young women
55 years old and non-Hispanic whites only. We recruited 421 breast cancer cases, consisting of 75 cases of ductal or lobular carcinomas in situ and 346 cases of invasive carcinomas. According to the American Joint Committee on Cancer (AJCC) staging classifications (35), there were 75 cases of stage 0, 159 cases of stage I, 98 cases of stage IIA, 43 cases of stage IIB, 22 cases of stage IIIA, 8 cases of stage IIIB, 12 cases of stage IIIC and 4 cases of stage IV. During the same time period, we randomly recruited 423 cancer-free female controls that were frequency-matched to the cases by age (±5 years) and smoking status (never versus ever) among hospital visitors who were either accompanying the patients but genetically unrelated with the cases included in this study or came to the hospital for our free skin cancer prevention screening and were found cancer-free. A simple questionnaire was used to collect information on demographic data, smoking and alcohol use for all cases and controls. Participants who had smoked at least 100 cigarettes in their lifetime were categorized as ever smokers and the rest were categorized as never smokers. Participants who had drunk alcoholic beverages at least once a week for 1 year or more in their lifetime were categorized as ever drinkers, and the rest were categorized as never drinkers. Additional information on the menopausal status, age at first full-term pregnancy, number of live births, use of hormone replacement or oral contraceptives, weight, height and education level were only available for the case subjects but not for the controls. Each subject was asked to donate 30 ml of blood after an informed consent was obtained. The study protocol was approved by the M. D. Anderson Cancer Center Institutional Review Board.
Selection of htSNPs
As shown in the GenBank SNP database (http://www.ncbi.nlm.nih.gov), there are >317 SNPs that have been identified in the NBS1 gene, of which 249 SNPs are also reported in the Environmental Genome Project (EGP) SNP database of the National Institute of Environmental Health Sciences (NIEHS) (http://egp.gs.washington.edu). On the basis of the available resequencing data, only 84 SNPs have a minor allele frequency
5% in 90 tested individuals, for which 54 sets of haplotypes were reconstructed with the PHASE program (version 2). The structure of LD block construction is shown in Figure 1 based on the NIEHS EGP-SNP database. A large block of high LD with limited haplotype diversity was marked with a triangle diagram. We used the BEST 1.0 program (36) and Stram's tagSNPs program (37,38) with the squared correlation Rh2 > 0.80 to select a minimal set of representative tagging SNPs (htSNPs) from the 84 SNPs included in the common [i.e. the frequency
5% (39)] haplotypes that account for 44.4% of the 54 haplotypes. Three htSNPs were selected and they are located in the high LD block (Figure 1) that has Rh2 = 1.00 with all other common haplotypes. Two of these three htSNPs are located in the 5' untranslated region of NBS1, at positions 626 (626G>A, 1418 nt to initiation transcription code ATG) and 924 (924T>C, 1120 nt), and one is located in intron 9, at position 30537 (30537G>C, rs1805787). Because the variants in the coding regions of the NBS1 gene may play roles in modifying the functions, another four SNPs with a frequency >5% were also selected as candidates from the Cancer Genome Anatomy Project SNP500Cancer Database, which are at the following positions: 3816 in exon 2 (3816G>A, Leu34Leu, rs1063045), 8360 in exon 5 (8360G>C, Glu185Gln, rs1805794), 31129 in exon 10 (31129C>T, Asp399Asp, rs709816) and 40419 in exon 13 (40419A>G, Pro672Pro, rs1061302).
|
On the basis of the genotype data of 90 subjects in the NIEHS SNP database, we calculated the D' values and correlation coefficient r2 with the 2LD program (40) and Statistical Analysis System (SAS)/Genetics software (Version 9.1, SAS Institute, Cary, NC) between any pair of these candidate SNPs and found that locus 626G>A was in complete LD with 8360G>C and 40419A>G (D' = 1.00, P = 0.00 and r2 = 1.00 for both) and locus 8360G>C was in complete LD with 40419A>G. As a result, genotyping for locus 8360G>C only would be sufficient for these three variants. We then selected five htSNPs [924T>C, 8360G>C, 30537G>C, 3816G>A and 31129C>T] and performed a pilot study to genotype 118 non-Hispanic white cancer-free controls to further confirm the LD of these five loci. We found that the SNP 3816G>A was in complete LD with 8360G>C (D' = 1.00, P = 0.00 and r2 = 1.00), and the SNP 31129C>T was in incomplete LD with 8360G>C (D' = 0.91, P = 0.00 and r2 = 0.67). As a result, the SNPs 3816G>A and 31129C>T were eliminated from further genotyping. Because 8360G>C is the mostly studied nsSNP, with the FORCE IN command of Stram's tagSNPs program we replaced 626G>A for 8360G>C, assuring the Rh2 = 1.00.
As shown in Figure 1, the three selected tagging SNPs (i.e. 924T>C, 8360G>C and 30537G>C) are located in the block with high LD and were further used to construct haplotypes for the analysis of association with risk of breast cancer.
Genotyping analysis
From each blood sample, the leukocyte cell pellet obtained from the buffy coat by centrifugation of 1 ml of whole blood was used for genomic DNA extraction with the QIAGEN DNA Blood Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer's instructions. DNA purity and concentrations were determined by spectrophotometric measurement of absorbance at 260 and 280 nm.
The NBS1 924T>C polymorphism was genotyped with the polymerase chain reactionrestriction fragment length polymorphism (PCRRFLP) method by using the following self-designed primers to amplify the SNP-containing fragment: 5'-CCAATGTTAACCCACCCTTTT-3' (forward) and 5'-AGGCTGGTTTACCAATCTCTG-3' (reverse). PCR was performed in 10 µl reaction systems containing 2 mM MgCl2, 0.04 mM deoxynucleotide triphosphates, 1.0 units of Taq polymerase and the manufacturer's buffer [20 mM TrisHCI (pH 8.4) and 50 mM KCl]. The PCR profile consisted of an initial melting step at 95°C for 5 min, 35 cycles of 95°C for 30 s, 61°C for 45 s and 72°C for 50 s and a final extension step of 72°C for 10 min. The amplified fragments were digested with BclI (New England BioLabs, Beverly, MA, USA) overnight, and the products were separated in 3% MetaPhor agarose gel. The 924TT genotype produced two bands (121 and 20 bp), whereas the 924CC genotype produced a single band (141 bp), and the heterozygotes displayed all three bands (141, 121 and 20 bp).
The NBS1 8360G>C polymorphism was also analyzed by the PCRRFLP method. We amplified the SNP-containing fragment by using the following self-designed primers: 5'-CGTCCAATTGTAAAGCCAGAA-3' (forward) and 5'-TCCTGAAACAAGCATTAAAGAGG-3' (reverse). PCR was performed in the 10 µl reaction systems. After an initial denaturation at 94°C for 5 min, there were 35 cycles at 95°C for 45 s, 54°C for 45 s and 72°C for 60 s, and then a final extension step of 72°C for 10 min. The amplified fragments were digested with HinfI (New England BioLabs) overnight, and the products were separated in 2% agarose gel. The 8360GG genotype produced two bands (125 and 49 bp), whereas the 8360CC genotype produced a single band (174 bp), and the heterozygotes displayed all three bands (174, 125 and 49 bp).
The NBS1 30537G>C polymorphism was also genotyped with the PCRRFLP method by using the following self-designed primers: 5'-TTGCTTGAATGACAGCCTGA-3' (forward) and 5'-TCAGCCAGATGGCAGACTC-3' (reverse). The PCR conditions were the same as that for the NBS1 924T>C polymorphism. The amplified fragments were digested with EarI (New England BioLabs) overnight, and the products were separated in 3% MetaPhor agarose gel. The homozygous 30537GG genotype produced a single DNA band (197 bp), whereas the 30537GC genotype produced two bands (172 and 25 bp), and the heterozygotes displayed all three bands (197, 172 and 25 bp).
The PCRs were conducted, and the results were evaluated without the knowledge of the subjects' case or control status. We randomly selected 96 samples (11.37% of 844 samples) for each of three SNPs for repeated assays, and the results were 100% concordant.
Statistical analysis
2-tests were used to compare the distribution of demographic variables and selected risk factors between cases and controls. The HardyWeinberg equilibrium (p2 + 2pq + q2 = 1, where p is the frequency of the variant allele and q =1 p) was tested by a goodness-of-fit
2-test to compare the observed genotype frequencies with the expected genotype frequencies in cancer-free controls. Unconditional logistical regression analysis was used to calculate odds ratios (ORs) and their 95% confidence intervals (CIs) with and without adjustment for age, ethnicity, smoking status and alcohol use. The logistic regression model was also used for the trend test. The 2LD program and Proc ALLELE in SAS/Genetics software were used to detect the LD of any pair of SNPs. Proc HAPLOTYPE in SAS/Genetics software using the expectation-maximization (EM) algorithm was applied to generate maximum likelihood estimates of haplotype frequencies (41). It was assumed that the HardyWeinberg equilibrium is applicable to the constructed haplotypes based on observed genotypes of NBS1. The PHASE program using Bayesian methods with Markov chain Monte Carlo algorithm was also used to construct haplotypes (42,43), and the results from PHASE and SAS were compared. Haplotype and diplotypes (haplotype pairs) were combined on the basis of haplotype information. The associations between the diplotypes and breast cancer risk were evaluated with logistic regression models. The joint effects of the NBS1 polymorphisms, smoking and drinking on breast cancer risk were calculated. All statistical tests were two-sided, a P-value of 0.05 being considered significant, by using SAS software (version 8.2; SAS Institute, Cary, NC).
| Results |
|---|
|
|
|---|
Characteristics of the study population
This analysis included 421 breast cancer cases and 423 cancer-free controls. Because we used frequency-matching, there were no significant differences in the distributions of age, smoking status and alcohol usage between the cases and controls (P = 0.545, 0.397 and 0.241, respectively). However, case subjects were slightly older (mean age ± standard deviation = 45.4 ± 6.6 years, range = 2255 years) than control subjects (mean age ± standard deviation = 44.9 ± 8.2 years, range = 2155 years), but the difference was not statistically significant (Student's t-test: P = 0.329). Additionally, multivariate logistic regression models were used to control these variables for any residual effect of confounding on the main effect of the selected NBS1 polymorphisms.
NBS1 genotype distributions and risk of breast cancer
The genotype and allele distributions of the NBS1 924T>C, 8360G>C and 30537G>C polymorphisms among the cases and controls are summarized in Table I. The observed genotype frequencies of these three SNPs were all in agreement with the HardyWeinberg equilibrium in the control subjects (P = 0.515, 0.603 and 0.212, respectively). The LD analyses in controls showed that the 924T>C locus and 8360G>C locus were relatively weak (D' = 0.166, P = 0.0024 and r2 = 0.0079), the 8360 C locus was not in LD with the 30537G>C locus (D' = 0.127, P = 0.319 and r2 = 0.0012) and the 924T>C locus was not in LD with the 30537G>C locus either (D' = 0.007, P = 0.977 and r2 = 0.0000), further confirming that these are independent tagging SNPs.
|
For the NBS1 924T>C polymorphism, compared with the 924TT genotype, only the homozygous 924CC carriers, but not the heterozygous 924TC carriers, had 4.55-fold elevated risk of breast cancer (95% CI = 1.5113.7). The NBS1 8360G>C genotype frequencies were 40.6% (GG), 47.0% (GC), and 12.4% (CC) in the cases and 48.9% (GG), 42.8% (CG) and 8.3% (CC) in the controls, and the difference was statistically significant (P = 0.023). The variant C allele frequency in the cases (35.9%) was also statistically higher than that of controls (29.7%) (P = 0.007). Logistic regression analysis showed that compared with the 8360GG genotype, the heterozygote 8360GC carriers had 1.33-fold elevated risk of breast cancer (95% CI = 1.001.78), and the homozygote 8360CC carriers had 1.83-fold elevated risk of breast cancer (95% CI = 1.142.94), and this trend was statistically significant (Ptrend = 0.005). For the NBS1 30537G>C polymorphism, the 30537CC variant homozygous carriers had a non-significantly increased risk (1.90-fold, 95% CI = 0.695.22), and the heterozygous 30537GC carriers had no risk for breast cancer at all, compared with the 30537GG genotype carriers (Table I). These data suggest that the 8360C polymorphism is a risk allele, since it causes amino acid change, whereas 924C and 30537C polymorphisms may be risk alleles or in linkage with other unknown risk alleles.
Distribution of NBS1 haplotypes/diplotypes and risk of breast cancer
To further analyze the combined effect of these three tagging SNPs, we generated the haplotypes on the basis of the observed genotypes. As shown in Table II, the PHASE and SAS Proc HAPLOTYPE programs both showed eight possible NBS1 haplotypes from the observed genotypes. The five common (
5%) haplotype alleles (TGG, TCG, TGC, CCG and TCC) represented 94.6% of the chromosomes in the controls. The frequency distribution of haplotype alleles between the cases and controls was not significantly different (P = 0.189). However, when we used the most common haplotype allele TGG as the reference, the CCG haplotype allele was associated with a 74% increased risk of breast cancer (OR = 1.74; 95% CI = 1.172.59). Because the frequencies of all minor variant alleles in the three selected htSNPs were associated with increased risk of breast cancer (Table I), and the frequencies of some haplotypes were <5%, we recombined all possible haplotypes into three groups on the basis of the number of minor variants within the haplotype alleles, assuming that the effects of the minor variant alleles would follow an additive model. Given that an individual's genotype of the gene is a mix of all possible SNPs, the analysis of diplotypes would provide the estimate for an actual effect of such a genotype on the risk of breast cancer. On the basis of the number of variants within the haplotype alleles, we recombined the haplotypes into three groups. The difference in frequency distributions of these haplotype groups among cases and controls was statistically significant (P = 0.012). Compared with the non-variant haplotype group Hap0, the risk for breast cancer was increased as the number of variant alleles increased in a doseresponse manner [OR = 1.11, 95% CI = 0.891.39 for Hap1 (i.e. having 1 variant allele) and OR = 1.51; 95% CI = 1.14-2.00 for Hap2 (i.e. having
2 variant alleles); Ptrend = 0.004] (Table II).
|
We further examined the distribution of the diplotypes based on the combination of the eight possible haplotypes. As shown in Table III, 25 possible diplotypes were generated. Using the most common TGG/TGG diplotype as the reference, only one diplotype CCG/CGG was associated with a 9.21-fold elevated risk of breast cancer (95% CI = 1.1077.3). When we divided these 25 diplotypes into three groups (i.e. no variant allele, 12 variant alleles and 36 variant alleles), the distributions between cases and controls were significantly different (P = 0.002). When the genotypes without any variant allele (i.e. any of 924C, 8360C and 30537C) were used as the reference, the cancer risk increased as the number of variant alleles increased in a doseresponse manner (adjusted OR = 1.07, 95% CI = 0.781.46 for 12 variant alleles; adjusted OR = 2.47; 95% CI = 1.484.14 for 36 variant alleles; Ptrend = 0.006) (Table III).
|
Stratification analysis of NBS1 diplotypes and risk of breast cancer
To take into account all three variants, we used the combined NBS1 diplotypes for stratification analysis by age, smoking status, alcohol use and tumor histological type. Compared with a zero variant allele, the trend of significantly increased cancer risk with an increased variant number of alleles was also obvious in subgroups of cases who were <45 years (adjusted OR = 1.30, 95% CI = 0.782.18 for 12 variant alleles; adjusted OR = 2.09, 95% CI = 1.014.35 for 36 variant alleles; Ptrend = 0.048), never smokers (adjusted OR = 1.09, 95% CI = 0.731.65 for 12 variant alleles; adjusted OR = 2.70, 95% CI = 1.375.30 for 36 variant alleles; Ptrend = 0.021) and never drinkers (adjusted OR = 1.32, 95% CI = 0.852.07 for 12 variant alleles; adjusted OR = 3.33, 95% CI = 1.587.04 for 36 variant alleles; Ptrend = 0.005) (Table IV).
|
However, there was no any difference between tumor histological types (i.e. ductal or lobular carcinomas in situ and invasive carcinomas) or among the AJCC stages in risk associated with the NBS1 polymorphism of 924T>C, 8360G>C, 30537G>C and its diplotypes (data not shown). Furthermore, there were no significant differences in the distribution of subjects' characteristics [i.e. menopausal status, age at first full-term pregnancy, number of live births, use of hormone replacement or oral contraceptives, education level and body mass index (BMI, kg/m2) by the NBS1 genotypes] (data not shown). Taken together, these data suggest that the NBS1 variant genotypes and diplotypes may be biomarkers for susceptibility rather than a disease biomarker. Finally, we analyzed the joint effects of the NBS1 polymorphism of 924T>C, 8360G>C, 30537G>C and its diplotypes with smoking and drinking exposure on breast cancer risk. We did not find any evidence of modification or interaction (data not shown).
| Discussion |
|---|
|
|
|---|
We have previously reported that young breast cancer patients tended to have suboptimal DNA repair capacity for removing BPDE-induced DNA adducts (7) or have higher levels of BPDE-induced chromosomal aberrations (8). In the present study, we were interested in the association between htSNPs of the NBS1 gene and risk for sporadic breast cancer. When we evaluated each of the selected htSNPs separately, the NBS1 8360C (185Gln) variant genotypes were found to be associated with increased risk in an allele-dose response manner, whereas only the homozygous 924CC genotype was associated with an increased risk. The homozygous 30537CC was associated with non-significantly increased risk. However, when the risk alleles (i.e. 924C, 8360C and 30537C) were combined either in the form of haplotypes or diplotypes, the risk for breast cancer appeared to increase as the number of risk alleles increased. Although the functional relevance of these selected htSNPs remains unknown, our findings are biologically plausible, because the htSNPs are either non-synonymous (Glu185Gln) or located in the promoter region.
The NBS1 8360G>C polymorphism is located in the area coding for breast cancer C-terminal (BRCT) domain (108196 amino acids) (9), which facilitates the interaction between NBS1 and BRCA1 (breast cancer 1), forming a BRCA1-associated genome surveillance complex (BASC), which is responsible for recognition and repair of aberrant DNA (4447). This might partially explain the possible effect of the 8360G>C polymorphism on breast cancer risk, that is, the amino acid change caused by the variant may interfere with such a proteinprotein interaction. However, this hypothesis needs to be further tested.
Several studies have investigated the association between polymorphisms of NBS1 8360G>C (Glu185Gln) and cancer risk, but the results are mixed. A recent large U.S. study of 894 AfricanAmerican breast cancer cases and 788 controls and 1417 whites breast cancer cases and 1234 controls, between 21 and 74 years, reported that Glu185Gln variant genotypes were not associated with risk of breast cancer (30). A similar result was reported in a large casecontrol study of 2205 breast cancer cases and 1826 controls, between 45 and 74 years, in Europe (31). We do not know what percentage of the subjects in these studies were 55 years or younger; it is possible that the analysis of mixed age groups may mask the moderate effect of Glu185Gln in young women. A combined study of 221 Finnish sporadic breast cancer cases and Polish familial cases and 319 controls also reported no association between the Glu185Gln polymorphism and breast cancer (28), as did a study of a Chinese population of 220 cases and 310 controls (29), but these studies were relatively small.
For other cancers, however, Medina et al. (48) report that the NBS1 185Gln variant genotypes (8360GC/CC) were associated with a p53 mutation in lung cancer and concluded that the genotypes may play a role in lung carcinogenesis. However, a study in a Chinese population reported that 8360CC (42.3% in controls) was the common genotype and that the 8360GG genotype was associated with increased risk of lung cancer (49). Other studies, including a casecontrol study of non-small cell lung cancer in Norway (50), a multi-center study of ovarian cancer from the United Kingdom, Denmark and the United States (51) and a study of non-Hodgkin's lymphoma, have not found any associations of NBS1 185Gln and cancer risk (52).
Compared with other studies, our present study is rather unique. We included relatively young (
55 years) non-Hispanic white women who were frequency-matched with controls on smoking status because genetic susceptibility is often characterized by an early age of onset and smoking is a potential risk factor for breast cancer at a young age. Therefore, a positive finding between the association of the 8360G>C polymorphism and breast cancer risk may be because of the inclusion of an at-risk population and frequency-matching, which may have enhanced the statistical power of our study. Furthermore, the findings from the present study are consistent with the findings from our two previously published pilot studies in which we found that increased risk of breast cancer in young women was significantly associated with a reduced capacity for repairing DNA damage and with chromosomal aberrations induced by a tobacco carcinogen (7,8).
However, it is likely that our results could be biased owing to possible selection biases in hospital-based casecontrol studies. Because the genotype frequencies of the NBS1 8360G>C polymorphism estimated from the hospital-based control subjects in our study were very close to the resequencing data in the Cancer Genome Anatomy Project SNP500 Cancer Database for Caucasians (http://snp500cancer.nci.nih.gov), the potential selection bias in genotype distribution is unlikely to be substantial, if any. Nevertheless, larger studies with young women are warranted to further test whether the NBS1 Glu185Gln variant is associated with breast cancer risk.
In this study, we also found that variant homozygote 924CC carriers had a 4.55-fold elevated risk of breast cancer compared with the homozygote 924TT carriers, whereas the 30537CC polymorphism was associated with a non-significantly increased risk of breast cancer compared with the homozygote 30537GG carriers. The significance of these findings, however, is limited because of the small number of cases in the subgroups. To the best of our knowledge, no reported study has demonstrated an association between NBS1 924T>C or 30537G>C polymorphisms and cancer risk. The NBS1 924T>C polymorphism locates in the 5' untranslated region (1120 nt) of the NBS1 gene, which has been shown by bioinformatics analysis that it is the transcription factor GATA-1 binding site. The activation domains of GATA-1 are capable of activating transcription in mammalian cells through GATA motifs (53). When this locus changes to 924C, it might cause a loss of the GATA-1 binding protein cis-element and thus may have an adverse effect on NBS1 gene expression.
The analyses of haplotypes and diplotypes with the combination of the three NBS1 variants showed that an elevated risk of breast cancer was associated with haplotypes or diplotypes in an allele-dose response manner. This suggests that NBS1 924T>C or 30537G>C polymorphisms also have a potential role in the etiology of breast cancer. Our analysis on the haplotype and diplotype provided additional information on these SNPs as markers of genetic susceptibility and a more efficient method for assessing the genetic susceptibility of a candidate gene than any of the polymorphisms. Further stratification analysis showed that the increased risk associated with variant diplotypes was more pronounced in the subgroup of never smokers and never drinkers, consistent with the concept that the susceptible individuals may have been exposed to a low level of carcinogens for a shorter period of time. However, these results may also be biased by the relatively small number of subjects in the subgroups. Further studies with large sample sizes are needed to verify these findings.
Our study has inherited shortcomings of a hospital-based casecontrol study. First of all, because our cases and controls were recruited from a hospital, they may not be representative of the general population and could have introduced selection bias. Also, because we restricted the subjects to non-Hispanic whites, it is uncertain whether there may be population stratification and whether these results are generalizable to other populations. A recent article reported the comparison of haplotype information between databases of HapMap and NIEHS; the authors indicated that the NIEHS provided more information with resequenced SNPs than that of HapMap, but it lacks the ethnic information of each individual (54). Furthermore, the tagging SNPs in our testing set of 118 non-Hispanic whites were also showing differences in LD compared with those of HapMap. Therefore, potential population stratification may have contributed to the differences observed in the different data sets. Finally, limited information on reproductive factors, particularly hormonal use, in this study population did not allow us to explore the interactions between genes and the environment.
In conclusion, in this hospital-based casecontrol study of sporadic breast cancer, we found a significant association between NBS1 polymorphisms and breast cancer risk in non-Hispanic whites. Furthermore, the haplotypes of NBS1 generated from the observed genotypes were associated with significantly increased cancer risk in an allele/dose response manner, especially in non-smokers and non-drinkers. This supports the notion that variant genotypes of the NBS1 gene may contribute to the risk of sporadic breast cancer, particularly in young women. However, because we used less stringent r2 values (i.e. 0.5 instead of 0.8) to define complete LD between SNPs in the SNP selection (55), we could have suffered some loss of our statistical power in our risk estimates. Furthermore, because the study design was retrospective, the observed associations may have been biased or simply due to chance. Therefore, larger, prospective studies are needed to verify these findings.
| Acknowledgments |
|---|
This study was supported in part by National Institutes of Health grants CA 108364 (L.E.W.), ES11740 (Q.W.) and CA 16672 (M. D. Anderson Cancer Center). We thank Margaret Lung, Maria Torres, and Leanel Fairly for their assistance in recruiting the subjects; Qiuling Shi for statistical support; Zhibin Hu, Luo Wang, Zhensheng Liu, Zhaozheng Guo, Yawei Qiao, Jianzhong He, Kejing Xu and Yinyan Li for their laboratory assistance; Joanne M. Sider for manuscript preparation; and Martha Morrison for scientific editing.
Conflict of Interest Statement: None declared.
| References |
|---|
|
|
|---|
- Jemal A., Siegel R., Ward E., Murray T., Xu J., Smigal C., Thun M.J. (2006) Cancer statistics, 2006. CA Cancer J. Clin. 56:106130.
[Abstract/Free Full Text] - Singletary S.E. (2003) Rating the risk factors for breast cancer. Ann. Surg. 4:474482.
- Thompson D. and Easton D. (2004) The genetic epidemiology of breast cancer genes. J. Mammary Gland Biol. Neoplasia 9:221236.[CrossRef][ISI][Medline]
- Dumitrescu R.G. and Cotarla I. (2005) Understanding breast cancer riskwhere do we stand in 2005? J. Cell Mol. Med. 9:208221.[ISI][Medline]
- Hu J.J., Smith T.R., Miller M.S., Lohman K., Case L.D. (2002) Genetic regulation of ionizing radiation sensitivity and breast cancer risk. Environ. Mol. Mutagen 39:208215.[CrossRef][ISI][Medline]
- Kennedy D.O., Agrawal M., Shen J., Terry M.B., Zhang F.F., Senie R.T., Motykiewicz G., Santella R.M. (2005) DNA repair capacity of lymphoblastoid cell lines from sisters discordant for breast cancer. J. Natl Cancer Inst. 97:127132.
[Abstract/Free Full Text] - Shi Q., Wang L.E., Bondy M.L., Brewster A., Singletary S.E., Wei Q. (2004) Reduced DNA repair of benzo[a]pyrene diol epoxide-induced adducts and common XPD polymorphisms in breast cancer patients. Carcinogenesis 9:16951700.
- Xiong P., Bondy M.L., Li D., Shen H., Wang L.E., Singletary S.E., Spitz M.R., Wei Q. (2001) Sensitivity to benzo(a)pyrene diol epoxide associated with risk of breast cancer in young women and modulation by glutathione S-transferase polymorphisms: a casecontrol analysis. Cancer Res. 61:84658469.
[Abstract/Free Full Text] - Kobayashi J., Antoccia A., Tauchi H., Matsuura S., Komatsu K. (2004) NBS1 and its functional role in the DNA damage response. DNA Repair 3:855861.[Medline]
- Matsuura S., Kobayashi J., Tauchi H., Komatsu K. (2004) Nijmegen breakage syndrome and DNA double strand break repair by NBS1 complex. Adv. Biophys. 38:6580.[CrossRef][ISI][Medline]
- Stephen P.J. (2002) Sensing and repairing DNA double-strand breaks. Carcinogenesis 23:687696.
[Abstract/Free Full Text] - Paull T.T. and Lee J.H. (2005) The Mre11/Rad50/NBS1 complex and its role as a DNA double-strand break sensor for ATM. Cell Cycle 4:737740.[ISI][Medline]
- Iijima K., Komatsu K., Matsuura S., Tauchi H. (2004) The Nijmegen breakage syndrome gene and its role in genome stability. Chromosoma 113:5361.[ISI][Medline]
- Stavridi E.S. and Halazonetis T.D. (2004) Nbs1 moving up in the world. Nat. Cell Biol. 7:648650.
- Stiff T., Reis C., Alderton G.K., Woodbine L., O'Driscoll M., Jeggo P.A. (2005) NBS1 is required for ATR-dependent phosphorylation events. EMBO J. 24:199208.[CrossRef][ISI][Medline]
- Howlett N.G., Scuric Z., D'Andrea A.D., Schiestl R.H. (2006) Impaired DNA double strand break repair in cells from Nijmegen breakage syndrome patients. DNA Repair (Amst) 5:251257.[Medline]
- Demuth I., Frappart P.O., Hildebrand G., et al. (2004) An inducible null mutant murine model of Nijmegen breakage syndrome proves the essential function of NBS1 in chromosomal stability and cell viability. Hum. Mol. Genet. 13:23852397.
[Abstract/Free Full Text] - Zhang Y., Lim C.U., Williams E.S., Zhou J., Zhang Q., Fox M.H., Bailey S.M., Liber H.L. (2005) NBS1 knockdown by small interfering RNA increases ionizing radiation mutagenesis and telomere association in human cells. Cancer Res. 65:55445553.
[Abstract/Free Full Text] - Varon R., Vissinga C., Platzer M., et al. (1998) Nibrin, a novel DNA double-strand break repair protein, is mutated in Nijimegen break syndrome. Cell 93:467476.[CrossRef][ISI][Medline]
- Carney J.P., Maser R.S., Olivares H., Davis E.M., Le Beau.M., Yates J.R. III, Hays L., Morgan W.F., Petrini J.H. (1998) The hMre11/hRad50 protein complex and Nijimegen break syndrome: linkage of double strand break repair to cellular DNA damage response. Cell 93:477486.[CrossRef][ISI][Medline]
- Matsuura S., Tauchi H., Nakamura A., et al. (1998) Positional cloning of the gene for Nijimegen break syndrome. Nat. Genet. 19:179181.[CrossRef][ISI][Medline]
- Weemaes C.M.R., Hustinx T.W.J., Scheres J.M.C., van Munster P.J.J., Bakkeren J.A.J.M., Taalman R.D.F.M. (1981) A new chromosomal instability disorder: the Nijimegen break syndrome. Acta Paediatr. Scand. 70:557564.[ISI][Medline]
- The International Nijimegen Break Syndrome Study Group. (2000) Nijimegen break syndrome. Arch. Dis. Child 82:400406.
[Abstract/Free Full Text] - Steffen J., Varon R., Mosor M., et al. (2004) Increased cancer risk of heterozygotes with NBS1 germline mutation in Poland. Int. J. Cancer 111:6771.[CrossRef][ISI][Medline]
- Buslov K.G., Iyevleva A.G., Chekmariova E.V., et al. (2005) NBS1 657del mutation may contribute to a limited fraction of breast cancer cases in Russia. Int. J. Cancer 114:585589.[CrossRef][ISI][Medline]
- Carlomagno F., Chang-Claude J., Dunning A.M., Ponder B.A. (1999) Determination of the frequency of the common 657Del5 Nijmegen breakage syndrome mutation in the German population: no association with risk of breast cancer. Genes Chromosomes Cancer 25:393395.[CrossRef][ISI][Medline]
- Tessitore A., Biordi L., Flati V., et al. (2003) New mutations and protein variants of NBS1 are identified in cancer cell lines. Genes Chromosomes Cancer 36:198204.[CrossRef][ISI][Medline]
- Forsti A., Angelini S., Festa F., et al. (2004) Single nucleotide polymorphisms in breast cancer. Oncol. Rep. 11:917922.[ISI][Medline]
- Zhang L., Zhang Z., Yan W. (2005) Single nucleotide polymorphisms for DNA repair genes in breast cancer patients. Clin. Chim. Acta 359:150155.[CrossRef][ISI][Medline]
- Millikan R.C., Player J.S., de Cotret A.R., Tse C.K., Keku T. (2005) Polymorphisms in DNA repair genes, medical exposure to ionizing radiation, and breast cancer risk. Cancer Epidemiol. Biomarkers Prev. 14:23262334.
[Abstract/Free Full Text] - Kuschel B., Auranen A., McBride S., et al. (2002) Variants in DNA double-strand break repair genes and breast cancer susceptibility. Human Mol. Genet. 12:13991407.
- Johnson G.C., Esposito L., Barratt B.J., et al. (2001) Haplotype tagging for the identification of common disease gene. Nat. Genet. 29:233237.[CrossRef][ISI][Medline]
- Bonnen P.E., Wang P.J., Kimmel M., Chakraborty R., Nelson D. (2002) Haplotype and linkage disequilibrium architecture for human cancer-associated genes. Genome Res. 12:18461853.
[Abstract/Free Full Text] - Rebbeck T.R., Ambrosone C.B., Bell D.A., Chanock S.J., Hayes R.B., Kadlubar F.F., Thomas D.C. (2004) SNPs, haplotypes, and cancer: applications in molecular epidemiology. Cancer Epidemiol. Biomarkers Prev. 13:681687.
[Free Full Text] - Greene F.L., Page D.L., Fleming I.D., Fritz A., Balch C.M., Haller D.G., Morrow M. (2002) AJCC Cancer Staging Manual. 6th Edition , Springer, NY.
- Sebastiani P., Lazarus R., Weiss S.T., Kunkel L.M., Kohane I.S., Ramoni M.F. (2003) Minimal haplotype tagging. Proc. Natl Acad. Sci. USA 17:99009905.
- Stram D.O., Haiman C.A., Hirschhorn J.N., Altshuler D., Kolonel L.N., Henderson B.E., Pike M.C. (2003) Choosing haplotype-tagging SNPs based on unphased genotype data using a preliminary sample of unrelated subjects with an example from the multiethnic cohort study. Hum. Hered. 55:2736.[CrossRef][ISI][Medline]
- Stram D.O. (2004) Tag SNP selection for association studies. Genet. Epidemiol. 27:365374.[CrossRef][ISI][Medline]
- Haiman C.A., Stram D.O., Pike M.C., Kolonel L.N., Burtt N.P., Altshuler D., Hirschhorn J., Henderson B.E. (2003) A comprehensive haplotype analysis of CYP19 and breast cancer risk: the Multiethnic Cohort. Hum. Mol. Genet. 12:26792692.
[Abstract/Free Full Text] - Zhao J.H. (2004) 2LD, gene counting and Hap: computer programs for linkage disequilibrium analysis. Bioinformatics 20:13251326.
[Abstract/Free Full Text] - Stram D.O., Leigh P.C., Bretsky P., Freedman M., Hirschhorn J.N., Altshuler D., Kolonel L.N., Henderson B.E., Thomas D.C. (2003) Modeling and E-M estimation of haplotype-specific relative risks from genotype data for a casecontrol study of unrelated individuals. Hum. Hered. 55:179190.[CrossRef][ISI][Medline]
- Stephens M., Smith N., Donnelly P. (2001) A new statistical method for haplotype reconstruction from population data. Am. J. Hum. Genet. 68:978989.[CrossRef][ISI][Medline]
- Stephens M. and Donnelly P. (2003) A comparison of Bayesian methods for haplotype reconstruction from population genotype data. Am. J. Hum. Genet. 73:11621169.[CrossRef][ISI][Medline]
- Wang Y., Cortez D., Yazdi P., Neff N., Elledge S.J., Qin J. (2000) BASC, a super complex of BRCA1-associated proteins involved in the recognition and repair of aberrant DNA structures. Genes Dev. 14:927939.
[Abstract/Free Full Text] - Futaki M. and Liu J.M. (2001) Chromosomal breakage syndromes and the BRCA1 genome surveillance complex. Trends Mol. Med. 7:560565.[CrossRef][ISI][Medline]
- Zhao S., Renthal W., Lee E.Y. (2002) Functional analysis of FHA and BRCT domains of NBS1 in chromatin association and DNA damage responses. Nucleic Acids Res. 30:48154822.
[Abstract/Free Full Text] - Kobayashi J. (2004) Molecular mechanism of the recruitment of NBS1/hMRE11/hRAD50 complex to DNA double-strand breaks: NBS1 binds to gamma-H2AX through FHA/BRCT domain. J. Radiat. Res. 45:473478.[Medline]
- Medina P.P., Ahrendt S.A., Pollan M., Fernandez P., Sidransky D., Sanchez-Cespedes M. (2003) Screening of homologous recombination gene polymorphisms in lung cancer patients reveals an association of NBS1-185Gln variant and p53 gene mutations. Cancer Epidemiol. Biomarkers Prev. 12:699704.
[Abstract/Free Full Text] - Lan Q., Shen M., Berndt S.I., et al. (2005) Smoky coal exposure, NBS1 polymorphisms, p53 protein accumulation, and lung cancer in Xuan Wei, China. Lung Cancer 49:317323.[CrossRef][ISI][Medline]
- Zienolddiny S., Campa D., Lind H., Ryberg D., Skaug V., Stangeland L., Phillips D.H., Canzian F., Haugen A. (2006) Polymorphisms of DNA repair genes and risk of non-small cell lung cancer. Carcinogenesis 27:560567.
[Abstract/Free Full Text] - Auranen A., Song H., Waterfall C., Dicioccio R.A., Kuschel B., Kjaer S.K., Hogdall E., Hogdall C., Stratton J., Whittemore A.S., Easton D.F., Ponder B.A., Novik K.L., Dunning A.M., Gayther S., Pharoah P.D. (2005) Polymorphisms in DNA repair genes and epithelial ovarian cancer risk. Int. J. Cancer 117:611618.[CrossRef][ISI][Medline]
- Cerosaletti K.M., Morrison V.A., Sabath D.E., Willerford D.M., Concannon P. (2002) Mutations and molecular variants of the NBS1 gene in non-Hodgkin lymphoma. Gene Chromosomes Cancer 35:282286.
- Merika M. and Orkin S.H. (1993) DNA-binding specificity of GATA family transcription factors. Mol. Cell. Biol. 13:39994010.
[Abstract/Free Full Text] - Taylor J.A., Xu Z.L., Kaplan N.L., Morris R.W. (2006) How well do HapMap haplotypes identify common haplotypes of genes? A comparison with haplotypes of 334 gene resequenced in the Environmental Genome Project. Cancer Epidemiol. Biomarkers Prev. 15:133137.
[Abstract/Free Full Text] - Carlson C.S., Eberle M.A., Rieder M.J., Smith J.D., Kruglyak L., Nickerson D.A. (2003) Additional SNPs and linkage-disequilibrium 95 analyses are necessary for whole-genome association studies in humans. Nat. Genet. 33:518521.[CrossRef][ISI][Medline]
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
J. Huang, M. A. Grotzer, T. Watanabe, E. Hewer, T. Pietsch, S. Rutkowski, and H. Ohgaki Mutations in the Nijmegen Breakage Syndrome Gene in Medulloblastomas Clin. Cancer Res., July 1, 2008; 14(13): 4053 - 4058. [Abstract] [Full Text] [PDF] |
||||
![]() |
H.-M. Hsu, H.-C. Wang, S.-T. Chen, G.-C. Hsu, C.-Y. Shen, and J.-C. Yu Breast Cancer Risk Is Associated with the Genes Encoding the DNA Double-Strand Break Repair Mre11/Rad50/Nbs1 Complex Cancer Epidemiol. Biomarkers Prev., October 1, 2007; 16(10): 2024 - 2032. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


