Widespread pyrazinamide-resistant Mycobacterium tuberculosis family in a low-incidence setting

Diagnostic Accuracy of Pyrazinamide Susceptibility Testing in Mycobacterium tuberculosis: A Systematic Review with Meta-Analysis

Introduction: Pyrazinamide (PZA) susceptibility testing plays a critical role in determining the appropriate treatment regimens for multidrug-resistant tuberculosis. We conducted a systematic review and meta-analysis to evaluate the diagnostic accuracy of sequencing PZA susceptibility tests against culture-based susceptibility testing methods as the reference standard. Methods: We searched the MEDLINE/PubMed, Embase, and Web of Science databases for the relevant records. The QUADAS-2 tool was used to assess the quality of the studies. Diagnostic accuracy measures (i.e., sensitivity and specificity) were pooled with a random-effects model. All statistical analyses were performed with Meta-DiSc (version 1.4, Cochrane Colloquium, Barcelona, Spain), STATA (version 14, Stata Corporation, College Station, TX), and RevMan (version 5.3, The Nordic Cochrane Centre, the Cochrane Collaboration, Copenhagen, Denmark) software. Results: A total of 72 articles, published between 2000 and 2019, comprising data for 8,701 isolates of Mycobacterium tuberculosis were included in the final analysis. The pooled sensitivity and specificity of the PZA sequencing test against all reference tests (the combination of BACTEC mycobacteria growth indicator tube 960 (MGIT 960), BACTEC 460, and proportion method) were 87% (95% CI: 85-88) and 94.7% (95% CI: 94-95). The positive likelihood ratio, negative likelihood ratio, diagnostic odds ratio, and the area under the curve estimates were found to be 12.0 (95% CI: 9.0-16.0), 0.17 (95% CI: 0.13-0.21), 106 (95% CI: 71-158), and 96%, respectively. Deek's test result indicated a low likelihood for publication bias (p = 0.01). Conclusions: Our analysis indicated that PZA sequencing may be used in combination with conventional tests due to the advantage of the time to result and in scenarios where culture tests are not feasible. Further work to improve molecular tests would benefit from the availability of standardized reference standards and improvements to the methodology.

Effects of Host, Sample, and in vitro Culture on Genomic Diversity of Pathogenic Mycobacteria

Mycobacterium tuberculosis (M. tb), an obligate human pathogen and the etiological agent of tuberculosis (TB), remains a major threat to global public health. Comparative genomics has been invaluable for monitoring the emergence and spread of TB and for gaining insight into adaptation of M. tb. Most genomic studies of M. tb are based on single bacterial isolates that have been cultured for several weeks in vitro. However, in its natural human host, M. tb comprises complex, in some cases massive bacterial populations that diversify over the course of infection and cannot be wholly represented by a single genome. Recently, enrichment via hybridization capture has been used as a rapid diagnostic tool for TB, circumventing culturing protocols and enabling the recovery of M. tb genomes directly from sputum. This method has further applicability to the study of M. tb adaptation, as it enables a higher resolution and more direct analysis of M. tb genetic diversity within hosts with TB. Here we analyzed genomic material from M. tb and Mycobacterium bovis populations captured directly from sputum and from cultured samples using metagenomic and Pool-Seq approaches. We identified effects of sampling, patient, and sample type on bacterial genetic diversity. Bacterial genetic diversity was more variable and on average higher in sputum than in culture samples, suggesting that manipulation in the laboratory reshapes the bacterial population. Using outlier analyses, we identified candidate bacterial genetic loci mediating adaptation to these distinct environments. The study of M. tb in its natural human host is a powerful tool for illuminating host pathogen interactions and understanding the bacterial genetic underpinnings of virulence.

Methods used in the spatial analysis of tuberculosis epidemiology: a systematic review

BACKGROUND

Tuberculosis (TB) transmission often occurs within a household or community, leading to heterogeneous spatial patterns. However, apparent spatial clustering of TB could reflect ongoing transmission or co-location of risk factors and can vary considerably depending on the type of data available, the analysis methods employed and the dynamics of the underlying population. Thus, we aimed to review methodological approaches used in the spatial analysis of TB burden.

METHODS

We conducted a systematic literature search of spatial studies of TB published in English using Medline, Embase, PsycInfo, Scopus and Web of Science databases with no date restriction from inception to 15 February 2017. The protocol for this systematic review was prospectively registered with PROSPERO ( CRD42016036655 ).

