Search

Mycobacterium brumae is a rapid-growing, non-pathogenic Mycobacterium species, originally isolated from environmental and human samples in Barcelona, Spain. Mycobacterium brumae is not pathogenic and it's in vitro phenotype and immunogenic properties have been well characterized. However, the knowledge of its underlying genetic composition is still incomplete. In this study, we first describe the 4 Mb genome of the M . brumae type strain ATCC 51384 T assembling PacBio reads, and second, we assess the low intraspecies variability by comparing the type strain with Illumina reads from three additional strains. Mycobacterium brumae genome is composed of a circular chromosome with a high GC content of 69.2% and containing 3,791 CDSs, 97 pseudogenes, one prophage and no CRISPR loci. Mycobacterium brumae has shown no pathogenic potential in in vivo experiments, and our genomic analysis confirms its phylogenetic position with other non-pathogenic and rapid growing mycobacteria. Accordingly, we determined the absence of virulence-related genes, such as ESX-1 locus and most PE/PPE genes, among others. Although the immunogenic potential of M . brumae was proved to be as high as Mycobacterium bovis BCG, the only mycobacteria licensed to treat cancer, the genomic content of M . tuberculosis T cell and B cell antigens in M . brumae genome is considerably lower than those antigens present in M . bovis BCG genome. Overall, this work provides relevant genomic data on one of the species of the mycobacterial genus with high therapeutic potential., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Renau-Mínguez, Herrero-Abadía, Ruiz-Rodriguez, Sentandreu, Torrents, Chiner-Oms, Torres-Puente, Comas, Julián and Coscolla.)

The S:A222V point mutation, within the G clade, was characteristic of the 20E (EU1) SARS-CoV-2 variant identified in Spain in early summer 2020. This mutation has since reappeared in the Delta subvariant AY.4.2, raising questions about its specific effect on viral infection. We report combined serological, functional, structural and computational studies characterizing the impact of this mutation. Our results reveal that S:A222V promotes an increased RBD opening and slightly increases ACE2 binding as compared to the parent S:D614G clade. Finally, S:A222V does not reduce sera neutralization capacity, suggesting it does not affect vaccine effectiveness., Competing Interests: NO authors have competing interests.

We have detected two mutations in the spike protein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) at amino acid positions 1163 and 1167 that appeared independently in multiple transmission clusters and different genetic backgrounds. Furthermore, both mutations appeared together in a cluster of 1,627 sequences belonging to clade 20E. This cluster is characterized by 12 additional single nucleotide polymorphisms but no deletions. The available structural information on the S protein in the pre- and postfusion conformations predicts that both mutations confer rigidity, which could potentially decrease viral fitness. Accordingly, we observed reduced infectivity of this spike genotype relative to the ancestral 20E sequence in vitro , and the levels of viral RNA in nasopharyngeal swabs were not significantly higher. Furthermore, the mutations did not impact thermal stability or antibody neutralization by sera from vaccinated individuals but moderately reduce neutralization by convalescent-phase sera from the early stages of the pandemic. Despite multiple successful appearances of the two spike mutations during the first year of SARS-CoV-2 evolution, the genotype with both mutations was displaced upon the expansion of the 20I (Alpha) variant. The midterm fate of the genotype investigated was consistent with the lack of advantage observed in the clinical and experimental data. IMPORTANCE We observed repeated, independent emergence of mutations in the SARS-CoV-2 spike involving amino acids 1163 and 1167, within the HR2 functional motif. Conclusions derived from evolutionary and genomic diversity analysis suggest that the co-occurrence of both mutations might pose an advantage for the virus and therefore a threat to effective control of the epidemic. However, biological characterization, including in vitro experiments and analysis of clinical data, indicated no clear benefit in terms of stability or infectivity. In agreement with this, continuous epidemiological surveillance conducted months after the first observations revealed that both mutations did not successfully outcompete other variants and stopped circulating 9 months after their initial detection. Additionally, we evaluated the potential of both mutations to escape neutralizing antibodies, finding that the presence of these two mutations on their own is not likely to confer antibody escape. Our results provide an example of how newly emerged spike mutations can be assessed to better understand the risk posed by new variants and indicate that some spike mutations confer no clear advantage to the virus despite independently emerging multiple times and are eventually displaced by fitter variants.

