Tuberculous Lymphadenitis in Ethiopia Predominantly Caused by Strains Belonging to the Delhi/CAS Lineage and Newly Identified Ethiopian Clades of the Mycobacterium tuberculosis Complex

Background

Recently, newly defined clades of Mycobacterium tuberculosis complex (MTBC) strains, namely Ethiopia 1–3 and Ethiopia H37Rv-like strains, and other clades associated with pulmonary TB (PTB) were identified in Ethiopia. In this study, we investigated whether these new strain types exhibit an increased ability to cause TB lymphadenitis (TBLN) and raised the question, if particular MTBC strains derived from TBLN patients in northern Ethiopia are genetically adapted to their local hosts and/or to the TBLN.

Methods

Genotyping of 196 MTBC strains isolated from TBLN patients was performed by spoligotyping and 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTR) typing. A statistical analysis was carried out to see possible associations between patient characteristics and phylogenetic MTBC strain classification.

Results

Among 196 isolates, the majority of strains belonged to the Delhi/CAS (38.8%) lineage, followed by Ethiopia 1 (9.7%), Ethiopia 3 (8.7%), Ethiopia H37RV-like (8.2%), Ethiopia 2 and Haarlem (7.7% each), URAL (3.6%), Uganda l and LAM (2% each), S-type (1.5%), X-type (1%), and 0.5% isolates of TUR, EAI, and Beijing genotype, respectively. Overall, 15 strains (7.7%) could not be allocated to a previously described phylogenetic lineage. The distribution of MTBC lineages is similar to that found in studies of PTB samples. The cluster rate (35%) in this study is significantly lower (P = 0.035) compared to 45% in the study of PTB in northwestern Ethiopia.

Conclusion

In the studied area, lymph node samples are dominated by Dehli/CAS genotype strains and strains of largely not yet defined clades based on MIRU-VNTR 24-loci nomenclature. We found no indication that strains of particular genotypes are specifically associated with TBLN. However, a detailed analysis of specific genetic variants of the locally contained Ethiopian clades by whole genome sequencing may reveal new insights into the host-pathogen co-evolution and specific features that are related to the local host immune system.

Whole-genome sequencing for prediction of Mycobacterium tuberculosis drug susceptibility and resistance: a retrospective cohort study

BACKGROUND

Diagnosing drug-resistance remains an obstacle to the elimination of tuberculosis. Phenotypic drug-susceptibility testing is slow and expensive, and commercial genotypic assays screen only common resistance-determining mutations. We used whole-genome sequencing to characterise common and rare mutations predicting drug resistance, or consistency with susceptibility, for all first-line and second-line drugs for tuberculosis.

METHODS

Between Sept 1, 2010, and Dec 1, 2013, we sequenced a training set of 2099 Mycobacterium tuberculosis genomes. For 23 candidate genes identified from the drug-resistance scientific literature, we algorithmically characterised genetic mutations as not conferring resistance (benign), resistance determinants, or uncharacterised. We then assessed the ability of these characterisations to predict phenotypic drug-susceptibility testing for an independent validation set of 1552 genomes. We sought mutations under similar selection pressure to those characterised as resistance determinants outside candidate genes to account for residual phenotypic resistance.

FINDINGS

We characterised 120 training-set mutations as resistance determining, and 772 as benign. With these mutations, we could predict 89·2% of the validation-set phenotypes with a mean 92·3% sensitivity (95% CI 90·7-93·7) and 98·4% specificity (98·1-98·7). 10·8% of validation-set phenotypes could not be predicted because uncharacterised mutations were present. With an in-silico comparison, characterised resistance determinants had higher sensitivity than the mutations from three line-probe assays (85·1% vs 81·6%). No additional resistance determinants were identified among mutations under selection pressure in non-candidate genes.

INTERPRETATION

A broad catalogue of genetic mutations enable data from whole-genome sequencing to be used clinically to predict drug resistance, drug susceptibility, or to identify drug phenotypes that cannot yet be genetically predicted. This approach could be integrated into routine diagnostic workflows, phasing out phenotypic drug-susceptibility testing while reporting drug resistance early.

FUNDING

Wellcome Trust, National Institute of Health Research, Medical Research Council, and the European Union.

Integration of published information into a resistance-associated mutation database for Mycobacterium tuberculosis

