Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

An outbreak of pulmonary tuberculosis and a follow-up investigation of latent tuberculosis in a high school in an eastern city in China, 2016-2019

Background

In October 2016, a senior high school student was diagnosed with sputum-smear positive [SS(+)] pulmonary tuberculosis (TB). We conducted an investigation of an outbreak in the school, including among students and teachers diagnosed with latent TB, who we followed until July 2019.

Methods

We defined latent TB infection (LTBI) as a tuberculin skin test (TST) induration of 15mm or larger; probable TB as a chest radiograph indicative of TB plus productive cough/hemoptysis for at least 2 weeks, or TST induration of 15mm or larger; and confirmed TB as two or more positive sputum smears or one positive sputum smear plus a chest radiograph indicative of TB or culture positive with M. tuberculosis. We conducted mycobacterial interspersed repetitive unit–variable number tandem repeat (MIRU-VNTR) typing based on 24 loci in the isolates.

Results

Between October 2016 and July 2019, we identified 52 cases, including nine probable, six confirmed, and 37 LTBI cases. The index case-student had attended school continuously despite having TB symptoms for almost three months before being diagnosed with TB. We obtained three isolates from classmates of the index case in 2016; all had identical MIRU-VNTR alleles with the index case. The LTBI rate was lower among students (7.41%, 30/405) than among teachers (26.92%, 7/26) (rate ratio [RR] = 0.28, 95% confidential interval [CI]: 0.13–0.57). Among the 17 students who had latent TB and refused prophylaxis in October 2016, 23.53% (4/17) became probable/confirmed cases by July 2019. None of the six teachers who also refused prophylaxis became probable or confirmed cases. Of the 176 students who were TST(-) in October 2016, 1.70% (3/176) became probable/confirmed cases, and among the 20 teachers who were TST(-), 1 became a probable case.

Conclusions

Delayed diagnosis of TB in the index patient may have contributed to the start of this outbreak; lack of post-exposure chemoprophylaxis facilitated spread of the outbreak. Post-exposure prophylaxis is strongly recommended for all TST-positive students; TST-negative students exposed to an SS(+) case should be followed up regularly so that prophylaxis can be started if LTBI is detected.

Predictive modeling to study the treatment-shortening potential of novel tuberculosis drug regimens, towards bundling of preclinical data

BACKGROUND

Given the persistently high global burden of tuberculosis (TB), effective and shorter treatment options are needed. Here, we explore the relationship between relapse and treatment length as well as inter-regimen differences for two novel anti-TB drug regimens using a mouse model of TB infection and mathematical modeling.

METHODS

Mycobacterium tuberculosis-infected mice were treated for up to 13 weeks with bedaquiline and pretomanid combined with moxifloxacin and pyrazinamide (BPaMZ) or linezolid (BPaL). Cure rates were evaluated 12 weeks after treatment completion. The standard regimen of isoniazid, rifampicin, pyrazinamide, and ethambutol (HRZE) was evaluated as a comparator.

RESULTS

Six weeks of BPaMZ was sufficient to cure all mice. In contrast, 13 weeks of BPaL and 24 weeks of HRZE did not achieve 100% cure rates. Based on mathematical model predictions, 95% probability of cure was predicted for BPaMZ, BPaL and HRZE to occur at 1.6, 4.3, and 7.9 months, respectively.

CONCLUSION

This study provides additional evidence for the treatment-shortening capacity of BPaMZ over BPaL and HRZE. To optimally utilize preclinical data for predicting clinical outcomes, and to overcome the limitations that hamper such extrapolation, we advocate bundling of available published preclinical data into mathematical models.

Human tuberculosis and Mycobacterium tuberculosis complex: A review on genetic diversity, pathogenesis and omics approaches in host biomarkers discovery

Mycobacterium tuberculosis complex (MTBC) refers to a group of mycobacteria encompassing nine members of closely related species that causes tuberculosis in animals and humans. Among the nine members, Mycobacterium tuberculosis (M. tuberculosis) remains the main causative agent for human tuberculosis that results in high mortality and morbidity globally. In general, MTBC species are low in diversity but exhibit distinctive biological differences and phenotypes among different MTBC lineages. MTBC species are likely to have evolved from a common ancestor through insertions/deletions processes resulting in species speciation with different degrees of pathogenicity. The pathogenesis of human tuberculosis is complex and remains poorly understood. It involves multi-interactions or evolutionary co-options between host factors and bacterial determinants for survival of the MTBC. Granuloma formation as a protection or survival mechanism in hosts by MTBC remains controversial. Additionally, MTBC species are capable of modulating host immune response and have adopted several mechanisms to evade from host immune attack in order to survive in humans. On the other hand, current diagnostic tools for human tuberculosis are inadequate and have several shortcomings. Numerous studies have suggested the potential of host biomarkers in early diagnosis of tuberculosis, in disease differentiation and in treatment monitoring. "Multi-omics" approaches provide holistic views to dissect the association of MTBC species with humans and offer great advantages in host biomarkers discovery. Thus, in this review, we seek to understand how the genetic variations in MTBC lead to species speciation with different pathogenicity. Furthermore, we also discuss how the host and bacterial players contribute to the pathogenesis of human tuberculosis. Lastly, we provide an overview of the journey of "omics" approaches in host biomarkers discovery in human tuberculosis and provide some interesting insights on the challenges and directions of "omics" approaches in host biomarkers innovation and clinical implementation.

Globally diverse Mycobacterium tuberculosis resistance acquisition: a retrospective geographical and temporal analysis of whole genome sequences

Background

Mycobacterium tuberculosis whole genome sequencing (WGS) data can provide insights into temporal and geographical trends in resistance acquisition and inform public health interventions. We aimed to use a large clinical collection of M tuberculosis WGS and resistance phenotype data to study how, when, and where resistance was acquired on a global scale.

Methods

We did a retrospective analysis of WGS data. We curated a set of clinical M tuberculosis isolates with high-quality sequencing and culture-based drug susceptibility data (spanning four lineages and 52 countries in Africa, Asia, the Americas, and Europe) using public databases and literature curation. For inclusion, sequence quality criteria and country of origin data were required. We constructed geographical and lineage specific M tuberculosis phylogenies and used Bayesian molecular dating with BEAST, version 1.10.4, to infer the most recent common susceptible ancestor age for 4869 instances of resistance to ten drugs.

Findings

Between Jan 1, 1987, and Sept 12, 2014, of 10 299 M tuberculosis clinical isolates, 8550 were curated, of which 6099 (71%) from 15 countries met criteria for molecular dating. The number of independent resistance acquisition events was lower than the number of resistant isolates across all countries, suggesting ongoing transmission of drug resistance. Ancestral age distributions supported the presence of old resistance, 20 years or more before, in most countries. A consistent order of resistance acquisition was observed globally starting with resistance to isoniazid, but resistance ancestral age varied by country. We found a direct correlation between gross domestic product per capita and resistance age (r ²=0·47; p=0·014). Amplification of fluoroquinolone and second-line injectable resistance among multidrug-resistant isolates is estimated to have occurred very recently (median ancestral age 4·7 years [IQR 1·9-9·8] before sample collection). We found the sensitivity of commercial molecular diagnostics for second-line resistance to vary significantly by country (p<0·0003).

Interpretation

Our results highlight that both resistance transmission and amplification are contributing to disease burden globally but vary by country. The observation that wealthier nations are more likely to have old resistance (most recent common susceptible ancestor >20 years before isolation) suggests that programmatic improvements can reduce resistance amplification, but that fit resistant strains can circulate for decades subsequently implies the need for continued surveillance.

Use of a whole genome sequencing-based approach for Mycobacterium tuberculosis surveillance in Europe in 2017-2019: an ECDC pilot study

Whole genome sequencing (WGS) can be used for molecular typing and characterisation of Mycobacterium tuberculosis complex (MTBC) strains. We evaluated the systematic use of a WGS-based approach for MTBC surveillance involving all European Union/European Economic Area (EU/EEA) countries and highlight the challenges and lessons learnt to be considered for the future development of a WGS-based surveillance system.

WGS and epidemiological data of patients with rifampicin-resistant (RR) and multidrug-resistant (MDR) tuberculosis (TB) were collected from EU/EEA countries between January 2017 and December 2019. WGS-based genetic relatedness analysis was performed using a standardised approach including both core genome multilocus sequence typing (cgMLST) and single nucleotide polymorphism (SNP)-based calculation of distances on all WGS data that fulfilled minimum quality criteria to ensure data comparability.

A total of 2218 RR/MDR-MTBC isolates were collected from 25 countries. Among these, 56 cross-border clusters with increased likelihood of recent transmission (≤5 SNPs distance) comprising 316 RR/MDR-MTBC isolates were identified. The cross-border clusters included between two and 30 resistant isolates from two to six countries, demonstrating different RR/MDR-TB transmission patterns in Western and Eastern EU countries.

