Patients with active tuberculosis often have different strains in the same sputum specimen

TBtypeR: Sensitive detection and sublineage classification of Mycobacterium tuberculosis complex mixed-strain infections

Tuberculosis (TB) is typically attributed to a single infecting strain of the Mycobacterium tuberculosis Complex (MTBC), however mixed-strain infections are more common than reported due to the limitations of conventional diagnostics. While whole genome sequencing (WGS) methods have improved the detection of mixed-strain infections, existing tools struggle to reliably identify mixed-strain infections with frequencies below 10%. TBtypeR, a new tool, addresses this challenge by comparing WGS data to a phylogenetic SNP panel of over 10,000 sites and 164 MTBC phylotypes and using a model based on the binomial distribution to classify sublineage mixtures at frequencies as low as 1%. Extensive benchmarking shows TBtypeR outperforms current methods. Application to a published dataset of 5000 WGS samples identified 305 mixed-strain infections, six-fold higher than previously reported. The TBtypeR R package and a Nextflow pipeline are available at github.com/bahlolab/TBtypeR, providing a powerful tool for studying TB epidemiology and mixed-strain infections.

Genetic stability of Mycobacterium smegmatis under the stress of first-line antitubercular agents

The sustained success of Mycobacterium tuberculosis as a pathogen arises from its ability to persist within macrophages for extended periods and its limited responsiveness to antibiotics. Furthermore, the high incidence of resistance to the few available antituberculosis drugs is a significant concern, especially since the driving forces of the emergence of drug resistance are not clear. Drug-resistant strains of Mycobacterium tuberculosis can emerge through de novo mutations, however, mycobacterial mutation rates are low. To unravel the effects of antibiotic pressure on genome stability, we determined the genetic variability, phenotypic tolerance, DNA repair system activation, and dNTP pool upon treatment with current antibiotics using Mycobacterium smegmatis. Whole-genome sequencing revealed no significant increase in mutation rates after prolonged exposure to first-line antibiotics. However, the phenotypic fluctuation assay indicated rapid adaptation to antibiotics mediated by non-genetic factors. The upregulation of DNA repair genes, measured using qPCR, suggests that genomic integrity may be maintained through the activation of specific DNA repair pathways. Our results, indicating that antibiotic exposure does not result in de novo adaptive mutagenesis under laboratory conditions, do not lend support to the model suggesting antibiotic resistance development through drug pressure-induced microevolution.

Effect of mixed Mycobacterium tuberculosis infection on rapid molecular diagnostics among patients starting MDR-TB treatment in Uganda

BACKGROUND

Mixed M. tuberculosis (MTB) infection occurs when one is infected with more than one clonally distinct MTB strain. This form of infection can assist MTB strains to acquire additional mutations, facilitate the spread of drug-resistant strains, and boost the rate of treatment failure. Hence, the presence of mixed MTB infection could affect the performance of some rapid molecular diagnostic tests such as Line Probe Assay (LPA) and GeneXpert MTB/RIF (Xpert) assays.

METHODS

This was a cross-sectional study that used sputum specimens collected from participants screened for STREAM 2 clinical trial between October 2017 and October 2019. Samples from 62 MTB smear-positive patients and rifampicin-resistant patients from peripheral health facilities were processed for Xpert and LPA as screening tests for eligibility in the trial. From November 2020, processed stored sputum samples were retrieved and genotyped to determine the presence of mixed-MTB strain infection using a standard 24-locus Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem-Repeat (MIRU-VNTR). Samples with at least 20/24 MIRU-VNTR loci amplified were considered for analysis. Agar proportional Drug Susceptibility Test (DST) was performed on culture isolates of samples that had discordant results between LPA and Xpert. The impact of the presence of mixed-MTB strain on Xpert and LPA test interpretation was analyzed.

RESULTS

A total of 53/62 (85%) samples had analyzable results from MIRU-VNTR. The overall prevalence of mixed-MTB infection was 5/53 (9.4%). The prevalence was highest among male's 3/31 (9.7%) and among middle-aged adults, 4/30 (33.3%). Lineage 4 of MTB contributed 3/5 (60.0%) of the mixed-MTB infection prevalence. Having mixed MTB strain infection increased the odds of false susceptible Xpert test results (OR 7.556, 95% CI 0.88-64.44) but not for LPA. Being HIV-positive (P = 0.04) independently predicted the presence of mixed MTB infection.

CONCLUSIONS

The presence of mixed-MTB strain infection may affect the performance of the GeneXpert test but not for LPA. For patients with high pre-test probability of rifampicin resistance, an alternative rapid method such as LPA should be considered.

Lineage classification and antitubercular drug resistance surveillance of Mycobacterium tuberculosis by whole-genome sequencing in Southern India

IMPORTANCE

Studies mapping genetic heterogeneity of clinical isolates of M. tuberculosis for determining their strain lineage and drug resistance by whole-genome sequencing are limited in high tuberculosis burden settings. We carried out whole-genome sequencing of 242 M. tuberculosis isolates from drug-sensitive and drug-resistant tuberculosis patients, identified and collected as part of the TB Portals Program, to have a comprehensive insight into the genetic diversity of M. tuberculosis in Southern India. We report several genetic variations in M. tuberculosis that may confer resistance to antitubercular drugs. Further wide-scale efforts are required to fully characterize M. tuberculosis genetic diversity at a population level in high tuberculosis burden settings for providing precise tuberculosis treatment.

Molecular investigations of Mycobacterium tuberculosis genotypes among baseline and follow-up strains circulating in four regions of Eswatini

BACKGROUND

The tuberculosis (TB) epidemic remains a major global health problem and Eswatini is not excluded. Our study investigated the circulating genotypes in Eswatini and compared them at baseline (start of treatment) and follow-up during TB treatment.

METHODS

Three hundred and ninety (n = 390) participants were prospectively enrolled from referral clinics and patients who met the inclusion criteria, were included in the study. A total of 103 participants provided specimens at baseline and follow-up within six months. Mycobacterium tuberculosis (M.tb) strains were detected by GeneXpert® MTB/RIF assay (Cephied, USA) and Ziehl -Neelsen (ZN) microscopy respectively at baseline and follow-up time-points respectively. The 206 collected specimens were decontaminated and cultured on BACTEC™ MGIT™ 960 Mycobacteria Culture System (Becton Dickinson, USA). Drug sensitivity testing was performed at both baseline and follow-up time points. Spoligotyping was performed on both baseline and follow-up strains after DNA extraction.

RESULTS

Resistance to at least one first line drug was detected higher at baseline compared to follow-up specimens with most of them developing into multidrug-resistant (MDR)-TB. A total of four lineages and twenty genotypes were detected. The distribution of the lineages varied among the different regions in Eswatini. The Euro-American lineage was the most prevalent with 46.12% (95/206) followed by the East Asian with 24.27% (50/206); Indo-Oceanic at 9.71% (20/206) and Central Asian at 1.94% (4/206). Furthermore, a high proportion of the Beijing genotype at 24.27% (50/206) and S genotype at 16.50% (34/206) were detected. The Beijing genotype was predominant in follow-up specimens collected from the Manzini region with 48.9% (23/47) (p = 0.001). A significant proportion of follow-up specimens developed MDR-TB (p = 0.001) with Beijing being the major genotype in most follow-up specimens (p < 0.000).

CONCLUSION

Eswatini has a high M.tb genotypic diversity. A significant proportion of the TB infected participants had the Beijing genotype associated with MDR-TB in follow-up specimens and thus indicate community wide transmission.

Mixed strain pathogen populations accelerate the evolution of antibiotic resistance in patients

Antibiotic resistance poses a global health threat, but the within-host drivers of resistance remain poorly understood. Pathogen populations are often assumed to be clonal within hosts, and resistance is thought to emerge due to selection for de novo variants. Here we show that mixed strain populations are common in the opportunistic pathogen P. aeruginosa. Crucially, resistance evolves rapidly in patients colonized by multiple strains through selection for pre-existing resistant strains. In contrast, resistance evolves sporadically in patients colonized by single strains due to selection for novel resistance mutations. However, strong trade-offs between resistance and growth rate occur in mixed strain populations, suggesting that within-host diversity can also drive the loss of resistance in the absence of antibiotic treatment. In summary, we show that the within-host diversity of pathogen populations plays a key role in shaping the emergence of resistance in response to treatment.

