MDR-TB Outbreak among HIV-Negative Tunisian Patients followed during 11 Years

Delineating the evolutionary pathway to multidrug-resistant outbreaks of a Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage

OBJECTIVES

We sought to capture the evolutionary itinerary of the Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, where it caused a major multidrug-resistant (MDR) tuberculosis outbreak in a context strictly negative for HIV infection.

METHODS

We combined whole genome sequencing and Bayesian approaches using a representative collection of drug-susceptible and drug-resistant L4.1.2.1/Haarlem clinical strains (n = 121) recovered from the outbreak region over 16 years.

RESULTS

In the absence of drug resistance, the L4.1.2.1/Haarlem sublineage showed a propensity for rapid transmission as witnessed by the high clustering (44.6%) and recent transmission rates (25%), as well as the reduced mean distance between genome pairs. The entire pool of L4.1.2.1/Haarlem MDR strains was found to be linked to either the aforementioned major outbreak (68 individuals, 2001-2016) or to a minor, newly uncovered outbreak (six cases, 2001-2011). Strikingly, the two outbreaks descended independently from a common ancestor that can be dated back to 1886.

CONCLUSIONS

Our data point to the intrinsic propensity for rapid transmission of the M. tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, linking the overall MDR tuberculosis epidemic to a single ancestor. These findings bring out the important role of the bacillus' genetic background in the emergence of successful MDR M. tuberculosis clones.

On the onset and dispersal of a major MDR TB clone among HIV-negative patients, Tunisia

BACKGROUND

To carry out a whole genome sequencing (WGS)-based investigation on the emergence and spread of the largest multidrug-resistant tuberculosis (MDR TB) outbreak that has been thriving among HIV-negative patients, Tunisia, since the early 2000s.

METHODS

We performed phylogeographic analyses and molecular dating based on a WGS dataset representing 68 unique Mycobacterium tuberculosis isolates, covering almost the entire MDR TB outbreak for the time period 2001-2016.

RESULTS

The data indicate that the ancestor of the MDR TB outbreak emerged in the region of Bizerte, as early as 1974 (95% CI 1951-1985), from where it spread to other regions by 1992 (95% CI 1980-1996). Analysis of a minimum spanning tree based on core genome Multilocus Sequence Typing (cgMLST) uncovered the early spill-over of the fitness-compensated MDR TB strain from the prison into the general population. Indeed, cases with history of incarceration were found to be directly or indirectly linked to up to 22 new outbreak cases (32.35%) among the non-imprisoned population. By around 2008, the MDR TB outbreak strain had acquired additional resistance, leading to an XDR phenotype.

CONCLUSIONS

WGS allowed refining our understanding of the emergence and evolution of the largest MDR TB outbreak in Tunisia, whose causative strain has been circulating silently for almost 26 years before. Our study lends further support to the critical role of prisons-related cases in the early spread of the outbreak among the general population. The shift to an XDR phenotype of such an epidemic clone prompts an urgent need to undertake drastic control measures.

Abattoir survey of bovine tuberculosis in tanta, centre of the Nile delta, with in silico analysis of gene mutations and protein-protein interactions of the involved mycobacteria

Bovine tuberculosis is a transboundary disease of high economic and public health burden worldwide. In this study, post-mortem examination of 750 cattle and buffalo in Tanta abattoir, Centre of the Nile Delta, revealed visible TB in 4% of animals and a true prevalence of 6.85% (95% CI: 5.3%-8.9%). Mycobacterial culture, histopathology and RT-PCR targeting all members of M. tuberculosis complex were performed, upon which 85%, 80% and 100% of each tested lesions were confirmed as TB, respectively. Mpb70-targeting PCR was conducted on ten RT-PCR positive samples for sequencing and identified nine Mycobacterium (M.) bovis strains and, interestingly, one M. tuberculosis (Mtb) strain from a buffalo. Bioinformatics tools were used for prediction of mutations, nucleotide polymorphisms, lineages, drug resistance and protein-protein interactions (PPI) of the sequenced strains. The Mtb strain was resistant to rifampicin, isoniazid and streptomycin, and to the best of our knowledge, this is the first report of multidrug resistant (MDR)-Mtb originating from buffaloes. Seven M. bovis strains were resistant to ethambutol and ethionamide. Such resistances were associated with KatG, rpoB, rpsL, embB and ethA genes mutations. Other mutations and nucleotide polymorphisms were also predicted, some are reported for the first time and require experimental work for validation. PPI revealed more interactions than what would be expected for a random set of proteins of similar size and had dense interactions between nodes that are biologically connected, as a group. Two M. bovis strains belonged to BOV AFRI lineage (Spoligotypes BOV 1; BOV 2) and eight strains belonged to East-Asian (Beijing) lineage. In conclusion, visible TB was prevalent in the study area, RT-PCR is the best to confirm the disease, MDR-Mtb is associated with buffalo TB, and mycobacteria of different lineages carry many resistance genes to chemotherapeutic agents used in treatment of human TB constituting a major public health risk.

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.

