- Research article
- Open Access
Predominance of Uganda genotype of Mycobacterium tuberculosis isolated from Ugandan patients with tuberculous lymphadenitis
BMC Research Notes volume 8, Article number: 398 (2015)
In Uganda, the emerging Uganda genotype of Mycobacterium tuberculosis is the most common cause of pulmonary tuberculosis (PTB), and accounts for up to 70 % of isolates. Extrapulmonary TB (EPTB) is less studied in Uganda.
Molecular characterization using deletion analysis and spoligotyping was performed on 121 M. tuberculosis isolates from lymph node fine needle biopsy aspirates of Ugandan patients with tuberculous lymphadenitis. The evolutionary relationships and worldwide distribution of the spoligotypes were analyzed.
Mycobacterium tuberculosis was the only cause of EPTB in this study. The T2 sublineage was the most predominant lineage and the Uganda genotype was the dominant genotype. There were 54 spoligotype patterns among the 121 study isolates. The dominant spoligotypes were shared international types (SIT) SIT420, SIT53, SIT 135, SIT 128 and SIT590 in descending order. All but SIT420 were previously reported in pulmonary TB in this setting. The phylogenetic analysis showed a long descendant branch of spoligotypes belonging to the T2-Uganda sublineage containing specifically SITs 135, 128 and 420.
In most cases, the spoligotypes were similar to those causing PTB, but the Uganda genotype was found to be less common in EPTB than previously reported for PTB in Uganda. The phylogenetic analysis and the study of the worldwide distribution of clustered spoligotypes indicate an ongoing evolution of the Uganda genotype, with the country of Uganda at the center of this evolution.
Tuberculosis (TB) is a leading infectious cause of morbidity and mortality, responsible for an estimated annual 8.2 million new cases and 1.4 million deaths globally . About one quarter of TB cases occur in Africa . Uganda with an incidence of 299 and a mortality of 84 cases per 100,000 per year, is ranked 16th among the 22 countries with the heaviest burden of TB . Although pulmonary TB (PTB) is the major type of TB, extrapulmonary TB (EPTB)  is associated with a high morbidity and mortality  and requires special attention . TB lymphadenitis is the commonest extrapulmonary manifestation of EPTB . Of 4072 patients seen in 2012 at the Mulago National Referral Hospital (MNRH) Tuberculosis Clinic, the largest TB treatment center in Uganda, 902 (22 %) had EPTB . EPTB also offers diagnostic challenges because of the paucibacillary nature of EPTB specimens, necessitating culture and molecular analysis for efficient diagnosis.
Mycobacterium tuberculosis has evolved into genetically diverse lineages and sublineages , and the genetic variety has ramifications including significant variation in virulence . The development of molecular techniques allows for the identification and tracking of individual strains of M. tuberculosis . In Uganda, the emerging Uganda genotype of M. tuberculosis is the prevalent cause of PTB, and accounts for up to 70 % of isolates . In a recent study of tuberculous lymphadenitis patients, we found that the Uganda genotype was also a major aetiological agent in EPTB .
Our objective was to analyze the genetic diversity of the clinical isolates from this study , and to compare them with the genotypic patterns of PTB isolates reported previously in this setting  in order to understand the epidemiological and phenotypic characteristics of M. tuberculosis strains causing EPTB.
The study was reviewed and approved by the Institutional Review Board of Makerere University School of Medicine and the Uganda National Council for Science and Technology. Written informed consent to obtain samples as well as to use isolates from the samples for studies was obtained from all enrolled study participants or their legal guardians.
Study design and population
Isolates of M. tuberculosis complex were studied from patients presenting with superficial persistent lymph node enlargement, at the Fine Needle Aspiration Clinic, at MNRH, Kampala. Altogether 121 lymph nodes isolates, each representing one patient, collected between February 2010 and July 2012 were included. EPTB patients were included with or without pulmonary clinical features.
The patients’ clinical, radiological and pathological features were described previously . Of the patients in the study, 58 (47.9 %) were males with a mean age of 30.3 years and 62 (52.1 %) were females with a mean age of 27.1 years. The gender of one patient was not recorded. Seventy-five (65.8 %) were HIV seropositive while 39 (34.2 %) were HIV negative. Seven (5.9 %) did not consent to HIV testing.
To confirm that all the isolates were mycobacteria, polymerase chain reaction (PCR) was performed using 16s reverse and 16s forward primers (Integrated DNA Technologies) to target the 16s rRNA region with a conserved sequence that is typical for genus Mycobacterium [13–15].
Furthermore, using an in-house PCR, we confirmed presence of M. tuberculosis complex members by amplification of the insertion sequence IS6110 with aid of reverse and forward IS6110 primers (Integrated DNA Technologies). Presence of bands on the gel of a size of approximately 500 bp signified positive results .
Region of difference analyses
The lineages and sublineages of M. tuberculosis isolates were determined by PCR based genomic deletion analysis for the presence or absence of specific regions of difference (RD), RD1, RD4, RD9, RD14, RD724, RD750 and RD105 using standard protocol . RD 9 ruled out other species and confirmed that the cases were M. tuberculosis, RD4 and RD 14 ruled out M. bovis, the RD724 deletion is characteristic of Uganda genotype while RD105 deletion is specific for Beijing strains. The analyses were performed using a previously described method .
All M. tuberculosis complex strains were assayed by spoligotyping using standard protocols  and following manufacturer’s instructions (reagents from Ocimum Biosolution, custom MasterMix from ABgene). The presence of spacers was visualized on film as black squares after incubation with streptavidin-peroxidase and ECL chemiluminescence detection reagents (RPN 2105 Amersham, GE Healthcare Bio-sciences). The spacer hybridization patterns were translated into binary and octal format as previously described . The 43-digit binary code was converted to 15-digit octal code (base 8, having the digits 0–7) .
The binary codes of the isolates were entered into the SITVIT2 database of the Pasteur Institute of Guadeloupe and assigned specific shared international spoligotype signatures (SIT) according to the SITVIT2 data base .
Identification of Uganda genotype
The evolutionary relationships between the spoligotypes of the 121 isolates was studied by (a) a Minimum Spanning Tree (MST) constructed on all isolates, (b) a spoligoforest tree drawn as a hierarchical layout and (c) a spoligoforest tree drawn using the Fruchterman–Reingold algorithm. Both spoligoforests were drawn using the spolTools software (http://www.emi.unsw.edu.au/spoltools/) [22, 23].
The worldwide distribution of clustered spoligotypes found in this study was further investigated using the SITVIT2 database, and was recorded for regions and countries representing ≥3 % of a given SIT as compared to their total number in the SITVIT2 database. The various macro-geographical regions and sub-regions were defined according to the specifications of the United Nations (http://en.wikipedia.org/wiki/ISO_3166-1_alpha-3).
Clinical-demographic data was entered into the computer using Epidata 3.1 software and then exported to SPSS version 21 for analysis. The spoligotyping data were digitized and analyzed with the BioNumerics software, version 5.0 (Applied Maths, Kortrijk, Belgium).
Statistical association between the dominant strain types and HIV sero-status were analyzed in a 2 × 2 contingency table using Fisher’s exact test. A p value of <0.05 was considered statistically significant. The independent T test was used to determine whether there was a statistically significant difference between the female and male mean ages.
Genus and species identification
All of the 121 isolates were confirmed to be M. tuberculosis complex based on 16S-rRNA as well as the presence of the insertion sequence IS6110 that is highly conserved in all members of M. tuberculosis complex.
All of the 121 isolates had spoligotypes characteristic of M. tuberculosis. The spoligotype octal code outcome results were compared to the SIT and lineages/sublineages described in the SITVIT2 database. Spoligotyping revealed 52 distinct spoligotype patterns, of which 32 patterns matched a pre-existing SIT in the SITVIT2 database, whereas 4 SITs were newly created either within the present study or after a match with an orphan in the database (Table 1). Sixteen patterns were true orphans (Additional file 1: Table S1).
The T2 sublineage was predominant comprising 55 (46 %) isolates, of which 46 were T2-Uganda, characterized by the absence of spacers 33–36 and spacer 40 and 43 , while 2 isolates had the controversial signature T2/LAM3. The T2 sublineage was followed in frequency by the T1 sublineage with 19 (16 %) isolates. The Central Asian Strain (CAS) family comprised 16 (18 %) isolates, of which CAS1-Dehli was the most common with 10 isolates followed by CAS1-Kili with 5 isolates. Other sublineages were EAI5 with 3 isolates, Beijing 3, LAM11-ZWE 3, MANU3 3 and LAM1 2. Seven spoligotypes could not be assigned to any sublineage.
The two predominant spoligotypes in our sample were SIT420 (T2-Uganda) with 13 (10.7 %) and SIT53 (T1) with 11 (19.8 %) isolates. Other significant spoligotypes identified in the sample were SIT128 (T2-Uganda) with 9 (7.4 %) and SIT135 (T2-Uganda) with 10 (8.3 %) isolates (Table 1). Sixteen isolates had no assigned SIT according to the SITVIT2 database and were true orphans (Additional file 1: Table S1). None of the predominant spoligotypes SIT420, SIT53 and SIT135 was significantly associated with HIV infection. No association was found between the major spoligotypes and HIV infection (Additional file 2: Table S2).
Of the 12 SIT420 patients tested for HIV, 8/74 (10.8 %) were HIV negative whereas 4/38 (10.5 %) were HIV positive (p = 1.00). One patient did not consent to HIV testing. Of ten SIT135 patients tested for HIV, 6/74 (8.1 %) were HIV negative while 4/38 (10.5 %) were HIV positive (p = 0.732). Ten SIT53 patients tested for HIV, 9/74 (12.2 %) were HIV negative and 1/38 (2.6 %) were HIV positive (p = 0.16). One SIT53 patient did not consent to HIV testing.
No correlation between a particular spoligotype and gender or age was found (Additional file 3: Table S3).
RD9 analysis confirmed that all isolates were M. tuberculosis, while RD4 and RD14 analysis ruled out M. bovis. RD105 analysis identified three Beijing isolates. RD724 analysis revealed 55 (46 %) isolates to be Uganda genotypes, that is all the 55 isolates with a T2 spoligotype signature, lacking spacers 33–36 and spacer 40 and/or 43.
Two of the isolates, designated SIT125, which were deleted for RD724, were initially assigned to the LAM lineage according to the SITVITWEB database, although the spoligotype signature was a typical T2. Their designation in the SITVIT2 database has now been changed to T2/LAM3.
Phylogenetic analysis showed the evolutionary relationships between the spoligotypes of the 121 isolates (Fig. 1). One may notice in Fig. 1b, a long descendant branch of spoligotypes belonging to T2-Uganda sublineage containing specifically SITs 135, 128 and 420, which supports an ongoing evolution of this sublineage in Uganda.
We also attempted to describe the worldwide distribution of predominant SITs encountered in this study. As shown in Table 2, we observed that the isolates of the Uganda genotype were less frequently observed in countries outside Uganda (0–9 countries) than isolates of other sublineages (6–13 countries).
In this study of the genotypic diversity of M. tuberculosis strains causing EPTB in Uganda, EPTB was found to be caused exclusively by M. tuberculosis. The Uganda genotype was the predominant genotype comprising 46 % of the isolates.
The Uganda genotype is defined by the RD724 deletion in combination with the lack of spacers 33–36 and spacer 40 and/or 43 by spoligotyping [20, 21]. Interestingly all isolates of the T2 sublineage were of the Uganda genotype, indicating that the RD724 deletion may be a marker of a majority of the T2 sublineage, a matter that needs to be further explored.
In a phylogenetic analysis SIT420 (T-2 Uganda), the most common spoligotype in this study was at the tip of a long descendant branch of spoligotypes belonging to the T2-Uganda sublineage, and this supports the ongoing evolution of this sublineage in Uganda. The fact that SIT420 is only reported in three other countries (USA, Tanzania and Belgium) in the SITVIT2 database supports the theory that Uganda is at the center of this evolution.
Following SIT420, SIT53 was the most prevalent spoligotype. SIT53 has been reported in several countries in sub-Saharan Africa. In a study performed in Cote d´Ivoire, SIT53 was the most predominant spoligotype in PTB retreated patients . In a study from South Africa, SIT53 was the most frequent spoligotype and was associated with mixed infections . However, in three previous studies in Uganda SIT53 was either rare [27, 28] or absent .
Although the Uganda genotype is the predominant genotype, it appears to be a less frequent aetiological agent in EPTB than in PTB. In a study of M. tuberculosis complex isolates from PTB done in the same setting , the Uganda genotype accounted for 70 % of isolates compared to 46 % seen in this current study. A more recent longitudinal study of sputum isolates in central Uganda also showed predominance of the Uganda genotype of M. tuberculosis at 63 % over a period from 1992 to 2009 , and in a study of sputum isolates done in Southwestern Uganda, 59.2 % of the isolates were of Uganda genotype . In another study done on PTB and EPTB subjects in a rural central Uganda district hospital, 23 different spoligotypes were detected comprising predominantly T2 with subtypes Uganda I and Uganda II , which are both characterized by lack of hybridization to spacer 40, while in addition, strains of genotype Uganda I also lack spacer 43 .
The combination of spoligotyping and deletion analysis revealed one incongruence: two isolates, designated SIT125 (signature being loss of spacers 1–24, 33–36 and 40) which were deleted for RD724, were originally assigned to the LAM lineage according to the SITVIT2 database. By spoligotyping, the Euro-American lineage consists of the sublineages T, Haarlem, LAM, S and X. The spoligotype patterns of isolates from different lineages may converge by deletion of single or contiguous spacers in the direct repeat (DR) locus . Although the phenomenon of homoplasy in the DR locus leading to convergent evolution to identical spoligotypes is considered a rare event , it argues against the use of spoligotyping for establishing phylogenetic relationships within the Euro-American lineage .
Indeed SIT125 with its characteristically ‘abridged’ profile could be linked either to T2 or LAM3 as well as an underlying S-lineage signature. It was initially classified as T2 in the SpolDB4 database and later reclassified as LAM3 in SITVITWEB, both being feasible based strictly on spoligotyping signatures. Although its definitive classification remains controversial, parallel data show its high phylogeographical specificity for Bulgaria . In that paper, a possible relatedness with the SIT34 spoligotype (prototype of the S family) was suggested. Furthermore, an absence of LAM-specific IS6110 insertion, which suggested that SIT125 (at least those from Bulgaria) were most probably not LAM. However, a study from Brazil on RDRio  showed that although 4/5 SIT125 strains were not LAM, 1/5 strain was LAM/RD-Rio. The SIT125 designation in the SITVIT2 database has now been changed to the controversial signature of T2/LAM3, while the exact phylogenetic position of SIT125 is debated and possible homoplasy in spoligotyping defined lineages is considered. For this reason, in this current study we have chosen to assign the two SIT125 isolates to the T2 lineage.
A previous PTB study from Uganda  indicated SIT128, SIT135, SIT52, SIT590 and SIT125 as the major SITs. SIT420, the predominant spoligotype in our study, was not identified in an urban area of Kampala  and was found only in low proportion in Southwestern Uganda  in patients with PTB. However, in a study from central Uganda  of both PTB and EPTB (lymph node isolates) SIT420 was the most prevalent with a cluster rate of 14.8 %, and mostly recovered from rural residents. In a study of sputum isolates from the neighboring country of Rwanda, SIT420 was the third most common spoligotype .
SIT52 constituted 3.3 % of the isolates compared to 7.6 % (26/334) seen in PTB in central Uganda  and 4.8 % (6/125) in Southwestern Uganda . In Ethiopia, SIT54, SIT53 and SIT149 were the most dominant spoligotypes in EPTB while T1 was the most dominant sublineage . Previous evidence from Italy however indicates that CAS lineages are associated with EPTB . A more recent study of EPTB conducted in India also showed spoligotypes belonging to the CAS family (57.27 %) were predominant .
There is evidence of genotypic linkage of strains isolated from patients from a distinctive geographical region regardless of whether the patients acquired the infection from their current locality or had reactivated disease from their native country, this suggests that M. tuberculosis lineages may be adapted to particular genetic, cultural or environmental characteristics of the host . Most of the genotypes in this study were similar to those previously found in PTB isolates in the same area , supporting the concept that EPTB and PTB are caused by the same type of strains and that mainly host factors determine whether PTB will transient through the lung to establish as chronic lymphatic disease . The quantitative differences in spoligotype patterns could partly be due to the fact that Mulago Hospital is a national referral hospital, which has the largest TB diagnostic and treatment facilities in the country.
The lower prevalence of the Uganda genotype in this current study of EPTB isolates, compared to studies of PTB isolates in the same location indicates that other genotypes may be more prone to cause EPTB. The fact that the predominant SIT420 and SIT53 in this study was not seen in the PTB isolates previously studied  is another intriguing difference between EPTB and PTB in this setting. The observations merit future prospective studies in which isolates from patients from the same rural and urban region with EPTB and PTB are compared head to head at the same time.
Evidence from previous studies indicate that pathogen characteristics also determine whether patients present with PTB or EPTB [39, 40]. Thus, strains of the Euro-American lineage (to which the Uganda genotype belongs) may be more likely to cause pulmonary disease and less capable of extrapulmonary dissemination .Studies in the United States indicate that Indo-Oceanic and East Indian lineages of MTB were associated with a higher frequency of exclusively EPTB . It is noteworthy that none of the predominant spoligotypes or sublineages was significantly associated with HIV infection.
In Uganda M. bovis control efforts are unsatisfactory which has increased the incidence of this infection in the population, and therefore the complete absence of M. bovis in the current study is worth noting. This absence of M. bovis was also seen in studies from Ethiopia . Even in pastoralists  no M. bovis was isolated from lymph node fine needle biopsies despite close contact between humans and livestock, including consumption of raw milk and meat. But in another study from Southern Ethiopia, six of 35 (17.1 %) PCR positive patients with tuberculous lymphadenitis were found to have M. bovis . Findings from a study in Tanzania of cervical adenitis patients confirmed that 10.8 % of the cases were due to M. bovis . Thus, there is increasing evidence that M. bovis is at least a minor cause of human lymphadenitis in African countries.
In conclusion, EPTB was caused exclusively by M. tuberculosis and the Uganda genotype was the predominant genotype, although at a lower frequency than in studies of PTB in the same location. The spoligotypes were in most cases similar to those causing PTB with the exception of SIT420 which predominated in EPTB but are mainly lacking in other studies of PTB. Exploring the diversity of M. tuberculosis is essential in appreciating the clinical manifestations and pathogenesis of the disease, as well as the development of new diagnostic strategies.
Mtb lineages are geographically restricted and may adapt to local human populations. Effective vaccine candidates may in future have to be evaluated against prevailing Mtb genotypes and host genetic backgrounds.
human immunodeficiency virus
mycobacteria other than tuberculosis
fine needle aspiration
mycobacterial growth indicator tube
Central Asian Strain
direct repeat locus
polymerase chain reaction
spoligotype international type
Geneva WHO. Global tuberculosis report. http://www.whoint/tb/publications/global_report/gtbr12_executivesummary.pdf. 2012.
WHO. Global tuberculosis report. WHO/HTM/TB/2010. 2010;7:15–32.
Sharma SK, Mohan A. Extrapulmonary tuberculosis. Indian J Med Res. 2004;120(4):316–53.
Kang’ombe CT, Harries AD, Ito K, Clark T, Nyirenda TE, Aldis W, et al. Long-term outcome in patients registered with tuberculosis in Zomba, Malawi: mortality at 7 years according to initial HIV status and type of TB. Int J Tuberc Lung Dis. 2004;8(7):829–36.
Fanning A. Tuberculosis: 6. Extrapulmonary disease. CMAJ. 1999;160(11):1597–603.
Handa U, Mundi I, Mohan S. Nodal tuberculosis revisited: a review. J Infect Dev Ctries. 2012;6(1):6–12.
Clinic MHT. Tuberculosis treatment centre. Mulago Hospital TB clinic patient register. 2012.
Brudey K, Driscoll JR, Rigouts L, Prodinger WM, Gori A, Al-Hajoj SA, et al. Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol. 2006;6:23 (epub 2006/03/08).
Gagneux S. Genetic diversity in Mycobacterium tuberculosis. Curr Top Microbiol Immunol. 2013;16:16.
Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, et al. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol. 1997;35(4):907–14 (epub 1997/04/01).
Asiimwe BB, Koivula T, Kallenius G, Huard RC, Ghebremichael S, Asiimwe J, et al. Mycobacterium tuberculosis Uganda genotype is the predominant cause of TB in Kampala, Uganda. Int J Tuberc Lung Dis Off J Int Union Against Tuberc Lung Dis. 2008;12(4):386–91 (epub 2008/03/29).
Wamala D, Asiimwe B, Kigozi E, Mboowa G, Joloba M, Kallenius G. Clinico-pathological features of tuberculosis due to Mycobacterium tuberculosis Uganda genotype in patients with tuberculous lymphadenitis: a cross sectional study. BMC Clin Pathol. 2014;14(1):1472–6890.
Huard RC, Lazzarini LC, Butler WR, van Soolingen D, Ho JL. PCR-based method to differentiate the subspecies of the Mycobacterium tuberculosis complex on the basis of genomic deletions. J Clin Microbiol. 2003;41(4):1637–50.
Springer B, Stockman L, Teschner K, Roberts GD, Bottger EC. Two-laboratory collaborative study on identification of mycobacteria: molecular versus phenotypic methods. J Clin Microbiol. 1996;34(2):296–303.
Turenne CY, Tschetter L, Wolfe J, Kabani A. Necessity of quality-controlled 16S rRNA gene sequence databases: identifying nontuberculous Mycobacterium species. J Clin Microbiol. 2001;39(10):3637–48.
Muhumuza J, Asiimwe B, Kayes S, Mugyenyi R, Whalen C, Mugerwa R, et al. Introduction of an in-house PCR for routine identification of M. tuberculosis in a low-income country. Int J Tuberc Lung Dis. 2006;10(11):1262–7.
Dale JW, Brittain D, Cataldi AA, Cousins D, Crawford JT, Driscoll J, et al. Spacer oligonucleotide typing of bacteria of the Mycobacterium tuberculosis complex: recommendations for standardised nomenclature. Int J Tuberc Lung Dis. 2001;5(3):216–9.
Demay C, Liens B, Burguière T, Hill V, Couvin D, Millet J, et al. SITVITWEB—a publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology. Infect Genet Evol. 2012;12(4):755–66.
Bazira J, Asiimwe B, Joloba M, Bwanga F, Matee M. Mycobacterium tuberculosis spoligotypes and drug susceptibility pattern of isolates from tuberculosis patients in South-Western Uganda. BMC Infect Dis. 2011;11(1):81.
Warren R, Streicher E, Sampson S, van der Spuy G, Richardson M, Nguyen D, et al. Microevolution of the direct repeat region of Mycobacterium tuberculosis: implications for interpretation of spoligotyping data. J Clin Microbiol. 2002;40(12):4457–65.
Sola C, Filliol I, Legrand E, Mokrousov I, Rastogi N. Mycobacterium tuberculosis phylogeny reconstruction based on combined numerical analysis with IS1081, IS6110, VNTR, and DR-based spoligotyping suggests the existence of two new phylogeographical clades. J Mol Evol. 2001;53(6):680–9.
Tang C, Reyes JF, Luciani F, Francis AR, Tanaka MM. spolTools: online utilities for analyzing spoligotypes of the Mycobacterium tuberculosis complex. Bioinformatics. 2008;24(20):2414–5 (epub 2008/08/20).
Reyes JF, Francis AR, Tanaka MM. Models of deletion for visualizing bacterial variation: an application to tuberculosis spoligotypes. BMC Bioinform. 2008;9:496 (epub 2008/11/28).
Asiimwe BB, Ghebremichael S, Kallenius G, Koivula T, Joloba ML. Mycobacterium tuberculosis spoligotypes and drug susceptibility pattern of isolates from tuberculosis patients in peri-urban Kampala, Uganda. BMC Infect Dis. 2008;8:101 (epub 2008/07/30).
Ouassa T, Borroni E, Loukou GY, Faye-Kette H, Kouakou J, Menan H, Cirillo DM. High prevalence of shared international type 53 among Mycobacterium tuberculosis complex strains in retreated patients from Côte d’Ivoire. PLoS ONE. 2012;7(9):e45363. doi:10.1371/journal.pone.0045363.
Stavrum RMM, Øvreås K, et al. High diversity of Mycobacterium tuberculosis genotypes in South Africa and preponderance of mixed infections among ST53 isolates. J Clin Microbiol. 2009;47(6):1848–56.
Muwonge A, Malama S, Johansen TB, Kankya C, Biffa D, Ssengooba W, et al. Molecular epidemiology, drug susceptibility and economic aspects of tuberculosis in Mubende district, Uganda. PLoS ONE. 2013;8(5):e64745.
Bazira J, Matte M, Asiimwe B, Joloba L. Genetic diversity of Mycobacterium tuberculosis in Mbarara, South Western Uganda. Afr Health Sci. 2010;10(4):306–11.
Wampande EM, Mupere E, Debanne SM, Asiimwe BB, Nsereko M, Mayanja H, et al. Long-term dominance of Mycobacterium tuberculosis Uganda family in peri-urban Kampala-Uganda is not associated with cavitary disease. BMC Infect Dis. 2013;13(1):484.
Rindi L, Lari N, Garzelli C. Large sequence polymorphisms of the Euro-American lineage of Mycobacterium tuberculosis: a phylogenetic reconstruction and evidence for convergent evolution in the DR locus. Infect Genet Evol. 2012;12(7):1551–7.
Valcheva V, Mokrousov I, Panaiotov S, Bachiiska E, Zozio T, Sola C, et al. Bulgarian specificity and controversial phylogeography of Mycobacterium tuberculosis spoligotype ST 125__BGR. FEMS Immunol Med Microbiol. 2010;59(1):90–9.
Gibson AL, Huard RC, Gey van Pittius NC, Lazzarini LC, Driscoll J, Kurepina N, et al. Application of sensitive and specific molecular methods to uncover global dissemination of the major RDRio sublineage of the Latin American-Mediterranean Mycobacterium tuberculosis spoligotype family. J Clin Microbiol. 2008;46(4):1259–67.
Gafirita J, Umubyeyi AN, Asiimwe BB. A first insight into the genotypic diversity of Mycobacterium tuberculosis from Rwanda. BMC Clin Pathol. 2012;12(20):1472–6890.
Garedew L, Mihret A, Abebe T, Ameni G. Molecular typing of mycobacteria isolated from extrapulmonary tuberculosis patients at Debre Birhan Referral Hospital, central Ethiopia. Scand J Infect Dis. 2013;45(7):512–8.
Lari N, Rindi L, Cristofani R, Rastogi N, Tortoli E, Garzelli C. Association of Mycobacterium tuberculosis complex isolates of BOVIS and Central Asian (CAS) genotypic lineages with extrapulmonary disease. Clin Microbiol Infect. 2009;15(6):538–43.
Sankar MM, Singh J, Diana SC, Singh S. Molecular characterization of Mycobacterium tuberculosis isolates from North Indian patients with extrapulmonary tuberculosis. Tuberculosis. 2013;93(1):75–83.
Hirsh AE, Tsolaki AG, DeRiemer K, Feldman MW, Small PM. Stable association between strains of Mycobacterium tuberculosis and their human host populations. Proc Natl Acad Sci USA. 2004;101(14):4871–6.
Behr MA, Waters WR. Is tuberculosis a lymphatic disease with a pulmonary portal? Lancet Infect Dis. 2013;14(3):250–5.
Coscolla M, Gagneux S. Does M. tuberculosis genomic diversity explain disease diversity? Drug Discov Today Dis Mech. 2010;7(1):e43–59.
Garcia de Viedma D, Lorenzo G, Cardona PJ, Rodriguez NA, Gordillo S, Serrano MJ, et al. Association between the infectivity of Mycobacterium tuberculosis strains and their efficiency for extrarespiratory infection. J Infect Dis. 2005;192(12):2059–65.
Caws M, Thwaites G, Dunstan S, Hawn TR, Thi Ngoc Lan N, Thuong NTT, et al. The influence of host and bacterial genotype on the development of disseminated disease with Mycobacterium tuberculosis. PLoS Pathog. 2008;4(3):e1000034.
Click E, Moonan P, Winston C, Cowan L, Oeltmann J. Relationship between Mycobacterium tuberculosis phylogenetic lineage and clinical site of tuberculosis. Clin Infect Dis. 2012;54(2):211–9.
Ameni G, Tadesse K, Hailu E, Deresse Y, Medhin G, Aseffa A, et al. Transmission of Mycobacterium tuberculosis between farmers and cattle in Central Ethiopia. PLoS ONE. 2013;8(10):e76891.
Kidane D, Olobo JO, Habte A, Negesse Y, Aseffa A, Abate G, et al. Identification of the causative organism of tuberculous lymphadenitis in Ethiopia by PCR. J Clin Microbiol. 2002;40(11):4230–4.
Cleaveland S, Shaw DJ, Mfinanga SG, Shirima G, Kazwala RR, Eblate E, et al. Mycobacterium bovis in rural Tanzania: risk factors for infection in human and cattle populations. Tuberculosis. 2007;87(1):30–43.
DW: Participated in the research conceptual development, drafting the manuscript, data acquisition and interpretation and intellectual content. EK: Participated in molecular biology studies and contributed to interpretation of data. MO: Participated in molecular biology studies and made contribution to interpretation of data. DC: Analyzed the data. MJ: Participated in the research conceptual development and design and contributed to the intellectual content of the manuscript. NR: Contributed to the intellectual content; GK: Participated in conceptual development, drafting and critically revising the manuscript. Contributed substantially to the intellectual content and gave the final approval of the manuscript. All authors read and approved the final manuscript.
The study was supported by funds from the Swedish International Development Cooperation Agency through Makerere University-Karolinska Institute Research collaboration, and the Swedish Heart-Lung Foundation. We thank Prof. Colleen Wright and Mr. Ian McDaniels who offered help in copy editing the manuscript, and Associate Professor Andrew Field for a final edit. We are grateful to Dr. Samuel Kyobe of Makerere University college of Health Sciences Molecular Laboratory and Ms. Carol Namaganda of Makerere College of Health Sciences Mycobacteriology Laboratory for their help in the laboratory work. Thanks go to Mr. Sam Nyole who helped in the recruitment of patients.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Table S1. A summary of genotyping, demographic and epidemiologic data on 121 M. tuberculosis strains isolated from Ugandan patient with tuberculous lymphadenitis.
Table S2. Correlation of the major spoligotypes with HIV status.
Table S3a. Correlation of the various lineages with sex. Table S3b. Correlation between the predominant spoligotypes and age.