Genomic Epidemiology of Tuberculosis

Integrating genomic and spatial analyses to describe tuberculosis transmission: a scoping review

Tuberculosis remains a leading cause of infection-related mortality, and efforts to reduce its incidence have been hindered by an incomplete understanding of local Mycobacterium tuberculosis transmission dynamics. Advances in pathogen sequencing and spatial analysis have created new opportunities to map M tuberculosis transmission patterns more precisely. In this scoping review, we searched for studies combining pathogen genetics and location data to analyse the spatial patterns of M tuberculosis transmission and identified 142 studies published between 1994 and 2024. Secular changes in genetic methods were observed, with genome sequencing approaches largely replacing lower-resolution genotyping methods since 2020. The included studies addressed four primary research questions: how are tuberculosis cases and M tuberculosis transmission clusters geographically distributed; do spatially concentrated M tuberculosis clusters exist, and where are these areas located; when spatial concentration occurs, what host, pathogen, or environmental factors contribute to these patterns; and do identifiable relationships exist between the spatial proximity of tuberculosis cases and the genetic similarity of the M tuberculosis isolates infecting these individuals? Collectively, in this Review, we examined the available study data, evaluated the analytical requirements for addressing these questions, and discussed opportunities and challenges for future research. We found that the integration of spatial and genomic data can inform a detailed understanding of local M tuberculosis transmission patterns, but improved study designs and new analytical methods to address gaps in sampling completeness and to integrate additional movement data are needed to fully realise the potential of these tools.

Genetic diversities and drug resistance in Mycobacterium bovis isolates from zoonotic tuberculosis using whole genome sequencing

BACKGROUND

Zoonotic human tuberculosis (TB) caused by Mycobacterium bovis (M. bovis) is as vital as Mycobacterium tuberculosis, however with scarce available information. We aimed to use whole-genome sequencing (WGS) technology to take a deep insight into the circulating genotypes of human M. bovis and the genomic characteristics underlying virulence and drug resistance.

METHODS

The study included smear positive Ziehl-Neelsen samples from patients with suspected tuberculosis. Samples were cultured on Lowenstein-Jensen media and suspected colonies of M. bovis were selected to undergo DNA extraction and WGS. Data was analysed using the Bacterial and Viral Bioinformatics Resource Center (BV-BRC), and online bioinformatics tools. A phylogenetic tree was constructed for our sequenced strains, in addition to a set of 59 previously sequenced M. bovis genomes from different hosts and countries.

RESULTS

Out of total 112 mycobacterial positive cultures, five M. bovis were isolated and underwent WGS. All sequenced strains belonged to Mycobacterium tuberculosis var bovis, spoligotype BOV_1; BOV_11. Resistance gene mutations were determined in 100% of strains to pyrazinamide (pncA and rpsA), isoniazid (KatG and ahpC), ethambutol (embB, embC, embR and ubiA), streptomycin (rpsl) and fluoroquinolones (gyrA and gyrB). Rifampin (rpoB and rpoC) and delamanid (fbiC) resistance genes were found in 80% of strains. The major represented virulence classes were the secretion system, cell surface components and regulation system. The phylogenetic analysis revealed close genetic relatedness of three sequenced M. bovis strains to previous reported cow strains from Egypt and human strains from France, as well as relatedness of one M. bovis strain to four human Algerian strains. One sequenced strain was related to one cow strain from Egypt and a human strain from South Africa.

CONCLUSIONS

All sequenced M. bovis isolates showed the same spoligotype, but diverse phylogeny. Resistance gene mutations were detected for anti-TB drugs including pyrazinamide, isoniazid, streptomycin, ethambutol, fluoroquinolones, cycloserine, rifampin and delamanid. The virulence profile comprised genes assigned mainly to secretion system, cell surface components and regulation system. Phylogenetic analysis revealed genetic relatedness between our isolates and previously sequenced bovine strains from Egypt as well as human strains from other nearby countries in the region.

Genomic analysis of Clostridium perfringens type D isolates from goat farms

