http://orcid.org/0000-0002-0449-818X
Figures
Abstract
Background
Human African Trypanosomiasis (HAT) is a neglected disease targeted for elimination as a public health problem by 2020. Elimination requires a better understanding of the epidemiology and clinical evolution of HAT. In addition to the classical clinical evolution of HAT, asymptomatic carriers and spontaneous cure have been reported in West Africa. A genetic component to human susceptibility to HAT has been suggested to explain these newly observed responses to infection. In order to test for genetic associations with infection response, genetic polymorphism in 17 genes were tested (APOL1, IL1B, IL4, IL4R, IL6, IL8, IL12B, IL12RB1, IL10, TNFA, INFG, MIF, HLA-G, HLA-A, HP, HPR and CFH).
Methodology
A case-control study was performed on 180 blood samples collected from 56 cases and 124 controls from Cameroon. DNA was extracted from blood samples. After quality control, 25 samples (24 controls and 1 case) were eliminated. The genotyping undertaken on 155 individuals including 55 cases and 100 controls were investigated at 96 loci (88 SNPs and 8 indels) located on 17 genes. Associations between these loci and HAT were estimated via a case-control association test.
Results
Analyses of 64 SNPs and 4 indels out of 96 identified in the selected genes reveal that the minor allele (T) of rs8062041 in haptoglobin (HP) appeared to be protective against HAT (p = 0.0002395, OR 0.359 (CI95 [0.204–0.6319])); indicating higher frequency in cases compared to controls. This minor allele with adjusted p value of 0.0163 is associated with a lower risk (protective effect) of developing sleeping sickness.
Author summary
Human African trypanosomiasis (HAT) or sleeping sickness is a neglected tropical disease targeted for elimination by 2020. This elimination requires a better understanding of the epidemiology and clinical evolution of this disease. Beside the classical clinical evolution, asymptomatic carriers, seropositive and spontaneous cure of infected persons have been reported in West Africa. Arguments in favor of human genetic susceptibility to HAT have been raised to explain this variability in clinical presentation. This study investigated the genetic polymorphism of 17 genes between controls and sleeping sickness patients in Southern Cameroon in order to improve our knowledge of human susceptibility to trypanosome infections. We identified single nucleotide polymorphisms and indels in 17 selected genes involved in immune responses and carried out a case-control candidate gene association study and demonstrated differences between variants associated with the disease. From these genes, only haptoglobin (HP) at the SNP rs8062041 was found to have polymorphisms which were strongly associated with trypanosomiasis. The minor allele (T) at this SNP position appeared to be protective against HAT (p = 0.0002395, OR 0.359 (CI95 [0.204–0.6319])) reducing the risk of developing disease approximately threefold. The haptoglobin related protein (HPR) is adjacent to HP and is a component of the Trypanolytic factor that kills trypanosomes. The HP and HPR locus is duplicated in some people. The rs8062041 variant may be associated with this duplication and it is possible that increased production of HPR is the cause of the protection associated with rs8062041. The results reported here will contribute to the knowledge of the role of human genetics in disease progression, and thus lead to the identification of novel biomarkers which could involve development of new diagnostics, treatments and intervention strategies.
Citation: Ofon E, Noyes H, Mulindwa J, Ilboudo H, Simuunza M, Ebo’o V, et al. (2017) A polymorphism in the haptoglobin, haptoglobin related protein locus is associated with risk of human sleeping sickness within Cameroonian populations. PLoS Negl Trop Dis 11(10): e0005979. https://doi.org/10.1371/journal.pntd.0005979
Editor: Jayne Raper, Hunter College, CUNY, UNITED STATES
Received: March 30, 2017; Accepted: September 20, 2017; Published: October 27, 2017
Copyright: © 2017 Ofon et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: We are grateful to the University of Dschang as well as the Wellcome Trust for funding this work (099310/Z/12/Z) as part of its support for the H3Africa initiative (http://h3africa.org/). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Human African Trypanosomiasis (HAT) or sleeping sickness is a lethal neglected tropical disease responsible for severe morbidity and economic losses in areas where it occurs [1]. HAT is caused by subspecies of Trypanosoma that are transmitted to humans through the bites of hematophagous flies of the genus Glossina, commonly known as tsetse flies. HAT exists in two forms: the acute form due to Trypanosoma brucei rhodesiense, which occurs in East Africa, and the chronic form due to T. b. gambiense, which is found in West and Central Africa. More than 98% of the currently reported cases belong to the chronic form. About 65 million people are estimated to be at risk of HAT and the current number of HAT cases is below 20,000 pa [2]. Control efforts undertaken during the last decades have reduced considerably the number of cases and 3,796 new cases were reported to WHO in 2014 [2]. With the success of these control efforts, HAT has been included in the WHO roadmap of neglected tropical diseases which are targeted for elimination as a public health problem by 2020. For effective control, it is important to gain a better understanding of the clinical evolution of the disease. Previously, HAT was classically considered to be fatal if untreated. During the last decades, a range of clinical presentations of T. b. gambiense HAT have been reported including asymptomatic carriers and spontaneous cure without treatment [3]. One hypothesis for the diversity of clinical outcomes that occur during infections due to T. b. gambiense is that it is due to human genetic variability. Previous investigations on genes such as HLA-A, HP, CFH, IL1B, IL12B, IL12RB1, IL4R and HPR [4, 5, 6, 7, 8, 9, 10, 11, 12] revealed associations between the polymorphisms in some of these genes with infectious diseases including HIV, viral hepatitis, malaria and tuberculosis. In HAT, polymorphisms in sequence or expression of genes involved in immune response such as APOL1, IL4, IL6, IL10, IL8, TNFA, HLA-G, MIF, HPR and INFG have been investigated for their association with the outcome of T. b. gambiense infections [13, 14, 15, 16, 17, 18, 19, 20, 21, 1]. These investigations found associations between some polymorphisms in genes and the risk of developing HAT. For instance, Courtin et al. [14–16] have shown that polymorphisms in IL6, IL10 and HLA-G were associated with a protective effect against HAT. In addition, a protective effect has been observed in vitro and in vivo of the G2 allele of APOL1 against infections due to T. b. rhodesiense [22]. In Guinea, the G1 allele of APOL1 was found to be associated with protection of asymptomatic individuals against development of active disease [22]. Despite these observations, the relationship between genetic polymorphisms and susceptibility to HAT is still not well understood and further investigations on populations of more HAT foci are required.
To improve our knowledge of the genetic determinants that could play important roles during infections due to T. b. gambiense, seventeen genes were selected in this study and their polymorphisms were investigated within populations of three sleeping sickness foci of the forest region of Cameroon.
Study area and study population
The study was conducted in three active sleeping sickness foci of the forest region of Southern Cameroon (Fig 1). The Cameroonian population is made up of more than 240 ethnic groups that can be grouped into Bantu (e.g.: Beti, Bassa, Bakundu, Maka, Douala, Pygmie), Semi Bantu (e.g.: Bamileke, Gbaya, Bamoun, Tika) and Sudano-Sao (e.g.: Fulbe, Mafa, Toupouri, Shoa-Arabs, Moundang, Massa, Mousgoum). The composition varies considerably between HAT foci and even within the same HAT focus. The three HAT foci where this study was undertaken were Bipindi and Campo in the Southern region and Fontem in the South-west region of Cameroon.
The Campo focus (2°82'00"N, 9°85'20"E) is located in the tropical forest and extends from the Atlantic coast along the Ntem river which delimits the Cameroon–Equatorial Guinea border. It is a hypo-endemic focus with no history of epidemic outbreaks [23] and a cumulative number of 98 cases were detected between 1998 and 2013. The main source of livelihood for the inhabitants of the Campo focus is agriculture, fishing and hunting. It is a cosmopolitan area with several ethnic groups including mainly the Iyassa, Kwasse, Maabi, Mvae and Ngoumba, most of whom are Bantu speaking. Other minor ethnic groups are semi Bantus and Sao-Sudanese and can be found at Campo for administrative and socioeconomic purposes.
The Bipindi HAT focus (3°82'00"N, 10°82'20"E) is located at about 75 km from the Atlantic coast in the South of Cameroon. It is an old HAT focus that has been known since 1920. During the last two decades, the Bipindi focus was among the most active HAT foci of Cameroon with around 83 HAT cases diagnosed from 1999 to 2011 [24]. About 95% of the inhabitants of the Bipindi HAT focus are Bantu speaking and belong to ethnic groups such as Ngoumba, Nti, Fan and Pygmies. The remaining 5% of inhabitants (semi Bantus and Sao-Sudanese) are there for administrative and socioeconomic purposes. The main livelihood for people in this focus is hunting, farming and seasonal harvesting of fruits.
The Fontem focus (5°40’00”N, 9°55’00”E) is located in the South-West Region of Cameroon where HAT has been known to occur since 1949 [25]. The Fontem focus was previously among the most active HAT foci of Cameroon [26], but in recent decades, it has become hypoendemic with about 8 patients detected among 16,000 persons examined between 1998 and 2007 [27]. In this focus, the Mundani, Bamoua and Bangwa are the major ethnics groups. Other minor ethnic groups such as Banyangue and Bamileke are also found.
Materials and methods
Sample collection
The blood samples were collected during medical surveys performed jointly with the national sleeping sickness control program of Cameroon. The sampling was done at Campo in 2014 and for Bipindi and Fontem, in 2015. During these surveys, all inhabitants were screened with CATT test [28] on whole blood. All inhabitants with positive CATT test were subjected to CATT dilution on plasma and each inhabitant positive on a CATT dilution ≥1/8 was subjected to parasitological examination (capillary tube centrifugation (CTC) [29] and minianion exchange centrifugation technique (mAECT)) [30]. For all inhabitants with CATT dilution ≥1/8 and negative for all parasitological tests, 90μl of plasma were spotted on a Whatman paper disc (divided in four equal parts with each bearing a spot of 30μl) that was sent to CIRDES in Burkina Faso for the trypanolysis test [31]. Beside the CTC and mAECT, lymph node aspiration followed by a microscopic examination was performed to search for trypanosomes in all individuals showing enlarged lymph nodes. A new HAT case was defined as an inhabitant in whom trypanosomes were seen by at least one parasitological method. Beside these new HAT cases, old HAT cases were also resampled. Old HAT cases were residents in whom trypanosomes had been previously seen on at least one parasitological test after passive or active case detection. Old HAT cases were only included in this study if the information regarding the clinical status, the CATT test and all parasitological tests were available in hospital records. A control was considered as any individual negative for the CATT test, all parasitological tests including CTC, mAECT and lymph node examination and when possible the trypanolysis test.
With these sampling criteria, 5ml of peripheral venous blood samples were collected from cases and controls into EDTA coated tubes. In the field, the tubes were mixed gently and stored at 4°C in an electric cooler before being transported to the laboratory.
Ethics statement
The protocol of this study was approved by the Ethical Committee of the Ministry of Public Health of Cameroon reference number N°2013/11/364/L/CNERSH/SP of 21 November 2013. The local administrative and traditional authorities of each HAT focus were also informed and gave their approval. Subsequently, the review board (LAMAS) of Laboratory of Microbiology and Anti-microbial substances of the Department of Biochemistry of the Faculty of Science of the University of Dschang gave their approval. All adult subjects provided informed consent, and a parent or guardian of any child participant below 18 years old provided informed consent on their behalf. Each informed consent was written because all individuals enrolled in this study gave their approval by signing an informed consent form and a Certificate of Confidentiality. During analyzes, data of each subject were anonymized.
DNA extraction
Blood samples were centrifuged at 5000rpm for 3 minutes and the buffy coat was collected. Genomic DNA was extracted from the Buffy-coat with the QIAamp DNA Blood Midi/Maxi kit (Qiagen) according to the manufacturer's instructions. The DNA was eluted with 200μl of elution buffer and stored at -20°C until use.
Power calculation
Power calculations were undertaken using the genetics analysis package gap in r [32]
Selection of candidate genes
The choices of candidate genes were based on previous observations. The cytokines IL4, IL6, IL10, IL8, INFG, TNFA, HP, HPR and MHC gene HLA-G were selected because they have been previously associated with HAT [14–16, 20, 33, 34, 35, 36, 37]. In addition, two genes for factors involved in the lysis of trypanosomes, APOL1 and haptoglobin-related protein (HPR) were also included [17, 18, 19]. Five further genes that had previously been reported to play an important role in the susceptibility to other infectious diseases were selected: Human Leukocytes Antigen A (HLA-A) [38, 4, 39], IL1B [10], Complement factor H (CFH) [6, 10], IL12B and IL12RB1 [5, 11] and Macrophage migration inhibitory factor (MIF) [40, 41, 42] genes were also included.
SNPs and INDELs selection
Most SNPs and indels for testing were selected after a Linkage scan (r² = 0.5) and quality control with Plink version 1.9 [43] using whole genome sequencing data. These data were obtained from a merged dataset between the African populations data from the 1000 Genomes Project combined with low fold coverage (8-10x) whole genome shotgun data generated from 230 residents living in regions (DRC, Guinea Conakry, Ivory Coast and Uganda, European Genome Archive A accession number) where trypanosomiasis is endemic [44]. The 88 SNPs and 8 indels loci were selected by two strategies: 1) by linkage scan of SNPs and indels (r2 < 0.5) across the gene; 2) by selection of SNPs and indels with published associations with HAT. Linked SNPs were identified for IL6, IL4, IL8, IFNG and HLA-G genes. For APOL1, HPR, HP, HLA-A, IL1B, IL12B, IL12RB1, IL4R, CFH, IL10, MIF and TNFA genes, individual published SNPs and indels were identified and selected based on literature searches.
Genotyping
Samples which had low DNA concentration or did not satisfy the quality control criteria were excluded prior to genotyping. Genotyping was performed by two commercial service providers: 1) “Plateforme Genome Transcriptome” at INRA of Bordeaux in France; 2) LGC Genomics Hoddesden, UK with approximately 1μg of genomic DNA per sample.
At INRA, genotyping was carried out with a Multiplex design (two sets of 40 SNPs or indels) using Assay Design Suite v2.0 (Agena Biosciences). For each SNP and indel, the genotyping was done with the iPLEX Gold genotyping kit (Agena Biosciences) for the Mass-Array iPLEX genotyping assay according to the manufacturer’s instructions. Products were detected on a Mass-Array mass spectrophotometer and data were obtained in real time with Mass-Array RT software (Agena Biosciences). SNP clustering and validation was carried out with Typer 4.0 software (Agena Biosciences). A summary of the candidate genes, and SNPs and indels is shown in the supplementary data S1 Table. Some SNPs and indels that failed genotyping at INRA and some additional SNPs and indels were genotyped at LGC Genomics, Hoddesden, UK where SNPs and indels were genotyped using the PCR based KASP assay [45].
Analysis
This was a case-control study where no familial controls were collected during sampling. The raw genotypic data were converted to PLINK format and quality control (QC) procedures implemented using the PLINK v1.9 package [43]. The Spearman Chi-square test was used to compare frequencies of observed and expected genotypes under Hardy–Weinberg equilibrium (HWE) and LD using R/Rstudio version 3.3.2 (2016-10-31)—‘Sincere Pumpkin Patch’ and Plink [43]. After quality control and filtering, poorly performing SNP loci with missing genotypes (≥10) and samples with missing loci (≥4) were removed. In addition, all loci with a MAF below 1% or a HWE P value < 1 × 10−4 were removed. SNP in linkage with adjacent SNP (r² > 0.5) were also pruned. These filters are as described by Anderson et al. [46] to minimize the influence of genotype-calling artifacts in a candidate gene study. The association between individual SNPs and indels within genes and HAT were tested using the Fisher exact test with Plink v1.9 software. Results were adjusted for multiple testing by Bonferroni correction. To show significant association during multiple tests, a single marker (SNP) must show, after Bonferroni correction, an alpha value (obtained P value before correction/number of SNPs analyzed) below 0.000746 (0.05/68). The Bonferroni correction assumes that each of the statistical tests is independent; however, this is not always true due to the possibility of linkage disequilibrium among the SNPs. In instances where the assumption is not true, the correction is often too strict, leading potentially to false negatives. A less stringent correction for multiple testing was also employed. The Benjamini-Hochberg false discovery rate (FDR) estimates the proportion of significant results (P < 0.05) when the Bonferroni correction considers them as false positives [47, 48]. FST is a measure of differences between populations. The analysis of FST was run to check for significant allele frequency difference between the cases and controls while Principal Component Analysis (PCA) was used to check for population stratification that might confound the analysis using Plink [43].
Results
Study design and population
This study was one of six studies of populations of HAT endemic areas in Cameroon, Cote d’Ivoire, Guinea, DRC, Malawi and Uganda. The studies were designed to have 80% power to detect odds ratios (OR) >2 for loci with disease allele frequencies of 0.15–0.65 and 100 cases and 100 controls with the 96 loci genotyped.
Overall, 216 individuals were included in this study: 56 (25.93%) HAT patients and 160 (74.07%) controls. The 216 individuals belonged to 22 different ethnic groups. The mean age (range) of HAT cases was 44.94 (15–82) years, while that of controls was 37.08 (9–86). The overall sex ratio (male/female) was 1.02 (109/107), with HAT cases being 0.75 (24/32) and controls 1.12 (84/75). Given that only 56 cases were available from our study area, the power of this study was reduced and it had 80% power to detect an OR >3 with disease allele frequencies of 0.1–0.45 with the 96 loci genotyped.
One hundred and eighty (56 HAT cases and 124 controls) of the 216 samples were sent for genotyping. After DNA quantification and quality control on each of these 180 samples, 25 were excluded from genotyping. 155 samples were genotyped: 55 (34.48%) HAT cases and 100 (64.52%) controls.
Genes and SNPs selected and genotyped
96 loci containing 88 SNPs and 8 indels were tested from 17 candidate genes. The number of SNPs and indels analyzed varied considerably (from 1 to 18) between genes (Table 1). The highest number of 18 SNPs and indels was observed for HLA-G and the lowest number of one SNP for IL10, IL1B and CFH. However, it is important to point out that for APOL1 (three SNPs), CFH, TNFA, HLA-A and IL10, the SNPs considered here are only those that have been already reported in the literature. Of the 88 SNPs and 8 indels used in this study, 24 SNPs and 4 indels (with 8 removed for MAF ≤1%, 7 for missing loci ≥10%, 5 with HWE P-values <1 x 10−4 and 8 for linkage at r² ≥ 0.5) of them were excluded during quality control which excluded one gene (HLA-A) completely. Four samples were also excluded during quality control due to missing individual data ≥4%. For subsequent analyses, 69.29% of loci including 64 SNPs and 4 indels from 16 genes and 151 (97.42%) samples will be considered for association analysis (Table 1).
The principal components (S1 Fig: supplementary data) and FST values (S2 Table: supplementary data) analysis showed that cases and controls were evenly dispersed (homogenous and samples did not cluster by phenotype); indicating that the population and subpopulation structure is not the driving force in our observations.
Genotyping results
Alleles for the sixteen genes and 64 SNPs and 4 indels analyzed were all in Hardy–Weinberg equilibrium (S2 and S3 Tables: supplementary data); suggesting random genetic exchange within the studied populations. The MAF varied considerably across SNP and indel (S2 Table: supplementary data) with the lowest MAF at rs11575934 in IL12RB1 (MAF = 0.0067) and the highest value at rs371194629 in HLA-G (MAF = 0.5).
The minor allele (T) of rs8062041 in HP appeared to be protective against HAT (p = 0.00024). An odds ratio (OR) of 0.359 (CI95 [0.20–0.63]) indicated low frequencies in cases compared to controls. This SNP is located in a copy number variation (CNV) essv41754 that spans both HP and HPR (Fig 2).
In addition, the minor alleles of IL4 and HLA-G also appeared protective (IL4: C rs2070874, uncorrected p = 0.047: and HLA-G: G rs1233330, uncorrected p = 0.011). The OR of 0.62 (CI95 [0.38–1.01]) for rs2070874-IL4 C and 0.2754 (CI95 [0.093–0.81]) for rs1233330 HLA-G also indicated low frequencies of the major allele in cases compared to controls. However, for HLA-G, the minor allele (A) of SNP rs17875389 had a higher frequency in cases than controls (p = 0.042). The OR of 2.29 (CI95 [0.97–5.39]) suggests that the A allele may increase the risk of developing HAT.
Of the 64 SNPs and 4 indels considered here, only four (SNP) of them belonging to three genes were associated with the development of HAT before Bonferroni correction (Table 2). After Bonferroni correction only one SNP (rs8062041 T/C) in HP was associated with HAT. The odds ratio of 0.359 suggests that the minor allele has a protective effect within the Cameroonian population with 95.3% (FDR) chance of this locus being associated with HAT (Table 2).
For the three remaining SNPs where the association was not significant after Bonferroni correction, our results show that the allele frequencies in cases and controls were not the same (Table 2) and that there is some possibility of an association with disease. FDR_BH is the probability of falsely rejecting the null hypothesis that allele frequencies are the same in cases and controls. rs1233330 and rs17875389 in HLA-G had FDR_BH values of 0.36 and 0.67 respectively; suggesting that there are 64% and 33% probabilities of an association between SNPs at these loci with HAT. For rs2070874 in IL4, the FDR_BH value of 0.722 suggests 27.8% chance of an association with HAT.
For the other genes (APOL1, IL4, IL6, IL10, IL8, TNFA and INFG) involved in immune response that have been previously investigated in HAT, our results revealed no statistical association with the disease within the Cameroonian population (S1 Table). No association was also observed with SNP and indel of APOL1 and all SNPs of HLA-A, IL1B, IL12B, CFH, IL12RB1, IL4R and MIF previously associated with the susceptibility to other infectious diseases (S2 Table: supplementary data).
Discussion
In this study we obtain good quality genotype data for a total of 64 SNPs and 4 indels in 16 genes to investigate associations with trypanosomiasis. Of these genes selected on the basis of their association with HAT or other infectious diseases, most were not statistically associated with HAT in Southern Cameroon.
The most important result of this study is the observation that the T allele of SNP rs8062041-HP with a p-value of 0.00024 (Bonferonni corrected p = 0.015) and an OR of 0.36 is associated with a lower risk (protective effect) of developing sleeping sickness. This SNP lies within intron 1–2 of HP of the CNV essv41754 that spans both HP and HPR transcripts (Fig 2). Although the biological significance of this CNV is not well understood, it is important to point out that HP and HPR have some biological similarities. Haptoglobin is involved in the scavenging of haem from lysed red blood cells. Trypanosome infections induce extensive lysis of red cells releasing haem which is scavenged by HP. In mice, the expression of the haptoglobin receptor (Cd163) on macrophages declines dramatically after infection with T. congolense [49] and is the earliest indicator of infection. HPR also binds haem but is not cleared from circulation after haemolysis. However, HPR is of particular interest because it plays a prominent role in the innate resistance of humans to most Trypanosoma species [50]. This innate resistance is linked to trypanosome lytic factors 1 and 2 (TLF1, TLF2) which are bound to a minor subclass of high-density lipoprotein (HDL) [51]. Both factors harbor APOL1, which is the trypanolytic component [52], and HPR which facilitates the uptake of APOL1 via trypanosome haptoglobin–hemoglobin receptors (HpHbR).
Interestingly, rs8062041-HP (T/C) is located on chromosome 16 (16q22.2) in the CNV essv41754 that spans both HP and HPR (Fig 2). Such genomic structural variants involving HP/HPR duplication have been reported with higher frequency in people of African descent [53, 7]. For instance, HP and HPR have been reported in 29 independent studies listed in the Database of Genome Variants [54]. T. b. gambiense protects itself against killing by APOL1 by reducing the abundance and affinity of the receptor for HPR [55]. If rs8062041, located in the CNV essv41754 spanning HP and HPR is correlated with CNV genotype, then an increase in HPR expression could drive increased uptake of APOL1 and parasite killing. As in other diseases such as heart disease, cancer, malaria and Crohn’s disease, polymorphism in HP could also have direct biological significance in HAT. Polymorphism in the haptoglobin gene may be associated with reduced cholesterol levels in the blood [56] and since cholesterol is specifically taken up by trypanosomes as a nutrient, any reduction in cholesterol might restrict parasite growth rate. There are numerous variants of HP, some of which may have arisen from gene conversion from HPR exons [56]. Single SNP tag these variants poorly (max r2 = 0.44), however SNP haplotypes can tag these variants efficiently (max r2 = 0.92) and are more strongly correlated with cholesterol levels than individual SNP [56]. High density genotyping of the HP/HPR locus will be required to understand the role of this locus in the response to trypanosome infection. Although rs8062041 is within HP, the known involvement of HPR in APOL1 mediated killing means that increased expression of HPR is another mechanism by which this SNP could be associated with the observed difference in likelihood of developing HAT. Our results showing an association between HAT and one SNP located within a CNV spanning HP and HPR duplication are not in line with results of Hardwick et al. [18] who observed no association with the HPR duplication allele and HAT in DRC. The difference between these results could be linked to the position of SNP within HP, the genetic diversity between the studied populations as well as the sampling methods. In our study, a case control approach was used while Hardwick et al. [18] used family-based sampling. Bresalier et al. [57] reported an association between polymorphisms at some HPR loci with an increasing risk of developing colon cancer. Similar associations were outlined by Tabak et al. [58] for HPR/APOL1 loci variations in hepatoma and leukemia. There are also examples of CNV mediating different susceptibilities to infectious diseases [59, 60].
Of the twelve SNPs of IL6 identified and investigated in our study, none of them revealed an association with HAT. However, with similar investigation on the same gene, Courtin et al. [15] showed a T allele of the IL6 (4339) SNP rs2069849 which was significantly (Bonferroni corrected p = 0.04) associated with a decreased risk of developing HAT in the DRC. This SNP was not genotyped in this study because it could not be multiplexed with the others in the panel. The discrepancy between our results and those of Courtin et al. [15] could be due to insufficient linkage between our marker SNP and rs2069849 and or genetic differences between the DRC and Cameroon populations. The study designs also differed; we used a case control approach while Courtin et al. [15] used a family-based design and our study was smaller.
It has been suggested that the G1 and G2 alleles of APOL1 which increase the risk of developing kidney disease are under selection because they confer resistance to HAT [17]. Our observations on APOL1 are consistent with Cooper et al. [22] who found no association with APOL1 G1 and G2 in a comparison of cases and active T.b.gambiense HAT.
Concerning HLA-G, our results showed a protective effect of developing HAT for the loci rs17875389 G/A (p = 0.0416 and OR of 2.291) and an increased risk effect for rs1233330 A/G (p = 0.01105 and OR of 0.2754). These results support those of Courtin et al. [14] who reported similar results for different SNPs of the same genes in the DRC.
The association with IL4 rs2070874 T/C (p = 0.00712 and OR of 0.6151) is the first time this has been observed in HAT although associations with IL4 have been observed in South American trypanosomiasis [61, 62]. The presence of IL4 in extravascular tissues promotes alternative activation of macrophages into M2 cells and inhibits classical activation of macrophages into M1 cells. This increase in repair macrophages (M2) is coupled with secretion of IL10 and TGFB that result in a diminution of pathological inflammation [63].
The results discussed above for IL4 and HLA-G are based on FDR_BH values should be used with caution because no association was found after correction for multiple testing. However, we were only able to collect a relatively small number of cases (56) for this study, despite conducting large-scale field surveys. Whilst our power calculations indicated that effects of the sizes observed could be detected with our relatively small number of samples, larger cohorts of well phenotyped cases and controls may be required to confirm these observations. Therefore, although the present data is only suggestive of an association, the finding of suggestive associations in multiple populations increases the probability that these are genuine associations with disease [64]. This challenge is precisely what the TrypanoGEN network, a consortium of partners in eight African and three European countries seeks to address. The network has collected from seven regions in six countries (Cameroon, Cote d’Ivoire, DRC, Malawi, Uganda, and Zambia) a total of 3301 samples from cases and controls to include in a genome-wide-association study [44] which will be used to test the hypotheses generated here.
Conclusion
The results of this study reveal an absence of association between HAT and several SNPs identified in genes previously associated with HAT within inhabitants of sleeping sickness foci of other African countries. An association between one SNP in HP and the susceptibility to HAT was revealed in inhabitants of sleeping sickness foci of Cameroon. Located within a CNV that spans both HP and HPR and given the known involvement of HPR in response to HAT, the association of rs8062041 with a CNV is the most plausible mechanism by which this SNP could be associated with protection against HAT. Our results reveal also that the association between host genetic determinants or gene polymorphisms and the susceptibility to T. b. gambiense infections may vary according to studied populations.
Supporting information
S1 Fig. Principal component analysis based on individual genotypes from the chronic HAT endemic area (foci) in Cameroon.
https://doi.org/10.1371/journal.pntd.0005979.s001
(TIF)
S1 Table. Candidate genes and SNPs identified and selected for this study.
https://doi.org/10.1371/journal.pntd.0005979.s002
(DOCX)
S2 Table. Fisher association analysis results.
https://doi.org/10.1371/journal.pntd.0005979.s003
(DOCX)
S3 Table. Fisher association of all loci genotype before quality control.
https://doi.org/10.1371/journal.pntd.0005979.s004
(DOCX)
Acknowledgments
We are grateful for the support and participation of the study populations included in this study as well as the technical team of the MINSANTE, Divisional Center for Diseases, PNLTHA. We also thank the anonymous reviewers for helpful comments on the manuscript.
Members of TrypanoGEN group: Hamidou Ilboudo, Harry Noyes, Julius Mulindwa, Magambo Phillip Kimuda, Mathurin Koffi, Justin Windingoudi Kabore, Bernadin Ahouty, Dieudonne Mumba Ngoyi, Olivier Fataki, Gustave Simo, Elvis Ofon, John Enyaru, John Chisi, Kelita Kamoto, Martin Simuunza, Vincent P. Alibu, Veerle Lejon, Vincent Jamonneau, Annette MacLeod, Mamadou Camara, Bruno Bucheton, Christiane Hertz-Fowler, Issa Sidibe, Enock Matovu.
References
- 1. Courtin D., Berthier D, Thevenon S., Dayo G-K., Garcia A., Bucheton B. (2008). Host genetics in African trypanosomiasis. Infect. Gen. Evol.8:229–238.
- 2.
World Health Organization, (2016). Human African Trypanosomiasis (HAT). Country Leadership and Collaboration on Neglected Tropical Diseases.
- 3. Jamonneau V., Ilboudo H., Kabore J., Kaba D., Koffi M., Solano P., Garcia A., Courtin D., Laveissiere C., Lingue K., Buscher P., Bucheton B. (2012). Untreated human infections by Trypanosoma brucei gambiense are not 100% fatal. PLoS Negl. Trop. Dis. 6, e1691. pmid:22720107
- 4. Blackwell J.M., Jamieson S.E. & Burgner D. (2009). HLA and Infectious Diseases. CLIN. MICROBIOL. REV., 22(2) 370–385 pmid:19366919
- 5. Carvalho F.M.C., Busser F.D., Freitas V.L.T., Furucho C.R., Sadahiro A., Kono A.S.G, Criado P.R., Moretti L.I, Sato P.K., Shikanai-Yasuda M.A.(2016). Polymorphisms on IFNG, IL12B and IL12RB1 genes and paracoccidioidomycosis in the Brazilian population. Infection, Genetics and Evolution 43:245–251. pmid:27223631
- 6. Francis P.J., Schultz D.W., Hamon S., Ott J., Weleber R.G., Klein M.L. (2007). Haplotypes in the Complement Factor H (CFH) Gene: Associations with Drusen and Advanced Age-Related Macular Degeneration. PLoS ONE 2(11): e1197. pmid:18043728
- 7. Kuhajda F.P., Katumuluwa A.I., and Pasternack G.R. (1989). Expression of haptoglobin-related protein and its potential role as a tumor antigen. Proc. Natl. Acad. Sci. 86;1188–1192. pmid:2465547
- 8. Minang J.T., Gyan B.A., Anchang J.K., Troye-Blomberg M., Perlmann H., Achidi E.A. (2004). Haptoglobin phenotypes and malaria infection in pregnant women at delivery in western Cameroon. Acta Trop. 90, 107–114. pmid:14739029
- 9. Quaye I.K.E., Ekuban F.A., Goka B.Q., Adabayeri V., Kurtzhals J.A.L., Gyan B., Ankrah N.A., Hviid L., Akanmori B.D. (2000). Haptoglobin 1–1 is associated with susceptibility to severe Plasmodium falciparum malaria. Trans. R. Soc. Trop. Med. Hyg. 94:216–219 pmid:10897372
- 10. Santos J.C., Ladeira M.S.P., Pedrazzoli J. Jr. & Ribeiro M.L. (2012). Relationship of IL-1 and TNF-á polymorphisms with Helicobacter pylori in gastric diseases in a Brazilian population. Braz. J. Med. Biol. Res., 45(9) 811–817. pmid:22714811
- 11. Sortica V.A., Cunha G.M., Ohnishi M.D.O., Souza J.M., Ribeiro-dos-Santos A.K.C., Santos N.P.C., Callegari-Jacques S.M., Santos S.E.B. and Hutz M.H. (2012). IL1B, IL4R, IL12RB1 and TNF gene polymorphisms are associated with Plasmodium vivax malaria in Brazil. Malaria Journal (2012), 11:409. pmid:23217179
- 12. Zhang D-F, Wang D., Li Y-Y., and Yao Y-G. (2014). Mapping genetic variants in the CFH gene for association with leprosy in Han Chinese. Genes and Immunity (2014) 15, 506–510; pmid:25030427
- 13. Bucheton B, MacLeod A, Jamonneau V. (2011). Human host determinants influencing the outcome of Trypanosoma brucei gambiense infections. Parasite Immunol 33: 438–447. pmid:21385185
- 14. Courtin D., Milet J., Sabbagh A., Massaro J.D., Castelli E.C., Jamonneau V., Bucheton B., Sese C., Favier B., Rouas-Freiss N., Moreau P., Donadi E.A., Garcia A., (2013). HLA-G 30 UTR-2 haplotype is associated with Human African trypanosomiasis susceptibility. Infection, Genetics and Evolution 17 (2013) 1–7. pmid:23541412
- 15. Courtin D., Milet J., Jamonneau V., Yeminanga C.S., Kumeso V.K., Bilengue C.M., Betard C., Garcia A., (2007). Association between human African trypanosomiasis and the IL6 gene in a Congolese population. Infect. Genet. Evol. 7, 60–68. pmid:16720107
- 16. Courtin D., Argiro L., Jamonneau V., N’Dri L., N’Guessan P., Abel L., Dessein A., Cot M., Laveissiere C., Garcia A., (2006). Interest of tumor necrosis factor-alpha 308 G/A and interleukin-10 592 C/A polymorphisms in human African trypanosomiasis. Infect. Genet. Evol. 6, 123–129. pmid:15894515
- 17. Genovese G., Friedman D.J., Ross M.D., Lecordier L., Uzureau P., Freedman B.I., Bowden D.W., Langefeld C.D., Oleksyk T.K., Uscinski Knob A.L., Bernhardy A.J., Hicks P.J., Nelson G.W., Vanhollebeke B., Winkler C.A., Kopp J.B., Pays E., Pollak M.R. (2010). Association of trypanolytic ApoL1 variants with kidney disease in African Americans. Science 329, 841–845. pmid:20647424
- 18. Hardwick J.R., Menard A., Sironi M., Milet J., Garcia A., C Sese C., Yang F., Fu B., David Courtin D.,. Hollox E.J. (2014). Haptoglobin (HP) and Haptoglobin-related protein (HPR) copy number variation, natural selection, and trypanosomiasis. Hum. Genet. 133:69–83. pmid:24005574
- 19. Ilboudo H., Berthier D., Camara M., Camara O., Kabore J., et al. (2012). APOL1 expression is induced by Trypanosoma brucei gambiense infection but is not associated with differential susceptibility to sleeping sickness. Infect Genet Evol 12: 1519–1523. pmid:22691369
- 20. Lejon V., Lardon J., Kenis G., Pinoges L., Legros D., Bisser S., N’Siesi X., Bosmans E., Buscher P. (2002). Interleukin (IL)-6, IL-8 and IL-10 in serum and CSF of Trypanosoma brucei gambiense sleeping sickness patients before and after treatment. Trans. R. Soc. Trop. Med. Hyg. 96, 329–333. pmid:12174791
- 21. Maclean L., Odiit M., Sternberg J.M. (2006). Intrathecal cytokine responses in Trypanosoma brucei rhodesiense sleeping sickness patients. Trans R Soc Trop Med Hyg 100: 270–275. pmid:16343570
- 22. Cooper A, Ilboudo H, Alibu VP, Ravel S, Enyaru J, Weir W, Noyes H, Capewell P, Camara M, Milet J, Jamonneau V, Camara O, Matovu E, Bucheton B, MacLeod A. APOL1 renal risk variants have contrasting resistance and susceptibility associations with African trypanosomiasis. Elife. 2017 May 24;6. pii: e25461. pmid:28537557
- 23. Penchenier L., Grébaut P., Ebo’o Eyenga V., Bodo J.M., Njiokou F., Binzouli J.J., Simarro P., Soula G., Herder S. (1999). Le foyer de la trypanosomiase humaine africaine de Campo (Cameroun) en 1998. Aspects épidémiologiques, état de l’endémie et comparaison des CATT 1.3 et CATT Latex dans le dépistage de masse. Bull. Soc. Path. Ex. 92, 185–190.
- 24. Simo G., Mbida J.A., Eyenga V.E., Asonganyi T., Njiokou F., Grébaut P. (2014). Challenges towards the elimination of Human African trypanosomiasis in the sleeping sickness focus of Campo in southern Cameroon. Parasites Vectors 7, 374. pmid:25129168
- 25. Assonganyi T., Bedifeh B.A., Ade S.S., Ngu J.L. (1994). An evaluation of the card agglutination test for trypanosomiasis (CATT) reagent in the Fontem sleeping Sickness focus, Cameroon. Afric. J. Med. Med. Sci. 23(1)[39–46].
- 26. Mélachio T.T., Simo G., Ravel S., De Meeûs T., Causse S., Solano P., Lutumba P., Asonganyi T., Njiokou F. (2011). Population genetics of Glossina palpalis palpalis from central African sleeping sickness foci. Parasites Vectors 4, 140. pmid:21767402
- 27. Simo G., Asonganyi T., Nkinin S.W., Njiokou F., Herder S. (2006). High prevalence of Trypanosoma brucei gambiense group 1 in pigs from the Fontem sleeping sickness focus in Cameroon. Veterinary Parasitology, 139: 57–66. pmid:16567049
- 28. Magnus E., Vervoort T., Van Meirvenne N. (1978). A card-agglutination test with stained trypanosomes (C.A.T.T.) for the serological diagnosis of T. b. gambiense trypanosomiasis. Ann. Soc. Belg. Méd. Trop. 58: 169–176. pmid:747425.
- 29.
WHO/World Health Organization, 1986. African trypanosomiasis: epidemiology and fight against. Technical Rapport Series, 739, 32pages.
- 30. Buscher P., Mumba N.D., Kabore J., Lejon V., Robays J., Jamonneau V., Bebronne N., Van der Veken W., Bieler S. (2009). Improved models of minianion exchange centrifugation technique (mAECT) and modified single centrifugation (MSC) for sleeping sickness diagnosis and staging. PLoS Negl. Trop. Dis. 3, e471. pmid:19936296
- 31. Jamonneau V., Bucheton B., Kabore J., Ilboudo H., Camara O., et al., (2010). Revisiting the immune trypanolysis test to optimise epidemiological surveillance and control of sleeping sickness in West Africa. PLoS Negl. Trop. Dis. 4: e917. pmid:21200417
- 32. Zhao J.H (2007). gap: Genetic Analysis Package. J Stat Softw [Internet]. 2007;23(8):1–18. Available from: http://www.jstatsoft.org/v23/i08
- 33. MacLean L., Odiit M., Sternberg J.M. (2001). Nitric oxide and cytokine synthesis in human African trypanosomiasis. J. Infect. Dis. 184, 1086–1090 pmid:11574928
- 34. MacLean L., Chisi J.E., Odiit M., Gibson W.C., Ferris V., Picozzi K., Sternberg J.M., (2004). Severity of Human African Trypanosomiasis inEast Africa is associated with geographic location, parasite genotype and host inflammatory cytokine response profile. Infect. Immun. 72, 7040–7044. pmid:15557627
- 35. Maclean L., Odiit M., Macleod A., Morrison L., Sweeney L., Cooper A., Kennedy P.E.G, and Sternberg J.M. (2007). Spatially and genetically distinct African Trypanosome virulence variants defined by host interferon-gamma response. J Infect Dis 196: 1620–1628 pmid:18008245
- 36. Vasilescu A., Heath S.C., Ivanova R., Hendel H., Do H., Mazoyer A., Khadivpour E., Goutalier F.X., Khalili K., Rappaport J., Lathrop G.M., Matsuda F., Zagury J.F. (2003). Genomic analysis of Th1-Th2 cytokine genes in an AIDS cohort: identification of IL4 and IL10 haplotypes associated with the disease progression. Genes Immun. 4, 441–449. pmid:12944981
- 37. Ilboudo H., Bras-Goncalves R., Camara M., Flori L., Camara O., Sakande H., Leno M., Petitdidier E., Jamonneau V., Bucheton B. (2014). Unravelling Human Trypanotolerance: IL8 is Associated with Infection Control whereas IL10 and TNFλ Are Associated with Subsequent Disease Development. PLoS Pathog 10(11): e1004469. pmid:25375156
- 38. Lajoie J., Hargrove J., Zijenah L.S., Humphrey J.H., Ward B.J., Roge M. (2006). Genetic variants in nonclassical major histocompatibility complex class I human leukocyte antigen (HLA)-E and HLA-G molecules are associated with susceptibility to heterosexual acquisition of HIV-1. J. Infect. Dis. 193, 298–301. pmid:16362895
- 39. Dias F.C, Castelli E.C., Collares C.V. A., Moreau P. and Donadi E.A. (2015). The role of HLA-G molecule and HLA-G gene polymorphisms in tumors, viral hepatitis, and parasitic diseases. pmid:25699038
- 40. Takahashi K., Koga K., Linge H.M., Zhang Y., Lin X., Metz C.N., Al-Abed Y., Ojamaa K., Miller E.J. (2009). "Macrophage CD74 contributes to MIF-induced pulmonary inflammation". Respir. Res. 10 (1): 33. PMC 2681459. pmid:19413900.
- 41. Das R., Mi-Sun K., Kim B.H., Jacob S.T., Subbianc S., Yaoa J., Lin L., Levy R., Murchison C.,. Burman W.J, Moore C.C., W. Scheld M., David J.R., Kaplan G., MacMicking J.D. and Bucala R. (2013). Macrophage migration inhibitory factor (MIF) is a critical mediator of the innate immune response to Mycobacterium tuberculosis. pmid:23882081
- 42. Jha A.N., Sundaravadivel P., Pati S.S., Patra P.K., Thangaraj K. (2013). Variations in ncRNA gene LOC284889 and MIF-794CATT repeats are associated with malaria susceptibility in Indian populations. Malar J 12: 345. pmid:24066864
- 43. Purcell S., Neale B., Todd-Brown K., Thomas L., Ferreira M.A.R, Bender D., Maller J., Sklar P., de Bakker P.I.W, Daly M.J. & Sham P.C. (2007). PLINK: a toolset for whole-genome association and population-based linkage analysis. Amer. J. Hum.Genet.81.
- 44. Ilboudo H, Noyes H, Mulindwa J, Kimuda MP, Koffi M, Kabore JW, et al. (2017). Introducing the TrypanoGEN biobank: A valuable resource for the elimination of human African trypanosomiasis. PLoS Negl Trop Dis 11(6): e0005438. pmid:28570558
- 45. Semagn K., Babu R., Hearne S., & Olsen M. (2014). Single nucleotide polymorphism genotyping using Kompetitive Allele Specific PCR (KASP): overview of the technology and its application in crop improvement. Molecular Breeding, 33(1), 1–14.
- 46. Anderson C.A., Pettersson F.H., Clarke G.M., Cardon L.R., Morris A.P., & Zondervan K.T. (2010). Data quality control in genetic-case control association studies. Nat. Protoc. 5:1564–1573. pmid:21085122
- 47. Benjamini Y., Hochberg Y. (1995): Controlling the false discovery rate: a practical and powerful approach to multiple testing. Journal of the Royal Statistical Society 1995:289–300.
- 48. Van Den Oord E.J.C.G. (2008): Controlling false discoveries in genetic studies. Am. J Med. Genet. B. Neuropsychiatr. Genet. 2008, 147B:637–644.
- 49. Noyes H.A., Alimohammadian M.H., Agaba M., Brass A., Fuchs H., Gailus-Durner V., et al. (2009). Mechanisms Controlling Anaemia in Trypanosoma congolense Infected Mice. PLoS ONE, 4(4), e5170. pmid:19365556
- 50. Pays E., Vanhollebeke B., (2008). Mutual self-defence. The trypanolytic factor story. Microbes Infect. 10, 985–989. pmid:18675374
- 51. Vanhollebeke B., Pays E., (2010). The trypanolytic factor of human serum: many ways to enter the parasite, a single way to kill. Mol. Microbiol. 76, 806–814. pmid:20398209
- 52. Vanhamme L., Paturiaux-Hanocq F., Poelvoorde P., Nolan D.P., Lins L., Van Den Abbeele J., Pays A., Tebabi P., Van Xong H., Jacquet A., Moguilevsky N., Dieu M., Kane J.P., De Baetselier P. Brasseur R. & Pays E. (2003). Apolipoprotein L-I is the trypanosome lytic factor of human serum. Nature 422, 83–87. pmid:12621437
- 53. Maeda N., McEvoy S.M., Harris H.F., Huisman T., Smithies O. (1986). Polymorphisms in the human haptoglobin gene cluster: chromosomes with multiple haptoglobin-related (Hpr) genes. Proc Natl Acad Sci 83:7395–7399. pmid:2876426
- 54. MacDonald J.R., Ziman R., Yuen R.K.C., Feuk L. & Scherer S.W. (2013). The Database of Genomic Variants: a curated collection of structural variation in the human genome. Nucleic Acids Research, 42(1), 986–992.
- 55. De Jesus E., Kieft R., Albright B., Stephens N. A. & Hajduk S. L. (2013). A Single Amino Acid Substitution in the Group 1 Trypanosoma brucei gambiense Haptoglobin-Hemoglobin Receptor Abolishes TLF-1 Binding. PLoS Pathogens, 9(4), e1003317. pmid:23637606
- 56. Boettger L.M., Rany M. Salem R.M., Handsaker R.E, Peloso G., Kathiresan S., Hirschhorn J. and McCarroll S.A. (2016). Recurring exon deletions in the haptoglobin (HP) gene associate with lower blood cholesterol levels. Nat Genet.; 48(4): 359–366. pmid:26901066
- 57. Bresalier R.S., Byrd J.C., Tessler D., Lebel J., Koomen J., Hawke D., Half E., Liu K.F., Mazurek N. (2004). A circulating ligand for galectin-3 is a haptoglobin-related glycoprotein elevated in individuals with colon cancer. Gastroent. 127:741–748.
- 58. Tabak S., Lev A., Valansi C., et al. (1997). "Transcriptionally active haptoglobin-related (Hpr) gene in hepatoma G2 and leukemia molt-4 cells." DNA. Cell Biol. 15 (11): 1001–7. pmid:8945641
- 59. Pelak K., Need A.C., Fellay J., Shianna K.V., Feng S., Urban T.J., Ge D., De Luca A., Martinez-Picado J., Wolinsky S.M., Martinson J., Jamieson B., Bream J., Martin M., Borrow P., McMichael A., Haynes B., Telenti A., Carrington M., Goldstein D., Alter G., Immunology NCfHAV. (2011). Copy number variation of KIR genes influences HIV-1 control. PLoS Biology 9:e1001208 pmid:22140359
- 60. Hardwick R.J., Amogne W., Mugusi S., Yimer G., Ngaimisi E., Habtewold A., Minzi O., Makonnen E., Janabi M., Machado L.R., Viskaduraki M., Mugusi F., Aderaye G., Lindquist L., Hollox E.J., Aklillu E. (2012). b-defensin Genomic Copy Number Is Associated With HIV Load and Immune Reconstitution in Sub-Saharan Africans. J Infect Dis 206:1012–1019. pmid:22837491
- 61. Arnez A.L.E., Venegas E.N., Ober C., Thompson E.E. (2011). Sequence variation in the IL4 gene and resistance to Trypanosoma cruzi infection in Bolivians. J Allergy Clin Immunol 127: 279–282, 282 e271-273. pmid:21211660
- 62. Florez O., Martin J., Gonzalez C.I. (2011). Interleukin 4, interleukin 4 receptor-alpha and interleukin 10 gene polymorphisms in Chagas disease. Parasite Immunol 33: 506–511. pmid:21729106
- 63. Hershey G.K., Friedrich M.F., Esswein L.A., Thomas M.L., Chatila T.A. (1997). The association of atopy with a gain-of-function mutation in the alpha subunit of the interleukin-4 receptor. N. Engl. J. Med. 337 (24): 1720–5. pmid:9392697
- 64. Kato C.D., Matovu E., Mugasa C.M., Nanteza A., Alibu V.P. (2016). The role of cytokines in the pathogenesis and staging of Trypanosoma brucei rhodesiense sleeping sickness. Allergy Asthma Clin Immunol. 2016;12:1. V. A, C. RD, E. JD, in Essentials of Glycobiolo. pmid:26734065