A Major Mycobacterium tuberculosis outbreak caused by one specific genotype in a low-incidence country: Exploring gene profile virulence explanations

Considering best practice standards for routine whole-genome sequencing for TB care and control

TB is a priority pathogen for the application of whole-genome sequencing (WGS) into routine public health practice. In low-incidence settings, a growing number of services have begun to incorporate routine WGS into standard practice. The increasing availability of real-time genomic information supports a variety of aspects of the public health response, including the detection of drug resistance, monitoring of laboratory and clinical practices, contact tracing investigations and active case finding. Optimal structures and approaches are needed to support the rapid translation of genomic information into practice and to evaluate outcomes and impact. In this consensus paper, we outline the elements needed to systemically incorporate routine WGS into the TB public health response, including the sustainability of services, multidisciplinary team models and monitoring and evaluation frameworks. If integrated in an efficient and thoughtful manner, routine WGS has the potential to significantly improve clinical TB care for individuals and the overall public health response.

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.

The occurrence rate of Haarlem and Beijing genotypes among Middle Eastern isolates of multi drug resistant Mycobacterium tuberculosis: A systematic review and meta-analysis

Antibiotic resistance is a serious problem that poses a major challenge to tuberculosis control worldwide. Many developing countries still struggle with this infection in term of various aspects as it remains a major health concern. A number of developing countries are located in the Middle East, one of the world's most important regions. The control of this infection remains largely suboptimal despite intensive research in the field, and the mechanisms that lead to its progression have not yet been fully understood. Therefore, TB control must be amended through the identification of new strategies. For this reason, monitoring genetic characterizations of TB strains by molecular typing methods in different geographical regions can be important to setting local programs and global strategies to control TB infection. It is important to know the genotype of Mycobacterium tuberculosis strains to evaluate the occurrence of outbreaks and the transmission of this disease. Beijing and Haarlem genotypes are the most prevalent and, in these families, there is greater association with drug resistance, resulting in more severe forms of TB and higher levels of treatment failure than in other families. The current study is planned to systematically conduct a review using a meta-analysis to show the prevalence of Beijing and Haarlem genotypes in the Middle Eastern MDR-TB cases. M. tuberculosis strains pose particular epidemiological and clinical concerns as they can endanger tuberculosis control programs.

Deciphering the Tangible Spatio-Temporal Spread of a 25-Year Tuberculosis Outbreak Boosted by Social Determinants

Outbreak strains of Mycobacterium tuberculosis are promising candidates as targets in the search for intrinsic determinants of transmissibility, as they are responsible for many cases with sustained transmission; however, the use of low-resolution typing methods and restricted geographical investigations represent flaws in assessing the success of long-lived outbreak strains. We can now address the nature of outbreak strains by combining large genomic data sets and phylodynamic approaches. We retrospectively sequenced the whole genome of representative samples assigned to an outbreak circulating in the Canary Islands (the GC strain) since 1993, which accounts for ~20% of local tuberculosis cases. We selected a panel of specific single nucleotide polymorphism (SNP) markers for an in-silico search for additional outbreak-related sequences within publicly available tuberculosis genomic data. Using this information, we inferred the origin, spread, and epidemiological parameters of the GC strain. Our approach allowed us to accurately trace the historical and more recent dispersion of the GC strain. We provide evidence of a highly successful nature within the Canarian archipelago but limited expansion abroad. Estimation of epidemiological parameters from genomic data disagree with a distinctive biology of the GC strain. With the increasing availability of genomic data allowing for the accurate inference of strain spread and critical epidemiological parameters, we can now revisit the link between Mycobacterium tuberculosis genotypes and transmission, as is routinely carried out for SARS-CoV-2 variants of concern. We demonstrate that social determinants rather than intrinsically higher bacterial transmissibility better explain the success of the GC strain. Importantly, our approach can be used to trace and characterize strains of interest worldwide. IMPORTANCE Infectious disease outbreaks represent a significant problem for public health. Tracing outbreak expansion and understanding the main factors behind emergence and persistence remain critical to effective disease control. Our study allows researchers and public health authorities to use Whole-Genome Sequencing-based methods to trace outbreaks, and shows how available epidemiological information helps to evaluate the factors underpinning outbreak persistence. Taking advantage of all the freely available information placed in public repositories, researchers can accurately establish the expansion of an outbreak beyond original boundaries, and determine the potential risk of a strain to inform health authorities which, in turn, can define target strategies to mitigate expansion and persistence. Finally, we show the need to evaluate strain transmissibility in different geographic contexts to unequivocally associate spread to local or pathogenic factors, an important lesson taken from genomic surveillance of SARS-CoV-2.