This pilot study shows that a WGS-based surveillance system is not only feasible but can efficiently elucidate the dynamics of in-country and cross-border RR/MDR-TB transmission across EU/EEA countries. Lessons learnt from this study highlight that the establishment of an EU/EEA centralised WGS-based surveillance system for TB will require strengthening of national integrated systems performing prospective WGS surveillance and the development of clear procedures to facilitate international collaboration for the investigation of cross-border clusters.

The implementation of a WGS-based surveillance system for monitoring the emergence of MDR-TB outbreaks in Europe is feasible and has the potential to provide supporting evidence to better elucidate cross-border transmission patternshttps://bit.ly/2ZTnPjk

Effect of mixed strain infections on clinical and epidemiological features of tuberculosis in Florida

Mixed infections with genetically distinct Mycobacterium tuberculosis (MTB) strains within a single host have been documented in different settings; however, this phenomenon is rarely considered in the management and care of new and relapse tuberculosis (T.B.) cases. This study aims to establish the epidemiological and clinical features of mixed infections among culture-confirmed T.B. patients enrolled in tuberculosis care at the Florida Department of Health (FDOH) and measure its association with T.B. mortality. We analyzed de-identified surveillance data of T.B. cases enrolled in T.B. care from April 2008 to January 2018. Mixed MTB infection was determined by the presence of more than one Copy Number Variant (CNV) in at least one locus, based on the genotype profile pattern of at least one isolate using 24-locus Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR). The prevalence of mixed MTB infections among the 4944 culture-confirmed TB cases included in this analysis was 2.6% (129). Increased odds of mixed infections were observed among middle-aged patients, 45-64 years (AOR = 2.38; 95% CI: 0.99, 5.69; p = 0.0513), older adults 65 years and above (AOR = 3.95; 95% CI: 1.63, 9.58; p = 0.0023) and patients with diabetes (OR = 1.77; 95% CI: 1.12, 2.80; p = 0.0150). There was no significant association between mixed infections and death. Our study provides insight into the epidemiological and clinical characteristics of patients with mixed MTB infections, which is essential in the management of T.B. patients.

Role of Immigration in Tuberculosis Transmission to Iran: A Systematic Review

Background:

Today, because of increasing immigration and the prevalence of drug-resistant tuberculosis in Iran, identifying intra-community cases is necessary in the country. It will be possible through the use of molecular epidemiologic methods. In this inquiry, in order to determine the role of immigrants in the transmission of specific strains to Iran, the studies have been examined which had been conducted based on molecular epidemiologic methods among Iranians and non-Iranians people.

Methods:

All studies from 1997 to the end of March 2017 were examined in three databases of PubMed, Scopus, and Google Scholar and finally, 16 studies were selected.

Results:

The common clustering rate between Iranians and non-Iranians was determined to be 19.8, and the intra-community recent transmission rate was from 0% to 49% with average of 18.1%. The rate of multidrug-resistant tuberculosis (MDR-TB) was 12.5%, which was higher among immigrants, especially Afghans, and a significant number of the strains were Beijing.

Conclusions:

The studies have shown that migrants, especially Afghans, are more effective in transmitting specific strains of tuberculosis to migratory areas. To control tuberculosis, it is necessary to register of immigrant's health information, while enter to the country, so that, by doing appropriate diagnostic tests, the curing the patients, the transmission of tuberculosis to the country would be prevented.

Robust detection of point mutations involved in multidrug-resistant Mycobacterium tuberculosis in the presence of co-occurrent resistance markers

Tuberculosis disease is a major global public health concern and the growing prevalence of drug-resistant Mycobacterium tuberculosis is making disease control more difficult. However, the increasing application of whole-genome sequencing as a diagnostic tool is leading to the profiling of drug resistance to inform clinical practice and treatment decision making. Computational approaches for identifying established and novel resistance-conferring mutations in genomic data include genome-wide association study (GWAS) methodologies, tests for convergent evolution and machine learning techniques. These methods may be confounded by extensive co-occurrent resistance, where statistical models for a drug include unrelated mutations known to be causing resistance to other drugs. Here, we introduce a novel ‘cannibalistic’ elimination algorithm (“Hungry, Hungry SNPos”) that attempts to remove these co-occurrent resistant variants. Using an M. tuberculosis genomic dataset for the virulent Beijing strain-type (n = 3,574) with phenotypic resistance data across five drugs (isoniazid, rifampicin, ethambutol, pyrazinamide, and streptomycin), we demonstrate that this new approach is considerably more robust than traditional methods and detects resistance-associated variants too rare to be likely picked up by correlation-based techniques like GWAS.

Tuberculosis is one of the deadliest infectious diseases, being responsible for more than one million deaths per year. The causing bacteria are becoming increasingly drug-resistant, which is hampering disease control. At the same time, an unprecedented amount of bacterial whole-genome sequencing is increasingly informing clinical practice. In order to detect the genetic alterations responsible for developing drug resistance and predict resistance status from genomic data, bio-statistical methods and machine learning models have been employed. However, due to strongly overlapping drug resistance phenotypes and genotypes in multidrug-resistant datasets, the results of these correlation-based approaches frequently also contain mutations related to resistance against other drugs. In the past, this issue has often been ignored or partially resolved by either restricting the input data or in post-analysis screening—with both strategies relying on prior information. Here we present a heuristic algorithm for finding resistance-associated variants and demonstrate that it is considerably more robust towards co-occurrent resistance compared to traditional techniques. The software is available at https://github.com/julibeg/HHS.

Successful expansion of Mycobacterium tuberculosis Latin American and Mediterranean sublineage (L4.3/LAM) in Tunisia mainly driven by a single, long-established clonal complex

OBJECTIVES

To explore the evolutionary history of Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in Tunisia, where it predominates.

METHODS

High-resolution genotyping of 252 L4.3/LAM clinical strains was undertaken, and whole-genome sequencing was performed on 31 representative isolates.

RESULTS

Genotyping data coupled with Bayesian analyses split the Tunisian L4.3/LAM strain collection into two divergent entities (65.07% vs 34.92%): a major subpopulation, dominated by a single clonal complex (CC), TUN4.3_CC1 (94.51%); and a minor subpopulation, dominated by TUN4.3_CC2 (42.04%). TUN4.3_CC1 is clearly thriving in Tunisia, accounting for 61.5% of the L4.3/LAM sublineage. TUN4.3_CC1 displayed higher mean allelic richness compared with TUN4.3_CC2 and predominated throughout the entire region, indicating a long-established history. The very low proportion of drug resistance among TUN4.3_CC1 isolates is indicative of their intrinsic ability to spread successfully in the host population. Genomic analyses further confirmed the clear genetic separation between the two main CCs (pairwise fixation index 0.56), and suggested the relatively ancient origin of TUN4.3_CC1. Consistent with its successful expansion, TUN4.3_CC1 showed reduced mean pairwise genetic distance between genomes.

CONCLUSIONS

These findings link the successful expansion of L4.3/LAM in Tunisia to a single long-established clone.

Emergence of multidrug-resistant Mycobacterium tuberculosis of the Beijing lineage in Portugal and Guinea-Bissau: a snapshot of moving clones by whole-genome sequencing

The Beijing genotype comprises a highly disseminated strain type that is frequently associated with multidrug resistant (MDR) tuberculosis (TB) and increased transmissibility but, countries such as Portugal and Guinea-Bissau fall outside the regions phylogeographically associated with this specific genotype. Nevertheless, recent data shows that this genotype might be gradually emerging in these two countries as an underlying cause of primary MDR-TB. Here, we describe the emergence of Mycobacterium tuberculosis Beijing strains associated with MDR-TB in Portugal and Guinea-Bissau demonstrating the presence of the well described superclusters 100-32 and 94-32 in Portugal and Guinea-Bissau, respectively. Genome-wide analysis and comparison with a global genomic dataset of M. tuberculosis Beijing strains, revealed the presence of two genomic clusters encompassing isolates from Portugal and Guinea-Bissau, GC1 (n = 121) and GC2 (n = 39), both of which bore SNP signatures compatible with the 100-32/B0/W148 and 94-32/Central Asia Outbreak clades, respectively. Moreover, GC2 encompasses a cross-border cluster between Portugal, Guinea-Bissau and Brazil thus supporting migration-associated introduction of MDR-TB and subsequent clonal expansion at the community-level. The comparison with global Beijing datasets demonstrates the global reach of the disease and its complex dissemination across multiple countries while in parallel there are clear microevolutionary trajectories towards extensively drug resistant TB.

Long-term molecular surveillance provides clues on a cattle origin for Mycobacterium bovis in Portugal