The coronavirus disease 2019 (COVID-19) pandemic has affected the world radically since 2020. Spain was one of the European countries with the highest incidence during the first wave. As a part of a consortium to monitor and study the evolution of the epidemic, we sequenced 2,170 samples, diagnosed mostly before lockdown measures. Here, we identified at least 500 introductions from multiple international sources and documented the early rise of two dominant Spanish epidemic clades (SECs), probably amplified by superspreading events. Both SECs were related closely to the initial Asian variants of SARS-CoV-2 and spread widely across Spain. We inferred a substantial reduction in the effective reproductive number of both SECs due to public-health interventions (R e  < 1), also reflected in the replacement of SECs by a new variant over the summer of 2020. In summary, we reveal a notable difference in the initial genetic makeup of SARS-CoV-2 in Spain compared with other European countries and show evidence to support the effectiveness of lockdown measures in controlling virus spread, even for the most successful genetic variants., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Pathogens of the Mycobacterium tuberculosis complex (MTBC) are considered to be monomorphic, with little gene content variation between strains. Nevertheless, several genotypic and phenotypic factors separate strains of the different MTBC lineages (L), especially L5 and L6 (traditionally termed Mycobacterium africanum ) strains, from each other. However, this genome variability and gene content, especially of L5 strains, has not been fully explored and may be important for pathobiology and current approaches for genomic analysis of MTBC strains, including transmission studies. By comparing the genomes of 355 L5 clinical strains (including 3 complete genomes and 352 Illumina whole-genome sequenced isolates) to each other and to H37Rv, we identified multiple genes that were differentially present or absent between H37Rv and L5 strains. Additionally, considerable gene content variability was found across L5 strains, including a split in the L5.3 sub-lineage into L5.3.1 and L5.3.2. These gene content differences had a small knock-on effect on transmission cluster estimation, with clustering rates influenced by the selected reference genome, and with potential overestimation of recent transmission when using H37Rv as the reference genome. We conclude that full capture of the gene diversity, especially high-resolution outbreak analysis, requires a variation of the single H37Rv-centric reference genome mapping approach currently used in most whole-genome sequencing data analysis pipelines. Moreover, the high within-lineage gene content variability suggests that the pan-genome of M. tuberculosis is at least several kilobases larger than previously thought, implying that a concatenated or reference-free genome assembly ( de novo ) approach may be needed for particular questions.

Background: Lineage 1 (L1) and 3 (L3) are two lineages of the Mycobacterium tuberculosis complex (MTBC) causing tuberculosis (TB) in humans. L1 and L3 are prevalent around the rim of the Indian Ocean, the region that accounts for most of the world's new TB cases. Despite their relevance for this region, L1 and L3 remain understudied. Methods: We analyzed 2,938 L1 and 2,030 L3 whole genome sequences originating from 69 countries. We reconstructed the evolutionary history of these two lineages and identified genes under positive selection. Results: We found a strongly asymmetric pattern of migration from South Asia toward neighboring regions, highlighting the historical role of South Asia in the dispersion of L1 and L3. Moreover, we found that several genes were under positive selection, including genes involved in virulence and resistance to antibiotics . For L1 we identified signatures of local adaptation at the esxH locus, a gene coding for a secreted effector that targets the human endosomal sorting complex, and is included in several vaccine candidates. Conclusions: Our study highlights the importance of genetic diversity in the MTBC, and sheds new light on two of the most important MTBC lineages affecting humans., Competing Interests: No competing interests were disclosed., (Copyright: © 2021 Menardo F et al.)

