Higher genome mutation rates of Beijing lineage of Mycobacterium tuberculosis during human infection

Mycobacterium tuberculosis is less likely to acquire pathogenic mutations during latent infection than during active disease

Most people infected with Mycobacterium tuberculosis (Mtb) are believed to be in a state of latent tuberculosis (TB) infection (LTBI). Although LTBI is asymptomatic and not infectious, there is a risk of developing active disease even decades after infection. Here, to characterize mutations acquired during LTBI, we collected and analyzed Mtb genomes from seven Japanese patient pairs, each pair consisting of two active TB patients whose starting dates of developing active disease were >3 years apart; one had a high suspicion of LTBI before developing active disease, whereas the other did not. Thereafter, we compared these genomes with those of longitudinal sample pairs within a host of chronic active TB infections combined with public data. The bacterial populations in patients with LTBI were genetically more homogeneous and accumulated single nucleotide polymorphisms (SNPs) slower than those from active disease. Moreover, the lower proportion of nonsynonymous SNPs indicated weaker selective pressures during LTBI than active disease. Finally, the different mutation spectrums indicated different mutators between LTBI and active disease. These results suggest that the likelihood of the acquisition of mutations responsible for antibiotic resistance and increased virulence was lower in the Mtb population from LTBI than active disease.IMPORTANCEControlling latent tuberculosis (TB) infection (LTBI) activation is an effective strategy for TB elimination, where understanding Mycobacterium tuberculosis (Mtb) dynamics within the host plays an important role. Previous studies on chronic active disease reported that Mtb accumulated genomic mutations within the host, possibly resulting in acquired drug resistance and increased virulence. However, several reports suggest that fewer mutations accumulate during LTBI than during the active disease, but the associated risk is largely unknown. Here, we analyzed the genomic dynamics of Mtb within the host during LTBI. Our results statistically suggest that Mtb accumulates mutations during LTBI, but most mutations are under low selective pressures, which induce mutations responsible for drug resistance and virulence. Thus, we propose that LTBI acts as a source for new TB disease rather than as a period for in-host genome evolution.

Delineating the evolutionary pathway to multidrug-resistant outbreaks of a Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage

OBJECTIVES

We sought to capture the evolutionary itinerary of the Mycobacterium tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, where it caused a major multidrug-resistant (MDR) tuberculosis outbreak in a context strictly negative for HIV infection.

METHODS

We combined whole genome sequencing and Bayesian approaches using a representative collection of drug-susceptible and drug-resistant L4.1.2.1/Haarlem clinical strains (n = 121) recovered from the outbreak region over 16 years.

RESULTS

In the absence of drug resistance, the L4.1.2.1/Haarlem sublineage showed a propensity for rapid transmission as witnessed by the high clustering (44.6%) and recent transmission rates (25%), as well as the reduced mean distance between genome pairs. The entire pool of L4.1.2.1/Haarlem MDR strains was found to be linked to either the aforementioned major outbreak (68 individuals, 2001-2016) or to a minor, newly uncovered outbreak (six cases, 2001-2011). Strikingly, the two outbreaks descended independently from a common ancestor that can be dated back to 1886.

CONCLUSIONS

Our data point to the intrinsic propensity for rapid transmission of the M. tuberculosis L4.1.2.1/Haarlem sublineage in northern Tunisia, linking the overall MDR tuberculosis epidemic to a single ancestor. These findings bring out the important role of the bacillus' genetic background in the emergence of successful MDR M. tuberculosis clones.

Mycobacterium tuberculosis Beijing in the State of Nuevo Leon, Mexico

Tuberculosis remains a serious threat to human health as an infectious disease in Mexico. Data about the genotypes of circulating Mycobacterium tuberculosis isolates (MTB) in the State of Nuevo Leon, Mexico are scarce. We aimed to determine the genotypes of circulating MTB belonging to the Beijing lineage recovered from patients in the State of Nuevo Leon, Mexico. A total of 406 MTB isolates from this state were genotyped using the spoligotyping method and 18-locus MIRU-VNTR. Lineage classification and MTB transmission analysis were performed. Based on the spoligotyping analysis, we found 24 strains belonging to the Beijing genotype that were characterized phylogenetically. The MIRUs showed greater discriminatory power than the standard RFLP-IS6110 method; therefore, the greatest allelic diversity among the Beijing strains was observed with MIRU10, MIRU31, MIRU39, MRU40, and MIRU 26. MVLA analysis showed a profile variation between Beijing and non-Beijing strains. The minimum spanning tree (MST) showed that 79% (19) of the strains are related. All Beijing strains exhibited the deletion of region TbD1, which is a characteristic of modern strains. The application of spoligotyping and MIRU-VNTR-18 methods together proved to be more sensitive, discriminatory, and rapid than the standard method for the epidemiological analysis of Mycobacterium Beijing isolates. This study is one of the first to describe the genomic diversity of M. Beijing in the State of Nuevo Leon, Mexico.

