Investigation of Association between Susceptibility to Leprosy and SNPs inside and near the <i>BCHE</i> Gene of Butyrylcholinesterase

Gomes, Henrique J. P.; Souza, Ricardo L. R.; Prevedello, Flávia Costa; Mira, Marcelo Távora; Chautard-Freire-Maia, Eleidi A.

doi:https://doi.org/10.1155/2012/184819

Journal of Tropical Medicine

On this page

Abstract Introduction Materials and Methods Results and Discussion References Copyright Related Articles

Research Article | Open Access

Volume 2012 | Article ID 184819 | https://doi.org/10.1155/2012/184819

Investigation of Association between Susceptibility to Leprosy and SNPs inside and near the BCHE Gene of Butyrylcholinesterase

Henrique J. P. Gomes,¹Ricardo L. R. Souza,¹Flávia Costa Prevedello,²Marcelo Távora Mira,²and Eleidi A. Chautard-Freire-Maia¹

Academic Editor: Ib Christian Bygbjerg

Received31 May 2011

Revised16 Oct 2011

Accepted13 Dec 2011

Published22 Feb 2012

Abstract

Leprosy is a chronic disease caused by Mycobacterium leprae and affects the skin and the peripheral nervous system. Butyrylcholinesterase is coded by the BCHE gene, and the atypical allele (70G; rs1799807) has been investigated as a leprosy risk factor, with conflicting results. The present study estimated the frequencies of variants of rs1799807 and of five additional SNPs at the BCHE gene or near it: rs1126680, rs1803274, rs2863381, rs4440084, and rs4387996. A total of 167 patients and 150 healthy controls were genotyped by TaqMan PCR. Significantly higher allelic (70G) and genotypic (70DG) frequencies in rs1799807 were found in the patient group, with odds ratio (OR) of 6.33 (1.40 to 28.53) for the heterozygote. This finding was replicated in a comparison of the cases against a control group of 361 blood donors. The present data suggest that the atypical BChE variant may predispose to leprosy per se.

1. Introduction

Leprosy is a granulomatous, chronic infectious disease that, in spite of its ancient origin, still affects thousands of people around the globe. It is caused by the Gram-negative bacteria Mycobacterium leprae and injures predominantly the skin and the peripheral nervous system, leading to sensory loss in the skin, muscle weakness, and, often, hand and foot permanent disabilities. The WHO classifies leprosy in two clinically distinct groups: paucibacillary leprosy, characterized by five or fewer skin lesions and negative bacilloscopy, and multibacillary leprosy, which includes patients with six or more skin lesions and positive bacilloscopy. Individuals with the paucibacillary form usually present a strong cell-mediated response (Th1 type), activating NK cells and macrophages that provide a more efficient response against the bacteria. Individuals with multibacillary leprosy show humoral response (Th2 type), which suppresses macrophages, increasing the chances of systemic infection [1, 2].

The genetic background of the host influences the type of damage caused by the immune responses against an etiologic agent. It is well known that leprosy is a complex disease with still unknown environmental and genetic risk factors. Many candidate genes concerned with susceptibility to infection per se, clinical manifestation, and reversal reaction have been studied to date, such as non-MHC genes VDR (12q13.11) [3]; NRAMP1 (2q35) [4, 5], TLR2 (4q32) [6], PARK2, and PACRG (6q25-q26) [7], and MHC genes DQB1, DQA1, DRB1 [8], MICA [9], TNFA [10], and C4B [11]. In the present study, the candidate is the BCHE gene that codes for butyrylcholinesterase.