Pan-genome association study of Mycobacterium tuberculosis lineage-4 revealed specific genes related to the high and low prevalence of the disease in patients from the North-Eastern area of Medellín, Colombia

Mycobacterium tuberculosis (Mtb) lineage 4 is responsible for the highest burden of tuberculosis (TB) worldwide. This lineage has been the most prevalent lineage in Colombia, especially in the North-Eastern (NE) area of Medellin, where it has been shown to have a high prevalence of LAM9 SIT42 and Haarlem1 SIT62 sublineages. There is evidence that regardless of environmental factors and host genetics, differences among sublineages of Mtb strains play an important role in the course of infection and disease. Nevertheless, the genetic basis of the success of a sublineage in a specific geographic area remains uncertain. We used a pan-genome-wide association study (pan-GWAS) of 47 Mtb strains isolated from NE Medellin between 2005 and 2008 to identify the genes responsible for the phenotypic differences among high and low prevalence sublineages. Our results allowed the identification of 12 variants in 11 genes, of which 4 genes showed the strongest association to low prevalence (mmpL12, PPE29, Rv1419, and Rv1762c). The first three have been described as necessary for invasion and intracellular survival. Polymorphisms identified in low prevalence isolates may suggest related to a fitness cost of Mtb, which might reflect a decrease in their capacity to be transmitted or to cause an active infection. These results contribute to understanding the success of some sublineages of lineage-4 in a specific geographical area.

Systems biology approach to functionally assess the Clostridioides difficile pangenome reveals genetic diversity with discriminatory power

SignificanceClostridioides difficile infections are the most common source of hospital-acquired infections and are responsible for an extensive burden on the health care system. Strains of the C. difficile species comprise diverse lineages and demonstrate genome variability, with advantageous trait acquisition driving the emergence of endemic lineages. Here, we present a systems biology analysis of C. difficile that evaluates strain-specific genotypes and phenotypes to investigate the overall diversity of the species. We develop a strain typing method based on similarity of accessory genomes to identify and contextualize genetic loci capable of discriminating between strain groups.

Impact of pathobiological diversity of Mycobacterium tuberculosis on clinical features and lethal outcome of tuberculosis

BACKGROUND

Mycobacterium tuberculosis population in Russia is dominated by the notorious Beijing genotype whose major variants are characterized by contrasting resistance and virulence properties. Here we studied how these strain features could impact the progression of pulmonary tuberculosis (TB) concerning clinical manifestation and lethal outcome.

RESULTS

The study sample included 548 M. tuberculosis isolates from 548 patients with newly diagnosed pulmonary TB in Omsk, West Siberia, Russia. Strains were subjected to drug susceptibility testing and genotyping to detect lineages, sublineages, and subtypes (within Beijing genotype). The Beijing genotype was detected in 370 (67.5%) of the studied strains. The strongest association with multidrug resistance (MDR) was found for epidemic cluster Beijing B0/W148 (modern sublineage) and two recently discovered MDR clusters 1071-32 and 14717-15 of the ancient Beijing sublineage. The group of patients infected with hypervirulent and highly lethal (in a mouse model) Beijing 14717-15 showed the highest rate of lethal outcome (58.3%) compared to Beijing B0/W148 (31.4%; P = 0.06), Beijing Central Asian/Russian (29.7%, P = 0.037), and non-Beijing (15.2%, P = 0.001). The 14717-15 cluster mostly included isolates from patients with infiltrative but not with fibrous-cavernous and disseminated TB. In contrast, a group infected with low virulent 1071-32-cluster had the highest rate of fibrous-cavernous TB, possibly reflecting the capacity of these strains for prolonged survival and chronicity of the TB process.