Define SNP thresholds for delineation of tuberculosis transmissions using whole-genome sequencing

For facilitating tuberculosis (TB) control, we used a whole-genome sequencing (WGS)-based approach to delineate transmission networks in a country with an intermediate burden of TB. A cluster was defined as Mycobacterium tuberculosis isolates with identical genotypes, and an outbreak was defined as clustered cases with epidemiological links (epi-links). To refine a cluster predefined using space oligonucleotide typing and mycobacterial interspersed repetitive unit variable tandem repeat typing, we analyzed one pansusceptible TB (C1) and three multidrug-resistant (MDR)-TB (C2-C4) clusters from different scenarios. Pansusceptible TB cluster (C1) consisting of 28 cases had ≤5 single nucleotide polymorphisms (SNPs) difference between their isolates. C1 was a definite outbreak, with cases attending the same junior high school in 2012. Three MDR-TB clusters (C2-C4) with distinct genotypes were identified, each consisting of 12-22 cases. Some of the cases had either ≤5 or ≤15 SNPs difference with clear or probable epi-links. Of note, even though WGS could effectively assist TB contact tracing, we still observed missing epi-links in some cases within the same cluster. Our results showed that thresholds of ≤5 and ≤15 SNPs difference between isolates were used to categorize definite and probable TB transmission, respectively. Furthermore, a higher SNP threshold might be required to define an MDR-TB outbreak. WGS still needs to be combined with classical epidemiological methods for improving outbreak investigations. Importantly, different SNP thresholds have to be applied to define outbreaks.

IMPORTANCE

TB is a chronic disease. Depending on host factors and TB burden, clusters of cases may continue to increase for several years. Conventional genotyping methods overestimate TB transmission, hampering precise detection of outbreaks and comprehensive surveillance. WGS can be used to obtain SNP information of M. tuberculosis to improve discriminative limitations of conventional methods and to strengthen delineation of transmission networks. It is important to define the country-specific SNP thresholds for investigation of transmission. This study demonstrated the use of thresholds of ≤5 and ≤15 SNPs difference between isolates to categorize definite and probable transmission, respectively. Different SNP thresholds should be applied while a higher cutoff was required to define an MDR-TB outbreak. The utilization of SNP thresholds proves to be crucial for guiding public health interventions, eliminating the need for unnecessary public health actions, and potentially uncovering undisclosed TB transmissions.

Unraveling tuberculosis patient cluster transmission chains: integrating WGS-based network with clinical and epidemiological insights

Background

Tuberculosis remains a global health threat, and the World Health Organization reports a limited reduction in disease incidence rates, including both new and relapse cases. Therefore, studies targeting tuberculosis transmission chains and recurrent episodes are crucial for developing the most effective control measures. Herein, multiple tuberculosis clusters were retrospectively investigated by integrating patients' epidemiological and clinical information with median-joining networks recreated based on whole genome sequencing (WGS) data of Mycobacterium tuberculosis isolates.

Methods

Epidemiologically linked tuberculosis patient clusters were identified during the source case investigation for pediatric tuberculosis patients. Only M. tuberculosis isolate DNA samples with previously determined spoligotypes identical within clusters were subjected to WGS and further median-joining network recreation. Relevant clinical and epidemiological data were obtained from patient medical records.

Results

