BioMed Research International

BioMed Research International / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 238532 |

Mulugeta Belay, Gobena Ameni, Gunnar Bjune, David Couvin, Nalin Rastogi, Fekadu Abebe, "Strain Diversity of Mycobacterium tuberculosis Isolates from Pulmonary Tuberculosis Patients in Afar Pastoral Region of Ethiopia", BioMed Research International, vol. 2014, Article ID 238532, 12 pages, 2014.

Strain Diversity of Mycobacterium tuberculosis Isolates from Pulmonary Tuberculosis Patients in Afar Pastoral Region of Ethiopia

Academic Editor: Tomasz Jagielski
Received24 Nov 2013
Accepted22 Jan 2014
Published06 Mar 2014


Data on genotypic diversity of Mycobacterium tuberculosis complex (MTBC) is important to understand its epidemiology, human adaptation, clinical phenotypes, and drug resistance. This study aimed to characterize MTBC clinical isolates circulating in a predominantly pastoralist area in Ethiopia, a country where tuberculosis is the second leading cause of mortality. Culture of sputum samples collected from a total of 325 pulmonary TB suspects was done to isolate MTBC. Spoligotyping was used to characterize 105 isolates from culture positive slopes and the result was compared with an international database. Forty-four spoligotype patterns were observed to correspond to 35 shared-types (SITs) containing 96 isolates and 9 orphan patterns; 27 SITs containing 83 isolates matched a preexisting shared-type in the database, whereas 8 SITs ( isolates) were newly created. A total of 19 SITs containing 80 isolates were clustered within this study (overall clustering of 76.19%). Three dominant lineages (T, CAS, and Manu) accounted for 76.19% of the isolates. SIT149/T3-ETH was one of the two most dominant sublineages. Unlike previous reports, we show that Manu lineage strains not only constitute a dominant lineage, but are also associated with HIV infection in Afar region of Ethiopia. The high level of clustering suggests the presence of recent transmission that should be further studied using additional genotyping markers.

1. Introduction

Tuberculosis (TB) is among the 10 major causes of death [1] and it is only outranked by HIV/AIDS among infectious causes of death worldwide [2]. Despite the global effort towards controlling TB for the last 20 years, the burden remains high with over 8.7 million incident cases and 1.4 million deaths in 2011 [3]. The HIV epidemics and the spread of MDR-TB are among the major obstacles for the successful control of TB and, in this regard, Africa is lagging behind in achieving to halve the 1990s mortality by 2015 [3].

Ethiopia is among the 22 high TB as well as the 27 high MDR-TB burden countries [3]. According to the Ministry of Health, TB is the second leading cause of death [4]. Apart from other factors, the HIV epidemics and the emergence of MDR-TB have contributed to the high TB burden in the country. According to the recent national TB prevalence survey, the prevalence of smear positive pulmonary TB is highest among pastoral communities [5].

Genotypic variations among M. tuberculosis strains as well as the existence of human genetic polymorphism linked to TB have resulted in the changing relationship between M. tuberculosis and the human host. Such changes have complicated TB control efforts [6]. Data on strain diversity of mycobacterial isolates is important to understand the transmission dynamics and phylogeographical distribution of dominant circulating strains of M. tuberculosis.

This study was part of a major project on molecular epidemiology, clinical epidemiology, and immunology of TB in a pastoral community and their livestock in Ethiopia. Some studies have reported the dominant strains among human isolates in central highlands of Ethiopia [711]; besides, strains isolated from goats and camels have been reported from Afar region recently [12, 13]. This study reports strains of M. tuberculosis isolates circulating in the Afar pastoral communities.

2. Methods

2.1. Study Area

The study area has been described elsewhere [14]. Briefly, the study was conducted in two public (Awash Health Centre and Dubti Hospital) health facilities in Afar region and three private (Selam hospital, Bati Hospital, and Amir Higher Clinic) health facilities in Dessie town. The reason for including these private health facilities is because they provide diagnostic services for a substantial number of TB patients coming from Afar region.

2.2. Study Design and Study Participants

A health facility-based cross-sectional study was conducted between September 2009 and March 2010. A total of 325 pulmonary TB suspects (≥18 years of age) who were residents of Afar region during the study period and came to the selected health facilities with cough lasting more than 2 weeks were included consecutively. TB patients who were already taking anti-TB drugs were excluded.

2.3. Data Collection

A semi-structured, pretested questionnaire was used to collect data on basic sociodemographic characteristics as well as clinically relevant symptoms. Three sputum samples (spot-morning-spot) were collected and smear microscopic examination for AFB was done at the respective health facilities according to the national TB and Leprosy control programme guideline [15]. The rest of the sputa were stored at +4°C and transported to Aklilu Lemma Institute of Pathobiology (ALIPB), Addis Ababa University, within 1 week. Participants were tested for HIV according to the national guideline [16].

2.4. Mycobacterial Culture

The three sputa from each participant were pooled at ALIPB laboratory; culture was done according to the WHO guideline [17]. Briefly, sputum samples were homogenized and decontaminated with equal volume of 4% NaOH and shaken for 15 minutes at room temperature. Subsequently, it was centrifuged at 3000 rpm for 15 minutes. After the supernatant was poured off, the sediment was neutralized with 2 N HCL and 2–4 loopfuls of the centrifuged sediment were inoculated into four slopes of LJ medium. Inspection of media was done every week for growth until 8 weeks. Based on colony morphology and smear microscopy, those with growth were identified and two colonies were transferred to eppendorf tube with 300 μL distilled water. Mycobacterial genomic DNA was extracted by heating isolates at 80°C for 60 minutes. Heat-killed isolates were stored at −20°C until spoligotyping was done.

