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

Pangenome and genomic signatures linked to the dominance of the lineage-4 of Mycobacterium tuberculosis isolated from extrapulmonary tuberculosis patients in western Ethiopia

BACKGROUND

The lineage 4 (L4) of Mycobacterium tuberculosis (MTB) is not only globally prevalent but also locally dominant, surpassing other lineages, with lineage 2 (L2) following in prevalence. Despite its widespread occurrence, factors influencing the expansion of L4 and its sub-lineages remain poorly understood both at local and global levels. Therefore, this study aimed to conduct a pan-genome and identify genomic signatures linked to the elevated prevalence of L4 sublineages among extrapulmonary TB (EPTB) patients in western Ethiopia.

METHODS

A cross-sectional study was conducted at an institutional level involving confirmed cases of extrapulmonary tuberculosis (EPTB) patients from August 5, 2018, to December 30, 2019. A total of 75 MTB genomes, classified under lineage 4 (L4), were used for conducting pan-genome and genome-wide association study (GWAS) analyses. After a quality check, variants were identified using MTBseq, and genomes were de novo assembled using SPAdes. Gene prediction and annotation were performed using Prokka. The pan-genome was constructed using GET_HOMOLOGUES, and its functional analysis was carried out with the Bacterial Pan-Genome Analysis tool (BPGA). For GWAS analysis, Scoary was employed with Benjamini-Hochberg correction, with a significance threshold set at p-value ≤ 0.05.

RESULTS

The analysis revealed a total of 3,270 core genes, predominantly associated with orthologous groups (COG) functions, notably in the categories of '[R] General function prediction only' and '[I] Lipid transport and metabolism'. Conversely, functions related to '[N] Cell motility' and '[Q] Secondary metabolites biosynthesis, transport, and catabolism' were primarily linked to unique and accessory genes. The pan-genome of MTB L4 was found to be open. Furthermore, the GWAS study identified genomic signatures linked to the prevalence of sublineages L4.6.3 and L4.2.2.2.

CONCLUSIONS

Apart from host and environmental factors, the sublineage of L4 employs distinct virulence factors for successful dissemination in western Ethiopia. Given that the functions of these newly identified genes are not well understood, it is advisable to experimentally validate their roles, particularly in the successful transmission of specific L4 sublineages over others.

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.

A first insight into tuberculosis transmission at the border of Ecuador and Colombia: a retrospective study of the population structure of Mycobacterium tuberculosis in Esmeraldas province

Objective

Tuberculosis (TB) is a major public health concern in Ecuador and Colombia, considering that both countries are high-burden TB settings. Molecular epidemiology is crucial to understand the transmission dynamics of Mycobacterium tuberculosis complex (MTBC) and to identify active transmission clusters of regional importance.

Methods

We studied the potential transmission of TB between Colombia and Ecuador through the analysis of the population structure of MTBC lineages circulating in the Ecuadorian province of Esmeraldas at the border with Colombia. A total of 105 MTBC strains were characterized by 24-loci MIRU-VNTR and spoligotyping.

Results

MTBC lineage 4 is only present in Esmeraldas; no MTBC strains belonging to Lineage 2-sublineage Beijing were found despite its presence in other provinces of Ecuador and, in Colombia. Genotyping results revealed a high degree of diversity for MTBC in Esmeraldas: Neither active transmission clusters within this province nor including MTBC strains from Colombia or other provinces of Ecuador were found.

Conclusion

Our data suggest that tuberculosis dynamics in this rural and isolated area may be not related to highly transmitted strains but could be influenced by other health determinants that favor TB relapse such as poverty and poor health system access. Further studies including a larger number of MTBC strains from Esmeraldas are necessary to test this hypothesis.

Mycobacterium tuberculosis complex molecular networks and their regulation: Implications of strain heterogeneity on epigenetic diversity and transcriptome regulation