Animal tuberculosis (TB), caused by Mycobacterium bovis, is maintained in Portugal in a multi-host system, with cattle, red deer and wild boar, playing a central role. However, the ecological processes driving transmission are not understood. The main aim of this study was thus to contribute to the reconstruction of the spatiotemporal history of animal TB and to refine knowledge on M. bovis population structure in order to inform novel intervention strategies. A collection of 948 M. bovis isolates obtained during long-term surveillance (2002–2016, 15 years) of cattle (n = 384), red deer (n = 303) and wild boar (n = 261), from the main TB hotspot areas, was characterized by spoligotyping and 8 to 12-loci MIRU-VNTR. Spoligotyping identified 64 profiles and MIRU-VNTR distinguished 2 to 36 subtypes within each spoligotype, enabling differentiation of mixed or clonal populations. Common genotypic profiles within and among livestock and wildlife in the same spatiotemporal context highlighted epidemiological links across hosts and regions, as for example the SB0119-M205 genotype shared by cattle in Beja district or SB0121-M34 shared by the three hosts in Castelo Branco and Beja districts. These genomic data, together with metadata, were integrated in a Bayesian inference framework, identifying five ancestral M. bovis populations. The phylogeographic segregation of M. bovis in specific areas of Portugal where the disease persists locally is postulated. Concurrently, robust statistics indicates an association of the most probable ancient population with cattle and Beja, providing a clue on the origin of animal TB epidemics. This relationship was further confirmed through a multinomial probability model that assessed the influence of host species on spatiotemporal clustering. Two significant clusters were identified, one that persisted between 2004 and 2010, in Beja district, with Barrancos county at the centre, overlapping the central TB core area of the Iberian Peninsula, and highlighting a significant higher risk associated to cattle. The second cluster was predominant in the 2012–2016 period, holding the county Rosmaninhal at the centre, in Castelo Branco district, for which wild boar contributed the most in relative risk. These results provide novel quantitative insights beyond empirical perceptions, that may inform adaptive TB control choices in different regions.

The value of the continuous genotyping of multi-drug resistant tuberculosis over 20 years in Spain

Molecular epidemiology of circulating clinical isolates is crucial to improve prevention strategies. The Spanish Working Group on multidrug resistant tuberculosis (MDR-TB) is a network that monitors the MDR-TB isolates in Spain since 1998. The aim of this study was to present the study of the MDR-TB and extensively drug-resistant tuberculosis (XDR-TB) patterns in Spain using the different recommended genotyping methods over time by a national coordinated system. Based on the proposed genotyping methods in the European Union until 2018, the preservation of one method, MIRU-VNTR, applied to selected clustered strains permitted to maintain our study open for 20 years. The distribution of demographic, clinical and epidemiological characteristics of clustered and non-clustered cases of MDR/XDR tuberculosis with proportion differences as assessed by Pearson’s chi-squared or Fisher’s exact test was compared. The differences in the quantitative variables using the Student's-t test and the Mann–Whitney U test were evaluated. The results obtained showed a total of 48.4% of the cases grouped in 77 clusters. Younger age groups, having a known TB case contact (10.2% vs 4.7%) and XDR-TB (16.5% vs 1.8%) were significantly associated with clustering. The largest cluster corresponded to a Mycobacterium bovis strain mainly spread during the nineties. A total of 68.4% of the clusters detected were distributed among the different Spanish regions and six clusters involving 104 cases were grouped in 17 and 18 years. Comparison of the genotypes obtained with those European genotypes included in The European Surveillance System (TESSy) showed that 87 cases had become part of 20 European clusters. The continuity of MDR strain genotyping in time has offered a widespread picture of the situation that allows better management of this public health problem. It also shows the advantage of maintaining one genotyping method over time, which allowed the comparison between ancient, present and future samples.

Circulating strains of Mycobacterium tuberculosis: 24 loci MIRU-VNTR analysis in Bangladesh

Bangladesh is among the high burden countries for tuberculosis (TB) and multidrug resistant TB (MDR-TB). As the genetic diversity and distinct phylogeographic distribution of Mycobacterium tuberculosis are responsible for regional differences in drug resistance, this cross sectional study was conducted to identify the circulating M. tuberculosis strains belonging to different lineages among pulmonary tuberculosis and, to investigate the contribution of distinct M. tuberculosis lineages to rifampicin resistant (RR) and rifampicin sensitive (RS) TB. A total of 40 RR and 20 RS isolates were enrolled in this study, all of which confirmed as M. tuberculosis by MPT 64 antigen detection. Furthermore, all isolates were genotyped by 24 loci Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR), thus comprising the first study to employ this approach in Bangladesh. Beijing was the predominant lineage (26.8%) followed by EAI (23.2%), Delhi/CAS (16.1%), H37Rv (8.9%), Haarlem (7.1%), LAM (5.4%), Cameroon (3.6%) and a NEW-1 (1.8%). Four (7.1%) isolates remained as unidentified. Beijing strains were the significantly predominant (36.8%; p = 0.0135) among the RR isolates in comparison with other strains whereas EAI was the predominant (38.8%) lineage among RS isolates. Also, approximately 13% RR isolates showed genotypic resistance against fluoroquinolones by LPA and, hence, classed as pre-XDR TB albeit no specific lineage was found associated with these latter strains. A low transmission rate (10.5%) and high genetic diversity was detected in this setting with all the clustered strains herein identified belonging to the Beijing lineage. This study highlights 24 loci MIRU-VNTR analysis as a powerful tool for genotyping of Mycobacterium tuberculosis in this setting as it shows a high discriminatory index (0.81).

Probable long-term prevalence for a predominant Mycobacterium tuberculosis clone of a Beijing genotype in Colon, Panama

Beijing genotype Mycobacterium tuberculosis strains associate with increased virulence, resistance and/or higher transmission rates. This study describes a specific Beijing strain predominantly identified in the Panamanian province of Colon with one of the highest incidences of tuberculosis in the country. Retrospective mycobacterial interspersed repetitive unit/variable number of tandem repeats analysis of 42 isolates collected between January and August 2018 allowed to identify a cluster (Beijing A) with 17 (40.5%) Beijing isolates. Subsequent prospective strain-specific PCR-based surveillance from September 2019 to March 2020 confirmed the predominance of the Beijing A strain (44.1%) in this province. Whole-genome sequencing revealed higher-than-expected diversity within the cluster, suggesting long-term prevalence of this strain and low number of cases caused by recent transmission. The Beijing A strain belongs to the Asian African 3 (Bmyc13, L2.2.5) branch of the modern Beijing sublineage, with their closest isolates corresponding to cases from Vietnam, probably introduced in Panama between 2000 and 2012.

Higher genome mutation rates of Beijing lineage of Mycobacterium tuberculosis during human infection

Mycobacterium tuberculosis (Mtb) strains of Beijing lineage have caused great concern because of their rapid emergence of drug resistance and worldwide spread. DNA mutation rates that reflect evolutional adaptation to host responses and the appearance of drug resistance have not been elucidated in human-infected Beijing strains. We tracked and obtained an original Mtb isolate of Beijing lineage from the 1999 tuberculosis outbreak in Japan, as well as five other isolates that spread in humans, and two isolates from the patient caused recurrence. Three isolates were from patients who developed TB within one year after infection (rapid-progressor, RP), and the other three isolates were from those who developed TB more than one year after infection (slow-progressor, SP). We sequenced genomes of these isolates and analyzed the propensity and rate of genomic mutations. Generation time versus mutation rate curves were significantly higher for RP. The ratio of oxidative versus non-oxidation damages induced mutations was higher in SP than RP, suggesting that persistent Mtb are exposed to oxidative stress in the latent state. Our data thus demonstrates that higher mutation rates of Mtb Beijing strains during human infection is likely to account for the higher adaptability and an emergence ratio of drug resistance.

Population Structure of Mycobacterium bovis in Germany: a Long-Term Study Using Whole-Genome Sequencing Combined with Conventional Molecular Typing Methods

Mycobacterium bovis is the primary cause of bovine tuberculosis (bTB) and infects a wide range of domestic animal and wildlife species and humans. In Germany, bTB still emerges sporadically in cattle herds, free-ranging wildlife, diverse captive animal species, and humans. In order to understand the underlying population structure and estimate the population size fluctuation through time, we analyzed 131 M. bovis strains from animals (n = 38) and humans (n = 93) in Germany from 1999 to 2017 by whole-genome sequencing (WGS), mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing, and spoligotyping. Based on WGS data analysis, 122 out of the 131 M. bovis strains were classified into 13 major clades, of which 6 contained strains from both human and animal cases and 7 only strains from human cases. Bayesian analyses suggest that the M. bovis population went through two sharp anticlimaxes, one in the middle of the 18th century and another one in the 1950s. WGS-based cluster analysis grouped 46 strains into 13 clusters ranging in size from 2 to 11 members and involving strains from distinct host types, e.g., only cattle and also mixed hosts. Animal strains of four clusters were obtained over a 9-year span, pointing toward autochthonous persistent bTB infection cycles. As expected, WGS had a higher discriminatory power than spoligotyping and MIRU-VNTR typing. In conclusion, our data confirm that WGS and suitable bioinformatics constitute the method of choice to implement prospective molecular epidemiological surveillance of M. bovis The population of M. bovis in Germany is diverse, with subtle, but existing, interactions between different host groups.