Association of Single-Nucleotide Polymorphisms in the VDR Gene with Tuberculosis and Infection of Beijing Genotype Mycobacterium tuberculosis

Background

The aim of the present study was to investigate the association between vitamin D receptor (VDR) gene polymorphism and tuberculosis susceptibility, as well as the potential interaction of host genetic factors with the heterogeneity of Mycobacterium tuberculosis in the population from Xinjiang, China.

Methods

From January 2019 to January 2020, we enrolled 221 tuberculosis patients as the case group and 363 staff with no clinical symptoms as the control group from four designated tuberculosis hospitals in southern Xinjiang, China. The polymorphisms of Fok I, Taq I, Apa I, Bsm I, rs3847987 and rs739837 in the VDR were detected by sequencing. M. tuberculosis isolates were collected from the case group and identified as Beijing or non-Beijing lineage by multiplex PCR. Propensity score (PS), univariate analysis and multivariable logistic regression models were used to perform the analysis.

Results

Our results showed that the allele and genotype frequencies of Fok I, Taq I, Apa I, Bsm I, rs3847987 and rs739837 in VDR were not correlated with tuberculosis susceptibility or lineages of M. tuberculosis. Two out of six loci of the VDR gene formed one haplotype block, and none of the haplotypes was found to correlate with tuberculosis susceptibility or lineages of M. tuberculosis infected.

Conclusion

Polymorphisms in the VDR gene may not indicate susceptibility to tuberculosis. There was also no evidence on the interaction between the VDR gene of host and the lineages of M. tuberculosis in the population from Xinjiang, China. Further studies are nonetheless required to prove our conclusions.

A Current Perspective on the Potential of Nanomedicine for Anti-Tuberculosis Therapy

Tuberculosis (TB) is one of the ten infectious diseases that cause the highest amount of human mortality and morbidity. This infection, which is caused by a single pathogen, Mycobacterium tuberculosis, kills over a million people every year. There is an emerging problem of antimicrobial resistance in TB that needs urgent treatment and management. Tuberculosis treatment is complicated by its complex drug regimen, its lengthy duration and the serious side-effects caused by the drugs required. There are a number of critical issues around drug delivery and subsequent intracellular bacterial clearance. Drugs have a short lifespan in systemic circulation, which limits their activity. Nanomedicine in TB is an emerging research area which offers the potential of effective drug delivery using nanoparticles and a reduction in drug doses and side-effects to improve patient compliance with the treatment and enhance their recovery. Here, we provide a minireview of anti-TB treatment, research progress on nanomedicine and the prospects for future applications in developing innovative therapies.

Check your assumptions: Further scrutiny of basic model frameworks of antimicrobial resistance

Since the mid-1990s, growing concerns over antimicrobial resistant (AMR) organisms has led to an increase in the use of mathematical models to explore the inter-host transmission of such infections. Previous work reviewing such models categorised them into generic frameworks based on their underlying assumptions. These assumptions dictated the coexistence between AMR and antimicrobial sensitive strains. We add to this work performing stability analyses of the frameworks, along with simulating them deterministically and stochastically. Stability analyses found that many of these assumptions lead to models having the same equilibria, but showed differences in the equilibria's stability between models. Deterministic simulations reveal that assuming replacement of one infecting strain by another leads to an unusual antimicrobial treatment threshold. Increasing beyond this threshold causes a discontinuous increase in disease burden. The cost of AMR to pathogen fitness (lowered transmission) dictates both the threshold of treatment that causes the discontinuous increase in disease burden and the size of that increase. It was also shown that Superinfection states can be biased against resident strains and so favour coexistence of both strains. Stochastic simulations demonstrated that differing scenario starting conditions can guide models to converge upon equilibria that they may not have under deterministic simulation. These findings highlight the importance of checking assumptions when modelling AMR and strain competition more widely.

The detection of mixed tuberculosis infections using culture filtrate and resuscitation promoting factor deficient filtrate

Tuberculosis (TB) infected individuals harbor a heterogenous population of differentially culturable tubercle bacilli (DCTB). Herein, we describe how DCTB assays using culture filtrate either containing or deficient in resuscitation promoting factors can uncover mixed infections. We demonstrate that Mycobacterium tuberculosis (Mtb) strain genotypes can be separated in DCTB assays based on their selective requirement for growth stimulatory factors. Beijing mixed infections appear to be associated with a higher bacterial load and reduced reliance on growth stimulatory factors. These data have important implications for identifying mixed infections and hetero-resistance, which in turn can affect selection of treatment regimen and establishment of transmission links.

Fluoroquinolone heteroresistance, antimicrobial tolerance, and lethality enhancement

With tuberculosis, the emergence of fluoroquinolone resistance erodes the ability of treatment to interrupt the progression of MDR-TB to XDR-TB. One way to reduce the emergence of resistance is to identify heteroresistant infections in which subpopulations of resistant mutants are likely to expand and make the infections fully resistant: treatment modification can be instituted to suppress mutant enrichment. Rapid DNA-based detection methods exploit the finding that fluoroquinolone-resistant substitutions occur largely in a few codons of DNA gyrase. A second approach for restricting the emergence of resistance involves understanding fluoroquinolone lethality through studies of antimicrobial tolerance, a condition in which bacteria fail to be killed even though their growth is blocked by lethal agents. Studies with Escherichia coli guide work with Mycobacterium tuberculosis. Lethal action, which is mechanistically distinct from blocking growth, is associated with a surge in respiration and reactive oxygen species (ROS). Mutations in carbohydrate metabolism that attenuate ROS accumulation create pan-tolerance to antimicrobials, disinfectants, and environmental stressors. These observations indicate the existence of a general death pathway with respect to stressors. M. tuberculosis displays a variation on the death pathway idea, as stress-induced ROS is generated by NADH-mediated reductive stress rather than by respiration. A third approach, which emerges from lethality studies, uses a small molecule, N-acetyl cysteine, to artificially increase respiration and additional ROS accumulation. That enhances moxifloxacin lethality with M. tuberculosis in culture, during infection of cultured macrophages, and with infection of mice. Addition of ROS stimulators to fluoroquinolone treatment of tuberculosis constitutes a new direction for suppressing the transition of MDR-TB to XDR-TB.

MGIT sensitivity testing and genotyping of drug resistant Mycobacterium tuberculosis isolates from Mizoram, Northeast India

PURPOSE

Tuberculosis, a crucial infectious disease is still a health concern globally. India is among the countries with high MDR-TB burden. Currently, sputum smear microscopy using Ziehl Neelsen stain and GeneXpert are the only diagnostic means in Mizoram. This study was done to characterize local tuberculosis strains circulating in Mizoram.

METHODS

Sputum was cultured using MGIT 960 and DST was performed for Streptomycin, Isoniazid, Rifampicin, Ethambutol and Pyrazinamide. GeneXpert test was done simultaneously. DNA was extracted using Trueprep AUTO v2, molbio diagnostics. Antibiotic Resistance Genes and LSP were amplified and sequenced.

RESULTS

Ser315Thr was the most common mutation in katG among MDR-TB isolates. GeneXpert probes A and D drop out upon sequencing showed L511P, H526Q and H526L mutation. The L511P and H526Q mutations were seen in new and treated cases. Discrepancy between MGIT 960 and GeneXpert were observed. LSP-PCR revealed that Indo-Oceanic, East-African Indian, Euro-American and Beijing lineages were found in Mizoram.

CONCLUSION

This study provides mutation information on the resistant genotypes detected with GeneXpert as well as MGIT 960. It also provides information on the lineages of Mycobacterium tuberculosis circulating in the state. Utilization of sequencing technologies is essential in diagnostic laboratories to rule out discrepant results and as a cautionary measure to prevent wrong diagnosis and treatment.

Application of Next Generation Sequencing for Diagnosis and Clinical Management of Drug-Resistant Tuberculosis: Updates on Recent Developments in the Field

