Journal of Tropical MedicineVolume 2012 (2012), Article ID 184819, 6 pagesdoi:10.1155/2012/184819

Research Article

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¹

¹Department of Genetics, Federal University of Paraná, P.O. Box 19071, 81531-980 Curitiba, PR, Brazil²Core for Advanced Molecular Investigation, Medical School, Pontifical Catholic University of Paraná, Imaculada Conceição, 1155, 80215-901 Curitiba, PR, Brazil

Received 31 May 2011; Revised 16 October 2011; Accepted 13 December 2011

Academic Editor: Ib Christian Bygbjerg

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.

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 ().

Table 1: Frequency distributions of allele and genotype data for the six studied SNPs in patients (: multibacillary; : paucibacillary) and contact controls, as indicated by the nucleotide base.

SNPa
Sample (n)
Genotypes (%)
Alleles (%)

rs2863381

TT
TC
CC
T
C

Patients (161)

M (111)

P (50)

Controls (143)

rs4440084

GG
GC
CC
G

Patients (157)

M (107)

P (50)

Controls (146)

rs1126680

GG
GA
AA
G
A

Patients (163)

M (113)

P (50)

Controls (142)

rs1799807b (D70G)

AA
AG
GG
A
G

Patients (165)

M (114)

P (51)

Controls (150)

rs1803274 (A539T)

GG
GA
AA
G
A

Patients (163)

M (113)

P (50)

Controls (149)

rs4387996

GG
GA
AA
G
A

Patients (162)

M (112)

P (50)

Controls (148)

aMost frequent alleles in Caucasians from the HapMap: rs2863381 (T), rs4440084 (G), rs1126680 (G), rs1799807 (A), rs1803274 (G), and rs4387996 (G).

bSignificant statistical differences between patients and controls (in bold) for genotype frequencies (; χ ² = 6.05 after the Yates correction) and for allele distributions (; χ ² = 5.90 after Yates correction).

Figure 1: Values of (%) obtained from Haploview for pairs of the studied SNPs in the total sample of patients plus contact controls: rs4387996 downstream the BCHE gene (1); rs1126680 (2), rs1799807 (3), and rs1803274 (4) in the BCHE gene; rs4440084 (5) and rs2863381 (6) upstream the BCHE gene.

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).

Table 2: Results from a forward stepwise multiple logistic regression in which leprosy (0 = contact control, 1 = patient) was the dependent variable.

Independent  variablesa
B ± S.E.
Wald
df

Odds ratio (95% C.I.)

Median age

13.03
1
0.000
2.48 (1.51 to 4.05)

SNP rs4387996
−
5.34
1
0.021
0.64 (0.44 to 0.93)

SNP rs1799807

4.34
1
0.037
5.18 (1.10 to 24.33)

Constant
−
5.07
1
0.024
0.11

aMedian age (≤51 = 1, >51 = 2); SNP rs4387996 ; SNP rs1799807 . Other independent variables were not significant: rs2863381; rs4440084; rs1126680; rs1803274 and sex.

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.

Table 3: Haplotypea frequencies compared between contact control (C) and patient (P) samples.

SNPsb

Samples (%)
; (Bonferroni’s correction)

rs2863381
rs4440084
rs1126680
rs1799807
rs1803274
rs4387996
C
P

T
G
G
A
G
G
24.9
25.1
n.s.

T
G
G
A
G
A
12.0
13.8
n.s.

T
C
G
A
G
G
11.6
18.5
5.28; 0.0216 (0.151)

T
C
G
A
G
A
9.5
3.4
8.62; 0.0033 (0.023)

C
C
G
A
G
G
8.9
3.8
6.69; 0.0097 (0.068)

C
C
G
A
G
A
8.1
6.5
n.s.

C
G
G
A
G
G
7.1
6.5
n.s.

aSelected on the basis of frequency higher than 5% at least in one of the samples. bMost frequent nucleotides in Caucasians from the HapMap: rs2863381 , rs4440084 , rs1126680 , rs1799807 , rs1803274 , and rs4387996 .

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.

M. T. Mira, “Genetic host resistance and susceptibility to leprosy,” Microbes and Infection, vol. 8, no. 4, pp. 1124–1131, 2006.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

B. Woolf, “On estimating the relation between blood group and disease,” Annals of Human Genetics, vol. 19, pp. 251–253, 1955.

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.

K. Kawashima and T. Fujii, “Extraneuronal cholinergic system in lymphocytes,” Pharmacology and Therapeutics, vol. 86, no. 1, pp. 29–48, 2000.

K. J. Tracey, “The inflammatory reflex,” Nature, vol. 420, no. 6917, pp. 853–859, 2002.

F. Herz and E. Kaplan, “A review: human erythrocyte acetylcholinesterase,” Pediatric Research, vol. 7, no. 4, pp. 204–214, 1973.

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.

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.