C. perfringens type D strains are the leading cause of enterotoxaemia in ruminants such as goats, sheep, and cattle. However, there has been no prior research on the genomic characteristics of C. perfringens type D strains from various regions in China. Here, we investigated the antibiotic resistance, genomic characteristics, and phylogenetic relationship of C. perfringens type D isolates recovered from goat farms in Shaanxi, Gansu, and Ningxia provinces. The antibiotic resistance test indicated that the isolates displayed high minimum inhibitory concentration (MIC) values to sulfafurazole, whereas the other antibiotics tested, such as penicillin, enrofloxacin, and florfenicol, worked well on them. Additionally, only tetracycline resistance genes [tetA(P) and tetB(P)] were identified from the isolates. A collective of 13 toxin genes, including etx and cpe were detected among the isolates. Sequence comparison revealed that the etx and cpe genes shared high sequence identities, and they could coexist on a pCW3-like plasmid, representing a potential risk to both animal breeding and public health. Phylogenetic analysis using core genome multi-locus sequence typing (cgMLST) and core genome single nucleotide polymorphisms (SNPs) revealed the close genetic relationship and potential regional/transregional transmission of the C. perfringens type D isolates in Shaanxi and Gansu provinces. Furthermore, pan-genomic analysis suggested the functional differences at the protein-coding gene level, although isolates from the same source shared a close genetic relationship. In conclusion, this study indicated the antibiotic resistance, virulence markers, potential transregional transmission, and genomic diversity of C. perfringens type D strains from various regions in China, which could provide references for the prevention of C. perfringens foodborne diseases and further research.

Genetic Characterization and Population Structure of Drug-Resistant Mycobacterium tuberculosis Isolated from Brazilian Patients Using Whole-Genome Sequencing

The present study aimed to determine the genetic diversity of isolates of Mycobacterium tuberculosis (Mtb) from presumed drug-resistant tuberculosis patients from several states of Brazil. The isolates had been submitted to conventional drug susceptibility testing for first- and second-line drugs. Multidrug-resistant (MDR-TB) (54.8%) was the most frequent phenotypic resistance profile, in addition to an important high frequency of pre-extensive resistance (p-XDR-TB) (9.2%). Using whole-genome sequencing (WGS), we characterized 298 Mtb isolates from Brazil. Besides the analysis of genotype distribution and possible correlations between molecular and clinical data, we determined the performance of an in-house WGS pipeline with other online pipelines for Mtb lineages and drug resistance profile definitions. Sub-lineage 4.3 (52%) was the most frequent genotype, and the genomic approach revealed a p-XDR-TB level of 22.5%. We detected twenty novel mutations in three resistance genes, and six of these were observed in eight phenotypically resistant isolates. A cluster analysis of 170 isolates showed that 43.5% of the TB patients belonged to 24 genomic clusters, suggesting considerable ongoing transmission of DR-TB, including two interstate transmissions. The in-house WGS pipeline showed the best overall performance in drug resistance prediction, presenting the best accuracy values for five of the nine drugs tested. Significant associations were observed between suffering from fatal disease and genotypic p-XDR-TB (p = 0.03) and either phenotypic (p = 0.006) or genotypic (p = 0.0007) ethambutol resistance. The use of WGS analysis improved our understanding of the population structure of MTBC in Brazil and the genetic and clinical data correlations and demonstrated its utility for surveillance efforts regarding the spread of DR-TB, hopefully helping to avoid the emergence of even more resistant strains and to reduce TB incidence and mortality rates.

Borderline rpoB mutations transmit at the same rate as common rpoB mutations in a tuberculosis cohort in Bangladesh

The spread of multidrug-resistant tuberculosis (MDR-TB) is a growing problem in many countries worldwide. Resistance to one of the primary first-line drugs, rifampicin, is caused by mutations in the Mycobacterium tuberculosis rpoB gene. So-called borderline rpoB mutations confer low-level resistance, in contrast to more common rpoB mutations which confer high-level resistance. While some borderline mutations show lower fitness in vitro than common mutations, their in vivo fitness is currently unknown. We used a dataset of 394 whole genome sequenced MDR-TB isolates from Bangladesh, representing around 44 % of notified MDR-TB cases over 6 years, to look at differences in transmission clustering between isolates with borderline rpoB mutations and those with common rpoB mutations. We found a relatively low percentage of transmission clustering in the dataset (34.8 %) but no difference in clustering between different types of rpoB mutations. Compensatory mutations in rpoA, rpoB, and rpoC were associated with higher levels of transmission clustering as were lineages two, three, and four relative to lineage one. Young people as well as patients with high sputum smear positive TB were more likely to be in a transmission cluster. Our findings show that although borderline rpoB mutations have lower in vitro growth potential this does not translate into lower transmission potential or in vivo fitness. Proper detection of these mutations is crucial to ensure they do not go unnoticed and spread MDR-TB within communities.

Molecular Epidemiology and Genetic Diversity of Multidrug-Resistant Mycobacterium tuberculosis Isolates in Bangladesh