2.5. Spoligotyping and Database Comparison

Spoligotyping of 105 isolates was performed using a commercially available kit from Ocimum Biosolutions Company, Iisselstein, The Netherlands, according to the company’s instructions and as described previously [9, 18]. Briefly, the direct repeat (DR) region was amplified by a Thermal Cycler PCR machine (VWR International) using oligonucleotides and primers derived from this region. The amplified product was hybridized to a set of 43 immobilized oligonucleotides, each corresponding to one of the unique spacer. Hybridized DNA was detected by chemiluminescence method (Amersham Biosciences, Little Chalfont, UK) and by exposure to X-ray film (Hyperfilm ECL, Amersham Biosciences), as specified by the manufacturer. The hybridization patterns were converted into binary and octal formats and compared with previously reported strains in the in-house SITVIT2 proprietary database of Institut Pasteur de la Guadeloupe, which is an updated version of the recently released SITVITWEB database [19]. In this database, a Spoligotype International Type (SIT) designates a pattern shared by 2 or more patient isolates, whereas “orphan” represents a pattern reported for a single isolate. At the time of this comparison, SITVIT2 contained genotyping data on more than 110,000 MTBC clinical isolates from 160 countries of patient origin, and >3500 SITs. Major phylogenetic clades were essentially assigned according to signatures provided earlier [19]; these included specific signatures for various MTBC members as well as rules defining major lineages/sublineages for M. tuberculosis sensu stricto. The latter included the Beijing clade, the Central-Asian (CAS) clade and its 2 sublineages, the East-African-Indian (EAI) clade with its 9 sublineages, the Haarlem clade and its 3 sublineages, the Latin-American-Mediterranean (LAM) clade and its 10 sublineages, the “Manu” family and its 3 sublineages, the IS6110-low banding X clade and its 3 sublineages, and the ill-defined T clade and its 5 sublineages. The high phylogeographical specificity of LAM10-CAM prototype SIT61 for Cameroon has led to its designation as the Cameroon lineage [20], whereas LAM7-TUR lineage was tentatively reclassified as the Turkey lineage [21].

We also studied the worldwide distribution of all major spoligotyping clusters in our study (clusters containing 3 or more isolates) by interrogating the database for their distribution in macrogeographical regions and subregions according to the United Nations (; for further details please refer to the footnote of Table 3). Lastly, the distribution of predominant spoligotype patterns in this study ( strains) was further compared to all other strains reported from Ethiopia in the international database ( strains excluding the present study) as well as the five African subregions ( strains), and the neighboring Western Asian countries ( strains; for further details please refer to the footnote of Table 4).

Minimum spanning trees (MSTs) illustrating evolutionary relationships between the M. tuberculosis spoligotypes and various sociodemographic variables were constructed using BioNumerics software version 6.6 (Applied Maths, Sint-Martens-Latem, Belgium). MST is an undirected network in which all of the samples are linked together with the fewest possible linkages between nearest neighbors. Lastly, Spoligoforest trees which represent another way to illustrate probable strain evolutionary relationships between spoligotypes were drawn [22, 23] using spolTools online utilities available through As opposed to the MSTs, the spoligoforest trees are directed and only evolve by loss of spacers. GraphViz software ( [24]) was used to color the strains based on their lineages on the spoligoforest trees.

2.6. Ethics Statement

The study has been ethically cleared by the Norwegian Ethics Committee and the Ethiopian National Ethics Committee. Written informed consent was obtained from each study participant.

2.7. Statistical Analysis

Chi-square test using STATA software version 12 was used to evaluate whether or not significant association may exist between lineages or SITs with sociodemographic and epidemiological characteristics. Pearson’s Chi-square test was used when more than 80% of data had an expected value greater than 5 and Fisher’s Exact Test for remaining data with smaller expected values (at least 20% of data having values less than 5). values less than 0.05 were considered as statistically significant.

3. Results

3.1. Demographic Characteristics of Study Participants

Among 325 pulmonary TB suspects, a total of 105 patients were confirmed to have culture positive pulmonary TB and included for further analysis. The median age of patients was 29 (IQR: 22.5–40.0) years and the majority were males with a male to female ratio of 1.92. Thirty-nine (37%) of culture confirmed TB patients were positive for acid fast bacilli on sputum smear examination. Besides, HIV result was available for 95 of the patients and 39% of them were found to be HIV co-infected. The majority (77.1%) of the study participants were new cases of pulmonary TB. Pastoralists represented 42.9% of the study participants. The proportion of HIV positives among ethnic Afar (8.7%) was significantly lower compared to other ethnic groups (28.0–47.5%) ().

3.2. Genetic Diversity and Distribution of Isolates

Overall, 103 (98.1%) isolates had spoligotypes characteristic of M. tuberculosis complex and 2 (1.9%) isolates had spoligotypes characteristic of M. bovis. A total of 9 isolates showed orphan patterns (Table 1), whereas 16 strains were unique making the total number of unclustered isolates to be 25 (23.81%). On the other hand, 19 SITs containing 80 isolates (76.19%) were clustered within this study (2 to 10 isolates per cluster). Besides, a total of 44 spoligotype patterns were observed and, therefore, the overall diversity of the isolates was 41.9%. A total of 27 SITs containing 83 isolates matched a preexisting shared-type in the SITVIT2 database, whereas 8 SITs ( isolates) did not match any of the previous isolates in the database and hence were newly created (Table 2).