Human butyrylcholinesterase (BChE; EC 3.1.1.8; OMIM 177400) is a serum enzyme produced by the liver that hydrolyses esters of choline and other esters. BChE has been associated with lipid metabolism and xenobiotics detoxification [12, 13]. More than 70 variants of the BCHE gene, located at chromosomal region 3q26.1-q26.2, have already been described [14]. The atypical variant 70G (209 , rs1799807) codes for a BChE resistant to the hydrolysis of the muscle relaxant succinylcholine and may cause prolonged apnea after its administration [15]. Thomas et al. [16], based on the observation that some patients with leprosy presented prolonged apnea after the administration of succinylcholine, searched for a relationship between leprosy and BChE and showed association of the atypical BChE variant (70G) with this disease. However, data on this association are controversial: significantly higher frequencies of 70G in patients when compared to controls were found again [17], whereas other studies failed to replicate this result [18, 19], possibly due to the lack of accuracy of the enzyme inhibition method of phenotyping. Here, we present results of an association study between leprosy and BCHE variants, including 70G, characterized by DNA genotyping.

2. Materials and Methods

2.1. Samples

A case-control approach was applied, including 167 nonconsanguineous patients and two independent control samples: the first was composed of 150 healthy individuals (contact controls) reporting previous history of exposure to leprosy through current or past cohabitation with leprosy affected individuals and no consanguinity with patients; the second control group was composed of blood donors recruited from the same geographic area and ethnic background as the cases. Patients were classified according to the WHO criteria [1]: 116 with multibacillary and 51 with paucibacillary leprosy. Patients and contact controls are from the Brazilian state of Paraná, being predominantly from Caucasian origin (%; %, resp.), representing the ethnic distribution of this state. The control group of male Euro-Brazilian blood donors (age range between 18 to 30 years) was described elsewhere [20] and is the only Brazilian population sample with BCHE SNPs data determined by DNA analysis.

2.2. Laboratory Methods

