Patients with suspected TBLN or pulmonary TB who came to hospitals or health centers in study sites and provided voluntary consent were recruited into the study during 2006–2010. Fine needle–aspirate samples and sputum samples were collected from 2,151 patients attending hospitals in Gondar, Woldiya, Ghimbi, Butajira, and Negelle, Ethiopia. In addition, sputum samples were collected from patients at hospitals in Fiche, Jinka, and Filtu and at health centers at 3 suburban sites in Addis Ababa (Holeta, Sululta, and Chancho). Samples were cultured on Löwenstein-Jensen medium supplemented with glycerol or pyruvate and on modified Middlebrook 7H11 medium optimized for culture of M. bovis.

We characterized isolates belonging to the M. tuberculosis complex by using multiplex PCR for large sequence polymorphisms (3,4), spoligotyping (5), and lineage-specific single-nucleotide polymorphism analysis (4,6). Isolates of selected spoligotypes were characterized by 24-loci mycobacterial interspersed repetitive unit–variable number tandem repeat (MIRU-VNTR) analysis (7). Four M. tuberculosis isolates from a group of 36 with unusual spoligotype patterns were further characterized by genome sequencing (Illumina Inc., San Diego, CA, USA). Sequencing reads were mapped to the inferred most recent common ancestor of the M. tuberculosis complex (6). A final alignment of 13,199 single-nucleotide polymorphism positions was generated and analyzed by using the neighbor-joining method with a Tamura-Nei evolutionary model (www.megasoftware.net/mega_papers.php). Nontuberculous mycobacteria were characterized by sequencing of the 16S rDNA gene.

Characteristics of 964 cultures positive for acid-fast bacilli are summarized in Table 1. Most of these isolates had an intact RD9 region, which identified them as M. tuberculosis. Only 4 (0.4%) of 964 isolates had undergone RD9 and RD4 deletions characteristic of M. bovis (3). The 4 M. bovis isolates were obtained from cases of pulmonary TB, 3 of which were from patients living in pastoralist communities in southern Ethiopia. The 10 nontuberculous mycobacterial isolates were identified as M. intracellulare, M. flavescens, and M. simiae; 2 of the isolates were from patients co-infected with M. tuberculosis.

Among the 954 isolates belonging to the M. tuberculosis complex, 671 (71%) belonged to lineage 4, which was the most common lineage in Ethiopia. However, lineage 3 was most prevalent in the northern sites of Gondar and Woldiya (122/257, 47%). Eleven strains belonging to lineage 1 were isolated in the southern region. Two isolates with a characteristic Beijing family spoligotype (spoligotype international type [SIT] 1) were identified as pseudo-Beijing strains belonging to lineage 3 (8). Thirty-six isolates with an unusual spoligotype pattern (missing spacers 4–24) and intact for the TbD1 region could not be assigned to known lineages. Genome sequencing identified these strains as members of a new lineage (lineage 7) localized between ancient lineage 1 and modern lineages 2, 3, and 4 of M. tuberculosis phylogeny (Figure 1). This new lineage 7 was prominent among strains collected in the Woldiya region (17/133 strains, 13%) (Table 1).

Lineage distribution was identical between the 2 disease forms at the national level; lineage 4 was isolated from 71% (442/662) and 70% (229/328) of pulmonary TB and TBLN patients, respectively, and lineage 3 was isolated from 25% (153/622) and 24% (79/328), respectively. The M. tuberculosis isolates encompassed 176 spoligotypes (Technical Appendix Table 1), of which 86 patterns were new to the international genotyping database 2 (SITVIT2) (www.pasteur-guadeloupe.fr:8081/SITVITDemo/) (9). A total of 11% (101/950) of the isolates represented single spoligotypes, and 62% (591/950) were included in 10 major spoligotype clusters (Table 2).

There was no difference in cluster distribution between pulmonary TB and lymph node TB isolates; 10% and 11%, respectively, were single types, and 64% and 60%, respectively, were included in dominant clusters. Two large clusters representative of lineage 4 (SIT 149) and lineage 3 (SIT 25) were further characterized by MIRU-VNTR typing (Technical Appendix Table 2) and network analysis (Figure 2). In each case, TBLN and pulmonary TB samples were dispersed throughout the network of spoligotype clusters.

All 4 M. bovis isolates from humans showed typical bovine spoligotype profiles lacking spacers 3, 9, 16, and 39–43. In addition, they lacked spacers 4–7 and had deletions of RDAf2, which are features that define strains of the African 2 clonal complex of M. bovis reported from TB-infected cattle in Ethiopia (10).

The frequency of M. bovis in persons in this study (0.4%) is similar to that found in other studies of human TB in Africa (11) and South and Central America (12), but much lower than that observed among selected populations in Tanzania (16%) (13), Ethiopia (17%) (14), and Mexico (28%) (15). These findings indicate that the overall contribution of M. bovis to human TB is minor but greater in specific areas. In Ethiopia, monitoring of zoonotic transmission is needed in urban areas with high rates of bovine TB associated with intensive farming of imported dairy cattle (R. Firdessa et al., unpub. data) and among pastoralist populations from which human M. bovis cases were identified in this study.

Zoonotic transmission of M. bovis can be excluded as the predominant cause of the high national incidence of TBLN in Ethiopia. Mapping of disease networks by spoligotyping and MIRU-VNTR analysis showed an integrated distribution of the 2 disease forms, which suggested that cases of TBLN arise from within the pulmonary TB transmission network, rather than from an external zoonotic source.

We identified a novel phylogenetic lineage of M. tuberculosis (lineage 7) in multiple sites and at a high frequency in Woldiya in the northeastern highlands of Ethiopia. Screening of the SITVIT2 database (9) and the US Centers for Disease Control and Prevention National Tuberculosis Genotyping Surveillance Network Database (L.S. Cowan, pers. comm.) identified 23 (0.03%) of ≈90,000 isolates as members of lineage 7; all were isolated from patients whose country of origin (when known) was in the Horn of Africa. Lineage 7 is of considerable evolutionary interest because it represents a phylogenetic branch intermediate between the ancient and modern lineages of M. tuberculosis (3,4,6).

Strain frequencies of respective spoligotype international types (SITs) of Mycobacterium tuberculosis isolated in this study, results of large and small sequence polymorphism typing, and typing results of 24-loci mycobacterial interspersed repetitive unit–variable number tandem repeats for isolates of spoligotypes SIT 25 and SIT 149, Ethiopia, 2006–2010.