Tuberculosis remains a major global public health challenge. Although incidence is decreasing, the proportion of drug-resistant cases is increasing. Technical and operational complexities prevent Mycobacterium tuberculosis drug susceptibility phenotyping in the vast majority of new and retreatment cases. The advent of molecular technologies provides an opportunity to obtain results rapidly as compared to phenotypic culture. However, correlations between genetic mutations and resistance to multiple drugs have not been systematically evaluated. Molecular testing of M. tuberculosis sampled from a typical patient continues to provide a partial picture of drug resistance. A database of phenotypic and genotypic testing results, especially where prospectively collected, could document statistically significant associations and may reveal new, predictive molecular patterns. We examine the feasibility of integrating existing molecular and phenotypic drug susceptibility data to identify associations observed across multiple studies and demonstrate potential for well-integrated M. tuberculosis mutation data to reveal actionable findings.

Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage

Mycobacterium tuberculosis strains of the Beijing lineage are globally distributed and are associated with the massive spread of multidrug-resistant (MDR) tuberculosis in Eurasia. Here we reconstructed the biogeographical structure and evolutionary history of this lineage by genetic analysis of 4,987 isolates from 99 countries and whole-genome sequencing of 110 representative isolates. We show that this lineage initially originated in the Far East, from where it radiated worldwide in several waves. We detected successive increases in population size for this pathogen over the last 200 years, practically coinciding with the Industrial Revolution, the First World War and HIV epidemics. Two MDR clones of this lineage started to spread throughout central Asia and Russia concomitantly with the collapse of the public health system in the former Soviet Union. Mutations identified in genes putatively under positive selection and associated with virulence might have favored the expansion of the most successful branches of the lineage.

Pre-Columbian mycobacterial genomes reveal seals as a source of New World human tuberculosis

Modern strains of Mycobacterium tuberculosis from the Americas are closely related to those from Europe, supporting the assumption that human tuberculosis was introduced post-contact. This notion, however, is incompatible with archaeological evidence of pre-contact tuberculosis in the New World. Comparative genomics of modern isolates suggests that M. tuberculosis attained its worldwide distribution following human dispersals out of Africa during the Pleistocene epoch, although this has yet to be confirmed with ancient calibration points. Here we present three 1,000-year-old mycobacterial genomes from Peruvian human skeletons, revealing that a member of the M. tuberculosis complex caused human disease before contact. The ancient strains are distinct from known human-adapted forms and are most closely related to those adapted to seals and sea lions. Two independent dating approaches suggest a most recent common ancestor for the M. tuberculosis complex less than 6,000 years ago, which supports a Holocene dispersal of the disease. Our results implicate sea mammals as having played a role in transmitting the disease to humans across the ocean.

Phylogenetic polymorphisms in antibiotic resistance genes of the Mycobacterium tuberculosis complex

OBJECTIVES

Sequence analysis of known antibiotic resistance genes of the Mycobacterium tuberculosis complex (MTBC) is increasingly being used to infer phenotypic resistance to a variety of antibiotics. However, a clear understanding of the genotype-phenotype relationship is required to interpret genotypic susceptibility results accurately. In this context, it is particularly important to distinguish phylogenetically informative neutral polymorphisms from true resistance-conferring mutations.

METHODS

Using a collection of 71 strains that encompasses all major MTBC genotypes, we mapped the genetic diversity in 18 genes that are known to be involved or were previously implicated in antibiotic resistance to eight current as well as two novel antibiotics. This included bedaquiline, capreomycin, ethambutol, fluoroquinolones, isoniazid, PA-824, para-aminosalicylic acid, prothionamide, rifampicin and streptomycin. Moreover, we included data from one of our prior studies that focused on two of the three known pyrazinamide resistance genes.

RESULTS

We found 58 phylogenetic polymorphisms that were markers for the genotypes M. tuberculosis Beijing, Haarlem, Latin American-Mediterranean (LAM), East African Indian (EAI), Delhi/Central Asian (CAS), Ghana, Turkey (Tur), Uganda I and II, Ural and X-type, as well as for Mycobacterium africanum genotypes West African I (WA I) and II (WA II), Mycobacterium bovis, Mycobacterium caprae, Mycobacterium pinnipedii, Mycobacterium microti and Mycobacterium canettii.

CONCLUSIONS

This study represents one of the most extensive overviews of phylogenetically informative polymorphisms in known resistance genes to date, and will serve as a resource for the design and interpretation of genotypic susceptibility assays.

