Relationship between Trypanosoma brucei rhodesiense genetic diversity and clinical spectrum among sleeping sickness patients in Uganda
BMC Research Notes volume 10, Article number: 518 (2017)
Human African trypanosomiasis (HAT) due to Trypanosoma brucei rhodesiense in East and southern Africa is reported to be clinically diverse. We tested the hypothesis that this clinical diversity is associated with a variation in trypanosome genotypes.
Trypanosome DNA isolated from HAT patients was genotyped using 7 microsatellite markers directly from blood spotted FTA cards following a whole genome amplification. All markers were polymorphic and identified 17 multi-locus genotypes with 56% of the isolates having replicate genotypes. We did not observe any significant clustering between isolates and bootstrap values across major tree nodes were insignificant. When genotypes were compared among patients with varying clinical presentation or outcome, replicate genotypes were observed at both extremes showing no significant association between genetic diversity and clinical outcome. Our study shows that T. b. rhodesiense isolates are homogeneous within a focus and that observed clinical diversity may not be associated with parasite genetic diversity. Other factors like host genetics and environmental factors might be involved in determining clinical diversity. Our study may be important in designing appropriate control measures that target the parasite.
Human African trypanosomiasis (HAT) or sleeping sickness is caused by extra-cellular protozoan parasites Trypanosoma brucei rhodesiense (East and southern Africa) and Trypanosoma brucei gambiense (West and central Africa). Previously, HAT due to T. b. rhodesiense has been classified as an acute disease with rapid progression [1, 2]. It is now clear that T. b. rhodesiense disease is clinically diverse, both within and across foci [3,4,5,6,7,8]. This observed clinical diversity has been partly attributed to a variation in trypanosome strains . A high degree of genetic diversity has been demonstrated among T. brucei stocks using highly polymorphic microsatellite markers [5, 10, 11]. Laboratory studies have shown that the parasite might be a strong driver for HAT associated pathology. When differences in pathology were compared in a mouse model using two distinct T. b. brucei strains (TREU-927 and STIB-247), mice infected with strain TREU-927 suffered a more severe clinical disease .
In a study comparing HAT clinical diversity among patients in Uganda and Malawi, patients in Uganda suffered a more severe disease . When parasite genotypes involved were analyzed, two sequence variants of the serum resistance antigen (SRA) gene were observed. However, these results were not conclusive since the SRA gene is not polymorphic to infer genetic diversity. When polymorphic microsatellite markers were used to compare disease in a restricted geographical area in Uganda (Tororo and Soroti), two distinct genotype clusters with varying disease response were observed . However, all these studies compared parasite genotypes between different HAT foci. With the observed clinical diversity within a localized HAT foci [3, 7] it is likely that the parasite might have a role to play.
We hypothesized that parasite genotypes within a single HAT focus are diverse. We genotyped trypanosomes directly from blood spotted FTA cards using 7 microsatellite markers following whole genome amplification. We further associated parasite genotypes to observed clinical diversity.
Materials and methods
Patients were recruited passively at Lwala hospital in northern Uganda from the year 2012–2014. Routine diagnosis of suspected HAT patients, was done by microscopic examination of wet and thick blood films . For inclusion in the study, patients needed to have a positive blood smear or those with a negative blood smear but with signs indicative of HAT with trypanosomes demonstrated in cerebrospinal fluid (CSF) or an elevation in CSF white blood cell count . Exclusion criteria were: terminally ill patients, children below 6 years and patients whose disease stage could not be ascertained. Physical examination was done by a medical officer and on each patient clinical data form demographic characteristics, signs and symptoms of HAT were recorded.
FTA card preparations and whole genome amplification
From each patient, approximately 200 µl of blood was spotted on the FTA card (Whatman). Whole genome amplification (WGA) was performed using the Ready-To-Go Genomiphi V3 DNA amplification kit (GE Healthcare, Sweden) as described previously [15, 16]. Briefly, 20 µl of cell lysis solution (400 mM KOH, 10 mM EDTA, 100 mM DTT) were added to 2-mm FTA discs. For WGA, 20 µl of the denatured cell lysate DNA was added to the Genomiphi V3 cake and samples incubated at 30 °C for 2 h followed by heating at 65 °C for 10 min with subsequent cooling at 4 °C. Three independent WGA reactions from the same sample were pooled and stored at − 20 °C until further use.
Polymerase chain reaction based genotyping and multi-locus genotype determination
We used 7 previously described microsatellite loci, Ch1/18, Ch2/5, Ch2/PLC, Ch4/M12C12, Ch3/5L5, Ch5/JS2 [5, 10, 17, 18] and M6C8 [19, 20] as shown in Additional file 1. All PCRs were performed in a final volume of 20 µl, containing: PCR buffer (50 mM tris-HCl (pH 9.0), 50 mM NaCl, 0.1 mg/ml BSA and 5 mM MgCl2), 200 µM of each dNTPs, 10 ng gDNA, 1 µm of forward and reverse primer and 1 unit of EconoTaq DNA polymerase (Lucigen, USA). PCR amplification conditions were, an initial denaturation at 95 °C for 3 min, followed by 45 cycles of 30 s at 95 °C, 30 s at 60–55 °C and a final elongation step for 20 min at 72 °C.
Allele size determination was done using a capillary based sequencer, the 3500xL Genetic Analyzer (Applied Biosystems) using GeneMapper Software v5.0 (Life Technologies). Multi-locus genotypes (MLGs) were defined by the specific combination of alleles across all loci (see Additional file 2).
Genetic and statistical analyses
We analyzed genetic data on allele frequencies, heterozygosity, allelic richness, and inbreeding coefficient (FIS) using GenAIEx v6.5 . A dendrogram based on MLGs was constructed using the neighbor-joining method in MEGA v6 . Data on clinical diversity was analyzed using IBM SPSS version 22. Numerical variables were summarized using medians. Univariate analysis to compare clinical diversity between those with the condition and those without was done using cross-tabulation with a Chi square or Fisher’s exact test.
Patient’s baseline characteristics and clinical diversity
A total of 25 HAT cases were recruited passively at Lwala hospital in Northern Uganda. A significant (P < 0.001) number of patients were diagnosed as late stage (18, 75%) compared to early stage patients (6, 25%). The male: female ratio for the patients was 2, with a median age of 23.8 years (Table 1). The observed clinical spectrum is presented in Table 1.
Microsatellite marker validation
We genotyped 25 trypanosome positive samples using 7 microsatellite markers. Marker Ch5/JS2 was monomorphic across all samples and was removed from further analysis. All markers except M6C8 exhibited an excess of heterozygosity (Ho > He, FIS < 0) with marker Ch4/M12C12 tending to heterozygosity fixation (FIS = − 1.00, Additional file 3). When allele frequencies were analyzed, marker Ch1/18 was the most polymorphic with 6 alleles in total, followed by Ch2/5 (4 alleles), Ch3/5L5/2 (3 alleles) the remaining markers had only 2 alleles each (Additional file 4).
Among the 25 genotyped samples, a total of 17 multi-locus genotypes (MLGs) were identified. Among the 17 MLGs, 44% (11/25) of the samples were associated with unique genotypes. Replicated genotypes were found in 56% (14/25) of the genotyped samples. MLG 4 occurred 4 times (16%) while all the other repeated MLGs were encountered only twice (8%). To demonstrate the relationship between T. b. rhodesiense isolates, we constructed a neighbor-joining tree using pairwise distance between the MLGs (Fig. 1). We did not observe any significant clustering across isolates and bootstrap values across major nodes were insignificant.
Correlation between clinical and genetic diversity
We carried out an evaluation to associate the scored MLGs with observed clinical diversity. However, due to the high heterogeneity of MLGs and the low power of the study, meaningful statistics could not be done. We compared the distribution of MLGs among patients with varying clinical presentation or outcome (Table 2). We did not observe any significant association between MLGs and clinical outcome. When genotypes were compared among patients with varying degree of neurological impairment as a measure of disease severity, MLGs 4, 13 and 14 were common in all extremes. Furthermore, when we compared clinical presentation between those with and without the clinical manifestation, repeated genotypes were observed at both extremes.
Sleeping sickness due to T. b. rhodesiense is described as clinically diverse, with varying disease severity and duration of illness across foci [4,5,6,7,8]. Similarly, a wide disease spectrum within the same HAT foci has been reported [3, 7] albeit, with reduced frequency. We tested the hypothesis that this variation in clinical diversity is associated with diversity in parasite genotypes.
In this study, we identified 17 multi-locus genotypes and among these, replicate genotypes were found in 56% of the isolates. With more than half of the samples constituted by repeated genotypes and the low bootstrap support for the phylogenetic tree, T. b. rhodesiense isolates in this focus appear homogenous with limited genetic diversity. Our findings are in agreement with a previous study in Soroti were T. b. rhodesiense isolates were homogeneous with replicate genotypes constituting 59% of the population . Indeed, when isolates from 2 related HAT foci in Uganda (Tororo and Soroti) were compared, the 2 populations could not be resolved with confidence due to the similar genotypes involved. This limited genetic diversity is consistent with previous studies describing T. b. rhodesiense stocks in Uganda as strictly clonal [10, 23].
In previous studies, a link between clinical diversity and parasite genetic diversity was proposed . We investigated this possibility using our data set. However, we did not observe any significant associations between genotypes and clinical outcome. In almost all cases replicate genotypes were found at all extremes of clinical outcome. Our results are in agreement with a previous T. b. gambiense study in Cote d’Ivoire that did not find any correlation between clinical diversity and parasite genotype . Similar to our study, the genotypes isolated from the latter study were homogeneous with limited genetic variation. Furthermore, when two zymodemes with varying clinical phenotypes (Zambezi and Busoga strains) were compared, the 2 strains were extremely similar . However, studies comparing T. b. rhodesiense isolates across different HAT foci associated clinical outcome to distinct parasite genotypes [5, 6, 8]. With the observed, genotypic homogeneity within a single HAT focus in this study, our results are expected. Due to the high polymorphic nature of microsatellite markers, a high number of MLGs appear even with minimal genetic differentiation. Indeed, with the high mutation rates within microsatellite loci [26, 27] related genes might be seen as slightly different  without necessarily changing virulence. Thus, although these markers have been shown to be perfect markers for inferring population structure, they pose a challenge in associating parasite diversity to clinical outcome.
Our study shows that T. b. rhodesiense isolates are homogeneous within a focus and that observed clinical diversity may not be related to parasite genetic diversity. It is likely that other factors like host genetic polymorphism and unknown environmental factors might be involved in driving clinical diversity. Due to the limited number of patients in this preliminary study, research studies recruiting large patient groups from different geographical locations are recommended for more conclusive results.
Limitations of the study
As a result of the low incidence of disease, we were unable to achieve a sample size greater than 25 patients. Further, this small sample size allowed us to compare only 18 late stage patients and 6 early stage patients. Consequently, this limited sample could not allow more conclusive statistics to be done. Indeed, correlation studies between observed genotypes and clinical presentation would have benefited from a larger sample size so as to cater for both disease extremes. To try and make the study more conclusive, we used a total of 7 markers so as to pave way for further investigations.
glasgow coma scale
human African trypanosomiasis
polymerase chain reaction
serum resistance antigen
white blood cell
whole genome amplification
Barrett MP, Burchmore RJ, Stich A, Lazzari JO, Frasch AC, Cazzulo JJ, et al. The trypanosomiases. Lancet. 2003;362(9394):1469–80.
Foulkes J. The six diseases WHO. Hum Trypanosomiasis Afr BMJ. 1981;283(6300):1172–4.
Kato CD, Nanteza A, Mugasa C, Edyelu A, Matovu E, Alibu VP. Clinical profiles, disease outcome and co-morbidities among T. b. rhodesiense sleeping sickness patients in Uganda. PLoS ONE. 2015;10(2):e0118370.
Kuepfer I, Hhary EP, Allan M, Edielu A, Burri C, Blum JA. Clinical presentation of T.b. rhodesiense sleeping sickness in second stage patients from Tanzania and Uganda. PLoS Negl Trop Dis. 2011;5(3):e968.
MacLean LM, Odiit M, MacLeod A, Morrison L, Sweeney L, Cooper A, et al. Spatially and genetically distinct African trypanosome virulence variants defined by host interferon-γ response. J Infect Dis. 2007;196(11):1620–8.
MacLean LM, Chisi JE, Odiit M, Gibson WC, Ferris V, Picozzi K, et al. Severity of human African trypanosomiasis in East Africa is associated with geographic location, parasite genotype, and host inflammatory cytokine response profile. Infect Immun. 2004;72(12):7040–4.
MacLean LM, Odiit M, Chisi JE, Kennedy PG, Sternberg JM. Focus-specific clinical profiles in human African trypanosomiasis caused by Trypanosoma brucei rhodesiense. PLoS Negl Trop Dis. 2010;4(12):e906.
Smith DH, Bailey JW. Human African trypanosomiasis in south-eastern Uganda: clinical diversity and isoenzyme profiles. Ann Trop Med Parasitol. 1997;91(7):851–6.
Hide G, Buchanan N, Welburn S, Maudlin I, Barry JD, Tait A. Trypanosoma brucei rhodesiense: characterisation of stocks from Zambia, Kenya, and Uganda using repetitive DNA probes. Exp Parasitol. 1991;72(4):430–9.
Duffy CW, MacLean L, Sweeney L, Cooper A, Turner CMR, Tait A, et al. Population genetics of Trypanosoma brucei rhodesiense: clonality and diversity within and between foci. PLoS Negl Trop Dis. 2013;7(11):1–11.
Echodu R, Sistrom M, Bateta R, Murilla G, Okedi L, Aksoy S, et al. Genetic diversity and population structure of Trypanosoma brucei in Uganda: implications for the epidemiology of sleeping sickness and Nagana. PLoS Negl Trop Dis. 2015;9(2):e0003353.
Morrison LJ, McLellan S, Sweeney L, Chan CN, MacLeod A, Tait A, et al. Role for parasite genetic diversity in differential host responses to Trypanosoma brucei infection. Infect Immun. 2010;78(3):1096–108.
Chappuis F, Loutan L, Simarro P, Lejon V, Büscher P. Options for field diagnosis of human African trypanosomiasis. Clin Microbiol Rev. 2005;18(1):133–46.
WHO. Control and surveillance of African trypanosomiasis: report of a WHO expert committee. WHO Tech Rep Ser. 2013;984:1–250.
Kato CD, Alibu VP, Nanteza A, Mugasa CM, Matovu E. Population genetic structure and temporal stability among Trypanosoma brucei rhodesiense isolates in Uganda. Parasites vectors. 2016;9(1):1.
Morrison LJ, Mccormack G, Sweeney L, Likeufack AC, Truc P, Turner CM, et al. Use of multiple displacement amplification to increase the detection and genotyping of Trypanosoma species samples immobilized on FTA filters. Am J Trop Med Hyg. 2007;76(6):1132–7.
MacLeod A, Tweedie A, McLellan S, Taylor S, Cooper A, Sweeney L, et al. Allelic segregation and independent assortment in T. brucei crosses: proof that the genetic system is Mendelian and involves meiosis. Mol Biochem Parasit. 2005;143(1):12–9.
Morrison LJ, Tait A, McCormack G, Sweeney L, Black A, Truc P, et al. Trypanosoma brucei gambiense type 1 populations from human patients are clonal and display geographical genetic differentiation. Infect Genet Evol. 2008;8(6):847–54.
Biteau N, Bringaud F, Gibson W, Truc P, Baltz T. Characterization of Trypanozoon isolates using a repeated coding sequence and microsatellite markers. Mol Biochem Parasit. 2000;105(2):187–202.
Koffi M, Solano P, Barnabé C, De Meeûs T, Bucheton B, Cuny G, et al. Genetic characterisation of Trypanosoma brucei sl using microsatellite typing: new perspectives for the molecular epidemiology of human African trypanosomosis. Infect Genet Evol. 2007;7(6):675–84.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research—an update. Bioinformatics. 2012;28(19):2537–9. https://doi.org/10.1093/bioinformatics/bts460.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9. https://doi.org/10.1093/molbev/mst197.
MacLeod A, Tweedie A, Welburn SC, Maudlin I, Turner CMR, Tait A. Minisatellite marker analysis of Trypanosoma brucei: reconciliation of clonal, panmictic, and epidemic population genetic structures. Proc Natl Acad Sci. 2000;97(24):13442–7.
Jamonneau V, Garcia A, Ravel S, Cuny G, Oury B, Solano P, et al. Genetic characterization of Trypanosoma brucei gambiense and clinical evolution of human African trypanosomiasis in Côte d’Ivoire. Tropical Med Int Health. 2002;7(7):610–21.
Goodhead I, Capewell P, Bailey JW, Beament T, Chance M, Kay S, et al. Whole-genome sequencing of Trypanosoma brucei reveals introgression between subspecies that is associated with virulence. MBio. 2013;4(4):e00197. https://doi.org/10.1128/mBio.00197-13.
Jin L, Macaubas C, Hallmayer J, Kimura A, Mignot E. Mutation rate varies among alleles at a microsatellite locus: phylogenetic evidence. Proc Natl Acad Sci. 1996;93(26):15285–8.
Schlötterer C. Evolutionary dynamics of microsatellite DNA. Chromosoma. 2000;109(6):365–71.
Rousset F. Equilibrium values of measures of population subdivision for stepwise mutation processes. Genetics. 1996;142(4):1357–62.
VPA and EM received the grant; CDK, VPA, AN, CMM and EM designed laboratory protocols; CDK, AN and CMM analyzed the data; EM, CDK and CMM wrote the paper; VPA and AN reviewed the paper. All authors read and approved the final manuscript.
We acknowledge writing support from the Makerere University–Uganda Virus Research Institute Centre of Excellence in Infection and Immunity Research and Training (MUII-Plus). We thank the Consortium for Advanced Research Training in Africa (CARTA) for the support in sample collection.
The authors declare that they have no competing interests.
Availability of data and materials
All data generated or analyzed during this study are included in this published article (and its Additional files).
Consent for publication
Ethics approval and consent to participate
The study protocol was reviewed by the Institutional Review Board (IRB) of the Vector Control Division (REF: VCD-IRC/021) and approved by the Uganda National Council for Science and Technology (REF: HS 1089) upon recommendation by the IRB. Before recruitment, patients were given written and verbal information about the project objectives in local language so as to give written informed consent. For patients below 16 years, written consent was given by the legal guardian, but they were also asked for assent. Samples used in the study were remnants from routine normal diagnostic procedures required to guide treatment as approved under the study protocol.
This work was supported within the framework of the consortium Afrique One “Ecosystem and Population Health: Expanding Frontiers in Health” (http://www.afriqueone.net/). Afrique One is funded by the Wellcome Trust (087535/Z/08/A). The funder had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Kato, C.D., Mugasa, C.M., Nanteza, A. et al. Relationship between Trypanosoma brucei rhodesiense genetic diversity and clinical spectrum among sleeping sickness patients in Uganda. BMC Res Notes 10, 518 (2017). https://doi.org/10.1186/s13104-017-2860-x
- Human African trypanosomiasis
- Sleeping sickness
- Clinical diversity
- Genetic diversity
- Multi-locus genotypes
- Microsatellite markers