Search

Abstract Background Culture contamination with environmental bacteria is a major challenge in tuberculosis (TB) laboratories in hot and humid climate zones. We studied the effect of cetylpyridinium chloride (CPC) preservation on culture results and performance of Xpert MTB/RIF. Methods Consecutive sputum samples from microscopy smear-positive TB patients were collected. Two-hundred samples were equally split in two aliquots, one aliquot was treated with CPC and stored at ambient temperature for 7 days. The second aliquot was immediately processed. Samples were decontaminated for 20, 15 or 10 min, and subsequently cultured on Löwenstein-Jensen medium. Furthermore, 50 samples were stored for 7, 14 and 21 days, and 100 CPC-pretreated samples tested by Xpert MTB/RIF. Results CPC pretreated samples showed a higher culture yield compared to non-treated sputum samples across all decontamination times: 94% vs. 73% at 10 min (p = 0.01), 94% vs. 64% at 15 min (p = 0.004), and 90% vs. 52% at 20 min (p

Abstract Background Bovine tuberculosis (BTB) is caused by Mycobacterium bovis, which belongs to the Mycobacterium tuberculosis complex. Mycobacterium bovis have been described to be responsible of most cases of bovine tuberculosis. Although M. tuberculosis, M. africanum and non-complex mycobacteria were isolated from cattle. In Morocco, so far, no molecular studies were conducted to characterize the strains responsible of BTB. The present study aims to characterize M. bovis in Morocco. The present study was conducted in slaughterhouses in Rabat and El Jadida. Samples were collected from 327 slaughtered animals with visible lesions suggesting BTB. Results A total of 225 isolates yielded cultures, 95% (n = 215) of them were acid-fast (AF). Sixty eight per cent of the AF positive samples were confirmed as tuberculous mycobacteria (n = 147), 99% of these (n = 146) having RD9 and among the latter, 98% (n = 143) positive while 2% (n = 3) negative for RD4 A total of 134 samples were analyzed by spoligotyping of which 14 were in cluster and with 41 different spoligotypes, ten of them were new patterns (23%). The most prevalent spoligotypes were SB0121, SB0265, and SB0120, and were already identified in many other countries, such as Algeria, Spain, Tunisia, the United States and Argentina. Conclusion The shared borders between Algeria and Morocco, in addition to the previous importation of cattle from Europe and the US could explain the similarities found in M. bovis spoligotypes. On the other hand, the desert of Morocco could be considered as an efficient barrier preventing the introduction of BTB to Morocco from West Central and East Africa. Our findings suggest a low level endemic transmission of BTB similar to other African countries. However, more research is needed for further knowledge about the transmission patterns of BTB in Morocco.

Abstract Background Combination therapy is one of the most effective tools for limiting the emergence of drug resistance in pathogens. Despite the widespread adoption of combination therapy across diseases, drug resistance rates continue to rise, leading to failing treatment regimens. The mechanisms underlying treatment failure are well studied, but the processes governing successful combination therapy are poorly understood. We address this question by studying the population dynamics of Mycobacterium tuberculosis within tuberculosis patients undergoing treatment with different combinations of antibiotics. Results By combining very deep whole genome sequencing (~1000-fold genome-wide coverage) with sequential sputum sampling, we were able to detect transient genetic diversity driven by the apparently continuous turnover of minor alleles, which could serve as the source of drug-resistant bacteria. However, we report that treatment efficacy has a clear impact on the population dynamics: sufficient drug pressure bears a clear signature of purifying selection leading to apparent genetic stability. In contrast, M. tuberculosis populations subject to less drug pressure show markedly different dynamics, including cases of acquisition of additional drug resistance. Conclusions Our findings show that for a pathogen like M. tuberculosis, which is well adapted to the human host, purifying selection constrains the evolutionary trajectory to resistance in effectively treated individuals. Nonetheless, we also report a continuous turnover of minor variants, which could give rise to the emergence of drug resistance in cases of drug pressure weakening. Monitoring bacterial population dynamics could therefore provide an informative metric for assessing the efficacy of novel drug combinations.

The molecular clock and its phylogenetic applications to genomic data have changed how we study and understand one of the major human pathogens, Mycobacterium tuberculosis (MTB), the etiologic agent of tuberculosis. Genome sequences of MTB strains sampled at different times are increasingly used to infer when a particular outbreak begun, when a drug-resistant clone appeared and expanded, or when a strain was introduced into a specific region. Despite the growing importance of the molecular clock in tuberculosis research, there is a lack of consensus as to whether MTB displays a clocklike behavior and about its rate of evolution. Here we performed a systematic study of the molecular clock of MTB on a large genomic data set (6,285 strains), covering different epidemiological settings and most of the known global diversity. We found that sampling times below 15-20 years were often insufficient to calibrate the clock of MTB. For data sets where such calibration was possible, we obtained a clock rate between 1x10-8 and 5x10-7 nucleotide changes per-site-per-year (0.04-2.2 SNPs per-genome-per-year), with substantial differences between clades. These estimates were not strongly dependent on the time of the calibration points as they changed only marginally when we used epidemiological isolates (sampled in the last 40 years) or three ancient DNA samples (about 1,000 years old) to calibrate the tree. Additionally, the uncertainty and the discrepancies in the results of different methods were sometimes large, highlighting the importance of using different methods, and of considering carefully their assumptions and limitations.

[This corrects the article DOI: 10.1371/journal.ppat.1006939.].

Transition bias, an overabundance of transitions relative to transversions, has been widely reported among studies of the rates and spectra of spontaneous mutations. However, demonstrating the role of transition bias in adaptive evolution remains challenging. In particular, it is unclear whether such biases direct the evolution of bacterial pathogens adapting to treatment. We addressed this challenge by analyzing adaptive antibiotic-resistance mutations in the major human pathogen Mycobacterium tuberculosis (MTB). We found strong evidence for transition bias in two independently curated data sets comprising 152 and 208 antibiotic-resistance mutations. This was true at the level of mutational paths (distinct adaptive DNA sequence changes) and events (individual instances of the adaptive DNA sequence changes) and across different genes and gene promoters conferring resistance to a diversity of antibiotics. It was also true for mutations that do not code for amino acid changes (in gene promoters and the 16S ribosomal RNA gene rrs) and for mutations that are synonymous to each other and are therefore likely to have similar fitness effects, suggesting that transition bias can be caused by a bias in mutation supply. These results point to a central role for transition bias in determining which mutations drive adaptive antibiotic resistance evolution in a key pathogen.

The Lineage 2–Beijing (L2–Beijing) sub-lineage of Mycobacterium tuberculosis has received much attention due to its high virulence, fast disease progression, and association with antibiotic resistance. Despite several reports of the recent emergence of L2–Beijing in Africa, no study has investigated the evolutionary history of this sub-lineage on the continent. In this study, we used whole genome sequences of 781 L2 clinical strains from 14 geographical regions globally distributed to investigate the origins and onward spread of this lineage in Africa. Our results reveal multiple introductions of L2–Beijing into Africa linked to independent bacterial populations from East- and Southeast Asia. Bayesian analyses further indicate that these introductions occurred during the past 300 years, with most of these events pre-dating the antibiotic era. Hence, the success of L2–Beijing in Africa is most likely due to its hypervirulence and high transmissibility rather than drug resistance.

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.

BACKGROUND:Bovine tuberculosis (bTB) caused by Mycobacterium bovis is a re-emerging problem in both livestock and humans. The association of some M. bovis strains with hyper-virulence, MDR-TB and disseminated disease makes it imperative to understand the biology of the pathogen. METHODS:Mycobacterium bovis (15) among 1755 M. tuberculosis complex (MTBC) isolated between 2012 and 2014 were characterized and analyzed for associated patient demography and other risk factors. Five of the M. bovis isolates were whole-genome sequenced and comparatively analyzed against a global collection of published M. bovis genomes. RESULTS:Mycobacterium bovis was isolated from 3/560(0.5%) females and 12/1195(1.0%) males with pulmonary TB. The average age of M. bovis infected cases was 46.8 years (7-72years). TB patients from the Northern region of Ghana (1.9%;4/212) had a higher rate of infection with M. bovis (OR = 2.7,p = 0.0968) compared to those from the Greater Accra region (0.7%;11/1543). Among TB patients with available HIV status, the odds of isolating M. bovis from HIV patients (2/119) was 3.3 higher relative to non-HIV patients (4/774). Direct contact with livestock or their unpasteurized products was significantly associated with bTB (p

INTRODUCTION:Recurrent tuberculosis (TB) is one of the main challenges in TB control. Genotyping based on Mycobacterial Interspersed Repetitive Units-Variable Tandem Repeats (MIRU-VNTR) has been widely used to differentiate between relapse and reinfection, which are the two main causes of recurrent TB. There is a lack of data regarding the causes of TB recurrence in Georgia, and while differentiating between relapse and reinfection plays a key role in defining appropriate interventions, the required genotyping methodologies have not been implemented. The objective of this study was to implement MIRU-VNTR genotyping at the National Center for Tuberculosis and Lung Diseases (NCTBLD) and differentiate between relapse and reinfection in multidrug resistant (MDR-) TB patients from Tbilisi, Georgia. METHODS:Recurrent MDR tuberculosis cases from 2014-2016 diagnosed at NCTLD were included in the study when bacterial samples from both episodes were available. Genotyping based on the MIRU-VNTR 24 loci was implemented and used for differentiating between relapse and reinfection. Paired samples showing the same MIRU-VNTR pattern or one locus difference were classified as relapse, while two and more loci differences were treated as reinfection. Exact logistic regression was used to identify predictors of recurrence. RESULTS:Thirty two MDR-TB patients (64 samples) were included and MIRU-VNTR 24 typing was performed on the corresponding paired samples. Of the 32 patients, 25 (83.3%) were identified as relapse while 5 (16.7%) were due to re-infection. Patients with a history of incarceration were significantly associated with TB reinfection (p< 0.05). CONCLUSION:Recurrent TB in MDR patients in Georgia are mainly caused by relapse, raising concerns on the efficacy of the TB control program. An association between incarceration and reinfection likely reflects high levels of ongoing TB transmission in prisons, indicating the need for better TB infection control measures in these settings. Our results add to the rationale for implementing genotypic surveillance of TB more broadly to support TB control in Georgia.

BackgroundHuman tuberculosis (TB) is caused by seven phylogenetic lineages of the Mycobacterium tuberculosis complex (MTBC), Lineage 1-7. Recent advances in rapid genotyping of MTBC based on single nucleotide polymorphisms (SNP), allow for phylogenetically robust strain classification, paving the way for defining genotype-phenotype relationships in clinical settings. Such studies have revealed that, in addition to host and environmental factors, strain variation in the MTBC influences the outcome of TB infection and disease. In Tanzania, such molecular epidemiological studies of TB however are scarce in spite of a high TB burden.Methods and findingsHere we used SNP-typing to characterize a nationwide collection of 2,039 MTBC clinical isolates representative of 1.6% of all new and retreatment TB cases notified in Tanzania during 2012 and 2013. Four lineages, namely Lineage 1-4 were identified within the study population. The distribution and frequency of these lineages varied across regions but overall, Lineage 4 was the most frequent (n = 866, 42.5%), followed by Lineage 3 (n = 681, 33.4%) and 1 (n = 336, 16.5%), with Lineage 2 being the least frequent (n = 92, 4.5%). We found Lineage 2 to be independently associated with female sex (adjusted odds ratio [aOR] 2.14; 95% confidence interval [95% CI] 1.31 - 3.50, p = 0.002) and retreatment cases (aOR 1.67; 95% CI 0.95 - 2.84, p = 0. 065) in the study population. We found no associations between MTBC lineage and patient age or HIV status. Our sublineage typing based on spacer oligotyping on a subset of Lineage 1, 3 and 4 strains revealed the presence of mainly EAI, CAS and LAM families. Finally, we detected low levels of multidrug resistant isolates among a subset of 144 retreatment cases.ConclusionsThis study provides novel insights into the MTBC lineages and the possible influence of pathogen-related factors on the TB epidemic in Tanzania.

BackgroundDiabetes Mellitus (DM) is a known risk factor for tuberculosis (TB) but little is known on TB-Diabetes Mellitus (TBDM) co-morbidity in Sub-Saharan Africa.MethodsConsecutive TB cases registered at a tertiary facility in Ghana were recruited from September 2012 to April 2016 and screened for DM using random blood glucose and glycated hemoglobin (HbA1c) level. TB patients were tested for other clinical parameters including HIV co-infection and TB lesion location. Mycobacterial isolates obtained from collected sputum samples were characterized by standard methods. Associations between TBDM patients' epidemiological as well as microbiological variables were assessed.ResultsThe prevalence of DM at time of diagnosis among 2990 enrolled TB cases was 9.4% (282/2990). TBDM cases were significantly associated with weight loss, poor appetite, night sweat and fatigue (pConclusionOur findings confirms the importance of screening for diabetes during TB diagnosis and highlights the association between genetic diversity and diabetes. in Ghana.

BackgroundLimited data existed exclusively describing Mycobacterium tuberculosis lineage 3 (MTB-L3), sub-lineages, and clinical manifestations in Kampala, Uganda. This study sought to elucidate the circulating MTB-L3 sub-lineages and their corresponding clinical phenotypes.MethodA total of 141 M. tuberculosis isolates were identified as M. tuberculosis lineage 3 using Single nucleotide polymorphism (SNP) marker analysis method. To ascertain the sub-lineages/sub-strains within the M. tuberculosis lineage 3, the direct repeat (DR) loci for all the isolates was examined for sub-lineage specific signatures as described in the SITVIT2 database. The infecting sub-strains were matched with patients' clinical and demographic characteristics to identify any possible association.ResultThe data showed 3 sub-lineages circulating with CAS 1 Delhi accounting for 55% (77/141), followed by CAS 1-Kili 16% (22/141) and CAS 2/CAS 8% (12/141). Remaining isolates 21% (30/141) were unclassifiable. To explore whether the sub-lineages differ in their ability to cause increased severe disease, we used extent of lung involvement as a proxy for severe disease. Multivariable analysis showed no association between M. tuberculosis lineage 3 sub-lineages with severe disease. The risk factors associated with severe disease include having a positive smear (OR = 9.384; CI 95% = 2.603-33.835), HIV (OR = 0.316; CI 95% = 0.114-0.876), lymphadenitis (OR = 0. 171; CI 95% = 0.034-0.856) and a BCG scar (OR = 0.295; CI 95% = 0.102-0.854).ConclusionIn Kampala, Uganda, there are three sub-lineages of M. tuberculosis lineage 3 that cause disease of comparable severity with CAS-Dehli as the most prevalent. Having HIV, lymphadenitis, a BCG scar and a smear negative status is associated with reduced severe disease.

Bacteria belonging to the genus Mycobacterium are predominantly responsible for pulmonary diseases; most notably Mycobacterium tuberculosis causes granulomatous pulmonary infections. Here we describe a novel slow growing mycobacterial species isolated from respiratory samples from five patients, four with underlying pulmonary disease. The isolates were characterized by biochemical and molecular techniques, including whole genome sequencing. Biochemical characteristics generally match those of M. marinum and M. ulcerans; however, the most striking difference of the new species is its ability to grow at 37°C. The new species was found to grow in human macrophages, but not amoebae, suggesting a pathogenic rather than an environmental lifestyle. Phylogenetic analysis reveals a deep-rooting relationship to M. marinum and M. ulcerans. A complete genome sequence was obtained through combining short and long-read sequencing, providing a genome of 5.6 Mb. The genome appears to be highly intact, syntenic with that of M. marinum, with very few insertion sequences. A vast array of virulence factors includes 283 PE/PPE surface-associated proteins, making up 10% of the coding capacity, and 22 non-ribosomal peptide synthase clusters. A comparison of six clinical isolates from the five patients shows that they differ by up to two single nucleotide polymorphisms, suggesting a common source of infection. Our findings are in accordance with the recognition of a new taxonomic entity. We propose the name M. basiliense, as all isolates were found in patients from the Basel area of Switzerland.

Abstract Background Vitamin D level is inversely associated with tuberculosis (TB) and diabetes (DM). Vitamin D could be a mediator in the association between TB and DM. We examined the associations between vitamin D, TB and DM. Methods Consecutive adults with TB and sex- and age-matched volunteers were included in a case-control study in Dar es Salaam, Tanzania. Glycemia and total vitamin D (25(OH)D) were measured at enrolment and after TB treatment in cases. The association between low 25(OH)D (

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.

Bacterial factors favoring the unprecedented multidrug-resistant tuberculosis (MDR-TB) epidemic in the former Soviet Union remain unclear. We utilized whole genome sequencing and Bayesian statistics to analyze the evolutionary history, temporal emergence of resistance and transmission networks of MDR Mycobacterium tuberculosis complex isolates from Karakalpakstan, Uzbekistan (2001–2006). One clade (termed Central Asian outbreak, CAO) dating back to 1974 (95% HPD 1969–1982) subsequently acquired resistance mediating mutations to eight anti-TB drugs. Introduction of standardized WHO-endorsed directly observed treatment, short-course in Karakalpakstan in 1998 likely selected for CAO-strains, comprising 75% of sampled MDR-TB isolates in 2005/2006. CAO-isolates were also identified in a published cohort from Russia (2008–2010). Similarly, the presence of mutations supposed to compensate bacterial fitness deficits was associated with transmission success and higher drug resistance rates. The genetic make-up of these MDR-strains threatens the success of both empirical and standardized MDR-TB therapies, including the newly WHO-endorsed short MDR-TB regimen in Uzbekistan.

Once considered a phenotypically monomorphic bacterium, there is a growing body of work demonstrating heterogeneity among Mycobacterium tuberculosis (Mtb) strains in clinically relevant characteristics, including virulence and response to antibiotics. However, the genetic and molecular basis for most phenotypic differences among Mtb strains remains unknown. To investigate the basis of strain variation in Mtb, we performed genome-wide transposon mutagenesis coupled with next-generation sequencing (TnSeq) for a panel of Mtb clinical isolates and the reference strain H37Rv to compare genetic requirements for in vitro growth across these strains. We developed an analytic approach to identify quantitative differences in genetic requirements between these genetically diverse strains, which vary in genomic structure and gene content. Using this methodology, we found differences between strains in their requirements for genes involved in fundamental cellular processes, including redox homeostasis and central carbon metabolism. Among the genes with differential requirements were katG, which encodes the activator of the first-line antitubercular agent isoniazid, and glcB, which encodes malate synthase, the target of a novel small-molecule inhibitor. Differences among strains in their requirement for katG and glcB predicted differences in their response to these antimicrobial agents. Importantly, these strain-specific differences in antibiotic response could not be predicted by genetic variants identified through whole genome sequencing or by gene expression analysis. Our results provide novel insight into the basis of variation among Mtb strains and demonstrate that TnSeq is a scalable method to predict clinically important phenotypic differences among Mtb strains.

Tuberculosis (TB) induces a systemic inflammatory state affecting iron homeostasis. Patients with TB often have additional comorbidities such as anemia which can result in poorer treat outcomes. We studied the contribution of anemia and the role of the iron regulatory hormone hepcidin among TB patients and household contacts.We analyzed serum samples from 102 TB cases and 98 controls without TB, matched by age/sex, for hepcidin, iron, and inflammation parameters. Five controls developed TB within 12 months. We used linear regression to assess associations.Anemia of chronic disease (ACD) was more frequent among cases than controls (59.8% vs. 26.1%), but iron-deficiency anemia more frequent in controls (10% vs. 1%). The median hepcidin level was higher in cases than controls (63.7 vs. 14.2 ng/mL), but coinfections with HIV, helminths, and respiratory pathogens did not show cumulative effects. Hepcidin was associated with more severe TB symptom scoring (coefficient 0.8, 95% confidence interval [CI] 0.5-1.2) and higher mycobacterial load (0.7, 95% CI 0.4-1.0). Hepcidin was higher in TB cases and controls who developed TB compared to controls without TB (p

The already enormous burden caused by tuberculosis (TB) will be further aggravated by the association of this disease with modern epidemics, as human immunodeficiency virus and diabetes. Furthermore, the increasingly aging population and the wider use of suppressive immune therapies hold the potential to enhance the incidence of TB. New preventive and therapeutic strategies based on recent advances on our understanding of TB are thus needed. In particular, understanding the intricate network of events modulating inflammation in TB will help to build more effective vaccines and host-directed therapies to stop TB. This review integrates the impact of host, pathogen, and extrinsic factors on inflammation and the almost scientifically unexplored complexity emerging from the interactions between these three factors. We highlight the exciting data showing a contribution of this troika for the clinical outcome of TB and the need of incorporating it when developing novel strategies to rewire the immune response in TB.

Combining different molecular typing methods for Mycobacterium tuberculosis complex (MTBC) can be a powerful tool for molecular epidemiology-based investigation of TB. However, the current standard method that provides high discriminatory power for such a combination, mycobacterial interspersed repetitive units-variable numbers of tandem repeats typing (MIRU-VNTR), is laborious, time-consuming and often too costly for many resource-limited laboratories. We aimed to evaluate a reduced set of loci for MIRU-VNTR typing in combination with spoligotyping and SNP-typing for routine molecular epidemiology of TB.Spoligotyping and SNP-typing, in combination with the 15 loci MIRU-VNTR typing, were first used to type clinical MTBC isolates (n = 158) from Madagascar. A step by step reduction of MIRU-VNTR loci number was then performed according to the Hunter and Gaston Discriminatory Index (HGDI) and to the Principal component analysis (PCA) correlation with the spoligotype profiles to evaluate the discrimination power inside the generated spoligotype clusters. The 15 MIRU-VNTR was used as reference and SNP-typing was used to determine the main MTBC lineages.Of the 158 clinical isolates studied, the SNP-typing classified 23 into Lineage 1 (14.6%), 31 into Lineage 2 (19.6%), 23 into Lineage 3 (14.6%) and 81 into Lineage 4 strains (51.3%). 37 different spoligotypes profiles were obtained, 15 of which were unique and 20 in clusters. 15-loci MIRU-VNTR typing revealed 144 different genotypes: 132 isolates had a unique MIRU-VNTR profile and 27 isolates were grouped into 12 clusters. After a stepwise reduction of the MIRU-VNTR loci number within each main spoligotype families, three different sets composed of 5 customised MIRU-VNTR loci had a similar discrimination level to the reference 15 loci MIRU-VNTR in lineage 1, lineage 2 and lineage 3. For lineage 4, a set of 4 and 3 MIRU-VNTR loci were proposed to subtype the Harleem and LAM spoligotype families, respectively. For the T spoligotype family, a set of 5 MIRU-VNTR loci was proposed.According to the lineages and the spoligotype families, the number of MIRU-VNTR loci can be reduced to get an optimal classification of MTBC. These customized sets of MIRU-VNTR loci reduce workload and save resources while maintaining optimal discriminatory power.

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.

Tuberculosis (TB) is caused by gram-positive bacteria known as the Mycobacterium tuberculosis complex (MTBC). MTBC include several human-associated lineages and several variants adapted to domestic and, more rarely, wild animal species. We report an M. tuberculosis strain isolated from a wild chimpanzee in Côte d’Ivoire that was shown by comparative genomic and phylogenomic analyses to belong to a new lineage of MTBC, closer to the human-associated lineage 6 (also known as M. africanum West Africa 2) than to the other classical animal-associated MTBC strains. These results show that the general view of the genetic diversity of MTBC is limited and support the possibility that other MTBC variants exist, particularly in wild mammals in Africa. Exploring this diversity is crucial to the understanding of the biology and evolutionary history of this widespread infectious disease.

Molecular typing of 964 specimens from patients in Ethiopia with lymph node or pulmonary tuberculosis showed a similar distribution of Mycobacterium tuberculosis strains between the 2 disease manifestations and a minimal role for M. bovis. We report a novel phylogenetic lineage of M. tuberculosis strongly associated with the Horn of Africa.

BACKGROUND:There is a perception that genomic differences in the species/lineages of the nine species making the Mycobacterium tuberculosis complex (MTBC) may affect the efficacy of distinct control tools in certain geographical areas. We therefore analyzed the prevalence and spatial distribution of MTBC species and lineages among isolates from pulmonary TB cases over an 8-year period, 2007-2014. METHODOLOGY:Mycobacterial species isolated by culture from consecutively recruited pulmonary tuberculosis patients presenting at selected district/sub-district health facilities were confirmed as MTBC by IS6110 and rpoß PCR and further assigned lineages and sub lineages by spoligotyping and large sequence polymorphism PCR (RDs 4, 9, 12, 702, 711) assays. Patient characteristics, residency, and risks were obtained with a structured questionnaire. We used SaTScan and ArcMap analyses to identify significantly clustered MTBC lineages within selected districts and spatial display, respectively. RESULTS:Among 2,551 isolates, 2,019 (79.1%), 516 (20.2%) and 16 (0.6%) were identified as M. tuberculosis sensu stricto (MTBss), M. africanum (Maf), 15 M. bovis and 1 M. caprae, respectively. The proportions of MTBss and Maf were fairly constant within the study period. Maf spoligotypes were dominated by Spoligotype International Type (SIT) 331 (25.42%), SIT 326 (15.25%) and SIT 181 (14.12%). We found M. bovis to be significantly higher in Northern Ghana (1.9% of 212) than Southern Ghana (0.5% of 2339) (p = 0.020). Using the purely spatial and space-time analysis, seven significant MTBC lineage clusters (p< 0.05) were identified. Notable among the clusters were Ghana and Cameroon sub-lineages found to be associated with north and south, respectively. CONCLUSION:This study demonstrated that overall, 79.1% of TB in Ghana is caused by MTBss and 20% by M. africanum. Unlike some West African Countries, we did not observe a decline of Maf prevalence in Ghana.

Mycobacterium africanum is a member of the Mycobacterium tuberculosis complex (MTBC) and an important cause of human tuberculosis in West Africa that is rarely observed elsewhere. Here we genotyped 613 MTBC clinical isolates from Ghana, and searched for associations between the different phylogenetic lineages of MTBC and patient variables. We found that 17.1% (105/613) of the MTBC isolates belonged to M. africanum, with the remaining belonging to M. tuberculosis sensu stricto. No M. bovis was identified in this sample. M. africanum was significantly more common in tuberculosis patients belonging to the Ewe ethnic group (adjusted odds ratio: 3.02; 95% confidence interval: 1.67-5.47, p

Objective/background: Mycobacterium africanum (MAF) remains an important TB causing pathogen in West Africa; however, little is known about its population structure and actual diversity which may have implications for diagnostics and vaccines. We carried out comparative genomics analysis of candidate Mycobacterium tuberculosis (MTB) and MAF using whole genome sequencing. Methods: Clinical MTB complex strains (n = 187) comprising L4 (n = 22), L5 (n = 126), and L6 (n = 39) isolated over 8 years from Ghana were whole genome sequenced. The reads were mapped onto a reference genome for phylogenetic and functional genomics analysis. A maximum likelihood tree with 100 bootstraps was constructed from the single nucleotide polymorphisms (SNPs) found using RAxML and clustered with hierBAPS. A total of 147 (18 L4, 36 L6, and 93 L5) of the genomes were de novo assembled and annotated for comparative pangenome analysis using Roary. Results: The population structure of MAF revealed at least five clusters of L5 as compared to three for L6. We also identified a group of three multi-drug-resistants (MDRs) within a single cluster of L5 strains from Southern Ghana isolated in 2013. Among the global collection of MTB complex, there were four Ghana-specific L5 clusters of which one (L5.1.1) had traits of clonal expansion. From the 5947pan genes extracted from the collection, 3215 (54.1%) were core to all the 147 genomes whereas 719 (12.1%) were found in single genomes. Most of the variable genes were PE-PGRS/PPE (1,281) duplicates of other genes (431). The genome degradation was more pronounced in Lineages 4 and 6 as compared to Lineage 5. We identified the absence of some unique genes among specific lineages and/or clades with possible clinical implications. For example, mpt64 and mlaD encoding respectively an immunogenic protein and a mammalian cell entry protein were missing from all L6 genomes. In addition, all L5 strains had an amino acid substitution I43N within the mpt64 gene. Analysis of SNPs within some genes encoding proteins for substrate metabolism, ion transport and secretory systems showed higher proportion of SNPs among L6 compared to L5 and L4. We also identified a number of lineage/sublineage specific SNPs and indels that may be utilized in rapid PCR based genotyping of MTB complex. Conclusion: This work emphasizes on the possibility that the mpt64-based rapid diagnostic kit would not be effective in MAF endemic settings. More mutations in ESAT-6 secretory system of MAF compared to MTB sensu stricto can affect efficacy of ESAT-6-based vaccines in the future.

ABSTRACT The Mycobacterium tuberculosis genome includes the large family of pe_pgrs genes, whose functions are unknown. Because of precedents in other pathogens in which gene families showing high sequence variation are involved in antigenic variation, a similar role has been proposed for the pe_pgrs genes. However, the impact of immune selection on pe_pgrs genes has not been examined. Here, we sequenced 27 pe_pgrs genes in 94 clinical strains from five phylogenetic lineages of the M. tuberculosis complex (MTBC). We found that pe_pgrs genes were overall more diverse than the remainder of the MTBC genome, but individual members of the family varied widely in their nucleotide diversity and insertion/deletion (indel) content: some were more, and others were much less, diverse than the genome average. Individual pe_pgrs genes also differed in the ratio of nonsynonymous to synonymous mutations, suggesting that different selection pressures act on individual pe_pgrs genes. Using bioinformatic methods, we tested whether sequence diversity in pe_pgrs genes might be selected by human T cell recognition, the major mechanism of adaptive immunity to MTBC. We found that the large majority of predicted human T cell epitopes were confined to the conserved PE domain and experimentally confirmed the antigenicity of this domain in tuberculosis patients. In contrast, despite being genetically diverse, the PGRS domains harbored few predicted T cell epitopes. These results indicate that human T cell recognition is not a significant force driving sequence diversity in pe_pgrs genes, which is consistent with the previously reported conservation of human T cell epitopes in the MTBC. IMPORTANCE Recognition of Mycobacterium tuberculosis antigens by T lymphocytes is known to be important for immune protection against tuberculosis, but it is unclear whether human T cell recognition drives antigenic variation in M. tuberculosis. We previously discovered that the known human T cell epitopes in the M. tuberculosis complex are highly conserved, but we hypothesized that undiscovered epitopes with naturally occurring sequence variants might exist. To test this hypothesis, we examined the pe_pgrs genes, a large family of genes that has been proposed to function in immune evasion by M. tuberculosis. We found that the pe_pgrs genes exhibit considerable sequence variation, but the regions containing T cell epitopes and the regions of variation are distinct. These findings confirm that the majority of human T cell epitopes of M. tuberculosis are highly conserved and indicate that selection forces other than T cell recognition drive sequence variation in the pe_pgrs genes.

BackgroundDifferent combinations of variable number of tandem repeat (VNTR) loci have been proposed for genotyping Mycobacterium tuberculosis complex (MTBC). Existing VNTR schemes show different discriminatory capacity among the six human MTBC lineages. Here, we evaluated the discriminatory power of a "customized MIRU12" loci format proposed previously by Comas et al. based on the standard 24 loci defined by Supply et al. for VNTR-typing of MTBC in Ghana.MethodOne hundred and fifty-eight MTBC isolates classified into Lineage 4 and Lineage 5 were used to compare a customized lineage-specific panel of 12 MIRU-VNTR loci ("customized MIRU-12") to the standard MIRU-15 genotyping scheme. The resolution power of each typing method was determined based on the Hunter-Gaston- Discriminatory Index (HGDI). A minimal set of customized MIRU-VNTR loci for typing Lineages 4 (Euro-American) and 5 (M. africanum West African 1) strains from Ghana was defined based on the cumulative HGDI.Results and conclusionAmong the 106 Lineage 4 strains, the customized MIRU-12 identified a total of 104 distinct genotypes consisting of 2 clusters of 2 isolates each (clustering rate 1.8%), and 102 unique strains while standard MIRU-15 yielded a total of 105 different genotypes, including 1 cluster of 2 isolates (clustering rate: 0.9%) and 104 singletons. Among, 52 Lineage 5 isolates, customized MIRU-12 genotyping defined 51 patterns with 1 cluster of 2 isolates (clustering rate: 0.9%) and 50 unique strains whereas MIRU-15 classified all 52 strains as unique. Cumulative HGDI values for customized MIRU-12 for Lineages 4 and 5 were 0.98 respectively whilst that of standard MIRU-15 was 0.99. A union of loci from the customised MIRU-12 and standard MIRU-15 revealed a set of customized eight highly discriminatory loci: 4052, 2163B, 40, 4165, 2165, 10,16 and 26 with a cumulative HGDI of 0.99 for genotyping Lineage 4 and 5 strains from Ghana.

The phylogeographic population structure of Mycobacterium tuberculosis suggests local adaptation to sympatric human populations. We hypothesized that HIV infection, which induces immunodeficiency, will alter the sympatric relationship between M. tuberculosis and its human host. To test this hypothesis, we performed a nine-year nation-wide molecular-epidemiological study of HIV-infected and HIV-negative patients with tuberculosis (TB) between 2000 and 2008 in Switzerland. We analyzed 518 TB patients of whom 112 (21.6%) were HIV-infected and 233 (45.0%) were born in Europe. We found that among European-born TB patients, recent transmission was more likely to occur in sympatric compared to allopatric host-pathogen combinations (adjusted odds ratio [OR] 7.5, 95% confidence interval [95% CI] 1.21-infinity, p = 0.03). HIV infection was significantly associated with TB caused by an allopatric (as opposed to sympatric) M. tuberculosis lineage (OR 7.0, 95% CI 2.5-19.1, p

Current tools available to study the molecular epidemiology of tuberculosis do not provide information about the directionality and sequence of transmission for tuberculosis cases occurring over a short period of time, such as during an outbreak. Recently, whole genome sequencing has been used to study molecular epidemiology of Mycobacterium tuberculosis over short time periods.To describe the microevolution of M. tuberculosis during an outbreak caused by one drug-susceptible strain. METHOD AND MEASUREMENTS: We included 9 patients with tuberculosis diagnosed during a period of 22 months, from a population-based study of the molecular epidemiology in San Francisco. Whole genome sequencing was performed using Illumina's sequencing by synthesis technology. A custom program written in Python was used to determine single nucleotide polymorphisms which were confirmed by PCR product Sanger sequencing.We obtained an average of 95.7% (94.1-96.9%) coverage for each isolate and an average fold read depth of 73 (1 to 250). We found 7 single nucleotide polymorphisms among the 9 isolates. The single nucleotide polymorphisms data confirmed all except one known epidemiological link. The outbreak strain resulted in 5 bacterial variants originating from the index case A1 with 0-2 mutations per transmission event that resulted in a secondary case.Whole genome sequencing analysis from a recent outbreak of tuberculosis enabled us to identify microevolutionary events observable during transmission, to determine 0-2 single nucleotide polymorphisms per transmission event that resulted in a secondary case, and to identify new epidemiologic links in the chain of transmission.

South Africa shows one of the highest global burdens of multidrug-resistant (MDR) and extensively drug-resistant (XDR) tuberculosis (TB). Since 2002, MDR-TB in South Africa has been treated by a standardized combination therapy, which until 2010 included ofloxacin, kanamycin, ethionamide, ethambutol and pyrazinamide. Since 2010, ethambutol has been replaced by cycloserine or terizidone. The effect of standardized treatment on the acquisition of XDR-TB is not currently known.We genetically characterized a random sample of 4,667 patient isolates of drug-sensitive, MDR and XDR-TB cases collected from three South African provinces, namely, the Western Cape, Eastern Cape and KwaZulu-Natal. Drug resistance patterns of a subset of isolates were analyzed for the presence of commonly observed resistance mutations.Our analyses revealed a strong association between distinct strain genotypes and the emergence of XDR-TB in three neighbouring provinces of South Africa. Strains predominant in XDR-TB increased in proportion by more than 20-fold from drug-sensitive to XDR-TB and accounted for up to 95% of the XDR-TB cases. A high degree of clustering for drug resistance mutation patterns was detected. For example, the largest cluster of XDR-TB associated strains in the Eastern Cape, affecting more than 40% of all MDR patients in this province, harboured identical mutations concurrently conferring resistance to isoniazid, rifampicin, pyrazinamide, ethambutol, streptomycin, ethionamide, kanamycin, amikacin and capreomycin.XDR-TB associated genotypes in South Africa probably were programmatically selected as a result of the standard treatment regimen being ineffective in preventing their transmission. Our findings call for an immediate adaptation of standard treatment regimens for M/XDR-TB in South Africa.

There is increasing evidence that strain variation in Mycobacterium tuberculosis complex (MTBC) might influence the outcome of tuberculosis infection and disease. To assess genotype-phenotype associations, phylogenetically robust molecular markers and appropriate genotyping tools are required. Most current genotyping methods for MTBC are based on mobile or repetitive DNA elements. Because these elements are prone to convergent evolution, the corresponding genotyping techniques are suboptimal for phylogenetic studies and strain classification. By contrast, single nucleotide polymorphisms (SNP) are ideal markers for classifying MTBC into phylogenetic lineages, as they exhibit very low degrees of homoplasy. In this study, we developed two complementary SNP-based genotyping methods to classify strains into the six main human-associated lineages of MTBC, the "Beijing" sublineage, and the clade comprising Mycobacterium bovis and Mycobacterium caprae. Phylogenetically informative SNPs were obtained from 22 MTBC whole-genome sequences. The first assay, referred to as MOL-PCR, is a ligation-dependent PCR with signal detection by fluorescent microspheres and a Luminex flow cytometer, which simultaneously interrogates eight SNPs. The second assay is based on six individual TaqMan real-time PCR assays for singleplex SNP-typing. We compared MOL-PCR and TaqMan results in two panels of clinical MTBC isolates. Both methods agreed fully when assigning 36 well-characterized strains into the main phylogenetic lineages. The sensitivity in allele-calling was 98.6% and 98.8% for MOL-PCR and TaqMan, respectively. Typing of an additional panel of 78 unknown clinical isolates revealed 99.2% and 100% sensitivity in allele-calling, respectively, and 100% agreement in lineage assignment between both methods. While MOL-PCR and TaqMan are both highly sensitive and specific, MOL-PCR is ideal for classification of isolates with no previous information, whereas TaqMan is faster for confirmation. Furthermore, both methods are rapid, flexible and comparably inexpensive.

Mycobacterium tuberculosis complex (MTBC) genomes contain 2 large gene families termed pe and ppe. The function of pe/ppe proteins remains enigmatic but studies suggest that they are secreted or cell surface associated and are involved in bacterial virulence. Previous studies have also shown that some pe/ppe genes are polymorphic, a finding that suggests involvement in antigenic variation. Using comparative sequence analysis of 18 publicly available MTBC whole genome sequences, we have performed alignments of 33 pe (excluding pe_pgrs) and 66 ppe genes in order to detect the frequency and nature of genetic variation. This work has been supplemented by whole gene sequencing of 14 pe/ppe (including 5 pe_pgrs) genes in a cohort of 40 diverse and well defined clinical isolates covering all the main lineages of the M. tuberculosis phylogenetic tree. We show that nsSNP's in pe (excluding pgrs) and ppe genes are 3.0 and 3.3 times higher than in non-pe/ppe genes respectively and that numerous other mutation types are also present at a high frequency. It has previously been shown that non-pe/ppe M. tuberculosis genes display a remarkably low level of purifying selection. Here, we also show that compared to these genes those of the pe/ppe families show a further reduction of selection pressure that suggests neutral evolution. This is inconsistent with the positive selection pressure of "classical" antigenic variation. Finally, by analyzing such a large number of genes we were able to detect large differences in mutation type and frequency between both individual genes and gene sub-families. The high variation rates and absence of selective constraints provides valuable insights into potential pe/ppe function. Since pe/ppe proteins are highly antigenic and have been studied as potential vaccine components these results should also prove informative for aspects of M. tuberculosis vaccine design.

Mycobacterium tuberculosis, the causative agent of most human tuberculosis, infects one third of the world's population and kills an estimated 1.7 million people a year. With the world-wide emergence of drug resistance, and the finding of more functional genetic diversity than previously expected, there is a renewed interest in understanding the forces driving genome evolution of this important pathogen. Genetic diversity in M. tuberculosis is dominated by single nucleotide polymorphisms and small scale gene deletion, with little or no evidence for large scale genome rearrangements seen in other bacteria. Recently, a single report described a large scale genome duplication that was suggested to be specific to the Beijing lineage. We report here multiple independent large-scale duplications of the same genomic region of M. tuberculosis detected through whole-genome sequencing. The duplications occur in strains belonging to both M. tuberculosis lineage 2 and 4, and are thus not limited to Beijing strains. The duplications occur in both drug-resistant and drug susceptible strains. The duplicated regions also have substantially different boundaries in different strains, indicating different originating duplication events. We further identify a smaller segmental duplication of a different genomic region of a lab strain of H37Rv. The presence of multiple independent duplications of the same genomic region suggests either instability in this region, a selective advantage conferred by the duplication, or both. The identified duplications suggest that large-scale gene duplication may be more common in M. tuberculosis than previously considered.

M. africanum West African 2 constitutes an ancient lineage of the M. tuberculosis complex that commonly causes human tuberculosis in West Africa and has an attenuated phenotype relative to M. tuberculosis.In search of candidate genes underlying these differences, the genome of M. africanum West African 2 was sequenced using classical capillary sequencing techniques. Our findings reveal a unique sequence, RD900, that was independently lost during the evolution of two important lineages within the complex: the "modern" M. tuberculosis group and the lineage leading to M. bovis. Closely related to M. bovis and other animal strains within the M. tuberculosis complex, M. africanum West African 2 shares an abundance of pseudogenes with M. bovis but also with M. africanum West African clade 1. Comparison with other strains of the M. tuberculosis complex revealed pseudogenes events in all the known lineages pointing toward ongoing genome erosion likely due to increased genetic drift and relaxed selection linked to serial transmission-bottlenecks and an intracellular lifestyle.The genomic differences identified between M. africanum West African 2 and the other strains of the Mycobacterium tuberculosis complex may explain its attenuated phenotype, and pave the way for targeted experiments to elucidate the phenotypic characteristic of M. africanum. Moreover, availability of the whole genome data allows for verification of conservation of targets used for the next generation of diagnostics and vaccines, in order to ensure similar efficacy in West Africa.

Tuberculosis (TB) is a major public health problem in Nepal. Strain variation in Mycobacterium tuberculosis may influence the outcome of TB infection and disease. To date, the phylogenetic diversity of M. tuberculosis in Nepal is unknown.We analyzed 261 M. tuberculosis isolates recovered from pulmonary TB patients recruited between August 2009 and August 2010 in Nepal. M. tuberculosis lineages were determined by single nucleotide polymorphisms (SNP) typing and spoligotyping. Drug resistance was determined by sequencing the hot spot regions of the relevant target genes. Overall, 164 (62.8%) TB patients were new, and 97 (37.2%) were previously treated. Any drug resistance was detected in 50 (19.2%) isolates, and 16 (6.1%) were multidrug-resistant. The most frequent M. tuberculosis lineage was Lineage 3 (CAS/Delhi) with 106 isolates (40.6%), followed by Lineage 2 (East-Asian lineage, includes Beijing genotype) with 84 isolates (32.2%), Lineage 4 (Euro-American lineage) with 41 (15.7%) isolates, and Lineage 1 (Indo-Oceanic lineage) with 30 isolates (11.5%). Based on spoligotyping, we found 45 different spoligotyping patterns that were previously described. The Beijing (83 isolates, 31.8%) and CAS spoligotype (52, 19.9%) were the dominant spoligotypes. A total of 36 (13.8%) isolates could not be assigned to any known spoligotyping pattern. Lineage 2 was associated with female sex (adjusted odds ratio [aOR] 2.58, 95% confidence interval [95% CI] 1.42-4.67, p = 0.002), and any drug resistance (aOR 2.79; 95% CI 1.43-5.45; p = 0.002). We found no evidence for an association of Lineage 2 with age or BCG vaccination status.We found a large genetic diversity of M. tuberculosis in Nepal with representation of all four major lineages. Lineages 3 and 2 were dominating. Lineage 2 was associated with clinical characteristics. This study fills an important gap on the map of the M. tuberculosis genetic diversity in the Asian region.

In Switzerland and other developed countries, the number of tuberculosis (TB) cases has been decreasing for decades, but HIV-infected patients and migrants remain risk groups. The aim of this study was to compare characteristics of TB in HIV-negative and HIV-infected patients diagnosed in Switzerland, and between coinfected patients enrolled and not enrolled in the national Swiss HIV Cohort Study (SHCS).All patients diagnosed with culture-confirmed TB in the SHCS and a random sample of culture-confirmed cases reported to the national TB registry 2000-2008 were included. Outcomes were assessed in HIV-infected patients and considered successful in case of cure or treatment completion. Ninety-three SHCS patients and 288 patients selected randomly from 4221 registered patients were analyzed. The registry sample included 10 (3.5%) coinfected patients not enrolled in the SHCS: the estimated number of HIV-infected patients not enrolled in the SHCS but reported to the registry 2000-2008 was 146 (95% CI 122-173). Coinfected patients were more likely to be from sub-Saharan Africa (51.5% versus 15.8%, P

BackgroundMycobacterium tuberculosis has a global population structure consisting of six main phylogenetic lineages associated with specific geographic regions and human populations. One particular M. tuberculosis genotype known as "Beijing" has repeatedly been associated with drug resistance and has been emerging in some parts of the world. "Beijing" strains are traditionally defined based on a characteristic spoligotyping pattern. We used three alternative genotyping techniques to revisit the phylogenetic classification of M. tuberculosis complex (MTBC) strains exhibiting the typical "Beijing" spoligotyping pattern.Methods and findingsMTBC strains were obtained from an ongoing molecular epidemiological study in Switzerland and Nepal. MTBC genotyping was performed based on SNPs, genomic deletions, and 24-loci MIRU-VNTR. We identified three MTBC strains from patients originating from Tibet, Portugal and Nepal which exhibited a spoligotyping patterns identical to the classical Beijing signature. However, based on three alternative molecular markers, these strains were assigned to Lineage 3 (also known as Delhi/CAS) rather than to Lineage 2 (also known as East-Asian lineage). Sequencing of the RD207 in one of these strains showed that the deletion responsible for this "Pseudo-Beijing" spoligotype was about 1,000 base pairs smaller than the usual deletion of RD207 in classical "Beijing" strains, which is consistent with an evolutionarily independent deletion event in the direct repeat (DR) region of MTBC.ConclusionsWe provide an example of convergent evolution in the DR locus of MTBC, and highlight the limitation of using spoligotypes for strain classification. Our results indicate that a proportion of "Beijing" strains may have been misclassified in the past. Markers that are more phylogenetically robust should be used when exploring strain-specific differences in experimental or clinical phenotypes.

ObjectiveThe aim of this study was to use spoligotyping and large sequence polymorphism (LSP) to study the population structure of M. tuberculosis complex (MTBC) isolates.MethodsMTBC isolates were identified using standard biochemical procedures, IS6110 PCR, and large sequence polymorphisms. Isolates were further typed using spoligotyping, and the phenotypic drug susceptibility patterns were determined by the proportion method.ResultOne hundred and sixty-two isolates were characterised by LSP typing. Of these, 130 (80.25%) were identified as Mycobacterium tuberculosis sensu stricto (MTBss), with the Cameroon sub-lineage being dominant (N = 59/130, 45.38%). Thirty-two (19.75%) isolates were classified as Mycobacterium africanum type 1, and of these 26 (81.25%) were identified as West-Africa I, and 6 (18.75%) as West-Africa II. Spoligotyping sub-lineages identified among the MTBss included Haarlem (N = 15, 11.53%), Ghana (N = 22, 16.92%), Beijing (4, 3.08%), EAI (4, 3.08%), Uganda I (4, 3.08%), LAM (2, 1.54%), X (N = 1, 0.77%) and S (2, 1.54%). Nine isolates had SIT numbers with no identified sub-lineages while 17 had no SIT numbers. MTBss isolates were more likely to be resistant to streptomycin (pConclusionThis study demonstrated that overall 36.4% of TB in South-Western Ghana is caused by the Cameroon sub-lineage of MTBC and 20% by M. africanum type 1, including both the West-Africa 1 and West-Africa 2 lineages. The diversity of MTBC in Ghana should be considered when evaluating new TB vaccines.

Mycobacterium africanum consists of two phylogenetically distinct lineages within the Mycobacterium tuberculosis complex, known as M. africanum West African 1 and M. africanum West African 2. These lineages are restricted to West Africa, where they cause up to half of human pulmonary tuberculosis. In this review we discuss the definition of M. africanum, describe the prevalence and restricted geographical distribution of M. africanum West African 1 and 2, review the occurrence of M. africanum in animals, and summarize the phenotypic differences described thus far between M. africanum and M. tuberculosis sensu stricto.

Recent studies suggest that M. tuberculosis lineage and host genetics interact to impact how active tuberculosis presents clinically. We determined the phylogenetic lineages of M. tuberculosis isolates from participants enrolled in the Tuberculosis Trials Consortium Study 28, conducted in Brazil, Canada, South Africa, Spain, Uganda and the United States, and secondarily explored the relationship between lineage, clinical presentation and response to treatment. Large sequence polymorphisms and single nucleotide polymorphisms were analyzed to determine lineage and sublineage of isolates. Of 306 isolates genotyped, 246 (80.4%) belonged to the Euro-American lineage, with sublineage 724 predominating at African sites (99/192, 51.5%), and the Euro-American strains other than 724 predominating at non-African sites (89/114, 78.1%). Uneven distribution of lineages across regions limited our ability to discern significant associations, nonetheless, in univariate analyses, Euro-American sublineage 724 was associated with more severe disease at baseline, and along with the East Asian lineage was associated with lower bacteriologic conversion after 8 weeks of treatment. Disease presentation and response to drug treatment varied by lineage, but these associations were no longer statistically significant after adjustment for other variables associated with week-8 culture status.

Because genetically monomorphic bacterial pathogens harbour little DNA sequence diversity, most current genotyping techniques used to study the epidemiology of these organisms are based on mobile or repetitive genetic elements. Molecular markers commonly used in these bacteria include Clustered Regulatory Short Palindromic Repeats (CRISPR) and Variable Number Tandem Repeats (VNTR). These methods are also increasingly being applied to phylogenetic and population genetic studies. Using the Mycobacterium tuberculosis complex (MTBC) as a model, we evaluated the phylogenetic accuracy of CRISPR- and VNTR-based genotyping, which in MTBC are known as spoligotyping and Mycobacterial Interspersed Repetitive Units (MIRU)-VNTR-typing, respectively. We used as a gold standard the complete DNA sequences of 89 coding genes from a global strain collection. Our results showed that phylogenetic trees derived from these multilocus sequence data were highly congruent and statistically robust, irrespective of the phylogenetic methods used. By contrast, corresponding phylogenies inferred from spoligotyping or 15-loci-MIRU-VNTR were incongruent with respect to the sequence-based trees. Although 24-loci-MIRU-VNTR performed better, it was still unable to detect all strain lineages. The DNA sequence data showed virtually no homoplasy, but the opposite was true for spoligotyping and MIRU-VNTR, which was consistent with high rates of convergent evolution and the low statistical support obtained for phylogenetic groupings defined by these markers. Our results also revealed that the discriminatory power of the standard 24 MIRU-VNTR loci varied by strain lineage. Taken together, our findings suggest strain lineages in MTBC should be defined based on phylogenetically robust markers such as single nucleotide polymorphisms or large sequence polymorphisms, and that for epidemiological purposes, MIRU-VNTR loci should be used in a lineage-dependent manner. Our findings have implications for strain typing in other genetically monomorphic bacteria.

Renewed efforts in tuberculosis (TB) research have led to important new insights into the biology and epidemiology of this devastating disease. Yet, in the face of the modern epidemics of HIV/AIDS, diabetes, and multidrug resistance--all of which contribute to susceptibility to TB--global control of the disease will remain a formidable challenge for years to come. New high-throughput genomics technologies are already contributing to studies of TB's epidemiology, comparative genomics, evolution, and host-pathogen interaction. We argue here, however, that new multidisciplinary approaches--especially the integration of epidemiology with systems biology in what we call "systems epidemiology"--will be required to eliminate TB.

BACKGROUND:Mycobacterium tuberculosis complex (MTBC), the causative agent of tuberculosis (TB), is characterized by low sequence diversity making this bacterium one of the classical examples of a genetically monomorphic pathogen. Because of this limited DNA sequence variation, routine genotyping of clinical MTBC isolates for epidemiological purposes relies on highly discriminatory DNA fingerprinting methods based on mobile and repetitive genetic elements. According to the standard view, isolates exhibiting the same fingerprinting pattern are considered direct progeny of the same bacterial clone, and most likely reflect ongoing transmission or disease relapse within individual patients. METHODOLOGY/PRINCIPAL FINDINGS:Here we further investigated this assumption and used massively parallel whole-genome sequencing to compare one drug-susceptible (K-1) and one multidrug resistant (MDR) isolate (K-2) of a rapidly spreading M. tuberculosis Beijing genotype clone from a high incidence region (Karakalpakstan, Uzbekistan). Both isolates shared the same IS6110 RFLP pattern and the same allele at 23 out of 24 MIRU-VNTR loci. We generated 23.9 million (K-1) and 33.0 million (K-2) paired 50 bp purity filtered reads corresponding to a mean coverage of 483.5 fold and 656.1 fold respectively. Compared with the laboratory strain H37Rv both Beijing isolates shared 1,209 SNPs. The two Beijing isolates differed by 130 SNPs and one large deletion. The susceptible isolate had 55 specific SNPs, while the MDR variant had 75 specific SNPs, including the five known resistance-conferring mutations. CONCLUSIONS:Our results suggest that M. tuberculosis isolates exhibiting identical DNA fingerprinting patterns can harbour substantial genomic diversity. Because this heterogeneity is not captured by traditional genotyping of MTBC, some aspects of the transmission dynamics of tuberculosis could be missed or misinterpreted. Furthermore, a valid differentiation between disease relapse and exogenous reinfection might be impossible using standard genotyping tools if the overall diversity of circulating clones is limited. These findings have important implications for clinical trials of new anti-tuberculosis drugs.

Mycobacterium tuberculosis infects one third of the human world population and kills someone every 15 seconds. For more than a century, scientists and clinicians have been distinguishing between the human- and animal-adapted members of the M. tuberculosis complex (MTBC). However, all human-adapted strains of MTBC have traditionally been considered to be essentially identical. We surveyed sequence diversity within a global collection of strains belonging to MTBC using seven megabase pairs of DNA sequence data. We show that the members of MTBC affecting humans are more genetically diverse than generally assumed, and that this diversity can be linked to human demographic and migratory events. We further demonstrate that these organisms are under extremely reduced purifying selection and that, as a result of increased genetic drift, much of this genetic diversity is likely to have functional consequences. Our findings suggest that the current increases in human population, urbanization, and global travel, combined with the population genetic characteristics of M. tuberculosis described here, could contribute to the emergence and spread of drug-resistant tuberculosis.

The factors that govern the development of tuberculosis disease are incompletely understood. We hypothesized that some strains of Mycobacterium tuberculosis (M. tuberculosis) are more capable of causing disseminated disease than others and may be associated with polymorphisms in host genes responsible for the innate immune response to infection. We compared the host and bacterial genotype in 187 Vietnamese adults with tuberculous meningitis (TBM) and 237 Vietnamese adults with uncomplicated pulmonary tuberculosis. The host genotype of tuberculosis cases was also compared with the genotype of 392 cord blood controls from the same population. Isolates of M. tuberculosis were genotyped by large sequence polymorphisms. The hosts were defined by polymorphisms in genes encoding Toll-interleukin 1 receptor domain containing adaptor protein (TIRAP) and Toll-like receptor-2 (TLR-2). We found a significant protective association between the Euro-American lineage of M. tuberculosis and pulmonary rather than meningeal tuberculosis (Odds ratio (OR) for causing TBM 0.395, 95% confidence intervals (C.I.) 0.193-0.806, P = 0.009), suggesting these strains are less capable of extra-pulmonary dissemination than others in the study population. We also found that individuals with the C allele of TLR-2 T597C allele were more likely to have tuberculosis caused by the East-Asian/Beijing genotype (OR = 1.57 [95% C.I. 1.15-2.15]) than other individuals. The study provides evidence that M. tuberculosis genotype influences clinical disease phenotype and demonstrates, for the first time, a significant interaction between host and bacterial genotypes and the development of tuberculosis.

Aims and objectives: Resistant-, multidrug resistant- (MDR) and even extensively drug resistant- (XDR) Mycobacterium tuberculosis complex (MTBC) strains have emerged worldwide and represent a serious challenge for global tuberculosis (TB) control. Interestingly, high rates of MDR TB have been associated with particular phylogenetic lineages of the MTBC such as the Beijing lineage and the strong clonal expansion of particular MDR strains in Eastern Europe. However, a systematic investigation of the evolution and cross border transmission of particular highly transmissible MDR strains in Eastern Europe has not been performed so far. Methods: To address this question, an in-depth analysis was performed of the population structure (based on 24 loci MIRU-VNTR-typing and spoligotyping) of strains from different settings in Eastern Europe, such as Uzbekistan, Turkmenistan, and Kazakhstan. Results: The results showed that two major clones of the Beijing lineage are the drivers of the MDR epidemic in Eastern Europe. In addition, retrospective investigation of an MDR strain collection from Germany using MIRU-VNTR and IS6110 RFLP patterns revealed that both clones were circulating for more than 15 years. In-depth analysis of genome characteristics confirmed the longitudinal transmission of the two major MDR clones and allowed the identification of particular beneficial combinations of mutations potentially compensating for the fitness defect of the MDR/XDR strains. Moreover, these allele combinations coincide with an enhanced transmissibility compared with other circulating MDR/XDR strain types and promote the selection of extremely drug resistant but also highly virulent MTBC strain variants. Conclusion: Taken together, these data argue strongly for an expansion of particular MDR Beijing clones in Eastern Europe that started at least 15 years ago. Evolution and selection over the past years has resulted in “hyper transmissible” MDR-TB strains that have the potential to spread rapidly and, thus, further complicate TB control.