Iso numberYearStrain Spoligotype description
Octal codeLineageSIT


Lineage designations for orphan patterns were done manually as expert-based interpretations using revised SpolDB4 rules.

SIT*Spoligotype descriptionOctal numberNb (%) in study% in study versus databaseLineage**Clustered versus unique patterns***

10238532.table.002a14777772774137711 (0.95)1.14EAI8-MDGUnique
20238532.table.002a26777776077607712 (1.9)0.24LAM1Clustered
25238532.table.002a37037777400031714 (3.81)0.64CAS1-DelhiClustered
26238532.table.002a47037777400037714 (3.81)0.26CAS1-DelhiClustered
34238532.table.002a57763777777607711 (0.95)0.12SUnique
37238532.table.002a677773777776077110 (9.52)2.04T3Clustered
41238532.table.002a77777774047607711 (0.95)0.21LAM7-TURUnique
50238532.table.002a87777777777207711 (0.95)0.03H3Unique
52238532.table.002a97777777777607316 (5.71)0.65T2Clustered
53238532.table.002a107777777777607719 (8.57)0.14T1Clustered
54238532.table.002a117777777777637719 (8.57)3.56Manu2Clustered
100238532.table.002a127777777777737711 (0.95)1.22Manu1Unique
119238532.table.002a137777767777607713 (2.86)0.28X1Clustered
121238532.table.002a147777777757207713 (2.86)5.77H3Clustered
149238532.table.002a1577700037776077110 (9.52)2.07T3-ETHClustered
241238532.table.002a167777777777604111 (0.95)3.7T1Unique
247238532.table.002a177037777400034712 (1.9)5.13CAS1-DelhiClustered
289238532.table.002a187037777400035712 (1.9)15.38CAS1-DelhiClustered
336238532.table.002a197777767777607312 (1.9)1.71X1Clustered
523238532.table.002a207777777777777713 (2.86)6.25Manu_ancestorClustered
804238532.table.002a214777777777607711 (0.95)3.85T1Unique
1264238532.table.002a227037777400000001 (0.95)4.35CAS1-DelhiUnique
1338238532.table.002a237777777777776711 (0.95)16.67UnknownUnique
1967238532.table.002a247037777400000312 (1.9)28.57CASClustered
2693238532.table.002a257037637400037711 (0.95)11.11CAS1-DelhiUnique
2781238532.table.002a267777777777677711 (0.95)11.11UnknownUnique
2820238532.table.002a277777777677607311 (0.95)14.29T2Unique
3896*238532.table.002a287777377777600001 (0.95)50T3Unique
3897*238532.table.002a297773477777636713 (2.86)75Manu2Clustered
3898*238532.table.002a306777770774137711 (0.95)50EAI2Unique
3899*238532.table.002a317037777000017711 (0.95)50CASUnique
3900*238532.table.002a327773477777637712 (1.9)66.67Manu2Clustered
3901*238532.table.002a330007737777776002 (1.9)66.67BOVClustered
3902*238532.table.002a347777774574137712 (1.9)66.67EAI6-BGD1Clustered
3903*238532.table.002a356777774077607711 (0.95)50LAM1Unique

A total of 27/35 SITs containing 83 isolates matched a preexisting shared-type in the database, whereas 8/35 SITs ( isolates) were newly created. A total of 19/35 SITs containing 80 isolates were clustered within this study (2 to 10 isolates per cluster) while 16/35 SITs containing 16 strains were unique (for total unique strains, one should add to this number the 9 orphan strains which brings the number of unclustered isolates in this study to 25/105 (23.81%) and clustered isolates to 80/105 (76.19%)). Note that SITs followed by an asterisk indicates “newly created” SITs due to 2 or more strains belonging to an identical new pattern within this study or after a match with an orphan in the database; SIT designations followed by number of strains: 3896* this study . DEU ; 3897* this study . USA ; 3898* this study . AUS ; 3899* this study . NPL ; 3900* this study . CHN ; 3901* this study . FXX ; 3902* this study . MYS ; 3903* this study . BRA .
**Lineage designations according to SITVIT2 using revised SpolDB4 rules; “unknown” designates patterns with signatures that do not belong to any of the major lineages described in the database.
***Clustered strains correspond to a similar spoligotype pattern shared by 2 or more strains “within this study”; as opposed to unique strains harboring a spoligotype pattern that does not match with another strain from this study. Unique strains matching a preexisting pattern in the SITVIT2 database are classified as SITs, whereas in case of no match, they are designated as “orphan” (see Table 1).

SIT (Lineage) Octal number spoligotype descriptionNb (%) in study% in study versus databaseDistribution in regions with ≥3% of a given SITs*Distribution in countries with ≥3% of a given SITs**

238532.table.003a10 (9.52)2.04AFRI-E 20.37, EURO-N 12.63, EURO-W 12.22, ASIA-W 11.2, AMER-N 9.16, ASIA-E 6.52, EURO-S 6.11, EURO-E 4.68, ASIA-S 3.87, AMER-S 3.67, AFRI-S 3.46ETH 17.72, USA 8.15, SWE 5.5, SAU 5.3, CHN 5.3, FXX 4.48, ITA 4.28, ZAF 3.46, DNK 3.46, FIN 3.06

238532.table.003b10 (9.52)2.07AFRI-E 58.06, EURO-N 21.69, ASIA-W 7.03, EURO-W 6.41, AMER-N 4.75ETH 58.06, DNK 11.78, SAU 6.61, USA 4.75, SWE 4.55, NLD 3.51

