Whole genome sequencing identifies circulating Beijing-lineage Mycobacterium tuberculosis strains in Guatemala and an associated urban outbreak

Endemic transmission of a Mycobacterium tuberculosis L2.2.M3 sublineage of the L2 lineage within Colon, Panama: A prospective study

Mycobacterium tuberculosis lineage 2 (L2) remains a globally significant lineage associated with increased drug resistance and rapid transmission. The L2 lineage exhibits a hotspot for genetic diversity and evolution in Panama, requiring an in-depth analysis. We conducted a prospective analysis of 274 Mycobacterium tuberculosis L2 isolates from Colon City between January 2021 and October 2023. Drug resistance was determined using GeneXpert and MTBDRplus-Genotype assays, strain lineage was determined by strain-specific PCR (ASO-PCR), and whole-genome sequencing was conducted for phylogenetic analysis. Sequencing data were analyzed using the mtb-call2 pipeline and TB-gen tools to predict drug resistance and sublineage, respectively. Genome-wide single-nucleotide polymorphisms (SNPs) were used for phylogenetic and evolutionary analyses. ASO-PCR results identified all 31.7 % (86/271) isolates as Modern L2.2. WGS analysis of 66 strains confirmed all isolates belonged to the L2.2.1 sublineage. Sixty-four strains were analyzed in depth, with 96.9 % (62/64) classified as pan-susceptible and 3.1 % (2/64) as rifampicin/pyrazinamide-resistant. The sublineage analysis based on SNPs using the TB-gen tool identified a SNP at position 1219683G > A, which genotyped all 64 strains as L2.2.M3 sublineage. Phylogenetic analysis revealed a correlation with geographical distribution compared to other Latin American L2 isolates. Transmission clusters (≤12 SNPs) were identified and used to determine recent transmission events or TB transmission clusters. These analyses also confirmed a relatively low evolutionary rate within Panama L2 isolates and a highly conserved common ancestor shared with L2 isolates from Peru, Colombia, and Guatemala. These findings suggest endemic transmission of the Mycobacterium tuberculosis L2.2.M3 sublineage in Colon, Panama. We recommend combining genomic information with epidemiological data to accurately track and identify the source hotspot for the L2.2.M3 sublineage and focus control measures.

An ancestral mycobacterial effector promotes dissemination of infection

The human pathogen Mycobacterium tuberculosis typically causes lung disease but can also disseminate to other tissues. We identified a M. tuberculosis (Mtb) outbreak presenting with unusually high rates of extrapulmonary dissemination and bone disease. We found that the causal strain carried an ancestral full-length version of the type VII-secreted effector EsxM rather than the truncated version present in other modern Mtb lineages. The ancestral EsxM variant exacerbated dissemination through enhancement of macrophage motility, increased egress of macrophages from established granulomas, and alterations in macrophage actin dynamics. Reconstitution of the ancestral version of EsxM in an attenuated modern strain of Mtb altered the migratory mode of infected macrophages, enhancing their motility. In a zebrafish model, full-length EsxM promoted bone disease. The presence of a derived nonsense variant in EsxM throughout the major Mtb lineages 2, 3, and 4 is consistent with a role for EsxM in regulating the extent of dissemination.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Probable long-term prevalence for a predominant Mycobacterium tuberculosis clone of a Beijing genotype in Colon, Panama

Beijing genotype Mycobacterium tuberculosis strains associate with increased virulence, resistance and/or higher transmission rates. This study describes a specific Beijing strain predominantly identified in the Panamanian province of Colon with one of the highest incidences of tuberculosis in the country. Retrospective mycobacterial interspersed repetitive unit/variable number of tandem repeats analysis of 42 isolates collected between January and August 2018 allowed to identify a cluster (Beijing A) with 17 (40.5%) Beijing isolates. Subsequent prospective strain-specific PCR-based surveillance from September 2019 to March 2020 confirmed the predominance of the Beijing A strain (44.1%) in this province. Whole-genome sequencing revealed higher-than-expected diversity within the cluster, suggesting long-term prevalence of this strain and low number of cases caused by recent transmission. The Beijing A strain belongs to the Asian African 3 (Bmyc13, L2.2.5) branch of the modern Beijing sublineage, with their closest isolates corresponding to cases from Vietnam, probably introduced in Panama between 2000 and 2012.

Increasing Globalization and the Movement of Antimicrobial Resistance between Countries