Human tuberculosis (TB) is caused by members of the Mycobacterium tuberculosis complex (MTBC). The MTBC comprises several human-adapted lineages known as M. tuberculosis sensu stricto , as well as two lineages (L5 and L6) traditionally referred to as Mycobacterium africanum . Strains of L5 and L6 are largely limited to West Africa for reasons unknown, and little is known of their genomic diversity, phylogeography and evolution. Here, we analysed the genomes of 350 L5 and 320 L6 strains, isolated from patients from 21 African countries, plus 5 related genomes that had not been classified into any of the known MTBC lineages. Our population genomic and phylogeographical analyses showed that the unclassified genomes belonged to a new group that we propose to name MTBC lineage 9 (L9). While the most likely ancestral distribution of L9 was predicted to be East Africa, the most likely ancestral distribution for both L5 and L6 was the Eastern part of West Africa. Moreover, we found important differences between L5 and L6 strains with respect to their phylogeographical substructure and genetic diversity. Finally, we could not confirm the previous association of drug-resistance markers with lineage and sublineages. Instead, our results indicate that the association of drug resistance with lineage is most likely driven by sample bias or geography. In conclusion, our study sheds new light onto the genomic diversity and evolutionary history of M. africanum , and highlights the need to consider the particularities of each MTBC lineage for understanding the ecology and epidemiology of TB in Africa and globally.

The Mycobacterium tuberculosis complex (MTBC) causes tuberculosis (TB) in humans and various other mammals. The human-adapted members of the MTBC comprise seven phylogenetic lineages that differ in their geographical distribution. There is growing evidence that this phylogeographic diversity modulates the outcome of TB infection and disease. For decades, TB research and development has focused on the two canonical MTBC laboratory strains H37Rv and Erdman, both of which belong to Lineage 4. Relying on only a few laboratory-adapted strains can be misleading as study results might not be directly transferrable to clinical settings where patients are infected with a diverse array of strains, including drug-resistant variants. Here, we argue for the need to expand TB research and development by incorporating the phylogenetic diversity of the MTBC. To facilitate such work, we have assembled a group of 20 genetically well-characterized clinical strains representing the seven known human-adapted MTBC lineages. With the "MTBC clinical strains reference set" we aim to provide a standardized resource for the TB community. We hope it will enable more direct comparisons between studies that explore the physiology of MTBC beyond the laboratory strains used thus far. We anticipate that detailed phenotypic analyses of this reference strain set will increase our understanding of TB biology and assist in the development of new control tools that are broadly effective., Competing Interests: The authors have declared that no competing interests exist.

Purpose: Differentiation of the Mycobacterium tuberculosis complex (MTBc) from non-tuberculous mycobacteria (NTM) is important for tuberculosis diagnosis and is a prerequisite for reliable phenotypic drug-resistance testing. We evaluated the performance of the rapid MPT64 antigen identification test for the detection of Mycobacterium africanum lineage 5 (MAF L5)., Methodology: Smear-positive tuberculosis patients' sputa were included prospectively. Culture was performed on Löwenstein-Jensen medium and, when positive, the MPT64 test and the classical para-nitro benzoic acid susceptibility and heat-labile catalase (PNB/catalase) identification tests were performed. The MPT64 test was repeated 14 days after an initially negative first testing. Direct spoligotyping was performed for MTBc lineage determination., Results: In total, 333 isolates were tested for all of the methods. Three hundred and twenty-two (96.7 %) were pure MTBc, by agreement between spoligotyping and PNB/catalase, and 11 were NTM or a mixture of MTBc/NTM. The MPT64 test conducted on day zero of culture-positivity correctly identified most of the pure MTBc isolates (93.2 %, 300/322), but it failed to detect 24 % of the L5 isolates (18/75) versus 2 % (4/202) of the L4 ones [OR=15.6 (5.3-45.8), P<0.0001], with improved sensitivity for L5 detection on repeat testing after 14 days. The L5-wide non-synonymous single-nucleotide polymorphism in the mpt64 gene may explain the poor performance of the MPT64 test for L5., Conclusion: The MPT64 test has a lower sensitivity for detecting L5 isolates of the MTBc, and can be considered as a first-screening test that should be confirmed by another identification method when it produces negative results in countries with L5. Given the microbiological bias in both the isolation and identification of MAF lineages, diagnostics with high sensitivity for direct testing on clinical material are preferable.