We investigated 18 clusters comprising 100 active tuberculosis patients 29 of whom were children at the time of diagnosis; nine patients experienced recurrent episodes. M. tuberculosis isolates of studied clusters belonged to Lineages 2 (sub-lineage 2.2.1) and 4 (sub-lineages 4.3.3, 4.1.2.1, 4.8, and 4.2.1), while sub-lineage 4.3.3 (LAM) was the most abundant. Isolates of six clusters were drug-resistant. Within clusters, the maximum genetic distance between closely related isolates was only 5-11 single nucleotide variants (SNVs). Recreated median-joining networks, integrated with patients' diagnoses, specimen collection dates, sputum smear microscopy, and epidemiological investigation results indicated transmission directions within clusters and long periods of latent infection. It also facilitated the identification of potential infection sources for pediatric patients and recurrent active tuberculosis episodes refuting the reactivation possibility despite the small genetic distance of ≤5 SNVs between isolates. However, unidentified active tuberculosis cases within the cluster, the variable mycobacterial mutation rate in dormant and active states, and low M. tuberculosis genetic variability inferred precise transmission chain delineation. In some cases, heterozygous SNVs with an allelic frequency of 10-73% proved valuable in identifying direct transmission events.

Conclusion

The complex approach of integrating tuberculosis cluster WGS-data-based median-joining networks with relevant epidemiological and clinical data proved valuable in delineating epidemiologically linked patient transmission chains and deciphering causes of recurrent tuberculosis episodes within clusters.

Actionable mechanisms of drug tolerance and resistance in Mycobacterium tuberculosis

The emergence of antimicrobial resistance (AMR) across bacterial pathogens presents a serious threat to global health. This threat is further exacerbated in tuberculosis (TB), mainly due to a protracted treatment regimen involving a combination of drugs. A diversity of factors contributes to the emergence of drug resistance in TB, which is caused by the pathogen Mycobacterium tuberculosis (Mtb). While the traditional genetic mutation-driven drug resistance mechanisms operate in Mtb, there are also several additional unique features of drug resistance in this pathogen. Research in the past decade has enriched our understanding of such unconventional factors as efflux pumps, bacterial heterogeneity, metabolic states, and host microenvironment. Given that the discovery of new antibiotics is outpaced by the emergence of drug resistance patterns displayed by the pathogen, newer strategies for combating drug resistance are desperately needed. In the context of TB, such approaches include targeting the efflux capability of the pathogen, modulating the host environment to prevent bacterial drug tolerance, and activating the host anti-mycobacterial pathways. In this review, we discuss the traditional mechanisms of drug resistance in Mtb, newer understandings and the shaping of a set of unconventional approaches to target both the emergence and treatment of drug resistance in TB.

Large-scale Pan Genomic Analysis of Mycobacterium tuberculosis Reveals Key Insights Into Molecular Evolutionary Rate of Specific Processes and Functions

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), an infectious disease that is a major killer worldwide. Due to selection pressure caused by the use of antibacterial drugs, Mtb is characterised by mutational events that have given rise to multi drug resistant (MDR) and extensively drug resistant (XDR) phenotypes. The rate at which mutations occur is an important factor in the study of molecular evolution, and it helps understand gene evolution. Within the same species, different protein-coding genes evolve at different rates. To estimate the rates of molecular evolution of protein-coding genes, a commonly used parameter is the ratio dN/dS, where dN is the rate of non-synonymous substitutions and dS is the rate of synonymous substitutions. Here, we determined the estimated rates of molecular evolution of select biological processes and molecular functions across 264 strains of Mtb. We also investigated the molecular evolutionary rates of core genes of Mtb by computing the dN/dS values, and estimated the pan genome of the 264 strains of Mtb. Our results show that the cellular amino acid metabolic process and the kinase activity function evolve at a significantly higher rate, while the carbohydrate metabolic process evolves at a significantly lower rate for M. tuberculosis. These high rates of evolution correlate well with Mtb physiology and pathogenicity. We further propose that the core genome of M. tuberculosis likely experiences varying rates of molecular evolution which may drive an interplay between core genome and accessory genome during M. tuberculosis evolution.

Prolonged survival of a patient with active MDR-TB HIV co-morbidity: insights from a Mycobacterium tuberculosis strain with a unique genomic deletion