Tuberculosis has been a public health crisis since the 1900, which has caused the highest mortalities due to a single bacterial infection worldwide, that was recently further complicated by the Coronavirus disease 2019 pandemic. The causative agent of Tuberculosis, Mycobacterium tuberculosis, belongs to a genetically well-characterized family of strains known as the Mycobacterium tuberculosis complex, which has complicated progress made towards eradicating Tuberculosis due to pathogen-specific phenotypic differences in the members of this complex. Mycobacterium tuberculosis complex strains are genetically diverse human- and animal-adapted pathogens belonging to 7 lineages (Indo-Oceanic, East-Asian, East-African Indian, Euro-American, M. africanum West Africa 1, M. africanum West Africa 2 and Ethopia), respectively and the recently identified Lineage 8 and M. africanum Lineage 9. Genomic studies have revealed that Mycobacterium tuberculosis complex members are ∼99 % similar, however, due to selective pressure and adaptation to human host, they are prone to mutations that have resulted in development of drug resistance and phenotypic heterogeneity that impact strain virulence. Furthermore, members of the Mycobacterium tuberculosis complex have preferred geographic locations and possess unique phenotypic characteristics that is linked to their pathogenicity. Due to the recent advances in development next generation sequencing platforms, several studies have revealed epigenetic changes in genomic regions combined with "unique" gene regulatory mechanisms through non-coding RNAs that are responsible for strain-specific behaviour on in vitro and in vivo infection models. The current review provides up to date epigenetic patterns, gene regulation through non-coding RNAs, together with implications of these mechanisms in down-stream proteome and metabolome, which may be responsible for "unique" responses to infection by members of the Mycobacterium tuberculosis complex. Understanding lineage-specific molecular mechanisms during infection may provide novel drug targets and disease control measures towards World Health organization END-TB strategy.

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.

Determination of genetic diversity of multidrug-resistant Mycobacterium tuberculosis strains in Turkey using 15 locus MIRU-VNTR and spoligotyping methods

Tuberculosis (TB) remains the leading cause of deaths from infectious disease worldwide. Nowadays, the tendency of Mycobacterium tuberculosis complex (MTBC) to spread between continents due to uncontrolled migration movements shows that TB is a global health problem. The number of studies for the detection of MTBC strains' epidemiological features in areas with TB spread risk using molecular-based methods such as spoligotyping and Mycobacterial Interspersed Repetitive Unit (MIRU) Variable Number Tandem Repeats (VNTR) at the clonal level is insufficient. In this study, it was aimed to determine the phylogenetic relationships of MTBC strains at the species level by spoligotyping and 15 locus MIRU-VNTR (MIRU-VNTR₁₅) molecular methods of 96 multidrug-resistant (MDR) MTBC strains isolated from sputum samples of patients with a preliminary diagnosis of pulmonary TB or suspected contact history those sent to National Tuberculosis Reference Laboratory from the centers that are members of the Tuberculosis Laboratory Surveillance Network. The phylogenetic relationship between 96 MDR-TB strains was investigated with the combination of bead-based spoligotyping and MIRU-VNTR₁₅ methods on the MAGPIX® Milliplex Map device. In this study, it was determined that the T1 family is more common in our country and LAM7-TUR family is less common than the Beijing family unlike other studies. It was determined that the strains in the same cluster had different locus profiles, and there was no transmission from the same clone in the clonal typing we performed with spoligotyping and MIRU-VNTR₁₅.

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.

Multidrug-Resistant Tuberculosis—Diagnostic Procedures and Treatment of Two Beijing-like TB Cases

The Beijing/W genotype is one of the major molecular families of Mycobacterium tuberculosis complex (MTBC), responsible for approximately 50% of tuberculosis (TB) cases in Far East Asia and at least 25% of TB cases globally. Studies have revealed that the Beijing genotype family is associated with a more severe clinical course of TB, increased ability to spread compared to other genotypes, and an unpredictable response to treatment. Based on the profile of spacers 35–43 in the Direct Repeat (DR) locus of the MTBC genome determined by spoligotyping, classical (typical) and modern (Beijing-like) clones can be identified within the Beijing family. While the modern and ancient Beijing strains appear to be closely related at the genetic level, there are marked differences in their drug resistance, as well as their ability to spread and cause disease. This paper presents two cases of drug-resistant tuberculosis caused by rare mycobacteria from the Beijing family: the Beijing 265 and Beijing 541 subtypes. The genotypes of isolated strains were linked with the clinical course of TB, and an attempt was made to initially assess whether the Beijing subtype can determine treatment outcomes in patients.