Although the number of multidrug-resistant (MDR) tuberculosis (TB) cases is high overall, a major gap exists in our understanding of the molecular characteristics and transmission dynamics of the MDR Mycobacterium tuberculosis isolates circulating in Bangladesh. The present study aims to characterize the MDR-TB isolates of Bangladesh and to investigate the mode of transmission. A total of 544 MDR-TB isolates were obtained from a nationwide drug-resistant TB surveillance study conducted between October 2011 and March 2017 covering all geographic divisions of Bangladesh. The isolates were characterized using TbD1 deletion analysis, spoligotyping, and mycobacterial interspersed repetitive-unit–variable-number tandem-repeat (MIRU-VNTR) typing. Deletion analysis showed that 440 (80.9%) isolates were the modern type, while the remainder were the ancestral type. The largest circulating lineage was the Beijing type, comprising 208 isolates (38.2%), followed by T, EAI, and LAM with 93 (17.1%), 58 (10.7%), and 52 (9.5%) isolates, respectively. Combined MIRU-VNTR and spoligotyping analysis demonstrated that the majority of the clustered isolates were of the Beijing and T1 lineages. The overall rate of recent transmission was estimated at 33.8%. In conclusion, the MDR M. tuberculosis isolates circulating in Bangladesh are mostly of the modern virulent type. The Beijing and T lineages are the predominant types and most of the transmission of MDR-TB can be attributed to them. The findings also suggest that, along with the remarkable transmission, the emergence of MDR-TB in Bangladesh is largely due to acquired resistance. Rapid and accurate diagnosis and successful treatment will be crucial for controlling MDR-TB in Bangladesh.

IMPORTANCE Multidrug-resistant TB is considered to be the major threat to tuberculosis control activities worldwide, including in Bangladesh. Despite the fact that the number of MDR-TB cases is high, a major gap exists in our understanding of the molecular epidemiology of the MDR-TB isolates in Bangladesh. In our study, we characterized and classified the MDR-TB isolates circulating in Bangladesh and investigated their mode of transmission. Our results demonstrated that the MDR M. tuberculosis isolates circulating in Bangladesh are mostly of the modern virulent type. The Beijing and T lineages are the predominant types and are implicated in the majority of MDR-TB transmission. Our findings reveal that, along with the remarkable transmission, the emergence of MDR-TB in Bangladesh is largely due to acquired resistance, which may be due to nonadherence to treatment or inadequate treatment of TB patients. Rapid diagnosis and adherence to an appropriate treatment regimen are therefore crucial to controlling MDR-TB in Bangladesh.

Genitourinary Tuberculosis: A Comprehensive Review of a Neglected Manifestation in Low-Endemic Countries

Genitourinary tuberculosis (GUTB) represents a disease often underestimated by urological specialists, particularly in settings such as the European one, where the pathology is less frequent. Similar to other uncommon diseases at these latitudes, GUTB is a neglected clinical problem. In this light, the aim of this review is to give a comprehensive overview of GUTB in order to provide a useful tool for urologists who seldomly manage this disease. A non-systematic review of genitourinary tuberculosis was performed on relevant articles published from January 1990 to July 2021 using PubMed, Scopus, and the Cochrane Central Register of Controlled Trials. GUTB represents up to a quarter of extrapulmonary tuberculosis (EPTB) cases. Diagnostic, therapeutic and surgical work-up have been deeply reviewed and summarized. The mass migration of refugees to Europe as well as the ease of international travel is gradually leading to an upsurge in urological diseases such as GUTB, which were previously only rarely encountered in some European countries. The poor TB knowledge of European urologists should be improved through medical education courses, webinars or telematic means.

The within-host evolution of antimicrobial resistance in Mycobacterium tuberculosis

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

Use of next generation sequencing technologies for the diagnosis and epidemiology of infectious diseases

For the first time, next generation sequencing technologies provide access to genomic information at a price and scale that allow their implementation in routine clinical practice and epidemiology. While there are still many obstacles to their implementation, there are also multiple examples of their major advantages compared with previous methods. Their main advantage is that a single determination allows epidemiological information on the causative microorganism to be obtained simultaneously, as well as its resistance profile, although these advantages vary according to the pathogen under study. This review discusses several examples of the clinical and epidemiological use of next generation sequencing applied to complete genomes and microbiomes and reflects on its future in clinical practice.

Molecular epidemiology of Mycobacterium tuberculosis in Brazil before the whole genome sequencing era: a literature review

Molecular-typing can help in unraveling epidemiological scenarios and improvement for disease control strategies. A literature review of Mycobacterium tuberculosis transmission in Brazil through genotyping on 56 studies published from 1996-2019 was performed. The clustering rate for mycobacterial interspersed repetitive units - variable tandem repeats (MIRU-VNTR) of 1,613 isolates were: 73%, 33% and 28% based on 12, 15 and 24-loci, respectively; while for RFLP-IS6110 were: 84% among prison population in Rio de Janeiro, 69% among multidrug-resistant isolates in Rio Grande do Sul, and 56.2% in general population in São Paulo. These findings could improve tuberculosis (TB) surveillance and set up a solid basis to build a database of Mycobacterium genomes.