Coinfection of HIV and multidrug-resistant tuberculosis (MDR-TB) presents significant challenges in terms of the treatment and prognosis of tuberculosis, leading to complexities in managing the disease and impacting the overall outcome for TB patients. This study presents a remarkable case of a patient with MDR-TB and HIV coinfection who survived for over 8 years, despite poor treatment adherence and comorbidities. Whole genome sequencing (WGS) of the infecting Mycobacterium tuberculosis (Mtb) strain revealed a unique genomic deletion, spanning 18 genes, including key genes involved in hypoxia response, intracellular survival, immunodominant antigens, and dormancy. This deletion, that we have called "Del-X," potentially exerts a profound influence on the bacterial physiology and its virulence. Only few similar deletions were detected in other non-related Mtb genomes worldwide. In vivo evolution analysis identified drug resistance and metabolic adaptation mutations and their temporal dynamics during the patient's treatment course.

Genomic, phenotypic and demographic characterization of Mycobacterium tuberculosis in Israel in 2021

According to World Health Organization WHO, Tuberculosis (TB) is the second cause of death from infectious disease worldwide. During 2021, 10.6 million people were infected with TB, and 1.6 million people died. TB is caused by pathogens belonging to the Mycobacterium tuberculosis complex (MTBC), mainly by Mycobacterium tuberculosis (M.tb). Members of this complex are acid-fast bacilli, which can cause intrapulmonary and extra pulmonary TB, and can be divided into various lineages, based on genomic markers. The main public health threat comes from drug resistant M.tb strains, which are responsible for about 25% of TB death and treatment failure worldwide. Treating drug resistant TB patients significantly raises the costs of TB treatment. This study provides an overview of the demographic and drug susceptibility characteristics of newly diagnosed TB patients in Israel in 2021. The State of Israel has a very low level of TB endemicity and is at a pre-elimination phase. Notably, only 11.7% of the newly diagnosed TB patients were born in Israel. In this report, of the 154 new laboratory-confirmed TB patients, 66.7% had pulmonary TB, while 16% had extrapulmonary TB. Males accounted for 52% of the patients, with the most prevalent age group being 21-40. Most patients were citizens of Israel (53.9%), while 37.7% had no Israeli citizenship. Among non-citizens, there was a predominance of males and patients aged 21-40. The susceptibility profile showed a high resistance rate to streptomycin (18.2%) and to a lower extent to isoniazid (13.6%), pyrazinamide (8.4%), rifampicin (7.8%), and ethambutol (3.2%). Only 2 cases of XDR-TB and 10 MDR-TB strains were detected in Israel in 2021, with both XDR strains and 5 out of 10 MDR strains belonging to the Beijing lineage. Most of Beijing isolates were resistant to at least one tested drug. Genomic sequencing of 134 out of 156 strains and bioinformatics analysis using the MTBseq program and WHO mutation catalogue shows a good match with only 9 discrepancies between phenotypic and genotypic susceptibility profiles in first line drugs. The most common lineage is Delhi-Cas (23%) followed by the Beijing lineage (17%). Most patients from the Delhi-Cas lineage were born in Africa, while patients with Beijing isolates were born in different countries. Minimum spanning tree analysis identified 15 clusters. The study highlights the need for ongoing surveillance of TB using molecular and phenotypic tools to further decreasing the spreading level of the disease and develop effective treatment strategies.

The Mycobacterium tuberculosis genome at 25 years: lessons and lingering questions

First achieved in 1998 by Cole et al., the complete genome sequence of Mycobacterium tuberculosis continues to provide an invaluable resource to understand tuberculosis (TB), the leading cause of global infectious disease mortality. At the 25-year anniversary of this accomplishment, we describe how insights gleaned from the M. tuberculosis genome have led to vital tools for TB research, epidemiology, and clinical practice. The increasing accessibility of whole-genome sequencing across research and clinical settings has improved our ability to predict antibacterial susceptibility, to track epidemics at the level of individual outbreaks and wider historical trends, to query the efficacy of the bacille Calmette-Guérin (BCG) vaccine, and to uncover targets for novel antitubercular therapeutics. Likewise, we discuss several recent efforts to extract further discoveries from this powerful resource.

Resistance patterns and transmission of mono- and polyresistant TB: clinical impact of WGS

Objectives

Rapidly diagnosing drug-resistant TB is crucial for improving treatment and transmission control. WGS is becoming increasingly accessible and has added value to the diagnosis and treatment of TB. The aim of the study was to perform WGS to determine the rate of false-positive results of phenotypic drug susceptibility testing (pDST) and characterize the molecular mechanisms of resistance and transmission of mono- and polyresistant Mycobacterium (M.) tuberculosis.

Methods

