Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage

Characterizing in vivo loss of virulence of an HN878 Mycobacterium tuberculosis isolate from a genetic duplication event

The increase of global cases of drug resistant (DR) Mycobacterium tuberculosis (M.tb) is a serious problem for the tuberculosis research community and the goals to END TB by 2030. Due to the need for advancing and screening next generation therapeutics and vaccines, we aimed to design preclinical DR models of Beijing lineage M.tb HN878 strain in different mouse backgrounds. We found escalating sensitivities of morbidity due to low dose aerosol challenge (50-100 bacilli) in CB6F1, C57BL/6 and SWR mice, respectively. We also observed that pulmonary bacterial burden at morbidity endpoints correlated inversely with survival over time between mouse strains. Interestingly, with in vitro passaging and in the process of selecting individual DR mutant colonies, we observed a significant decrease in in vivo HN878 strain virulence, which correlated with the acquisition of a large genetic duplication. We confirmed that low passage infection stocks with no or low prevalence of the duplication, including stocks directly acquired from the BEI resources biorepository, retained virulence, measured by morbidity over time. These data help confirm previous reports and emphasize the importance of monitoring virulence and stock fidelity.

A multi-omics investigation into the mechanisms of hyper-virulence in Mycobacterium tuberculosis

Clinical manifestations of tuberculosis range from asymptomatic infection to a life-threatening disease such as tuberculous meningitis (TBM). Recent studies showed that the spectrum of disease severity could be related to genetic diversity among clinical strains of Mycobacterium tuberculosis (Mtb). Certain strains are reported to preferentially invade the central nervous system, thus earning the label “hypervirulent strains”.However, specific genetic mutations that accounted for enhanced mycobacterial virulence are still unknown. We previously identified a set of 17 mutations in a hypervirulent Mtb strain that was from TBM patient and exhibited significantly better intracellular survivability. These mutations were also commonly shared by a cluster of globally circulating hyper-virulent strains. Here, we aimed to validate the impact of these hypervirulent-specific mutations on the dysregulation of gene networks associated with virulence in Mtb via multi-omic analysis. We surveyed transcriptomic and proteomic differences between the hyper-virulent and low-virulent strains using RNA-sequencing and label-free quantitative LC-MS/MS approach, respectively. We identified 25 genes consistently differentially expressed between the strains at both transcript and protein level, regardless the strains were growing in a nutrient-rich or a physiologically relevant multi-stress condition (acidic pH, limited nutrients, nitrosative stress, and hypoxia). Based on integrated genomic-transcriptomic and proteomic comparisons, the hypervirulent-specific mutations in FadE5 (g. 295,746 C >T), Rv0178 (p. asp150glu), higB (p. asp30glu), and pip (IS6110-insertion) were linked to deregulated expression of the respective genes and their functionally downstream regulons. The result validated the connections between mutations, gene expression, and mycobacterial pathogenicity, and identified new possible virulence-associated pathways in Mtb.

High Throughput Sequencing for Clinical Tuberculosis: An Overview

High throughput sequencing (HTS) can identify the presence of Mycobacterium tuberculosis DNA in a clinical sample while also providing information on drug susceptibility. Multiple studies have provided a context for exploring the clinical application of HTS for TB diagnosis. The workflow challenges, strengths and limitations of the various sequencing platforms, and tools used for analysis are presented to provide a framework for further innovations in the field.

Whole-Genome Sequencing for Resistance Prediction and Transmission Analysis of Mycobacterium tuberculosis Complex Strains from Namibia

Namibia is among 30 countries with a high burden of tuberculosis (TB), with an estimated incidence of 460 per 100,000 population and around 800 new multidrug-resistant (MDR) TB cases per year. Still, data on the transmission and evolution of drug-resistant Mycobacterium tuberculosis complex (Mtbc) strains are not available. Whole-genome sequencing data of 136 rifampicin-resistant (RIFr) Mtbc strains obtained from 2016 to 2018 were used for phylogenetic classification, resistance prediction, and cluster analysis and linked with phenotypic drug susceptibility testing (pDST) data. Roughly 50% of the strains investigated were resistant to all first-line drugs. Furthermore, 13% of the MDR Mtbc strains were already pre-extensively drug resistant (pre-XDR). The cluster rates were high, at 74.6% among MDR and 85% among pre-XDR strains. A significant proportion of strains had borderline resistance-conferring mutations, e.g., inhA promoter mutations or rpoB L430P. Accordingly, 25% of the RIFr strains tested susceptible by pDST. Finally, we determined a potentially new bedaquiline resistance mutation (Rv0678 D88G) occurring in two independent clusters. High rates of resistance to first-line drugs in line with emerging pre-XDR and likely bedaquiline resistance linked with the ongoing recent transmission of MDR Mtbc clones underline the urgent need for the implementation of interventions that allow rapid diagnostics to break MDR TB transmission chains in the country. A borderline RIFr mutation in the dominant outbreak strain causing discrepancies between phenotypic and genotypic resistance testing results may require breakpoint adjustments but also may allow individualized regimens with high-dose treatment. IMPORTANCE The transmission of drug-resistant tuberculosis (TB) is a major problem for global TB control. Using genome sequencing, we showed that 13% of the multidrug-resistant (MDR) M. tuberculosis complex strains from Namibia are already pre-extensively drug resistant (pre-XDR), which is substantial in an African setting. Our data also indicate that the ongoing transmission of MDR and pre-XDR strains contributes significantly to the problem. In contrast to other settings with higher rates of drug resistance, we found a high proportion of strains having so-called borderline low-level resistance mutations, e.g., inhA promoter mutations or rpoB L430P. This led to the misclassification of 25% of the rifampicin-resistant strains as susceptible by phenotypic drug susceptibility testing. This observation potentially allows individualized regimens with high-dose treatment as a potential option for patients with few treatment options. We also found a potentially new bedaquiline resistance mutation in rv0678.

First insights into the phylogenetic diversity of Mycobacterium tuberculosis in Kuwait and evaluation of REBA MTB-MDR assay for rapid detection of MDR-TB

Early detection of Mycobacterium tuberculosis (Mtb) in clinical specimens, its susceptibility to anti-TB drugs and disruption of infection transmission to new hosts are essential components for global tuberculosis (TB) control efforts. This study investigated major Mtb genotypes circulating in Kuwait and evaluated the performance of REBA MTB-MDR (REBA) test in comparison to GenoType MTBDRplus (gMTBDR⁺) assay for rapid detection of resistance of Mtb to isoniazid and rifampicin (MDR-TB). M. tuberculosis isolates (n = 256) originating predominantly from expatriate patients during a 6-month period were tested by spoligotyping and a dendrogram was created by UPGMA using MIRU-VNTRplus software. Phenotypic drug susceptibility testing (DST) was performed by MGIT 960 system. Genotypic DST for isoniazid and rifampicin was done by REBA and gMTBDR⁺ assays. Spoligotyping assigned 188 (73.4%) isolates to specific spoligotype international type (SIT) while 68 isolates exhibited orphan patterns. All major M. tuberculosis lineages were detected and EAI, CAS and Beijing families were predominant. Phylogenetic tree showed 131 patterns with 105 isolates exhibiting a unique pattern while 151 isolates clustered in 26 patterns. Fifteen isolates were resistant to one/more drugs. REBA and gMTBDR⁺ detected isoniazid resistance in 11/12 and 10/12 and rifampicin resistance in 4/5 and 4/5 resistant isolates, respectively. The diversity of SIT patterns are highly suggestive of infection of most expatriate patients with unique Mtb strains, likely acquired in their native countries before their arrival in Kuwait. Both, REBA and gMTBDR⁺ assays performed similarly for detection of resistance of Mtb to isoniazid and rifampicin for rapid detection of MDR-TB.

Transcontinental spread and evolution of Mycobacterium tuberculosis W148 European/Russian clade toward extensively drug resistant tuberculosis

Transmission-driven multi-/extensively drug resistant (M/XDR) tuberculosis (TB) is the largest single contributor to human mortality due to antimicrobial resistance. A few major clades of the Mycobacterium tuberculosis complex belonging to lineage 2, responsible for high prevalence of MDR-TB in Eurasia, show outstanding transnational distributions. Here, we determined factors underlying the emergence and epidemic spread of the W148 clade by genome sequencing and Bayesian demogenetic analyses of 720 isolates from 23 countries. We dated a common ancestor around 1963 and identified two successive epidemic expansions in the late 1980s and late 1990s, coinciding with major socio-economic changes in the post-Soviet Era. These population expansions favored accumulation of resistance mutations to up to 11 anti-TB drugs, with MDR evolving toward additional resistances to fluoroquinolones and second-line injectable drugs within 20 years on average. Timescaled haplotypic density analysis revealed that widespread acquisition of compensatory mutations was associated with transmission success of XDR strains. Virtually all W148 strains harbored a hypervirulence-associated ppe38 gene locus, and incipient recurrent emergence of prpR mutation-mediated drug tolerance was detected. The outstanding genetic arsenal of this geographically widespread M/XDR strain clade represents a “perfect storm” that jeopardizes the successful introduction of new anti-M/XDR-TB antibiotic regimens.

An outbreak of the multidrug-resistant Mycobacterium tuberculosis lineage W148 has spread widely across Russia, Central Asia and Europe. Here, the authors use whole genome sequences of ~700 isolates of this lineage collected over ~20 years to analyze its spread, evolution of drug resistance, and impact of compensatory mutations.

High fluoroquinolone resistance proportions among multidrug-resistant tuberculosis driven by dominant L2 Mycobacterium tuberculosis clones in the Mumbai Metropolitan Region