The World Health Organization’s End TB Strategy prioritizes universal access to an early diagnosis and comprehensive drug susceptibility testing (DST) for all individuals with tuberculosis (TB) as a key component of integrated, patient-centered TB care. Next generation whole genome sequencing (WGS) and its associated technology has demonstrated exceptional potential for reliable and comprehensive resistance prediction for Mycobacterium tuberculosis isolates, allowing for accurate clinical decisions. This review presents a descriptive analysis of research describing the potential of WGS to accelerate delivery of individualized care, recent advances in sputum-based WGS technology and the role of targeted sequencing for resistance detection. We provide an update on recent research describing the mechanisms of resistance to new and repurposed drugs and the dynamics of mixed infections and its potential implication on TB diagnosis and treatment. Whilst the studies reviewed here have greatly improved our understanding of recent advances in this arena, it highlights significant challenges that remain. The wide-spread introduction of new drugs in the absence of standardized DST has led to rapid emergence of drug resistance. This review highlights apparent gaps in our knowledge of the mechanisms contributing to resistance for these new drugs and challenges that limit the clinical utility of next generation sequencing techniques. It is recommended that a combination of genotypic and phenotypic techniques is warranted to monitor treatment response, curb emerging resistance and further dissemination of drug resistance.

Genotypic and phenotypic drug-resistance detection and prevalence of heteroresistance in patients with isoniazid- and multidrug-resistant tuberculosis in Ethiopia

Objective

To assess the agreement between genotypic and phenotypic methods for detecting drug resistance, and examine the prevalence of heteroresistance among isoniazid (INH)- and multidrug/rifampicin-resistant (MDR/RR) TB.

Method

In total, 127 Mycobacterium tuberculosis (Mtb) isolates, including 65 MDR/RR and 62 INH resistant, were used. First-line drug susceptibility testing (DST) was performed using the LJ method to determine the percentage of resistant bacteria. All drug-resistant isolates underwent testing with LPA. Heteroresistance was defined as simultaneous detection of wild-type and resistance-conferring mutations using LPA.

Result

The sensitivity of LPA (compared with LJ DST) was 96% for any INH-resistant TB and 94% for any RR TB. The prevalence of heteroresistance among the 123. Mtb isolates was 9.8%. The percentage of resistant bacteria ranged from 1% to 10% for heteroresistant TB. Rifampicin heteroresistance was detected in 1.6% of MDR TB patients. INH heteroresistance was detected in 1.6% and 16.7% of MDR and INH-resistant TB patients, respectively. The proportion of INH heteroresistance was significantly higher (p = 0.030) in persons living with HIV.

Conclusion

Some phenotypic drug resistances were not captured by LPA. The prevalence and percentage of resistant bacteria among heteroresistant TB highlight the importance of LPA for early detection of heteroresistant TB.

Detection of Mycobacterium tuberculosis multiple strains in sputum samples from patients with pulmonary tuberculosis in south western Uganda using MIRU-VNTR

Infections with multiple strains of Mycobacterium tuberculosis are now widely recognized as a common occurrence. Identification of patients infected with multiple strains provides both insight into the disease dynamics and the epidemiology of tuberculosis. Analysis of Mycobacterial Interspersed Repetitive Unit-Variable-Number Tandem Repeats (MIRU-VNTR) has been shown to be highly sensitive in detecting multiple M. tuberculosis strains even in sputum. The goal of this study was to identify cases of multiple M. tuberculosis strain infections among patients diagnosed with pulmonary tuberculosis in Southwestern Uganda and assessment of factors associated with multiple strain infections. DNA extracted directly from 78 sputum samples, each from an individual patient, was analyzed using the standard 24 loci MIRU-VNTR typing. Five (6.4%) of the 78 patients were infected with multiple strains of M. tuberculosis with all of them being the newly diagnosed cases while two-thirds of them were co-infected with HIV. Exact regression analysis projected that the natives were more likely to harbor multiple strains (OR; 0.981, 95% CI 0–7.926) as well as those with a high microbial load (OR; 0.390, 95% CI 0–3.8167). Despite these findings being not statistically significant due to the small sample size, this points to a critical component of disease dynamics that has clinical implications and emphasizes a need for a study using a larger cohort. It is also essential to study the potential factors associated with higher risk of exposure to newly diagnosed and HIV positive patients at the community level. In addition, our ability to detect multiple M. tuberculosis strains using the standard 24 loci MIRU-VNTR typing especially with allelic diversity in loci 2059 and 3171, which are excluded from the 15-locus MIRU-VNTR, lead us to recommend the use of this genotyping technique, especially in areas with tuberculosis endemicity similar to this study.

Transmission patterns of a Mycobacterium avium subsp. paratuberculosis clone within a single heard investigated by Whole Genome Sequencing

Mycobacterium avium subsp. paratuberculosis (MAP) is characterized by a low genomic rate of mutation. Current subtyping tools, such as Mini-Micro-satellite analyses, do to have not sufficient discriminatory power to disclose MAP's evolution on small spatial and temporal scales. The aim of the study was to investigate the population structure of MAP inside a single dairy herd using whole genome sequencing (WGS) approaches. For this purpose, the genomes of 43 field isolates, recovered from the faeces of 36 cows of the same dairy herd from 2012 to 2016, were sequenced by WGS. The isolates' genomes showed a low number (43) of polymorphic sites (SNPs), confirming the clonal origin of the herd infection. However, despite the limited genomic diversity found in WGS, the phylogenetic analysis was discriminatory enough to detect the presence of different genomic clades and sub-clades inside the herd population. In addition, the phylodynamic reconstruction showed the existence of an ancestor clade from which the other clades and sub-clades originated. Moreover, by reconstructing the putative within-herd transmission networks using WGS data, we demonstrated that: (i) in a herd where MAP is endemic, multiple isolates recovered from a single animal and differing from each other by few (three/four) SNPs can originate from different transmission or passive shedding events and not from intra-host evolution; and (ii) variability of minisatellites coupled with a few microsatellites does not represent reliable tracers of within-herd infection chains. Our findings show that WGS, coupled with relevant epidemiological information, represents a valuable tool to work out fine epidemiological and micro-evolutionary relationships such as those at herd-level scale.

Blood transcriptomics reveal the evolution and resolution of the immune response in tuberculosis

Blood transcriptomics have revealed major characteristics of the immune response in active TB, but the signature early after infection is unknown. In a unique clinically and temporally well-defined cohort of household contacts of active TB patients that progressed to TB, we define minimal changes in gene expression in incipient TB increasing in subclinical and clinical TB. While increasing with time, changes in gene expression were highest at 30 d before diagnosis, with heterogeneity in the response in household TB contacts and in a published cohort of TB progressors as they progressed to TB, at a bulk cohort level and in individual progressors. Blood signatures from patients before and during anti-TB treatment robustly monitored the treatment response distinguishing early and late responders. Blood transcriptomics thus reveal the evolution and resolution of the immune response in TB, which may help in clinical management of the disease.

High Rates of Recurrent Tuberculosis Disease: A Population-level Cohort Study

Background

Retreatment tuberculosis (TB) disease is common in high-prevalence settings. The risk of repeated episodes of recurrent TB is unknown. We calculated the rate of recurrent TB per subsequent episode by matching individual treatment episodes over a period of 13 years.

Methods

All recorded TB episodes in Cape Town between 2003 and 2016 were matched by probabilistic linkage of personal identifiers. Among individuals with a first episode notified in Cape Town and who completed their prior treatment successfully we estimated the recurrence rate stratified by subsequent episode and HIV status. We adjusted person-time to background mortality by age, sex, and HIV status.

Results

A total of 292 915 TB episodes among 263 848 individuals were included. The rate of recurrent TB was 16.4 per 1000 person-years (95% CI, 16.2–16.6), and increased per subsequent episode (8.4-fold increase, from 14.6 to 122.7 per 1000 from episode 2 to 6, respectively). These increases were similar stratified by HIV status. Rates among HIV positives were higher than among HIV negatives for episodes 2 and 3 (2- and 1.5-fold higher, respectively), and the same thereafter.

Conclusions

TB recurrence rates were high and increased per subsequent episode, independent of HIV status. This suggests that HIV infection is insufficient to explain the high burden of recurrence; it is more likely due to a high annual risk of infection combined with an increased risk of infection or progression to disease associated with a previous TB episode. The very high recurrence rates would justify increased TB surveillance of patients with >1 episode.

Social network analysis and whole-genome sequencing to evaluate disease transmission in a large, dynamic population: A study of avian mycobacteriosis in zoo birds