DNA was extracted from whole blood by salting out [21]. Tag SNPs were selected from the HapMap Project website (http://www.hapmap.org), each representing a different bin of linkage disequilibrium and having a minimum allele frequency (MAF) of 5% in the Caucasian samples. Following this criterion, two SNPs upstream (rs2863381 and rs4440084), three intragenic (rs1126680, G–116A; rs1799807, D70G; rs1803274, A539T), and one downstream the BCHE gene (rs4387996) were selected for genotyping.

Genotyping was performed by fluorescence-based allelic discrimination as implemented in the ABI 7500 TaqMan platform. A 10 μL solution containing 1.5 μL of DNA solution (20 ng/μL), 0.3 μL of the TaqMan genotyping kit, 5.0 μL of TaqMan universal PCR Master Mix, and 3.2 μL of ultrapure water was submitted to PCR (45 cycles at 50°C for 2 min, 95°C for 10 min, 95°C for 15 secs of denaturation, and 60°C for 1 min of hybridization and extension).

SNP rs1803274 (A539T) was also genotyped by PCR-SSCA for quality control. This SNP was amplified by PCR (35 cycles at 94°C, 48°C, and 72°C for 1 min each; final extension at 72°C for 10 min): 1 μL of DNA (about 100 ng), 9 μL of PCR super MIX (Invitrogen), and 10 pmol of each primer (P45 and P43 [22]). A volume of 10 μL of PCR product was added to 10 μL of the denaturation solution (95% formamide, 0.25% bromophenol blue, 0.25% xylene cyanole, 10 mM EDTA, and 10 mM NaOH), submitted to 94°C (5 min), and kept chilled on ice until application in a polyacrylamide gel (8%), followed by silver nitrate staining [23].

2.3. Statistical Analysis

Frequency distributions, means ± standard errors, -tests, and tests were calculated using Statistica for Windows (Statsoft, Inc., 2000; http://www.statsoft.com). Arlequin 3.0 [24] was used to test for the Hardy-Weinberg equilibrium (HWE) and for estimating haplotype frequencies by maximum likelihood, using an expectation-maximization (EM) algorithm for genotypic data with unknown gametic phase. Linkage disequilibrium (LD) was estimated by the parameter as implemented in Haploview 4.2 [25]. Odds ratios were also calculated [26]. SPSS 13.0 for Windows (IBM; SPSS Inc., 1989–2004; http://www.spss.com) was used for forward (Wald) stepwise multiple logistic regression analyses. Bonferroni’s correction for multiple testing was used to adjust values, when necessary.

3. Results and Discussion

All markers were in HWE in patient and contact control samples. Figure 1 shows pairwise estimates of LD for all six studied SNPs in the samples of patients and contact controls; highest value was 26%, confirming independence of all tag SNPs genotyped in this study. Table 1 shows frequency distributions of allele and genotype data for all tag SNPs. Significant statistical differences were found between patients and contact controls for genotypic ( = 6.05; ) and allelic ( = 5.90; ) frequency distributions of SNP rs1799807: % of heterozygotes (70DG) and % of the 70G allele in patients; % of heterozygotes and % of the 70G allele in contact controls. The odds ratio (OR) was 6.33 (CI_(95%) from 1.40 to 28.53) for the 70DG heterozygote. Pauci and multibacillary patients did not differ on 70G or 70DG frequencies, suggesting an impact of this variant over the early events of interaction between host and pathogen that result in infection per se, that is, the disease independent of its clinical form. Based on previous reports [16, 17] of association of the atypical variant with higher risk of susceptibility to leprosy, our initial hypothesis was that the frequency of this variant would be higher in patients than in controls. In view of this, Bonferroni’s correction was not used for this case. Statistically significant difference between patients and contact controls was also found for allele frequencies of the SNP rs4387996 (Table 1) located downstream the BCHE gene ( = 5.07; ). However, this significance did not resist Bonferroni’s correction ().

Patients and contact controls differ in sex composition ( = 5.20; ) and mean age ( = 5.10; = 5.72 ×10⁻⁵): the M/F ratio and mean age of the patients were 1.3 and years and of the controls were 0.8 and years, respectively. Since both gender and age are well-known risk factors for leprosy, these covariables were included, as independent variables, in a forward stepwise multiple logistic regression analysis in which leprosy was the dependent variable (Table 2). The strategy revealed three significant independent variables: median age, SNP rs4387996, and SNP rs1799807. In the presence of the age effect, association analysis revealed that the allele of rs4387996 is protective, whereas the risk of the 70DG genotype for acquiring leprosy is 5.18 times higher than that of the 70DD genotype, confirming our previous analyses. The rs4387996 association was significant in both additive (Table 2, ) and dominant models for allele A (GG versus GA+AA; ; data not shown).

To confirm the findings for the 70DG genotype, data from the case group was compared with genetic information obtained previously for a sample of 361 Euro-Brazilian blood donors of Curitiba [20], which presented 3.60% of heterozygotes and 1.80% of 70G frequency. Again, statistically significant differences for genotype ( = 4.45; ) and allele ( = 4.36; ) frequencies were found. Estimate of OR for this second comparison was 2.29 (CI_(95%) from 1.04 to 5.05). The blood donor sample is appropriated for comparison, given the similar Euro-Brazilian origin of the cases, and the result clearly reinforces our hypothesis. Of note, the OR found in the comparison of cases against blood donor controls is lower than the one found for cases against contact controls. This may be due to the selection criterion for the contact control that may have led to a sample with a higher frequency of individuals presenting innate resistance to leprosy, as they have more likely been exposed than the blood donor controls and yet did not develop the disease. Furthermore, the contact controls report no consanguinity with any leprosy patient. Unfortunately, DNA data of the blood donor control group were not available for marker rs4387996, also associated with leprosy in the primary comparison.

To further advance on the understanding of the impact of BCHE variants over leprosy susceptibility, we performed a haplotypic analysis, involving all haplotypes with frequency > 5% in at least one of the samples (Table 3). Significantly higher frequency of the [T; C; G; A; G; A] haplotype was found in contact controls when compared to cases ( after Bonferroni’s correction, OR = 0.34; CI_(95%) from 0.16 to 0.72). Although genotypes were of unknown phase and haplotype frequencies had to be estimated [24], the frequencies of this haplotype are expected to be rather accurate, considering that the EM algorithm usually presents very reliable performances for haplotype estimation [27]. Furthermore, both the total number of gametes (contact controls: 268; patients: 297) and the estimated haplotype numbers for contact controls () and patients () reinforce the accuracy of this frequency inference.

This protective haplotype differs from the most frequent one observed among cases ([T; G; G; A; G; G], 25.07%) only at positions corresponding to SNPs rs4440084 and rs4387996. These two variants (C; A) are found in the cis conformation in 10.92% and in 18.86% of haplotypes in patients and contact controls, respectively ( = 8.93; ). The OR value for this comparison is 0.49 (CI_(95%) 0.30–0.79) and may indicate a protective effect that is two times higher than that of the haplotypes-bearing variants and in these respective SNPs. Considering that the rs4440084 SNP was not significant in the regression analysis (Table 2), it is possible that the protective effect of this haplotype is only due to the variant of SNP rs4387996. However, the SNP independent variables in the regression analyses were the genotypes, as haplotypes were not known individually.

If the present data referring to the atypical variant (70G) is not due to linkage disequilibrium with other causal genetic variants, it is possible that the atypical variant of BChE may predispose to leprosy, although the physiological mechanism remains unclear. One possible hypothesis would be a relation of 70G with the acetylcholine- (ACh-) mediated immune response. Activated leucocytes release small quantities of ACh [28], which acts like autocrine and paracrine signaling factors, mediating the contact with target cells. Moreover, when ACh is produced, it acts as part of an anti-inflammatory parasympathetic reaction, inhibiting the nuclear action of NF-κB and stopping the production of proinflammatory citokynes as TNF, IL-1, IL-6, IL-8, and HMGB1, leading to the interruption of the inflammatory process [29]. Therefore, high levels of ACh are likely to cause a decrease in the inflammatory response. It is possible to speculate that this decrease in the inflammatory process can lead to a less competent host response to M. leprae. Acetylcholinesterase (AChE) activity is inhibited by excess of ACh [30]. Once AChE is inhibited, BChE becomes responsible for ACh hydrolysis [31]. The usual (70D) BChE variant is able to accomplish this function, but the atypical variant (70G), which has less affinity to the ACh substrate, could interfere in this regulation, failing in reducing ACh levels, consequently preventing the inflammatory response and contributing to the acquisition of leprosy.

Considering that SNP rs4387996 independently associated with leprosy in our population sample is intergenic, it is difficult to infer a functional role. However, one possible hypothesis is that the variants (G, A) are in linkage disequilibrium with variants of the SLITRK3 gene located downstream the BCHE gene but in the same LD bin. The Slitrk3 protein encoded by this gene is a neuronal transmembrane protein that controls axons and dendrites outgrowth [32]. The association among BCHE variants and the SLITRK3 gene is unclear, but it is possible that they may be markers of the influence of this gene in the neuronal development which, somehow, can protect the nerves from M. leprae.

The present study reveals two independent signals of association between leprosy and variants of the BCHE gene, including the known functional variant 70G (rs1799807), establishing BCHE as a strong candidate gene for the control of leprosy susceptibility. It is not clear whether the association of the 70G variant with leprosy is really functional or is due to linkage disequilibrium with a variant of another, unknown gene. However, it is valid to suppose that the functional impact of this SNP that can lead to a less efficient BChE may predispose to the disease by failing to degrade a supposed substrate, such as acetylcholine, whenever it should act upon it. This fact could disturb the immune response when mediated by acetylcholine, reducing the inflammatory process and increasing the susceptibility to leprosy.

References

WHO, “Chemotherapy of leprosy for control programmes,” Technical Report Series 675, Geneva, Switzerland, 1982.
View at: Google Scholar
M. T. Mira, “Genetic host resistance and susceptibility to leprosy,” Microbes and Infection, vol. 8, no. 4, pp. 1124–1131, 2006.
View at: Publisher Site | Google Scholar
S. Roy, A. Frodsham, B. Saha, S. K. Hazra, C. G. N. Mascie-Taylor, and A. V. S. Hill, “Association of vitamin D receptor genotype with leprosy type,” Journal of Infectious Diseases, vol. 179, no. 1, pp. 187–191, 1999.
View at: Publisher Site | Google Scholar
L. Abel, F. O. Sánchez, J. Oberti et al., “Susceptibility to leprosy is linked to the human NRAMP1 gene,” Journal of Infectious Diseases, vol. 177, no. 1, pp. 133–145, 1998.
View at: Google Scholar
A. Alcaïs, F. O. Sanchez, N. V. Thuc et al., “Granulomatous reaction to intradermal injection of lepromin (Mitsuda reaction) is linked to the human NRAMP1 gene in Vietnamese leprosy sibships,” Journal of Infectious Diseases, vol. 181, no. 1, pp. 302–308, 2000.
View at: Publisher Site | Google Scholar
T. J. Kang and G. T. Chae, “Detection of Toll-like receptor 2 (TLR2) mutation in the lepromatous leprosy patients,” FEMS Immunology and Medical Microbiology, vol. 31, no. 1, pp. 53–58, 2001.
View at: Publisher Site | Google Scholar
M. T. Mira, A. Alcaïs, H. Van Thuc et al., “Susceptibility to leprosy is associated with PARK2 and PACRG,” Nature, vol. 427, no. 6975, pp. 636–640, 2004.
View at: Publisher Site | Google Scholar
M. A. Shaw, I. J. Donaldson, A. Collins et al., “Association and linkage of leprosy phenotypes with HLA class II and tumour necrosis factor genes,” Genes and Immunity, vol. 2, no. 4, pp. 196–204, 2001.
View at: Publisher Site | Google Scholar
L. M. Wang, A. Kimura, M. Satoh, and S. Mineshita, “HLA linked with leprosy in southern China; HLA-linked resistance alleles to leprosy,” International Journal of Leprosy and Other Mycobacterial Diseases, vol. 67, no. 4, pp. 403–408, 1999.
View at: Google Scholar
S. Roy, W. McGuire, C. G. N. Mascie-Taylor et al., “Tumor necrosis factor promoter polymorphism and susceptibility to lepromatous leprosy,” Journal of Infectious Diseases, vol. 176, no. 2, pp. 530–532, 1997.
View at: Google Scholar
I. J. T. De Mesias, J. Santamaria, M. Brenden, A. Reis, and G. Mauff, “Association of C4B deficiency (C4B*Q0) with erythema nodosum in leprosy,” Clinical and Experimental Immunology, vol. 92, no. 2, pp. 284–287, 1993.
View at: Google Scholar
B. Li, E. G. Duysen, and O. Lockridge, “The butyrylcholinesterase knockout mouse is obese on a high-fat diet,” Chemico-Biological Interactions, vol. 175, no. 1–3, pp. 88–91, 2008.
View at: Publisher Site | Google Scholar
E. G. Duysen, B. Li, S. Darvesh, and O. Lockridge, “Sensitivity of butyrylcholinesterase knockout mice to (- -)-huperzine A and donepezil suggests humans with butyrylcholinesterase deficiency may not tolerate these Alzheimer's disease drugs and indicates butyrylcholinesterase function in neurotransmission,” Toxicology, vol. 233, no. 1–3, pp. 60–69, 2007.
View at: Publisher Site | Google Scholar
R. L. R. Souza, L. R. Mikami, R. O. B. Maegawa, and E. A. Chautard-Freire-Maia, “Four new mutations in the BCHE gene of human butyrylcholinesterase in a Brazilian blood donor sample,” Molecular Genetics and Metabolism, vol. 84, no. 4, pp. 349–353, 2005.
View at: Publisher Site | Google Scholar
M. C. McGuire, C. P. Nogueira, C. F. Bartels et al., “Identification of the structural mutation responsible for the dibucaine-resistant (atypical) variant form of human serum cholinesterase,” Proceedings of the National Academy of Sciences of the United States of America, vol. 86, no. 3, pp. 953–957, 1989.
View at: Google Scholar
M. Thomas, C. K. Job, and P. V. Kurian, “Susceptibility to leprosy and serum atypical pseudocholinesterase,” International Journal of Leprosy, vol. 44, no. 3, pp. 315–318, 1976.
View at: Google Scholar
M. Thomas and C. K. Job, “Serum atypical pseudocholinesterase and genetic factors in leprosy,” British Medical Journal, vol. 3, no. 823, pp. 390–391, 1972.
View at: Google Scholar
D. P. Agarwal, H. W. Goedde, and W. Schloot, “A note on atypical serum cholinesterase and genetic factors in leprosy,” Human Heredity, vol. 23, no. 4, pp. 370–373, 1973.
View at: Google Scholar
M. Whittaker, R. F. Lowe, and B. P. B. Ellis, “Serum cholinesterase variants in African leprosy patients resident in Rhodesia,” Human Heredity, vol. 26, no. 5, pp. 372–379, 1976.
View at: Google Scholar
L. R. Mikami, S. Wieseler, R. L. R. Souza et al., “Five new naturally occurring mutations of the BCHE gene and frequencies of 12 butyrylcholinesterase alleles in a Brazilian population,” Pharmacogenetics and Genomics, vol. 18, no. 3, pp. 213–218, 2008.
View at: Publisher Site | Google Scholar
S. A. Miller, D. D. Dykes, and H. F. Polesky, “A simple salting out procedure for extracting DNA from human nucleated cells,” Nucleic Acids Research, vol. 16, no. 3, p. 1215, 1988.
View at: Publisher Site | Google Scholar
T. Hohler, M. Hundt, C. Rittner, P. M. Schneider, and K. H. Meyer zum Buschenfelde, “Cholinesterase variants: rapid characterisation by PCR/SSCP and evidence for molecular homogeneity,” Journal of Medical Genetics, vol. 32, no. 2, pp. 109–112, 1995.
View at: Google Scholar
B. Budowle, R. Chakraborty, A. M. Giusti, A. J. Eisenberg, and R. C. Allen, “Analysis of the VNTR locus DIS80 by the PCR followed by right-resolution PAGE,” American Journal of Human Genetics, vol. 48, pp. 137–144, 1991.
View at: Google Scholar
L. Excoffier, G. Laval, and S. Schneider, “Arlequin ver. 3.0: an integrated software package for population genetics data analysis,” Evolutionary Bioinformatics Online, vol. 1, pp. 47–50, 2005.
View at: Google Scholar
J. C. Barrett, B. Fry, J. Maller, and M. J. Daly, “Haploview: analysis and visualization of LD and haplotype maps,” Bioinformatics, vol. 21, no. 2, pp. 263–265, 2005.
View at: Publisher Site | Google Scholar
B. Woolf, “On estimating the relation between blood group and disease,” Annals of Human Genetics, vol. 19, pp. 251–253, 1955.
View at: Google Scholar
E. C. Castelli, C. T. Mendes-Junior, L. C. Veiga-Castelli, N. F. Pereira, M. L. Petzl-Erler, and E. A. Donadi, “Evaluation of computational methods for the reconstruction of HLA haplotypes,” Tissue Antigens, vol. 76, no. 6, pp. 459–466, 2010.
View at: Publisher Site | Google Scholar
K. Kawashima and T. Fujii, “Extraneuronal cholinergic system in lymphocytes,” Pharmacology and Therapeutics, vol. 86, no. 1, pp. 29–48, 2000.
View at: Publisher Site | Google Scholar
K. J. Tracey, “The inflammatory reflex,” Nature, vol. 420, no. 6917, pp. 853–859, 2002.
View at: Publisher Site | Google Scholar
F. Herz and E. Kaplan, “A review: human erythrocyte acetylcholinesterase,” Pediatric Research, vol. 7, no. 4, pp. 204–214, 1973.
View at: Google Scholar
B. Li, J. A. Stribley, A. Ticu et al., “Abundant tissue butyrylcholinesterase and its possible function in the acetylcholinesterase knockout mouse,” Journal of Neurochemistry, vol. 75, no. 3, pp. 1320–1331, 2000.
View at: Google Scholar
J. Aruga and K. Mikoshiba, “Identification and characterization of Slitrk, a novel neuronal transmembrane protein family controlling neurite outgrowth,” Molecular and Cellular Neuroscience, vol. 24, no. 1, pp. 117–129, 2003.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2012 Henrique J. P. Gomes et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1333

Downloads

1060

Citations