Whole-Genome Sequencing of Drug-Resistant Mycobacterium tuberculosis Strains, Tunisia, 2012-2016

To investigate transmission of drug-resistant strains of Mycobacterium tuberculosis in Tunisia, we performed whole-genome sequencing on 46 multidrug-resistant strains isolated during 2012-2016. Core-genome multilocus sequence typing grouped 30 strains (65.2%) into 3 clusters, indicating extensive recent transmission and Haarlem clone predominance. Whole-genome sequencing might help public health services undertake appropriate control actions.

Automated IS6110-based fingerprinting of Mycobacterium tuberculosis: Reaching unprecedented discriminatory power and versatility

BACKGROUND

Several technical hurdles and limitations have restricted the use of IS6110 restriction fragment length polymorphism (IS6110 RFLP), the most effective typing method for detecting recent tuberculosis (TB) transmission events. This has prompted us to conceive an alternative modality, IS6110-5'3'FP, a plasmid-based cloning approach coupled to a single PCR amplification of differentially labeled 5' and 3' IS6110 polymorphic ends and their automated fractionation on a capillary sequencer. The potential of IS6110-5'3'FP to be used as an alternative to IS6110 RFLP has been previously demonstrated, yet further technical improvements are still required for optimal discriminatory power and versatility.

OBJECTIVES

Here we introduced critical amendments to the original IS6110-5'3'FP protocol and compared its performance to that of 24-loci multiple interspersed repetitive unit-variable number tandem repeats (MIRU-VNTR), the current standard method for TB transmission analyses.

METHODS

IS6110-5'3'FP protocol modifications involved: (i) the generation of smaller-sized polymorphic fragments for efficient cloning and PCR amplification, (ii) omission of the plasmid amplification step in E. coli for shorter turnaround times, (iii) the use of more stable fluorophores for increased sensitivity, (iv) automated subtraction of background fluorescent signals, and (v) the automated conversion of fluorescent peaks into binary data.

RESULTS

In doing so, the overall turnaround time of IS6110-5'3'FP was reduced to 4 hours. The new protocol allowed detecting almost all 5' and 3' IS6110 polymorphic fragments of any given strain, including IS6110 high-copy number Beijing strains. IS6110-5'3'FP proved much more discriminative than 24-loci MIRU-VNTR, particularly with strains of the M. tuberculosis lineage 4.

CONCLUSIONS

The IS6110-5'3'FP protocol described herein reached the optimal discriminatory potential of IS6110 fingerprinting and proved more accurate than 24-loci MIRU-VNTR in estimating recent TB transmission. The method, which is highly cost-effective, was rendered versatile enough to prompt its evaluation as an automatized solution for a TB integrated molecular surveillance.

Genotyping and spatial analysis of pulmonary tuberculosis and diabetes cases in the state of Veracruz, Mexico

Background

Genotyping and georeferencing in tuberculosis (TB) have been used to characterize the distribution of the disease and occurrence of transmission within specific groups and communities.

Objective

The objective of this study was to test the hypothesis that diabetes mellitus (DM) and pulmonary TB may occur in spatial and molecular aggregations.

Material and methods

Retrospective cohort study of patients with pulmonary TB. The study area included 12 municipalities in the Sanitary Jurisdiction of Orizaba, Veracruz, México. Patients with acid-fast bacilli in sputum smears and/or Mycobacterium tuberculosis in sputum cultures were recruited from 1995 to 2010. Clinical (standardized questionnaire, physical examination, chest X-ray, blood glucose test and HIV test), microbiological, epidemiological, and molecular evaluations were carried out. Patients were considered “genotype-clustered” if two or more isolates from different patients were identified within 12 months of each other and had six or more IS6110 bands in an identical pattern, or < 6 bands with identical IS6110 RFLP patterns and spoligotype with the same spacer oligonucleotides. Residential and health care centers addresses were georeferenced. We used a Jeep hand GPS. The coordinates were transferred from the GPS files to ArcGIS using ArcMap 9.3. We evaluated global spatial aggregation of patients in IS6110-RFLP/ spoligotype clusters using global Moran´s I. Since global distribution was not random, we evaluated “hotspots” using Getis-Ord Gi* statistic. Using bivariate and multivariate analysis we analyzed sociodemographic, behavioral, clinic and bacteriological conditions associated with “hotspots”. We used STATA® v13.1 for all statistical analysis.