Background: The threat of antimicrobial resistance continues to grow worldwide, exacerbated by poor antibiotic stewardship practices, limited development of new antimicrobial agents, and increasing globalization. Methods: This review covers previously published studies examining how human movement contributes to the global spread of antimicrobial resistance, including between low- and middle-income and high-income countries. Results: The emergence of resistance in one country or part of the world can become a worldwide event quickly. Human movement, including travel, medical tourism, military service, and migration, results in the globalization of resistant bacterial strains. Conclusions: Increased surveillance, whole-genome sequencing, focused infection control, and effective stewardship practices are needed to maintain the efficacy of antibiotics.

Genetic profiling of Mycobacterium tuberculosis revealed "modern" Beijing strains linked to MDR-TB from Southwestern Colombia

Beijing strains of Mycobacterium tuberculosis (lineage 2) have been associated with drug-resistance and transmission of tuberculosis worldwide. Most of the Beijing strains identified in the Colombian Pacific coast have exhibited a multidrug resistant (MDR) phenotype. We sought to evaluate the clonality and sublineage of Beijing strains circulating in Southwestern Colombia. Thirty-seven Beijing strains were identified through spoligotyping out of 311 clinical isolates collected in 9 years from 2002-2010. Further analysis by MIRU-VNTR 24 loci was conducted for the Beijing strains. For sublineage classification, deletions of RD105, RD207, and RD131 and point mutations at fbpB, mutT2, and acs were evaluated. Drug-resistance associated mutations to first- and second-line anti-TB drugs were also evaluated. Additionally, two Beijing strains were Illumina-whole genome sequenced (one MDR and one drug-susceptible). Among the 37 Beijing strains characterized, 36 belonged to the SIT190 type from which 28 were MDR, four pre-extensively drug resistant (XDR) TB, and four XDR-TB. The remaining strain was SIT1 and drug susceptible. MIRU-VNTR analysis allowed the identification of three Beijing clusters and two unique strains. Beijing strains were confirmed as "modern" sublineage. The mutations rpoB S531L and katG S315T were the most common among MDR strains. Moreover, the two strains evaluated by whole genome sequencing (WGS) shared most of the genetic features with the sublineage 2.2.1 "modern" Beijing previously characterized from Asian strains. WGS analysis of the MDR strain revealed the presence of eight SNPs previously reported in other MDR "Beijing-like" strains from Colombia. The presence of "modern" Beijing strains in Southwestern Colombia, most of them with MDR phenotype, suggests a different origin of this M. tuberculosis sublineage compared to other Beijing strains found in neighboring South American countries. This work may serve as a genetic baseline to study the evolution and spread of M. tuberculosis Beijing strains in Colombia, which play an important role in the propagation of MDR-TB.

Circulation of M. tuberculosis Beijing genotype in Latin America and the Caribbean

Lineage 2 (East Asian), which includes the Beijing genotype, is one of the most prevalent lineages of Mycobacterium tuberculosis (Mtb) throughout the world. The Beijing family is associated to hypervirulence and drug-resistant tuberculosis. The study of this genotype's circulation in Latin America is crucial for achieving total control of TB, the goal established by the World Health Organization, for the American sub-continent, before 2035. In this sense, the present work presents an overview of the status of the Beijing genotype for this region, with a bibliographical review, and data analysis of MIRU-VNTRs for available Beijing isolates. Certain countries present a prevalent trend of <5%, suggesting low transmissibility for the region, with the exception of Cuba (17.2%), Perú (16%) and Colombia (5%). Minimum Spanning Tree analysis, obtained from MIRU-VNTR data, shows distribution of specific clonal complex strains in each country. From this data, in most countries, we found that molecular epidemiology has not been a tool used for the control of TB, suggesting that the Beijing genotype may be underestimated in Latin America. It is recommended that countries with the highest incidence of the Beijing genotype use effective control strategies and increased care, as a requirement for public health systems.

Mycobacterial Evolution Intersects With Host Tolerance

Over the past 200 years, tuberculosis (TB) has caused more deaths than any other infectious disease, likely infecting more people than it has at any other time in human history. Mycobacterium tuberculosis (Mtb), the etiologic agent of TB, is an obligate human pathogen that has evolved through the millennia to become an archetypal human-adapted pathogen. This review focuses on the evolutionary framework by which Mtb emerged as a specialized human pathogen and applies this perspective to the emergence of specific lineages that drive global TB burden. We consider how evolutionary pressures, including transmission dynamics, host tolerance, and human population patterns, may have shaped the evolution of diverse mycobacterial genomes.

