Relationship between Crohn’s disease, infection with Mycobacterium avium subspecies paratuberculosis and SLC11A1 gene polymorphisms in Sardinian patients

Abstract

AIM: To study the association between Crohn’s disease (CD), Mycobacterium avium subspecies paratuberculosis (MAP), and genetic factors by examining the role of natural resistance-associated macrophage protein 1 (NRAMP1) gene polymorphisms (now SLC11A1) in Sardinian patients with CD and controls.

METHODS: Thirty-seven CD patients and 34 controls with no inflammatory bowel disease (IBD) were recruited at the University of Sassari after giving written consent. Six SCL11A1 polymorphisms previously reported to be the most significantly associated with IBD were searched. M. paratuberculosis was identified by IS900 PCR and sequencing. Logistic regression was used to calculate odds ratios (OR) for the associations among CD, presence of MAP, and 6 loci described above.

RESULTS: For the first time, a strong association was observed between polymorphisms at NRAMP1 locus 823C/T and CD. While CD was strongly associated with both NRAMP1 and MAP, NRAMP1 polymorphisms and MAP themselves were not correlated.

CONCLUSION: Combined with previous work on the NOD2/CARD15 gene, it is clear that the interplay of genetic, infectious, and immunologic factors in the etiology of CD is complex.

Key Words: Mycobacterium avium subspecies paratuberculosis; Crohn’s disease; SCL11A1 polymorphisms

INTRODUCTION

Natural resistance-associated macrophage protein 1 (NRAMP1), now strictly referred to as SLC11A1 (Solute carrier 11a1) and the gene which encodes for it is recognized as having a role in the susceptibility of men and animals to a number of mycobacterial infections. In human beings, the NRAMP1 gene is located on human chromosome region 2q35. It is composed of 15 exons and covers at least 16 kb of DNA. It encodes an integral membrane protein of 550 amino acids that is expressed exclusively in the lysosomal compartment of monocytes and macrophages[1]. The promoter region of the NRAMP1 gene possesses a polymorphism within a possible enhancer element containing a Z-DNA-forming dinucleotide repeat[2,3].

It has been proposed that NRAMP1 polymorphisms play a role in susceptibility to mycobacterial infections[4,5]. To date, studies have yielded contradictory results[2]. In a West African population, NRAMP1 variants have been associated with susceptibility to Mycobacterium tuberculosis[6]. Malik et al[7] reported that variants of NRAMP1 are associated with tuberculosis (TB) in children. NRAMP1 gene polymorphisms have been linked with genetic susceptibility to infection with M. tuberculosis and progression of TB into severe clinical forms in eastern China[8]. Abe et al[9] reported that genetic variation in the NRAMP1 gene is associated with tuberculous cavitation of the lungs in Japanese patients. Kim et al[10] found that NRAMP1 polymorphisms are associated with tuberculous pleurisy.

Various studies have been carried out to establish a connection between NRAMP1 variants and diseases caused by mycobacteria other than the TB bacilli, such as Mycobacterium ulcerans causing Buruli ulcers[11] and Mycobacterium leprae[12]. A recent study[5] showed that the NRAMP1 promoter region polymorphism are positively associated with leprosy but not with the Mitsuda reaction (intradermal injection of lepromin), and that variants of the NRAMP1 gene favour microbial survival inside the macrophages by blocking the efficient transport of iron. It was reported that the NRAMP1 gene may have a role in some autoimmune diseases, including rheumatoid arthritis (RA)[13]. It was also reported that NRAMP1 823 C/C prevents the development of rheumatoid nodules in RA patients[14], type 1 diabetes[15] and multiple sclerosis[16]. Variants of the NRAMP1 gene have been associated with improved response to Bacillus Calmette-Guerin immunotherapy for superficial bladder cancer[17]. Hofmeister et al showed that NRAMP1 variants are specifically associated with Crohn’s disease (CD)[18].