Results

From 1995 to 2010, 1,370 patients >20 years were diagnosed with pulmonary TB; 33% had DM. The proportion of isolates that were genotyped was 80.7% (n = 1105), of which 31% (n = 342) were grouped in 91 genotype clusters with 2 to 23 patients each; 65.9% of total clusters were small (2 members) involving 35.08% of patients. Twenty three (22.7) percent of cases were classified as recent transmission. Moran`s I indicated that distribution of patients in IS6110-RFLP/spoligotype clusters was not random (Moran`s I = 0.035468, Z value = 7.0, p = 0.00). Local spatial analysis showed statistically significant spatial aggregation of patients in IS6110-RFLP/spoligotype clusters identifying “hotspots” and “coldspots”. GI* statistic showed that the hotspot for spatial clustering was located in Camerino Z. Mendoza municipality; 14.6% (50/342) of patients in genotype clusters were located in a hotspot; of these, 60% (30/50) lived with DM. Using logistic regression the statistically significant variables associated with hotspots were: DM [adjusted Odds Ratio (aOR) 7.04, 95% Confidence interval (CI) 3.03–16.38] and attending the health center in Camerino Z. Mendoza (aOR18.04, 95% CI 7.35–44.28).

Conclusions

The combination of molecular and epidemiological information with geospatial data allowed us to identify the concurrence of molecular clustering and spatial aggregation of patients with DM and TB. This information may be highly useful for TB control programs.

Molecular clustering of patients with diabetes and pulmonary tuberculosis: A systematic review and meta-analysis

Introduction

Many studies have explored the relationship between diabetes mellitus (DM) and tuberculosis (TB) demonstrating increased risk of TB among patients with DM and poor prognosis of patients suffering from the association of DM/TB. Owing to a paucity of studies addressing this question, it remains unclear whether patients with DM and TB are more likely than TB patients without DM to be grouped into molecular clusters defined according to the genotype of the infecting Mycobacterium tuberculosis bacillus. That is, whether there is convincing molecular epidemiological evidence for TB transmission among DM patients. Objective: We performed a systematic review and meta-analysis to quantitatively evaluate the propensity for patients with DM and pulmonary TB (PTB) to cluster according to the genotype of the infecting M. tuberculosis bacillus.

Materials and methods

We conducted a systematic search in MEDLINE and LILACS from 1990 to June, 2016 with the following combinations of key words “tuberculosis AND transmission” OR “tuberculosis diabetes mellitus” OR “Mycobacterium tuberculosis molecular epidemiology” OR “RFLP-IS6110” OR “Spoligotyping” OR “MIRU-VNTR”. Studies were included if they met the following criteria: (i) studies based on populations from defined geographical areas; (ii) use of genotyping by IS6110- restriction fragment length polymorphism (RFLP) analysis and spoligotyping or mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTR) or other amplification methods to identify molecular clustering; (iii) genotyping and analysis of 50 or more cases of PTB; (iv) study duration of 11 months or more; (v) identification of quantitative risk factors for molecular clustering including DM; (vi) > 60% coverage of the study population; and (vii) patients with PTB confirmed bacteriologically. The exclusion criteria were: (i) Extrapulmonary TB; (ii) TB caused by nontuberculous mycobacteria; (iii) patients with PTB and HIV; (iv) pediatric PTB patients; (v) TB in closed environments (e.g. prisons, elderly homes, etc.); (vi) diabetes insipidus and (vii) outbreak reports. Hartung-Knapp-Sidik-Jonkman method was used to estimate the odds ratio (OR) of the association between DM with molecular clustering of cases with TB. In order to evaluate the degree of heterogeneity a statistical Q test was done. The publication bias was examined with Begg and Egger tests. Review Manager 5.3.5 CMA v.3 and Biostat and Software package R were used.

Results

Selection criteria were met by six articles which included 4076 patients with PTB of which 13% had DM. Twenty seven percent of the cases were clustered. The majority of cases (48%) were reported in a study in China with 31% clustering. The highest incidence of TB occurred in two studies from China. The global OR for molecular clustering was 0.84 (IC 95% 0.40–1.72). The heterogeneity between studies was moderate (I² = 55%, p = 0.05), although there was no publication bias (Beggs test p = 0.353 and Eggers p = 0.429).

Conclusion

There were very few studies meeting our selection criteria. The wide confidence interval indicates that there is not enough evidence to draw conclusions about the association. Clustering of patients with DM in TB transmission chains should be investigated in areas where both diseases are prevalent and focus on specific contexts.