Microbe hunting in the modern era: reflecting on a decade of microbial genomic epidemiology

Since the first recognition that infectious microbes serve as the causes of many human diseases, physicians and scientists have sought to understand and control their spread. For the past 150+ years, these ‘microbe hunters’ have learned to combine epidemiological information with knowledge of the infectious agent(s). In this essay, I reflect on the evolution of microbe hunting, beginning with the history of pre-germ theory epidemiological studies, through the microbiological and molecular eras. Now in the genomic age, modern-day microbe hunters are combining pathogen whole-genome sequencing with epidemiological data to enhance epidemiological investigations, advance our understanding of the natural history of pathogens and drivers of disease, and ultimately reshape our plans and priorities for global disease control and eradication. Indeed, as we have seen during the ongoing Covid-19 pandemic, the role of microbe hunters is now more important than ever. Despite the advances already made by microbial genomic epidemiology, the field is still maturing, with many more exciting developments on the horizon.

The Biotin Biosynthetic Pathway in Mycobacterium tuberculosis is a Validated Target for the Development of Antibacterial Agents

Mycobacterium tuberculosis, responsible for Tuberculosis (TB), remains the leading cause of mortality among infectious diseases worldwide from a single infectious agent, with an estimated 1.7 million deaths in 2016. Biotin is an essential cofactor in M. tuberculosis that is required for lipid biosynthesis and gluconeogenesis. M. tuberculosis relies on de novo biotin biosynthesis to obtain this vital cofactor since it cannot scavenge sufficient biotin from a mammalian host. The biotin biosynthetic pathway in M. tuberculosis has been well studied and rigorously genetically validated providing a solid foundation for medicinal chemistry efforts. This review examines the mechanism and structure of the enzymes involved in biotin biosynthesis and ligation, summarizes the reported genetic validation studies of the pathway, and then analyzes the most promising inhibitors and natural products obtained from structure-based drug design and phenotypic screening.

Genetic Diversity in Mycobacterium tuberculosis Clinical Isolates and Resulting Outcomes of Tuberculosis Infection and Disease

Tuberculosis claims more human lives than any other bacterial infectious disease and represents a clear and present danger to global health as new tools for vaccination, treatment, and interruption of transmission have been slow to emerge. Additionally, tuberculosis presents with notable clinical heterogeneity, which complicates diagnosis, treatment, and the establishment of nonrelapsing cure. How this heterogeneity is driven by the diversity ofclinical isolates of the causative agent, Mycobacterium tuberculosis, has recently garnered attention. Herein, we review advances in the understanding of how naturally occurring variation in clinical isolates affects transmissibility, pathogenesis, immune modulation, and drug resistance. We also summarize how specific changes in transcriptional responses can modulate infection or disease outcome, together with strain-specific effects on gene essentiality. Further understanding of how this diversity of M. tuberculosis isolates affects disease and treatment outcomes will enable the development of more effective therapeutic options and vaccines for this dreaded disease.

Preclinical Progress of Subunit and Live Attenuated Mycobacterium tuberculosis Vaccines: A Review following the First in Human Efficacy Trial

Tuberculosis (TB) is the global leading cause of death from an infectious agent with approximately 10 million new cases of TB and 1.45 million deaths in 2018. Bacille Calmette-Guérin (BCG) remains the only approved vaccine for Mycobacterium tuberculosis (M. tb, causative agent of TB), however clinical studies have shown BCG has variable effectiveness ranging from 0–80% in adults. With 1.7 billion people latently infected, it is becoming clear that vaccine regimens aimed at both post-exposure and pre-exposure to M. tb will be crucial to end the TB epidemic. The two main strategies to improve or replace BCG are subunit and live attenuated vaccines. However, following the failure of the MVA85A phase IIb trial in 2013, more varied and innovative approaches are being developed. These include recombinant BCG strains, genetically attenuated M. tb and naturally attenuated mycobacteria strains, novel methods of immunogenic antigen discovery including for hypervirulent M. tb strains, improved antigen recognition and delivery strategies, and broader selection of viral vectors. This article reviews preclinical vaccine work in the last 5 years with focus on those tested against M. tb challenge in relevant animal models.

Diagnostic Potential of a PPE Protein Derived from Mycobacterium tuberculosis Beijing/K Strain

PURPOSE

The prevalence of Mycobacterium tuberculosis (M. tb) and the status of M. bovis BCG vaccination may affect host immune responses to M. tb antigens. Understanding of the predominant local M. tb strain and immune signatures induced by its strain-specific antigens may contribute to an improved diagnosis of tuberculosis (TB). The aim of this study was to determine immune responses to M. tb antigen which was identified from the hyper-virulent Beijing/K strain in South Korea.

MATERIALS AND METHODS

Pulmonary TB patients (n=52) and healthy subjects (n=92) including individuals with latent TB infection (n=31) were recruited, and QuantiFERON-TB Gold In-Tube tests were performed. The Beijing/K-antigen specific immune signatures were examined by diluted whole blood assays and multiplex bead arrays in a setting where nationwide BCG vaccination is employed.

RESULTS

Statistical analyses demonstrated that three [C-X-C motif chemokine (CXCL10), interleukin (IL)-6, interferon (IFN)-α] of 17 cytokines/chemokines distinguished active cases from healthy controls following stimulation with the Beijing/K-specific antigen. IFN-α also differentiated between active diseases and latent TB infection (p<0.01), and the detection rate of TB was dramatically increased in combination with IL-6 and CXCL10 at the highest levels of specificity (95-100%).

CONCLUSION

Our data indicate that immune signatures to the M. tb Beijing/K-specific antigen can provide useful information for improved TB diagnostics. The antigen may be developed as a diagnostic marker or a vaccine candidate, particularly in regions where the M. tb Beijing/K strain is endemic.

Insertion and deletion evolution reflects antibiotics selection pressure in a Mycobacterium tuberculosis outbreak

In genome evolution, genetic variants are the source of diversity, which natural selection acts upon. Treatment of human tuberculosis (TB) induces a strong selection pressure for the emergence of antibiotic resistance-conferring variants in the infecting Mycobacterium tuberculosis (MTB) strains. MTB evolution in response to treatment has been intensively studied and mainly attributed to point substitutions. However, the frequency and contribution of insertions and deletions (indels) to MTB genome evolution remains poorly understood. Here, we analyzed a multi-drug resistant MTB outbreak for the presence of high-quality indels and substitutions. We find that indels are significantly enriched in genes conferring antibiotic resistance. Furthermore, we show that indels are inherited during the outbreak and follow a molecular clock with an evolutionary rate of 5.37e-9 indels/site/year, which is 23 times lower than the substitution rate. Inherited indels may co-occur with substitutions in genes along related biological pathways; examples are iron storage and resistance to second-line antibiotics. This suggests that epistatic interactions between indels and substitutions affect antibiotic resistance and compensatory evolution in MTB.

Extensive Homoplasy but No Evidence of Convergent Evolution of Repeat Numbers at MIRU Loci in Modern Mycobacterium tuberculosis Lineages

More human deaths have been attributable to Mycobacterium tuberculosis than any other pathogen, and the epidemic is sustained by ongoing transmission. Various typing schemes have been developed to identify strain-specific differences and track transmission dynamics in affected communities, with recent introduction of whole genome sequencing providing the most accurate assessment. Mycobacterial interspersed repetitive unit (MIRU) typing is a family of variable number tandem repeat schemes that have been widely used to study the molecular epidemiology of M. tuberculosis. MIRU typing was used in most well-resourced settings to perform routine molecular epidemiology. Instances of MIRU homoplasy have been observed in comparison with sequence-based phylogenies, limiting its discriminatory value. A fundamental question is whether the observed homoplasy arises purely through stochastic processes, or whether there is evidence of natural selection. We compared repeat numbers at 24 MIRU loci with a whole genome sequence-based phylogeny of 245 isolates representing three modern M. tuberculosis lineages. This analysis demonstrated extensive homoplasy of repeat numbers, but did not detect any evidence of natural selection of repeat numbers, at least since the ancestral branching of the three modern lineages of M. tuberculosis. In addition, we observed good sensitivity but poor specificity and positive predictive values of MIRU-24 to detect clusters of recent transmission, as defined by whole-genome single nucleotide polymorphism analysis. These findings provide mechanistic insight, and support a transition away from VNTR-based typing toward sequence-based typing schemes for both research and public health purposes.