In a previous study, we found that CD is associated with polymorphisms of the NOD2 gene and the pathogen MAP in Sardinian patients[19]. NOD2 protein is an intracellular protein that activates NFkB upon binding to microbial peptidoglycan. MAP is the etiological agent of Johne’s disease, a granulomatous enteritis of ruminants and other monogastric animals[20]. Although the association between MAP and CD (a human equivalent of paratuberculosis) has been postulated for a long time, only recent studies have made improvements in isolation and genomic techniques have allowed this link to be firmly established in a number of different populations[21-24].

As a part of our ongoing studies with NOD2/CARD15, we also looked at the role of other genes, especially NRAMP1, which may be associated with the survival of intracellular pathogens such as MAP. The role of NRAMP1 in regulating microbial survival inside phagosomes is related to iron transport, although the mechanism has not yet been completely elucidated[5,25].

MATERIALS AND METHODS

Patients were recruited at the University of Sassari after giving written consent. Using a case-control design, we analyzed 37 CD patients and 34 controls with no inflammatory bowel disease (IBD).

We searched for polymorphisms in 6 loci previously reported to be the most significantly associated with IBD[13]. The SCL11A1 polymorphisms that we analyzed included a (GT)n microsatellite in the promoter region, (-)237 C/T; 469 + 14G/C, INT4; a non-conservative base substitution at codon 543 (D543N); 823 C/T; and a 4-bp TGTG deletion locating 55 nucleotides downstream of the last codon in exon 15 (1729 + del55del4).

Amplification was performed as previously reported[13] using 100 ng of template genomic DNA previously extracted from intestinal tissues. The primer sequences that we used were previously reported[13]. Detection of MAP also was performed as previously reported[24].

Logistic regression was used to calculate odds ratios (OR) for the associations among CD, presence of MAP, and the 6 loci described above. Saturated and reduced models were computed and their comparability was assessed using the likelihood ratio test. The reduced model was assessed for goodness of fit using the covariate patterns as groups. In addition, the actual and model-predicted probabilities of CD according to levels of the independent variables in the final model were calculated.

RESULTS

Table 1 shows the full and reduced logistic regression models computed for this study. The full model showed significant associations between CD and the presence of MAP, the 1729 + 55del4 deletion polymorphism and the 823 C/T CT polymorphism. No significant associations were found between CD and the INT4, D534R, GT (n) and (-) 237 C/T loci. These non-significant loci were dropped in the reduced model and likelihood ratio testing showed no difference between the full and reduced models (likelihood ratio test: P = 0.436). Furthermore, the reduced model showed good fit (Pearson goodness-of-fit test: P = 0.407) and predictive power. We tested for significant interactions between MAP infection and each of the genetic loci, none of which was statistically significant. We therefore selected the reduced main effect model as our final model. No association was found between MAP and the 823 C/T or 1729 + 55del4 polymorphisms.

Table 1 Full and reduced logistic regression main effects models for the association between Crohn's disease, MAP infection, and the presence of a number of genetic polymorphisms of the NRAMP1 gene.

		Full model		Reduced model2
		Odds Ratio	P	Odds Ratio	P
MAP	Negative	1.01		1.01
	Positive	45.5	< 0.001	42.4	< 0.001
1729 + 55del4	+/+	1.01		1.01
	+/DEL, DEL/DEL	5.6	0.014	5.6	0.011
823 C/T	CC	1.01		1.01
	CT	75.6	0.001	50.8	0.001
D534R	GG	1.01
	GC	1.1	0.913
(-)237C/T	CC	1.01
	CT, TT	0.5	0.350
INT4	GG	1.01
	GC, CC	2.5	0.224
GT (n)	Allele1	1.01
	Allele2	1.2	0.807
	Allele3	0.6	0.673

¹Reference category;

²Model statistics: Likelihood ratio test for equivalence with full model: P = 0.436; Pearson goodness-of-fit test: P = 0.407; Area under Receiver-Operating Characteristic curve: 0.886.

