Identifying Mixed Mycobacterium tuberculosis Infection and Laboratory Cross-Contamination during Mycobacterial Sequencing Programs

Application of Metagenomic Next-Generation Sequencing in Mycobacterium tuberculosis Infection

The fight against Mycobacterium tuberculosis (MTB) has been going on for thousands of years, while it still poses a threat to human health. In addition to routine detections, metagenomic next-generation sequencing (mNGS) has begun to show presence as a comprehensive and hypothesis-free test. It can not only detect MTB without isolating specific pathogens but also suggest the co-infection pathogens or underlying tumor simultaneously, which is of benefit to assist in comprehensive clinical diagnosis. It also shows the potential to detect multiple drug resistance sites for precise treatment. However, considering the cost performance compared with conventional assays (especially Xpert MTB/RIF), mNGS seems to be overqualified for patients with mild and typical symptoms. Technology optimization of sequencing and analyzing should be conducted to improve the positive rate and broaden the applicable fields.

Population structure, biogeography and transmissibility of Mycobacterium tuberculosis

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

Rapid and accurate SNP genotyping of clonal bacterial pathogens with BioHansel

Hierarchical genotyping approaches can provide insights into the source, geography and temporal distribution of bacterial pathogens. Multiple hierarchical SNP genotyping schemes have previously been developed so that new isolates can rapidly be placed within pre-computed population structures, without the need to rebuild phylogenetic trees for the entire dataset. This classification approach has, however, seen limited uptake in routine public health settings due to analytical complexity and the lack of standardized tools that provide clear and easy ways to interpret results. The BioHansel tool was developed to provide an organism-agnostic tool for hierarchical SNP-based genotyping. The tool identifies split k-mers that distinguish predefined lineages in whole genome sequencing (WGS) data using SNP-based genotyping schemes. BioHansel uses the Aho-Corasick algorithm to type isolates from assembled genomes or raw read sequence data in a matter of seconds, with limited computational resources. This makes BioHansel ideal for use by public health agencies that rely on WGS methods for surveillance of bacterial pathogens. Genotyping results are evaluated using a quality assurance module which identifies problematic samples, such as low-quality or contaminated datasets. Using existing hierarchical SNP schemes for Mycobacterium tuberculosis and Salmonella Typhi, we compare the genotyping results obtained with the k-mer-based tools BioHansel and SKA, with those of the organism-specific tools TBProfiler and genotyphi, which use gold-standard reference-mapping approaches. We show that the genotyping results are fully concordant across these different methods, and that the k-mer-based tools are significantly faster. We also test the ability of the BioHansel quality assurance module to detect intra-lineage contamination and demonstrate that it is effective, even in populations with low genetic diversity. We demonstrate the scalability of the tool using a dataset of ~8100 S. Typhi public genomes and provide the aggregated results of geographical distributions as part of the tool’s output. BioHansel is an open source Python 3 application available on PyPI and Conda repositories and as a Galaxy tool from the public Galaxy Toolshed. In a public health context, BioHansel enables rapid and high-resolution classification of bacterial pathogens with low genetic diversity.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Clinical and public health utility of Mycobacterium tuberculosis whole genome sequencing

The World Health Organization (WHO) estimates that around 10 million people develop tuberculosis (TB) every year, with 1.5 million deaths attributed to TB in 2019 (World Health Organization, 2020). The majority of the disease burden occurs in low-income countries, where access to diagnostics and tailored treatment remains problematic. The current COVID-19 pandemic further threatens to impact global TB control by diverting resources, reducing notifications and hence significantly increasing deaths attributable to TB (World Health Organization, 2020). Whole genome sequencing (WGS) is becoming increasingly accessible, and has particular value in the diagnosis and management of TB disease (Cabibbe et al., 2018; Meehan et al., 2019). Not only does it have the potential to give more rapid and complete information on drug-resistance, but the high discriminatory power it offers allows detection of clusters and transmission pathways, as well as likely contamination events, mixed infections and to differentiate between re-infection and relapse with much greater confidence than previous typing methods.

In-host population dynamics of Mycobacterium tuberculosis complex during active disease