Cross-sectional Whole-genome Sequencing and Epidemiological Study of Multidrug-resistant Mycobacterium tuberculosis in China

Background

The increase in multidrug-resistant tuberculosis (MDR-TB) severely hampers tuberculosis prevention and control in China, a country with the second highest MDR-TB burden globally. The first nationwide drug-resistant tuberculosis surveillance program provides an opportunity to comprehensively investigate the epidemiological/drug-resistance characteristics, potential drug-resistance mutations, and effective population changes of Chinese MDR-TB.

Methods

We sequenced 357 MDR strains from 4600 representative tuberculosis-positive sputum samples collected during the survey (70 counties in 31 provinces). Drug-susceptibility testing was performed using 18 anti-tuberculosis drugs, representing the most comprehensive drug-resistance profile to date. We used 3 statistical and 1 machine-learning methods to identify drug-resistance genes/single-nucleotide polymorphisms (SNPs). We used Bayesian skyline analysis to investigate changes in effective population size.

Results

Epidemiological/drug-resistance characteristics showed different MDR profiles, co-resistance patterns, preferred drug combination/use, and recommended regimens among 7 Chinese administrative regions. These factors not only reflected the serious multidrug co-resistance and drug misuse but they were also potentially significant in facilitating the development of appropriate regimens for MDR-TB treatment in China. Further investigation identified 86 drug-resistance genes/intergenic regions/SNPs (58 new), providing potential targets for MDR-TB diagnosis and treatment. In addition, the effective population of Chinese MDR-TB displayed a strong expansion during 1993–2000, reflecting socioeconomic transition within the country. The phenomenon of expansion was restrained after 2000, likely attributable to the advances in diagnosis/treatment technologies and government support.

Conclusions

Our findings provide an important reference and improved understanding of MDR-TB in China, which are potentially significant in achieving the goal of precision medicine with respect to MDR-TB prevention and treatment.

Spoligotype and Drug Susceptibility Profiles of Mycobacterium tuberculosis Complex Isolates in Golestan Province, North Iran

Introduction

Despite the moderate incidence of tuberculosis (TB) in many parts of Iran, Golestan province had a permanently higher TB incidence rate than the national average. Moreover, Golestan province receives immigrants, mainly from TB-endemic areas of Iran and neighbor countries. Here, we aimed to characterize the circulating Mycobacterium tuberculosis complex (MTBC) isolates in terms of the spoligotype and drug resistance patterns, across Golestan province.

Materials and Methods

A set of 166 MTBC isolates was collected during July 2014 to July 2015 and subjected to drug susceptibility testing for first- and second-line anti-TB drugs and spoligotyping.

Results

Of 166 MTBC isolates, 139 (83.7%) isolates were assigned to 28 spoligotype international types (SITs). The most frequent SITs were SIT127/Ural-2 (n=25, 15.1%), followed by SIT1/Beijing (n=21, 12.7%) and SIT3427/Ural-2 (n=18, 10.8%). The set of 18 isolates (10.8%) showed resistance to at least one drug, which mainly belonged to SIT1/Beijing (n=7, 38.9%), orphan patterns (n=4, 22.2%) and SIT357/CAS1-Delhi (n=3, 16.7%). In addition, four isolates (2.4%) were resistant to pyrazinamide. The analysis of mutation corresponded to resistance to rifampin and isoniazid showed that two isolates had Ser531Leu substitution in rpoB, four isolates had Ser315Thr substitution in katG and one isolate had [C(−15)T] in inhA locus.

Conclusion

High diversity in spoligotypes of the MTBC isolates and lack of dominant genotype might be due to residence of immigrants in this region and consequent reactivation of latent infection. In addition, due to the presence of extensively drug-resistant (XDR) isolates in Golestan province, it is important to conduct future studies to determine transmission pattern of drug-resistant isolates in this region.

Targeted surveillance strategies for efficient detection of novel antibiotic resistance variants

Genotype-based diagnostics for antibiotic resistance represent a promising alternative to empiric therapy, reducing inappropriate antibiotic use. However, because such assays infer resistance based on known genetic markers, their utility will wane with the emergence of novel resistance. Maintenance of these diagnostics will therefore require surveillance to ensure early detection of novel resistance variants, but efficient strategies to do so remain undefined. We evaluate the efficiency of targeted sampling approaches informed by patient and pathogen characteristics in detecting antibiotic resistance and diagnostic escape variants in Neisseria gonorrhoeae, a pathogen associated with a high burden of disease and antibiotic resistance and the development of genotype-based diagnostics. We show that patient characteristic-informed sampling is not a reliable strategy for efficient variant detection. In contrast, sampling informed by pathogen characteristics, such as genomic diversity and genomic background, is significantly more efficient than random sampling in identifying genetic variants associated with resistance and diagnostic escape.

Local adaptive evolution of two distinct clades of Beijing and T families of Mycobacterium tuberculosis in Chongqing: a Bayesian population structure and phylogenetic study

Background

Beijing sub-pedigree 2 (BSP2) and T sub-lineage 6 (TSL6) are two clades belonging to Beijing and T family of Mycobacterium tuberculosis (MTB), respectively, defined by Bayesian population structure analysis based on 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTR). Globally, over 99% of BSP2 and 89% of TSL6 isolates were distributed in Chongqing, suggesting their possible local adaptive evolution. The objective of this paper is to explore whether BSP2 and TSL6 originated by their local adaptive evolution from the specific isolates of Beijing and T families in Chongqing.

Methods

The genotyping data of 16 090 MTB isolates were collected from laboratory collection, published literatures and SITVIT database before subjected to Bayesian population structure analysis based on 24-loci MIRU-VNTR. Spacer Oligonucleotide Forest (Spoligoforest) and 24-loci MIRU-VNTR-based minimum spanning tree (MST) were used to explore their phylogenetic pathways, with Bayesian demographic analysis for exploring the recent demographic change of TSL6.

Results

Phylogenetic analysis suggested that BSP2 and TSL6 in Chongqing may evolve from BSP4 and TSL5, respectively, which were locally predominant in Tibet and Jiangsu, respectively. Spoligoforest showed that Beijing and T families were genetically distant, while the convergence of the MIRU-VNTR pattern of BSP2 and TSL6 was revealed by WebLogo. The demographic analysis concluded that the recent demographic change of TSL6 might take 111.25 years.

Conclusions

BSP2 and TSL6 clades might originate from BSP4 and TSL5, respectively, by their local adaptive evolution in Chongqing. Our study suggests MIRU-VNTR be combined with other robust markers for a more comprehensive genotyping approach, especially for families of clades with the same MIRU-VNTR pattern.

Applied phyloepidemiology: Detecting drivers of pathogen transmission from genomic signatures using density measures

Understanding the driving forces of an epidemic is key to inform intervention strategies against it. Correlating measures of the epidemic success of a pathogen with ancillary parameters such as its drug resistance profile provides a flexible tool to identify such driving forces. The recently described time‐scaled haplotypic density (THD) method facilitates the inference of a pathogen's epidemic success from genetic data. Contrary to demogenetic approaches that define success in an aggregated fashion, the THD computes an independent index of success for each isolate in a collection. Modeling this index using multivariate regression, thus, allows us to control for various sources of bias and to identify independent predictors of success. We illustrate the use of THD to address key questions regarding three exemplary epidemics of multidrug‐resistant (MDR) bacterial lineages, namely Mycobacterium tuberculosis Beijing, Salmonella Typhi H58, and Staphylococcus aureus ST8 (including ST8‐USA300 MRSA), based on previously published, international genetic datasets. In each case, THD analysis allowed to identify the impact, or lack thereof, of various factors on the epidemic success, independent of confounding by population structure and geographic distribution. Our results suggest that rifampicin resistance drives the MDR Beijing epidemic and that fluoroquinolone resistance drives the S. aureus ST8/USA300 epidemic, in line with previous evidence of a lack of resistance‐associated fitness cost in these pathogens. Conversely, fluoroquinolone resistance measurably hampered the success of S. Typhi H58 and non‐H58. These findings illustrate how THD can help leverage the massive genomic datasets generated by molecular epidemiology studies to address new questions. THD implementation for the R platform is available at https://github.com/rasigadelab/thd.

Evidence for the Effect of Vaccination on Host-Pathogen Interactions in a Murine Model of Pulmonary Tuberculosis by Mycobacterium tuberculosis