WGS was performed on 53 monoresistant and 25 polyresistant M. tuberculosis isolates characterized by pDST. Sequencing data were bioinformatically processed to infer mutations encoding resistance and determine the origin of resistance and phylogenetic relationship between isolates studied.

Results

The data showed the variable sensitivity and specificity of WGS in comparison with pDST as the gold standard: isoniazid 92.7% and 92.3%; streptomycin 41.9% and 100.0%; pyrazinamide 15% and 94.8%; and ethambutol 75.0% and 98.6%, respectively. We found novel mutations encoding resistance to streptomycin (in gidB) and pyrazinamide (in kefB). Most isolates belonged to lineage 4 (80.1%) and the overall clustering rate was 11.5%. We observed lineage-specific gene variations encoding resistance to streptomycin and pyrazinamide.

Conclusions

This study highlights the clinical potential of WGS in ruling out false-positive drug resistance following phenotypic or genetic drug testing, and recommend this technology together with the WHO catalogue in designing an optimal individualized treatment regimen and preventing the development of MDR TB. Our results suggest that resistance is primarily developed through spontaneous mutations or selective pressure.

Characteristics of subtype III-A CRISPR-Cas system in Mycobacterium tuberculosis: An overview

CRISPR-Cas systems are the only RNA- guided adaptive immunity pathways that trigger the detection and destruction of invasive phages and plasmids in bacteria and archaea. Due to its prevalence and mystery, the Class 1 CRISPR-Cas system has lately been the subject of several studies. This review highlights the specificity of CRISPR-Cas system III-A in Mycobacterium tuberculosis, the tuberculosis-causing pathogen, for over twenty years. We discuss the difference between the several subtypes of Type III and their defence mechanisms. The anti-CRISPRs (Acrs) recently described, the critical role of Reverse transcriptase (RT) and housekeeping nuclease for type III CRISPR-Cas systems, and the use of this cutting-edge technology, its impact on the search for novel anti-tuberculosis drugs.

Exploration of the diversity of multi-drug resistant Mycobacterium tuberculosis complex in Lagos, Nigeria using WGS: Distribution of lineages, drug resistance patterns and genetic mutations

Multidrug-resistant (MDR) or extensively drug-resistant (XDR) Tuberculosis (TB) is a major challenge to global TB control. Therefore, accurate tracing of in-country MDR-TB transmission are crucial for the development of optimal TB management strategies. This study aimed to investigate the diversity of MTBC in Nigeria. The lineage and drug-resistance patterns of the clinical MTBC isolates of TB patients in Southwestern region of Nigeria were determined using the WGS approach. The phenotypic DST of the isolates was determined for nine anti-TB drugs. The sequencing achieved average genome coverage of 65.99X. The most represented lineages were L4 (n = 52, 83%), L1 (n = 8, 12%), L2 (n = 2, 3%) and L5 (n = 1, 2%), suggesting a diversified MTB population. In term of detection of M/XDR-TB, while mutations in katG and rpoB genes are the strong predictors for the presence of M/XDR-TB, the current study also found the lack of good genetic markers for drug resistance amongst the MTBC in Nigeria which may pose greater problems on local tuberculosis management efforts. This high-resolution molecular epidemiological data provides valuable insights into the mechanistic for M/XDR TB in Lagos, Nigeria.

Drug Resistance in Tuberculous Lymphadenitis: Molecular Characterization

Background

Drug-resistant tuberculosis (TB) epidemic in high-TB-incidence countries, particularly Ethiopia, remains a significant challenge. As a result, we investigated the drug resistance, common gene mutation, and molecular characterization of mycobacterial isolates from patients with suspected tuberculous lymphadenitis (TBLN). Methodology. A cross-sectional study of 218 FNA samples from TBLN patients inoculated on Lowenstein-Jensen media was carried out. The culture isolates were identified as MTB by polymerase chain reaction (PCR) and the difference-9 (RD9) test region. In addition, the GenoType MTBDRplus assay tested the first and second-line MTB drugs, and the spoligotyping strain-dependent polymorphism test was determined.

Results