The fast evolution of pathogenic viruses has allowed for the development of phylodynamic approaches that extract information about the epidemiological characteristics of viral genomes. Thanks to advances in whole genome sequencing, they can be applied to slowly evolving bacterial pathogens like Mycobacterium tuberculosis. In this study, we investigate and compare the epidemiological dynamics underlying two M. tuberculosis outbreaks using phylodynamic methods. Specifically, we (i) test if the outbreak data sets contain enough genetic variation to estimate short-term evolutionary rates and (ii) reconstruct epidemiological parameters such as the effective reproduction number. The first outbreak occurred in the Swiss city of Bern (1987-2012) and was caused by a drug-susceptible strain belonging to the phylogenetic M. tuberculosis Lineage 4. The second outbreak was caused by a multidrug-resistant (MDR) strain of Lineage 2, imported from the Wat Tham Krabok (WTK) refugee camp in Thailand into California. There is little temporal signal in the Bern data set and moderate temporal signal in the WTK data set. Thanks to its high sampling proportion (90%) the Bern outbreak allows robust estimation of epidemiological parameters despite the poor temporal signal. Conversely, there is much uncertainty in the epidemiological estimates concerning the sparsely sampled (9%) WTK outbreak. Our results suggest that both outbreaks peaked around 1990, although they were only recognized as outbreaks in 1993 (Bern) and 2004 (WTK). Furthermore, individuals were infected for a significantly longer period (around 9 years) in the WTK outbreak than in the Bern outbreak (4-5 years). Our work highlights both the limitations and opportunities of phylodynamic analysis of outbreaks involving slowly evolving pathogens: (i) estimation of the evolutionary rate is difficult on outbreak time scales and (ii) a high sampling proportion allows quantification of the age of the outbreak based on the sampling times, and thus allows for robust estimation of epidemiological parameters., (Copyright © 2018 The Authors. Published by Elsevier B.V. All rights reserved.)

Tuberculosis (TB) affects humans and other animals and is caused by bacteria from the Mycobacterium tuberculosis complex (MTBC). Previous studies have shown that there are at least nine members of the MTBC infecting animals other than humans; these have also been referred to as ecotypes. However, the ecology and the evolution of these animal-adapted MTBC ecotypes are poorly understood. Here we screened 12,886 publicly available MTBC genomes and newly sequenced 17 animal-adapted MTBC strains, gathering a total of 529 genomes of animal-adapted MTBC strains. Phylogenomic and comparative analyses confirm that the animal-adapted MTBC members are paraphyletic with some members more closely related to the human-adapted Mycobacterium africanum Lineage 6 than to other animal-adapted strains. Furthermore, we identified four main animal-adapted MTBC clades that might correspond to four main host shifts; two of these clades are hypothesized to reflect independent cattle domestication events. Contrary to what would be expected from an obligate pathogen, MTBC nucleotide diversity was not positively correlated with host phylogenetic distances, suggesting that host tropism in the animal-adapted MTBC seems to be driven by contact rates and demographic aspects of the host population rather by than host relatedness. By combining phylogenomics with ecological data, we propose an evolutionary scenario in which the ancestor of Lineage 6 and all animal-adapted MTBC ecotypes was a generalist pathogen that subsequently adapted to different host species. This study provides a new phylogenetic framework to better understand the evolution of the different ecotypes of the MTBC and guide future work aimed at elucidating the molecular mechanisms underlying host range.