CONCLUSIONS

The group of patients infected with hypervirulent and highly lethal in murine model 14717-15 cluster had the highest proportion of the lethal outcome (58.3%) compared to the groups infected with Beijing B0/W148 (31.4%) and non-Beijing (15.2%) isolates. This study carried out in the TB high-burden area highlights that not only drug resistance but also strain virulence should be considered in the implementation of personalized TB treatment.

Tuberculosis Drug Susceptibility, Treatment, and Outcomes for Belarusian HIV-Positive Patients with Tuberculosis: Results from a National and International Laboratory

BACKGROUND

To cure drug-resistant (DR) tuberculosis (TB), the antituberculous treatment should be guided by Mycobacterium tuberculosis drug-susceptibility testing (DST). In this study, we compared conventional DST performed in Minsk, Belarus, a TB DR high-burden country, with extensive geno- and phenotypic analyses performed at the WHO TB Supranational Reference Laboratory in Copenhagen, Denmark, for TB/HIV coinfected patients. Subsequently, DST results were related to treatment regimen and outcome.

METHODS

Thirty TB/HIV coinfected patients from Minsk were included and descriptive statistics applied.

RESULTS

Based on results from Minsk, 10 (33%) TB/HIV patients had drug-sensitive TB. Two (7%) had isoniazid monoresistant TB, 8 (27%) had multidrug-resistant (MDR) TB, 5 (17%) preextensive drug-resistant (preXDR) TB, and 5 (17%) had extensive drug-resistant (XDR) TB. For the first-line drugs rifampicin and isoniazid, there was DST agreement between Minsk and Copenhagen for 90% patients. For the second-line anti-TB drugs, discrepancies were more pronounced. For 14 (47%) patients, there were disagreements for at least one drug, and 4 (13%) patients were classified as having MDR-TB in Minsk but were classified as having preXDR-TB based on DST results in Copenhagen. Initially, all patients received standard anti-TB treatment with rifampicin, isoniazid, pyrazinamide, and ethambutol. However, this was only suitable for 40% of the patients based on DST. On average, DR-TB patients were changed to 4 (IQR 3-5) active drugs after 1.5 months (IQR 1-2). After treatment adjustment, the treatment duration was 8 months (IQR 2-11). Four (22%) patients with DR-TB received treatment for >18 months. In total, sixteen (53%) patients died during 24 months of follow-up.

CONCLUSIONS

We found high concordance for rifampicin and isoniazid DST between the Minsk and Copenhagen laboratories, whereas discrepancies for second-line drugs were more pronounced. For patients with DR-TB, treatment was often insufficient and relevant adjustments delayed. This example from Minsk, Belarus, underlines two crucial points in the management of DR-TB: the urgent need for implementation of rapid molecular DSTs and availability of second-line drugs in all DR-TB high-burden settings. Carefully designed individualized treatment regimens in accordance with DST patterns will likely improve patients' outcome and reduce transmission with drug-resistant Mycobacterium tuberculosis strains.

Recurrences of multidrug-resistant tuberculosis: Strains involved, within-host diversity, and fine-tuned allocation of reinfections