Background

Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains are a serious health problem in India, also contributing to one-fourth of the global MDR tuberculosis (TB) burden. About 36% of the MDR MTBC strains are reported fluoroquinolone (FQ) resistant leading to high pre-extensively drug-resistant (pre-XDR) and XDR-TB (further resistance against bedaquiline and/or linezolid) rates. Still, factors driving the MDR/pre-XDR epidemic in India are not well defined.

Methods

In a retrospective study, we analyzed 1852 consecutive MTBC strains obtained from patients from a tertiary care hospital laboratory in Mumbai by whole genome sequencing (WGS). Univariate and multivariate statistics was used to investigate factors associated with pre-XDR. Core genome multi locus sequence typing, time scaled haplotypic density (THD) method and homoplasy analysis were used to analyze epidemiological success, and positive selection in different strain groups, respectively.

Results

In total, 1016 MTBC strains were MDR, out of which 703 (69.2%) were pre-XDR and 45 (4.4%) were XDR. Cluster rates were high among MDR (57.8%) and pre-XDR/XDR (79%) strains with three dominant L2 (Beijing) strain clusters (Cl 1–3) representing half of the pre-XDR and 40% of the XDR-TB cases. L2 strains were associated with pre-XDR/XDR-TB (P < 0.001) and, particularly Cl 1–3 strains, had high first-line and FQ resistance rates (81.6–90.6%). Epidemic success analysis using THD showed that L2 strains outperformed L1, L3, and L4 strains in short- and long-term time scales. More importantly, L2 MDR and MDR + strains had higher THD success indices than their not-MDR counterparts. Overall, compensatory mutation rates were highest in L2 strains and positive selection was detected in genes of L2 strains associated with drug tolerance (prpB and ppsA) and virulence (Rv2828c). Compensatory mutations in L2 strains were associated with a threefold increase of THD indices, suggesting improved transmissibility.

Conclusions

Our data indicate a drastic increase of FQ resistance, as well as emerging bedaquiline resistance which endangers the success of newly endorsed MDR-TB treatment regimens. Rapid changes in treatment and control strategies are required to contain transmission of highly successful pre-XDR L2 strains in the Mumbai Metropolitan region but presumably also India-wide.

Tuberculosis in Poland: Epidemiological and Molecular Analysis during the COVID-19 Pandemic

The COVID-19 pandemic may have a negative impact on the proper implementation of TB control programmes and may increase TB incidence rates in the near future. The aim of this study was to perform an epidemiological and molecular analysis of Mycobacterium tuberculosis strains cultured from tuberculosis patients in Poland in 2020 and to compare the results of monitoring drug-resistant tuberculosis in Poland with previous studies in 2012 and 2016. The analysis was based on questionnaires and strains sent by regional laboratories during the 12 months of 2020. Molecular analysis was performed by spoligotyping 20% of the strains sensitive to the four primary antimycobacterial drugs and all of the drug-resistant strains. The number of strains sent for analysis dropped threefold, from 4136 in 2012 to 1383 in 2020. The incidence of tuberculosis among men was higher than among women. There was an increase in strains’ resistance to antimycobacterial drugs in both newly diagnosed patients, from 4.4% in 2012 to 6.1% in 2020, and previously treated patients, from 11.7% to 12.3%. Four-year resistance increased to 1% and 2.1%, respectively. The spoligotype SIT1 was the most abundant among the resistant strains (17%), and SIT53 (13.9%) was the most common among susceptible strains.

Diversity and novel mutations in membrane transporters of Mycobacterium tuberculosis

Mycobacterium tuberculosis (MTB), the causative agent of tuberculosis (TB), encodes a family of membrane proteins belonging to Resistance-Nodulation-Cell Division (RND) permeases also called multidrug resistance pumps. Mycobacterial membrane protein Large (MmpL) transporters represent a subclass of RND transporters known to participate in exporting of lipid components across the cell envelope. These proteins perform an essential role in MTB survival; however, there are no data regarding mutations in MmpL, polyketide synthase (PKS) and acyl-CoA dehydrogenase FadE proteins from Khyber Pakhtunkhwa, Pakistan. This study aimed to screen mutations in transmembrane transporter proteins including MmpL, PKS and Fad through whole-genome sequencing (WGS) in local isolates of Khyber Pakhtunkhwa province, Pakistan. Fourteen samples were collected from TB patients and drug susceptibility testing was performed. However, only three samples were completely sequenced. Moreover, 209 whole-genome sequences of the same geography were also retrieved from NCBI GenBank to analyze the diversity of mutations in MmpL, PKS and Fad proteins. Among the 212 WGS (Accession ID: PRJNA629298, PRJNA629388, and ERR2510337-ERR2510345, ERR2510546-ERR2510645), numerous mutations in Fad (n = 756), PKS (n = 479), and MmpL (n = 306) have been detected. Some novel mutations were also detected in MmpL, PKS and acyl-CoA dehydrogenase Fad. Novel mutations including Asn576Ser in MmpL8, Val943Gly in MmpL9 and Asn145Asp have been detected in MmpL3. The presence of a large number of mutations in the MTB membrane may have functional consequences on proteins. However, further experimental studies are needed to elucidate the variants' effect on MmpL, PKS and FadE functions.

A convolutional neural network highlights mutations relevant to antimicrobial resistance in Mycobacterium tuberculosis

Long diagnostic wait times hinder international efforts to address antibiotic resistance in M. tuberculosis. Pathogen whole genome sequencing, coupled with statistical and machine learning models, offers a promising solution. However, generalizability and clinical adoption have been limited by a lack of interpretability, especially in deep learning methods. Here, we present two deep convolutional neural networks that predict antibiotic resistance phenotypes of M. tuberculosis isolates: a multi-drug CNN (MD-CNN), that predicts resistance to 13 antibiotics based on 18 genomic loci, with AUCs 82.6-99.5% and higher sensitivity than state-of-the-art methods; and a set of 13 single-drug CNNs (SD-CNN) with AUCs 80.1-97.1% and higher specificity than the previous state-of-the-art. Using saliency methods to evaluate the contribution of input sequence features to the SD-CNN predictions, we identify 18 sites in the genome not previously associated with resistance. The CNN models permit functional variant discovery, biologically meaningful interpretation, and clinical applicability.

Mycobacterium tuberculosis Lineages Associated with Mutations and Drug Resistance in Isolates from India

Current knowledge on resistance-conferring determinants in Mycobacterium tuberculosis is biased toward globally dominant lineages 2 and 4. In contrast, lineages 1 and 3 are predominant in India. In this study, we performed whole-genome sequencing of 498 MDR M. tuberculosis isolates from India to determine the prevalence of drug resistance mutations and to understand the genomic diversity. A retrospective collection of 498 M. tuberculosis isolates submitted to the National Institute for Research in Tuberculosis for phenotypic susceptibility testing between 2014 to 2016 were sequenced. Genotypic resistance prediction was performed using known resistance-conferring determinants. Genotypic and phenotypic results for 12 antituberculosis drugs were compared, and sequence data were explored to characterize lineages and their association with drug resistance. Four lineages were identified although lineage 1 predominated (43%). The sensitivity of prediction for isoniazid and rifampicin was 92% and 98%, respectively. We observed lineage-specific variations in the proportion of isolates with resistance-conferring mutations, with drug resistance more common in lineages 2 and 3. Disputed mutations (codons 430, 435, 445, and 452) in the rpoB gene were more common in isolates other than lineage 2. Phylogenetic analysis and pairwise SNP difference revealed high genetic relatedness of lineage 2 isolates. WGS based resistance prediction has huge potential, but knowledge of regional and national diversity is essential to achieve high accuracy for resistance prediction. IMPORTANCE Current knowledge on resistance-conferring determinants in Mycobacterium tuberculosis is biased toward globally dominant lineages 2 and 4. In contrast, lineages 1 and 3 are predominant in India. We performed whole-genome sequencing of 498 MDR M. tuberculosis isolates from India to determine the prevalence of drug resistance mutations and to understand genomic diversity. Four lineages were identified although lineage 1 predominated (43%). The sensitivity of prediction for isoniazid and rifampicin was 92% and 98%, respectively. We observed lineage-specific variations in the proportion of isolates with resistance-conferring mutations, with drug resistance more common in lineages 2 and 3. Disputed mutations (codons 430, 435, 445, and 452) in the rpoB gene were more common in isolates other than lineage 2. Phylogenetic analysis and pairwise SNP difference revealed high genetic relatedness of lineage 2 isolates. WGS based resistance prediction has huge potential, but knowledge of regional and national diversity is essential to achieve high accuracy for resistance prediction.

Origin and Global Expansion of Mycobacterium tuberculosis Complex Lineage 3

Simple Summary

Tuberculosis still causes 1.5 million deaths annually and is mainly caused by Mycobacterium tuberculosis complex strains belonging to three evolutionary modern lineages (Lineages 2–4). While Lineage 2 and Lineage 4 virtually conquered the world, Lineage 3 is particularly successful in Northern and Eastern Africa, as well as in Southern Asia, the suspected evolutionary origin of these strains. Here, we sought to understand how Lineage 3 strains came to the African continent. To this end, we performed routine genotyping to characterize over 2500 clinical isolates from 38 countries. We then selected a representative collection of 373 isolates for a whole-genome analysis and a modeling approach to infer the geographic origin of different sublineages. In fact, the origin of Lineage 3 could be located in India, and we found evidence for independent introductions of four distinct sublineages into North/East Africa, in line with known ancient exchanges and migrations between both world regions. Our study illustrates that the evolutionary history of humans and their pathogens are closely connected and further provides a systematic understanding of the genomic diversity of Lineage 3, which could be important for the development of new tuberculosis vaccines or new therapeutics.

Abstract