Background: Large sequence datasets are difficult to visualize and handle. Additionally, they often do not represent a random subset of the natural diversity, but the result of uncoordinated and convenience sampling. Consequently, they can suffer from redundancy and sampling biases., Results: Here we present Treemmer, a simple tool to evaluate the redundancy of phylogenetic trees and reduce their complexity by eliminating leaves that contribute the least to the tree diversity., Conclusions: Treemmer can reduce the size of datasets with different phylogenetic structures and levels of redundancy while maintaining a sub-sample that is representative of the original diversity. Additionally, it is possible to fine-tune the behavior of Treemmer including any kind of meta-information, making Treemmer particularly useful for empirical studies.

Zoonotic tuberculosis is a risk for human health, especially when animals are in close contact with humans. Mycobacterium tuberculosis was cultured from several organs, including lung tissue and gastric mucosa, of three captive elephants euthanized in a Swiss zoo. The elephants presented weight loss, weakness and exercise intolerance. Molecular characterization of the M. tuberculosis isolates by spoligotyping revealed an identical profile, suggesting a single source of infection. Multilocus variable-number of tandem-repeat analysis (MLVA) elucidated two divergent populations of bacteria and mixed infection in one elephant, suggesting either different transmission chains or prolonged infection over time. A total of eight M. tuberculosis isolates were subjected to whole-genome sequence (WGS) analysis, confirming a single source of infection and indicating the route of transmission between the three animals. Our findings also show that the methods currently used for epidemiological investigations of M. tuberculosis infections should be carefully applied on isolates from elephants. Moreover the importance of multiple sampling and analysis of within-host mycobacterial clonal populations for investigations of transmission is demonstrated.

Mycobacterium africanum lineage (L) 6 is an important pathogen in West Africa, causing up to 40% of pulmonary tuberculosis (TB). The biology underlying the clinical differences between M. africanum and M. tuberculosis sensu stricto remains poorly understood. We performed ex vivo expression of 2179 genes of the most geographically dispersed cause of human TB, M. tuberculosis L4 and the geographically restricted, M. africanum L6 directly from sputa of 11 HIV-negative TB patients from The Gambia who had not started treatment. The DosR regulon was the most significantly decreased category in L6 relative to L4. Further, we identified nonsynonymous mutations in major DosR regulon genes of 44 L6 genomes of TB patients from The Gambia and Ghana. Using Lebek's test, we assessed differences in oxygen requirements for growth. L4 grew only at the aerobic surface while L6 grew throughout the medium. In the host, the DosR regulon is critical for M. tuberculosis in adaptation to oxygen limitation. However, M. africanum L6 appears to have adapted to growth under hypoxic conditions or to different biological niches. The observed under expression of DosR in L6 fits with the genomic changes in DosR genes, microaerobic growth and the association with extrapulmonary disease., (Copyright © 2017 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Antibiotic resistance is one of the biggest threats to human health globally. Alarmingly, multidrug-resistant and extensively drug-resistant Mycobacterium tuberculosis have now spread worldwide. Some key antituberculosis antibiotics are prodrugs, for which resistance mechanisms are mainly driven by mutations in the bacterial enzymatic pathway required for their bioactivation. We have developed drug-like molecules that activate a cryptic alternative bioactivation pathway of ethionamide in M. tuberculosis , circumventing the classic activation pathway in which resistance mutations have now been observed. The first-of-its-kind molecule, named SMARt-420 (Small Molecule Aborting Resistance), not only fully reverses ethionamide-acquired resistance and clears ethionamide-resistant infection in mice, it also increases the basal sensitivity of bacteria to ethionamide., (Copyright © 2017, American Association for the Advancement of Science.)