The global control of Tuberculosis remains elusive, and Bacillus Calmette-Guérin (BCG) -the most widely used vaccine in history-has proven insufficient for reversing this epidemic. Several authors have suggested that the mass presence of vaccinated hosts might have affected the Mycobacterium tuberculosis (MTB) population structure, and this could in turn be reflected in a prevalence of strains with higher ability to circumvent BCG-induced immunity, such as the recent Beijing genotype. The effect of vaccination on vaccine-escape variants has been well-documented in several bacterial pathogens; however the effect of the interaction between MTB strains and vaccinated hosts has never been previously described. In this study we show for the first time the interaction between MTB Beijing-genotype strains and BCG-vaccinated hosts. Using a well-controlled murine model of progressive pulmonary tuberculosis, we vaccinated BALB/c mice with two different sub-strains of BCG (BCG-Phipps and BCG-Vietnam). Following vaccination, the mice were infected with either one of three selected MTB strains. Strains were selected based on lineage, and included two Beijing-family clinical isolates (strains 46 and 48) and a well-characterized laboratory strain (H37Rv). Two months after infection, mice were euthanized and the bacteria extracted from their lungs. We characterized the genomic composite of the bacteria before and after exposure to vaccinated hosts, and also characterized the local response to the bacteria by sequencing the lung transcriptome in animals during the infection. Results from this study show that the interaction within the lungs of the vaccinated hosts results in the selection of higher-virulence bacteria, specifically for the Beijing genotype strains 46 and 48. After exposure to the BCG-induced immune response, strains 46 and 48 acquire genomic mutations associated with several virulence factors. As a result, the bacteria collected from these vaccinated hosts have an increased ability for immune evasion, as shown in both the host transcriptome and the histopathology studies, and replicates far more efficiently compared to bacteria collected from unvaccinated hosts or to the original-stock strain. Further research is warranted to ascertain the pathways associated with the genomic alterations. However, our results highlight novel host-pathogen interactions induced by exposure of MTB to BCG vaccinated hosts.

Epidemiology characteristics of the clonal complexes of Mycobacterium tuberculosis Lineage 4 in China

Mycobacterium tuberculosis (M. tuberculosis) Lineage 4 (L4) is frequently prevailing in Western regions of China, where the tuberculosis incidence rate is high. However, the epidemiology characteristics of M. tuberculosis L4 in China remain poorly understood. Here, the 15-loci Variable number of tandem repeats (VNTR) patterns of 975 L4 isolates from a National Survey of Tuberculosis in China were used to construct a Minimum Spanning Tree (MST), which divided the 975 isolates into 5 major clonal complexes (CC; named CC1 to CC5). We found that the CCs of M. tuberculosis L4 were nationally distributed, geographically restricted, and different in epidemiology characteristics. For example, CC1 was mainly concentrated in East and Central China and significantly related to the farmer occupation and income of an individual (>4200 yuan) (p < .05); CC5 was mainly distributed in Southwest China and was associated with ethnic minorities. Notably, using whole genome sequencing (WGS) data of 141 strains that matched our samples, we found that both CC1 and CC5 were mapped to the sublineage L4.5. Nevertheless, due to the difference of geographical distribution, the epidemiology characteristics of these CCs were largely different. We found that income and occupation significantly contributed to the odds of infection by CC1 to CC5. Consequently, our findings revealed the epidemiology characteristics of the CCs of M. tuberculosis L4, and will help in the formulation of more effective intervention measures in line with regional specifications and patient characteristics in China.

Epidemic of venereal treponematosis in wild monkeys: a paradigm for syphilis origin

Treponema pallidum infections have been primarily known as slightly contagious mucocutaneous infections called yaws (tropical Africa and America) and bejel (subtropical North Africa). T. pallidum emerged as a highly infectious venereal syphilis agent in South America, probably about 500 years ago, and because of its venereal transmission, it quickly caused a worldwide pandemic. The disease manifests as lesions, including a chancre; then antibodies become detectable when or slightly after the chancre appears, and before the development of a rash and other systemic manifestations. Venereal diseases are poorly known in monkeys. During fieldwork in Senegal, we discovered an epizootic outbreak of venereal disease that we explored. We detected a venereal form of T. pallidum subsp. pertenue infection in green monkeys (Chlorocebus sabaeus), then observed an epizootic outbreak in Senegal and its spread among baboons a year later. Comparative analysis of T. pallidum genomes from the monkeys' chancres and other Treponema genomes showed an acceleration of the number of single nucleotide polymorphisms, comparable to that observed in syphilis. Identified T. pallidum clones seem to be epizootic through the acceleration of their mutation rate, which is linked to their larger diffusion.

Genotyping, Assessment of Virulence and Antibacterial Resistance of the Rostov Strain of Mycobacterium tuberculosis Attributed to the Central Asia Outbreak Clade

The Central Asia Outbreak (CAO) clade is a growing public health problem for Central Asian countries. Members of the clade belong to the narrow branch of the Mycobacterium tuberculosis Beijing genotype and are characterized by multidrug resistance and increased transmissibility. The Rostov strain of M. tuberculosis isolated in Russia and attributed to the CAO clade based on PCR-assay and whole genome sequencing and the laboratory strain H37Rv were selected to evaluate the virulence on C57Bl/6 mice models by intravenous injection. All mice infected with the Rostov strain succumbed to death within a 48-day period, while more than half of the mice infected by the H37Rv strain survived within a 90-day period. Mice weight analysis revealed irreversible and severe depletion of animals infected with the Rostov strain compared to H37Rv. The histological investigation of lung and liver tissues of mice on the 30th day after injection of mycobacterial bacilli showed that the pattern of pathological changes generated by two strains were different. Moreover, bacterial load in the liver and lungs was higher for the Rostov strain infection. In conclusion, our data demonstrate that the drug-resistant Rostov strain exhibits a highly virulent phenotype which can be partly explained by the CAO-specific mutations.

Live attenuated TB vaccines representing the three modern Mycobacterium tuberculosis lineages reveal that the Euro-American genetic background confers optimal vaccine potential

Background

Human tuberculosis (TB) is caused by a plethora of Mycobacterium tuberculosis complex (MTBC) strains belonging to seven phylogenetic branches. Lineages 2, 3 and 4 are considered “modern” branches of the MTBC responsible for the majority of worldwide TB. Since the current BCG vaccine confers variable protection against pulmonary TB, new candidates are investigated. MTBVAC is the unique live attenuated vaccine based on M. tuberculosis in human clinical trials.

Methods

MTBVAC was originally constructed by unmarked phoP and fadD26 deletions in a clinical isolate belonging to L4. Here we construct new vaccines based on isogenic gene deletions in clinical isolates of the L2 and L3 modern lineages. These three vaccine candidates were characterized at molecular level and also in animal experiments of protection and safety.

Findings

Safety studies in immunocompromised mice showed that MTBVAC-L2 was less attenuated than BCG Pasteur, while the original MTBVAC was found even more attenuated than BCG and MTBVAC-L3 showed an intermediate phenotype. The three MTBVAC candidates showed similar or superior protection compared to BCG in immunocompetent mice vaccinated with each MTBVAC candidate and challenged with three representative strains of the modern lineages.

Interpretation

MTBVAC vaccines, based on double phoP and fadD26 deletions, protect against TB independently of the phylogenetic linage used as template strain for their construction. Nevertheless, lineage L4 confers the best safety profile.

Funding

European Commission (TBVAC2020, H2020-PHC-643381), Spanish Ministry of Science (RTI2018-097625-B-I00), Instituto de Salud Carlos III (PI18/0336), Gobierno de Aragón/Fondo Social Europeo and the French National Research Council (ANR-10-LABX-62-IBEID, ANR-16-CE35-0009, ANR-16-CE15-0003).

Niche specialization and spread of Staphylococcus capitis involved in neonatal sepsis

The multidrug-resistant Staphylococcus capitis NRCS-A clone is responsible for sepsis in preterm infants in neonatal intensive care units (NICUs) worldwide. Here, to retrace the spread of this clone and to identify drivers of its specific success, we investigated a representative collection of 250 S. capitis isolates from adults and newborns. Bayesian analyses confirmed the spread of the NRCS-A clone and enabled us to date its emergence in the late 1960s and its expansion during the 1980s, coinciding with the establishment of NICUs and the increasing use of vancomycin in these units, respectively. This dynamic was accompanied by the acquisition of mutations in antimicrobial resistance- and bacteriocin-encoding genes. Furthermore, combined statistical tools and a genome-wide association study convergently point to vancomycin resistance as a major driver of NRCS-A success. We also identified another S. capitis subclade (alpha clade) that emerged independently, showing parallel evolution towards NICU specialization and non-susceptibility to vancomycin, indicating convergent evolution in NICU-associated pathogens. These findings illustrate how the broad use of antibiotics can repeatedly lead initially commensal drug-susceptible bacteria to evolve into multidrug-resistant clones that are able to successfully spread worldwide and become pathogenic for highly vulnerable patients.