Mycobacterium tuberculosis complex (MTBC) Lineage 3 (L3) strains are abundant in world regions with the highest tuberculosis burden. To investigate the population structure and the global diversity of this major lineage, we analyzed a dataset comprising 2682 L3 strains from 38 countries over 5 continents, by employing 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats genotyping (MIRU-VNTR) and drug susceptibility testing. We further combined whole-genome sequencing (WGS) and phylogeographic analysis for 373 strains representing the global L3 genetic diversity. Ancestral state reconstruction confirmed that the origin of L3 strains is located in Southern Asia and further revealed multiple independent introduction events into North-East and East Africa. This study provides a systematic understanding of the global diversity of L3 strains and reports phylogenetic variations that could inform clinical trials which evaluate the effectivity of new drugs/regimens or vaccine candidates.

Prediction of drug resistance profile of multidrug-resistant Mycobacterium tuberculosis (MDR-MTB) isolates from newly diagnosed case by whole genome sequencing (WGS): a study from a high tuberculosis burden country

OBJECTIVES

Our aim was to assess the ability of the Whole-genome sequencing (WGS) in predicting drug resistance profile of multidrug-resistant mycobacterium tuberculosis (MDR-MTB) from newly diagnosed cases in China.

METHODS

We validated the Phenotypic drug Sensitivity Test (pDST) for 12 anti-tuberculosis drugs using the Bactec MGIT 960 system. We described the characteristics of the isolates enrolled and compared the pDST results with resistance profiles predicted by WGS.

RESULTS

The pDST showed that of the 43 isolates enrolled, 25.6% were sensitive to rifabutin (RFB); 97.7%、97.7%、93.0% and 93.0% were sensitive to cycloserine (Cs), amikacin/kanamycin (Ak/Km), para-aminosalicylic acid (Pas) and ethionamide Eto), respectively; 18.6% were resistant to fluoroquinolones (FQs) or second-line injections. Genotype DST determined by WGS of Ak/Km、Eto and RFP reached high consistency to 97.7% compared with pDST, followed by moxifloxacin (Mfx) 95.3%, levofloxaci (Lfx) and Pas 93%, streptomycin (Sm) 90.3%. The genotype DST of RFB and EMB showed low consistency with the pDST of 67.2 and 79.1%. WGS also detected 27.9% isolates of pyrazinamide(PZA)-related drug-resistant mutation. No mutations associated with linezolid (Lzd), bedaquiline (Bdq) and clofazimine (Cfz) were detectd.

CONCLUSIONS

WGS has the potential to infer resistance profiles without time-consuming phenotypic methods, which could be provide a basis to formulate reasonable treatment in high TB burden areas.

Molecular characteristics of Mycobacterium tuberculosis drug-resistant isolates from HIV- and HIV+ tuberculosis patients in Russia

Background

High burden of drug-resistant (DR) tuberculosis (TB) is a significant threat to national TB control programs all over the world and in the Russian Federation. Different Mycobacterium tuberculosis (MTB) genotypes are hypothesized to have specific characteristics affecting TB control programs. For example, Beijing strains are supposed to have higher mutation rates compared to strains of other genotypes and subsequently higher capability to develop drug-resistance.

Results

Clinical MTB isolates from HIV- and HIV+ patients from four regions of Russia were analyzed for genotypes and mutations conferring resistance to Isoniazid, Rifampicin, Ethambutol, aminoglycosides, and fluoroquinolones. Analysis of genotypes and polymorphism of genomic loci according to the HIV status of the patients – sources of MTB isolates were performed. Studied MTB isolates from HIV- TB patients belonged to 15 genotypes and from HIV + TB patients – to 6 genotypes. Beijing clinical isolates dominated in HIV- (64,7%) and HIV+ (74,4%) groups. Other isolates were of LAM (including LAM1 and LAM9), Ural, and 4 minor groups of genotypes (including 5 subclones T). The spectrum of genotypes in the HIV- group was broader than in the HIV+ group. PR of B0/W148 Beijing was significantly lower than of other Beijing genotypes in susceptible and MDR-XDR isolates. Rates of isolates belonging to non-Beijing genotypes were higher than Beijing in susceptible isolates from HIV- patients.

Conclusions

Beijing genotype isolates prevailed in clinical isolates of all drug susceptibility profiles both from HIV- and HIV+ patients, although B0/W148 Beijing genotype did not dominate in this study. Genome loci and mutations polymorphisms were more pronounced in clinical isolates from HIV- patients, than from HIV+.

Evolution and spread of a highly drug resistant strain of Mycobacterium tuberculosis in Papua New Guinea

BACKGROUND

Molecular mechanisms determining the transmission and prevalence of drug resistant tuberculosis (DR-TB) in Papua New Guinea (PNG) are poorly understood. We used genomic and drug susceptibility data to explore the evolutionary history, temporal acquisition of resistance and transmission dynamics of DR-TB across PNG.

METHODS

We performed whole genome sequencing on isolates from Central Public Health Laboratory, PNG, collected 2017-2019. Data analysis was done on a composite dataset that also included 100 genomes previously sequenced from Daru, PNG (2012-2015).

RESULTS

Sampled isolates represented 14 of the 22 PNG provinces, the majority (66/94; 70%) came from the National Capital District (NCD). In the composite dataset, 91% of strains were Beijing 2.2.1.1, identified in 13 provinces. Phylogenetic tree of Beijing strains revealed two clades, Daru dominant clade (A) and NCD dominant clade (B). Multi-drug resistance (MDR) was repeatedly and independently acquired, with the first MDR cases in both clades noted to have emerged in the early 1990s, while fluoroquinolone resistance emerged in 2009 (95% highest posterior density 2000-2016). We identified the presence of a frameshift mutation within Rv0678 (p.Asp47fs) which has been suggested to confer resistance to bedaquiline, despite no known exposure to the drug. Overall genomic clustering was significantly associated with rpoC compensatory and inhA promoter mutations (p < 0.001), with high percentage of most genomic clusters (12/14) identified in NCD, reflecting its role as a potential national amplifier.

CONCLUSIONS

The acquisition and evolution of drug resistance among the major clades of Beijing strain threaten the success of DR-TB treatment in PNG. With continued transmission of this strain in PNG, genotypic drug resistance surveillance using whole genome sequencing is essential for improved public health response to outbreaks. With occurrence of resistance to newer drugs such as bedaquiline, knowledge of full drug resistance profiles will be important for optimal treatment selection.

Whole genome sequencing of multidrug-resistant Mycobacterium tuberculosis isolates collected in the Czech Republic, 2005-2020

The emergence and spread of resistant tuberculosis (TB) pose a threat to public health, so it is necessary to diagnose the drug-resistant forms in a clinically short time frame and closely monitor their transmission. In this study, we carried out a first whole genome sequencing (WGS)-based analysis of multidrug resistant (MDR) M. tuberculosis strains to explore the phylogenetic lineages diversity, drug resistance mechanisms, and ongoing transmission chains within the country. In total, 65 isolates phenotypically resistant to at least rifampicin and isoniazid collected in the Czech Republic in 2005-2020 were enrolled for further analysis. The agreement of the results obtained by WGS with phenotypic drug susceptibility testing (pDST) in the determination of resistance to isoniazid, rifampicin, pyrazinamide, streptomycin, second-line injectables and fluoroquinolones was more than 80%. Phylogenetic analysis of WGS data revealed that the majority of MDR M. tuberculosis isolates were the Beijing lineage 2.2.1 (n = 46/65; 70.8%), while the remaining strains belonged to Euro-American lineage. Cluster analysis with a predefined cut-off distance of less than 12 single nucleotide polymorphisms between isolates showed 19 isolates in 6 clusters (clustering rate 29.2%), located mainly in the region of the capital city of Prague. This study highlights the utility of WGS as a high-resolution approach in the diagnosis, characterization of resistance patterns, and molecular-epidemiological analysis of resistant TB in the country.

Multiplexed Strain Phenotyping Defines Consequences of Genetic Diversity in Mycobacterium tuberculosis for Infection and Vaccination Outcomes

There is growing evidence that genetic diversity in Mycobacterium tuberculosis, the causative agent of tuberculosis, contributes to the outcomes of infection and public health interventions, such as vaccination. Epidemiological studies suggest that among the phylogeographic lineages of M. tuberculosis, strains belonging to a sublineage of Lineage 2 (mL2) are associated with concerning clinical features, including hypervirulence, treatment failure, and vaccine escape. The global expansion and increasing prevalence of this sublineage has been attributed to the selective advantage conferred by these characteristics, yet confounding host and environmental factors make it difficult to identify the bacterial determinants driving these associations in human studies. Here, we developed a molecular barcoding strategy to facilitate high-throughput, experimental phenotyping of M. tuberculosis clinical isolates. This approach allowed us to characterize growth dynamics for a panel of genetically diverse M. tuberculosis strains during infection and after vaccination in the mouse model. We found that mL2 strains exhibit distinct growth dynamics in vivo and are resistant to the immune protection conferred by Bacillus Calmette-Guerin (BCG) vaccination. The latter finding corroborates epidemiological observations and demonstrates that mycobacterial features contribute to vaccine efficacy. To investigate the genetic and biological basis of mL2 strains’ distinctive phenotypes, we performed variant analysis, transcriptional studies, and genome-wide transposon sequencing. We identified functional genetic changes across multiple stress and host response pathways in a representative mL2 strain that are associated with variants in regulatory genes. These adaptive changes may underlie the distinct clinical characteristics and epidemiological success of this lineage.

IMPORTANCE Tuberculosis, caused by the bacterium Mycobacterium tuberculosis, is a remarkably heterogeneous disease, a feature that complicates clinical care and public health interventions. The contributions of pathogen genetic diversity to this heterogeneity are uncertain, in part due to the challenges of experimentally manipulating M. tuberculosis, a slow-growing, biosafety level 3 organism. To overcome these challenges, we applied a molecular barcoding strategy to a panel of M. tuberculosis clinical isolates. This novel application of barcoding permitted the high-throughput characterization of M. tuberculosis strain growth dynamics and vaccine resistance in the mouse model of infection. Integrating these results with genomic analyses, we uncover bacterial pathways that contribute to infection outcomes, suggesting targets for improved therapeutics and vaccines.