Close Related Drug-Resistance Beijing Isolates of Mycobacterium tuberculosis Reveal a Different Transcriptomic Signature in a Murine Disease Progression Model

Mycobacterium tuberculosis (MTB) lineage 2/Beijing is associated with high virulence and drug resistance worldwide. In Colombia, the Beijing genotype has circulated since 1997, predominantly on the pacific coast, with the Beijing-Like SIT-190 being more prevalent. This genotype conforms to a drug-resistant cluster and shows a fatal outcome in patients. To better understand virulence determinants, we performed a transcriptomic analysis with a Beijing-Like SIT-190 isolate (BL-323), and Beijing-Classic SIT-1 isolate (BC-391) in progressive tuberculosis (TB) murine model. Bacterial RNA was extracted from mice lungs on days 3, 14, 28, and 60. On average, 0.6% of the total reads mapped against MTB genomes and of those, 90% against coding genes. The strains were independently associated as determined by hierarchical cluster and multidimensional scaling analysis. Gene ontology showed that in strain BL-323 enriched functions were related to host immune response and hypoxia, while proteolysis and protein folding were enriched in the BC-391 strain. Altogether, our results suggested a differential bacterial transcriptional program when evaluating these two closely related strains. The data presented here could potentially impact the control of this emerging, highly virulent, and drug-resistant genotype.

Snapshot of Mycobacterium tuberculosis Phylogenetics from an Indian State of Arunachal Pradesh Bordering China

Uncontrolled transmission of Mycobacterium tuberculosis (M. tuberculosis, MTB) drug resistant strains is a challenge to control efforts of the global tuberculosis program. Due to increasing multi-drug resistant (MDR) cases in Arunachal Pradesh, a northeastern state of India, the tracking and tracing of these resistant MTB strains is crucial for infection control and spread of drug resistance. This study aims to correlate the phenotypic DST, genomic DST (gDST) and phylogenetic analysis of MDR-MTB strains in the region. Of the total 200 samples 22 (11%) patients suspected of MDR-TB and 160 (80%) previously treated MDR-TB cases, 125 (62.5%) were identified as MTB. MGIT-960 SIRE DST detected 71/125 (56.8%) isolates as MDR/RR-MTB of which 22 (30.9%) were detected resistant to second-line drugs. Whole-genome sequencing of 65 isolates and their gDST found Ser315Thr mutation in katG (35/45; 77.8%) and Ser531Leu mutation in rpoB (21/41; 51.2%) associated with drug resistance. SNP barcoding categorized the dataset with Lineage2 (41; 63.1%) being predominant followed by Lineage3 (10; 15.4%), Lineage1 (8; 12.3%) and Lineage4 (6; 9.2%) respectively. Phylogenetic assignment by cgMLST gave insights of two Beijing sub-lineages viz; 2.2.1 (SNP difference < 19) and 2.2.1.2 (SNP difference < 9) associated with recent ongoing transmission in Arunachal Pradesh. This study provides insights in identifying two virulent Beijing sub-lineages (sub-lineage 2.2.1 and 2.2.1.2) with ongoing transmission of TB drug resistance in Arunachal Pradesh.

Snapshot of Mycobacterium tuberculosis Phylogenetics from an Indian State of Arunachal Pradesh Bordering China