Figure 1 shows the actual and final model-predicted probabilities of CD by MAP infection and the two statistically significant loci, 823 C/T and 1729 + 55del4. The close similarity between the actual and predicted probabilities reflected the good fit of our final model to the data. Strong independent effects of MAP infection and each of the loci on the probability of CD were observed (Figure 1). Even with no mutant polymorphism, MAP infection was highly predictive of CD. Having the CT polymorphism at the 823 locus, the probability of CD was greatly increased both among patients infected with MAP and among patients not infected with MAP. The effect of deletions at the 1729 locus was more moderate. It should be noted that these probabilities were sensitive to the ratio of cases to controls in the study, which were shown here to illustrate the particularly strong association among CD, MAP infection and the 823 C/T locus as well as the close fit between the experimental data and our chosen statistical model.

Figure 1 Actual probability of CD and predicted probability of CD showing strong independent effects of MAP infection and each of the loci on the probability of CD.

DISCUSSION

For the first time, a strong association has been observed between CD and polymorphisms at the 823C/T and 1729 + 55del4 loci in the NRAMP1 gene in Sardinians.

Although the combination of MAP infection and 823 CT mutation could perfectly predict CD, generalizations were limited by the small sample size in our study. Only 8 patients in this study simultaneously had MAP infection, 823 CT mutation, and CD. It would be worthwhile to follow up this finding in a larger study to determine if these two factors are a useful prognostic index for the development of CD.

It was reported that polymorphisms in the NOD2 gene and the presence of MAP are strongly associated with CD[19]. Both NOD2 polymorphisms and MAP infection had a strong independent association with CD, and are also associated with one another, suggesting that susceptibility to MAP infection may be influenced by NOD2 polymorphisms[26]. NOD2 protein is involved in the activation of NFkB factor and failure of its priming signal causes failure of pathogen clearance, possibly explaining the abnormal adaptive immune response to pathogens[27].

Our findings agree and disagree with some previous genetic studies of various NRAMP1 loci and CD. Our finding which is lacking of association between GT (n) alleles and CD is consistent with the failure to find an association between any of these alleles and IBD in Americans[28]. Stokkers et al[29] found that mutations at the 823C/T or (-) 274C/T loci mutations are not associated with CD, which is in agreement with our findings at (-) 274C/T locus in this study. However, we have found a very strong association between CD and the 823CT polymorphism.

Although previous studies have suggested that NRAMP1 mutations may favour microbial survival, the fact that this study failed to find any association between NRAMP1 polymorphisms and MAP infection does not support that viewpoint. Interestingly, although a previous study[6] showed allele 2 of this polymorphism is associated with mycobacterial infections, we did not find such an association between this allele and MAP infection.

The lack of association between NRAMP1 polymorphisms and MAP infection is particularly curious because NRAMP1 is thought to directly impact microbial survival inside of macrophages. NRAMP1 is thought to either deny needed iron in intraphagosomal microbes, or to increase transphagosomal iron needed to create bactericidal hydroxyl radicals via the Haber-Weiss/Fenton reaction[25].

The 823CT and 1729 + 55del4 polymorphisms have been associated with autoimmune disease[14], but not with another mycobacterial disease, tuberculosis[30], which probably does not have an autoimmune etiology. In the present study, these polymorphisms were found to be associated with CD which probably does have an autoimmune component, but not associated with MAP infection. These findings, combined with the strong association between MAP infection and CD, suggest that NRAMP1 polymorphisms may not cause CD by affecting the survival of MAP in intestinal tissue, but rather work together with MAP infection to cause CD via an autoimmune mechanism. In other words, one possible interpretation may be that these data support the aetiology of CD with both infectious and autoimmune components. At least, one study has elucidated a mechanism by which NRAMP1 affects the survival of intracellular pathogens and is simultaneously involved in autoimmune responses[25].

A further step of the study is to test the derived model on a large cohort of patients and controls.

In conclusion, CD may be strongly influenced by both NRAMP1 and MAP, although NRAMP1 polymorphisms and MAP infection are not themselves correlated. Combined with previous work on the NOD2/CARD15 gene, a complex interplay of genetic, immunologic and infectious factors may play a role in the aetiology of CD.

ACKNOWLEDGMENTS

The authors are thankful to Professor Giovanni Fadda for his guidance and support. The authors are also thankful to the International Society for Genomic and Evolutionary Microbiology (ISOGEM) for supporting and endorsing the study.