The TbD1 Locus Mediates a Hypoxia-Induced Copper Response in Mycobacterium bovis

The Mycobacterium tuberculosis complex (MTBC) contains the causative agents of tuberculosis (TB) in mammals. The archetypal members of the MTBC, Mycobacterium tuberculosis and Mycobacterium bovis, cause human tuberculosis and bovine tuberculosis, respectively. Although M. tuberculosis and M. bovis share over 99.9% genome identity, they show distinct host adaptation for humans and animals; hence, while the molecular basis of host adaptation is encoded in their genomes, the mechanistic basis of host tropism is still unclear. Exploration of the in vitro phenotypic consequences of known genetic difference between M. bovis and M. tuberculosis offers one route to explore genotype–phenotype links that may play a role in host adaptation. The TbD1 (“Mycobacterium tuberculosis deletion 1 region”) locus encompasses the mmpS6 and mmpL6 genes. TbD1 is absent in M. tuberculosis “modern” lineages (Lineages 2, 3, and 4) but present in “ancestral” M. tuberculosis (Lineages 1 and 7), Mycobacterium africanum lineages (Lineages 5 and 6), newly identified M. tuberculosis lineages (Lineages 8 and 9), and animal adapted strains, such as M. bovis. The function of TbD1 has previously been investigated in M. tuberculosis, where conflicting data has emerged on the role of TbD1 in sensitivity to oxidative stress, while the underlying mechanistic basis of such a phenotype is unclear. In this study, we aimed to shed further light on the role of the TbD1 locus by exploring its function in M. bovis. Toward this, we constructed an M. bovis TbD1 knockout (ΔTbD1) strain and conducted comparative transcriptomics to define global gene expression profiles of M. bovis wild-type (WT) and the ΔTbD1 strains under in vitro culture conditions (rolling and standing cultures). This analysis revealed differential induction of a hypoxia-driven copper response in WT and ΔTbD1 strains. In vitro phenotypic assays demonstrated that the deletion of TbD1 sensitized M. bovis to H2O2 and hypoxia-specific copper toxicity. Our study provides new information on the function of the TbD1 locus in M. bovis and its role in stress responses in the MTBC.

The aldehyde hypothesis: metabolic intermediates as antimicrobial effectors

There are many reactive intermediates found in metabolic pathways. Could these potentially toxic molecules be exploited for an organism's benefit? We propose that during certain microbial infections, the production of inherently reactive aldehydes by an infected host is a previously unappreciated innate immune defence mechanism. While there has been a significant focus on the effects of aldehydes on mammalian physiology, the idea that they might be exploited or purposefully induced to kill pathogens is new. Given that aldehydes are made as parts of metabolic programmes that accompany immune cell activation by the cytokine interferon-gamma (IFN-γ) during infections, we hypothesize that aldehydes are among the arsenal of IFN-γ-inducible effectors needed for pathogen control.

Distribution and identification of Mycobacterium tuberculosis lineage in Kashgar prefecture

Objectives

Kashgar prefecture is an important transportation and trade hub with a high incidence of tuberculosis. The following study analyzed the composition and differences in Mycobacterium tuberculosis (M.tb) lineage and specific tags to distinguish the lineage of the M.tb in Kashgar prefecture, thus providing a basis for the classification and diagnosis of tuberculosis in this area.

Methods

Whole-genome sequencing (WGS) of 161 M.tb clinical strains was performed. The phylogenetic tree was constructed using Maximum Likelihood (ML) based on single nucleotide polymorphisms (SNPs) and verified through principal component analysis (PCA). The composition structure of M.tb in different regions was analyzed by combining geographic information.

Results

M.tb clinical strains were composed of lineage 2 (73/161, 45.34%), lineage 3 (52/161, 32.30%) and lineage 4 (36/161, 22.36%). Moreover, the 3 lineages were subdivided into 11 sublineages, among which lineage 2 included lineage 2.2.2/Asia Ancestral 1 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 2 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 3 (18/73, 24.66%), and lineage 2.2.1-Modern Beijing (39/73, 53.42%). Lineage 3 included lineage 3.2 (14/52, 26.92%) and lineage 3.3 (38/52, 73.08%), while lineage 4 included lineage 4.1 (3/36, 8.33%), lineage 4.2 (2/36, 5.66%), lineage 4.4.2 (1/36, 2.78%), lineage 4.5 (28/36, 77.78%) and lineage 4.8 (2/36, 5.66%), all of which were consistent with the PCA results. One hundred thirty-six markers were proposed for discriminating known circulating strains. Reconstruction of a phylogenetic tree using the 136 SNPs resulted in a tree with the same number of delineated clades. Based on geographical location analysis, the composition of Lineage 2 in Kashgar prefecture (45.34%) was lower compared to other regions in China (54.35%-90.27%), while the composition of Lineage 3 (32.30%) was much higher than in other regions of China (0.92%-2.01%), but lower compared to the bordering Pakistan (70.40%).

Conclusion

Three lineages were identified in M.tb clinical strains from Kashgar prefecture, with 136 branch-specific SNP. Kashgar borders with countries that have a high incidence of tuberculosis, such as Pakistan and India, which results in a large difference between the M.tb lineage and sublineage distribution in this region and other provinces of China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07307-4.

Transborder molecular analysis of drug-resistant tuberculosis in Mongolia and Eastern Siberia, Russia

Eastern Siberia (Russia) and Mongolia are borderline regions in Asia with a high incidence of tuberculosis (TB). In this study, we investigated the transborder transmission of Mycobacterium tuberculosis with a focus on endemic and epidemic clones and drug resistance. M. tuberculosis isolates (287 from Mongolia and 754 from Russia) were collected using cross-sectional population-based surveys between 2010 and 2016. The isolates were genotyped using 24 variable number of tandem repeat (VNTR) loci and by testing of the key markers to discriminate within the Beijing genotype. All isolates were divided into 427 MIRU-types that were assigned to Lineage 2 (Beijing) and Lineage 4 (Ural, Haarlem, Latin American-Mediterranean [LAM], S, unclassified). The Beijing genotype was dominant in both countries (69% in Russia, 75% in Mongolia). However, the Beijing isolates differed significantly between the countries, in terms of the identified subtypes. LAM was the most common non-Beijing genotype (11.1% in Mongolia and 14.9% in Russia) and LAM isolates mostly belonged to the LAM-RUS branch in both countries. The MDR rate was higher in Russia than in Mongolia among newly diagnosed patients: 29.4% versus 5.6% (p < 0.001). In Mongolia, the MDR rate was similar in Beijing (29.7%) and non-Beijing (27.5%) genotypes. In Russia, a higher MDR rate was observed in (i) Beijing compared to non-Beijing (48.7% versus 38.3%, p = 0.03) and (ii) Beijing B0/W148 compared to Beijing Central Asian/Russian (63.4% versus 37.3%, p<0.001). In conclusion, the M. tuberculosis population structure in Mongolia was shaped by mainly historical interaction with China (dominance of the Beijing genotype) and Northern Eurasia (presence of the LAM-RUS branch). In contrast, the transborder transmission of M. tuberculosis since the 1990s between Mongolia and its neighbors has been negligible, and the adverse trends of MDR-TB in Russia did not impact the current situation in Mongolia. This article is protected by copyright. All rights reserved.

Molecular Epidemiology of Mycobacterium tuberculosis Complex Strains in Urban and Slum Settings of Nairobi, Kenya

Kenya is a country with a high tuberculosis (TB) burden. However, knowledge on the genetic diversity of Mycobacterium tuberculosis complex (MTBC) strains and their transmission dynamics is sparsely available. Hence, we used whole-genome sequencing (WGS) to depict the genetic diversity, molecular markers of drug resistance, and possible transmission clusters among MTBC strains in urban and slum settings of Nairobi. We analyzed 385 clinical MTBC isolates collected between 2010 and 2015 in combination with patients’ demographics. We showed that the MTBC population mainly comprises strains of four lineages (L1–L4). The two dominating lineages were L4 with 55.8% (n = 215) and L3 with 25.7% (n = 99) of all strains, respectively. Genome-based cluster analysis showed that 30.4% (117/385) of the strains were clustered using a ≤5 single-nucleotide polymorphism (SNP) threshold as a surrogate marker for direct patient-to-patient MTBC transmission. Moreover, 5.2% (20/385) of the strains were multidrug-resistant (MDR), and 50.0% (n = 10) were part of a genome-based cluster (i.e., direct MDR MTBC transmission). Notably, 30.0% (6/20) of the MDR strains were resistant to all first-line drugs and are part of one molecular cluster. Moreover, TB patients in urban living setting had 3.8 times the odds of being infected with a drug-resistant strain as compared to patients from slums (p-value = 0.002). Our results show that L4 strains are the main causative agent of TB in Nairobi and MDR strain transmission is an emerging concern in urban settings. This emphasizes the need for more focused infection control measures and contact tracing of patients with MDR TB to break the transmission chains.

A Cluster-based Method to Quantify Individual Heterogeneity in Tuberculosis Transmission

BACKGROUND

Recent evidence suggests transmission of Mycobacterium tuberculosis (Mtb) may be characterized by extreme individual heterogeneity in secondary cases (i.e., few cases account for the majority of transmission). Such heterogeneity implies outbreaks are rarer but more extensive and has profound implications in infectious disease control. However, discrete person-to-person transmission events in tuberculosis (TB) are often unobserved, precluding our ability to directly quantify individual heterogeneity in TB epidemiology.