238532.table.003c9 (8.57)0.14EURO-W 15.71, AMER-N 13.54, AMER-S 12.08, EURO-S 9.45, EURO-N 7.51, ASIA-W 7.34, AFRI-S 4.99, AFRI-E 4.67, ASIA-E 4.25, AFRI-N 3.53, EURO-E 3.28, AMER-C 3.25USA 13.25, FXX 7.91, ITA 5.36, BRA 5.15, ZAF 4.88, TUR 3.49, AUT 3.44, CHN 3.1

238532.table.003d9 (8.57)3.56ASIA-W 16.6, ASIA-S 14.63, ASIA-E 14.63, AFRI-N 14.23, AMER-N 7.91, AFRI-S 6.32, AFRI-E 5.14, ASIA-N 3.95, AMER-S 3.56CHN 14.63, IND 13.83, EGY 13.44, USA 7.91, SAU 7.91, ZAF 6.32, IRQ 5.53, RUS 3.95, ETH 3.56

238532.table.003e6 (5.71)0.65EURO-W 19.63, ASIA-E 15.29, EURO-N 13.88, AMER-N 12.04, EURO-S 5.32, ASIA-W 5.32, AFRI-E 4.99, EURO-E 4.34, AFRI-M 3.58, AMER-C 3.47CHN 11.82, USA 11.61, SWE 9.76, FXX 9.0, BEL 4.88, ITA 3.47, JPN 3.25, ETH 3.25

238532.table.003f4 (3.81)0.64ASIA-W 27.05, AFRI-E 24.96, AFRI-N 13.53, ASIA-S 12.4, AMER-N 8.37, EURO-N 5.8, EURO-W 4.35ETH 23.83, SAU 15.14, SDN 11.11, USA 8.37, IND 7.89, OMN 6.76, IRQ 4.67

238532.table.003g4 (3.81)0.26ASIA-S 50.0, AMER-N 14.47, ASIA-W 14.07, EURO-W 5.61, AFRI-E 4.95, EURO-N 4.36, EURO-S 3.3IND 30.19, USA 14.47, PAK 9.51, OMN 7.2, BGD 5.48, SAU 5.15, IRN 4.1, ITA 3.17, ETH 3.1

238532.table.003h3 (2.86)0.28AMER-N 58.87, AFRI-S 13.22, AMER-C 9.89, AMER-S 3.97, EURO-N 3.51USA 58.23, ZAF 13.22, MEX 9.61

238532.table.003i3 (2.86)5.77ASIA-E 38.46, AFRI-E 19.23, EURO-W 11.54, EURO-E 7.69, EURO-N 5.77, ASIA-W 5.77, AMER-N 5.77, EURO-S 3.85TWN 36.54, ETH 19.23, FXX 11.54, CZE 7.69, USA 5.77, SAU 5.77

238532.table.003j3 (2.86)6.25ASIA-E 27.08, ASIA-SE 16.67, AMER-N 10.42, EURO-W 8.33, ASIA-W 8.33, AFRI-W 8.33, AFRI-E 6.25, AFRI-S 4.17USA 10.42, SAU 8.33, JPN 8.33, CHN 8.33, NGA 6.25, MYS 6.25, KOR 6.25, ETH 6.25, ZAF 4.17, VNM 4.17, TWN 4.17, THA 4.17, BEL 4.17

238532.table.003k3 (2.86)75AFRI-E 75.0, AMER-N 25.0ETH 75.0, USA 25.0