RESULTS

We identified 168 eligible studies with spatial methods used to describe the spatial distribution (n = 154), spatial clusters (n = 73), predictors of spatial patterns (n = 64), the role of congregate settings (n = 3) and the household (n = 2) on TB transmission. Molecular techniques combined with geospatial methods were used by 25 studies to compare the role of transmission to reactivation as a driver of TB spatial distribution, finding that geospatial hotspots are not necessarily areas of recent transmission. Almost all studies used notification data for spatial analysis (161 of 168), although none accounted for undetected cases. The most common data visualisation technique was notification rate mapping, and the use of smoothing techniques was uncommon. Spatial clusters were identified using a range of methods, with the most commonly employed being Kulldorff's spatial scan statistic followed by local Moran's I and Getis and Ord's local Gi(d) tests. In the 11 papers that compared two such methods using a single dataset, the clustering patterns identified were often inconsistent. Classical regression models that did not account for spatial dependence were commonly used to predict spatial TB risk. In all included studies, TB showed a heterogeneous spatial pattern at each geographic resolution level examined.

CONCLUSIONS

A range of spatial analysis methodologies has been employed in divergent contexts, with all studies demonstrating significant heterogeneity in spatial TB distribution. Future studies are needed to define the optimal method for each context and should account for unreported cases when using notification data where possible. Future studies combining genotypic and geospatial techniques with epidemiologically linked cases have the potential to provide further insights and improve TB control.

Global expansion of Mycobacterium tuberculosis lineage 4 shaped by colonial migration and local adaptation

On the basis of population genomic and phylogeographic analyses of 1669 Mycobacterium tuberculosis lineage 4 (L4) genomes, we find that dispersal of L4 has been completely dominated by historical migrations out of Europe. We demonstrate an intimate temporal relationship between European colonial expansion into Africa and the Americas and the spread of L4 tuberculosis (TB). Markedly, in the age of antibiotics, mutations conferring antimicrobial resistance overwhelmingly emerged locally (at the level of nations), with minimal cross-border transmission of resistance. The latter finding was found to reflect the relatively recent emergence of these mutations, as a similar degree of local restriction was observed for susceptible variants emerging on comparable time scales. The restricted international transmission of drug-resistant TB suggests that containment efforts at the level of individual countries could be successful.

Phenotypic and genotypic drug sensitivity of Mycobacterium tuberculosis complex isolated from South Omo Zone, Southern Ethiopia

Background

Knowledge of drug-sensitivity patterns of Mycobacterium tuberculosis complex (MTBC) strains isolated from patients is an important aspect of TB control strategy. This study was conducted to evaluate the drug sensitivity of MTBC isolates in South Omo, southern Ethiopia.

Materials and methods

A total of 161 MTBC isolates (153 from new cases and eight re-treatment TB cases) were isolated using Lowenstein Jensen medium of which 126 isolates were able to be tested for drug sensitivity by BACTEC™MGIT™ 960 system, while all the 161 isolates were tested by GenoType^® MTBDRplus VER 2.0. Descriptive statistics and logistic regression were used to express and present results.

Results

On the basis of MGIT 960 system, the prevalence of mono-resistance was 9.2% (11/119) in the new cases, although neither poly-resistance nor multidrug resistance (MDR) was recorded in these cases. On the basis of GenoType MTBDRplus assay, two of the 153 isolates (1.3%) of the new cases were mono-resistant for rifampicin (RIF) and one of these isolates had known rpoB gene mutation (H526D). One of the eight (12.5%) isolates obtained from the re-treatment cases was MDR with rpoB gene mutation (D516V) and katG gene mutation (S315T2). Taking MGIT 960 system as a gold standard, the sensitivities of the MTBDRplus assay were 33.3%, 100% and 100% for detection of resistance to isoniazid, RIF and MDR, respectively. On the other hand, its specificities were 99.2%, 100% and 100% for detection of resistance to RIF, isoniazid and MDR, respectively.

Conclusion

The magnitude of drug resistance was relatively low in the new TB cases of South Omo as compared to the reports from the other regions of the country. This is encouraging and hence the TB Control Program in the Zone should strengthen its program so that the emergence of drug resistance is inhibited.

Signatures of Selection at Drug Resistance Loci in Mycobacterium tuberculosis