METHODS

We used a modified negative binomial branching process model to quantify the extent of individual heterogeneity using only observed transmission cluster size distribution data (i.e., the simple sum of all cases in a transmission chain) without knowledge of individual-level transmission events. The negative binomial parameter k quantifies the extent of individual heterogeneity (generally, indicates extensive heterogeneity, and as transmission becomes more homogenous). We validated the robustness of the inference procedure considering common limitations affecting cluster size data. Finally, we demonstrate the epidemiologic utility of this method by applying it to aggregate US molecular surveillance data from the US Centers for Disease Control and Prevention.

RESULTS

The cluster-based method reliably inferred k using TB transmission cluster data despite a high degree of bias introduced into the model. We found that the TB transmission in the United States was characterized by a high propensity for extensive outbreaks (; 95% confidence interval = 0.09, 0.10).

CONCLUSIONS

The proposed method can accurately quantify critical parameters that govern TB transmission using simple, more easily obtainable cluster data to improve our understanding of TB epidemiology.

Application and Challenge of 3rd Generation Sequencing for Clinical Bacterial Studies

Over the past 25 years, the powerful combination of genome sequencing and bioinformatics analysis has played a crucial role in interpreting information encoded in bacterial genomes. High-throughput sequencing technologies have paved the way towards understanding an increasingly wide range of biological questions. This revolution has enabled advances in areas ranging from genome composition to how proteins interact with nucleic acids. This has created unprecedented opportunities through the integration of genomic data into clinics for the diagnosis of genetic traits associated with disease. Since then, these technologies have continued to evolve, and recently, long-read sequencing has overcome previous limitations in terms of accuracy, thus expanding its applications in genomics, transcriptomics and metagenomics. In this review, we describe a brief history of the bacterial genome sequencing revolution and its application in public health and molecular epidemiology. We present a chronology that encompasses the various technological developments: whole-genome shotgun sequencing, high-throughput sequencing, long-read sequencing. We mainly discuss the application of next-generation sequencing to decipher bacterial genomes. Secondly, we highlight how long-read sequencing technologies go beyond the limitations of traditional short-read sequencing. We intend to provide a description of the guiding principles of the 3rd generation sequencing applications and ongoing improvements in the field of microbial medical research.

Connection between two historical tuberculosis outbreak sites in Japan, Honshu, by a new ancestral Mycobacterium tuberculosis L2 sublineage

By gathering 680 publicly available Sequence Read Archives from isolates of Mycobacterium tuberculosis complex (MTBC) including 190 belonging to the lineage 2 Beijing , and using an in-house bioinformatical pipeline, the TB-Annotator , that analyses more than 50 000 characters, we describe herein a new L2 sublineage from 20 isolates found in the Tochigi province, (Japan), that we designate as asia ancestral 5 (AAnc5). These isolates harbour a number of specific criteria (42 SNPs) and their intra-cluster pairwise distance suggests historical and not epidemiological transmission. These isolates harbour a mutation in rpoC , and do not fulfil, any of the modern Beijing lineage criteria, nor any of the other ancestral Beijing lineages described so far. Asia ancestral 5 isolates do not possess mutT2 58 and ogt 12 characteristics of modern Beijing , but possess ancestral Beijing SNPs characteristics. By looking into the literature, we found a reference isolate ID381, described in Kobe and Osaka belonging to the ‘G3’ group, sharing 36 out of the 42 specific SNPs found in AAnc5. We also assessed the intermediate position of the asia ancestral 4 (AAnc4) sublineage recently described in Thailand and propose an improved classification of the L2 that now includes AAnc4 and AAnc5. By increasing the recruitment into TB-Annotator to around 3000 genomes (including 642 belonging to L2), we confirmed our results and discovered additional historical ancestral L2 branches that remain to be investigated in more detail. We also present, in addition, some anthropological and historical data from Chinese and Japan history of tuberculosis, as well as from Korea, that could support our results on L2 evolution. This study shows that the reconstruction of the early history of tuberculosis in Asia is likely to reveal complex patterns since its emergence.

OUP accepted manuscript

Background

Rifampicin- or multidrug-resistant (RR/MDR) Mycobacterium tuberculosis complex (MTBC) strains account for considerable morbidity and mortality globally. WGS-based prediction of drug resistance may guide clinical decisions, especially for the design of RR/MDR-TB therapies.

Methods

We compared WGS-based drug resistance-predictive mutations for 42 MTBC isolates from MDR-TB patients in Tanzania with the MICs of 14 antibiotics measured in the Sensititre™ MycoTB assay. An isolate was phenotypically categorized as resistant if it had an MIC above the epidemiological-cut-off (ECOFF) value, or as susceptible if it had an MIC below or equal to the ECOFF.

Results

Overall, genotypically non-wild-type MTBC isolates with high-level resistance mutations (gNWT-R) correlated with isolates with MIC values above the ECOFF. For instance, the median MIC value (mg/L) for rifampicin-gNWT-R strains was >4.0 (IQR 4.0–4.0) compared with 0.5 (IQR 0.38–0.50) in genotypically wild-type (gWT-S, P < 0.001); isoniazid-gNWT-R >4.0 (IQR 2.0–4.0) compared with 0.25 (IQR 0.12–1.00) among gWT-S (P = 0.001); ethionamide-gNWT-R 15.0 (IQR 10.0–20.0) compared with 2.50 (IQR; 2.50–5.00) among gWT-S (P < 0.001). WGS correctly predicted resistance in 95% (36/38) and 100% (38/38) of the rifampicin-resistant isolates with ECOFFs >0.5 and >0.125 mg/L, respectively. No known resistance-conferring mutations were present in genes associated with resistance to fluoroquinolones, aminoglycosides, capreomycin, bedaquiline, delamanid, linezolid, clofazimine, cycloserine, or p-amino salicylic acid.

Conclusions

WGS-based drug resistance prediction worked well to rule-in phenotypic drug resistance and the absence of second-line drug resistance-mediating mutations has the potential to guide the design of RR/MDR-TB regimens in the future.

Putative extensive and pre-extensive drug resistant-tuberculosis associated with unusual genotypes on the Thailand-Myanmar border

Extensive drug-resistant tuberculosis (XDR-TB) is highly life threatening and its diagnosis is usually difficult and time-consuming. Here we present the first two cases of XDR and pre-XDR-TB diagnosed in 2018 on the Thailand-Myanmar border, more specifically in Tak province. Rapid detection of XDR-TB was performed by loop-mediated isothermal amplification (LAMP), Xpert MTB/RIF, and line probe assays. Mutation analyses targeting rpoB, katG, inhA, gyrA and rrs genes showed an association with drug-resistant phenotypes, except for rifampicin resistance. Spoligotyping revealed uncommon Beijing and T2 genotypes and the analysis of M. tuberculosis interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) showed the presence of more polymorphisms. This report highlights the importance of the early detection of drug-resistant tuberculosis by molecular tests followed by phenotyping assays. Based on the up-to-date definition of XDR- and pre-XDR-TB, the susceptibility testing for bedaquiline and linezolid is required and the two reported cases may correspond to putative XDR-TB.

Genetic diversity and drug susceptibility of Mycobacterium tuberculosis in a city with a high prevalence of drug resistant tuberculosis from Southeast of Mexico

Background

Mexico is on the top five countries with the highest number of TB cases in America continent, nevertheless, information about genotypes circulating is practically unknown. Considering the above this study aims to characterize the genetic diversity of TB in the city of Veracruz, México.

Methods

A cross-sectional study was conducted among positive smear samples from patients living in Veracruz City, samples were cultured, and first-line drug profiles determined. Genotyping was made by spoligotyping and MIRU-VNTR 24 loci. Associations of lineages, clusters, and variables were also analyzed.

Results

Among the 202 isolates analyzed resistance to at least one drug was observed in 60 (30%) isolates and 41(20%) were multidrug-resistant. Three major lineages were identified: L4/Euro-American (88%), L1/Indo-Oceanic (9%), and L2/East Asian (3%). The Euro-American lineage included more than six sublineages, the most abundant were: H (32%), T (23%), LAM (18%), and X (12%). 140 isolates (70%) were placed in 42 SITs patterns.

Conclusions

These results provide the first baseline data on the genetic structure of TB in the city of Veracruz. Sublineages H, X and LAM were predominant; however, it was founded an important diversity of genotypes that could contribute to the dispersion of TB and explain the high prevalence. This information might be useful for the development of further interventions to reduce impact of TB.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-021-06904-z.

Revised nomenclature and SNP barcode for Mycobacterium tuberculosis lineage 2

Mycobacterium tuberculosis (Mtb) lineage 2 (L2) strains are present globally, contributing to a widespread tuberculosis (TB) burden, particularly in Asia where both prevalence of TB and numbers of drug resistant TB are highest. The increasing availability of whole-genome sequencing (WGS) data worldwide provides an opportunity to improve our understanding of the global genetic diversity of Mtb L2 and its association with the disease epidemiology and pathogenesis. However, existing L2 sublineage classification schemes leave >20 % of the Modern Beijing isolates unclassified. Here, we present a revised SNP-based classification scheme of L2 in a genomic framework based on phylogenetic analysis of >4000 L2 isolates from 34 countries in Asia, Eastern Europe, Oceania and Africa. Our scheme consists of over 30 genotypes, many of which have not been described before. In particular, we propose six main genotypes of Modern Beijing strains, denoted L2.2.M1–L2.2.M6. We also provide SNP markers for genotyping L2 strains from WGS data. This fine-scale genotyping scheme, which can classify >98 % of the studied isolates, serves as a basis for more effective monitoring and reporting of transmission and outbreaks, as well as improving genotype-phenotype associations such as disease severity and drug resistance. This article contains data hosted by Microreact.