The use of whole-genome sequencing in cluster investigation of a multidrug-resistant tuberculosis outbreak

We used whole-genome sequencing (WGS) to delineate transmission networks and investigate the benefits of WGS during cluster investigation.We included clustered cases of multidrug-resistant (MDR) tuberculosis (TB)/extensively drug-resistant (XDR) TB linked by mycobacterial interspersed repetitive unit variable tandem repeat (MIRU-VNTR) strain typing or epidemiological information in the national cluster B1006, notified between 2007 and 2013 in the UK. We excluded from further investigation cases whose isolates differed by greater than 12 single nucleotide polymorphisms (SNPs). Data relating to patients' social networks were collected.27 cases were investigated and 22 had WGS, eight of which (36%) were excluded as their isolates differed by more than 12 SNPs to other cases. 18 cases were ruled into the transmission network based on genomic and epidemiological information. Evidence of transmission was inconclusive in seven out of 18 cases (39%) in the transmission network following WGS and epidemiological investigation.This investigation of a drug-resistant TB cluster illustrates the opportunities and limitations of WGS in understanding transmission in a setting with a high proportion of migrant cases. The use of WGS should be combined with classical epidemiological methods. However, not every cluster will be solvable, regardless of the quality of genomic data.

Annotated Genome Sequences of 16 Lineage 4 Mycobacterium tuberculosis Strains from Guatemala

Whole-genome sequencing has resulted in new insights into the phylogeography of Mycobacterium tuberculosis. However, only limited genomic data are available from M. tuberculosis strains in Guatemala. Here we report 16 complete genomes of clinical strains belonging to the Euro-American lineage 4, the most common lineage found in Guatemala and Central America.

Recent developments in genomics, bioinformatics and drug discovery to combat emerging drug-resistant tuberculosis

Emergence of drug-resistant tuberculosis (DR-TB) is a big challenge in TB control. The delay in diagnosis of DR-TB leads to its increased transmission, and therefore prevalence. Recent developments in genomics have enabled whole genome sequencing (WGS) of Mycobacterium tuberculosis (M. tuberculosis) from 3-day-old liquid culture and directly from uncultured sputa, while new bioinformatics tools facilitate to determine DR mutations rapidly from the resulting sequences. The present drug discovery and development pipeline is filled with candidate drugs which have shown efficacy against DR-TB. Furthermore, some of the FDA-approved drugs are being evaluated for repurposing, and this approach appears promising as several drugs are reported to enhance efficacy of the standard TB drugs, reduce drug tolerance, or modulate the host immune response to control the growth of intracellular M. tuberculosis. Recent developments in genomics and bioinformatics along with new drug discovery collectively have the potential to result in synergistic impact leading to the development of a rapid protocol to determine the drug resistance profile of the infecting strain so as to provide personalized medicine. Hence, in this review, we discuss recent developments in WGS, bioinformatics and drug discovery to perceive how they would transform the management of tuberculosis in a timely manner.

Detection and characterization of drug-resistant conferring genes in Mycobacterium tuberculosis complex strains: A prospective study in two distant regions of Ghana

We spoligotyped and screened 1490 clinical Mycobacterium tuberculosis complex strains from Northern and Greater Accra regions of Ghana against INH and RIF using the microplate alamar blue phenotypic assay. Specific drug resistance associated genetic elements of drug resistant strains were analyzed for mutations. A total of 111 (7.5%), 10 (0.7%) and 40 (2.6%) were mono-resistant to INH, RIF, and MDR, respectively. We found the Ghana spoligotype to be associated with drug resistance (INH: 22.1%; p = 0.0000, RIF: 6.2%; p = 0.0103, MDR: 4.6%; p = 0.0240) as compared to the Cameroon spoligotype (INH: 6.7%, RIF: 2.4%, MDR: 1.6%). The propensity for an isolate to harbour katG S315T mutation was higher in M. tuberculosis (75.8%) than Mycobacterium africanum (51.7%) (p = 0.0000) whereas the opposite was true for inhApro mutations; MAF (48.3%) compared to MTBSS (26.7%) (p = 0.0419). We identified possible novel compensatory INH resistance mutations in inhA (G204D) and ahpCpro (-88G/A and -142G/A) and a novel ndh mutation K32R. We detected two possible rpoC mutations (G332R and V483G), which occurred independently with rpoB S450L, respectively. The study provides the first evidence that associate the Ghana spoligotype with DR-TB and calls for further genome analyses for proper classification of this spoligotype and to explore for fitness implications and mechanisms underlying this observation.