Tuberculosis (TB) is the leading cause of death by an infectious disease, and global TB control efforts are increasingly threatened by drug resistance in Mycobacterium tuberculosis. Unlike most bacteria, where lateral gene transfer is an important mechanism of resistance acquisition, resistant M. tuberculosis arises solely by de novo chromosomal mutation. Using whole-genome sequencing data from two natural populations of M. tuberculosis, we characterized the population genetics of known drug resistance loci using measures of diversity, population differentiation, and convergent evolution. We found resistant subpopulations to be less diverse than susceptible subpopulations, consistent with ongoing transmission of resistant M. tuberculosis. A subset of resistance genes ("sloppy targets") were characterized by high diversity and multiple rare variants; we posit that a large genetic target for resistance and relaxation of purifying selection contribute to high diversity at these loci. For "tight targets" of selection, the path to resistance appeared narrower, evidenced by single favored mutations that arose numerous times in the phylogeny and segregated at markedly different frequencies in resistant and susceptible subpopulations. These results suggest that diverse genetic architectures underlie drug resistance in M. tuberculosis and that combined approaches are needed to identify causal mutations. Extrapolating from patterns observed for well-characterized genes, we identified novel candidate variants involved in resistance. The approach outlined here can be extended to identify resistance variants for new drugs, to investigate the genetic architecture of resistance, and when phenotypic data are available, to find candidate genetic loci underlying other positively selected traits in clonal bacteria. IMPORTANCEMycobacterium tuberculosis, the causative agent of tuberculosis (TB), is a significant burden on global health. Antibiotic treatment imposes strong selective pressure on M. tuberculosis populations. Identifying the mutations that cause drug resistance in M. tuberculosis is important for guiding TB treatment and halting the spread of drug resistance. Whole-genome sequencing (WGS) of M. tuberculosis isolates can be used to identify novel mutations mediating drug resistance and to predict resistance patterns faster than traditional methods of drug susceptibility testing. We have used WGS from natural populations of drug-resistant M. tuberculosis to characterize effects of selection for advantageous mutations on patterns of diversity at genes involved in drug resistance. The methods developed here can be used to identify novel advantageous mutations, including new resistance loci, in M. tuberculosis and other clonal pathogens.

Genotypic characterization of drug resistant Mycobacterium tuberculosis in Quebec, 2002-2012

Background

The increasing emergence of drug-resistant tuberculosis presents a threat to the effective control of tuberculosis (TB). Rapid detection of drug-resistance is more important than ever to address this scourge. The purpose of this study was to genotypically characterize the first-line antitubercular drug-resistant isolates collected over 11 years in Quebec.

Results

The main mutations found in our resistant strains collection (n = 225) include: the S315T substitution in katG (50.2 %), the -15 C/T mutation in the inhA promoter (29 %); the S531L substitution in rpoB (43 %); the deletion 8 bp 446 / + R140S in pncA (72.9 %); the M306I (35.7 %) and M306V (21.4 %) substitutions in embB. Ten of the mutations in katG and 4 mutations identified in pncA were previously undescribed.

Conclusion

Screening of mutations conferring resistance to first-line antituberculous drugs using DNA-sequencing approach seems to be feasible and would drastically shorten the time to determine the resistance profile compared to the proportion method.

Outbreak of pyrazinamide-monoresistant tuberculosis identified using genotype cluster and social media analysis

SETTING

Monoresistance to pyrazinamide (PZA) has infrequently been associated with Mycobacterium tuberculosis.

OBJECTIVE

To report an outbreak of PZA-monoresistant M. tuberculosis in Virginia involving two genotype clusters from December 2004 to August 2010.

RESULTS

Thirty cases were identified involving a predominantly young, US-born population with histories of substance use and incarceration and a large proportion of children aged <15 years (n= 6, 20%); of these, 23 cases (77%) were culture-confirmed as M. tuberculosis complex. DNA fingerprinting and molecular analysis of the PZA resistance gene, pncA, demonstrated a clonal strain that was not M. bovis. Genotypic data provided the initial link between seemingly unrelated cases, and helped reveal a historic genotype cluster of cases from 2004. Further genotype cluster and contact investigation procedures, including the novel use of the social networking website Facebook.com, revealed additional links between the 2004 and 2009 genotype clusters and described an ongoing, extensive outbreak necessitating an enhanced screening and treatment protocol for contacts.