Practical approach to detection and surveillance of emerging highly resistant Mycobacterium tuberculosis Beijing 1071-32-cluster

Ancient sublineage of the Mycobacterium tuberculosis Beijing genotype is endemic and prevalent in East Asia and rare in other world regions. While these strains are mainly drug susceptible, we recently identified a novel clonal group Beijing 1071-32 within this sublineage emerging in Siberia, Russia and present in other Russian regions. This cluster included only multi/extensive drug resistant (MDR/XDR) isolates. Based on the phylogenetic analysis of the available WGS data, we identified three synonymous SNPs in the genes Rv0144, Rv0373c, and Rv0334 that were specific for the Beijing 1071-32-cluster and developed a real-time PCR assay for their detection. Analysis of the 2375 genetically diverse M. tuberculosis isolates collected between 1996 and 2020 in different locations (European and Asian parts of Russia, former Soviet Union countries, Albania, Greece, China, Vietnam, Japan and Brazil), confirmed 100% specificity and sensitivity of this real-time PCR assay. Moreover, the epidemiological importance of this strain and the newly developed screening assay is further stressed by the fact that all identified Beijing 1071-32 isolates were found to exhibit MDR genotypic profiles with concomitant resistance to additional first-line drugs due to a characteristic signature of six mutations in rpoB450, rpoC485, katG315, katG335, rpsL43 and embB497. In conclusion, this study provides a set of three concordant SNPs for the detection and screening of Beijing 1071-32 isolates along with a validated real-time PCR assay easily deployable across multiple settings for the epidemiological tracking of this significant MDR cluster.

Genomic Profiling of Mycobacterium tuberculosis Strains, Myanmar

Multidrug resistance is a major threat to global elimination of tuberculosis (TB). We performed phenotypic drug-susceptibility testing and whole-genome sequencing for 309 isolates from 342 consecutive patients who were given a diagnosis of TB in Yangon, Myanmar, during July 2016‒June 2018. We identified isolates by using the GeneXpert platform to evaluate drug-resistance profiles. A total of 191 (62%) of 309 isolates had rifampin resistance; 168 (88%) of these rifampin-resistant isolates were not genomically related, indicating the repeated emergence of resistance in the population, rather than extensive local transmission. We did not detect resistance mutations to new oral drugs, including bedaquiline and pretomanid. The current GeneXpert MTB/RIF system needs to be modified by using the newly launched Xpert MTB/XDR cartridge or line-probe assay. Introducing new oral drugs to replace those currently used in treatment regimens for multidrug-resistant TB will also be useful for treating TB in Myanmar.

Population structure, biogeography and transmissibility of Mycobacterium tuberculosis

Mycobacterium tuberculosis is a clonal pathogen proposed to have co-evolved with its human host for millennia, yet our understanding of its genomic diversity and biogeography remains incomplete. Here we use a combination of phylogenetics and dimensionality reduction to reevaluate the population structure of M. tuberculosis, providing an in-depth analysis of the ancient Indo-Oceanic Lineage 1 and the modern Central Asian Lineage 3, and expanding our understanding of Lineages 2 and 4. We assess sub-lineages using genomic sequences from 4939 pan-susceptible strains, and find 30 new genetically distinct clades that we validate in a dataset of 4645 independent isolates. We find a consistent geographically restricted or unrestricted pattern for 20 groups, including three groups of Lineage 1. The distribution of terminal branch lengths across the M. tuberculosis phylogeny supports the hypothesis of a higher transmissibility of Lineages 2 and 4, in comparison with Lineages 3 and 1, on a global scale. We define an expanded barcode of 95 single nucleotide substitutions that allows rapid identification of 69 M. tuberculosis sub-lineages and 26 additional internal groups. Our results paint a higher resolution picture of the M. tuberculosis phylogeny and biogeography.

Genetic Profile of Linezolid-Resistant M. tuberculosis Clinical Strains from Moscow

Background: Linezolid, bedaquiline, and newer fluoroquinolones are currently placed as priority Group A drugs for the treatment of drug-resistant tuberculosis. The number of reported linezolid-resistant clinical strains is still low, and the correlation of molecular determinants with phenotype is not perfect. Methods: We determined the linezolid MICs for clinical isolates from the Moscow region and identified mutations in rplC and rrl genes. Results: All 16 linezolid-resistant isolates had previously reported mutations in the rplC or rrl loci, and 13 of them bore a RplC C154R substitution. Detection of this substitution in a heteroresistant state was not successful, probably, due to the more stable DNA secondary structure of the mutated fragment, which precludes its amplification in mixes with the wild-type DNA. Strains with an rplC mutation had higher linezolid MIC compared to isolates with rrl mutations. Conclusions: Linezolid resistance mostly emerged during treatment with the latest regimen. Three primary cases with linezolid resistance question the possible transmission of totally drug-resistant tuberculosis in the Moscow region, which demands further investigation.

Whole-genome comparative analysis at the lineage/sublineage level discloses relationships between Mycobacterium tuberculosis genotype and clinical phenotype

Background

Human tuberculosis (TB) caused by members of the Mycobacterium tuberculosis complex (MTBC) is the main cause of death among infectious diseases worldwide. Pulmonary TB (PTB) is the most common clinical phenotype of the disease, but some patients develop an extrapulmonary (EPTB) phenotype in which any organ or tissue can be affected. MTBC species include nine phylogenetic lineages, with some appearing globally and others being geographically restricted. EPTB can or not have pulmonary involvement, challenging its diagnosis when lungs are not implicated, thus causing an inadequate treatment. Finding evidence of a specific M. tuberculosis genetic background associated with EPTB is epidemiologically relevant due to the virulent and multidrug-resistant strains isolated from such cases. Until now, the studies conducted to establish associations between M. tuberculosis lineages and PTB/EPTB phenotypes have shown inconsistent results, which are attributed to the strain predominance from specific M. tuberculosis lineages/sublineages in the samples analyzed and the use of low-resolution phylogenetic tools that have impaired sublineage discrimination abilities. The present work elucidates the relationships between the MTBC strain lineages/sublineages and the clinical phenotypes of the disease as well as the antibiotic resistance of the strains.

Methods

To avoid biases, we retrieved the raw genomic reads (RGRs) of all (n = 245) the M. tuberculosis strains worldwide causing EPTB available in databases and an equally representative sample of the RGRs (n = 245) of PTB strains. A multiple alignment was constructed, and a robust maximum likelihood phylogeny based on single-nucleotide polymorphisms was generated, allowing effective strain lineage/sublineage assignment.

Results

A significant Odds Ratio (OR range: 1.8–8.1) association was found between EPTB and the 1.1.1, 1.2.1, 4.1.2.1 and ancestral Beijing sublineages. Additionally, a significant association between PTB with 4.3.1, 4.3.3, and 4.5 and Asian African 2 and Europe/Russia B0/W148 modern Beijing sublineages was found. We also observed a significant association of Lineage 3 strains with multidrug resistance (OR 3.8; 95% CI [1.1–13.6]), as well as between modern Beijing sublineages and antibiotic resistance (OR 4.3; 3.8–8.6). In this work, it was found that intralineage diversity can drive differences in the immune response that triggers the PTB/EPTB phenotype.

How Molecular Epidemiology Can Affect Tuberculosis Control in the Middle East Countries: A Systematic Review and Meta-Analysis

BACKGROUND

Nowadays, due to the incidence of specific strains of Mycobacterium tuberculosis and also an increase in the rate of drug resistant-TB, the mortality rate has been increased by this disease. The identification of common strains in the region, as well as the sources of transmission, is essential to control the disease, which is possible by using molecular epidemiology.

OBJECTIVES

In this survey, the studies utilizing the spoligotyping method in Muslim Middle East countries are reviewed to determine their role in the control of TB.

METHODS

All studies conducted from 2005 to June 2016 were considered systematically in three electronic databases out of which 23 studies were finally selected.

RESULTS

The average rate of clustering was 84% and the rate of recent transmission varied from 21.7% to 92.4%. The incidence of Beijing strains has been found to be rising in the abovementioned countries. In Iran and Saudi Arabia known as immigration and labour-hosting countries, respectively, rapid transmittable and drug-resistant Beijing strains were higher than those in other Muslim Middle East countries.

CONCLUSION

Considering the incidence of highly virulent strains, due to the increase in immigration and people infected with HIV, tuberculosis, especially drug-resistant form, careful monitoring is needed.

The rise and the fall of a Pseudomonas aeruginosa endemic lineage in a hospital

The biological features that allow a pathogen to survive in the hospital environment are mostly unknown. The extinction of bacterial epidemics in hospitals is mostly attributed to changes in medical practice, including infection control, but the role of bacterial adaptation has never been documented. We analysed a collection of Pseudomonas aeruginosa isolates belonging to the Besançon Epidemic Strain (BES), responsible for a 12year nosocomial outbreak, using a genotype-to-phenotype approach. Bayesian analysis estimated the emergence of the clone in the hospital 5 years before its opening, during the creation of its water distribution network made of copper. BES survived better than the reference strains PAO1 and PA14 in a copper solution due to a genomic island containing 13 metal-resistance genes and was specifically able to proliferate in the ubiquitous amoeba Vermamoeba vermiformis. Mutations affecting amino-acid metabolism, antibiotic resistance, lipopolysaccharide biosynthesis, and regulation were enriched during the spread of BES. Seven distinct regulatory mutations attenuated the overexpression of the genes encoding the efflux pump MexAB-OprM over time. The fitness of BES decreased over time in correlation with its genome size. Overall, the resistance to inhibitors and predators presumably aided the proliferation and propagation of BES in the plumbing system of the hospital. The pathogen further spread among patients via multiple routes of contamination. The decreased prevalence of patients infected by BES mirrored the parallel and convergent genomic evolution and reduction that affected bacterial fitness. Along with infection control measures, this may have participated in the extinction of BES in the hospital setting.