Whole genome sequencing reveals complex evolution patterns of multidrug-resistant Mycobacterium tuberculosis Beijing strains in patients

Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains represent a major threat for tuberculosis (TB) control. Treatment of MDR-TB patients is long and less effective, resulting in a significant number of treatment failures. The development of further resistances leads to extensively drug-resistant (XDR) variants. However, data on the individual reasons for treatment failure, e.g. an induced mutational burst, and on the evolution of bacteria in the patient are only sparsely available. To address this question, we investigated the intra-patient evolution of serial MTBC isolates obtained from three MDR-TB patients undergoing longitudinal treatment, finally leading to XDR-TB. Sequential isolates displayed identical IS6110 fingerprint patterns, suggesting the absence of exogenous re-infection. We utilized whole genome sequencing (WGS) to screen for variations in three isolates from Patient A and four isolates from Patient B and C, respectively. Acquired polymorphisms were subsequently validated in up to 15 serial isolates by Sanger sequencing. We determined eight (Patient A) and nine (Patient B) polymorphisms, which occurred in a stepwise manner during the course of the therapy and were linked to resistance or a potential compensatory mechanism. For both patients, our analysis revealed the long-term co-existence of clonal subpopulations that displayed different drug resistance allele combinations. Out of these, the most resistant clone was fixed in the population. In contrast, baseline and follow-up isolates of Patient C were distinguished each by eleven unique polymorphisms, indicating an exogenous re-infection with an XDR strain not detected by IS6110 RFLP typing. Our study demonstrates that intra-patient microevolution of MDR-MTBC strains under longitudinal treatment is more complex than previously anticipated. However, a mutator phenotype was not detected. The presence of different subpopulations might confound phenotypic and molecular drug resistance tests. Furthermore, high resolution WGS analysis is necessary to accurately detect exogenous re-infection as classical genotyping lacks discriminatory power in high incidence settings.

Novel Mycobacterium tuberculosis complex isolate from a wild chimpanzee

Tuberculosis (TB) is caused by gram-positive bacteria known as the Mycobacterium tuberculosis complex (MTBC). MTBC include several human-associated lineages and several variants adapted to domestic and, more rarely, wild animal species. We report an M. tuberculosis strain isolated from a wild chimpanzee in Côte d’Ivoire that was shown by comparative genomic and phylogenomic analyses to belong to a new lineage of MTBC, closer to the human-associated lineage 6 (also known as M. africanum West Africa 2) than to the other classical animal-associated MTBC strains. These results show that the general view of the genetic diversity of MTBC is limited and support the possibility that other MTBC variants exist, particularly in wild mammals in Africa. Exploring this diversity is crucial to the understanding of the biology and evolutionary history of this widespread infectious disease.

Molecular epidemiology and transmission dynamics of Mycobacterium tuberculosis in Northwest Ethiopia: new phylogenetic lineages found in Northwest Ethiopia

BACKGROUND

Although Ethiopia ranks seventh among the world's 22 high-burden tuberculosis (TB) countries, little is known about strain diversity and transmission. In this study, we present the first in-depth analysis of the population structure and transmission dynamics of Mycobacterium tuberculosis strains from Northwest Ethiopia.

METHODS

In the present study, 244 M. tuberculosis isolates where analysed by mycobacterial interspersed repetitive unit - variable number tandem repeat 24-loci typing and spoligotyping methods to determine phylogenetic lineages and perform cluster analysis. Clusters of strains with identical genotyping patterns were considered as an indicator for the recent transmission.

RESULTS

Of 244 isolates, 59.0% were classified into nine previously described lineages: Dehli/CAS (38.9%), Haarlem (8.6%), Ural (3.3%), LAM (3.3%), TUR (2.0%), X-type (1.2%), S-type (0.8%), Beijing (0.4%) and Uganda II (0.4%). Interestingly, 31.6% of the strains were grouped into four new lineages and were named as Ethiopia_3 (13.1%), Ethiopia_1 (7.8%), Ethiopia_H37Rv like (7.0%) and Ethiopia_2 (3.7%) lineages. The remaining 9.4% of the isolates could not be assigned to the known or new lineages. Overall, 45.1% of the isolates were grouped in clusters, indicating a high rate of recent transmission.

CONCLUSIONS

This study confirms a highly diverse M. tuberculosis population structure, the presence of new phylogenetic lineages and a predominance of the Dehli/CAS lineage in Northwest Ethiopia. The high rate of recent transmission indicates defects of the TB control program in Northwest Ethiopia. This emphasizes the importance of strengthening laboratory diagnosis of TB, intensified case finding and treatment of TB patients to interrupt the chain of transmission.

First molecular epidemiology study of Mycobacterium tuberculosis in Kiribati

Tuberculosis incidence rates in Kiribati are among the highest in the Western Pacific Region, however the genetic diversity of circulating Mycobacterium tuberculosis complex strains (MTBC) and transmission dynamics are unknown. Here, we analysed MTBC strains isolated from culture positive pulmonary tuberculosis (TB) cases from the main TB referral centre between November 2007 and October 2009. Strain genotyping (IS6110 typing, spoligotyping, 24-loci MIRU-VNTR and SNP typing) was performed and demographic information collected. Among 73 MTBC strains analysed, we identified seven phylogenetic lineages, dominated by Beijing strains (49%). Beijing strains were further differentiated in two main branches, Beijing-A (n = 8) and -B (n = 28), that show distinct genotyping patterns and are characterized by specific deletion profiles (Beijing A: only RD105, RD207 deleted; Beijing B: RD150 and RD181 additionally deleted). Many Kiribati strains (59% based on IS6110 typing of all strains) occurred in clusters, suggesting ongoing local transmission. Beijing-B strains and over-crowded living conditions were associated with strain clustering (likely recent transmission), however little evidence of anti-tuberculous drug resistance was observed. We suggest enhanced case finding amongst close contacts and continued supervised treatment of all identified cases using standard first-line drugs to reduce TB burden in Kiribati. Beijing strains can be subdivided in different principle branches that might be associated with differential spreading patterns in the population.

Whole genome sequencing versus traditional genotyping for investigation of a Mycobacterium tuberculosis outbreak: a longitudinal molecular epidemiological study

BACKGROUND

Understanding Mycobacterium tuberculosis (Mtb) transmission is essential to guide efficient tuberculosis control strategies. Traditional strain typing lacks sufficient discriminatory power to resolve large outbreaks. Here, we tested the potential of using next generation genome sequencing for identification of outbreak-related transmission chains.

METHODS AND FINDINGS

During long-term (1997 to 2010) prospective population-based molecular epidemiological surveillance comprising a total of 2,301 patients, we identified a large outbreak caused by an Mtb strain of the Haarlem lineage. The main performance outcome measure of whole genome sequencing (WGS) analyses was the degree of correlation of the WGS analyses with contact tracing data and the spatio-temporal distribution of the outbreak cases. WGS analyses of the 86 isolates revealed 85 single nucleotide polymorphisms (SNPs), subdividing the outbreak into seven genome clusters (two to 24 isolates each), plus 36 unique SNP profiles. WGS results showed that the first outbreak isolates detected in 1997 were falsely clustered by classical genotyping. In 1998, one clone (termed "Hamburg clone") started expanding, apparently independently from differences in the social environment of early cases. Genome-based clustering patterns were in better accordance with contact tracing data and the geographical distribution of the cases than clustering patterns based on classical genotyping. A maximum of three SNPs were identified in eight confirmed human-to-human transmission chains, involving 31 patients. We estimated the Mtb genome evolutionary rate at 0.4 mutations per genome per year. This rate suggests that Mtb grows in its natural host with a doubling time of approximately 22 h (400 generations per year). Based on the genome variation discovered, emergence of the Hamburg clone was dated back to a period between 1993 and 1997, hence shortly before the discovery of the outbreak through epidemiological surveillance.

CONCLUSIONS

Our findings suggest that WGS is superior to conventional genotyping for Mtb pathogen tracing and investigating micro-epidemics. WGS provides a measure of Mtb genome evolution over time in its natural host context.

Drug resistance-conferring mutations in Mycobacterium tuberculosis from Madang, Papua New Guinea

BACKGROUND

Monitoring drug resistance in Mycobacterium tuberculosis is essential to curb the spread of tuberculosis (TB). Unfortunately, drug susceptibility testing is currently not available in Papua New Guinea (PNG) and that impairs TB control in this country. We report for the first time M. tuberculosis mutations associated with resistance to first and second-line anti-TB drugs in Madang, PNG. A molecular cluster analysis was performed to identify M. tuberculosis transmission in that region.

RESULTS

Phenotypic drug susceptibility tests showed 15.7% resistance to at least one drug and 5.2% multidrug resistant (MDR) TB. Rifampicin resistant strains had the rpoB mutations D516F, D516Y or S531L; Isoniazid resistant strains had the mutations katG S315T or inhA promoter C15T; Streptomycin resistant strains had the mutations rpsL K43R, K88Q, K88R), rrs A514C or gidB V77G. The molecular cluster analysis indicated evidence for transmission of resistant strain.

CONCLUSIONS

We observed a substantial rate of MDR-TB in the Madang area of PNG associated with mutations in specific genes. A close monitoring of drug resistance is therefore urgently required, particularly in the presence of drug-resistant M. tuberculosis transmission. In the absence of phenotypic drug susceptibility testing in PNG, molecular assays for drug resistance monitoring would be of advantage.