Among the 50 culture-positive isolates, 14% (7/50) had drug resistance caused by a gene mutation. Out of these, 4 (8%) isolates were mono-resistant to isoniazid drug, which is caused by a gene mutation in katG in the region of interrogated at codon 315 in the amino acid sequence of S315T1, and 3 (6%) isolates were resistant to both rifampicin and isoniazid drugs. The mutation was observed for katG (at codon 315 with a change in the sequence of amino acid S315T) and rpoB (at codon 530-533 with a change in the sequence of amino acid S531L (S450L)) genes. The most prevalent spoligotypes were orphan and SIT53 strains.

Conclusion

The predominance of INH mono-resistance poses a critical risk for the potential development of MDR-TB, as INH mono-resistance is a typical pathway to the occurrence of MDR-TB. The orphan and SIT53 (T) strains were the most common in the study area, and a drug-resistant strain caused by a common gene mutation could indicate the transmission of clonal-resistant strains in the community.

Whole-genome sequencing-based analyses of drug-resistant Mycobacterium tuberculosis from Taiwan

Drug-resistant tuberculosis (DR-TB) posed challenges to global TB control. Whole-genome sequencing (WGS) is recommended for predicting drug resistance to guide DR-TB treatment and management. Nevertheless, data are lacking in Taiwan. Phenotypic drug susceptibility testing (DST) of 12 anti-TB drugs was performed for 200 Mycobacterium tuberculosis isolates. WGS was performed using the Illumina platform. Drug resistance profiles and lineages were predicted in silico using the Total Genotyping Solution for TB (TGS-TB). Using the phenotypic DST results as a reference, WGS-based prediction demonstrated high concordance rates of isoniazid (95.0%), rifampicin (RIF) (98.0%), pyrazinamide (98.5%) and fluoroquinolones (FQs) (99.5%) and 96.0% to 99.5% for second-line injectable drugs (SLIDs); whereas, lower concordance rates of ethambutol (87.5%), streptomycin (88.0%) and ethionamide (84.0%). Furthermore, minimum inhibitory concentrations confirmed that RIF rpoB S450L, FQs gyrA D94G and SLIDs rrs a1401g conferred high resistance levels. Besides, we identified lineage-associated mutations in lineage 1 (rpoB H445Y and fabG1 c-15t) and predominant lineage 2 (rpoB S450L and rpsL K43R). The WGS-based prediction of drug resistance is highly concordant with phenotypic DST results and can provide comprehensive genetic information to guide DR-TB precision therapies in Taiwan.

Genetic Diversity and Primary Drug Resistance of Mycobacterium tuberculosis Beijing Genotype Strains in Northwestern Russia

The Beijing genotype is the main family of Mycobacterium tuberculosis in Russia. We analyzed its diversity and drug resistance in provinces across Northwestern Russia to identify the epidemiologically relevant Beijing strains. The study collection included 497 isolates from newly-diagnosed tuberculosis (TB) patients. Bacterial isolates were subjected to drug-susceptibility testing and genotyping. The Beijing genotype was detected in 57.5% (286/497); 50% of the Beijing strains were multidrug-resistant (MDR). Central Asian/Russian and B0/W148 groups included 176 and 77 isolates, respectively. MDR was more frequent among B0/W148 strains compared to Central Asian/Russian strains (85.7% vs. 40.3%, p < 0.0001). Typing of 24 minisatellite loci of Beijing strains revealed 82 profiles; 230 isolates were in 23 clusters. The largest Central Asian/Russian types were 94-32 (n = 75), 1065-32 (n = 17), and 95-32 (n = 12). B0/W148 types were 100-32 (n = 59) and 4737-32 (n = 5). MDR was more frequent in types 1065-32 (88.2%), 100-32 (83.1%), and 4737-32 (100%). In contrast, type 9391-32 (n = 9) included only drug-susceptible strains. To conclude, M. tuberculosis Beijing genotype is dominant in Northwestern Russia, and an active transmission of overwhelmingly MDR B0/W148 types explains the reported increase of MDR-TB. The presence of MDR-associated minor variants (type 1071-32/ancient Beijing and Central Asia Outbreak strain) in some of the studied provinces also requires attention.

Evaluation of Five User-Friendly Whole Genome Sequencing Software for Mycobacterium tuberculosis in Clinical Application

BACKGROUND

Whole genome sequencing (WGS) is an increasingly useful tool for tuberculosis (TB) diagnosis and disease management. In this study, we evaluated the utility of user-friendly WGS tools in reporting resistance profiles and identifying lineages of clinical TB isolates from South Korea.

METHODS

