Whole genome sequencing (WGS) of Mycobacterium tuberculosis has rapidly evolved from a research tool to a clinical application for the diagnosis and management of tuberculosis and in public health surveillance. This evolution has been facilitated by the dramatic drop in costs, advances in technology, and concerted efforts to translate sequencing data into actionable information. There is however a risk that, in the absence of a consensus and international standards, the widespread use of WGS technology may result in data and processes that lack harmonisation, comparability and validation. In this review, we outline the current landscape of WGS pipelines and applications and set out best practices for M. tuberculosis WGS, including standards for bioinformatics pipelines, curated repository of resistance-causing variants, phylogenetic analyses, quality control processes, and standardised reporting. 1. Introduction Mycobacterium tuberculosis complex (Mtbc) pathogens are collectively the top infectious disease killer globally, causing 10 million new tuberculosis (TB) cases annually 1. Increasingly, 95 new TB cases are already resistant to rifampicin and isoniazid (termed multidrug resistance; 96 MDR-TB), the key first line drugs 1. Tackling the spread and drug resistance burden of this pathogen requires concerted global effort in prevention, diagnosis, treatment and surveillance.

Mycobacterium tuberculosis requires a large number of secreted and exported proteins for its virulence, immune modulation and nutrient uptake. Most of these proteins are transported by the different type VII secretion systems. The most recently evolved type VII secretion system, ESX-5, secretes dozens of substrates belonging to the PE and PPE families, which are named for conserved proline and glutamic acid residues close to the amino terminus. However, the role of these proteins remains largely elusive . Here, we show that mutations of ppe38 completely block the secretion of two large subsets of ESX-5 substrates, that is, PPE-MPTR and PE_PGRS, together comprising>80 proteins. Importantly, hypervirulent clinical M. tuberculosis strains of the Beijing lineage have such a mutation and a concomitant loss of secretion . Restoration of PPE38-dependent secretion partially reverted the hypervirulence phenotype of a Beijing strain, and deletion of ppe38 in moderately virulent M. tuberculosis increased virulence. This indicates that these ESX-5 substrates have an important role in virulence attenuation. Phylogenetic analysis revealed that deletion of ppe38 occurred at the branching point of the 'modern' Beijing sublineage and is shared by Beijing outbreak strains worldwide, suggesting that this deletion may have contributed to their success and global distribution.

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.

Mycobacterium africanum consists of Lineages L5 and L6 of the Mycobacterium tuberculosis complex (MTBC) and causes human tuberculosis in specific regions of Western Africa, but is generally not transmitted in other parts of the world. Since M. africanum is evolutionarily closely placed between the globally dispersed Mycobacterium tuberculosis and animal-adapted MTBC-members, these lineages provide valuable insight into M. tuberculosis evolution. Here, we have collected 15 M. africanum L5 strains isolated in France over 4 decades. Illumina sequencing and phylogenomic analysis revealed a previously underappreciated diversity within L5, which consists of distinct sublineages. L5 strains caused relatively high levels of extrapulmonary tuberculosis and included multi- and extensively drug-resistant isolates, especially in the newly defined sublineage L5.2. The specific L5 sublineages also exhibit distinct phenotypic characteristics related to in vitro growth, protein secretion and in vivo immunogenicity. In particular, we identified a PE_PGRS and PPE-MPTR secretion defect specific for sublineage L5.2, which was independent of PPE38. Furthermore, L5 isolates were able to efficiently secrete and induce immune responses against ESX-1 substrates contrary to previous predictions. These phenotypes of Type VII protein secretion and immunogenicity provide valuable information to better link genome sequences to phenotypic traits and thereby understand the evolution of the MTBC.

BackgroundWhole genome sequencing has revolutionised the interrogation of mycobacterial genomes. Recent studies have reported conflicting findings on the genomic stability of Mycobacterium tuberculosis during the evolution of drug resistance. In an age where whole genome sequencing is increasingly relied upon for defining the structure of bacterial genomes, it is important to investigate the reliability of next generation sequencing to identify clonal variants present in a minor percentage of the population. This study aimed to define a reliable cut-off for identification of low frequency sequence variants and to subsequently investigate genetic heterogeneity and the evolution of drug resistance in M. tuberculosis. MethodsGenomic DNA was isolated from single colonies from 14 rifampicin mono-resistant M. tuberculosis isolates, as well as the primary cultures and follow up MDR cultures from two of these patients. The whole genomes of the M. tuberculosis isolates were sequenced using either the Illumina MiSeq or Illumina HiSeq platforms. Sequences were analysed with an in-house pipeline.ResultsUsing next-generation sequencing in combination with Sanger sequencing and statistical analysis we defined a read frequency cut-off of 30 % to identify low frequency M. tuberculosis variants with high confidence. Using this cut-off we demonstrated a high rate of genetic diversity between single colonies isolated from one population, showing that by using the current sequencing technology, single colonies are not a true reflection of the genetic diversity within a whole population and vice versa. We further showed that numerous heterogeneous variants emerge and then disappear during the evolution of isoniazid resistance within individual patients. Our findings allowed us to formulate a model for the selective bottleneck which occurs during the course of infection, acting as a genomic purification event.ConclusionsOur study demonstrated true levels of genetic diversity within an M. tuberculosis population and showed that genetic diversity may be re-defined when a selective pressure, such as drug exposure, is imposed on M. tuberculosis populations during the course of infection. This suggests that the genome of M. tuberculosis is more dynamic than previously thought, suggesting preparedness to respond to a changing environment.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2067-2) contains supplementary material, which is available to authorized users.

No abstract

No abstract

Although currently available model organisms such as Mycobacterium smegmatis and Mycobacterium bovis Bacillus Calmette-Guérin (BCG) have significantly contributed to our understanding of tuberculosis (TB) biology, these models have limitations such as differences in genome size, growth rates and virulence. However, attenuated Mycobacterium tuberculosis strains may provide more representative, safer models to study M. tuberculosis biology. For example, the M. tuberculosis Δ leuD Δ panCD double auxotroph, has undergone rigorous in vitro and in vivo safety testing. Like other auxotrophic strains, this has subsequently been approved for use in biosafety level (BSL) 2 facilities. Auxotrophic strains have been assessed as models for drug-resistant M. tuberculosis and for studying latent TB. These offer the potential as safe and useful models, but it is important to understand how well these recapitulate salient features of non-attenuated M. tuberculosis. We therefore performed a comprehensive comparison of M. tuberculosis H37Rv and M. tuberculosis Δ leuD Δ panCD . These strains demonstrated similar in vitro and intra-macrophage replication rates, similar responses to anti-TB agents and whole genome sequence conservation. Shotgun proteomics analysis suggested that M. tuberculosis Δ leuD Δ panCD has a heightened stress response that leads to reduced bacterial replication during exposure to acid stress, which has been verified using a dual-fluorescent replication reporter assay. Importantly, infection of human peripheral blood mononuclear cells with the 2 strains elicited comparable cytokine production, demonstrating the suitability of M. tuberculosis Δ leuD Δ panCD for immunological assays. We provide comprehensive evidence to support the judicious use of M. tuberculosis Δ leuD Δ panCD as a safe and suitable model organism for M. tuberculosis research, without the need for a BSL3 facility.