Uncontrolled transmission of Mycobacterium tuberculosis (M. tuberculosis, MTB) drug resistant strains is a challenge to control efforts of the global tuberculosis program. Due to increasing multi-drug resistant (MDR) cases in Arunachal Pradesh, a northeastern state of India, the tracking and tracing of these resistant MTB strains is crucial for infection control and spread of drug resistance. This study aims to correlate the phenotypic DST, genomic DST (gDST) and phylogenetic analysis of MDR-MTB strains in the region. Of the total 200 samples 22 (11%) patients suspected of MDR-TB and 160 (80%) previously treated MDR-TB cases, 125 (62.5%) were identified as MTB. MGIT-960 SIRE DST detected 71/125 (56.8%) isolates as MDR/RR-MTB of which 22 (30.9%) were detected resistant to second-line drugs. Whole-genome sequencing of 65 isolates and their gDST found Ser315Thr mutation in katG (35/45; 77.8%) and Ser531Leu mutation in rpoB (21/41; 51.2%) associated with drug resistance. SNP barcoding categorized the dataset with Lineage2 (41; 63.1%) being predominant followed by Lineage3 (10; 15.4%), Lineage1 (8; 12.3%) and Lineage4 (6; 9.2%) respectively. Phylogenetic assignment by cgMLST gave insights of two Beijing sub-lineages viz; 2.2.1 (SNP difference < 19) and 2.2.1.2 (SNP difference < 9) associated with recent ongoing transmission in Arunachal Pradesh. This study provides insights in identifying two virulent Beijing sub-lineages (sub-lineage 2.2.1 and 2.2.1.2) with ongoing transmission of TB drug resistance in Arunachal Pradesh.

Snapshot of Mycobacterium tuberculosis Phylogenetics from an Indian State of Arunachal Pradesh Bordering China

Uncontrolled transmission of Mycobacterium tuberculosis (M. tuberculosis, MTB) drug resistant strains is a challenge to control efforts of the global tuberculosis program. Due to increasing multi-drug resistant (MDR) cases in Arunachal Pradesh, a northeastern state of India, the tracking and tracing of these resistant MTB strains is crucial for infection control and spread of drug resistance. This study aims to correlate the phenotypic DST, genomic DST (gDST) and phylogenetic analysis of MDR-MTB strains in the region. Of the total 200 samples 22 (11%) patients suspected of MDR-TB and 160 (80%) previously treated MDR-TB cases, 125 (62.5%) were identified as MTB. MGIT-960 SIRE DST detected 71/125 (56.8%) isolates as MDR/RR-MTB of which 22 (30.9%) were detected resistant to second-line drugs. Whole-genome sequencing of 65 isolates and their gDST found Ser315Thr mutation in katG (35/45; 77.8%) and Ser531Leu mutation in rpoB (21/41; 51.2%) associated with drug resistance. SNP barcoding categorized the dataset with Lineage2 (41; 63.1%) being predominant followed by Lineage3 (10; 15.4%), Lineage1 (8; 12.3%) and Lineage4 (6; 9.2%) respectively. Phylogenetic assignment by cgMLST gave insights of two Beijing sub-lineages viz; 2.2.1 (SNP difference < 19) and 2.2.1.2 (SNP difference < 9) associated with recent ongoing transmission in Arunachal Pradesh. This study provides insights in identifying two virulent Beijing sub-lineages (sub-lineage 2.2.1 and 2.2.1.2) with ongoing transmission of TB drug resistance in Arunachal Pradesh.

Beijing genotype of Mycobacterium tuberculosis is associated with extensively drug-resistant tuberculosis: A global analysis

We found that the frequency of Beijing genotype among XDR-TB strains was high. The data in this study would help guide the TB control program, and we however need further investigation to confirm the reliability of the present findings.

The Relevance of Genomic Epidemiology for Control of Tuberculosis in West Africa

Tuberculosis (TB), an airborne infectious disease caused by Mycobacterium tuberculosis complex (MTBC), remains a global health problem. West Africa has a unique epidemiology of TB that is characterized by medium- to high-prevalence. Moreover, the geographical restriction of M. africanum to the sub-region makes West Africa have an extra burden to deal with a two-in-one pathogen. The region is also burdened with low case detection, late reporting, poor treatment adherence leading to development of drug resistance and relapse. Sporadic studies conducted within the subregion report higher burden of drug resistant TB (DRTB) than previously thought. The need for more sensitive and robust tools for routine surveillance as well as to understand the mechanisms of DRTB and transmission dynamics for the design of effective control tools, cannot be overemphasized. The advancement in molecular biology tools including traditional fingerprinting and next generation sequencing (NGS) technologies offer reliable tools for genomic epidemiology. Genomic epidemiology provides in-depth insight of the nature of pathogens, circulating strains and their spread as well as prompt detection of the emergence of new strains. It also offers the opportunity to monitor treatment and evaluate interventions. Furthermore, genomic epidemiology can be used to understand potential emergence and spread of drug resistant strains and resistance mechanisms allowing the design of simple but rapid tools. In this review, we will describe the local epidemiology of MTBC, highlight past and current investigations toward understanding their biology and spread as well as discuss the relevance of genomic epidemiology studies to TB control in West Africa.