Increased male susceptibility to Mycobacterium tuberculosis infection is associated with smaller B cell follicles in the lungs

Tuberculosis prevalence is significantly higher among men than women. We have previously revealed an increased susceptibility of male C57BL/6 mice towards Mycobacterium tuberculosis (Mtb) H37Rv. In the current study, we confirm the male bias for infection with the Beijing strain HN878. Males succumbed to HN878 infection significantly earlier than females. In both models, premature death of males was associated with smaller B cell follicles in the lungs. Analysis of homeostatic chemokines and their receptors revealed differences between H37Rv and HN878 infected animals, indicating different immune requirements for follicle formation in both models. However, expression of IL-23, which is involved in long-term containment of Mtb and lymphoid follicle formation, was reduced in male compared to female lungs in both models. Our study reveals sex differences in the formation of B cell follicles in the Mtb infected lung and we propose that impaired follicle formation is responsible for accelerated disease progression in males.

Phylogenetically informative mutations in genes implicated in antibiotic resistance in Mycobacterium tuberculosis complex

BACKGROUND

A comprehensive understanding of the pre-existing genetic variation in genes associated with antibiotic resistance in the Mycobacterium tuberculosis complex (MTBC) is needed to accurately interpret whole-genome sequencing data for genotypic drug susceptibility testing (DST).

METHODS

We investigated mutations in 92 genes implicated in resistance to 21 anti-tuberculosis drugs using the genomes of 405 phylogenetically diverse MTBC strains. The role of phylogenetically informative mutations was assessed by routine phenotypic DST data for the first-line drugs isoniazid, rifampicin, ethambutol, and pyrazinamide from a separate collection of over 7000 clinical strains. Selected mutations/strains were further investigated by minimum inhibitory concentration (MIC) testing.

RESULTS

Out of 547 phylogenetically informative mutations identified, 138 were classified as not correlating with resistance to first-line drugs. MIC testing did not reveal a discernible impact of a Rv1979c deletion shared by M. africanum lineage 5 strains on resistance to clofazimine. Finally, we found molecular evidence that some MTBC subgroups may be hyper-susceptible to bedaquiline and clofazimine by different loss-of-function mutations affecting a drug efflux pump subunit (MmpL5).

CONCLUSIONS

Our findings underline that the genetic diversity in MTBC has to be studied more systematically to inform the design of clinical trials and to define sound epidemiologic cut-off values (ECOFFs) for new and repurposed anti-tuberculosis drugs. In that regard, our comprehensive variant catalogue provides a solid basis for the interpretation of mutations in genotypic as well as in phenotypic DST assays.

Differential determinants of virulence in two Mycobacterium tuberculosis Colombian clinical isolates of the LAM09 family

The heterogeneity of the clinical outcome of Mycobacterium tuberculosis (Mtb) infection may be due in part to different strategies used by circulating strains to cause disease. This heterogeneity is one of the main limitations to eradicate tuberculosis disease. In this study, we have compared the transcriptional response of two closely related Colombian clinical isolates (UT127 and UT205) of the LAM family under two axenic media conditions. These clinical isolates are phenotypically different at the level of cell death, cytokine production, growth kinetics upon in vitro infection of human tissue macrophages, and membrane vesicle secretion upon culture in synthetic medium. Using RNA-seq, we have identified different pathways that account for two different strategies to cope with the stressful condition of a carbon-poor media such as Sauton’s. We showed that the clinical isolate UT205 focus mainly in the activation of virulence systems such as the ESX-1, synthesis of diacyl-trehalose, polyacyl-trehalose, and sulfolipids, while UT127 concentrates its efforts mainly in the survival mode by the activation of the DNA replication, cell division, and lipid biosynthesis. This is an example of two Mtb isolates that belong to the same family and lineage, and even though they have a very similar genome, its transcriptional regulation showed important differences. This results in summary highlight the necessity to reach a better understanding of the heterogeneity in the behavior of these circulating Mtb strains which may help us to design better treatments and vaccines and to identify new targets for drugs.

Whole Genome Sequencing Results Associated with Minimum Inhibitory Concentrations of 14 Anti-Tuberculosis Drugs among Rifampicin-Resistant Isolates of Mycobacterium Tuberculosis from Iran

Accurate and timely detection of drug resistance can minimize the risk of further resistance development and lead to effective treatment. The aim of this study was to determine the resistance to first/second-line anti-tuberculosis drugs in rifampicin/multidrug-resistant Mycobacterium tuberculosis (RR/MDR-MTB) isolates. Molecular epidemiology of strains was determined using whole genome sequencing (WGS)-based genotyping. A total of 35 RR/MDR-MTB isolates were subjected to drug susceptibility testing against first/second-line drugs using 7H9 Middlebrook in broth microdilution method. Illumina technology was used for paired-end WGS applying a Maxwell 16 Cell DNA Purification kit and the NextSeq platform. Data analysis and single nucleotide polymorphism calling were performed using MTBseq pipeline. The genome-based resistance to each drug among the resistant phenotypes was as follows: rifampicin (97.1%), isoniazid (96.6%), ethambutol (100%), levofloxacin (83.3%), moxifloxacin (83.3%), amikacin (100%), kanamycin (100%), capreomycin (100%), prothionamide (100%), D-cycloserine (11.1%), clofazimine (20%), bedaquiline (0.0%), and delamanid (44.4%). There was no linezolid-resistant phenotype, and a bedaquiline-resistant strain was wild type for related genes. The Beijing, Euro-American, and Delhi-CAS were the most populated lineage/sublineages. Drug resistance-associated mutations were mostly linked to minimum inhibitory concentration results. However, the role of well-known drug-resistant genes for D-cycloserine, clofazimine, bedaquiline, and delamanid was found to be more controversial.

A landscape of genomic alterations at the root of a near-untreatable tuberculosis epidemic

BACKGROUND

Atypical Beijing genotype Mycobacterium tuberculosis strains are widespread in South Africa and have acquired resistance to up to 13 drugs on multiple occasions. It is puzzling that these strains have retained fitness and transmissibility despite the potential fitness cost associated with drug resistance mutations.

METHODS

We conducted Illumina sequencing of 211 Beijing genotype M. tuberculosis isolates to facilitate the detection of genomic features that may promote acquisition of drug resistance and restore fitness in highly resistant atypical Beijing forms. Phylogenetic and comparative genomic analysis was done to determine changes that are unique to the resistant strains that also transmit well. Minimum inhibitory concentration (MIC) determination for streptomycin and bedaquiline was done for a limited number of isolates to demonstrate a difference in MIC between isolates with and without certain variants.

RESULTS

Phylogenetic analysis confirmed that two clades of atypical Beijing strains have independently developed resistance to virtually all the potent drugs included in standard (pre-bedaquiline) drug-resistant TB treatment regimens. We show that undetected drug resistance in a progenitor strain was likely instrumental in this resistance acquisition. In this cohort, ethionamide (ethA A381P) resistance would be missed in first-line drug-susceptible isolates, and streptomycin (gidB L79S) resistance may be missed due to an MIC close to the critical concentration. Subsequent inadequate treatment historically led to amplification of resistance and facilitated spread of the strains. Bedaquiline resistance was found in a small number of isolates, despite lack of exposure to the drug. The highly resistant clades also carry inhA promoter mutations, which arose after ethA and katG mutations. In these isolates, inhA promoter mutations do not alter drug resistance, suggesting a possible alternative role.

CONCLUSION

The presence of the ethA mutation in otherwise susceptible isolates from ethionamide-naïve patients demonstrates that known exposure is not an adequate indicator of drug susceptibility. Similarly, it is demonstrated that bedaquiline resistance can occur without exposure to the drug. Inappropriate treatment regimens, due to missed resistance, leads to amplification of resistance, and transmission. We put these results into the context of current WHO treatment regimens, underscoring the risks of treatment without knowledge of the full drug resistance profile.

The Genetic Background Modulates the Evolution of Fluoroquinolone-Resistance in Mycobacterium tuberculosis