TB DEPOT (Data Exploration Portal): A multi-domain tuberculosis data analysis resource

The NIAID TB Portals Program (TBPP) established a unique and growing database repository of socioeconomic, geographic, clinical, laboratory, radiological, and genomic data from patient cases of drug-resistant tuberculosis (DR-TB). Currently, there are 2,428 total cases from nine country sites (Azerbaijan, Belarus, Moldova, Georgia, Romania, China, India, Kazakhstan, and South Africa), 1,611 (66%) of which are multidrug- or extensively-drug resistant and 1,185 (49%), 863 (36%), and 952 (39%) of which contain X-ray, computed tomography (CT) scan, and genomic data, respectively. We introduce the Data Exploration Portal (TB DEPOT, https://depot.tbportals.niaid.nih.gov) to visualize and analyze these multi-domain data. The TB DEPOT leverages the TBPP integration of clinical, socioeconomic, genomic, and imaging data into standardized formats and enables user-driven, repeatable, and reproducible analyses. It furthers the TBPP goals to provide a web-enabled analytics platform to countries with a high burden of multidrug-resistant TB (MDR-TB) but limited IT resources and inaccessible data, and enables the reusability of data, in conformity with the NIH's Findable, Accessible, Interoperable, and Reusable (FAIR) principles. TB DEPOT provides access to "analysis-ready" data and the ability to generate and test complex clinically-oriented hypotheses instantaneously with minimal statistical background and data processing skills. TB DEPOT is also promising for enhancing medical training and furnishing well annotated, hard to find, MDR-TB patient cases. TB DEPOT, as part of TBPP, further fosters collaborative research efforts to better understand drug-resistant tuberculosis and aid in the development of novel diagnostics and personalized treatment regimens.

The relationship between transmission time and clustering methods in Mycobacterium tuberculosis epidemiology

BACKGROUND

Tracking recent transmission is a vital part of controlling widespread pathogens such as Mycobacterium tuberculosis. Multiple methods with specific performance characteristics exist for detecting recent transmission chains, usually by clustering strains based on genotype similarities. With such a large variety of methods available, informed selection of an appropriate approach for determining transmissions within a given setting/time period is difficult.

METHODS

This study combines whole genome sequence (WGS) data derived from 324 isolates collected 2005-2010 in Kinshasa, Democratic Republic of Congo (DRC), a high endemic setting, with phylodynamics to unveil the timing of transmission events posited by a variety of standard genotyping methods. Clustering data based on Spoligotyping, 24-loci MIRU-VNTR typing, WGS based SNP (Single Nucleotide Polymorphism) and core genome multi locus sequence typing (cgMLST) typing were evaluated.

FINDINGS

Our results suggest that clusters based on Spoligotyping could encompass transmission events that occurred almost 200 years prior to sampling while 24-loci-MIRU-VNTR often represented three decades of transmission. Instead, WGS based genotyping applying low SNP or cgMLST allele thresholds allows for determination of recent transmission events, e.g. in timespans of up to 10 years for a 5 SNP/allele cut-off.

INTERPRETATION

With the rapid uptake of WGS methods in surveillance and outbreak tracking, the findings obtained in this study can guide the selection of appropriate clustering methods for uncovering relevant transmission chains within a given time-period. For high resolution cluster analyses, WGS-SNP and cgMLST based analyses have similar clustering/timing characteristics even for data obtained from a high incidence setting.

Mycobacterium tuberculosis Acquires Limited Genetic Diversity in Prolonged Infections, Reactivations and Transmissions Involving Multiple Hosts

Background:Mycobacterium tuberculosis (MTB) has limited ability to acquire variability. Analysis of its microevolution might help us to evaluate the pathways followed to acquire greater infective success. Whole-genome sequencing (WGS) in the analysis of the transmission of MTB has elucidated the magnitude of variability in MTB. Analysis of transmission currently depends on the identification of clusters, according to the threshold of variability (<5 SNPs) between isolates. Objective: We evaluated whether the acquisition of variability in MTB, was more frequent in situations which could favor it, namely intrapatient, prolonged infections or reactivations and interpatient transmissions involving multiple sequential hosts. Methods: We used WGS to analyze the accumulation of variability in sequential isolates from prolonged infections or translations from latency to reactivation. We then measured microevolution in transmission clusters with prolonged transmission time, high number of involved cases, simultaneous involvement of latency and active transmission. Results: Intrapatient and interpatient acquisition of variability was limited, within the ranges expected according to the thresholds of variability proposed, even though bursts of variability were observed. Conclusions: The thresholds of variability proposed for MTB seem to be valid in most circumstances, including those theoretically favoring acquisition of variability. Our data point to multifactorial modulation of microevolution, although further studies are necessary to elucidate the factors underlying this modulation.