Molecular epidemiology and drug-resistance of tuberculosis in Luodian revealed by whole genome sequencing

In this study, we aimed to investigate the molecular epidemiology and drug-resistance profiles of tuberculosis (TB) in Luodian, an area with highest TB incidence and limited healthcare resources in Guizhou, China. The passive case finding strategy was used to identify suspected pulmonary TB with symptoms, and individuals with positive Mycobacterium tuberculosis (MTB) culture were enrolled from May 22, 2018 to April 21, 2019. All the 107 cases except three came from nine towns, including 55.1% from Longping and Bianyang. The phylogeny tree showed that 53.3% of strains were Lineage 2 (Beijing genotype), while 46.7% were Lineage 4 (Euro-American genotype). Among Lineage 2 strains, 66.7% were of "modern" Beijing type. Seven clusters with genomic distance within 12 SNPs were identified. The clusters included 14 strains, accounting for a clustering rate of 13.1%. The distance separating the clustered cases was between 2.1 and 71.0 km (Km), with an average paired distance of 21.8 Km (interquartile range, 2.8-38.0 Km). Based on the gene mutations associated with drug-resistance, we predicted that 4.8% of strains were resistant to isoniazid, 3.7% to rifampicin, and 3.7% to streptomycin; only one strain (0.9%) had multidrug resistance (MDR). This study found low drug-resistance rates in Luodian, and the sub-lineage of the "modern" Beijing branch has recent expansion in Luodian. This work may also serve as a genomic baseline to assess the evolution and spread of MTB in Guizhou.

Magnitude of Phenotypic and MTBDRplus Line Probe Assay First-Line Anti-Tuberculosis Drug Resistance Among Tuberculosis Patients; Northwest Ethiopia

Background

Mycobacterium tuberculosis (Mtb) drug resistance is a key challenge in ending TB.

Objective

The study aimed to determine anti-TB drug resistance and compare the discordance between phenotypic and genotypic drug-susceptibility testing (DST).

Methods

Prospective enrollment and sputum collection from patients suspected of active pulmonary TB from May 2018 to December 2019 at the University of Gondar Hospital. Phenotypic DST study for streptomycin, isoniazid, rifampin, and ethambutol was done by MGIT 360 SIRE Kit. Genotypic resistance for isoniazid and rifampin was performed by MTBDRplus v2 line probe assay (LPA) and compared to phenotypic drug resistance.

Results

A total of 376 patients, median age 32 years, and 53.7% male were enrolled. Mtb was isolated from 126 patients. 106/126 (84%) patients were newly diagnosed with TB and 20 patients with prior TB treatment. Seventy (66.0%) were susceptible to all anti‐TB drugs tested. Twenty-five (19.8%) of the isolates were resistant to isoniazid, 12 (9.5%) to rifampicin and six (5%) were multidrug resistant. Among previously treated TB patients, 4 (20.0%) and 5 (25.0%) were mono-resistant and poly-resistant, respectively. The sensitivity and specificity of LPA resistance for isoniazid were 94.4% and 100%, and for rifampin was 75.0% and 100%, respectively.

Conclusion

The frequency of mono- and poly-drug resistance among both newly diagnosed and previously treated TB patients was high to the rest of the nation. MTBDRplus showed excellent concordance for isoniazid and rifampin. We concluded that DST should be performed for all patients to improve management and decrease spread of drug-resistant Mtb strains in the community.