The transcontinental spread of multidrug-resistant (MDR) tuberculosis is poorly characterized in molecular epidemiologic studies. We used genomic sequencing to understand the establishment and dispersion of MDR Mycobacterium tuberculosis within a group of immigrants to the United States. We used a genomic epidemiology approach to study a genotypically matched (by spoligotype, IS6110 restriction fragment length polymorphism, and mycobacterial interspersed repetitive units-variable number of tandem repeat signature) lineage 2/Beijing MDR strain implicated in an outbreak of tuberculosis among refugees in Thailand and consecutive cases within California. All 46 MDR M. tuberculosis genomes from both Thailand and California were highly related, with a median difference of 10 single-nucleotide polymorphisms (SNPs). The Wat Tham Krabok (WTK) strain is a new sequence type distinguished from all known Beijing strains by 55 SNPs and a genomic deletion (Rv1267c) associated with increased fitness. Sequence data revealed a highly prevalent MDR strain that included several closely related but distinct allelic variants within Thailand, rather than the occurrence of a single outbreak. In California, sequencing data supported multiple independent introductions of WTK with subsequent transmission and reactivation within the state, as well as a potential super spreader with a prolonged infectious period. Twenty-seven drug resistance-conferring mutations and 4 putative compensatory mutations were found within WTK strains. Genomic sequencing has substantial epidemiologic value in both low- and high-burden settings in understanding transmission chains of highly prevalent MDR strains.

The Beijing family of Mycobacterium tuberculosis strains is part of lineage 2 (also known as the East Asian lineage). In clinical studies, we have observed that isolates from the sublineage RD207 of lineage 2 were more readily transmitted among humans. To investigate the basis for this difference, we tested representative strains with the characteristic Beijing spoligotype from four of the five sublineages of lineage 2 in the guinea pig model and subjected these strains to comparative whole-genome sequencing. The results of these studies showed that all of the clinical strains were capable of growing and causing lung pathology in guinea pigs after low-dose aerosol exposure. Differences between the abilities of the four sublineages to grow in the lungs of these animals were not overt, but members of RD207 were significantly more pathogenic, resulting in severe lung damage. The RD207 strains also induced much higher levels of markers associated with regulatory T cells and showed a significant loss of activated T cells in the lungs over the course of the infections. Whole-genome sequencing of the strains revealed mutations specific for RD207 which may explain this difference. Based on these data, we hypothesize that the sublineages of M. tuberculosis are associated with distinct pathological and clinical phenotypes and that these differences influence the transmissibility of particular M. tuberculosis strains in human populations.

The addition of glucose to soil has long been used to study the metabolic activity of microbes in soil; however, the response of the microbial ecophysiology remains poorly characterized. To address this, we sequenced the metagenomes and metatranscriptomes of glucose-amended soil microbial communities in a laboratory incubation.