- Research article
- Open Access
Genetic variation in the NBS1, MRE11, RAD50 and BLM genes and susceptibility to non-Hodgkin lymphoma
BMC Medical Genetics volume 10, Article number: 117 (2009)
Translocations are hallmarks of non-Hodgkin lymphoma (NHL) genomes. Because lymphoid cell development processes require the creation and repair of double stranded breaks, it is not surprising that disruption of this type of DNA repair can cause cancer. The members of the MRE11-RAD50-NBS1 (MRN) complex and BLM have central roles in maintenance of DNA integrity. Severe mutations in any of these genes cause genetic disorders, some of which are characterized by increased risk of lymphoma.
We surveyed the genetic variation in these genes in constitutional DNA of NHL patients by means of gene re-sequencing, then conducted genetic association tests for susceptibility to NHL in a population-based collection of 797 NHL cases and 793 controls.
114 SNPs were discovered in our sequenced samples, 61% of which were novel and not previously reported in dbSNP. Although four variants, two in RAD50 and two in NBS1, showed association results suggestive of an effect on NHL, they were not significant after correction for multiple tests.
These results suggest an influence of RAD50 and NBS1 on susceptibility to diffuse large B-cell lymphoma and marginal zone lymphoma. Larger association and functional studies could confirm such a role.
Non-Hodgkin lymphoma (NHL) is a heterogeneous group of hematological malignancies that in aggregate constitutes the 5th highest cause of cancer mortality in the United States  and Canada . NHL subtypes vary in presentation, survival expectation, morbidity and responses to treatment. Chromosomal translocations are so characteristic of NHL that many genes now known to be important in the development of cancer, such as BCL2 , were originally discovered due to their position at recurrent translocation breakpoints in NHL tumours.
During development and differentiation, the DNA of B- and T-cells is subject to double stranded breaks necessary for the rearrangement of immunoglobulin genes. Genes functioning in double-stranded break repair are involved in successfully controlling and repairing these breaks, thus protecting the genome from molecular events that could lead to cancer. This study examined four genes with key roles in maintaining genome stability: the MRN complex, MRE11, RAD50 and NBS1, and the Bloom syndrome gene (BLM). We have previously shown association with NHL of a genetic variant in H2AX, which encodes a histone involved in signalling the presence of double stranded breaks . The MRN complex forms foci at sites of double stranded breaks induced by ionizing radiation or immunoglobulin rearrangements during B- and T-cell development, sensing DNA damage and initiating DNA repair [5–7].
The chromosome instability syndromes (reviewed in ) form a group of rare autosomal recessive diseases characterized by an increased risk of cancer. This group includes ataxia-telangiectasia (AT, OMIM 208900), Nijmegen breakage syndrome (NBS, OMIM 251260), Bloom syndrome (OMIM 210900) and Fanconi anemia (OMIM 227650). NBS includes an increased risk of lymphoid malignancies , particularly B-cell lymphoma [10, 11]. Some patients with an NBS-like phenotype have mutations in RAD50 . Hypomorphic mutations in MRE11 result in an AT-like disorder (AT-LD). NBS and AT-LD share many features, including immunodeficiency and genome instability caused by failure of timely activation of cell cycle checkpoint pathways [13–16].
Mutations in NBS1 cause aplastic anemia and acute lymphoblastic leukemia [17, 18]. RAD50 variants have also been associated with an increased risk of sporadic , but not necessarily familial breast cancer [19, 20]. MRE11 inactivation has been identified in colorectal cancer cell lines and primary tumours , suggesting that inactivation of the MRN complex could be a frequent event in cancers.
Bloom syndrome is also marked by a predisposition to cancer, particularly lymphoma and leukemia in young patients . Although homozygous loss of Blm in mice leads to embryonic lethality, heterozygotes show increased risk of neoplasia, with augmented T-cell tumourigenesis . This haploinsufficiency is supported by the increased risk of cancer in BLM heterozygotes of Ashkenazi Jewish descent , although there is some controversy regarding this finding . This illustrates BLM's role in response to DNA damage , particularly during DNA replication stress .
While both Nbs1  and Mre11  null mutants are inviable in vertebrates, the hypermorphic Rad50 Smutation causes hematopoietic stem cell failure so that mice that do not die of lymphoma die of bone marrow attrition , highlighting the delicate balance the MRN complex exerts on cell survival. This is illustrated by the dosage sensitivity to this mutation and the bidirectional phenotypic rescue in Rad50 S/S Atm -/- mice , leading the authors to speculate that while mutations that cause gross chromosomal instability would have a wide array of outcomes, less severe mutations would primarily affect tissues developed from a limited number of precursor stem cells. Since the hematopoietic system is such a system, this reinforces the need to look for variants in genes already known to be associated with severe genetic disorders, with the rationale that varying degrees of mutation severity affect the spectrum of possible effects.
To systematically investigate the role of NBS1, MRE11, RAD50 and BLM in susceptibility to NHL, we carried out re-sequencing of these four genes to establish the spectrum of genetic variation in NHL cases, and genotyped 797 NHL cases and 793 controls. Just as total inactivation of a gene and attenuation of its activity lead to different phenotypes in mice, we expected that subtle variation in DNA repair genes could be pertinent to NHL risk in the general population, while complete inactivation of these genes leads to rare and severe syndromes.
The methodology has been described previously [32, 33]. Informed consent was obtained as approved by the joint University of British Columbia/British Columbia Cancer Agency Research Ethics Board. All HIV-negative NHL cases diagnosed in British Columbia from March 2000 to February 2004, residing in the Greater Vancouver Regional District and greater Victoria (Capital Regional District), aged 20 to 79 were invited to participate. Cases were reviewed and coded using the World Health Organization classification by an experienced lymphoma pathologist (RDG). Population controls were identified from the Client Registry of the British Columbia Ministry of Health and were frequency matched to cases by sex, age, and area of residence in a 1:1 ratio. 828 cases and 848 controls completed at least part of a study questionnaire; however, only those subjects with DNA available were used in this study. Table 1 summarizes the characteristics of the 797 cases and 793 controls available for analysis.
DNA extraction and sequencing
Genomic DNA was extracted from whole blood (in 10% of cases from a mouthwash or saliva sample) using the PureGene DNA isolation kit (Gentra Systems) following manufacturer's instructions. DNA was then quantified using PicoGreen (Molecular Probes) in a Victor2 fluorescence plate reader (Perkin-Elmer).
The genomic sequences for all genes were downloaded from the UCSC genome browser . All coding and non-coding exons were sequenced, as well as 1000 base pairs upstream of transcription start. Conserved non-coding sequence regions (CNS regions) were identified using the VISTA genome browser . The six most highly conserved CNS regions with at least 100 base pairs of at least 70% identity with the mouse and rat homolog were also sequenced.
Primers were selected for all amplicons using Primer3 . The -21M13F (TGTAAAACGACGGCCAGT) forward or M13R (CAGGAAACAGCTATGAC) extensions were added to the 5' ends of the forward and reverse PCR primers, respectively, to allow uniform sequencing conditions. PCR and sequencing reactions were carried out as previously described . Primers and conditions used in PCR reactions are listed in Additional file 1. The quality of sequencing reads was assessed using Phred [38, 39], potential variants identified by Polyphred version 5  and all sequences assembled with reference sequences using Phrap  and viewed in Consed version 12 .
Haplotypes of variants with minor allele frequency (MAF) >5% in the sequence data were inferred using PHASE v2.1.1 [43, 44]. Four tagSNPs were selected for each gene using TagSNP, version 1.1 . Three additional SNPs of potential functional relevance in NBS1 were also tested.
TaqMan® was used for all genotyping. Assays were designed using the Assays-by-Design service (Applied Biosystems). Primers and probes used are listed in Additional file 2. 10 ng of each sample was aliquoted in 384-well plates and the DNA dried down at room temperature. TaqMan reactions were carried out in 5 uL volumes as per the manufacturer's protocols. Fluorescence data was obtained in the ABI PRISM 7900 HT, after 10 min at 95°C, followed by 40 cycles of 92°C for 15 s and 60°C for 1 min. The SDS2.2 software (Applied Biosystems) was used to assign genotypes to individual samples.
Statistical analyses were carried out as described previously . Briefly, all controls were tested for deviation from Hardy-Weinberg equilibrium. Odds ratios (OR) and 95% confidence intervals were estimated using logistic regression. These analyses were conducted using SPSS version 15, with adjustment for sex, age group (categories: 20-49, 50-59, 60-69, 70+), residence (Vancouver or Victoria), and for ethnicity (Caucasian, Asian, South Asian, Mixed, Unknown/Refused) when all cases and all controls were analyzed together. Heterozygotes and rare homozygotes were combined for analysis when the number of rare homozygotes was less than five. Tests were not performed when the sum of the number of heterozygotes and rare homozygotes was less than five for cases or controls. Tests for trend were conducted when there were at least five samples in each genotype category for both cases and controls. Multiple testing correction was carried out by the false discovery rate (FDR) method . Because we tested nineteen markers, the p-value of the most significant marker must be below the threshold of 0.0026 to be considered significant. The haplotypes inferred were analyzed as categorical variables and assessed for risk effect using R version 2.1.1 . Haplotypes with frequency <4.5% were combined into a "rare" category.
Re-sequencing for variant discovery
We sequenced DNA samples from 87 NHL cases to survey the germline genetic variation in the NBS1, MRE11, RAD50 and BLM genes in NHL patients in our population. By using recent methods  the number of unseen variants using data from deep sequencing projects (such as ENCODE ) can be estimated. Using such methods, the sequencing of 174 chromosomes in our population is expected to have revealed 99.99% of SNPs with a MAF of 1% or more, and 76% of SNPs with a MAF of 0.5% or more. Samples were derived from 74 cases with B-cell NHL and 13 with T-cell NHL (see Additional file 3). The number of amplicons bi-directionally sequenced for each gene is shown in Table 2. In total, 63 amplicons were used. On average, 91.6% of sample-amplicon combinations produced good quality reads in both directions, and 96.6% of sample-amplicon combinations produced good quality reads in at least one direction.
Re-sequencing revealed 114 variants (Additional file): 12 small deletions or insertions (10.5%), 73 (64%) transitions and 29 (25.4%) transversions. Twenty-nine variants (25.4%) were in coding regions, with 17 (58.6%) non-synonymous mutations, 4 of which were ranked as "probably" or "possibly damaging" by PolyPhen . Only one of these, BLM_X13_(2603)_C/T, was observed more than once, with a MAF of 5.6%. Fifty-five (48%) variants were "singletons", meaning the minor allele was only observed once in this data set of 87 samples, or 174 chromosomes. Forty-one (36%) variants were "common", with MAF of at least 5%. 59% of variants were previously described in dbSNP (build 128) ; their rs numbers are included in Additional file 4. Of the common polymorphisms (MAF ≥5%), 14% were novel.
Overall, sequence variation was found at 34 of 12,352 nucleotides in coding regions (or 8 of 3,805 nucleotides in RAD50, 13 of 2,265 nucleotides in NBS1, 2 of 2,127 nucleotides in MRE11, and 11 of 4,155 nucleotides in BLM) and at 95 of 17,257 nucleotides in non-coding regions (or 21 of 5,226 nucleotides in RAD50, 32 of 3,777 nucleotides in NBS1, 20 of 3,872 nucleotides in MRE11, and 22 of 4,382 nucleotides in BLM). The K a/K s value for these four genes together is 0.6 (or 0.56 for RAD50, 0.75 for NBS1, 0.50 for MRE11, and 0.54 for BLM), indicating moderate negative selection.
Linkage Disequilibrium (LD) calculations were performed in sequence data using Haploview v4.0 ; singletons were excluded from these calculations. r 2 values for pairwise combinations of SNPs in each gene are shown in Additional files 5, 6, 7 &8.
Haplotypes were inferred using the 41 variants that were observed more than once in the sequence data, using PHASE v2.1.1 [43, 44]. The number of haplotypes inferred for each gene is indicated in Table 2. Haplotype tagging SNPs (tagSNPs) were selected using TagSNP version 1.1 . Nineteen variants were chosen for genotyping and are indicated in bold in Additional file 4.
The 19 tagSNPs were genotyped in 797 cases and 793 controls, with an average genotype call rate of 97.6%. Their respective MAFs, as calculated using all 1590 samples, are in Additional file 2. The concordance of genotypes (in the 87 samples that were sequenced) between the independent methods of sequencing and TaqMan genotyping was complete; no discrepancies were found. As a quality assurance measure, we also genotyped the 19 SNPs in DNA samples from five three-generation CEPH families (purchased from Coriell Cell Repositories, NJ, USA) and confirmed that the alleles segregated according to Mendelian inheritance.
NHL association tests
We compared all European ancestry controls against all European ancestry NHL cases, all B-cell NHL, all T-cell NHL and major subtypes individually. One of the variants, MRE11_5UP_(-1456)_C/T, was excluded from analysis due to deviation from Hardy-Weinberg equilibrium in controls. Results for the two most common subtypes - diffuse large B-cell lymphoma (DLBCL) and follicular lymphoma (FL) - and results suggestive of association with Marginal Zone lymphoma/Mucosa-Associated Lymphoid Tissue (MZ/MALT) are shown in Table 3; see Additional file 9 for all results. RAD50_IVS22(+24)_A/G showed a possible association with DLBCL that was strong enough to influence the overall NHL analysis (p-trend of 0.022 for DLBCL). Another example is RAD50_IVS7(-38)_C/T in MZL, with an OR of 3.39 (95% CI: 1.48-7.75, p = 0.004).
Analyses for all NHL were performed separately for the Asian and South Asian cases (see Additional file 10). NBS1_3UTR_(+273)_G/A (rs1063053) gave an OR of 5.3 (95% CI = 1.023 - 27.579, p = 0.004) in samples of South Asian ethnicity.
Combined analyses of all samples from all ethnicities were also performed, adjusting for ethnicity in the model (data not shown); some SNPs (usually the same as in the European ancestry only analysis) again showed results suggestive of association but failed to reach p < 0.05 upon correction for multiple testing. The ethnic diversity of our study population could mask a real signal and so we focused on the European subpopulation.
The haplotypes inferred from individual SNP genotypes were also tested for association with NHL using R version 2.1.1 (data not shown). No haplotype was more significantly associated with NHL than the individual SNPs forming that haplotype.
RAD50, NBS1, MRE11 and BLM were re-sequenced in 87 NHL cases to characterize the variation in these genes in NHL cases in our population. All genes had similar numbers of variants and similar nucleotide diversity, albeit slightly greater for NBS1 (Table 2). All four genes showed evidence of negative selection, as indicated by a K a/K s value of less than one (0.56 for all four genes combined), which we would expect for genes involved in such a conserved and critical process as DNA repair. The most variable gene, NBS1, also showed the lowest conservation.
Two SNPs in RAD50 were suggestive of association with specific NHL subtypes (Table 3). RAD50_IVS7(-38)_C/T was suggestive of association with MZ/MALT (p = 0.004). The low frequency of this allele (MAF 2.6%), and the low incidence of MZL/MALT (12% of our cases) make it difficult to conclusively implicate this marker in a single study. Interestingly, MZL lymphomas usually develop in tissue subjected to chronic antigenic stimulation, for example gastric MALT lymphoma which arises as a result of chronic Helicobacter pylori infection. Such tissue, with persistent and accelerated cell lymphoid cell proliferation, may be uniquely susceptible to neoplastic transformation associated with faulty DNA repair. Our results may serve to highlight specific mechanistic hypotheses for further testing in other association studies, or for in vitro functional studies. Mechanisms of tumourigenesis, and the basis for NHL susceptibility, may differ between NHL subtypes. Observations such as ours, if replicable, will help us understand the basis for the diversity of NHL types.
We did not find that variants in NBS1 conferred an increased risk of lymphoma, as in most other studies [53–57], although there remain some contradictory positive reports [58–61]. In contrast, non-synonymous mutations in NBS1 have been shown to be associated with acute lymphoblastic leukemia in German  and Polish  children. A study by Rollinson et al  of haplotypic variation in NHL found no increased risk associated with haplotypes of NBS1 and RAD50; however, they observed the variant rs601341 in MRE11 to have a protective effect on FL and a protective effect of an MRE11 haplotype on DLBCL. We did not sequence the part of intron 18 where rs601341 is located and so did not explicitly test this SNP. The difference between our results and those of Rollinson et al. could be the result of a SNP-specific effect, and/or the different populations studied.
Although there have been other studies of susceptibility to NHL looking at the genes addressed in this study, most have relied on the genotyping of rare variants discovered in studies of the rare recessive syndromes discussed above. Genotyping was generally done using single-strand conformation polymorphisms [17, 53, 54, 56, 58, 61, 62] or by TaqMan . One study  used public databases to collect the information on the SNPs in the regions of interest. However, sequencing of germline DNA of patients with sporadic lymphoma to systematically identify genetic variants had not been previously done. Our systematic characterization of these genes provides valuable information on the variation found in these genes in individuals with NHL. Previous systematic investigations of another double-stranded break repair gene, ATM, by our group did not reveal any association between common variants in ATM and NHL or its subtypes . In contrast, a common SNP in the promoter region of H2AX showed a protective effect on NHL and on FL in particular .
Limitations of our study include the histological heterogeneity of NHL, which is composed of many subtypes, many of which are rare. Identification of genetic susceptibility factors that differ between subtypes will be limited by the lack of availability of adequate sample numbers for less common subtypes. The clinical diversity of NHL enabled us to make the strongest conclusions only for DLBCL and FL. Our sample is also ethnically heterogeneous, and so has reduced power to detect genetic factors that are present only in specific ethnic groups. Future replication of results in the context of large international consortia, such as the InterLymph Consortium , will help to overcome such limitations.
While the genes in this study were not significantly associated with NHL independently, it is possible that they could modify NHL risk in combination with other variants. Larger studies would be required to detect such gene-gene interactions. Our observation of possible associations of SNPs in RAD50 with DLBCL and MZ/MALT lymphomas may contribute to the refinement of biological hypotheses for confirmation in larger association studies and functional studies. Mechanisms of tumourigenesis, and the basis for NHL susceptibility, probably differ between NHL subtypes. Specific observations such as these will help us understand the etiological basis for the diversity of NHL types.
diffuse large B-cell lymphoma
- MRN complex:
single nucleotide polymorphism
Nijmegen breakage syndrome
confidence interval, CNS: conserved non-coding sequence, PCR: polymerase chain reaction, MZ: marginal zone lymphoma, MALT: mucosa-associated lymphoid tissue
Amercian_Cancer_Society_document_"Non-Hodgkin_Lymphoma". 2008, [http://www.cancer.org/docroot/LRN/LRN_0.asp?dt=32]
Canadian_Cancer_Statistics. 2008, [http://www.cancer.ca]
Tsujimoto Y, Finger LR, Yunis J, Nowell PC, Croce CM: Cloning of the chromosome breakpoint of neoplastic B cells with the t(14;18) chromosome translocation. Science (New York, NY). 1984, 226 (4678): 1097-1099.
Novik KL, Spinelli JJ, Macarthur AC, Shumansky K, Sipahimalani P, Leach S, Lai A, Connors JM, Gascoyne RD, Gallagher RP, et al: Genetic variation in H2AFX contributes to risk of non-Hodgkin lymphoma. Cancer Epidemiol Biomarkers Prev. 2007, 16 (6): 1098-1106. 10.1158/1055-9965.EPI-06-0639.
De la Torre C, Pincheira J, Lopez-Saez JF: Human syndromes with genomic instability and multiprotein machines that repair DNA double-strand breaks. Histology and histopathology. 2003, 18 (1): 225-243.
Carney JP, Maser RS, Olivares H, Davis EM, Le Beau M, Yates JR, Hays L, Morgan WF, Petrini JH: The hMre11/hRad50 protein complex and Nijmegen breakage syndrome: linkage of double-strand break repair to the cellular DNA damage response. Cell. 1998, 93 (3): 477-486. 10.1016/S0092-8674(00)81175-7.
Hopfner KP, Karcher A, Craig L, Woo TT, Carney JP, Tainer JA: Structural biochemistry and interaction architecture of the DNA double-strand break repair Mre11 nuclease and Rad50-ATPase. Cell. 2001, 105 (4): 473-485. 10.1016/S0092-8674(01)00335-X.
Thompson LH, Schild D: Recombinational DNA repair and human disease. Mutation research. 2002, 509 (1-2): 49-78.
Burgt van der I, Chrzanowska KH, Smeets D, Weemaes C: Nijmegen breakage syndrome. Journal of medical genetics. 1996, 33 (2): 153-156. 10.1136/jmg.33.2.153.
Seemanova E, Passarge E, Beneskova D, Houstek J, Kasal P, Sevcikova M: Familial microcephaly with normal intelligence, immunodeficiency, and risk for lymphoreticular malignancies: a new autosomal recessive disorder. American journal of medical genetics. 1985, 20 (4): 639-648. 10.1002/ajmg.1320200410.
Chrzanowska KH, Kleijer WJ, Krajewska-Walasek M, Bialecka M, Gutkowska A, Goryluk-Kozakiewicz B, Michalkiewicz J, Stachowski J, Gregorek H, Lyson-Wojciechowska G, et al: Eleven Polish patients with microcephaly, immunodeficiency, and chromosomal instability: the Nijmegen breakage syndrome. American journal of medical genetics. 1995, 57 (3): 462-471. 10.1002/ajmg.1320570321.
Heikkinen K, Rapakko K, Karppinen SM, Erkko H, Knuutila S, Lundan T, Mannermaa A, Borresen-Dale AL, Borg A, Barkardottir RB, et al: RAD50 and NBS1 are breast cancer susceptibility genes associated with genomic instability. Carcinogenesis. 2006, 27 (8): 1593-1599. 10.1093/carcin/bgi360.
Ranganathan V, Heine WF, Ciccone DN, Rudolph KL, Wu X, Chang S, Hai H, Ahearn IM, Livingston DM, Resnick I, et al: Rescue of a telomere length defect of Nijmegen breakage syndrome cells requires NBS and telomerase catalytic subunit. Curr Biol. 2001, 11 (12): 962-966. 10.1016/S0960-9822(01)00267-6.
Tauchi H, Matsuura S, Kobayashi J, Sakamoto S, Komatsu K: Nijmegen breakage syndrome gene, NBS1, and molecular links to factors for genome stability. Oncogene. 2002, 21 (58): 8967-8980. 10.1038/sj.onc.1206136.
Grenon M, Gilbert C, Lowndes NF: Checkpoint activation in response to double-strand breaks requires the Mre11/Rad50/Xrs2 complex. Nature cell biology. 2001, 3 (9): 844-847. 10.1038/ncb0901-844.
Uziel T, Lerenthal Y, Moyal L, Andegeko Y, Mittelman L, Shiloh Y: Requirement of the MRN complex for ATM activation by DNA damage. The EMBO journal. 2003, 22 (20): 5612-5621. 10.1093/emboj/cdg541.
Varon R, Reis A, Henze G, von Einsiedel HG, Sperling K, Seeger K: Mutations in the Nijmegen Breakage Syndrome gene (NBS1) in childhood acute lymphoblastic leukemia (ALL). Cancer research. 2001, 61 (9): 3570-3572.
Shimada H, Shimizu K, Mimaki S, Sakiyama T, Mori T, Shimasaki N, Yokota J, Nakachi K, Ohta T, Ohki M: First case of aplastic anemia in a Japanese child with a homozygous missense mutation in the NBS1 gene (I171V) associated with genomic instability. Human genetics. 2004, 115 (5): 372-376. 10.1007/s00439-004-1155-1.
Heikkinen K, Karppinen SM, Soini Y, Makinen M, Winqvist R: Mutation screening of Mre11 complex genes: indication of RAD50 involvement in breast and ovarian cancer susceptibility. Journal of medical genetics. 2003, 40 (12): e131-10.1136/jmg.40.12.e131.
Tommiska J, Seal S, Renwick A, Barfoot R, Baskcomb L, Jayatilake H, Bartkova J, Tallila J, Kaare M, Tamminen A, et al: Evaluation of RAD50 in familial breast cancer predisposition. International journal of cancer. 2006, 118 (11): 2911-2916. 10.1002/ijc.21738.
Giannini G, Ristori E, Cerignoli F, Rinaldi C, Zani M, Viel A, Ottini L, Crescenzi M, Martinotti S, Bignami M, et al: Human MRE11 is inactivated in mismatch repair-deficient cancers. EMBO reports. 2002, 3 (3): 248-254. 10.1093/embo-reports/kvf044.
German J: Bloom's syndrome. Dermatologic clinics. 1995, 13 (1): 7-18.
Goss KH, Risinger MA, Kordich JJ, Sanz MM, Straughen JE, Slovek LE, Capobianco AJ, German J, Boivin GP, Groden J: Enhanced tumor formation in mice heterozygous for Blm mutation. Science (New York, NY). 2002, 297 (5589): 2051-2053.
Gruber SB, Ellis NA, Scott KK, Almog R, Kolachana P, Bonner JD, Kirchhoff T, Tomsho LP, Nafa K, Pierce H, et al: BLM heterozygosity and the risk of colorectal cancer. Science (New York, NY). 2002, 297 (5589): 2013-
Cleary SP, Zhang W, Di Nicola N, Aronson M, Aube J, Steinman A, Haddad R, Redston M, Gallinger S, Narod SA, et al: Heterozygosity for the BLM(Ash) mutation and cancer risk. Cancer research. 2003, 63 (8): 1769-1771.
Wang W: Emergence of a DNA-damage response network consisting of Fanconi anaemia and BRCA proteins. Nature reviews. 2007, 8 (10): 735-748.
Shimura T, Torres MJ, Martin MM, Rao VA, Pommier Y, Katsura M, Miyagawa K, Aladjem MI: Bloom's syndrome helicase and Mus81 are required to induce transient double-strand DNA breaks in response to DNA replication stress. Journal of molecular biology. 2008, 375 (4): 1152-1164. 10.1016/j.jmb.2007.11.006.
Dumon-Jones V, Frappart PO, Tong WM, Sajithlal G, Hulla W, Schmid G, Herceg Z, Digweed M, Wang ZQ: Nbn heterozygosity renders mice susceptible to tumor formation and ionizing radiation-induced tumorigenesis. Cancer research. 2003, 63 (21): 7263-7269.
Xiao Y, Weaver DT: Conditional gene targeted deletion by Cre recombinase demonstrates the requirement for the double-strand break repair Mre11 protein in murine embryonic stem cells. Nucleic acids research. 1997, 25 (15): 2985-2991. 10.1093/nar/25.15.2985.
Bender CF, Sikes ML, Sullivan R, Huye LE, Le Beau MM, Roth DB, Mirzoeva OK, Oltz EM, Petrini JH: Cancer predisposition and hematopoietic failure in Rad50(S/S) mice. Genes & development. 2002, 16 (17): 2237-2251. 10.1101/gad.1007902.
Morales M, Theunissen JW, Kim CF, Kitagawa R, Kastan MB, Petrini JH: The Rad50S allele promotes ATM-dependent DNA damage responses and suppresses ATM deficiency: implications for the Mre11 complex as a DNA damage sensor. Genes & development. 2005, 19 (24): 3043-3054. 10.1101/gad.1373705.
Sipahimalani P, Spinelli JJ, MacArthur AC, Lai A, Leach SR, Janoo-Gilani RT, Palmquist DL, Connors JM, Gascoyne RD, Gallagher RP, et al: A systematic evaluation of the ataxia telangiectasia mutated gene does not show an association with non-Hodgkin lymphoma. International journal of cancer. 2007, 121 (9): 1967-1975. 10.1002/ijc.22888.
Spinelli JJ, Ng CH, Weber JP, Connors JM, Gascoyne RD, Lai AS, Brooks-Wilson AR, Le ND, Berry BR, Gallagher RP: Organochlorines and risk of non-Hodgkin lymphoma. International journal of cancer. 2007, 121 (12): 2767-2775. 10.1002/ijc.23005.
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D: The human genome browser at UCSC. Genome research. 2002, 12 (6): 996-1006.
Frazer KA, Pachter L, Poliakov A, Rubin EM, Dubchak I: VISTA: computational tools for comparative genomics. Nucleic acids research. 2004, W273-279. 10.1093/nar/gkh458. 32 Web Server
Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods in molecular biology (Clifton, NJ). 2000, 132: 365-386.
Brooks-Wilson AR, Kaurah P, Suriano G, Leach S, Senz J, Grehan N, Butterfield YS, Jeyes J, Schinas J, Bacani J, et al: Germline E-cadherin mutations in hereditary diffuse gastric cancer: assessment of 42 new families and review of genetic screening criteria. Journal of medical genetics. 2004, 41 (7): 508-517. 10.1136/jmg.2004.018275.
Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome research. 1998, 8 (3): 175-185.
Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome research. 1998, 8 (3): 186-194.
Stephens M, Sloan JS, Robertson PD, Scheet P, Nickerson DA: Automating sequence-based detection and genotyping of SNPs from diploid samples. Nature genetics. 2006, 38 (3): 375-381. 10.1038/ng1746.
Green P: Phrap. 1996, [http://www.phrap.org/phredphrap/phrap.html]
Gordon D: Viewing and editing assembled sequences using Consed. Current protocols in bioinformatics/editoral board, Andreas D Baxevanis [et al]. 2003, Chapter 11 (Unit11): 12-
Stephens M, Donnelly P: A comparison of bayesian methods for haplotype reconstruction from population genotype data. American journal of human genetics. 2003, 73 (5): 1162-1169. 10.1086/379378.
Stephens M, Smith NJ, Donnelly P: A new statistical method for haplotype reconstruction from population data. American journal of human genetics. 2001, 68 (4): 978-989. 10.1086/319501.
Stram DO, Haiman CA, Hirschhorn JN, Altshuler D, Kolonel LN, Henderson BE, Pike MC: Choosing haplotype-tagging SNPS based on unphased genotype data using a preliminary sample of unrelated subjects with an example from the Multiethnic Cohort Study. Human heredity. 2003, 55 (1): 27-36. 10.1159/000071807.
Benjamini Y, Y H: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995, 57: 289-300.
Bates D, Chambers J, Dalgaard P, Falcon S, Gentleman R, Hornik K, Iacus S, Ihaka R, Leisch F, Lumley T, et al: R. The R Foundation for Statistical Computing, 2.1.1
Ionita-Laza I, Lange C, N ML: Estimating the number of unseen variants in the human genome. Proceedings of the National Academy of Sciences of the United States of America. 2009, 106 (13): 5008-5013. 10.1073/pnas.0807815106.
Encode. 2009, [http://www.hapmap.org/downloads/encode1.html.en]
Ramensky V, Bork P, Sunyaev S: Human non-synonymous SNPs: server and survey. Nucleic acids research. 2002, 30 (17): 3894-3900. 10.1093/nar/gkf493.
Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics (Oxford, England). 2005, 21 (2): 263-265. 10.1093/bioinformatics/bth457.
Cerosaletti KM, Morrison VA, Sabath DE, Willerford DM, Concannon P: Mutations and molecular variants of the NBS1 gene in non-Hodgkin lymphoma. Genes, chromosomes & cancer. 2002, 35 (3): 282-286. 10.1002/gcc.10114.
Chrzanowska KH, Piekutowska-Abramczuk D, Popowska E, Gladkowska-Dura M, Maldyk J, Syczewska M, Krajewska-Walasek M, Goryluk-Kozakiewicz B, Bubala H, Gadomski A, et al: Carrier frequency of mutation 657del5 in the NBS1 gene in a population of Polish pediatric patients with sporadic lymphoid malignancies. International journal of cancer. 2006, 118 (5): 1269-1274. 10.1002/ijc.21439.
Rischewski J, Bismarck P, Kabisch H, Janka-Schaub G, Obser T, Schneppenheim R: The common deletion 657del5 in the Nibrin gene is not a major risk factor for B or T cell non-Hodgkin lymphoma in a pediatric population. Leukemia. 2000, 14 (8): 1528-1529. 10.1038/sj.leu.2401836.
Soucek P, Gut I, Trneny M, Skovlund E, Grenaker Alnaes G, Kristensen T, Borresen-Dale AL, Kristensen VN: Multiplex single-tube screening for mutations in the Nijmegen Breakage Syndrome (NBS1) gene in Hodgkin's and non-Hodgkin's lymphoma patients of Slavic origin. Eur J Hum Genet. 2003, 11 (5): 416-419. 10.1038/sj.ejhg.5200972.
Stanulla M, Stumm M, Dieckvoss BO, Seidemann K, Schemmel V, Muller Brechlin A, Schrappe M, Welte K, Reiter A: No evidence for a major role of heterozygous deletion 657del5 within the NBS1 gene in the pathogenesis of non-Hodgkin's lymphoma of childhood and adolescence. British journal of haematology. 2000, 109 (1): 117-120. 10.1046/j.1365-2141.2000.01973.x.
Resnick IB, Kondratenko I, Pashanov E, Maschan AA, Karachunsky A, Togoev O, Timakov A, Polyakov A, Tverskaya S, Evgrafov O, et al: 657del5 mutation in the gene for Nijmegen breakage syndrome (NBS1) in a cohort of Russian children with lymphoid tissue malignancies and controls. Am J Med Genet A. 2003, 120A (2): 174-179. 10.1002/ajmg.a.20188.
Seemanova E, Hoch J, Herzogova J, Kawaciuk I, Janda J, Kohoutova M, Seeman P, Varon R, Sperling K: Mutations in tumor suppressor gene NBS1 in adult patients with malignancies. Casopis lekaru ceskych. 2006, 145 (3): 201-203.
Seemanova E, Jarolim P, Seeman P, Varon R, Sperling K: Increased risk of malignancies in heterozygotes in families of patients with Nijmegen breakage syndrome. Casopis lekaru ceskych. 2006, 145 (2): 138-143.
Steffen J, Varon R, Mosor M, Maneva G, Maurer M, Stumm M, Nowakowska D, Rubach M, Kosakowska E, Ruka W, et al: Increased cancer risk of heterozygotes with NBS1 germline mutations in Poland. International journal of cancer. 2004, 111 (1): 67-71. 10.1002/ijc.20239.
Mosor M, Ziolkowska I, Pernak-Schwarz M, Januszkiewicz-Lewandowska D, Nowak J: Association of the heterozygous germline I171V mutation of the NBS1 gene with childhood acute lymphoblastic leukemia. Leukemia. 2006, 20 (8): 1454-1456. 10.1038/sj.leu.2404285.
Rollinson S, Kesby H, Morgan GJ: Haplotypic variation in MRE11, RAD50 and NBS1 and risk of non-Hodgkin's lymphoma. Leukemia & lymphoma. 2006, 47 (12): 2567-2583. 10.1080/10428190600909743.
Rothman N, Skibola CF, Wang SS, Morgan G, Lan Q, Smith MT, Spinelli JJ, Willett E, De Sanjose S, Cocco P, et al: Genetic variation in TNF and IL10 and risk of non-Hodgkin lymphoma: a report from the InterLymph Consortium. The lancet oncology. 2006, 7 (1): 27-38. 10.1016/S1470-2045(05)70434-4.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/10/117/prepub
This work was supported by the Canadian Cancer Society through grants from the National Cancer Institute of Canada, and by grants from the Canadian Institutes of Health Research (CIHR). It was initiated with funds from the Chan Sisters Foundation and the BC Cancer Foundation. AB-W is a Senior Scholar of the Michael Smith Foundation for Health Research (MSFHR). JMS was supported by scholarships from MSFHR, the National Sciences and Engineering Research Council of Canada, and CIHR. ACM was supported by a research unit infrastructure grant from MSFHR. We thank all the participants of the study for making this research possible. We thank Rozmin Janoo-Gilani, Kuldip Bagga, Agnes Bauzon, Betty Hall, Lina Hsu, Pat Ostrow, Lynne Tse, Anthony Tung, Tim Lee and Zenaida Abanto for their assistance.
The authors declare that they have no competing interests.
JMS carried out some experiments, participated in the coordination and interpretation of the analysis, and drafted the manuscript. AM performed the statistical analysis and assisted in interpretation. SL carried out experiments. AL coordinated sample collection and interpreted phenotypes. JC and RDG interpreted the medical relevance and pathology, respectively. RPG and JJS participated in the design of the study and JJS advised on the statistical analysis. ABW participated in the study design, coordinated the experiments and helped to draft the manuscript. All authors read and approved the manuscript.
Electronic supplementary material
Additional file 1: PCR primers & conditions. For each exon, the primers and temperatures used for each PCR reaction, as well as the product size. (PDF 41 KB)
Additional file 2: Genotyping assays: primers, probes, and minor allele frequencies. Primers and probes used in genotyping assays, with genotyped minor allele frequencies. (PDF 33 KB)
Additional file 3: Samples sequenced. Composition (gender, age and pathology) of the samples used in the re-sequencing phase. (PDF 28 KB)
Additional file 4: SNPs discovered by re-sequencing. Full list of all SNPs found by re-sequencing and their properties (chromosomal position, dbSNP number if known, flanking sequence, nucleotide/codon/amino acid change, and minor allele frequency). (PDF 68 KB)
Additional file 9: Regression analysis in Caucasian samples. Table containing statistical analysis results in Caucasian samples, for overall NHL and all subtypes examined, with the following sections: Additional file 9a: Regression analysis results for RAD50 SNPs in Caucasian samples. Additional file 9b: Regression analysis results for NBS1 SNPs in Caucasian samples. Additional file 9c: Regression analysis results for MRE11 SNPs in Caucasian samples. Additional file 9d: Regression analysis results for BLM SNPs in Caucasian samples. (PDF 53 KB)
Additional file 10: Regression analysis for overall NHL in Asians and South-East Asian samples. Table containing statistical analysis results for overall NHL only. (PDF 20 KB)
About this article
Cite this article
Schuetz, J.M., MacArthur, A.C., Leach, S. et al. Genetic variation in the NBS1, MRE11, RAD50 and BLM genes and susceptibility to non-Hodgkin lymphoma. BMC Med Genet 10, 117 (2009). https://doi.org/10.1186/1471-2350-10-117
- Minor Allele Frequency
- Follicular Lymphoma
- Nijmegen Breakage Syndrome
- Bloom Syndrome
- Rare Homozygote