This study combined a social network analysis and whole-genome sequencing (WGS) to test for general patterns of contagious spread of a mycobacterial infection for which pathways of disease acquisition are not well understood. Our population included 275 cases diagnosed with avian mycobacteriosis that were nested in a source population of 16,430 birds at San Diego Zoo Wildlife Alliance facilities from 1992 through mid-2014. Mycobacteria species were determined using conventional methods and whole genome sequencing (WGS). Mycobacterium avium avium (MAA) and Mycobacterium genavense were the most common species of mycobacteria identified and were present in different proportions across bird taxa. A social network for the birds was constructed from the source population to identify directly and indirectly connected cases during time periods relevant to disease transmission. Associations between network connectivity and genetic similarity of mycobacteria (as determined by clusters of genotypes separated by few single nucleotide polymorphisms, or SNPs) were then evaluated in observed and randomly generated network permutations. Findings showed that some genotypes clustered along pathways of bird connectivity, while others were dispersed throughout the network. The proportion of directly connected birds having a similar mycobacterial genotype was 0.36 and significant (p<0.05). This proportion was higher (0.58) and significant for MAA but not for M. genavense. Evaluations of SNP distributions also showed genotypes of MAA were more related in connected birds than expected by chance; however, no significant patterns of genetic relatedness were identified for M. genavense, although data were sparse. Integrating the WGS analysis of mycobacteria with a social network analysis of their host birds revealed significant genetic clustering along pathways of connectivity, namely for MAA. These findings are consistent with a contagious process occurring in some, but not all, case clusters.

Screening for Mycobacterium tuberculosis Infection Using Beijing/K Strain-Specific Peptides in a School Outbreak Cohort

Background

The Beijing strain of Mycobacterium tuberculosis (M. tb) has been most frequently isolated from TB patients in South Korea, and the hyper-virulent Beijing/K genotype is associated with TB outbreaks. To examine the diagnostic potential of Beijing/K-specific peptides, we performed IFN-γ release assays (IGRA) using a MTBK antigen tube containing Beijing/K MTBK_24800, ESAT-6, and CFP-10 peptides in a cohort studied during a school TB outbreak.

Methods

A total of 758 contacts were investigated for M. tb infection, and 43 contacts with latent TB infection (LTBI) and 25 active TB patients were enrolled based on serial screening with QuantiFERON-TB Gold In-Tube tests followed by clinical examinations. Blood collected in MTBK antigen tubes was utilized for IGRA and multiplex cytokine bead arrays. Immune responses were retested in 24 patients after TB treatment, and disease progression was investigated in subjects with LTBI.

Results

Total proportions of active disease and LTBI during the outbreak were 3.7% (28/758) and 9.2% (70/758), respectively. All clinical isolates had a Beijing/K M. tb genotype. IFN-γ responses to the MTBK antigen identified M. tb infection and distinguished between active disease and LTBI. After anti-TB treatment, IFN-γ responses to the MTBK antigen were significantly reduced, and strong TNF-α responses at diagnosis were dramatically decreased.

Conclusions

MTBK antigen-specific IFN-γ has diagnostic potential for differentiating M. tb infection from healthy controls, and between active TB and LTBI as well. In addition, TNF-α is a promising marker for monitoring therapeutic responses. These data provide informative readouts for TB diagnostics and vaccine studies in regions where the Beijing/K strain is endemic.

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

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

Methods for Detecting Mycobacterial Mixed Strain Infections-A Systematic Review