Recurrent tuberculosis occurs due to exogenous reinfection or reactivation/persistence. We analysed 90 sequential MDR Mtb isolates obtained in Argentina from 27 patients with previously diagnosed MDR-TB that recurred in 2018 (1-10 years, 2-10 isolates per patient). Three long-term predominant strains were responsible for 63% of all MDR-TB recurrences. Most of the remaining patients were infected by strains different from each other. Reactivation/persistence of the same strain caused all but one recurrence, which was due to a reinfection with a predominant strain. One of the prevalent strains showed marked stability in the recurrences, while in another strain higher SNP-based diversity was observed. Comparisons of intra- versus inter-patient SNP distances identified two possible reinfections with closely related variants circulating in the community. Our results show a complex scenario of MDR-TB infections in settings with predominant MDR Mtb strains.

An Overview of Genetic Information of Latent Mycobacterium tuberculosis

Mycobacterium tuberculosis has infected more than two billion individuals worldwide, of whom 5%–10% have clinically active disease and 90%–95% remain in the latent stage with a reservoir of viable bacteria in the macrophages for extended periods of time. The tubercle bacilli at this stage are usually called dormant, non-viable, and/or non-culturable microorganisms. The patients with latent bacilli will not have clinical pictures and are not infectious. The infections in about 2%–23% of the patients with latent status become reactivated for various reasons such as cancer, human immunodeficiency virus infection, diabetes, and/or aging. Many studies have examined the mechanisms involved in the latent state of Mycobacterium and showed that latency modified the expression of many genes. Therefore, several mechanisms will change in this bacterium. Hence, this study aimed to briefly examine the genes involved in the latent state as well as the changes that are caused by Mycobacterium tuberculosis. The study also evaluated the relationship between the functions of these genes.

Molecular epidemiology of leprosy: An update

Molecular epidemiology investigations are notoriously challenging in the leprosy field mainly because the inherent characteristics of the disease as well as its yet uncultivated causative agents, Mycobacterium leprae and M. lepromatosis. Despite significant developments in understanding the biology of leprosy bacilli through genomic approaches, the exact mechanisms of transmission is still unclear and the factors underlying pathological variation of the disease in different patients remain as major gaps in our knowledge about leprosy. Despite these difficulties, the last two decades have seen the development of genotyping procedures based on PCR-sequencing of target loci as well as by the genome-wide analysis of an increasing number of geographically diverse isolates of leprosy bacilli. This has provided a foundation for molecular epidemiology studies that are bringing a better understanding of strain evolution associated with ancient human migrations, and phylogeographical insights about the spread of disease globally. This review discusses the advantages and drawbacks of the main tools available for molecular epidemiological investigations of leprosy and summarizes various methods ranging from PCR-based genotyping to genome-typing techniques. We also describe their main applications in analyzing the short-range and long-range transmission of the disease. Finally, we summarise the current gaps and challenges that remain in the field of molecular epidemiology of leprosy.

Genomic signatures of drug resistance in highly resistant Mycobacterium tuberculosis strains of the early ancient sublineage of Beijing genotype in Russia

The Mycobacterium tuberculosis Beijing genotype is a clinically and epidemiologically important lineage that is subdivided into ancient/ancestral and modern strains. In our previous study in western Siberia, we identified variable number of tandem repeats (VNTR)-based clusters within the early ancient sublineage of the Beijing genotype characterized by an unexpectedly high rate of extensive drug resistance (XDR). In the current study, next generation sequencing data were analysed to gain insight into genomic signatures underlying drug resistance of these strains. A total of 184 genomes of the Beijing early ancient sublineage from Russia (16), China (15), Japan (36), Korea (25), Vietnam (18), Thailand (73), and the USA (1) were used for phylogenetic analysis. The drug-resistant profile was deduced genotypically. The Russian isolates were distributed into two clusters and were all drug resistant, mainly pre-XDR and XDR. The largest of these clusters included only Russian isolates from remote locations in both Asian and European parts of the country. All its isolates had a quadruple drug resistance (to isoniazid, rifampin, ethambutol and streptomycin) due to the 6-mutation signature (KatG Ser315Thr, KatG Ile335Val, RpoB Ser450Leu, RpoC Asp485Asn, EmbB Gln497Arg, and RpsL Lys43Arg). In most samples, it was complemented with additional and different pncA, gyrA and rrs mutations leading to the pre-XDR/XDR genotype. Phylogenomic analysis indicates a distant origin of this Russian resistant cluster in the early 1970s but location and circumstances are yet to be clarified.