Forty clinical samples from TB patients showing discrepancies between their rapid molecular and conventional drug susceptibility tests were used in this study. Among these clinical isolates, 37 strains were successfully evaluated via WGS software, using the GenTB, TB Profiler, PhyResSE, CASTB, and Mykrobe.

RESULTS

More accurate and faster susceptibility results could be obtained with isoniazid than with rifampin. Using the phenotypic test as the gold standard, the isoniazid concordance rate between phenotypic drug susceptibility test (DST) and WGS (GenTB: 45.9%, TB profiler: 40.5%, PhyResSE: 40.5%, CASTB: 48.6%, and Mykrobe: 43.2%) was much higher than between phenotypic DST and rapid molecular genotypic DST (18.9%) among the 37 strains. In contrast, the rifampin concordance rate between phenotypic DST and WGS and that between phenotypic DST and rapid molecular genotypic DST was similar (81.1-89.2%). We also found novel mutations associated with INH in katG and ahpC gene region, not covered by the line probe assay. In addition, lineage analysis identified 81.1% of these samples as L2 East Asian lineage strains, and 18.9% as L4 Euro-American lineage strains.

CONCLUSION

WGS may play a pivotal role in TB diagnosis and the detection of drug resistance, genetic diversity, and transmission dynamics in the near future because of its accuracy, speed, and extensibility.

Evolution of Mycobacterium tuberculosis drug resistance in the genomic era

Mycobacterium tuberculosis has acquired drug resistance to all drugs that have been used against it, including those only recently introduced into clinical practice. Compared to other bacteria, it has a well conserved genome due to its role as an obligate human pathogen that has adapted to a niche over five to ten thousand years. These features facilitate reconstruction and dating of M. tuberculosis phylogenies, giving key insights into how resistance has been acquired and spread globally. Resistance to each new drug has occurred within five to ten years of clinical use and has occurred even more rapidly with recently introduced drugs. In most cases, resistance-conferring mutations come with a fitness cost, but this can be overcome by compensatory mutations which restore fitness to that of wild-type bacteria. It is likely that M. tuberculosis acquires drug resistance while maintaining limited genomic variability due the generation of low frequency within-host variation, combined with ongoing purifying selection causing loss of variants without a clear fitness advantage. However, variants that do confer an advantage, such as drug resistance, can increase in prevalence amongst all bacteria within a host and become the dominant clone. These resistant strains can then be transmitted leading to primary drug resistant infection in a new host. As many countries move towards genomic methods for diagnosis of M. tuberculosis infection and drug resistance, it is important to be aware of the implications for the evolution of resistance. Currently, understanding of resistance-conferring mutations is incomplete, and some targeted genetic diagnostics create their own selective pressures. We discuss an example where a rifampicin resistance-conferring mutation which was not routinely covered by standard testing became dominant. Finally, resistance to new drugs such as bedaquiline and delamanid is caused by individually rare mutations occurring across a large mutational genomic target that have been detected over a short time, and do not provide statistical power for genotype-phenotype correlation – in contrast to longer-established drugs that form the backbone of drug-sensitive antituberculosis therapy. Therefore, we need a different approach to identify resistance-conferring mutations of new drugs before their resistance becomes widespread, abrogating their usefulness.

Whole genome sequencing analysis to evaluate the influence of T2DM on polymorphisms associated with drug resistance in M. tuberculosis

Background

Type 2 diabetes mellitus (T2DM) has been associated with treatment failure, and the development of drug resistance in tuberculosis (TB). Also, whole-genome sequencing has provided a better understanding and allowed the growth of knowledge about polymorphisms in genes associated with drug resistance. Considering the above, this study analyzes genome sequences to evaluate the influence of type 2 diabetes mellitus in the development of mutations related to tuberculosis drug resistance. M. tuberculosis isolates from individuals with (n = 74), and without (n = 74) type 2 diabetes mellitus was recovered from online repositories, and further analyzed.

Results