Tuberculosis (TB) is a leading cause of death globally. Understanding the population dynamics of TB’s causative agent Mycobacterium tuberculosis complex (Mtbc) in-host is vital for understanding the efficacy of antibiotic treatment. We use longitudinally collected clinical Mtbc isolates that underwent Whole-Genome Sequencing from the sputa of 200 patients to investigate Mtbc diversity during the course of active TB disease after excluding 107 cases suspected of reinfection, mixed infection or contamination. Of the 178/200 patients with persistent clonal infection >2 months, 27 developed new resistance mutations between sampling with 20/27 occurring in patients with pre-existing resistance. Low abundance resistance variants at a purity of ≥19% in the first isolate predict fixation in the subsequent sample. We identify significant in-host variation in 27 genes, including antibiotic resistance genes, metabolic genes and genes known to modulate host innate immunity and confirm several to be under positive selection by assessing phylogenetic convergence across a genetically diverse sample of 20,352 isolates.

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.

Discordant line probe genotypic testing vs culture-based drug susceptibility phenotypic testing in TB endemic KwaZulu-Natal: Impact on bedside clinical decision making

The recommendations for Mycobacterium tuberculosis drug susceptibility testing include both phenotypic and genotypic methods. This concurrent use of differing testing platforms has created an emerging challenge of discordant results, creating a diagnostic dilemma for the laboratorians as well as attending clinicians. We undertook a retrospective study to determine the prevalence of discordant results between the MTBDRplus line probe assay and solid culture-based drug susceptibility testing for rifampicin and isoniazid. The analysis was conducted for the period January 2013 and December 2015 at the Inkosi Albert Luthuli Central Hospital. Rifampicin and isoniazid resistance testing data were "paired" on 8273 isolates for culture-based drug susceptibility testing and line probe assay. The latter method showed high sensitivity and specificity of 93% and 95% respectively for isoniazid testing. For rifampicin testing, sensitivity and specificity were 95% and 75%. Overall, discordance was 14.6% for rifampicin and 7.2% for isoniazid. This report is not intended to determine superiority of one method over another. It is merely to show that discordance does exist between different methods of testing. Given the burden of HIV and Tuberculosis in Sub-Saharan Africa, these findings have clinical significance and huge public health implications. Clinicians should understand the limitations of phenotypic testing methods.

Mycobacterium bovis: From Genotyping to Genome Sequencing

Mycobacterium bovis is the main pathogen of bovine, zoonotic, and wildlife tuberculosis. Despite the existence of programs for bovine tuberculosis (bTB) control in many regions, the disease remains a challenge for the veterinary and public health sectors, especially in developing countries and in high-income nations with wildlife reservoirs. Current bTB control programs are mostly based on test-and-slaughter, movement restrictions, and post-mortem inspection measures. In certain settings, contact tracing and surveillance has benefited from M. bovis genotyping techniques. More recently, whole-genome sequencing (WGS) has become the preferential technique to inform outbreak response through contact tracing and source identification for many infectious diseases. As the cost per genome decreases, the application of WGS to bTB control programs is inevitable moving forward. However, there are technical challenges in data analyses and interpretation that hinder the implementation of M. bovis WGS as a molecular epidemiology tool. Therefore, the aim of this review is to describe M. bovis genotyping techniques and discuss current standards and challenges of the use of M. bovis WGS for transmission investigation, surveillance, and global lineages distribution. We compiled a series of associated research gaps to be explored with the ultimate goal of implementing M. bovis WGS in a standardized manner in bTB control programs.

Application of MIRU-VNTR on smear slides: a shortcut for detection of polyclonal infections in tuberculosis patients

Mixed (polyclonal) infections are one of the main problems in tuberculosis (TB) management. The best available method for detecting polyclonal infections in TB is mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR). According to multiple studies, MIRU-VNTR method can be applied to detect TB-related polyclonal infections in sputum samples or cultures. Setup of MIRU-VNTR on smear slides can be an efficient approach, regardless of the limitations of cultures and sputum samples in many laboratories. The present study aimed at investigating the diagnostic potential of MIRU-VNTR on smear slides in detecting mixed infections. Ziehl-Neelsen-stained microscopic slides were prepared from 14 clinical specimens. For amplifying 24 MIRU-VNTR loci, PCR assay was performed on the smear slides, clinical specimens, and cultures. Based on the 24-locus MIRU-VNTR analysis, polyclonal infections were reported in 42.85% of smear slides, while the corresponding rate was estimated at 57.1% (8/14) in the clinical samples. In the corresponding cultures, the rate of mixed infection was 7.14% (1/14). Use of smear slides can be a safe option for transferring clinical specimens between environmental and reference laboratories. Considering their significant impact on TB treatment, it is essential to diagnose mixed infections in low-resource countries with a high prevalence of mixed infections. The present findings show that direct MIRU-VNTR on smear slides can be conveniently used for the detection of mixed infections.