CONCLUSIONS

This outbreak demonstrates how tuberculosis can spread through a young, vulnerable population. The use of genotypic data and the novel incorporation of social media investigations were critical to understanding the settings and context of infectivity.

Gene sequencing for routine verification of pyrazinamide resistance in Mycobacterium tuberculosis: a role for pncA but not rpsA

Pyrazinamide (PZA) is an important component of first-line therapy for the treatment of tuberculosis. Here, we evaluate targeted gene sequencing as a supplement to phenotypic PZA susceptibility testing of Mycobacterium tuberculosis. Routine sequencing of pncA, but not rpsA, is effective for verification of PZA susceptibility results.

How close is close enough? Exploring matching criteria in the estimation of recent transmission of tuberculosis

If Mycobacterium tuberculosis isolates from 2 people have the same genotype, transmission may have occurred between them. Genotyping based on the insertion sequence IS6110 uses identical restriction fragment length polymorphisms ("fingerprints") to infer transmission. However, once transmission has occurred, the genotypes may mutate, resulting in divergent fingerprints. Estimation of the proportion of tuberculosis (TB) cases due to recent transmission includes 3 approaches to determine if genotypes match: exact matching (assumes no fingerprint change); band-addition, band-loss, band-shift matching (ad hoc attempt to account for fingerprint changes); and genetic distance (directly accounts for fingerprint changes). Via simulation study, the authors varied the fingerprint change rate, level of recent transmission, and background genetic heterogeneity and estimated sensitivity, specificity, and bias of the recent transmission index by matching method. For exact matching, specificity was always high, but sensitivity decreased as the change rate increased. For band-addition, band-loss, band-shift matching, specificity decreased as genetic diversity decreased, and sensitivity remained high as the change rate increased. Genetic distance offered a compromise between the 2. Results from this study suggest that interpretation of the recent transmission index and the resulting necessary public health interventions will vary according to how researchers account for spontaneous mutation when estimating transmission from genotyping data.

Pharmacokinetics-pharmacodynamics of pyrazinamide in a novel in vitro model of tuberculosis for sterilizing effect: a paradigm for faster assessment of new antituberculosis drugs

There are currently renewed efforts to develop drugs that could shorten the duration of antituberculosis therapy. This is best achieved by optimizing the sterilizing effect. However, the current pathway for the development of new molecules with the potential to have a sterilizing effect is inefficient. We designed an in vitro pharmacokinetic-pharmacodynamic model in which Mycobacterium tuberculosis replicating slowly at pH 5.8 was exposed to pyrazinamide by use of the concentration-time profiles encountered in patients. The sterilizing effect rates and the time to the emergence of drug resistance were examined. Daily pyrazinamide dosing for 28 days accurately achieved (i) the pyrazinamide pharmacokinetic parameters, (ii) the lack of early bactericidal activity, (iii) a sterilizing effect rate of 0.10 log(10) CFU/ml per day starting on day 6 of therapy, and (iv) a time to the emergence of resistance of the from 2 to 3 weeks of monotherapy encountered in patients with tuberculosis. Next, dose-scheduling studies were performed. The sterilizing effect was linked to the pyrazinamide ratio of the area under the concentration-time curve from 0 to 24 h (AUC(0-24)) to the MIC (r(2) = 0.80 to 0.90), with 90% of the maximal effect being achieved by an AUC(0-24)/MIC of 209.08. Resistance suppression was associated with the percentage of time that the concentration persisted above the MIC (r(2) = 0.73 to 0.91). Monte Carlo simulations of 10,000 patients demonstrated that the currently recommended pyrazinamide doses (15 to 30 mg/kg of body weight/day) achieved the AUC(0-24)/MIC of 209.08 in the epithelial lining fluid of only 15.1 to 53.3% of patients. Doses of >60 mg/kg per day performed better. Our vitro model for the sterilizing effect, together with Monte Carlo simulations, can be used for the faster identification of the clinical doses that are needed to achieve a sterilizing effect and that can then be studied in clinical trials.

Geography and genealogy of the human host harbouring a distinctive drug-resistant strain of tuberculosis