Tuberculosis is caused by different groups of bacteria belonging to the Mycobacterium tuberculosis complex (MTBC). The combined action of human factors, environmental conditions and bacterial virulence determine the extent and form of human disease. MTBC virulence is a composite of different clinical phenotypes such as transmission rate and disease severity among others. Clinical phenotypes are also influenced by cellular and immunological phenotypes. MTBC phenotypes are determined by the genotype, therefore finding genotypes responsible for clinical phenotypes would allow discovering MTBC virulence factors. Different MTBC strains display different cellular and clinical phenotypes. Strains from Lineage 5 and Lineage 6 are metabolically different, grow slower, and are less virulent. Also, at least certain groups of Lineage 2 and Lineage 4 strains are more virulent in terms of disease severity and human-to-human transmission. Because phenotypic differences are ultimately caused by genotypic differences, different genomic loci have been related to various cellular and clinical phenotypes. However, defining the impact of specific bacterial genomic loci on virulence when other bacterial determinants, human and environmental factors are also impacting the phenotype would contribute to a better knowledge of tuberculosis virulence and ultimately benefit tuberculosis control.

Molecular epidemiological assessments, drug treatment optimization, and development of immunological interventions all depend on understanding pathogen adaptation and genetic variation, which differ for specific pathogens. Mycobacterium tuberculosis is an exceptionally successful human pathogen, yet beyond knowledge that this bacterium has low overall genomic variation but acquires drug resistance mutations, little is known of the factors that drive its population genomic characteristics. Here, we compared the genetic diversity of the bacteria that established infection to the bacterial populations obtained from infected tissues during murine M. tuberculosis pulmonary infection and human disseminated M. bovis BCG infection. We found that new mutations accumulate during in vitro culture, but that in vivo, purifying selection against new mutations dominates, indicating that M. tuberculosis follows a dominant lineage model of evolution. Comparing bacterial populations passaged in T cell-deficient and immunocompetent mice, we found that the presence of T cells is associated with an increase in the diversity of the M. tuberculosis genome. Together, our findings put M. tuberculosis genetic evolution in a new perspective and clarify the impact of T cells on sequence diversity of M. tuberculosis., Competing Interests: The authors have declared that no competing interests exist.

Immigrants from regions with a high incidence of tuberculosis (TB) are a risk group for TB in low-incidence countries such as Switzerland. In a previous analysis of a nationwide collection of 520 Mycobacterium tuberculosis isolates from 2000 to 2008, we identified 35 clusters comprising 90 patients based on standard genotyping (24-locus mycobacterial interspersed repetitive-unit-variable-number tandem-repeat [MIRU-VNTR] typing and spoligotyping). Here, we used whole-genome sequencing (WGS) to revisit these transmission clusters. Genome-based transmission clusters were defined as isolate pairs separated by ≤12 single nucleotide polymorphisms (SNPs). WGS confirmed 17/35 (49%) MIRU-VNTR typing clusters; the other 18 clusters contained pairs separated by >12 SNPs. Most transmission clusters (3/4) of Swiss-born patients were confirmed by WGS, as opposed to 25% (4/16) of the clusters involving only foreign-born patients. The overall clustering proportion was 17% (90 patients; 95% confidence interval [CI], 14 to 21%) by standard genotyping but only 8% (43 patients; 95% CI, 6 to 11%) by WGS. The clustering proportion was 17% (67/401; 95% CI, 13 to 21%) by standard genotyping and 7% (26/401; 95% CI, 4 to 9%) by WGS among foreign-born patients and 19% (23/119; 95% CI, 13 to 28%) and 14% (17/119; 95% CI, 9 to 22%), respectively, among Swiss-born patients. Using weighted logistic regression, we found weak evidence of an association between birth origin and transmission (adjusted odds ratio of 2.2 and 95% CI of 0.9 to 5.5 comparing Swiss-born patients to others). In conclusion, standard genotyping overestimated recent TB transmission in Switzerland compared to WGS, particularly among immigrants from regions with a high TB incidence, where genetically closely related strains often predominate. We recommend the use of WGS to identify transmission clusters in settings with a low incidence of TB., (Copyright © 2016, American Society for Microbiology. All Rights Reserved.)

