MUTYH gene polymorphisms as risk factors for rheumatoid arthritis
Yung-Jen Kung1, Kun Shi Tsai2, Chung-Ming Huang3,4, Hui-Ju Lin5, Ter-Hsin Chen6, Yu-An Hsu7, Ching-Yao Chang8,9, Yong-San Huang1, and Lei Wan9, 10, 11*
1Department of Veterinary Medicine, National Chung Hsing University; 2 Department
of Biological Science and Technology, National Chiao Tung University, Hsinchu, Taiwan; 3Graduate Institute of Integrated medicine, China Medical University, Taichung, Taiwan; 4Division of Immunology and Rheumatology, Department of Internal Medicine, China Medical University Hospital, Taichung, Taiwan; 5Department of Ophthalmology, China Medical University Hospital, Taichung,
Taiwan; 6Graduate Institute of Veterinary Pathobiology, National Chung Hsing University, Taichung, Taiwan; 7Department of Life Science, National Tsing Hua University, HsinChu, Taiwan; 8 Department of Medical Research, China Medical University Hospital, Taichung, Taiwan; 9Department of Biotechnology, Asia University, Taichung, Taiwan; 10School of Chinese Medicine, China Medical University, Taichung, Taiwan; 11Department of Gynecology, China Medical University Hospital, Taichung, Taiwan
Yung-Jen Kung and Kun Shi Tsai contributed equally to the work and therefore should be considered equivalent authors.
* Reprints and correspondence to: Lei Wan, PhD.
School of Chinese Medicine, China Medical University, Taichung, Taiwan No. 91, Hsueh-Shih Road, Taichung 40402, Taiwan.
Telephone: +886-4-22053366 ext. 3326
Abstract
Objectives: MUTYH glycosylase directly interacts with various proteins involved in DNA repair pathways and may be associated with the risk of
autoimmune diseases such as Rheumatoid arthritis (RA). Therefore, the association between polymorphisms in the MUTYH gene and RA was evaluated.
Methods: We recruited 192 RA patients and 192 healthy subjects in Taiwan. The 4 MUTYH polymorphisms (rs3219463, rs3219476, rs3219489, and rs3219493) were detected by a polymerase chain reaction-restriction fragment length polymorphism technique. Haplotype analysis was performed using unphased genotype data using the Bayesian method. The genotype and allelic frequency distributions of the
polymorphisms in both RA patients and healthy patients were compared by the chi-square test.
Results: Comparison of the genotype/allele frequencies between individuals with RA and the control groups revealed significant differences in 2 MUTYH gene
polymorphisms, rs3219463 and rs3219476. The haplotype frequencies were estimated among the 4 identified polymorphisms. After we performed a haplotype-specific analysis, the haplotypes Ht6-GTGC and Ht8-GGCG had lower presenting rates in RA patients than in the control groups. The association between the clinical features of
RA patients with genotypes rs3219463 and rs3219476 were analyzed. The genotype frequency of rs3219463 G/- was significantly increased among patients with
immunoglobulin M rheumatoid factors, whereas that of rs3219476 was not. There was no significant difference between Ht6-GTGC and Ht8-GGCG among RA patients.
Conclusion: We demonstrated that the rs3219463 and rs3219476
polymorphisms in RA patients from a Taiwan Chinese population were associated with disease susceptibility. These data indicate that the MUTYH gene may play a role in the progression of RA.
1. Introduction
Rheumatoid arthritis (RA) is a chronic autoimmune disease that can lead to
deformities and severe disability because of irreversible damage to the tendons, joints, and bones. Although RA is a common type of arthritis with a prevalence of 1%, the pathogenesis is still unknown[1, 2]. RA can be considered to be a complex genetic disease and many “disease-risk” genes or “disease-protective” genes may be involved[3, 4].
Numerous studies have reported an association between HLA-DR gene
polymorphisms[4] and RA, and CD4 T cells play important roles in the development of RA[5]. In addition, the chronic inflammatory process of RA is mediated through the cytokines network, which induces the expression of some destructive enzymes, such as matrix metalloproteinases in the cellular synovial tissue of joints[1]. Although there is no conclusive evidence, there are hints that genes and enzymes of DNA repair may be involved in some types of autoimmune diseases. It had been suggested that the T cells in RA patients fail to producesufficient transcripts and
proteins of the DNA repair kinaseataxia telangiectasia mutated (ATM) gene[6]. Therefore, we are interested in whether there are other mechanisms of DNA repair related to the etiology of RA. It has been reported that
p53 mutations have been observed in micro-dissected RA synovial tissue[7]. This strongly suggests that somatic DNA mutations induced by oxidative stress enhance the transformation of cells involved in the invasive excess growth of the RA joints and that a loss of genomic stability is one of the key mechanisms related to the severity of RA. Additionally, because inflammatory responses are the main characteristics of RA[4], we have focused on the deficient DNA repair enzymes related to severe or chronic inflammation.
Inflammation can elicit a variety of cellular signals and induce the formation of reactive oxygen species (ROS), which threaten the integrity of cellular DNA. The oxidized base 7,8-dihydroxy-8-oxoguanine (8-oxo-G) is one of the mutagenic products of oxidative DNA damage; however, misincorporated 8-oxo-G can be excised by the base excision repair (BER) pathway[8]. MUTYH glycosylase, which is encoded by MUTYH (mutY Homolog (Escherichia coli)), is involved in this pathway. When incorrect pairing, oxo-G·A, occurs, MUTYH glycosylase intercepts the 8-oxo-G·A base pair and removes the inappropriate A base and provides an opportunity for creation of an 8-oxo-G·C substrate for repair by 8-oxo-G glycosylase (OGG1)[9-11]. MUTYH glycosylase directly interacts with various proteins involved in other DNA repair pathways[12], and defective MUTYH function is associated with various
forms of cancer, including lung cancer[13], gastric cancer[14], head and neck squamous cell carcinomas[15], and colorectal cancer (CRC)[16]. Because there is evidence which suggests that ROS-induced DNA damage contributes to tumor development in a mouse model of chronic inflammation[9], we sought to determine whether mutations in MUTYH are risk factors of RA.
We hypothesized that genetic variations in MUTYH (in its promoter region, introns, or exons) may be associated with the risk of
autoimmune diseases such as RA. Therefore, we analyzed the genetic and allelic frequencies of 4 polymorphisms in MUTYH and the haplotypes among RA patients and healthy controls in the Taiwan Chinese population. We also compared the genotypes and haplotypes among RA patients with various clinical variables to investigate whether a relationship exists between MUTYH polymorphisms and the clinical manifestation of RA.
2. Methods 2.1. Patients
The study subjects including a total of 192 patients with RA and 192 healthy subjects were recruited from China Medical University Hospital in Taiwan. Patients with RA according to the revised America College of Rheumatology criteria were enrolled[17]. Nephelometry was used to detect rheumatoid factor (RF). Values 30 IU/ml were ≧ defined as positive. The presence or history of extra-articular manifestations in patients with RA was recorded[18]. The healthy control from the general population was selected from health examination. All individuals’ samples were collected by venipuncture for genomic DNA isolation. Informed consent was from all participants and was approved by the local Ethics Committee.
2.2. Genomic DNA extraction and genotyping
Genomic DNA was prepared from peripheral blood according to standard protocols of the DNA extraction kit (Qiagen, Valencia, CA, USA). The four MUTYH single nucleotide polymorphisms (SNPs) (rs3219463, rs3219476, rs3219489, and
rs3219493) were detected by using a polymerase chain reaction-restriction fragment length polymorphism technique (PCR-RFLP). PCRs for MUTYH gene
polymorphisms were carried out in a 50-μL reaction mixture containing 50 ng of genomic DNA, 2 to 6 pmole of each primers, 1× Taq polymerase buffer (1.5 mM
MgCl2), and 0.25 U of AmpliTaq DNA polymerase (Applied Biosystems). The primers, PCR conditions, and restriction enzyme cutting sites used to determine MUTYH polymorphisms were listed in Table 1.
2.3. Haplotype analysis
Haplotypes were inferred from unphased genotype data using the Bayesian statistical method available in the software program Phase 2.1[19, 20]. All four SNPs were analyzed with the Phase 2.1 software and the population data were divided into 12 groups, from Ht1 to Ht12. A linkage disequilibrium (LD) map was constructed using Haploview software (version 4.2;
http://www.broadinstitute.org/scientific-community/science/programs/medical-and-population-genetics/haploview/haploview)
to estimate the correlation coefficient between those SNPs.
2.4. Statistical analysis
The chi-square test was utilized to evaluate the Hardy-Weinberg equilibrium in genotypic distributions. Also, the genotype frequency and allelic frequency distributions of the polymorphisms in both RA patients and healthy subjects were compared by the chi-square method (2 × 3 table for genotype frequency; 2 × 2 table for allele frequency) by using the SPSS Version 10.0 software (SPSS Inc., Chicago, IL, USA). Results were considered significant when the P values were < 0.05. The
Bonferroni method was used to correct for multiple testing. The odd ratios (OR) were also calculated from the genomic and allelic frequencies, with a 95% confidence interval (95% CI).
3. Results
3.1. Demographic characteristics of control subjects and RA patients.
A case and control association study was performed. As shown in the supplementary table1, there are no significant differences in gender and age between control subjects and RA patients. Within the 192 RA patients, 151 people were female and the mean age was 53.4 years. In the healthy controls, 151 people were female and the mean age was 50.5.
3.2. Allele and genotype frequency of MUTYH polymorphisms
Of the four tested SNPs, three SNPs was found to be in Hardy–Weinberg equilibrium (P value > 0.05); however, the genotype distribution of the rs3219476 was not
compatible with Hardy–Weinberg equilibrium. The genotype, allele, dominant and recessive distributions of the 4 polymorphisms are shown in Table 2.
Comparison of the genotypes between individuals with RA and the control groups revealed significant differences in the frequency of 2 MUTYH gene polymorphisms, rs3219463 and rs3219476. The P values were corrected for multiple comparison by Bonferroni correction, represented as c P value . For the rs3219463 G/A
polymorphism, the c P value was 0.00024 and the odds ratio (OR) was 1.16 (95% confidence interval [CI] 0.62–2.17) for the heterozygous mutant G/A, and 0.41 (95%
CI 0.21–0.82) for the homozygous mutant G/G. There was also a significant difference for the polymorphism rs3219476 T/G ( c P = 0.0036): the OR was 2.12 (95% CI 1.26–3.59) for the heterozygous mutant T/G and 3.81 (95% CI 1.81–8.04) for the homozygous mutant T/T. However, there was no significant difference between the RA patients and control individuals for rs3219489 or rs3219493. The differences in allele frequencies of these polymorphisms between individuals with and without RA were similar to the results of the genotypes frequencies (Table 2). There were significant differences in allele frequencies for the rs3219463 ( c P = 0.0064, OR 0.63, 95% CI 0.47–0.84) and rs3219476 polymorphism ( c P = 0.0096, OR 1.56, 95% CI 1.17–2.07), while the allele frequencies of the rs3219489 and rs3219493
polymorphisms were not significantly different ( c P = 0.318 and 2.684, respectively). The dominant and recessive distributions of these two RA-related SNPs also
presented significant differences of the frequencies, except for dominant distribution of rs319463. These results demonstrate a significant difference in the allele, genotype, dominant and recessive distribution for the rs3219463 and rs3219476 polymorphisms between the RA and control groups.
3.3. Distribution of MUTYH haplotypes
LD map was constructed (Figure 1) and haplotype frequencies differences were estimated for the 4 identified polymorphisms. The pattern of LD in the MUTYH locus was measured by D’ score among the 4 SNPs. Linkage disequilibrium indicated that rs3219493 and rs3219489 had a high-LD score (D’ > 80), and so did rs3219493 and rs3219476. However, rs3219463 and rs3219476 had a low-LD score (D’ = 61).
Of the 12 observed haplotypes, haplotypes with frequency greater than 5% were presented (Table 3) and were further analyzed. After performing a haplotype-specific analysis, the haplotypes Ht6-GTGC and Ht8-GGCG demonstrated significantly lower frequencies in RA patients than in the control groups (Ht6-GTGC, P = 0.0076, OR 0.41, 95% CI 0.21–0.81; Ht8-GGCG, P = 0.0070, OR 0.32, 95% CI 0.13–0.76), while after Bonferroni correction, the differences were not significant. There were no differences in haplotype frequency between the RA and control groups for the other haplotypes.
3.4. MUTYH polymorphisms and clinical features of RA
The association between the clinical features of patients with RA with the genotypes rs3219463 and rs3219476 were analyzed (Table 4). Among the patients with RA, the genotype frequency of rs3219463 G/- was significantly greater in patients with immunoglobulin M rheumatoid factors (RF) (P = 0.019, OR 2.70, 95% CI 1.155–
6.312). However, a comparison of the genotype frequency of rs3219476 and the clinical features revealed no significant difference among RA patients. There was no significant association between the haplotypes Ht6-GTGC and Ht8-GGCG and the clinical features among RA patients (Table 5).
4. Discussions
In the present study, we found two SNPs of the MUTYH gene were associated with RA. The genetic/allelic frequencies of rs3219463 (in the promoter region) and rs3219476 (in an intron) showed significant differences between RA and control patients. Our results provide the evidence for genetic association coordinated by these polymorphisms with the clinical features of RA. There have been some reports that showed the relationship between the DNA repair pathway and autoimmune diseases. For example, dysfunctional T cells show chronic inflammatory immune responses in the synovium of RA patients, and the DNA repair kinase ATM may be involved [6] . Furthermore, it has been postulated that SNPs in the X-ray repair cross complementing gene 1 ( XRCC1 ), one of the BER proteins, may correlate with RA [21] . To our knowledge, this study is the first investigation to demonstrate that gene
polymorphisms in MUTYH , which encodes a glycosylase member of BER, may contribute to the pathogenesis of RA.
ROS-associated base damage can be repaired by the BER pathway, which is initiated by many damage-specific glycosylases that cleave the glycosylic bond between the damaged base and the sugar
phosphate backbone[8]. MUTYH, which removes adenine paired with 8-oxo-G, is important in suppressing G·C to T·A transversions and in preventing CRC in people and mice. In murine studies, MUTYH-deficient mice were more susceptible to spontaneous tumors [22] and oxidative stress-induced intestinal tumors. MUTYH-associated polyposis (MAP), an autosomal recessive disorder, is associated with biallelic germline mutations in
MUTYH[11, 23] and was first reported in 2002. Since then, many aspects of MAP have been investigated. The possibility that some mutations are nonpathogenic polymorphisms that are found
coincidentally in patients with colorectal adenomas and carcinomas cannot be ruled out[8].
RA, a chronic inflammatory disease with tissue-destructive potential, is now
recognized as a complex genetic disorder. Although numerous studies have reported an association between RA and HLA-DRB1*04, cytokines, and CD4 cells, other genes may contribute to disease susceptibility[4]. Disease risk genes may increase the risk of pathological reactions, which suggests a gene–dose effect. One effect involves a snow-balling mechanism that leads to increased levels of
cellular damage and death resulting in more inflammation that induces the production of more ROS[9]. Several SNPs of various genes were found related with the risk of RA among Taiwanese or Chinese population. For example, an SNP (rs7135617) located in the intron region of ORAI1 gene was associated with a risk of RA in the Taiwanese population [24] . The polymorphism of FCGR2B was related with early-onset of RA and the polymorphism of FCGR3A may influence RF production [25, 26] . Furthermore, there were many studies discussing the IL-family were associated with
pathogenesis of RA in Chinese population, such as IL-4, IL-6 , and IL-18 [27, 28] . The question remains whether defects in DNA damage repairrenders individuals susceptible to RA or whether this is a consequenceof the disease. The accumulation of DNA damage may also have broader implicationson impairing diverse cellular functions.
In haplotype studies, 2 haplotypes of the MUTYH gene showed differences between the RA and control groups (Table 3). The haplotypes Ht6-GTGC and Ht8-GGCG had lower presenting rates in RA patients than in the control groups; these appeared to be
“protective” haplotypes. The rs3219463 polymorphism positioned in the MUTYH gene promoter region; the rs3219476 and rs3219493 polymorphisms positioned in the MUTYH gene intron region; and rs3219489 polymorphism positioned in the MUTYH gene exon region (Gln338His). We found that investigation of the haplotype of the MUTYH gene gave an additional sign for determining whether this gene is associated with RA.
Previous studies have shown that patients with erosive RF-positive joint disease are characterized by the genotype HLA-DRB1*4/RA non-associated allele[4]. In this study, we found that the G/- frequency of rs3219463 was significantly increased in patients with RF (Table 3). Therefore, we postulated that rs3219463 may be a risk factor for RF in RA patients. How the RF factor contributes to the disease process is unknown, but it has been known to facilitate antigen presentation. Immune complexes of RF and immunoglobulin G may contribute to disease activity by activating
complements and stimulating cytokine synthesis[1]. From a clinical perspective, the dissection of the phenotypic and genotypic variants of RA is critical for further exploration of targeted therapeutic treatment.
There were two major limitations in this study. Firstly, because of the
non-compatibility in Hardy-Weinberg equilibrium of rs3219476 genotype distributions, we minimized errors in genotyping and validated our findings by repeating the genotyping analysis several times and obtaining consistent results. Therefore, occurrence of genotyping errors in this study was kept to a minimum. Moreover, the T/G genotype compared with TT genotype of rs3219476 has been revealed as a protective genotype of cholangiocarcinoma among individuals of Han Chinese living in Jiangsu province in China [29] . Based on the genotype frequencies in the controls, the
P value is 0.0859, which is not in Hardy-Weinberg equilibrium. Compared to our
data, we regard the occurrence frequency of rs3219476 is consistent. Furthermore, we
tried to trace the ancestry background of controls and patients, and found that study subjects were a mixed population of Minnan, Hakka, and Canton descendants. According to the paper published by Pan et al.[30], they mentioned that the SNP profiles in the major histocompatibility complex (MHC) region (6p21.3) showed no significant difference among Minnan descendants, Hakka descendants and mixed population of Minnan and Hakka descendants, which indicate the homogeneity of the population. Thus population stratification should not produce in this study[30].
many comparisons revealed no significant differences after Bonferroni correction. Besides, neither the genotype of rs3219463 and rs3219476 nor the haplotypes were significantly associated in patients with extra-articular RA erosion. However, we had evaluated the power in this sample size. For rs3219463 and rs3219476, the statistical powers were 72.1% and 69.9%, respectively. And these two SNPs represent low linkage disequilibrium (Figure 1). Thus, the MUTYH polymorphisms did associated with susceptibility of RA. Even so, our data may be biased by the relatively small number of subjects and larger scale studies are needed to confirm our findings.
In conclusion, our findings demonstrate that the rs3219463 and rs3219476 polymorphisms in the MUTYH gene in patients with RA in a Taiwanese population were associated with disease
susceptibility. Although the exact function of the MUTYH polymorphisms cannot be confirmed in this study, we suggest that genetic variations in the MUTYH gene may affect BER efficiency and play a role in the progression of RA.
Acknowledgements
This study was supported by a grant from the National Science Council (NSC101-2320-B-039-038 and NSC102-2320-B-039-043-MY3), Taipei, Taiwan, grants from
the China Medical University Hospital (DMR-103-064 and CMU102-ASIA-09), Taichung, Taiwan, and a grant in part by Taiwan Ministry of Health and Welfare Clinical Trial and Research Center of Excellence (DOH102-TD-B-111-004).
Conflicts of interest
Reference
1. M. Feldmann, F. M. Brennan and R. N. Maini, "Rheumatoid arthritis," Cell, vol. 85, no. 3, pp. 307-310, 1996.
2. M. Ziff, "Rheumatoid arthritis--its present and future," J Rheumatol, vol. 17, no. 2, pp. 127-133, 1990.
3. A. H. Lynn, C. K. Kwoh, C. M. Venglish, C. E. Aston and A. Chakravarti, "Genetic epidemiology of rheumatoid arthritis," Am J Hum Genet, vol. 57, no. 1, pp. 150-159, 1995.
4. C. M. Weyand, P. A. Klimiuk and J. J. Goronzy, "Heterogeneity of rheumatoid arthritis: from phenotypes to genotypes," Springer Semin Immunopathol, vol. 20, no. 1-2, pp. 5-22, 1998.
5. M. R. Ehrenstein, J. G. Evans, A. Singh, S. Moore, G. Warnes, D. A. Isenberg and C. Mauri, "Compromised function of regulatory T cells in rheumatoid arthritis and reversal by anti-TNFalpha therapy," J Exp Med, vol. 200, no. 3, pp. 277-285, 2004. 6. L. Shao, H. Fujii, I. Colmegna, H. Oishi, J. J. Goronzy and C. M. Weyand,
"Deficiency of the DNA repair enzyme ATM in rheumatoid arthritis," J Exp Med, vol. 206, no. 6, pp. 1435-1449, 2009.
7. Y. Yamanishi, D. L. Boyle, S. Rosengren, D. R. Green, N. J. Zvaifler and G. S. Firestein, "Regional analysis of p53 mutations in rheumatoid arthritis synovium," Proc Natl Acad Sci U S A, vol. 99, no. 15, pp. 10025-10030, 2002.
8. S. S. David, V. L. O'Shea and S. Kundu, "Base-excision repair of oxidative DNA damage," Nature, vol. 447, no. 7147, pp. 941-950, 2007.
9. L. B. Meira, J. M. Bugni, S. L. Green, C. W. Lee, B. Pang, D. Borenshtein, B. H. Rickman, A. B. Rogers, C. A. Moroski-Erkul, J. L. McFaline, D. B. Schauer, P. C. Dedon, J. G. Fox and L. D. Samson, "DNA damage induced by chronic inflammation contributes to colon carcinogenesis in mice," J Clin Invest, vol. 118, no. 7, pp. 2516-2525, 2008.
10. R. C. Niessen, R. H. Sijmons, J. Ou, S. G. Olthof, J. Osinga, M. J. Ligtenberg, F. B. Hogervorst, M. M. Weiss, C. M. Tops, F. J. Hes, G. H. de Bock, C. H. Buys, J. H. Kleibeuker and R. M. Hofstra, "MUTYH and the mismatch repair system: partners in crime?," Hum Genet, vol. 119, no. 1-2, pp. 206-211, 2006.
11. M. L. Poulsen and M. L. Bisgaard, "MUTYH Associated Polyposis (MAP)," Curr Genomics, vol. 9, no. 6, pp. 420-435, 2008.
12. A. R. Parker and J. R. Eshleman, "Human MutY: gene structure, protein functions and interactions, and role in carcinogenesis," Cell Mol Life Sci, vol. 60, no. 10, pp. 2064-2083, 2003.
13. A. Miyaishi, K. Osawa, Y. Osawa, N. Inoue, K. Yoshida, M. Kasahara, A. Tsutou, Y. Tabuchi, K. Sakamoto, N. Tsubota and J. Takahashi, "MUTYH Gln324His gene
polymorphism and genetic susceptibility for lung cancer in a Japanese population," J Exp Clin Cancer Res, vol. 28, pp. 10, 2009.
14. Y. Zhang, X. Liu, Y. Fan, J. Ding, A. Xu, X. Zhou, X. Hu, M. Zhu, X. Zhang, S. Li, J. Wu, H. Cao, J. Li and Y. Wang, "Germline mutations and polymorphic variants in MMR, E-cadherin and MYH genes associated with familial gastric cancer in Jiangsu of China," Int J Cancer, vol. 119, no. 11, pp. 2592-2596, 2006.
15. T. Sliwinski, L. Markiewicz, P. Rusin, W. Pietruszewska, J. Olszewski, A. Morawiec-Sztandera, W. Mlynarski and I. Majsterek, "Polymorphisms of the DNA base excision repair gene MUTYH in head and neck cancer," Exp Oncol, vol. 31, no. 1, pp. 57-59, 2009.
16. H. Chen, L. Xu, Q. Qi, Y. Yao, M. Zhu and Y. Wang, "A haplotype variation affecting the mitochondrial transportation of hMYH protein could be a risk factor for colorectal cancer in Chinese," BMC Cancer, vol. 8, pp. 269, 2008.
17. F. C. Arnett, S. M. Edworthy, D. A. Bloch, D. J. McShane, J. F. Fries, N. S. Cooper, L. A. Healey, S. R. Kaplan, M. H. Liang, H. S. Luthra and et al., "The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis," Arthritis Rheum, vol. 31, no. 3, pp. 315-324, 1988.
18. J. H. Yen, J. R. Chen, W. J. Tsai, J. J. Tsai and H. W. Liu, "HLA-DRB1
genotyping in patients with rheumatoid arthritis in Taiwan," J Rheumatol, vol. 22, no. 8, pp. 1450-1454, 1995.
19. M. Stephens and P. Donnelly, "A comparison of bayesian methods for haplotype reconstruction from population genotype data," Am J Hum Genet, vol. 73, no. 5, pp. 1162-1169, 2003.
20. M. Stephens, N. J. Smith and P. Donnelly, "A new statistical method for
haplotype reconstruction from population data," Am J Hum Genet, vol. 68, no. 4, pp. 978-989, 2001.
21. A. Koyama, Y. Kubota, T. Shimamura and S. Horiuchi, "Possible association of the X-ray cross complementing gene 1 (XRCC1) Arg280His polymorphism as a risk for rheumatoid arthritis," Rheumatol Int, vol. 26, no. 8, pp. 749-751, 2006.
22. K. Sakamoto, Y. Tominaga, K. Yamauchi, Y. Nakatsu, K. Sakumi, K. Yoshiyama, A. Egashira, S. Kura, T. Yao, M. Tsuneyoshi, H. Maki, Y. Nakabeppu and T.
Tsuzuki, "MUTYH-null mice are susceptible to spontaneous and oxidative stress induced intestinal tumorigenesis," Cancer Res, vol. 67, no. 14, pp. 6599-6604, 2007. 23. J. P. Cheadle and J. R. Sampson, "MUTYH-associated polyposis--from defect in base excision repair to clinical genetic testing," DNA Repair (Amst), vol. 6, no. 3, pp. 274-279, 2007.
24. J. H. Yen, C. M. Chang, Y. W. Hsu, C. H. Lee, M. S. Wu, D. Y. Hwang, B. K. Chen, H. T. Liao, M. T. Wu and W. C. Chang, "A polymorphism of ORAI1
rs7135617, is associated with susceptibility to rheumatoid arthritis," Mediators Inflamm, vol. 2014, pp. 834831, 2014.
25. J. Y. Chen, C. M. Wang, C. C. Ma, L. A. Hsu, H. H. Ho, Y. J. Wu, S. N. Kuo and J. Wu, "A transmembrane polymorphism in FcgammaRIIb (FCGR2B) is associated with the production of anti-cyclic citrullinated peptide autoantibodies in Taiwanese RA," Genes Immun, vol. 9, no. 8, pp. 680-688, 2008.
26. J. Y. Chen, C. M. Wang, J. M. Wu, H. H. Ho and S. F. Luo, "Association of rheumatoid factor production with FcgammaRIIIa polymorphism in Taiwanese rheumatoid arthritis," Clin Exp Immunol, vol. 144, no. 1, pp. 10-16, 2006. 27. S. Chen, F. Jiang, J. Ren, J. Liu and W. Meng, "Association of IL-18
polymorphisms with rheumatoid arthritis and systemic lupus erythematosus in Asian populations: a meta-analysis," BMC Med Genet, vol. 13, pp. 107, 2012.
28. X. Li, W. Chai, M. Ni, M. Xu, Z. Lian, L. Shi, Y. Bai and Y. Wang, "The effects of gene polymorphisms in interleukin-4 and interleukin-6 on the susceptibility of rheumatoid arthritis in a Chinese population," Biomed Res Int, vol. 2014, pp. 265435, 2014.
29. S. H. You, X. Wang, S. Huang, M. Wang, G. Z. Ji, J. R. Xia and Z. N. Fan, "MYH rs3219476 and rs3219472 polymorphisms and risk of cholangiocarcinoma," Mol Med Rep, vol. 7, no. 1, pp. 347-351, 2013.
30. W. H. Pan, C. S. Fann, J. Y. Wu, Y. T. Hung, M. S. Ho, T. H. Tai, Y. J. Chen, C. J. Liao, M. L. Yang, A. T. Cheng and Y. T. Chen, "Han Chinese cell and genome bank in Taiwan: purpose, design and ethical considerations," Hum Hered, vol. 61, no. 1, pp. 27-30, 2006.
Figure 1: Four SNPs on the LD map pf MUTYH gene. LD coefficinets (D’) between four SNPs of MUTYH gene were graphically indicated.