Fluoroquinolones (FQ) form the backbone in experimental treatment regimens against drug-susceptible tuberculosis. However, little is known on whether the genetic variation present in natural populations of Mycobacterium tuberculosis (Mtb) affects the evolution of FQ-resistance (FQ-R). To investigate this question, we used nine genetically distinct drug-susceptible clinical isolates of Mtb and measured their frequency of resistance to the FQ ofloxacin (OFX) in vitro. We found that the Mtb genetic background led to differences in the frequency of OFX-resistance (OFX-R) that spanned two orders of magnitude and substantially modulated the observed mutational profiles for OFX-R. Further, in vitro assays showed that the genetic background also influenced the minimum inhibitory concentration and the fitness effect conferred by a given OFX-R mutation. To test the clinical relevance of our in vitro work, we surveyed the mutational profile for FQ-R in publicly available genomic sequences from clinical Mtb isolates, and found substantial Mtb lineage-dependent variability. Comparison of the clinical and the in vitro mutational profiles for FQ-R showed that 51% and 39% of the variability in the clinical frequency of FQ-R gyrA mutation events in Lineage 2 and Lineage 4 strains, respectively, can be attributed to how Mtb evolves FQ-R in vitro. As the Mtb genetic background strongly influenced the evolution of FQ-R in vitro, we conclude that the genetic background of Mtb also impacts the evolution of FQ-R in the clinic.

TbD1 deletion as a driver of the evolutionary success of modern epidemic Mycobacterium tuberculosis lineages

Mycobacterium tuberculosis (Mtb) strains are classified into different phylogenetic lineages (L), three of which (L2/L3/L4) emerged from a common progenitor after the loss of the MmpS6/MmpL6-encoding Mtb-specific deletion 1 region (TbD1). These TbD1-deleted “modern” lineages are responsible for globally-spread tuberculosis epidemics, whereas TbD1-intact “ancestral” lineages tend to be restricted to specific geographical areas, such as South India and South East Asia (L1) or East Africa (L7). By constructing and characterizing a panel of recombinant TbD1-knock-in and knock-out strains and comparison with clinical isolates, here we show that deletion of TbD1 confers to Mtb a significant increase in resistance to oxidative stress and hypoxia, which correlates with enhanced virulence in selected cellular, guinea pig and C3HeB/FeJ mouse infection models, the latter two mirroring in part the development of hypoxic granulomas in human disease progression. Our results suggest that loss of TbD1 at the origin of the L2/L3/L4 Mtb lineages was a key driver for their global epidemic spread and outstanding evolutionary success.

Mycobacterium tuberculosis (Mtb) modern strains emerged from a common progenitor after the loss of Mtb-specific deletion 1 region (TbD1). Here, the authors show that deletion of TbD1 correlates with enhanced Mtb virulence in animal models, mirroring the development of hypoxic granulomas in human disease progression.

The epidemic of multidrug resistant tuberculosis in China in historical and phylogenetic perspectives

OBJECTIVES

For the past decade, the epidemic of multidrug resistance tuberculosis (MDR-TB) stays high in China. We investigated the possible driving forces behind the epidemics from phylogenetic and historical perspectives.

METHODS

420 representative strains were selected from the first national drug resistance survey based on their genotypes, drug susceptibility patterns and geographic information. We reconstructed the phylogeny by whole genome sequencing and compared it to the global phylogeny including MDR outbreaks reported in other settings. We estimated the historical trajectory of population dynamics by Bayesian Skygrid plot for all strains and MDR-TB alone. Integrating geographic information and mutations in drug resistance related genes, we investigated the spatial scale of transmission, recent selection of drug resistant mutant, and mechanism for fitness restoration.

RESULTS

Three new subgroups within Beijing clade are described for the first time, but none of the MDR-TB outbreak strains reported in other high MDR-TB burden settings is identified. The overall epidemics experienced two successive phases of expansion at different rates between 1660s and 1950s, followed by a sharp decline till today. Four fifths of the clustered MDR-TB strains suggest transmission of DR strains and nearly half suggest recent selection of (additional) mutations in rpoB. Among all identified transmission events, about one fifth occurred between far distant locations. Possible intergenic and intragenic compensatory mutations both presented in our dataset at comparable frequencies.

CONCLUSIONS

MDR-TB epidemic in China is not yet driven by the spread of a few highly successful clonal expansions but by repeated emergence of smaller and currently less successful clusters. However, internal migration and undertreatment could escalate MDR-TB epidemic. To prevent generating of drug resistance and restoration of fitness as well as to stop transmission of MDR-TB at early stage, national TB control program needs to strengthen management of floating populations and promote universal drug susceptibility testing in China.

Diverse Clinical Isolates of Mycobacterium tuberculosis Develop Macrophage-Induced Rifampin Tolerance

The Mycobacterium tuberculosis lineage 4 strains CDC1551 and H37Rv develop tolerance to multiple antibiotics upon macrophage residence. To determine whether macrophage-induced tolerance is a general feature of clinical M. tuberculosis isolates, we assessed macrophage-induced drug tolerance in strains from lineages 1–3, representing the other predominant M. tuberculosis strains responsible for tuberculosis globally. All 3 lineages developed isoniazid tolerance. While lineage 1, 3, and 4 strains developed rifampin tolerance, lineage 2 Beijing strains did not. Their failure to develop tolerance may be explained by their harboring of a loss-of-function mutation in the Rv1258c efflux pump that is linked to macrophage-induced rifampicin tolerance.

Tuberculosis in Liberia: high multidrug-resistance burden, transmission and diversity modelled by multiple importation events

Tuberculosis (TB) surveillance is scarce in most African countries, even though it is the continent with the greatest disease incidence according to the World Health Organization. Liberia is within the 30 countries with the highest TB burden, probably as a consequence of the long civil war and the recent Ebola outbreak, both crippling the health system and depreciating the TB prevention and control programmes. Due to difficulties working in the country, there is a lack of resistance surveys and bacillus characterization. Here, we use genome sequencing of Mycobacteriumtuberculosis clinical isolates to fill this gap. Our results highlight that the bacillus population structure is dominated by lineage 4 strains that harbour an outstanding genetic diversity, higher than in the rest of Africa as a whole. Coalescent analyses demonstrate that strains currently circulating in Liberia were introduced several times beginning in the early year 600 CE until very recently coinciding with migratory movements associated with the civil war and Ebola epidemics. A higher multidrug-resistant (MDR)-TB frequency (23.5 %) than current estimates was obtained together with non-catalogued drug-resistance mutations. Additionally, 39 % of strains were in genomic clusters revealing that ongoing transmission is a major contribution to the TB burden in the country. Our report emphasizes the importance of TB surveillance and control in African countries where bacillus diversity, MDR-TB prevalence and transmission are coalescing to jeopardize TB control programmes.

New insights into the mycobacterial PE and PPE proteins provide a framework for future research

The PE and PPE proteins of Mycobacterium tuberculosis have been studied with great interest since their discovery. Named after the conserved proline (P) and glutamic acid (E) residues in their N-terminal domains, these proteins are postulated to perform wide-ranging roles in virulence and immune modulation. However, technical challenges in studying these proteins and their encoding genes have hampered the elucidation of molecular mechanisms and leave many open questions regarding the biological functions mediated by these proteins. Here, I review the shared and unique characteristics of PE and PPE proteins from a molecular perspective linking this information to their functions in mycobacterial virulence. I discuss how the different subgroups (PE_PGRS, PPE-PPW, PPE-SVP and PPE-MPTR) are defined and why this classification of paramount importance to understand the PE and PPE proteins as individuals and or groups. The goal of this MicroReview is to summarize and structure the existing information on this gene family into a simplified framework of thinking about PE and PPE proteins and genes. Thereby, I hope to provide helpful starting points in studying these genes and proteins for researchers with different backgrounds. This has particular implications for the design and monitoring of novel vaccine candidates and in understanding the evolution of the M. tuberculosis complex.

Mycobacterium tuberculosis clinical isolates of the Beijing and East-African Indian lineage induce fundamentally different host responses in mice compared to H37Rv

Substantial differences exist in virulence among Mycobacterium tuberculosis strains in preclinical TB models. In this study we show how virulence affects host responses in mice during the first four weeks of infection with a mycobacterial strain belonging to the Beijing, East-African-Indian or Euro-American lineage. BALB/c mice were infected with clinical isolates of the Beijing-1585 strain or the East-African Indian (EAI)-1627 strain and host responses were compared to mice infected with the non-clinical H37Rv strain of the Euro-American lineage. We found that H37Rv induced a ‘classical’ T-cell influx with high IFN-γ levels, while Beijing-1585 and EAI-1627 induced an influx of B-cells into the lungs together with elevated pulmonary IL-4 protein levels. Myeloid cells in the lungs appeared functionally impaired upon infection with Beijing-1585 and EAI-1627 with reduced iNOS and IL-12 expression levels compared to H37Rv infection. This impairment might be related to significantly reduced expression in the bone marrow of IFN-γ, TNF-α and IFN-β in mice infected with Beijing-1585 and EAI-1627, which could be detected from the third day post infection onwards. Our findings suggest that increased virulence of two clinical isolates compared to H37Rv is associated with a fundamentally different systemic immune response, which already can be detected early during infection.