Pathogens that evade adaptive immunity typically exhibit antigenic variation. By contrast, it appears that although the chronic human tuberculosis (TB)-causing pathogen Mycobacterium tuberculosis needs to counter host T cell responses, its T cell epitopes are hyperconserved. Here we present an extensive analysis of the T cell epitopes of M. tuberculosis. We combined population genomics with experimental immunology to determine the number and identity of T cell epitope sequence variants in 216 phylogenetically diverse strains of M. tuberculosis. Antigen conservation is indeed a hallmark of M. tuberculosis. However, our analysis revealed a set of seven variable antigens that were immunogenic in subjects with active TB. These findings suggest that M. tuberculosis uses mechanisms other than antigenic variation to evade T cells. T cell epitopes that exhibit sequence variation may not be subject to the same evasion mechanisms, and hence vaccines that include such variable epitopes may be more efficacious.

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 Author 2015. Published by Oxford University Press on behalf of the Infectious Diseases Society of America. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

The DNA replication machinery is an important target for antibiotic development in increasingly drug-resistant bacteria, including Mycobacterium tuberculosis. Although blocking DNA replication leads to cell death, disrupting the processes used to ensure replication fidelity can accelerate mutation and the evolution of drug resistance. In Escherichia coli, the proofreading subunit of the replisome, the ɛ exonuclease, is essential for high-fidelity DNA replication; however, we find that the corresponding subunit is completely dispensable in M. tuberculosis. Rather, the mycobacterial replicative polymerase DnaE1 itself encodes an editing function that proofreads DNA replication, mediated by an intrinsic 3'-5' exonuclease activity within its PHP domain. Inactivation of the DnaE1 PHP domain increases the mutation rate by more than 3,000-fold. Moreover, phylogenetic analysis of DNA replication proofreading in the bacterial kingdom suggests that E. coli is a phylogenetic outlier and that PHP domain-mediated proofreading is widely conserved and indeed may be the ancestral prokaryotic proofreader.

Although small molecules shed from pathogens are widely used to diagnose infection, such tests have not been widely implemented for tuberculosis. Here we show that the recently identified compound, 1-tuberculosinyladenosine (1-TbAd), accumulates to comprise >1% of all Mycobacterium tuberculosis lipid. In vitro and in vivo, two isomers of TbAd were detected that might serve as infection markers. Using mass spectrometry and nuclear magnetic resonance, we established the structure of the previously unknown molecule, N(6)-tuberculosinyladenosine (N(6)-TbAd). Its biosynthesis involves enzymatic production of 1-TbAd by Rv3378c followed by conversion to N(6)-TbAd via the Dimroth rearrangement. Intact biosynthetic genes are observed only within M. tuberculosis complex bacteria, and TbAd was not detected among other medically important pathogens, environmental bacteria, and vaccine strains. With no substantially similar known molecules in nature, the discovery and in vivo detection of two abundant terpene nucleosides support their development as specific diagnostic markers of tuberculosis., (Copyright © 2015 Elsevier Ltd. All rights reserved.)