Whole genome sequencing of clinical samples reveals extensively drug resistant tuberculosis (XDR TB) strains from the Beijing lineage in Nigeria, West Africa

Multi-drug (MDR) and extensively drug-resistant (XDR) tuberculosis (TB) continues to be a global public health problem especially in high TB burden countries like Nigeria. Many of these cases are undetected and go on to infect high risk individuals. Clinical samples from positive rifampicin resistant Xpert®MTB/Rif assay were subjected to direct whole genome sequencing and bioinformatics analysis to identify the full antibiotics resistance and lineage profile. We report two (2) XDR TB samples also belonging to the East-Asian/Beijing family of lineage 2 Mycobacterium tuberculosis complex from clinical samples in Nigeria. Our findings further reveal the presence of mutations that confer resistance to first-line drugs (rifampicin, isoniazid, ethambutol and pyrazanimide), second-line injectables (capreomycin, streptomycin, kanamycin and/or amikacin) and at least one of the fluoroquinolones (ofloxacin, moxifloxacin, levofloxacin and/or ciprofloxacin) in both samples. The genomic sequence data from this study not only provide the first evidence of XDR TB in Nigeria and West Africa, but also emphasize the importance of WGS in accurately detecting MDR and XDR TB, to ensure adequate and proper management treatment regimens for affected individuals. This will greatly aid in preventing the spread of drug resistance TB in high burden countries.

MDR and Pre-XDR Clinical Mycobacterium tuberculosis Beijing Strains: Assessment of Virulence and Host Cytokine Response in Mice Infectious Model

Mycobacterium tuberculosis Beijing genotype associated with drug resistance is a growing public health problem worldwide. The aim of this study was the assessment of virulence for C57BL/6 mice after infection by clinical M. tuberculosis strains 267/47 and 120/26, which belong to the modern sublineages B0/W148 and Central Asia outbreak of the Beijing genotype, respectively. The sublineages were identified by the analysis of the strains' whole-genomes. The strains 267/47 and 120/26 were characterized as agents of pre-extensively drug-resistant (pre-XDR) and multidrug-resistant (MDR) tuberculosis, respectively. Both clinical strains were slow-growing in 7H9 broth compared to the M. tuberculosis H37Rv strain. The survival rates of C57BL/6 mice infected by 267/47, 120/26, and H37Rv on the 150th day postinfection were 10%, 40%, and 70%, respectively. Mycobacterial load in the lungs, spleen, and liver was higher and histopathological changes were more expressed for mice infected by the 267/47 strain compared to those infected by the 120/26 and H37Rv strains. The cytokine response in the lungs of C57BL/6 mice after infection with the 267/47, 120/26, and H37Rv strains was different. Notably, proinflammatory cytokine genes Il-1α, Il-6, Il-7, and Il-17, as well as anti-inflammatory genes Il-6 and Il-13, were downregulated after an infection caused by the 267/47 strain compared to those after infection with the H37Rv strain.

A genome epidemiological study of mycobacterium tuberculosis in subpopulations with high and low incidence rate in Guangxi, South China

BACKGROUND

Tuberculosis (TB) is caused by a bacterium called Mycobacterium tuberculosis (Mtb). China is the third in top 8 high TB burden countries and Guangxi is one of the high incidence areas in South China. Determine bacterial factors that affected TB incidence rate is a step toward Ending the TB epidemic.

RESULTS

Genomes of M. tuberculosis cultures from a relatively high and low incidence region in Guangxi have been sequenced. 347 of 358(96.9%) were identified as M. tuberculosis. All the strains belong to Lineage 2 and Lineage 4, except for one in Lineage 1. We found that the genetic structure of the M. tuberculosis population in each county varies enormously. Low incidence rate regions have a lower prevalence of Beijing genotypes than other regions. Four isolates which harbored mutT4-48 also had mutT2-58 mutations. It is suggested that strains from the ancestors of modern Beijing lineage is circulating in Guangxi. Strains of modern Beijing lineage (OR=2.04) were more likely to acquire drug resistances than Lineage 4. Most of the lineage differentiation SNPs are related to cell wall biosynthetic pathways.

CONCLUSIONS

These results provided a higher resolution to better understand the history of transmission of M. tuberculosis from/to South China. And the incidence rate of tuberculosis might be affected by bacterial population structure shaped by demographic history. Our findings also support the hypothesis that Modern Beijing lineage originated in South China.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

Tuberculosis (TB) has been responsible for the greatest number of human deaths due to an infectious disease in general, and due to antimicrobial resistance (AMR) in particular. The etiological agents of human TB are a closely-related group of human-adapted bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Understanding how MTBC populations evolve within-host may allow for improved TB treatment and control strategies. In this review, we highlight recent works that have shed light on how AMR evolves in MTBC populations within individual patients. We discuss the role of heteroresistance in AMR evolution, and review the bacterial, patient and environmental factors that likely modulate the magnitude of heteroresistance within-host. We further highlight recent works on the dynamics of MTBC genetic diversity within-host, and discuss how spatial substructures in patients' lungs, spatiotemporal heterogeneity in antimicrobial concentrations and phenotypic drug tolerance likely modulates the dynamics of MTBC genetic diversity in patients during treatment. We note the general characteristics that are shared between how the MTBC and other bacterial pathogens evolve in humans, and highlight the characteristics unique to the MTBC.

Accounting for the Biological Complexity of Pathogenic Fungi in Phylogenetic Dating

In the study of pathogen evolution, temporal dating of phylogenies provides information on when species and lineages may have diverged in the past. When combined with spatial and epidemiological data in phylodynamic models, these dated phylogenies can also help infer where and when outbreaks occurred, how pathogens may have spread to new geographic locations and/or niches, and how virulence or drug resistance has developed over time. Although widely applied to viruses and, increasingly, to bacterial pathogen outbreaks, phylogenetic dating is yet to be widely used in the study of pathogenic fungi. Fungi are complex organisms with several biological processes that could present issues with appropriate inference of phylogenies, clock rates, and divergence times, including high levels of recombination and slower mutation rates although with potentially high levels of mutation rate variation. Here, we discuss some of the key methodological challenges in accurate phylogeny reconstruction for fungi in the context of the temporal analyses conducted to date and make recommendations for future dating studies to aid development of a best practices roadmap in light of the increasing threat of fungal outbreaks and antifungal drug resistance worldwide.

Implementation of whole genome sequencing for tuberculosis diagnostics in a low-middle income, high MDR-TB burden country

Whole genome sequencing (WGS) is revolutionary for diagnostics of TB and its mutations associated with drug-resistances, but its uptake in low- and middle-income countries is hindered by concerns of implementation feasibility. Here, we provide a proof of concept for its successful implementation in such a setting. WGS was implemented in the Kyrgyz Republic. We estimated needs of up to 55 TB-WGS per week and chose the MiSeq platform (Illumina, USA) because of its capacity of up to 60 TB-WGS per week. The project's timeline was completed in 93-weeks. Costs of large equipment and accompanying costs were 222,065 USD and 8462 USD, respectively. The first 174 WGS costed 277 USD per sequence, but this was skewed by training inefficiencies. Based on real prices and presuming optimal utilization of WGS capacities, WGS costs could drop to 167 and 141 USD per WGS using MiSeq Reagent Kits v2 (500-cycles) and v3 (600-cycles), respectively. Five trainings were required to prepare the staff for autonomous WGS which cost 48,250 USD. External assessment confirmed excellent performance of WGS by the Kyrgyz laboratory in an interlaboratory comparison of 30 M. tuberculosis genomes showing complete agreeance of results.

Molecular epidemiology and drug susceptibility profiles of Mycobacterium tuberculosis complex isolates from Northern Ghana

OBJECTIVE

We conducted a cross-sectional study in the five administrative regions of Northern Ghana to determine the diversity of Mycobacterium tuberculosis complex (MTBC) sub/lineages and their susceptibility to isoniazid (INH) and rifampicin (RIF).

METHODS

Sputum specimens were collected and cultured from 566 pulmonary tuberculosis patients reporting to 17 health facilities from 2015 to 2019. Mycobacterial isolates obtained from solid cultures were confirmed as members of the MTBC by PCR amplification of IS6110 and rpoß and assigned lineages and sub-lineages using spoligotyping.

RESULTS

Of 294 mycobacterial isolates recovered, MTBC species identified were: M. tuberculosis sensu stricto (Mtbss) 241 (82.0%), M. africanum 41 (13.9%) and M. bovis four (1.4%) with eight (2.7%) unidentified. The human-adapted lineages (L) identified (N=279) were L1 (8/279, 2.9%), L2 (15/279, 5.4%), L3 (7/279, 2.5%), L4 (208/279, 74.5%), L5 (13/279, 4.7%) and L6 (28/279, 10.0%) with three unidentified lineages. Among the 208 L4, the dominant sub-lineages in the region were the Cameroon 120/208 (57.7%) and Ghana 50/208 (24.0%). We found 4.4% (13/294) and 0.7% (2/294) of the patients infected with MTBC isolates resistant to INH only and RIF only, respectively, with 2.4% (7/294) being infected with MDR strains. Whereas L6 was associated with the elderly, we identified that the Ghana sub-lineage of L4 was associated with both INH and MDR (p<0.05), making them important TB pathogens in Northern Ghana and a growing public health concern.

Genotypic diversity of multi- and pre-extremely drug-resistant Mycobacterium tuberculosis isolates from Morocco