For a strain of Mycobacterium tuberculosis mono-resistant to pyrazinamide (PZA), we report the geographic distribution within Quebec of the 77 cases diagnosed during 1990-2000. Known as the Quebec mutation (or the pncA deletion), the strain is rare in urban areas and showed an unexpected concentration in Mauricie, one of the 16 health districts of the province, with a cluster of 10 cases situated in a rural area of 35-km radius. The cases occurred among people >50 (98%), of French Canadian origins (90%), and are understood to have arisen by reactivation. The rarity in Montreal and smaller cities is explained by the youthfulness of massive postwar migrations. To reach back into the history of settlement, we examined genealogies: 92,429 ancestral marriages for 32 of the 77 PZA-resistant isolates and 226,535 for a set of 85 controls with isolates of more diverse mycobacterial strains. Genealogical analysis showed no salient common ancestor for the cases, and kinship among them was no greater than observed in control samples from the same regions. But it identified an unsuspected geographical region as the site of ancestral concentrations prior to 1840, for both resistant strains and controls. The following scenario is proposed for the resistant strain: endemic in a specific geographical region by 1800, it dispersed with families moving into regions opened to settlement in the 1840s and 1850s, among them Mauricie, where dispersion was intensified by seasonal mobility of labour in logging, milling and marketing timber. In high-incidence areas, it is difficult to distinguish cases of reactivation from recent infections, but the low-incidence context allows us to observe a 200-year trajectory of a distinctive drug-resistant strain of M. tuberculosis.

Contact investigation and genotyping to identify tuberculosis transmission to children

BACKGROUND

Tuberculosis (TB) in young children is an indicator of ongoing community transmission. We examined contact investigations related to pediatric TB, yield for source case identifications and genotypes for relevant Mycobacterium tuberculosis isolates in a low-incidence setting.

METHODS

We reviewed public health data for all patients with TB aged <18 years reported to Montreal authorities during 1996 to 2000. M. tuberculosis isolates from patients of all ages were subjected to IS6110-based genotyping, supplemented by spoligotyping, to compare isolates from children and adults during the same years.

RESULTS

Sixty-six patients aged <18 years were diagnosed with active TB from 1996 to 2000. Mean age was 11.1 years (standard deviation 6.7 years). Twenty-five children (38%) were Canadian-born, all with at least one foreign-born parent. Nineteen children were diagnosed after contact investigations of known adult cases; 8 underwent no contact investigation. For the remaining 39 children, a total of 616 contacts were identified. The median number of contacts per child was 9 (interquartile range, 6-10). Four hundred eighty-one contacts (78%) underwent tuberculin testing; 188 (39%) were reactors and 186 (39%) began treatment of latent TB. Investigations uncovered 4 probable source cases, all involving parents or other relatives. M. tuberculosis genotyping for 38 children identified up to 14 additional possible source cases; in only one was a possible epidemiologic link evident from public health records.

CONCLUSIONS

Among largely foreign-born children with active TB, contact investigations were extensive and often identified latent tuberculosis infection--but rarely source cases. However, genotyping suggested substantial, previously unrecognized transmission to children despite low overall incidence.

Spoligotype signatures in the Mycobacterium tuberculosis complex

Evolution of the direct repeat region in Mycobacterium tuberculosis has created unique spoligotype signatures specifically associated with IS6110-defined strain families. Spoligotyping signatures may enable the analysis of the strain population structure in different settings and will enable the rapid identification of strain families that acquire drug resistance or escape protective immunity in drug and vaccine trials.

Sensitivities and specificities of spoligotyping and mycobacterial interspersed repetitive unit-variable-number tandem repeat typing methods for studying molecular epidemiology of tuberculosis

The development of PCR-based genotyping modalities (spoligotyping and mycobacterial interspersed repetitive unit-variable-number tandem repeat [MIRU-VNTR] typing) offers promise for real-time molecular epidemiological studies of tuberculosis (TB). However, the utility of these methods depends on their capacity to appropriately classify isolates. To determine the operating parameters of spoligotyping and MIRU-VNTR typing, we have compared results generated by these newer tests to the standard typing method, IS6110 restriction fragment length polymorphism, in analyses restricted to high-copy-number IS6110 isolates. Sensitivities of the newer tests were estimated as the percentages of isolates with identical IS6110 fingerprints that had identical spoligotypes and MIRU-VNTR types. The specificities of these tests were estimated as the percentages of isolates with unique IS6110 fingerprints that had unique spoligotypes and MIRU-VNTR types. The sensitivity of MIRU-VNTR typing was 52% (95% confidence interval [CI], 31 to 72%), and the sensitivity of spoligotyping was 83% (95% CI, 63 to 95%). The specificity of MIRU-VNTR typing was 56% (95% CI, 51 to 62%), and the specificity of spoligotyping was 40% (95% CI, 35 to 46%). The proportion of isolates estimated to be due to recent transmission was 4% by identical IS6110 patterns, 19% by near-identical IS6110 patterns, 33% by MIRU-VNTR typing, and 53% by spoligotyping. The low calculated specificities of spoligotyping and MIRU-VNTR typing led to misclassification of cases, inflated estimates of TB transmission, and low positive predictive values, suggesting that these techniques have unsuitable operating parameters for population-based molecular epidemiology studies.