Background: Whole-genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single-nucleotide polymorphism (SNP) typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland., Methods: On the basis of genome sequences of 3 historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1642 patient isolates and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected., Results: We identified 68 patients associated with the outbreak strain. Most received a tuberculosis diagnosis in 1991-1995, but cases were observed until 2011. Two thirds were homeless and/or substance abusers. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into 3 subclusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs., Conclusions: Strain-specific SNP genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real time and at high resolution., (© The Author 2014. Published by Oxford University Press on behalf of the Infectious Diseases Society of America. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Mycolic acids (MAs) are α-alkyl, β-hydroxy long-chain fatty acids found in abundance in the cell envelope of the Mycobacterium tuberculosis complex (MTBC). MAs form an efficient permeability barrier, modulate host innate immune responses, and are the targets of several anti-tuberculosis drugs. Using mass spectrometry, we measured the relative abundance of 80 MA species across 36 clinical isolates of MTBC covering four major phylogenetic lineages. We found significant variations in the MA patterns between different MTBC strains and lineages. MA patterns of "ancient" lineages contrasted those from "modern" lineages, with a lower representation of alpha-mycolates among Lineage 6 strains and an inversion of the methoxy: keto-mycolates ratio in Lineage 1 strains. By interrogating the whole genome sequences of these MTBC strains, we identified relevant single-nucleotide polymorphisms that may sustain the lineage-specific MA patterns. Our results show that the strain genetic background influences MA metabolism and suggests that strain diversity should be considered in the development of new anti-tuberculosis drugs that target MA synthesis., (© 2014 The Authors. MicrobiologyOpen published by John Wiley & Sons Ltd.)

The causative agent of human tuberculosis, Mycobacterium tuberculosis complex (MTBC), comprises seven phylogenetically distinct lineages associated with different geographical regions. Here we review the latest findings on the nature and amount of genomic diversity within and between MTBC lineages. We then review recent evidence for the effect of this genomic diversity on mycobacterial phenotypes measured experimentally and in clinical settings. We conclude that overall, the most geographically widespread Lineage 2 (includes Beijing) and Lineage 4 (also known as Euro-American) are more virulent than other lineages that are more geographically restricted. This increased virulence is associated with delayed or reduced pro-inflammatory host immune responses, greater severity of disease, and enhanced transmission. Future work should focus on the interaction between MTBC and human genetic diversity, as well as on the environmental factors that modulate these interactions., (Copyright © 2014 The Authors. Published by Elsevier Ltd.. All rights reserved.)

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.

Background: High-throughput DNA sequencing produces vast amounts of data, with millions of short reads that usually have to be mapped to a reference genome or newly assembled. Both reference-based mapping and de novo assembly are computationally intensive, generating large intermediary data files, and thus require bioinformatics skills that are often lacking in the laboratories producing the data. Moreover, many research and practical applications in microbiology require only a small fraction of the whole genome data., Results: We developed KvarQ, a new tool that directly scans fastq files of bacterial genome sequences for known variants, such as single nucleotide polymorphisms (SNP), bypassing the need of mapping all sequencing reads to a reference genome and de novo assembly. Instead, KvarQ loads "testsuites" that define specific SNPs or short regions of interest in a reference genome, and directly synthesizes the relevant results based on the occurrence of these markers in the fastq files. KvarQ has a versatile command line interface and a graphical user interface. KvarQ currently ships with two "testsuites" for Mycobacterium tuberculosis, but new "testsuites" for other organisms can easily be created and distributed. In this article, we demonstrate how KvarQ can be used to successfully detect all main drug resistance mutations and phylogenetic markers in 880 bacterial whole genome sequences. The average scanning time per genome sequence was two minutes. The variant calls of a subset of these genomes were validated with a standard bioinformatics pipeline and revealed >99% congruency., Conclusion: KvarQ is a user-friendly tool that directly extracts relevant information from fastq files. This enables researchers and laboratory technicians with limited bioinformatics expertise to scan and analyze raw sequencing data in a matter of minutes. KvarQ is open-source, and pre-compiled packages with a graphical user interface are available at http://www.swisstph.ch/kvarq.

Tuberculosis caused 20% of all human deaths in the Western world between the seventeenth and nineteenth centuries and remains a cause of high mortality in developing countries. In analogy to other crowd diseases, the origin of human tuberculosis has been associated with the Neolithic Demographic Transition, but recent studies point to a much earlier origin. We analyzed the whole genomes of 259 M. tuberculosis complex (MTBC) strains and used this data set to characterize global diversity and to reconstruct the evolutionary history of this pathogen. Coalescent analyses indicate that MTBC emerged about 70,000 years ago, accompanied migrations of anatomically modern humans out of Africa and expanded as a consequence of increases in human population density during the Neolithic period. This long coevolutionary history is consistent with MTBC displaying characteristics indicative of adaptation to both low and high host densities.