The results showed the presence of 431 SNPs with similar proportions between diabetics, and non-diabetics individuals (48% vs. 52%), but with no significant relationship. A greater number of mutations associated with rifampicin resistance was observed in the T2DM-TB individuals (23.2% vs. 16%), and the exclusive presence of rpoBQ432L, rpoBQ432P, rpoBS441L, and rpoBH445L variants. While these variants are not private to T2DM-TB cases they are globally rare highlighting a potential role of T2DM. The phylogenetic analysis showed 12 sublineages, being 4.1.1.3, and 4.1.2.1 the most prevalent in T2DM-TB individuals but not differing from those most prevalent in their geographic location. Four clonal complexes were found, however, no significant relationship with T2DM was observed. Samples size and potential sampling biases prevented us to look for significant associations.

Conclusions

The occurrence of globally rare rifampicin variants identified only in isolates from individuals with T2DM could be due to the hyperglycemic environment within the host. Therefore, further studies about the dynamics of SNPs’ generation associated with antibiotic resistance in patients with diabetes mellitus are necessary.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08709-z.

Spoligotype Variation of Mycobacterium tuberculosis Strains Prevailing in Korea

Tuberculosis (TB) is an ongoing global health problem, including in South Korea. To manage TB efficiently, it is necessary to understand the epidemiology, transmission route, and characteristics of prevailing Mycobacterium tuberculosis strains. In this study, we investigated microevolutions over time in the spoligotype patterns of M. tuberculosis isolated from TB patients in Korea. We collected 1,055 clinical M. tuberculosis isolates from 16 provinces in Korea from 1994 to 2006 and analyzed them by spoligotyping. We observed 26 subfamilies, including two large predominant families: a Beijing family (72.7%) and the T family (19.1%). Specifically, the abundance of spoligotype SIT269 from the Beijing-like subfamily significantly increased in the 2000s relative to the 1990s in Korea. This study provides an overview of the M. tuberculosis genotype trends over time in Korea. These data also indicate that we should consider the influence of the newly growing SIT269 subtype identified in the Beijing family.

Comparative Genomic Analysis of Mycobacteriaceae Reveals Horizontal Gene Transfer-Mediated Evolution of the CRISPR-Cas System in the Mycobacterium tuberculosis Complex

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) genes are conserved genetic elements in many prokaryotes, including Mycobacterium tuberculosis, the causative agent of tuberculosis. Although knowledge of CRISPR locus variability has been utilized in M. tuberculosis strain genotyping, its evolutionary path in Mycobacteriaceae is not well understood. In this study, we have performed a comparative analysis of 141 mycobacterial genomes and identified the exclusive presence of the CRISPR-Cas type III-A system in M. tuberculosis complex (MTBC). Our global phylogenetic analysis of CRISPR repeats and Cas10 proteins offers evidence of horizontal gene transfer (HGT) of the CRISPR-Cas module in the last common ancestor of MTBC and Mycobacterium canettii from a Streptococcus-like environmental bacterium. Additionally, our results show that the variation of CRISPR-Cas organization in M. tuberculosis lineages, especially in the Beijing sublineage of lineage 2, is due to the transposition of insertion sequence IS6110 The direct repeat (DR) region of the CRISPR-Cas locus acts as a hot spot for IS6110 insertion. We show in M. tuberculosis H37Rv that the repeat at the 5' end of CRISPR1 of the forward strand is an atypical repeat made up partly of IS-terminal inverted repeat and partly CRISPR DR. By tracing an undetectable spacer sequence in the DR region, the two CRISPR loci could theoretically be joined to reconstruct the ancestral single CRISPR-Cas locus organization, as seen in M. canettii This study retracing the evolutionary events of HGT and IS6110-driven genomic deletions helps us to better understand the strain-specific variations in M. tuberculosis lineages.IMPORTANCE Comparative genomic analysis of prokaryotes has led to a better understanding of the biology of several pathogenic microorganisms. One such clinically important pathogen is M. tuberculosis, the leading cause of bacterial infection worldwide. Recent evidence on the functionality of the CRISPR-Cas system in M. tuberculosis has brought back focus on these conserved genetic elements, present in many prokaryotes. Our study advances understanding of mycobacterial CRISPR-Cas origin and its diversity among the different species. We provide phylogenetic evidence of acquisition of CRISPR-Cas type III-A in the last common ancestor shared between MTBC and M. canettii, by HGT-mediated events. The most likely source of HGT was an environmental Firmicutes bacterium. Genomic mapping of the CRISPR loci showed the IS6110 transposition-driven variations in M. tuberculosis strains. Thus, this study offers insights into events related to the evolution of CRISPR-Cas in M. tuberculosis lineages.