Recurrent tuberculosis in patients infected with the predominant Mycobacterium tuberculosis outbreak strain in Denmark. New insights gained through whole genome sequencing

OBJECTIVES

Recurrent tuberculosis (TB) is defined by more than one TB episode per patient and is caused by reinfection with a new M. tuberculosis (Mtb) strain or relapse with the previous strain. In Denmark, a major TB outbreak caused by one specific Mtb genotype "DKC2" is ongoing. Of the 892 patients infected with DKC2 between 1992 and 2014, 32 had recurrent TB with 67 TB episodes in total.

METHODS

The 32 cases were evaluated in terms of number of single-nucleotide polymorphism (SNP) differences and time between episodes derived from whole-genome sequencing data.

RESULTS

For four TB cases, the subsequent episodes could be confirmed as relapse and for one case as reinfection. Eight cases with SNP differences <6, theoretically indicating relapse, could be classified as likely reinfections based on phylogenetic analysis in combination with geographical data. Subsequent TB episodes for the remaining 19 cases could not be classified as relapse or reinfection even though they all had a SNP difference of <6 SNPs.

CONCLUSIONS

In newer studies, investigating recurrent TB with the use of WGS, the number of SNPs has been used to distinguish between relapse and reinfection. The algorithm proposed for this is not valid in the Danish TB outbreak setting as our findings challenge the interpretation of few SNP differences as representing relapse. However, when including phylogenetic analysis and geographical data in the analysis, classification of 13 of the 32 cases were possible.

Current topics of molecular mycobacteriology

The 2nd St. Petersburg Symposium on Tuberculosis and Mycobacteria: Molecular Approach, was held in St. Petersburg, Russia on 5-6 December 2018. A special issue of Infection, Genetics and Evolution will publish articles based on the selected presentations. In this paper, I will discuss some of the hot topics of molecular mycobacteriology highlighted at this meeting that I had the pleasure to organize and honor to chair. The symposium addressed interrelated fundamental and applied issues of modern mycobacteriology such as molecular evolution and phylogenomics, host-microbe interactions and pathogenesis, coevolution of M. tuberculosis with humans, new genomic and postgenomic technologies. Molecular methods for TB diagnostics and drug resistance detection are supported by WHO and whole genome/next generation sequencing presents a comprehensive approach. At the same time, cost and implementation of new methods for direct analysis of clinical samples and/or in low-resource settings remain a great challenge. A due attention was also given to the medically important nontuberculous mycobacteria. Assessment of spectrum of the circulating mycobacterial species in the Russian Federation and the countries of the European Union was presented and the underlying reasons of the observed diversity were discussed. To conclude, the symposium became a multidisciplinary event that was useful to promote networking and exchange of knowledge and experience. The next (third) symposium was planned to be organized in 2021.

Complete Genome Sequence of Mycobacterium tuberculosis DKC2, the Predominant Danish Outbreak Strain

The largest clonal outbreak of Mycobacterium tuberculosis infection in Scandinavia has been monitored by the International Reference Laboratory of Mycobacteriology (IRLM) since 1992. Here, we present the complete genome sequence of M. tuberculosis strain DKC2 substrain PP1, a representative isolate collected in 1993 from a Danish patient with pulmonary tuberculosis.

The largest clonal outbreak of Mycobacterium tuberculosis infection in Scandinavia has been monitored by the International Reference Laboratory of Mycobacteriology (IRLM) since 1992. Here, we present the complete genome sequence of M. tuberculosis strain DKC2 substrain PP1, a representative isolate collected in 1993 from a Danish patient with pulmonary tuberculosis.