ABSTRACTThe Beijing family ofMycobacterium tuberculosisstrains 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 ofM. tuberculosisare associated with distinct pathological and clinical phenotypes and that these differences influence the transmissibility of particularM. tuberculosisstrains in human populations.

The organisms in aerosol microenvironments, especially densely populated urban areas, are relevant to maintenance of public health and detection of potential epidemic or biothreat agents. To examine aerosolized microorganisms in this environment, we performed sequencing on the material from an urban aerosol surveillance program. Whole metagenome sequencing was applied to DNA extracted from air filters obtained during periods from each of the four seasons. The composition of bacteria, plants, fungi, invertebrates, and viruses demonstrated distinct temporal shifts. Bacillus thuringiensis serovar kurstaki was detected in samples known to be exposed to aerosolized spores, illustrating the potential utility of this approach for identification of intentionally introduced microbial agents. Together, these data demonstrate the temporally dependent metagenomic complexity of urban aerosols and the potential of genomic analytical techniques for biosurveillance and monitoring of threats to public health.Electronic supplementary materialThe online version of this article (doi:10.1007/s00248-014-0517-z) contains supplementary material, which is available to authorized users.

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.

Summary: PIQA is a quality analysis pipeline designed to examine genomic reads produced by Next Generation Sequencing technology (Illumina G1 Genome Analyzer). A short statistical summary, as well as tile-by-tile and cycle-by-cycle graphical representation of clusters density, quality scores and nucleotide frequencies allow easy identification of various technical problems including defective tiles, mistakes in sample/library preparations and abnormalities in the frequencies of appearance of sequenced genomic reads. PIQA is written in the R statistical programming language and is compatible with bustard, fastq and scarf Illumina G1 Genome Analyzer data formats.Availability: The PIQA pipeline, installation instructions and examples are available at the supplementary web site (http://bioinfo.uh.edu/PIQA).Contact: yfofanov@bioinfo.uh.edu