Research Article | Open Access
Anwar Nuru, Gezahegne Mamo, Adane Worku, Aschalew Admasu, Girmay Medhin, Rembert Pieper, Gobena Ameni, "Genetic Diversity of Mycobacterium tuberculosis Complex Isolated from Tuberculosis Patients in Bahir Dar City and Its Surroundings, Northwest Ethiopia", BioMed Research International, vol. 2015, Article ID 174732, 9 pages, 2015. https://doi.org/10.1155/2015/174732
Genetic Diversity of Mycobacterium tuberculosis Complex Isolated from Tuberculosis Patients in Bahir Dar City and Its Surroundings, Northwest Ethiopia
The knowledge of the diversity of strains of Mycobacterium tuberculosis complex (MTBC) species in a specific geographical region can contribute to the control of tuberculosis (TB). This study was conducted to identify the MTBC isolates to the species and spoligotype international type (SIT) level by spoligotyping. A total of 168 MTBC isolates were recovered from TB patients, spoligotyped, and their patterns were compared with those of the strains registered in the SITVIT2 database. Of 168 isolates spoligotyped, 89 patterns were identified. Ninety-eight isolates were clustered into 19 strain groups with clustering percentage of 58.3%. Forty-four strains matched the preexisting SITs in the SITVIT2 database. The dominant strains were SIT289, SIT134, and SIT3411, comprising 16.7% (28/168), 7.14% (12/168), and 4.76% (8/168) of the isolates, respectively. Euro-American (51.2%), East-African-Indian (34.5%), and M. africanum (9.52%) were the major lineages identified. Two strains of M. bovis were isolated from TB lymphadenitis cases. The high percentage of clustered strains of M. tuberculosis could suggest that a small number of lineages of M. tuberculosis are causing the disease in the area while isolation of M. bovis could suggest its zoonotic potential. Additionally, identification of M. africanum requires further confirmation by tools with a better discriminatory power.
The incidence of TB has continued to increase in many parts of the world . An estimated 9 million people developed TB in 2013, which is 6% greater than that reported in 2012. A quarter of these cases occurred in the African region . About 1.5 million deaths were attributed to TB globally, of which approximately 75% occurred in Africa and Southeast Asia . The 22 high TB burden countries (HBCs) collectively accounted for 82% of all the estimated TB incidences worldwide . Ethiopia is among the top ten HBCs with an estimated incidence rate of 224 patients per 100,000 population , which is much greater than 126/100,000 reported globally in 2013 . Moreover, TB is the second most common cause of hospital deaths in the country . TB lymphadenitis in cervical lymph nodes is also common in Ethiopia and accounted for 33% of all new TB cases , which is greater than the global average of 15% . According to the 2013 Heal TB Project Report of 2014, the incidence of TB in the Amhara Region (present study area) for one year (October 2012 to September 2013) was estimated to be 247 per 100,000 population (unpublished, URL: http://pdf.usaid.gov/pdf%20docs/pa00jn8p.pdf), which is higher than the national incidence rate of Ethiopia during the same year. However, little or no information is available on the type of MTBC species and strains causing the disease in the study area.
Thus, identification of strains circulating in a certain geographic region using molecular tools can contribute to the TB control program of that region. Various molecular epidemiology methods allow identifying mycobacterial strains and tracking the transmission of TB in different geographical regions . For example, spoligotyping is widely used for identification of M. tuberculosis . Hence, the purpose of this study was to identify the MTBC isolates to the species and SIT level by the use of spoligotyping.
2. Materials and Methods
2.1. Study Area
The study was conducted in Bahir Dar city and its surrounding districts. Bahir Dar is a city in the highland in northwest of Ethiopia and is the capital city of the Amhara Regional Administrative State. The city is located at geographical coordinates of 11°38′ north in latitude and 37°15′ east in longitude. It has an elevation of 1830 meters above sea level and is characterized by hot and humid weather with an average temperature of 29°C. The population size of the city and its surrounding is 221,991, of which 180,174 (81.2%) are residing in the Bahir Dar city . To seek a better economic situation, rural inhabitants, particularly those living in a radius of 60 km around the city, recently migrated and caused the rapid rise in population in Bahir Dar city . The migration led to homelessness and poverty, as the anticipated job opportunities were not realized. In addition, crowded living conditions, lack of ventilation in temporary housing, malnutrition, and lack of education facilitated the spread of TB .
2.2. Study Subjects
Smear positive pulmonary TB (PTB) and extrapulmonary TB (EPTB) patients, who were diagnosed as TB cases between September 2012 and January 2014 at Felegehiwot Referral and GAMBY General Hospitals, were included in this study. High TB patient flow, existence of better diagnostic facilities, and skilled human resource were the major reasons for selecting the specified health facilities. The average TB case flow in these two study hospitals over four years (2010–2013) was 321 as assessed from the respective hospital records. TB patients who visited these health facilities during the study period were enrolled in the study, excluding those below 18 years of age and those who had started treatment prior to launching the study. Children under the age of 18 years were not included as (1) the study was not a pure epidemiological study and its main objective was the identification of strains circulating in the study area, (2) it was not easy to obtain consent from a family member or a guardian for children below 18 years of age, (3) collection of sputum samples can be difficult in children, and (4) systematic differences of MTBC strains comparing the adult population with children were not expected.
2.3. Sample and Data Collection
A structured questionnaire was used to collect data from all study subjects. These data included patient origin, age, sex, household size, TB category, clinical presentation, and family history of TB infections. Clinical examination of patients suspected to be infected with TB was performed by the attending physicians. Sputum samples submitted for the routine Ziehl Neelsen staining for diagnostic purpose were used for bacteriological examination. Similarly, fine needle aspirates (FNA) collected by a pathologist for the routine diagnosis of TB were used for mycobacterial culture.
2.4. Mycobacterium Culture
Isolation of mycobacteria was made on Lowenstein-Jensen (LJ) medium using the procedure described by the National TB and Leprosy Control Programme Guideline  that was adopted from WHO guideline . Both sputum and FNA samples were cultured at the Bahir Dar Regional Health Research Laboratory Centre. Briefly, sputum or FNA samples were homogenized and decontaminated with an equal volume of 4% NaOH by centrifugation at 3000 rpm for 15 minutes at room temperature. The supernatant was decanted while the sediment was neutralized with 2 N HCl using phenol red as an indicator. Neutralization was achieved when the color of the solution changed from purple to yellow. About 100 µL of the suspension was inoculated on four sterile LJ medium slopes (two were supplemented with pyruvate and the other two with glycerol) and then incubated at 37°C with weekly examination for growth. Specimens without colonies at eighth week after culturing were considered as negative. Specimens with growth of colonies were examined by Ziehl Neelsen microscopy. AFB positive colonies were harvested and resuspended in 200 µL sterile distilled water. Thereafter, they were inactivated by heating at 80°C for 45 minutes in a water bath and transported to Aklilu Lemma Institute of Pathobiology (ALIPB), Addis Ababa University, for spoligotyping.
All of the 168 isolates were characterized by spoligotyping as previously described by Kamerbeek et al.  following the instruction supplied with the spoligotyping kit (Ocimum Biosolutions Company, Ijsselstein, The Netherlands). Briefly, the direct repeat (DR) region of the isolate was amplified using DRa and DRb primers. The amplified biotinylated products were hybridized with a set of 43 oligonucleotides covalently bound to a membrane (Animal and Plant Health Agency, Great Britain). Known strains of M. bovis and M. tuberculosis H37Rv were used as positive controls, whereas Qiagen water (Qiagen Company, Germany) was used as a negative control. Hybridized DNA was detected by the enhanced chemiluminescence method. Images were captured by exposure to X-ray film (Hyperfilm ECL, Amersham) as specified by the manufacturer’s instruction. The presence and absence of spacers were visualized on film as black and white squares, respectively.
2.6. Comparison of Experimental Data with the SITVIT2 Database
The spoligotype patterns were converted into binary and octal formats and entered into the open source spoligotype database available at the website http://www.pasteur-guadeloupe.fr:8081/SITVIT_ONLINE/tools.jsp. The shared international spoligotype (SIT) number and lineages/sublineages were retrieved from the database. The results were compared with the existing designations in the SITVIT2 database of Institute Pasteur de la Guadeloupe. Two or more mycobacterial isolates sharing a spoligotype pattern in the study were identified as a cluster, whereas single spoligotype patterns in the study were recognized as unique. Strains matching a preexisting pattern in the SITVIT2 database were identified with the SIT number, whereas strains for which SIT numbers were not found from the database were considered as orphan strains. In addition, the online tool “Run TB-Lineage” (http://tbinsight.cs.rpi.edu/run_tb_lineage.html) was also used to predict the major lineages to which the strains belong by a conformal Bayesian network (CBN) analysis.
2.7. Statistical Analysis
The statistical analysis was performed using STATA software version 12 . Descriptive statistics were used to depict the demographic and clinical variables. Chi-square or Fisher’s exact tests were used to evaluate the association of clusters and major lineages with selected patient characteristics. values of less than 0.05 were considered statistically significant.
2.8. Ethical Considerations
Ethical clearance was obtained from Ethical Review Board (Ref. number IRB/05-02/2013) of the Aklilu Lemma Institute of Pathobiology, Addis Ababa University. In addition, permission was obtained from the Research Committee of Bureau of Health, Amhara National Regional State, Ethiopia.
3.1. Demographic and Clinical Characteristics of the Study Subjects
Data generated from 168 subjects were used in the analysis of the demographic and clinical results. Among the study participants, 52.4% were female, 73.8% were in age range of 18–39 years, 84.5% were new cases, 27.4% had a history of TB pertaining to one of their family members, and 67.9% were EPTB patients. Surprisingly, all EPTB cases were identified as TB lymphadenitis (TBLN), of which 67 (60.4%) and 18 (16.2%) were TBLN in cervical and axillary lymph nodes, respectively. Of the 168 isolates, 33.9% and 25.6% originated from South Gondar and West Gojjam, respectively (Table 1). Nonetheless, the sociodemographic and clinical characteristics of the patients did not affect the clustering rates and distribution of the lineages of MTBC strains (Table 1).
|EA: Euro-American; EAI: East-African Indian; IO: Indo-Oceanic; MA: M. africanum; MB: M. bovis; CBN: conformal Bayesian network; values were presented at 95% confidence interval and considered statistically significant.|
3.2. Spoligotyping Patterns of Mycobacterium tuberculosis Complex Strains
A total of 168 MTBC isolates were spoligotyped, and 89 (53%) different spoligotype patterns (strains) were identified. Clustering of isolates into strains was observed, and a total of 98 isolates were grouped in 19 (58.3%) different clusters of strains. The dominant strains were SIT289, SIT134, and SIT3411, each consisting of 28 (16.7%), 12 (7.14%), and 8 (4.76%) isolates, respectively. These strains contributed 28.6% (48/168) of all isolates with known spoligotype patterns. Out of the 89 spoligotype patterns (strains), 44 strains associated with 122 isolates matched the preexisting patterns in the SITVIT2 database while the remaining 45 spoligotype patterns associated with 46 isolates were not registered in the international spoligotype SITVIT2 database and thus designated as orphan strains.
Classification of MTBC strains showed the occurrence of the following lineages: Euro-American (86/168; 51.2%), East-African Indian (58/168; 34.5%), M. africanum (16/168; 9.52%), Indo-Oceanic (6/168; 3.57%), and M. bovis (2/168; 1.19%). Two of the spoligotypes (i.e., SIT982 and SIT665) were M. bovis and both were isolated from EPTB cases. Detailed spoligotyping results and their corresponding SITs/orphan strains and lineages are summarized in Tables 2 and 3.
|CBN: conformal Bayesian network; unknown: designates patterns with signatures that do not belong to any of the major lineages/sublineages described in the SITVIT2 database. The 168 isolates were grouped into 89 different spoligotype patterns (strains). Of the total 89 strains, 44 strains (patterns) have already been registered in the SITVIT2 database (Table 2), while the remaining 45 patters were orphans and presented in Table 3. The dominant strains were SIT289 (28 isolates), SIT134 (12 isolates), and SIT3411 (8 isolates) (Table 2). Furthermore, the 168 isolates were grouped into five different lineages including the Euro-American, East-African Indian, M. africanum, Indo-Oceanic, and M. bovis lineages in the order of decreasing percentage.|
|CBN: conformal Bayesian network; unknown: designates patterns with signatures that do not belong to any of the major lineages/sublineages described in the SITVIT2 database. Forty-five of the total 89 strains were identified as orphan strains in the present study and shown in Table 3. The orphan strains belonged to four lineages including Euro-American, East-African Indian, M. africanum, and Indo-Oceanic lineages. |
3.3. Distribution of Strains and Lineages in the Study Area
The majority of MTBC strains were identified from the South Gondar Zone (57/168; 33.9%) followed by the West Gojjam Zone (43/168; 25.6%), each with a strain-clustering rate of 17.3% (Table 1). The distribution of the three dominant strains (SIT289, SIT134, and SIT3411) in the area is depicted in Figure 1.
The Euro-American (EA), East-African-Indian (EAI), Indo-Oceanic (IO), and M. africanum (MA) lineages were identified in all study zones, whereas M. bovis (MB) was recorded only from patients with TBLN located in South Gondar. Figure 2 depicts the distribution of the major lineages.
In the present study, MTBC species were isolated from 168 TB patients from Bahir Dar city itself and the surrounding districts who visited health institutions in Bahir Dar city. The isolates were identified at strain and lineage levels on the basis of spoligotyping. Identification at a higher level of resolution by using variable nucleotide tandem repeat (VNTR) typing, facilitated by mycobacterial interspersed repetitive units- (MIRU-) VNTR, is desirable although this method is not currently available in Ethiopia.
Spoligotyping of 168 mycobacterial isolates revealed 89 distinct patterns, which corresponded to 53% of genotype diversity. The high diversity of spoligotypes strains that we observed in this study was consistent with the 59% reported by Tessema et al.  but higher than the percentages reported earlier by other studies in Ethiopia [3, 15–18]. This finding suggests the circulation of genetically variable strains in the study area, which could be the result of significant migration of infected people to Bahir Dar city and its surroundings from other regions of the country. In addition, the long period of MTBC clonal evolution may contribute to the diversity of strains . Ninety-eight mycobacterial isolates were grouped into 19 clusters with an overall clustering percentage of 58.3%. The clustering rate observed in this study was slightly higher than those reported previously in Ethiopia [14, 18, 20]. On the other hand, it was lower than those reported by several other national studies [15, 16, 21] and international studies (e.g., South Africa  and Malawi ). The observed differences in clustering rates might be related to differences of sociocultural origin, sanitation, and population density. High level of strain clustering could suggest recent and ongoing TB transmission .
The prevalent strains identified in this study were SIT289, SIT134, and SIT3411. All three strains seem to be specific for the Bahir Dar city and its surroundings since they were not reported previously from other sites in Ethiopia [3, 14, 17, 18]. However, in the SITVIT database, SIT289 has only been reported from Brazil and Europe (mainly France, French Guiana, Martinique, Italy, Netherlands, and Sweden), while SIT134 has been reported from Central Asia and Middle East (Bangladesh, India, Pakistan, and Saudi Arabia), Australia, Netherlands, and United States of America . A considerable number (46/168) of orphan strains were also recorded in this study. This is nearly identical to the average reported by Tessema et al. , Belay et al. , and Mihret et al.  in Ethiopia. The existence of mixed infections may also complicate spoligotyping results [6, 26], and hence higher resolution molecular tools should be applied toreveal thus far undefined mixed spoligotyping signatures.
Five different major lineages, namely, Euro-American, East-African-Indian, Indo-Oceanic, M. africanum, and M. bovis, were identified in this study. Euro-American was the dominant lineage, and more than half (51.5%) of the overall strains belonged to this lineage. This finding agreed with the results of previous studies in Ethiopia [3, 16–18] and Morocco . The high proportion of new MTBC lineages is supposed to be related to their successful geographical spread as compared to ancient lineages [28, 29]. Even though Euro-American was identified as the prevalent major lineage, the CAS1_DELHI sublineage (consisting of 46 isolates) in the East-African Indian lineage appeared to have had a high transmission rate in our study population. This lineage is localized in South Asia, preferentially India, countries of the Middle East, and several other regions, including Africa . It can be hypothesized that East-African Indian ancestral strains spread back from Asia to Africa through India as a result of human migration .
Screening of the SITVIT2 database also identified 9.52% (16/168) of the isolates as members of M. africanum. The clustering rate was 1.19% (2/168), indicative of a low rate of recent human-to-human transmission. Since it was not reported previously in Ethiopia [3, 14–18], isolation and identification of M. africanum in this study represent a novel finding. Further studies are needed to explore evolutionary aspects that may have contributed to the spread of M. africanum in the study population. Two strains of M. bovis (SIT982 and SIT665) were identified in this study. This finding was interesting and could implicate the public health importance of M. bovis in northwestern Ethiopia.
Molecular characterization of MTBC isolates from TB patients in Bahir Dar city and its surroundings was performed using spoligotyping. The high percentage of clustered strains of M. tuberculosis could suggest that a small number of lineages of M. tuberculosis are causing the disease in the area and isolation of M. bovis could suggest its zoonotic potential in the study area. Meanwhile, identification of M. africanum requires confirmation by molecular tools with a better discriminatory power than spoligotyping.
Conflict of Interests
The authors declared that they have no competing interests.
Anwar Nuru was involved in the design, data collection and laboratory work, and statistical analysis and interpretation of data and drafted the paper; Gezahegne Mamo was involved in the design, interpretation, and critical revision of the paper; Adane Worku was involved in the spoligotyping of the isolates and interpretation; Aschalew Admasu was involved in data collection, TB culture, and interpretation; Girmay Medhin was involved in the design, statistical analysis, and critical revision of the paper; Rembert Pieper was involved in critical revision of the paper; Gobena Ameni was involved in the design, interpretation, and critical revision of the paper. All authors have read and approved the paper for submission.
This study was jointly funded by the National Institute of Health (NIH, USA) through its H3Africa Consortium Program (Grant Ref. no. U01HG007472-01), Addis Ababa University through its Thematic Research Program, and University of Gondar. The authors are grateful to Dr. Getachew Mulualem and the other technical personnel for collecting samples as well as the institutions which facilitated the conducting of this study. The authors would like to acknowledge David Couvin for comparing the experimental strains with SITVIT2 database and Eyaya Belay for his support in doing the GIS maps.
- WHO, Global Tuberculosis Report, WHO, Geneva, Switzerland, 2014.
- 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, 15th edition, 2012.
- R. Firdessa, S. Berg, E. Hailu et al., “Mycobacterial lineages causing pulmonary and extrapulmonary Tuberculosis, Ethiopia,” Emerging Infectious Diseases, vol. 19, no. 3, pp. 460–463, 2013.
- WHO, Global Tuberculosis Report, WHO, Geneva, Switzerland, 2011.
- D. García De Viedma, I. Mokrousov, and N. Rastogi, “Innovations in the molecular epidemiology of tuberculosis,” Enfermedades Infecciosas y Microbiología Clínica, vol. 29, no. 1, pp. 8–13, 2011.
- K. Brudey, J. R. Driscoll, L. Rigouts et al., “Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology,” BMC Microbiology, vol. 6, article 23, 2006.
- Central Statistics Agency, Summary and Statistical Report of the 2007 Population and Housing Census: Population Size by Age and Sex, Central Statistics Agency, Addis Ababa, Ethiopia, 2007.
- M. Ali, “Status of homeless population in urban Ethiopia: a case study of Amhara region,” International Journal of Management and Social Sciences Research, vol. 3, no. 1, pp. 61–68, 2014.
- C. W. Schmidt, “Linking TB and the environment: an overlooked mitigation strategy,” Environmental Health Perspectives, vol. 116, no. 11, pp. A478–A485, 2008.
- Federal Ministry of Health, Tuberculosis, Leprosy and TB/HIV Prevention and Control Programme, Federal Ministry of Health, Addis Ababa, Ethiopia, 4th edition, 2008.
- WHO, Laboratory Services in Tuberculosis Control: Culture. Part III, World Health Organization (WHO), Geneva, Switzerland, 1998.
- 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.
- StataCorp LP, Stata Statistical Software: Release 12, StataCorp LP, College Station, Tex, USA, 2011.
- 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, no. 1, article 131, 2013.
- M. Belay, G. Ameni, G. Bjune, D. Couvin, N. Rastogi, and F. 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.
- A. Mihret, Y. Bekele, M. Aytenew et al., “Modern lineages of Mycobacterium tuberculosis in Addis Ababa, Ethiopia: implications for the tuberculosis control programe,” African Journal of Health Sciences, vol. 12, no. 3, pp. 339–344, 2012.
- 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.
- T. Debebe, A. Admassu, G. Mamo, and G. Ameni, “Molecular characterization of Mycobacterium tuberculosis isolated from pulmonary tuberculosis patients in Felege Hiwot Referral Hospital, northwest Ethiopia,” Journal of Microbiology, Immunology and Infection, vol. 47, no. 4, pp. 333–338, 2014.
- 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.
- J. Bruchfeld, G. Aderaye, I. B. Palme et al., “Molecular epidemiology and drug resistance of Mycobacterium tuberculosis isolates from ethiopian pulmonary tuberculosis patients with and without human immunodeficiency virus infection,” Journal of Clinical Microbiology, vol. 40, no. 5, pp. 1636–1643, 2002.
- B. Diriba, T. Berkessa, G. Mamo, Y. Tedla, and G. Ameni, “Spoligotyping of multidrug-resistant Mycobacterium tuberculosis isolates in Ethiopia,” International Journal of Tuberculosis and Lung Disease, vol. 17, no. 2, pp. 246–250, 2013.
- S. Verver, R. M. Warren, Z. Munch et al., “Transmission of tuberculosis in a high incidence urban community in South Africa,” International Journal of Epidemiology, vol. 33, no. 2, pp. 351–357, 2004.
- J. R. Glynn, A. C. Crampin, M. D. Yates et al., “The importance of recent infection with Mycobacterium tuberculosis in an area with high HIV prevalence: a long-term molecular epidemiological study in Northern Malawi,” Journal of Infectious Diseases, vol. 192, no. 3, pp. 480–487, 2005.
- P. J. Easterbrook, A. Gibson, S. Murad et al., “High rates of clustering of strains causing tuberculosis in Harare, Zimbabwe: a molecular epidemiological study,” Journal of Clinical Microbiology, vol. 42, no. 10, pp. 4536–4544, 2004.
- 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.
- I. C. Shamputa, L. Rigouts, L. A. Eyongeta et al., “Genotypic and phenotypic heterogeneity among Mycobacterium tuberculosis isolates from pulmonary tuberculosis patients,” Journal of Clinical Microbiology, vol. 42, no. 12, pp. 5528–5536, 2004.
- O. Lahlou, J. Millet, I. Chaoui et al., “The genotypic population structure of Mycobacterium tuberculosis complex from Moroccan patients reveals a predominance of Euro-American lineages,” PLoS ONE, vol. 7, no. 10, Article ID e47113, 2012.
- S. Gagneux, “Host-pathogen coevolution in human tuberculosis,” Philosophical Transactions of the Royal Society B: Biological Sciences, vol. 367, no. 1590, pp. 850–859, 2012.
- M. Coscolla and S. Gagneux, “Consequences of genomic diversity in Mycobacterium tuberculosis,” Seminars in Immunology, vol. 26, no. 6, pp. 431–444, 2014.
- F. Cruciani, P. Santolamazza, P. Shen et al., “A back migration from Asia to sub-Saharan Africa is supported by high-resolution analysis of human Y-chromosome haplotypes,” American Journal of Human Genetics, vol. 70, no. 5, pp. 1197–1214, 2002.
Copyright © 2015 Anwar Nuru 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.