Mixed strain infection (MSI) refers to the concurrent infection of a susceptible host with multiple strains of a single pathogenic species. Known to occur in humans and animals, MSIs deserve special consideration when studying transmission dynamics, evolution, and treatment of mycobacterial diseases, notably tuberculosis in humans and paratuberculosis (or Johne's disease) in ruminants. Therefore, a systematic review was conducted to examine how MSIs are defined in the literature, how widespread the phenomenon is across the host species spectrum, and to document common methods used to detect such infections. Our search strategy identified 121 articles reporting MSIs in both humans and animals, the majority (78.5%) of which involved members of the Mycobacterium tuberculosis complex, while only a few (21.5%) examined non-tuberculous mycobacteria (NTM). In addition, MSIs exist across various host species, but most reports focused on humans due to the extensive amount of work done on tuberculosis. We reviewed the strain typing methods that allowed for MSI detection and found a few that were commonly employed but were associated with specific challenges. Our review notes the need for standardization, as some highly discriminatory methods are not adapted to distinguish between microevolution of one strain and concurrent infection with multiple strains. Further research is also warranted to examine the prevalence of NTM MSIs in both humans and animals. In addition, it is envisioned that the accurate identification and a better understanding of the distribution of MSIs in the future will lead to important information on the epidemiology and pathophysiology of mycobacterial diseases.

Detection of Beijing strains of MDR M. tuberculosis and their association with drug resistance mutations in katG, rpoB, and embB genes

Background

Molecular epidemiological studies of Mycobacterium tuberculosis (MTB) are the core of current research to find out the association of the M. tuberculosis genotypes with its outbreak and transmission. The high prevalence of the Beijing genotype strain among multidrug resistance (MDR) TB has already been reported in various studies around India. The overall objective of this study was to detect the prevalence of Beijing genotype strains of MDR M. tuberculosis and their association with the clinical characteristics of TB patients.

Methods

In this study 381 M. tuberculosis clinical isolates were obtained from sputum samples from 2008 to 2014. The multiplex-PCR and Spoligotyping (n = 131) methods were used to investigate the prevalence of the Beijing genotype strain by targeting the Rv2820 gene and their association with drug resistance and clinical characteristics of TB patients. The drug susceptibility testing of first-line anti-TB drugs was performed by using the proportion method and MGIT960. A collection of isolates having Beijing and non-Beijing strains were also characterized to see if Beijing genotype strains had a higher rate of mutations at codons 516, 526 and 531 of the 81-bp region of the rpoB gene, codon 315 of the katG gene, and codon 306 of the embB gene.

Results

The sensitivities and specificities of multiplex-PCR assay compared to that of standard Spoligotyping was detected to be 100%. Further, we observe that the multi drug-resistance was significantly associated with Beijing genotype strains (p = 0.03) and a strong correlation between Beijing genotype strains and specific resistance mutations at the katG315, rpoB531, and embB306 codons (p = < 0.0001, < 0.0001 & 0.0014 respectively) was also found.

Conclusions

This rapid, simple, and cost-effective multiplex PCR assay can effectively be used for monitoring the prevalence of Beijing genotype strains in low resource settings. Findings of this study may provide a scientific basis for the development of new diagnostic tools for detection and effective management of DR-TB in countries with a higher incidence rate of Beijing genotype strains.

Investigating Non-sterilizing Cure in TB Patients at the End of Successful Anti-TB Therapy

Mycobacterium tuberculosis (Mtb) is extremely recalcitrant to antimicrobial chemotherapy requiring 6 months to treat drug-sensitive tuberculosis (TB). Despite this, 4-10% of cured patients will develop recurrent disease within 12 months after completing therapy. Reasons for relapse in cured TB patients remains speculative, attributed to both pathogen and host factors. Populations of dormant bacilli are hypothesized to cause relapse in initially cured TB patients however, development of tests to convincingly demonstrate their presence at the end of anti-TB treatment has been challenging. Previous studies have indicated the utility of culture filtrate supplemented media (CFSM) to detect differentially culturable tubercle bacilli (DCTB). Here, we show that 3/22 of clinically cured patients retained DCTB in induced sputum and bronchoalveolar lavage fluid (BALF), with one DCTB positive patient relapsing within the first year of completing therapy. We also show a correlation of DCTB status with "unresolved" end of treatment FDG PET-CT imaging. Additionally, 19 end of treatment induced sputum samples from patients not undergoing bronchoscopy were assessed for DCTB, identifying a further relapse case with DCTB. We further show that induced sputum is a less reliable source for the DCTB assay at the end of treatment, limiting the utility of this assay in a clinical setting. We next investigated the host proteome at the site of disease (BALF) using multiplexed proteomic analysis and compared these to active TB cases to identify host-specific factors indicative of cure. Distinct signatures stratified active from cured TB patients into distinct groups, with a DCTB positive, subsequently relapsing, end of treatment patient showing a proteomic signature closer to active TB disease than cure. This exploratory study offers evidence of live Mtb, undetectable with conventional culture methods, at the end of clinically successful treatment and putative host protein biomarkers of active disease and cure. These findings have implications for the assessment of true sterilizing cure in TB patients and opens new avenues for targeted approaches to monitor treatment response.

Towards standardisation: comparison of five whole genome sequencing (WGS) analysis pipelines for detection of epidemiologically linked tuberculosis cases

Background

Whole genome sequencing (WGS) is a reliable tool for studying tuberculosis (TB) transmission. WGS data are usually processed by custom-built analysis pipelines with little standardisation between them.

Aim

To compare the impact of variability of several WGS analysis pipelines used internationally to detect epidemiologically linked TB cases.

Methods

From the Netherlands, 535 Mycobacterium tuberculosis complex (MTBC) strains from 2016 were included. Epidemiological information obtained from municipal health services was available for all mycobacterial interspersed repeat unit-variable number of tandem repeat (MIRU-VNTR) clustered cases. WGS data was analysed using five different pipelines: one core genome multilocus sequence typing (cgMLST) approach and four single nucleotide polymorphism (SNP)-based pipelines developed in Oxford, United Kingdom; Borstel, Germany; Bilthoven, the Netherlands and Copenhagen, Denmark. WGS clusters were defined using a maximum pairwise distance of 12 SNPs/alleles.

Results

The cgMLST approach and Oxford pipeline clustered all epidemiologically linked cases, however, in the other three SNP-based pipelines one epidemiological link was missed due to insufficient coverage. In general, the genetic distances varied between pipelines, reflecting different clustering rates: the cgMLST approach clustered 92 cases, followed by 84, 83, 83 and 82 cases in the SNP-based pipelines from Copenhagen, Oxford, Borstel and Bilthoven respectively.

Conclusion

Concordance in ruling out epidemiological links was high between pipelines, which is an important step in the international validation of WGS data analysis. To increase accuracy in identifying TB transmission clusters, standardisation of crucial WGS criteria and creation of a reference database of representative MTBC sequences would be advisable.

Mycobacterium tuberculosis polyclonal infections through treatment and recurrence

Background

Mixed/polyclonal infections due to different genotypes are reported in Tuberculosis. The current study was designed to understand the fate of mixed infections during the course of treatment and follow-up and its role in disease pathogenesis.

Methods

Sputum samples were collected on 0,1,2,3,6,12 and 24 months from 157 treatment-naïve patients, cultures subjected to Drug-Susceptibility-testing (MGIT 960), spoligotyping, MIRU-VNTR and SNP genotyping. All isolated colonies on thin layer agar (7H11) were subjected to spoligotyping.

Findings

One thirty three baseline cultures were positive (133/157, 84.7%), 43(32.3%) had mixture of genotypes. Twenty-four of these patients (55.8%) showed change in genotype while six showed different drug-susceptibility patterns while on treatment. Twenty-three (53.5%) patients with polyclonal infections showed resistance to at least one drug compared to 10/90 (11.1%) monoclonal infections (P<0.0001). Eight patients had recurrent TB, two with a new genotype and two with altered phenotypic DST.

Conclusions

The coexistence of different genotypes and change of genotypes during the same disease episode, while on treatment, confirms constancy of polyclonal infections. The composition of the mixture of genotypes and the relative predominance may be missed by culture due to its limit of detection. Polyclonal infections in TB could be a rule rather than exception and challenges the age-old dogma of reactivation/reinfection.

Inoculum size and traits of the infecting clinical strain define the protection level against Mycobacterium tuberculosis infection in a rabbit model

Host protective immunity against pathogenic Mycobacterium tuberculosis (Mtb) infection is variable and poorly understood. Both prior Mtb infection and BCG vaccination have been reported to confer some protection against subsequent infection and/or disease. However, the immune correlates of host protection with or without BCG vaccination remain poorly understood. Similarly, the host response to concomitant infection with mixed Mtb strains is unclear. In this study, we used the rabbit model to examine the host response to various infectious doses of virulent Mtb HN878 with and without prior BCG vaccination, as well as simultaneous infection with a mixture of two Mtb clinical isolates HN878 and CDC1551. We demonstrate that both the ability of host immunity to control pulmonary Mtb infection and the protective efficacy of BCG vaccination against the progression of Mtb infection to disease is dependent on the infectious inoculum. The host response to infection with mixed Mtb strains mirrors the differential responses seen during infection with each of the strains alone. The protective response mounted against a hyperimmunogenic Mtb strain has a limited impact on the control of disease caused by a hypervirulent strain. This preclinical study will aid in predicting the success of any vaccination strategy and in optimizing TB vaccines.

QuantTB - a method to classify mixed Mycobacterium tuberculosis infections within whole genome sequencing data

BACKGROUND

Mixed infections of Mycobacterium tuberculosis and antibiotic heteroresistance continue to complicate tuberculosis (TB) diagnosis and treatment. Detection of mixed infections has been limited to molecular genotyping techniques, which lack the sensitivity and resolution to accurately estimate the multiplicity of TB infections. In contrast, whole genome sequencing offers sensitive views of the genetic differences between strains of M. tuberculosis within a sample. Although metagenomic tools exist to classify strains in a metagenomic sample, most tools have been developed for more divergent species, and therefore cannot provide the sensitivity required to disentangle strains within closely related bacterial species such as M. tuberculosis. Here we present QuantTB, a method to identify and quantify individual M. tuberculosis strains in whole genome sequencing data. QuantTB uses SNP markers to determine the combination of strains that best explain the allelic variation observed in a sample. QuantTB outputs a list of identified strains, their corresponding relative abundances, and a list of drugs for which resistance-conferring mutations (or heteroresistance) have been predicted within the sample.

RESULTS

We show that QuantTB has a high degree of resolution and is capable of differentiating communities differing by less than 25 SNPs and identifying strains down to 1× coverage. Using simulated data, we found QuantTB outperformed other metagenomic strain identification tools at detecting strains and quantifying strain multiplicity. In a real-world scenario, using a dataset of 50 paired clinical isolates from a study of patients with either reinfections or relapses, we found that QuantTB could detect mixed infections and reinfections at rates concordant with a manually curated approach.

CONCLUSION

QuantTB can determine infection multiplicity, identify hetero-resistance patterns, enable differentiation between relapse and re-infection, and clarify transmission events across seemingly unrelated patients - even in low-coverage (1×) samples. QuantTB outperforms existing tools and promises to serve as a valuable resource for both clinicians and researchers working with clinical TB samples.

Usefulness of 3'- 5' IS6110-RFLP genotyping and spoligotyping of Mycobacterium tuberculosis isolated in a tertiary hospital: a retrospective study detecting unsuspected epidemiological events

A drug resistance survey involving Mycobacterium tuberculosis isolated from patients of a tertiary Hospital in the Rio de Janeiro city (RJ), Brazil, between the years 1996 and 1998 revealed a high frequency of isoniazid (H^R) resistance. These isolates were revisited and genotyped. Patients came from different RJ neighborhoods and municipalities, and 70% were outpatients. Applying the 3’ and 5’ IS 6110 -RFLP and the Spoligotype genotyping methods, the clonal structure of this population was investigated obtaining a snapshot of past epidemiological events. The 3’ clusters were subsequently 5’ IS 6110 -RFLP typed. Spoligotyping was analyzed in the SITVIT2 database. Epidemiological relationships were investigated. The major lineage was T (54.4%), and SIT 53/T1 and SIT 535/T1 were the most frequent. The T1 sublineage comprises 12.8% of resistant strains and SIT 535 were assigned for 31.8% of them. Orphan patterns corresponded to 12% and 73.3% and belonged to the T lineage. One pattern was unlisted in the SITVIT2. The 5’ IS 6110 -RFLP did not confirm 3/12 of the 3’ IS 6110 -RFLP clusters. A combination of all methods decreased the number of clusters to three. Nosocomial transmission was associated with one cluster involving a hospital cupbearer. This event was suspected in a multidrug resistant-TB inpatient caregiver who harbored a mixed infection. The 3’ IS 6110 clusters were associated with H^R (p=0.046). These genotypic retrospective data may reflect a fraction of more extensive recent transmission in different communities that may be corroborated by the concentration of H^R patients, and may serve as a database for further evolutionary and characterization evaluation of circulating strains and together with epidemiological data favors a more effective transmission control.

Effect of mixed infection on TB dynamics

Poor living conditions, overcrowding and strain diversity are some of the factors that influence mixed infection in Tuberculosis (TB) at the population level. We formulate a mathematical model for mixed infection in TB using nonlinear ordinary differential equations where such factors were represented as probabilities of acquiring mixed infection. A qualitative analysis of the model shows that it exhibits multiple endemic equilibria and backward bifurcation for certain parameter values. The reactivation rate and transmission rate of individuals with mixed infection were of importance as well as the probabilities for latent individuals to acquire mixed infection. We calculate the prevalence of mixed infection from the model and the effect of mixed infection on TB incidence, TB prevalence and Mycobacterium tuberculosis (MTB) infection rate. Numerical simulations show that mixed infection may explain high TB incidences in areas which have a high strain diversity, poor living conditions and are overcrowded even without HIV.

Bacterial Genotyping of Central Nervous System Tuberculosis in South Africa: Heterogenic Mycobacterium tuberculosis Infection and Predominance of Lineage 4

Tuberculous meningitis (TBM), the most severe extrapulmonary manifestation of tuberculosis, is caused by the pathogen Mycobacterium tuberculosis The M. tuberculosis complex includes seven lineages, all described to harbor a unique geographical dissemination pattern and clinical presentation. In this study, we set out to determine whether a certain M. tuberculosis lineage demonstrated tropism to cause TBM in patients from Cape Town, South Africa. DNA was extracted from formalin-fixed paraffin-embedded central nervous system (CNS) tissue from a unique neuropathological cohort of 83 TBM patients, collected between 1975 and 2012. M. tuberculosis lineages 1, 2, 3, and 4 were determined using an allele-specific PCR and Sanger sequencing. Of the 83 patient specimens tested, bacterial characterization could be performed on 46 specimens (55%). M. tuberculosis lineage 4 was present in 26 patient specimens (56%), and non-lineage 4 was identified in 10 cases (22%). Moreover, genomic heterogeneity was detected in the CNS specimens of 7 adults and 3 children. We could show that infection of the CNS is not restricted to a single M. tuberculosis lineage and that even young children with rapid progression of disease can harbor more than one M. tuberculosis lineage in the CNS.

Deciphering drug resistance in Mycobacterium tuberculosis using whole-genome sequencing: progress, promise, and challenges

Tuberculosis (TB) is a global infectious threat that is intensified by an increasing incidence of highly drug-resistant disease. Whole-genome sequencing (WGS) studies of Mycobacterium tuberculosis, the causative agent of TB, have greatly increased our understanding of this pathogen. Since the first M. tuberculosis genome was published in 1998, WGS has provided a more complete account of the genomic features that cause resistance in populations of M. tuberculosis, has helped to fill gaps in our knowledge of how both classical and new antitubercular drugs work, and has identified specific mutations that allow M. tuberculosis to escape the effects of these drugs. WGS studies have also revealed how resistance evolves both within an individual patient and within patient populations, including the important roles of de novo acquisition of resistance and clonal spread. These findings have informed decisions about which drug-resistance mutations should be included on extended diagnostic panels. From its origins as a basic science technique, WGS of M. tuberculosis is becoming part of the modern clinical microbiology laboratory, promising rapid and improved detection of drug resistance, and detailed and real-time epidemiology of TB outbreaks. We review the successes and highlight the challenges that remain in applying WGS to improve the control of drug-resistant TB through monitoring its evolution and spread, and to inform more rapid and effective diagnostic and therapeutic strategies.

Mycobacterium tuberculosis Next-Generation Whole Genome Sequencing: Opportunities and Challenges

Mycobacterium tuberculosis drug resistance is a threat to global tuberculosis (TB) control. Comprehensive and timely drug susceptibility determination is critical to inform appropriate treatment of drug-resistant tuberculosis (DR-TB). Phenotypic drug susceptibility testing (DST) is the gold standard for M. tuberculosis drug resistance determination. M. tuberculosis whole genome sequencing (WGS) has the potential to be a one-stop method for both comprehensive DST and epidemiological investigations. We discuss in this review the tremendous opportunities that next-generation WGS presents in terms of understanding the molecular epidemiology of tuberculosis and mechanisms of drug resistance. The potential clinical value and public health impact in the areas of DST for patient management and tracing of transmission chains for timely public health intervention are also discussed. We present the current challenges for the implementation of WGS in low and middle-income settings. WGS analysis has already been adapted routinely in laboratories to inform patient management and public health interventions in low burden high-income settings such as the United Kingdom. We predict that the technology will be adapted similarly in high burden settings where the impact on the epidemic will be greatest.

Mixed Mycobacterium tuberculosis-Strain Infections Are Associated With Poor Treatment Outcomes Among Patients With Newly Diagnosed Tuberculosis, Independent of Pretreatment Heteroresistance

Background

Heteroresistant Mycobacterium tuberculosis infections (defined as concomitant infection with drug-resistant and drug-susceptible strains) may explain the higher risk of poor tuberculosis treatment outcomes observed among patients with mixed-strain M. tuberculosis infections. We investigated the clinical effect of mixed-strain infections while controlling for pretreatment heteroresistance in a population-based sample of patients with tuberculosis starting first-line tuberculosis therapy in Botswana.

Methods

We performed 24-locus mycobacterial interspersed repetitive unit-variable number tandem-repeat analysis and targeted deep sequencing on baseline primary cultured isolates to detect mixed infections and heteroresistance, respectively. Drug-sensitive, micro-heteroresistant, macro-heteroresistant, and fixed-resistant infections were defined as infections in which the frequency of resistance was <0.1%, 0.1%-4%, 5%-94%, and ≥95%, respectively, in resistance-conferring domains of the inhA promoter, the katG gene, and the rpoB gene.

Results

Of the 260 patients with tuberculosis included in the study, 25 (9.6%) had mixed infections and 30 (11.5%) had poor treatment outcomes. Micro-heteroresistance, macro-heteroresistance, and fixed resistance were found among 11 (4.2%), 2 (0.8%), and 11 (4.2%), respectively, for isoniazid and 21 (8.1%), 0 (0%), and 10 (3.8%), respectively, for rifampicin. In multivariable analysis, mixed infections but not heteroresistant infections independently predicted poor treatment outcomes.

Conclusions

Among patients starting first-line tuberculosis therapy in Botswana, mixed infections were associated with poor tuberculosis treatment outcomes, independent of heteroresistance.

The value of transcriptomics in advancing knowledge of the immune response and diagnosis in tuberculosis

Blood transcriptomics analysis of tuberculosis has revealed an interferon-inducible gene signature that diminishes in expression after successful treatment; this promises improved diagnostics and treatment monitoring, which are essential for the eradication of tuberculosis. Sensitive radiography revealing lung abnormalities and blood transcriptomics have demonstrated heterogeneity in patients with active tuberculosis and exposed asymptomatic people with latent tuberculosis, suggestive of a continuum of infection and immune states. Here we describe the immune response to infection with Mycobacterium tuberculosis revealed through the use of transcriptomics, as well as differences among clinical phenotypes of infection that might provide information on temporal changes in host immunity associated with evolving infection. We also review the diverse blood transcriptional signatures, composed of small sets of genes, that have been proposed for the diagnosis of tuberculosis and the identification of at-risk asymptomatic people and suggest novel approaches for the development of such biomarkers for clinical use.

The Evolution of Genotyping Strategies To Detect, Analyze, and Control Transmission of Tuberculosis

The introduction of genotypic tools to analyze Mycobacterium tuberculosis isolates has transformed our knowledge of the transmission dynamics of this pathogen. We discuss the development of the laboratory methods that have been applied in recent years to study the epidemiology of M. tuberculosis. This review integrates two approaches: on the one hand, it considers how genotyping techniques have evolved over the years; and on the other, it looks at how the way we think these techniques should be applied has changed. We begin by examining the application of fingerprinting tools to suspected outbreaks only, before moving on to universal genotyping schemes, and finally we describe the latest real-time strategies used in molecular epidemiology. We also analyze refined approaches to obtaining epidemiological data from patients and to increasing the discriminatory power of genotyping by techniques based on genomic characterization. Finally, we review the development of integrative solutions to reconcile the speed of PCR-based methods with the high discriminatory power of whole-genome sequencing in easily implemented formats adapted to low-resource settings. Our analysis of future considerations highlights the need to bring together the three key elements of high-quality surveillance of transmission in tuberculosis, namely, speed, precision, and ease of implementation.

Identifying mixed Mycobacterium tuberculosis infections from whole genome sequence data

BACKGROUND

Mixed, polyclonal Mycobacterium tuberculosis infection occurs in natural populations. Developing an effective method for detecting such cases is important in measuring the success of treatment and reconstruction of transmission between patients. Using whole genome sequence (WGS) data, we assess two methods for detecting mixed infection: (i) a combination of the number of heterozygous sites and the proportion of heterozygous sites to total SNPs, and (ii) Bayesian model-based clustering of allele frequencies from sequencing reads at heterozygous sites.

RESULTS

In silico and in vitro artificially mixed and known pure M. tuberculosis samples were analysed to determine the specificity and sensitivity of each method. We found that both approaches were effective in distinguishing between pure strains and mixed infection where there was relatively high (> 10%) proportion of a minor strain in the mixture. A large dataset of clinical isolates (n = 1963) from the Karonga Prevention Study in Northern Malawi was tested to examine correlations with patient characteristics and outcomes with mixed infection. The frequency of mixed infection in the population was found to be around 10%, with an association with year of diagnosis, but no association with age, sex, HIV status or previous tuberculosis.

CONCLUSIONS

Mixed Mycobacterium tuberculosis infection was identified in silico using whole genome sequence data. The methods presented here can be applied to population-wide analyses of tuberculosis to estimate the frequency of mixed infection, and to identify individual cases of mixed infections. These cases are important when considering the evolution and transmission of the disease, and in patient treatment.

Geospatial distribution of Mycobacterium tuberculosis genotypes in Africa

Objective

To investigate the distribution of Mycobacterium tuberculosis genotypes across Africa.

Methods

The SITVIT2 global repository and PUBMED were searched for spoligotype and published genotype data respectively, of M. tuberculosis from Africa. M. tuberculosis lineages in Africa were described and compared across regions and with those from 7 European and 6 South-Asian countries. Further analysis of the major lineages and sub-lineages using Principal Component analysis (PCA) and hierarchical cluster analysis were done to describe clustering by geographical regions. Evolutionary relationships were assessed using phylogenetic tree analysis.

Results

A total of 14727 isolates from 35 African countries were included in the analysis and of these 13607 were assigned to one of 10 major lineages, whilst 1120 were unknown. There were differences in geographical distribution of major lineages and their sub-lineages with regional clustering. Southern African countries were grouped based on high prevalence of LAM11-ZWE strains; strains which have an origin in Portugal. The grouping of North African countries was due to the high percentage of LAM9 strains, which have an origin in the Eastern Mediterranean region. East African countries were grouped based on Central Asian (CAS) and East-African Indian (EAI) strain lineage possibly reflecting historic sea trade with Asia, while West African Countries were grouped based on Cameroon lineage of unknown origin. A high percentage of the Haarlem lineage isolates were observed in the Central African Republic, Guinea, Gambia and Tunisia, however, a mixed distribution prevented close clustering.

Conclusions

This study highlighted that the TB epidemic in Africa is driven by regional epidemics characterized by genetically distinct lineages of M. tuberculosis. M. tuberculosis in these regions may have been introduced from either Europe or Asia and has spread through pastoralism, mining and war. The vast array of genotypes and their associated phenotypes should be considered when designing future vaccines, diagnostics and anti-TB drugs.

Backward bifurcation and hysteresis in models of recurrent tuberculosis

An epidemiological model is presented that provides a comprehensive description of the transmission pathways involved for recurrent tuberculosis (TB), whereby cured individuals can become reinfected. Our main goal is to determine conditions that lead to the appearance of a backward bifurcation. This occurs when an asymptotically stable infection free equilibrium concurrently exists with a stable non-trivial equilibria even though the basic reproduction number R₀ is less than unity. Although, some 10-30% cases of TB are recurrent, the role of recurrent TB as far as the formation of backward bifurcation is concerned, is rarely if ever studied. The model used here incorporates progressive primary infection, exogenous reinfection, endogenous reactivation and recurrent TB as transmission mechanisms that contribute to TB progression. Unlike other studies of TB dynamics that make use of frequency dependent transmission rates, our analysis provides exact backward bifurcation threshold conditions without resorting to commonly applied approximations and simplifying assumptions. Exploration of the model through analytical and numerical analysis reveal that recurrent TB is sometimes capable of triggering hysteresis effects which allow TB to persist when R₀ < 1 even though there is no backward bifurcation. Furthermore, recurrent TB can independently induce backward bifurcation phenomena if it exceeds a certain threshold.

Mixed infections in tuberculosis: The missing part in a puzzle

The mixed strains infection phenomenon is a major problem posing serious challenges in control of tuberculosis (TB). In patients with mixed infection, several different strains of Mycobacterium tuberculosis can be isolated simultaneously. Although different genotyping methods and various molecular approaches can be employed for detection of mixed infection in clinical samples, the MIRU-VNTR technique is more sensitive with higher discriminative power than many widely used techniques. Furthermore, the recent introduction of whole genome sequencing (WGS) promises to reveal more details about mixed infection with high resolution. WGS has been used for detection of mixed infection with high sensitivity and discriminatory, but the technology is currently limited to developed countries. Mixed infection may involve strains with different susceptibility patterns, which may alter the treatment outcome. In this report, we review the current concepts of mixed strains infection and also infection involving strains with a different susceptibility pattern in TB. We evaluate the importance of identifying mixed infection for diagnosis as well as treatment and highlight the accuracy and clinical utility of direct genotyping of clinical specimens.

Prediction of Local Transmission of Mycobacterium tuberculosis Isolates of a Predominantly Beijing Lineage by Use of a Variable-Number Tandem-Repeat Typing Method Incorporating a Consensus Set of Hypervariable Loci

Strain genotyping based on the variable-number tandem repeat (VNTR) is widely applied for identifying the transmission of Mycobacterium tuberculosis A consensus set of four hypervariable loci (1982, 3232, 3820, and 4120) has been proposed to improve the discrimination of Beijing lineage strains. Herein, we evaluated the utility of these four hypervariable loci for tracing local tuberculosis transmission in 981 cases over a 14-month period in Japan (2010 to 2011). We used six different VNTR systems, with or without the four hypervariable loci. Patient ages and weighted standard distances (a measure of the dispersion of genotype-clustered cases) were used as proxies for estimating local tuberculosis transmission. The highest levels of isolate discrimination were achieved with VNTR systems that incorporated the four hypervariable loci (i.e., the Japan Anti-Tuberculosis Association [JATA]18-VNTR, mycobacterial interspersed repetitive unit [MIRU]28-VNTR, and 24_Beijing-VNTR). The clustering rates by JATA12-VNTR, MIRU15-VNTR, JATA15-VNTR, JATA18-VNTR, MIRU28-VNTR, and 24_Beijing-VNTR systems were 52.2%, 51.0%, 39.0%, 24.1%, 23.1%, and 22.0%, respectively. As the discriminative power increased, the median weighted standard distances of the clusters tended to decrease (from 311 to 80 km, P < 0.001, Jonckheere-Terpstra trend test). Concurrently, the median ages of patients in the clusters tended to decrease (from 68 to 60 years, P < 0.001, Jonckheere-Terpstra trend test). These findings suggest that strain typing using the four hypervariable loci improves the prediction of active local tuberculosis transmission. The four-locus set can therefore contribute to the targeted control of tuberculosis in settings with high prevalence of Beijing lineage strains.

Evaluation of spoligotyping, SNPs and customised MIRU-VNTR combination for genotyping Mycobacterium tuberculosis clinical isolates in Madagascar

BACKGROUND

Combining different molecular typing methods for Mycobacterium tuberculosis complex (MTBC) can be a powerful tool for molecular epidemiology-based investigation of TB. However, the current standard method that provides high discriminatory power for such a combination, mycobacterial interspersed repetitive units-variable numbers of tandem repeats typing (MIRU-VNTR), is laborious, time-consuming and often too costly for many resource-limited laboratories. We aimed to evaluate a reduced set of loci for MIRU-VNTR typing in combination with spoligotyping and SNP-typing for routine molecular epidemiology of TB.

METHOD

Spoligotyping and SNP-typing, in combination with the 15 loci MIRU-VNTR typing, were first used to type clinical MTBC isolates (n = 158) from Madagascar. A step by step reduction of MIRU-VNTR loci number was then performed according to the Hunter and Gaston Discriminatory Index (HGDI) and to the Principal component analysis (PCA) correlation with the spoligotype profiles to evaluate the discrimination power inside the generated spoligotype clusters. The 15 MIRU-VNTR was used as reference and SNP-typing was used to determine the main MTBC lineages.

RESULTS

Of the 158 clinical isolates studied, the SNP-typing classified 23 into Lineage 1 (14.6%), 31 into Lineage 2 (19.6%), 23 into Lineage 3 (14.6%) and 81 into Lineage 4 strains (51.3%). 37 different spoligotypes profiles were obtained, 15 of which were unique and 20 in clusters. 15-loci MIRU-VNTR typing revealed 144 different genotypes: 132 isolates had a unique MIRU-VNTR profile and 27 isolates were grouped into 12 clusters. After a stepwise reduction of the MIRU-VNTR loci number within each main spoligotype families, three different sets composed of 5 customised MIRU-VNTR loci had a similar discrimination level to the reference 15 loci MIRU-VNTR in lineage 1, lineage 2 and lineage 3. For lineage 4, a set of 4 and 3 MIRU-VNTR loci were proposed to subtype the Harleem and LAM spoligotype families, respectively. For the T spoligotype family, a set of 5 MIRU-VNTR loci was proposed.

CONCLUSION

According to the lineages and the spoligotype families, the number of MIRU-VNTR loci can be reduced to get an optimal classification of MTBC. These customized sets of MIRU-VNTR loci reduce workload and save resources while maintaining optimal discriminatory power.

Mycobacterium tuberculosis Whole Genome Sequences From Southern India Suggest Novel Resistance Mechanisms and the Need for Region-Specific Diagnostics

BACKGROUND.

India is home to 25% of all tuberculosis cases and the second highest number of multidrug resistant cases worldwide. However, little is known about the genetic diversity and resistance determinants of Indian Mycobacterium tuberculosis, particularly for the primary lineages found in India, lineages 1 and 3.

METHODS.

We whole genome sequenced 223 randomly selected M. tuberculosis strains from 196 patients within the Tiruvallur and Madurai districts of Tamil Nadu in Southern India. Using comparative genomics, we examined genetic diversity, transmission patterns, and evolution of resistance.

RESULTS.

Genomic analyses revealed (11) prevalence of strains from lineages 1 and 3, (11) recent transmission of strains among patients from the same treatment centers, (11) emergence of drug resistance within patients over time, (11) resistance gained in an order typical of strains from different lineages and geographies, (11) underperformance of known resistance-conferring mutations to explain phenotypic resistance in Indian strains relative to studies focused on other geographies, and (11) the possibility that resistance arose through mutations not previously implicated in resistance, or through infections with multiple strains that confound genotype-based prediction of resistance.

CONCLUSIONS.

In addition to substantially expanding the genomic perspectives of lineages 1 and 3, sequencing and analysis of M. tuberculosis whole genomes from Southern India highlight challenges of infection control and rapid diagnosis of resistant tuberculosis using current technologies. Further studies are needed to fully explore the complement of diversity and resistance determinants within endemic M. tuberculosis populations.

Relapse, re-infection and mixed infections in tuberculosis disease

Tuberculosis (TB) disease can be characterized by genotypic and phenotypic complexity in Mycobacterium tuberculosis bacilli within a single patient. This microbiological heterogeneity has become an area of intense study due its perceived importance in drug tolerance, drug resistance and as a surrogate measure of transmission rates. This review presents a descriptive analysis of research describing the prevalence of mixed-strain TB infections in geographically distinct locations. Despite significant variation in disease burden and a rampant human immunodeficiency virus (HIV)-TB co-epidemic, there was no difference in the prevalence range of mixed infections reported in African countries when compared to the rest of the world. The occurrence of recurrent TB was associated with a higher prevalence of mixed-strain infections, but this difference was not reported as statistically significant. These interpretations were limited by differences in the design and overall size of the studies assessed. Factors such as sputum quality, culture media, number of repeated culture steps, molecular typing methods and HIV-infection status can affect the detection of mixed-strain infection. It is recommended that future clinical studies should focus on settings with varying TB burdens, with a common sample processing protocol to gain further insight into these phenomena and develop novel transmission blocking strategies.

Rapid and simple identification of Beijing genotype strain of Mycobacterium tuberculosis using a loop-mediated isothermal amplification assay

Beijing genotype strains of Mycobacterium tuberculosis are geographically widespread and pose a notorious public health problem, these strains causing outbreaks of multidrug-resistant tuberculosis (TB); some studies have reported an association with drug resistance. Because the prevalence of Beijing strain has a substantial impact on TB control programs, the availability of a rapid and reliable method for detecting these strains is important for epidemiological monitoring of their circulation. The main methods currently used to identify Beijing genotype strains are IS6110 DNA fingerprinting, spoligotyping and PCR to detect specific deletions such as region of difference (RD)207. More recently, multiplex PCR assay using a Beijing-specific single nucleotide polymorphism (SNP) has been developed for detecting Beijing lineage strains. However, these methods are time-consuming and technically demanding. In the present study, a loop-mediated isothermal amplification (LAMP) assay that allows specific identification of Beijing genotype strain was developed. This Beijing genotype strain-identifying LAMP assay was performed 214 clinical isolates and the results compared with those of conventional PCR that targeted RD207 and Rv0679c-targreting multiplex PCR for Beijing lineage identification. LAMP assay showed 100% sensitivity and specificity compared with RD207-PCR. Furthermore, the sensitivity and specificity were 99.3% and 100%, respectively, compared with Rv0679c-multiplex PCR. This LAMP assay could be used routinely in local laboratories to monitor the prevalence of the Beijing genotype strain and thereby used to help control the spread of these potentially highly virulent and drug resistant strains.

Whole-genome analysis of mycobacteria from birds at the San Diego Zoo

Methods

Mycobacteria isolated from more than 100 birds diagnosed with avian mycobacteriosis at the San Diego Zoo and its Safari Park were cultured postmortem and had their whole genomes sequenced. Computational workflows were developed and applied to identify the mycobacterial species in each DNA sample, to find single-nucleotide polymorphisms (SNPs) between samples of the same species, to further differentiate SNPs between as many as three different genotypes within a single sample, and to identify which samples are closely clustered genomically.

Results

Nine species of mycobacteria were found in 123 samples from 105 birds. The most common species were Mycobacterium avium and Mycobacterium genavense, which were in 49 and 48 birds, respectively. Most birds contained only a single mycobacterial species, but two birds contained a mixture of two species. The M. avium samples represent diverse strains of M. avium avium and M. avium hominissuis, with many pairs of samples differing by hundreds or thousands of SNPs across their common genome. By contrast, the M. genavense samples are much closer genomically; samples from 46 of 48 birds differ from each other by less than 110 SNPs. Some birds contained two, three, or even four genotypes of the same bacterial species. Such infections were found in 4 of 49 birds (8%) with M. avium and in 11 of 48 birds (23%) with M. genavense. Most were mixed infections, in which the bird was infected by multiple mycobacterial strains, but three infections with two genotypes differing by ≤ 10 SNPs were likely the result of within-host evolution. The samples from 31 birds with M. avium can be grouped into nine clusters within which any sample is ≤ 12 SNPs from at least one other sample in the cluster. Similarly, the samples from 40 birds with M. genavense can be grouped into ten such clusters. Information about these genomic clusters is being used in an ongoing, companion study of mycobacterial transmission to help inform management of bird collections.

Clonal Complexity in Mycobacterium tuberculosis Can Hamper Diagnostic Procedures

Clonal complexity is increasingly accepted in Mycobacterium tuberculosis infection, including mixed infections by ≥2 strains, which usually occur in settings with a high burden of tuberculosis and/or a high risk of overexposure to infected patients. Mixed infections can hamper diagnostic procedures; obtaining an accurate antibiogram is difficult when the susceptibility patterns of the strains differ. Here, we show how mixed infections can also prove challenging for other diagnostic procedures, even outside settings where mixed infections are traditionally expected. We show how an unnoticed mixed infection in an HIV-positive patient diagnosed in Madrid, Spain, with differences in the representativeness of the coinfecting strains in different sputum samples, markedly complicated the resolution of a laboratory cross-contamination false positivity alert.