Novel DNA chip based on a modified DigiTag2 assay for high-throughput species identification and genotyping of Mycobacterium tuberculosis complex isolates

The MAGMA pipeline for comprehensive genomic analyses of clinical Mycobacterium tuberculosis samples

BACKGROUND

Whole genome sequencing (WGS) holds great potential for the management and control of tuberculosis. Accurate analysis of samples with low mycobacterial burden, which are characterized by low (<20x) coverage and high (>40%) levels of contamination, is challenging. We created the MAGMA (Maximum Accessible Genome for Mtb Analysis) bioinformatics pipeline for analysis of clinical Mtb samples.

METHODS AND RESULTS

High accuracy variant calling is achieved by using a long seedlength during read mapping to filter out contaminants, variant quality score recalibration with machine learning to identify genuine genomic variants, and joint variant calling for low Mtb coverage genomes. MAGMA automatically generates a standardized and comprehensive output of drug resistance information and resistance classification based on the WHO catalogue of Mtb mutations. MAGMA automatically generates phylogenetic trees with drug resistance annotations and trees that visualize the presence of clusters. Drug resistance and phylogeny outputs from sequencing data of 79 primary liquid cultures were compared between the MAGMA and MTBseq pipelines. The MTBseq pipeline reported only a proportion of the variants in candidate drug resistance genes that were reported by MAGMA. Notable differences were in structural variants, variants in highly conserved rrs and rrl genes, and variants in candidate resistance genes for bedaquiline, clofazmine, and delamanid. Phylogeny results were similar between pipelines but only MAGMA visualized clusters.

CONCLUSION

The MAGMA pipeline could facilitate the integration of WGS into clinical care as it generates clinically relevant data on drug resistance and phylogeny in an automated, standardized, and reproducible manner.

A Bioinformatics Whole-Genome Sequencing Workflow for Clinical Mycobacterium tuberculosis Complex Isolate Analysis, Validated Using a Reference Collection Extensively Characterized with Conventional Methods and In Silico Approaches

The use of whole-genome sequencing (WGS) for routine typing of bacterial isolates has increased substantially in recent years. For Mycobacterium tuberculosis (MTB), in particular, WGS has the benefit of drastically reducing the time required to generate results compared to most conventional phenotypic methods. Consequently, a multitude of solutions for analyzing WGS MTB data have been developed, but their successful integration in clinical and national reference laboratories is hindered by the requirement for their validation, for which a consensus framework is still largely absent. We developed a bioinformatics workflow for (Illumina) WGS-based routine typing of MTB complex (MTBC) member isolates allowing complete characterization, including (sub)species confirmation and identification (16S, csb/RD, hsp65), single nucleotide polymorphism (SNP)-based antimicrobial resistance (AMR) prediction, and pathogen typing (spoligotyping, SNP barcoding, and core genome multilocus sequence typing). Workflow performance was validated on a per-assay basis using a collection of 238 in-house-sequenced MTBC isolates, extensively characterized with conventional molecular biology-based approaches supplemented with public data. For SNP-based AMR prediction, results from molecular genotyping methods were supplemented with in silico modified data sets, allowing us to greatly increase the set of evaluated mutations. The workflow demonstrated very high performance with performance metrics of >99% for all assays, except for spoligotyping, where sensitivity dropped to ∼90%. The validation framework for our WGS-based bioinformatics workflow can aid in the standardization of bioinformatics tools by the MTB community and other SNP-based applications regardless of the targeted pathogen(s). The bioinformatics workflow is available for academic and nonprofit use through the Galaxy instance of our institute at https://galaxy.sciensano.be.

Cas12a/Guide RNA-Based Platform for Rapid and Accurate Identification of Major Mycobacterium Species

Mycobacterium tuberculosis infection and nontuberculous mycobacteria (NTM) infections exhibit similar clinical symptoms; however, the therapies for these two types of infections are different. Therefore, the rapid and accurate identification of M. tuberculosis and NTM species is very important for the control of tuberculosis and NTM infections. In the present study, a Cas12a/guide RNA (gRNA)-based platform was developed to identify M. tuberculosis and most NTM species. By designing species-specific gRNA probes targeting the rpoB sequence, a Cas12a/gRNA-based platform successfully identified M. tuberculosis and six major NTM species (Mycobacterium abscessus, Mycobacterium intracellulare, Mycobacterium avium, Mycobacterium kansasii, Mycobacterium gordonae, and Mycobacterium fortuitum) without cross-reactivity. In a blind assessment, a total of 72 out of 73 clinical Mycobacterium isolates were correctly identified, which is consistent with previous rpoB sequencing results. These results suggest that the Cas12a/gRNA-based platform is a promising tool for the rapid, accurate, and cost-effective identification of both M. tuberculosis and NTM species.

Mycobacterium bovis and Other Uncommon Members of the Mycobacterium tuberculosis Complex