Worldwide distribution (analysis made on August 18, 2013) is reported for regions with more than 3% of given SITs as compared to their total number in the SITVIT2 database. The definition of macrogeographical regions and subregions ( is according to the United Nations; regions: AFRI (Africa), AMER (Americas), ASIA (Asia), EURO (Europe), and OCE (Oceania), subdivided in: E (Eastern), M (Middle), C (Central), N (Northern), S (Southern), SE (South-Eastern), and W (Western). Furthermore, CARIB (Caribbean) belongs to Americas, while Oceania is subdivided in 4 subregions, AUST (Australasia), MEL (Melanesia), MIC (Micronesia), and POLY (Polynesia). Note that in our classification scheme, Russia has been attributed a new subregion by itself (Northern Asia) instead of including it among the rest of Eastern Europe. It reflects its geographical localization as well as due to the similarity of specific TB genotypes circulating in Russia (a majority of Beijing genotypes) with those prevalent in Central, Eastern, and Southeastern Asia.
**The 3-letter country codes are according to; countrywide distribution is only shown for SITs with ≥3% of given SITs as compared to their total number in the SITVIT2 database.

Octal code (SIT) lineageEthiopia valueDistribution in five African subregions and neighboring Western Asia in SITVIT2b
This study SITVIT2 AFRI-E -valueAFRI-N valueAFRI-W valueAFRI-M valueAFRI-S valueASIA-W value
(%) (%) (%) (%) (%) (%) (%) (%)

37 (T3) 10/105 77/1507 0.07290/4138 0.0003*7/2145 <0.0001*6/1985 <0.0001*0/808 <0.0001*17/4580 <0.0001*55/4790 <0.0001*
149 (T3-ETH) 10/105 271/1507 0.056271/4138 0.261/2145 <0.0001*0/1985 <0.0001*0/808 <0.0001*1/4580 <0.0001*34/4790 <0.0001*
53 (T1) 9/105 154/1507 0.62283/4138 0.52221/2145 0.60109/1985 0.2124/808 0.006*312/4580 0.51459/4790 0.75
777777777760771 (8.57)(10.22)(6.84)(10.30)(5.49)(2.97)(6.81)(9.58)
54 (Manu2) 9/105 0/1507 <0.0001*4/4138 <0.0001*36/2145 <0.0001*3/1985 <0.0001*0/808 <0.0001*16/4580 <0.0001*42/4790 <0.0001*
52 (T2) 6/105 24/1507 0.0036*40/4138 0.001*6/2145 <0.0001*6/1985 <0.0001*33/808 0.4421/4580 <0.0001*49/4790 0.0013*
25 (CAS1-Delhi) 4/105 144/1507 0.066151/4138 0.7984/2145 1.02/1985 <0.0001*1/808 0.0009*3/4580 <0.0001*168/4790 0.79
26 (CAS1-Delhi) 4/105 43/1507 0.5571/4138 0.1216/2145 0.013*0/1985 <0.0001*0/808 0.0002*6/4580 <0.0001*213/4790 1.0
119 (X1) 3/105 1/1507 0.001*4/4138 0.0005*5/2145 0.005*2/1985 0.0012*0/808 0.0016*143/4580 1.09/4790 0.002*
121 (H3) 3/105 7/15070.024*7/4138 0.002*0/2145 0.0001*0/1985 0.0001*0/808 0.0016*0/4580 <0.0001*3/4790 0.0002*
777777775720771(2.86) (0.46)(0.17)(0.00)(0.00)(0.00)(0.00)(0.06)
523 (Manu_ancestor) 3/105 0/1507 0.0003*0/4138 <0.0001*1/2145 0.0004*4/1985 0.004*0/808 0.0016*2/4580 0.0001*4/4790 0.0003*
3897 (Manu2) 3/105 0/15070.0003*0/4138 <0.0001*0/2145 0.0001*0/1985 0.0001*0/808 0.0016*0/4580 <0.0001*0/4790 <0.0001*
777347777763671(2.86) (0.00)(0.00)(0.00)(0.00)(0.00)(0.00)(0.00)

Asterisk denotes statistically significant differences as compared to the total number of strains for a given pattern in this study ( strains). For statistical comparison of differences observed, note that Pearson’s Chi-square test was used when more than 80% of data had an expected value greater than 5, and Fisher’s Exact Test for remaining data with smaller values (at least 20% of data having values less than 5).
bThe definition of subregions is according to the United Nations ( The table summarizes distribution in this study versus all strains reported from Ethiopia in the international database ( , excluding the present study), the five African subregions , and neighboring Western Asian countries (ASIA-W, ). AFRI-E: Burundi, Comoros, Ethiopia, Kenya, Madagascar, Malawi, Mauritius, Mozambique, Reunion, Tanzania, Uganda, Zambia, and Zimbabwe; AFRI-N: Algeria, Egypt, Libyan Arab Jamahiriya, Morocco, Sudan, and Tunisia; AFRI-W: Benin, Burkina Faso, Gambia, Guinea, Guinea-Bissau, Ivory Coast, Mali, Mauritania, Nigeria, Senegal, and Sierra Leone; AFRI-M: Cameroon, Central African Republic, Chad, and Democratic Republic of Congo; AFRI-S: Botswana, Namibia, South Africa, and Swaziland; ASIA-W: Armenia, Azerbaijan, Georgia, Iraq, Israel and Palestinian territories, Oman, Saudi Arabia, Turkey, and Yemen (SITVIT2 database comparison made on August 21st 2013 by David Couvin and Nalin Rastogi).

The majority (76.19%) of the isolates belonged to the following three dominant lineages: T (, 40.95%), Manu (, 19.05%), and CAS (, 16.19%). The other minor lineages in this study include X (, 5.71%), Haarlem (, 3.81%), EAI (, 3.81%), LAM (, 2.86%), BOV (, 1.90%), AFRI (, 0.95%), Turkey (, 0.95%), and S (, 0.95%) (Table 2). Three isolates (2.86%) had patterns with signatures that do not belong to any of the major lineages described in the database and hence labelled as “unknown” (Tables 1 and 2). No Beijing strains were identified in this study.

The dominant SITs in this study included SIT 149 of the T3-ETH lineage, SIT 37 of the T3 lineage, SIT 53 of the T1 lineage, and SIT 54 of the Manu2 lineage all together accounting for 36.2% of the isolates (Table 3). On comparing the dominant isolates from this study with the isolates in the SITVIT2 database from Ethiopia, the five African regions, and neighboring Western Asia, statistically significant differences were observed (Table 4). Significantly higher proportions of T2 (), X1 () and H () lineages/sublineages were found among isolates in this study compared to their representations in the SITVIT2 database from Ethiopia. Notably, the Manu lineage was significantly higher compared to previous isolates from Africa and neighboring Western Asia (). Besides, this lineage has not been reported to SITVIT2 database from Ethiopia so far. On the other hand, there was no significant difference in the distribution of T1, T3, T3-ETH, and CAS1-Delhi lineages compared to previous isolates from Ethiopia ().

A minimum spanning tree (MST) illustrating evolutionary relationships among the strains has been constructed (Figure 1). In this evolutionary tree, MANU strains were closest to the central node of the unrooted tree which is represented by SIT53, whereas CAS strains are further away from the central node of the tree. On spoligoforest trees (Figure 2), SIT149/T3-ETH and SIT37/T3 are the biggest nodes (each with 10 isolates), followed by SIT53/T1 and SIT54/Manu2 (each with 9 isolates); SIT52/T2 (), SIT25, and SIT26/CAS1-Delhi (each with 4 isolates) are the other dominant isolates in our study.

On investigating the association between strains and HIV co-infection, strains belonging to Manu2 lineage were significantly associated with HIV infection compared to strains belonging to CAS lineage (Fisher’s Exact Test, ). However, strain clustering was not significantly associated with HIV infection and other sociodemographic characteristics of the study participants.

4. Discussion

Analysis of strain diversity and comparison with international database is key to have an insight into the global distribution of M. tuberculosis strains. In this regard, we describe the strain diversity of M. tuberculosis complex among clinical isolates of pulmonary TB patients from a predominantly pastoralist area in Northeast Ethiopia and compared the strains with SITVIT2 database strains.

Among 105 isolates, the majority (98.1%) of the isolates were identified as M. tuberculosis complex with only 2 (1.9%) isolates being M. bovis. This finding is in agreement with previous reports from Ethiopia [25, 26]. Unlike a previous study from Ethiopia [27] which reported a significant contribution (17%) of bovine TB among TB lymphadenitis patients, the contribution of M. bovis to human TB seems to be minimal and mainly restricted to the pastoral communities as supported by current evidence [25, 26].

In our study, the overall clustering was 76.19% indicating a high rate of recent transmission in the study area. Although this high level of spoligotyping based clustering should be further studied using additional genotyping markers, it is comparable with a previous study from Addis Ababa [8]. The recent national TB prevalence survey also argues for an ongoing TB epidemic mainly affecting the young probably indicating a defect in the TB control programme [28]. Besides, we have previously reported [29] that TB patients in Afar region suffered long diagnostic delays contributing to continued transmission with strain circulation.

Three dominant lineages were identified in this study and in agreement with previous reports [8, 9], the ill-defined T-lineage is the most dominant lineage accounting for 40.95% of the isolates. Among the T-lineages, T3-ETH is the most common lineage accounting for 25.6% of the T-lineage in this study. Similarly, previous studies from Ethiopia [8, 9, 11] reported a high proportion of T3-ETH among isolates indicating that this lineage is one of the dominant lineages in the country. CAS-Delhi lineage is the other dominant lineage in this study accounting for 16.19% of the isolates. In agreement with our finding, this lineage has been reported in previous studies as one of the dominant lineages in Ethiopia [79, 11] as well as in neighboring Djibouti [30]. EAI lineage was also isolated in smaller proportion in the study area. Both CAS and EAI are mainly prevalent in Middle East and Central Asia (Table 3). Two hypotheses could explain the presence of these lineages in Ethiopia as well: (i) it could have emerged in Ethiopia and migrated to Middle East and Central Asia, a hypothesis in agreement with the suggestion that East Africa is the cradle of M. tuberculosis complex species [31]; this is supported by a recent evidence [32] that MTBC co-evolved with the modern human host and migrated from Africa (particularly East Africa) to Asia and other parts of the world; (ii) alternatively, it could be due to migration of this lineage from Middle East and Central Asia to Ethiopia due to the recent human migrations from these areas to Ethiopia as suggested by a previous study [7].

Interestingly, Manu is one of the dominant lineages in our study contributing to 19.05% of the isolates. In agreement with this, investigators from Egypt reported a high proportion (27%) of Manu among their isolates [33]. In Ethiopia, only one study [9] from the Central highlands reported SIT 54 contributing to 8.33% of the isolates which is similar to our finding for SIT 54. However, Manu has not been reported from Ethiopia in other previous studies [7, 8, 11] as well as in the SITVIT2 database in which SIT 54 (Manu2) has been mainly reported from South and East Asia, Middle East including Egypt, and America. Furthermore, a recent study on a relatively large number of isolates mainly from Djibuti [34] did not find Manu lineage. The diverse distribution of strains in the study areas included in Ethiopia so far implies the need to map the distribution of Mycobacterium tuberculosis in the country. Interestingly, Manu2 lineage was found to be significantly associated with HIV infection although the reason for such association is not apparent from the current study.

Regarding the evolutionary relationship among strains, as opposed to CAS strains, the Manu strains in the MST tree are the ones closest to the dominant T lineage as reported before [33]. Besides, the T1 strain represented by its prototype SIT53 constitutes the central node of this unrooted tree and SIT 52 of T2 and SIT 37 of T3 are the other big nodes closest to the central node (Figure 1(a)). SIT149 (T3-ETH), which is the commonest lineage in Ethiopia, is further from the central node.

In conclusion, this study identified the presence of high level of clustering suggesting the presence of recent transmission in the study area. Nonetheless, one limitation of the present study is the fact that Mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) [35, 36] that allows splitting certain spoligotyping defined clusters in smaller subclusters [37] could not be used in the present study. Considering a recent study showing that some of the Manu lineage strains could result from mixed or polyclonal infections [38], and seeing the higher than usual proportion of Manu lineage strains in our study, it might be important to systematically perform MIRU typing of all spoligotyping clusters in future studies in the Afar pastoral region. We suggest that a nationwide study aiming to map the population structure of M. tuberculosis complex should be planned using spoligotyping, MIRU-VNTRs, and SNP in Ethiopia.

Conflict of Interests

The authors declared that they have no competing interests.

Authors’ Contribution

Mulugeta Belay was involved in the design, data collection and laboratory work, statistical analysis, and interpretation of data and drafted the manuscript; Fekadu Abebe was involved in the design, interpretation, and critical revision of the paper; Gunnar Bjune was involved in the design, interpretation, and critical revision of the paper; Gobena Ameni was involved in data collection, spoligotyping of the isolates, and critical revision of the paper. Nalin Rastogi and David Couvin analyzed the strains and compared them with SITVIT2 database; besides, both were involved in the critical revision of the paper. All authors have read and approved the paper for submission.


This study was supported by the Norwegian Programme for Development, Research and Education (NUFU) (Project no. NUFUPRO.2007/10198), and the Norwegian Research Council (Project no. 196397/S50). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the paper. The authors’ sincere gratitude goes to Mr. Adane Worku at Aklilu Lemma Institute of Pathobiology for carrying out spoligotyping. Laboratory technicians working at the respective health facilities deserve our sincere gratitude for collecting samples. The authors would like to acknowledge participants for consenting to participate in the study.


  1. D. Bleed, C. Dye, and M. C. Raviglione, “Dynamics and control of the global tuberculosis epidemic,” Current Opinion in Pulmonary Medicine, vol. 6, no. 3, pp. 174–179, 2000. View at: Google Scholar
  2. T. R. Frieden, T. R. Sterling, S. S. Munsiff, C. J. Watt, and C. Dye, “Tuberculosis,” The Lancet, vol. 362, no. 9387, pp. 887–899, 2003. View at: Publisher Site | Google Scholar
  3. WHO, Global Tuberculosis Report, WHO, Geneva, Switzerland, 2012.
  4. Federal Ministry of Health, Guidelines For Clinical and Programmatic Management of TB, Leprosy and TB/HIV in Ethiopia, Federal Ministry of Health, Addis Ababa, Ethiopia, 12th edition, 2012.
  5. Z. Alebachew, “Ethiopian National TB Prevalence survey 2010–2011: Preliminary Result,” View at: Google Scholar
  6. C. Dye and B. G. Williams, “The population dynamics and control of tuberculosis,” Science, vol. 328, no. 5980, pp. 856–861, 2010. View at: Publisher Site | Google Scholar
  7. B. Tessema, J. Beer, M. Merker et al., “Molecular epidemiology and transmission dynamics of Mycobacterium tuberculosis in Northwest Ethiopia: new phylogenetic lineages found in Northwest Ethiopia,” BMC Infectious Diseases, vol. 13, article 131, 2013. View at: Google Scholar
  8. A. Mihret, Y. Bekele, A. G. Loxton et al., “Diversity of Mycobacterium tuberculosis isolates from new pulmonary tuberculosis cases in Addis Ababa, Ethiopia,” Tuberculosis Research and Treatment, vol. 2012, Article ID 892079, 7 pages, 2012. View at: Publisher Site | Google Scholar
  9. L. Garedew, A. Mihret, G. Mamo et al., “Strain diversity of mycobacteria isolated from pulmonary tuberculosis patients at Debre Birhan Hospital, Ethiopia,” International Journal of Tuberculosis and Lung Disease, vol. 17, no. 8, pp. 1076–1081, 2013. View at: Google Scholar
  10. A. Deribew, G. Abebe, L. Apers et al., “Prevalence of pulmonary TB and spoligotype pattern of Mycobacterium tuberculosis among TB suspects in a rural community in Southwest Ethiopia,” BMC Infectious Diseases, vol. 12, article 54, 2012. View at: Publisher Site | Google Scholar
  11. M. Agonafir, E. Lemma, D. Wolde-Meskel et al., “Phenotypic and genotypic analysis of multidrug-resistant tuberculosis in Ethiopia,” International Journal of Tuberculosis and Lung Disease, vol. 14, no. 10, pp. 1259–1265, 2010. View at: Google Scholar
  12. G. Mamo, F. Abebe, Y. Worku et al., “Tuberculosis in goats and sheep in afar pastoral region of ethiopia and isolation of Mycobacterium tuberculosis from goat,” Veterinary Medicine International, vol. 2012, Article ID 869146, 8 pages, 2012. View at: Publisher Site | Google Scholar
  13. G. Mamo, G. Bayleyegn, T. S. Tessema et al., “Pathology of camel tuberculosis and molecular characterization of its causative agents in pastoral regions of Ethiopia,” PLoS ONE, vol. 6, no. 1, Article ID e15862, 2011. View at: Publisher Site | Google Scholar
  14. M. Belay, G. Bjune, G. Ameni, M. Abebe, and F. Abebe, “Serodiagnostic performance of rESAT-6-CFP-10 in the diagnosis of pulmonary tuberculosis in Ethiopia,” Mycobacterial Diseases, vol. 1, no. 2, article 103, 2011. View at: Google Scholar
  15. Federal Ministry of Health, Tuberculosis, Leprosy and TB/HIV Prevention and Control Programme, Federal Ministry of Health, Addis Ababa, Ethiopia, 4th edition, 2008.
  16. HAPCO, Guidelines for HIV Counselling and Testing in Ethiopia, HAPCO, Addis Ababa, Ethiopia, 2007.
  17. WHO, Laboratory Services in Tuberculosis Control: Culture. Part III, WHO, Geneva, Switzerland, 1998.
  18. J. Kamerbeek, L. Schouls, A. Kolk et al., “Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology,” Journal of Clinical Microbiology, vol. 35, no. 4, pp. 907–914, 1997. View at: Google Scholar
  19. C. Demay, B. Liens, T. Burguière et al., “SITVITWEB—a publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology,” Infection, Genetics and Evolution, vol. 12, no. 4, pp. 755–766, 2012. View at: Publisher Site | Google Scholar
  20. F. Koro Koro, Y. Kamdem Simo, F. F. Piam et al., “Population dynamics of tuberculous Bacilli in Cameroon as assessed by spoligotyping,” Journal of Clinical Microbiology, vol. 51, no. 1, pp. 299–302, 2013. View at: Google Scholar
  21. O. Kisa, G. Tarhan, S. Gunal et al., “Distribution of spoligotyping defined genotypic lineages among drug-resistant Mycobacterium tuberculosis complex clinical isolates in Ankara, Turkey,” PLoS ONE, vol. 7, no. 1, Article ID e30331, 2012. View at: Publisher Site | Google Scholar
  22. J. F. Reyes, A. R. Francis, and M. M. Tanaka, “Models of deletion for visualizing bacterial variation: an application to tuberculosis spoligotypes,” BMC Bioinformatics, vol. 9, article 496, 2008. View at: Publisher Site | Google Scholar
  23. C. Tang, J. F. Reyes, F. Luciani, A. R. Francis, and M. M. Tanaka, “spolTools: online utilities for analyzing spoligotypes of the Mycobacterium tuberculosis complex,” Bioinformatics, vol. 24, no. 20, pp. 2414–2415, 2008. View at: Publisher Site | Google Scholar
  24. J. G. E. Ellson, E. Koutsofios, S. C. North, and G. Woodhull, “Graphviz and Dynagraph—Static and Dynamic Graph Drawing Tools,” 2002, View at: Google Scholar
  25. R. Firdessa, S. Berg, E. Hailu et al., “Mycobacterial lineages causing pulmonary and extrapulmonary tuberculosis,” Emerging Infectious Diseases, vol. 19, no. 3, pp. 460–463, 2013. View at: Google Scholar
  26. B. Gumi, E. Schelling, S. Berg et al., “Zoonotic transmission of tuberculosis between pastoralists and their livestock in south-east Ethiopia,” EcoHealth, vol. 9, no. 2, pp. 139–149, 2012. View at: Publisher Site | Google Scholar
  27. D. Kidane, J. O. Olobo, A. Habte et al., “Identification of the causative organism of tuberculous lymphadenitis in Ethiopia by PCR,” Journal of Clinical Microbiology, vol. 40, no. 11, pp. 4230–4234, 2002. View at: Publisher Site | Google Scholar
  28. A. Kebede, Z. Alebachew, F. Tsegaye et al., “Background document 6: The First Ethiopian National Population-based TB Prevalence Survey,” View at: Google Scholar
  29. M. Belay, G. Bjune, G. Ameni, and F. Abebe, “Diagnostic and treatment delay among Tuberculosis patients in Afar Region, Ethiopia: a cross-sectional study,” BMC Public Health, vol. 12, article 369, 2012. View at: Google Scholar
  30. S. Godreuil, F. Renaud, M. Choisy et al., “Highly structured genetic diversity of the Mycobacterium tuberculosis population in Djibouti,” Clinical Microbiology and Infection, vol. 16, no. 7, pp. 1023–1026, 2010. View at: Publisher Site | Google Scholar
  31. G. M Cristina, S. Brisse, R. Brosch et al., “Ancient origin and gene mosaicism of the progenitor of Mycobacterium tuberculosis,” PLoS Pathogens, vol. 1, no. 1, article e5, 2005. View at: Publisher Site | Google Scholar
  32. I. Comas, M. Coscolla, T. Luo et al., “Out-of-Africa migration and Neolithic coexpansion of Mycobacterium tuberculosis with modern humans,” Nat Genet, vol. 45, no. 10, pp. 1176–1182, 2013. View at: Google Scholar
  33. Z. H. Helal, M. S. E.-D. Ashour, S. A. Eissa et al., “Unexpectedly high proportion of ancestral manu genotype Mycobacterium tuberculosis strains cultured from tuberculosis patients in Egypt,” Journal of Clinical Microbiology, vol. 47, no. 9, pp. 2794–2801, 2009. View at: Publisher Site | Google Scholar
  34. Y. Blouin, Y. Hauck, C. Soler et al., “Significance of the identification in the horn of Africa of an exceptionally deep branching Mycobacterium tuberculosis clade,” PLoS ONE, vol. 7, no. 12, Article ID e52841, 2012. View at: Google Scholar
  35. P. Supply, S. Lesjean, E. Savine, K. Kremer, D. van Soolingen, and C. Locht, “Automated high-throughput genotyping for study of global epidemiology of Mycobacterium tuberculosis based on mycobacterial interspersed repetitive units,” Journal of Clinical Microbiology, vol. 39, no. 10, pp. 3563–3571, 2001. View at: Publisher Site | Google Scholar
  36. P. Supply, C. Allix, S. Lesjean et al., “Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis,” Journal of Clinical Microbiology, vol. 44, no. 12, pp. 4498–4510, 2006. View at: Publisher Site | Google Scholar
  37. C. Sola, I. Filliol, E. Legrand et al., “Genotyping of the Mycobacterium tuberculosis complex using MIRUs: association with VNTR and spoligotyping for molecular epidemiology and evolutionary genetics,” Infection, Genetics and Evolution, vol. 3, no. 2, pp. 125–133, 2003. View at: Publisher Site | Google Scholar
  38. L. C. O. Lazzarini, J. Rosenfeld, R. C. Huard et al., “Mycobacterium tuberculosis spoligotypes that may derive from mixed strain infections are revealed by a novel computational approach,” Infection, Genetics and Evolution, vol. 12, no. 4, pp. 798–806, 2012. View at: Publisher Site | Google Scholar

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

More related articles

 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles