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.

Evaluation of Mycobacterium tuberculosis typing methods in a 4-year study in Schleswig-Holstein, Northern Germany.

In order to evaluate the discriminatory power of different methods for genotyping of Mycobacterium tuberculosis complex (MTBC) isolates, we compared the performance of (i) IS6110 DNA fingerprint typing, (ii) spoligotyping, and (iii) 24-loci mycobacterial interspersed repetitive units-variable number of tandem repeats (MIRU-VNTR) typing in a long-term study on the epidemiology of tuberculosis (TB) in Schleswig-Holstein, the northernmost federal state of Germany. In total, we analyzed 277 MTBC isolates collected from patients between the years 2006 and 2010. The collection comprised a broad spectrum of 13 different genotypes, among which strains of the Haarlem genotype (31%) were most prominent, followed by strains belonging to the Delhi and Beijing lineages (7% and 6%, respectively). On the basis of IS6110 restriction fragment length polymorphism (RFLP) and spoligotyping analyses, 211 isolates had unique patterns (76%) and 66 isolates (24%) were in 20 clusters. MIRU-VNTR combined with spoligotyping analyses revealed 202 isolates with unique patterns (73%) and 75 isolates in 18 clusters (27%). Overall, there was 93.1% concordance between the typing results obtained; 198 strains were identified as unique, and 60 isolates were clustered by both typing combinations (including all 31 isolates with confirmed epidemiological links). Of the remaining 19 isolates with discrepant results, 15 were falsely clustered by MIRU-VNTR (six Beijing genotype strains) and four were clustered by IS6110 RFLP (low IS6110 copy number) only. In conclusion, in the study population investigated, a minority of isolates, especially of the Beijing genotype, clustered by standard 24-loci MIRU-VNTR and without an obvious epidemiological link may require second-line typing by IS6110 RFLP or hypervariable MIRU-VNTR loci.