Genomic characterization of an endemic Mycobacterium tuberculosis strain: evolutionary and epidemiologic implications

In a study of 302 Mycobacterium tuberculosis clinical isolates from the low-incidence Canadian-born population of Quebec, we characterized a large endemic strain family by using genomic deletions. The DS6(Quebec) deleted region (11.4 kb) defined a strain family of 143 isolates encompassing two subgroups: one characterized by pyrazinamide (PZA) susceptibility and the other marked by a PZA-monoresistant phenotype. A second deletion (8 bp) in the pncA gene was shared by all 76 isolates with the PZA resistance phenotype, whereas a third DRv0961 deletion (970 bp) defined a further subset of 15 isolates. From their deletion profiles, we derived a most parsimonious evolutionary scenario and compared multiple standard genotyping modalities (using IS6110 restriction fragment length polymorphism [RFLP], spoligotyping, and mycobacterial interspersed repetitive units [MIRU]) across the deletion-based subgroups. The use of a single genotyping modality yielded an unexpectedly high proportion of clustered isolates for a high IS6110 copy strain (27% by IS6110 RFLP, 61% by MIRU, and 77% by spoligotyping). By combining all three modalities, only 14% were genotypically clustered overall, a result more congruent with the epidemiologic profile of reactivation tuberculosis, as suggested by the older age (mean age, 60 years), rural setting, and low proportion of epidemiologic links. These results provide insight into the evolution of genotypes in endemic strains and the potential for false clustering in molecular epidemiologic studies.

Genetic characterization of the Guinea-Bissau family of Mycobacterium tuberculosis complex strains

In a previous study of Mycobacterium tuberculosis complex isolates from Guinea-Bissau in West Africa, we identified a unique group of strains, designated here as the Guinea-Bissau family of strains, which, although genotypically closely related, phenotypically demonstrated a considerable heterogeneity. We conducted here a detailed genotypic analysis of a subset (n = 35) of these isolates. Based on the data obtained, and by comparison of known corresponding genes in mycobacteria outside the M. tuberculosis complex, we propose that the Guinea-Bissau strains belong to a unique branch of the M. tuberculosis complex tree in between classical M. tuberculosis and classical M. bovis. These observations are discussed in their significance in M. tuberculosis complex classification.

Tuberculosis in the Inuit Community of Quebec, Canada

In low-incidence countries targeting tuberculosis (TB) elimination, TB remains a problem of a few high-risk groups. In Canada, Aboriginals, and particularly the Arctic Inuit communities, have witnessed dramatic decreases in TB during the 1960s to 1970s, but rates remain at least 10 to 20 times higher than the national average. We are describing the results of an integrated traditional and molecular epidemiology study of all culture-positive Mycobacterium tuberculosis cases in the Arctic Inuit communities of Quebec from 1990 until 2000. The demographic characteristics of the 46 TB cases included in the study were most notable for a bimodal age distribution (48% under 25 years). Genotyping analysis using multiple modalities (IS6110 restriction fragment length polymorphism, spoligotype, mycobacterial interspersed repetitive units-variable number tandem repeats) showed that 76% (35/46) of TB cases were clustered (six clusters, median size four cases) and estimated that at least 62.5% of TB cases were due to ongoing transmission. By integrating the epidemiologic and genotyping data, we observed that the genotyping clustering results were concordant with recognized epidemiologic links but most notably identified previously unrecognized intervillage transmission. This study demonstrates significant ongoing transmission in a geographically isolated, low-density population. In a resource-rich country such as Canada, these communities illustrate some of the persistent challenges of TB control and elimination.