Since its discovery by Theobald Smith, Mycobacterium bovis has been a human pathogen closely related to animal disease. At present, M. bovis tuberculosis is still a problem of importance in many countries and is considered the main cause of zoonotic tuberculosis throughout the world. Recent development of molecular epidemiological tools has helped us to improve our knowledge about transmission patterns of this organism, which causes a disease indistinguishable from that caused by Mycobacterium tuberculosis. Diagnosis and treatment of this mycobacterium are similar to those for conventional tuberculosis, with the important exceptions of constitutive resistance to pyrazinamide and the fact that multidrug-resistant and extremely drug-resistant M. bovis strains have been described. Among other members of this complex, Mycobacterium africanum is the cause of many cases of tuberculosis in West Africa and can be found in other areas mainly in association with immigration. M. bovis BCG is the currently available vaccine for tuberculosis, but it can cause disease in some patients. Other members of the M. tuberculosis complex are mainly animal pathogens with only exceptional cases of human disease, and there are even some strains, like "Mycobacterium canettii," which is a rare human pathogen that could have an important role in the knowledge of the evolution of tuberculosis in the history.

Molecular Epidemiology of Mycobacterium bovis in Humans and Cattle

Bovine tuberculosis (bTB), caused by Mycobacterium bovis (M. bovis), is a serious re-emerging disease in both animals and humans. The evolution of the Multi- and Extensively drug-resistant M. bovis strains (MDR-TB and XDR-TB) represents a global threat to public health. Worldwide, the disease is responsible for great economic losses in the veterinary field, serious threat to the ecosystem, and about 3.1% of human TB cases, up to 16% in Tanzania. Only thorough investigation to understand the pathogen's epidemiology can help in controlling the disease and minimizing its threat. For this purpose, various tools have been developed for use in advanced molecular epidemiological studies of bTB, either alone or in combination with standard conventional epidemiological approaches. These techniques enable the analysis of the intra- and inter-species transmission dynamics of bTB. The delivered data can reveal detailed insights into the source of infection, correlations among human and bovine isolates, strain diversity and evolution, spread, geographical localization, host preference, tracing of certain virulence factors such as antibiotic resistance genes, and finally the risk factors for the maintenance and spread of M. bovis. They also allow for the determination of epidemic and endemic strains. This, in turn, has a significant diagnostic impact and helps in vaccine development for bTB eradication programs. The present review discusses many topics including the aetiology, epidemiology and importance of M. bovis, the prevalence of bTB in humans and animals in various countries, the molecular epidemiology of M. bovis, and finally applied molecular epidemiological techniques.

Genetic structure and relationships of 16 Asian and European cattle populations using DigiTag2 assay

In this study, we genotyped 117 autosomal single nucleotide polymorphisms using a DigiTag2 assay to assess the genetic diversity, structure and relationships of 16 Eurasian cattle populations, including nine cattle breeds and seven native cattle. Phylogenetic and principal component analyses showed that Bos taurus and Bos indicus populations were clearly distinguished, whereas Japanese Shorthorn and Japanese Polled clustered with European populations. Furthermore, STRUCTURE analysis demonstrated the distinct separation between Bos taurus and Bos indicus (K=2), and between European and Asian populations (K=3). In addition, Japanese Holstein exhibited an admixture pattern with Asian and European cattle (K=3‐5). Mongolian (K=13‐16) and Japanese Black (K=14‐16) populations exhibited admixture patterns with different ancestries. Bos indicus populations exhibited a uniform genetic structure at K=2‐11, thereby suggesting that there are close genetic relationships among Bos indicus populations. However, the Bhutan and Bangladesh populations formed a cluster distinct from the other Bos indicus populations at K=12‐16. In conclusion, our study could sufficiently explain the genetic construction of Asian cattle populations, including: (i) the close genetic relationships among Bos indicus populations; (ii) the genetic influences of European breeds on Japanese breeds; (iii) the genetic admixture in Japanese Holstein, Mongolian and Japanese Black cattle; and (iv) the genetic subpopulations in Southeast Asia.

Genetic diversity and dynamic distribution of Mycobacterium tuberculosis isolates causing pulmonary and extrapulmonary tuberculosis in Thailand

This study examined the genetic diversity and dynamicity of circulating Mycobacterium tuberculosis strains in Thailand using nearly neutral molecular markers. The single nucleotide polymorphism (SNP)-based genotypes of 1,414 culture-positive M. tuberculosis isolates from 1,282 pulmonary tuberculosis (PTB) and 132 extrapulmonary TB (EPTB) patients collected from 1995 to 2011 were characterized. Among the eight SNP cluster groups (SCG), SCG2 (44.1%), which included the Beijing (BJ) genotype, and SCG1 (39.4%), an East African Indian genotype, were dominant. Comparisons between the genotypes of M. tuberculosis isolates causing PTB and EPTB in HIV-negative cases revealed similar prevalence trends although genetic diversity was higher in the PTB patients. The identification of 10 reported sequence types (STs) and three novel STs was hypothesized to indicate preferential expansion of the SCG2 genotype, especially the modern BJ ST10 (15.6%) and ancestral BJ ST19 (13.1%). An association between SCG2 and SCG1 genotypes and particular patient age groups implies the existence of different genetic advantages among the bacterial populations. The results revealed that increasing numbers of young patients were infected with M. tuberculosis SCGs 2 and 5, which contrasts with the reduction of the SCG1 genotype. Our results indicate the selection and dissemination of potent M. tuberculosis genotypes in this population. The determination of heterogeneity and dynamic population changes of circulating M. tuberculosis strains in countries using the Mycobacterium bovis BCG (bacillus Calmette-Guérin) vaccine are beneficial for vaccine development and control strategies.