In Morocco, the prevalence of multidrug resistant tuberculosis (MDR-TB) continues to increase especially within previously treated cases; these MDR cases may evolve to extensively drug resistant tuberculosis (XDR-TB) raising major concern to TB control programs. From an epidemiological window, scarce informations are available about the genetic diversity of Mycobacterium tuberculosis (MTB) strains fueling these forms of resistance. The aim of this study was to assess to genetic diversity of MDR-MTB strains. Hence, this prospective study was conducted on patients diagnosed with MDR-TB at Pasteur Institute of Casablanca from 2010 to 2013. A total of 70 MDR-MTB isolates were genotyped by spoligotyping and 15-loci MIRU-VNTR methods. Spoligotyping generated four orphan patterns, five unique profiles whereas 61 strains were grouped in nine clusters (2 to 25 strains per cluster), the clustering rates being 87.1%. Subtyping by 15 loci MIRU-VNTR splitted all clusters already established by spoligotyping and generated 70 unique profiles not recognized in SITVIT2 database; clustering rate was equal to zero. HGDI analysis of 15 loci MIRU demonstrated that eight out of 15 loci were highly discriminant. Of note, all pre-XDR strains belongs to many clades, meaning that there no association between gyrA mutants and particular clade. Overall, the data generated by this study (i) describe the population structure of MDR MTBC in Morocco which is highly homogenous, (ii) confirm that TB in Morocco is almost exclusively transmitted by modern and evolutionary lineages with high level of biodiversity seen by MIRU, and (iii) validate the use of optimized 15-loci MIRU-VNTR format for future investigations in Morocco.

Phenotypic and genotypic features of the Mycobacterium tuberculosis lineage 1 subgroup in central Vietnam

Mycobacterium tuberculosis (Mtb) has different features depending on different geographic areas. We collected Mtb strains from patients with smear-positive pulmonary tuberculosis in Da Nang, central Vietnam. Using a whole genome sequencing platform, including genome assembly complemented by long-read-sequencing data, genomic characteristics were studied. Of 181 Mtb isolates, predominant Vietnamese EAI4_VNM and EAI4-like spoligotypes (31.5%), ZERO strains (5.0%), and part of EAI5 (11.1%) were included in a lineage-1 (L1) sublineage, i.e., L1.1.1.1. These strains were found less often in younger people, and they genetically clustered less frequently than other modern strains. Patients infected with ZERO strains demonstrated less lung infiltration. A region in RD2bcg spanning six loci, i.e., PE_PGRS35, cfp21, Rv1985c, Rv1986, Rv1987, and erm(37), was deleted in EAI4_VNM, EAI4-like, and ZERO strains, whereas another 118 bp deletion in furA was specific only to ZERO strains. L1.1.1.1-sublineage-specific deletions in PE_PGRS4 and PE_PGRS22 were also identified. RD900, seen in ancestral lineages, was present in majority of the L1 members. All strains without IS6110 (5.0%) had the ZERO spoligo-pattern. Distinctive features of the ancestral L1 strains provide a basis for investigation of the modern versus ancestral Mtb lineages and allow consideration of countermeasures against this heterogeneous pathogen.

Characterization of genetic diversity and clonal complexes by whole genome sequencing of Mycobacterium tuberculosis isolates from Jalisco, Mexico

Whole genome sequencing (WGS) analysis in tuberculosis allows the prediction of drug-resistant phenotypes, identification of lineages, and to better understanding of the epidemiology and transmission chains. Nevertheless the procedure has been scarcely assessed in Mexico, in this work we analyze by WGS isolates of Mycobacterium tuberculosis circulating in Jalisco, Mexico. Lineage and phylogenetic characterization, drug resistant prediction, "in silico" spoligotyping determination, were provided by WGS in 32 M. tuberculosis clinical isolates. Lineage 4 (L4), with 28 isolates (87%) and eleven sublineages was dominant. Forty SNPs and INDELs were found in genes related to first-, and second-line drugs. Eleven isolates were sensitive, seven (22%) were predicted to be resistant to isoniazid, two resistant to rifampicin (6%) and two (6%) were multidrug-resistant tuberuclosis. Spoligotyping shows that SIT 53 (19%) and SIT 119 (16%) were dominant. Four clonal transmission complexes were found. This is the first molecular epidemiological description of TB isolates circulating in western Mexico, achieved through WGS. L4 was dominant and included a high diversity of sublineages. It was possible to track the transmission route of two clonal complexes. The WGS demonstrated to be of great utility and with further implications for clinical and epidemiological study of TB in the region.

Multiple Merger Genealogies in Outbreaks of Mycobacterium tuberculosis

The Kingman coalescent and its developments are often considered among the most important advances in population genetics of the last decades. Demographic inference based on coalescent theory has been used to reconstruct the population dynamics and evolutionary history of several species, including Mycobacterium tuberculosis (MTB), an important human pathogen causing tuberculosis. One key assumption of the Kingman coalescent is that the number of descendants of different individuals does not vary strongly, and violating this assumption could lead to severe biases caused by model misspecification. Individual lineages of MTB are expected to vary strongly in reproductive success because 1) MTB is potentially under constant selection due to the pressure of the host immune system and of antibiotic treatment, 2) MTB undergoes repeated population bottlenecks when it transmits from one host to the next, and 3) some hosts show much higher transmission rates compared with the average (superspreaders).

Here, we used an approximate Bayesian computation approach to test whether multiple-merger coalescents (MMC), a class of models that allow for large variation in reproductive success among lineages, are more appropriate models to study MTB populations. We considered 11 publicly available whole-genome sequence data sets sampled from local MTB populations and outbreaks and found that MMC had a better fit compared with the Kingman coalescent for 10 of the 11 data sets. These results indicate that the null model for analyzing MTB outbreaks should be reassessed and that past findings based on the Kingman coalescent need to be revisited.

Current advances in drug delivery of nanoparticles for respiratory disease treatment

Cases of respiratory diseases have been increasing around the world, affecting the health and quality of life of millions of people every year. Chronic respiratory diseases (CRDs) and acute respiratory infections (ARIs) are responsible for many hospital admissions and deaths, requiring sophisticated treatments that facilitate the delivery of therapeutics to specific target sites with controlled release. In this context, different nanoparticles (NPs) have been explored to match this demand, such as lipid, liposome, protein, carbon-based, polymeric, metallic, oxide, and magnetic NPs. The use of NPs as drug delivery systems can improve the efficacy of commercial drugs due to their advantages related to sustained drug release, targeting effects, and patient compliance. The current review presents an updated summary of recent advances regarding the use of NPs as drug delivery systems to treat diseases related to the respiratory tract, such as CRDs and ARIs. The latest applications presented in the literature were considered, and the opportunities and challenges of NPs in the drug delivery field are discussed.

Factors influencing the effectiveness of the Gudair vaccine for controlling Johne's disease in sheep flocks in Australia

Ovine Johne's disease is a chronic debilitating disease of sheep caused by Mycobacterium avium subsp. paratuberculosis (Mptb) which results in diarrhoea, emaciation and mortalities in infected animals. Vaccination with Gudair® has been a key strategy for controlling the disease in Australia since its approval in 2002. Previous research conducted in Australia has demonstrated that the vaccine is quite effective in reducing sheep mortalities. While some farms have also been successful in reducing the prevalence of the disease in their flocks to undetectable levels, sheep in other flocks continue to shed Mptb in faeces even after an ongoing vaccination program . This study was conducted to investigate management, husbandry and biosecurity factors associated with paratuberculosis infection in Gudair® vaccinated sheep flocks in Australia. We enrolled 64 sheep farmers and interviewed them to obtain information about their management and biosecurity practices. Pooled faecal samples were collected from sheep at each farm and cultured to create two outcome variables: Mptb positive (yes/no) and disease prevalence level (nil, < 1 %, ≥ 1 %). Binary and ordinal logistic regression analyses were conducted to evaluate the association of management, husbandry and biosecurity factors with these outcome variables. Farms were more likely to have Mptb positive sheep and a higher disease prevalence in their flocks if they: (a) provided supplementary feed on the ground (instead of in a trough); (b) had a greater number of neighbours with sheep; and (c) had introduced rams from a greater number of sources. The results suggest the effectiveness of Gudair® vaccination to control OJD can be improved if sheep producers maintain other risk management strategies and biosecurity practices. Extension agencies should advise farmers not to relax their biosecurity practices and to purchase rams from only low-risk sources, even if they are continuing to vaccinate their flocks.

Improving tuberculosis surveillance by detecting international transmission using publicly available whole genome sequencing data

IntroductionImproving the surveillance of tuberculosis (TB) is especially important for multidrug-resistant (MDR) and extensively drug-resistant (XDR) TB. The large amount of publicly available whole genome sequencing (WGS) data for TB gives us the chance to re-use data and to perform additional analyses at a large scale.AimWe assessed the usefulness of raw WGS data of global MDR/XDR Mycobacterium tuberculosis isolates available from public repositories to improve TB surveillance.MethodsWe extracted raw WGS data and the related metadata of M. tuberculosis isolates available from the Sequence Read Archive. We compared this public dataset with WGS data and metadata of 131 MDR- and XDR M. tuberculosis isolates from Germany in 2012 and 2013.ResultsWe aggregated a dataset that included 1,081 MDR and 250 XDR isolates among which we identified 133 molecular clusters. In 16 clusters, the isolates were from at least two different countries. For example, Cluster 2 included 56 MDR/XDR isolates from Moldova, Georgia and Germany. When comparing the WGS data from Germany with the public dataset, we found that 11 clusters contained at least one isolate from Germany and at least one isolate from another country. We could, therefore, connect TB cases despite missing epidemiological information.ConclusionWe demonstrated the added value of using WGS raw data from public repositories to contribute to TB surveillance. Comparing the German with the public dataset, we identified potential international transmission events. Thus, using this approach might support the interpretation of national surveillance results in an international context.