Population structure and genetic diversity of Mycobacterium tuberculosis in Ecuador

Molecular epidemiology of Mycobacterium tuberculosis in Ecuador: Recent advances and future challenges

Tuberculosis (TB) is one of the three leading causes of death from a single infectious agent, Mycobacterium tuberculosis (MTB), together with COVID-19 and HIV/AIDS. This disease places a heavy burden on countries with low socio-economic development and aggravates existing inequalities. For the year 2021, estimations for Ecuador were 8500 TB cases, of which 370 were associated to multiple drug resistance (TB-MDR), and 1160 deaths. In the same year, Ecuador notified 5973 total cases, 401 of them were TB-MDR, pointing out an under diagnosis problem. The few molecular epidemiology studies available conclude that L4 is the most prevalent MTB lineage in Ecuador (with LAM as the main L4 sublineage), but L2-Beijing family is also present at low prevalence. Nevertheless, with less than 1 % MTB isolates genetically characterized by either MIRU-VNTR, spolygotyping or WGS to date, molecular epidemiology research must me improved to assist the TB surveillance and control program in Ecuador.

Prevalence, drug resistance, and genotypic diversity of the RDRio subfamily of Mycobacterium tuberculosis in Ecuador: a retrospective analysis for years 2012-2016

Introduction

A major sublineage within the Mycobacterium tuberculosis (MTB) LAM family characterized by a new in-frame fusion gene Rv3346c/55c was discovered in Rio de Janeiro (Brazil) in 2007, called RDRio, associated to drug resistance. The few studies about prevalence of MTB RDRio strains in Latin America reported values ranging from 3% in Chile to 69.8% in Venezuela, although no information is available for countries like Ecuador.

Methods

A total of 814 MTB isolates from years 2012 to 2016 were screened by multiplex PCR for RDRio identification, followed by 24-loci MIRU-VNTR and spoligotyping.

Results

A total number of 17 MTB RDRio strains were identified, representing an overall prevalence of 2.09% among MTB strains in Ecuador. While 10.9% of the MTB isolates included in the study were multidrug resistance (MDR), 29.4% (5/17) of the RDRio strains were MDR.

Discussion

This is the first report of the prevalence of MTB RDRio in Ecuador, where a strong association with MDR was found, but also a very low prevalence compared to other countries in Latin America. It is important to improve molecular epidemiology tools as a part of MTB surveillance programs in Latin America to track the transmission of potentially dangerous MTB stains associated to MDR TB like MTB RDRio.

Rapid and accurate identification and differentiation of Mycobacterium tuberculosis and non-tuberculous mycobacteria using PCR kits available in a high-burden setting

Infections caused by mycobacteria, including Mycobacterium tuberculosis complex (MTBC) and non-tuberculous mycobacteria (NTM), are a major public health issue worldwide. An accurate diagnosis of mycobacterial species is a challenge for surveillance and treatment, particularly in high-burden settings usually associated with low- and middle-income countries. In this study, we analyzed the clinical performance of two commercial PCR kits designed for the identification and differentiation of MTBC and NTM, available in a high-burden setting such as Ecuador. A total of 109 mycobacteria isolates were included in the study, 59 of which were previously characterized as M. tuberculosis and the other 59 as NTM. Both kits displayed great clinical performance for the identification of M. tuberculosis, with 100% sensitivity. On the other hand, for NTM, one of the kits displayed a good clinical performance with a sensitivity of 94.9% (CI 95%: 89-100%), while the second kit had a reduced sensitivity of 77.1% (CI 95%: 65-89%). In conclusion, one of the kits is a fast and reliable tool for the identification and discrimination of MTBC and NTM from clinical isolates.

Population structure of Mycobacterium tuberculosis in El Oro: A first insight into Ecuador-Peru tuberculosis transmission

BACKGROUND

Tuberculosis (TB) is a major public health concern in Ecuador and Peru, both settings of high burden of drug resistance TB. Molecular epidemiology tools are important to understand the transmission dynamics of Mycobacterium tuberculosis Complex (MTBC) and to track active transmission clusters of regional importance. This study is the first to address the transmission of TB between Peru and Ecuador through the population structure of MTBC lineages circulating in the Ecuadorian border province of "El Oro".

METHODS

A total number of 56 MTBC strains from this province for years 2012-2015 were included in the study and analyzed by 24-loci MIRU-VNTR and spoligotyping.

RESULTS

Genotyping revealed a high degree of diversity for MTBC in "El Oro", without active transmission clusters. MTBC L4 was predominant, with less than 2% of strains belonging to MTBC L2-Beijing.

CONCLUSIONS

These results may suggest that TB dynamics in this rural and semi-urban area would not be linked to highly transmitted strains like MTBC L2-Beijing from Peru, but related to TB relapse; although further studies with larger MTBC cultures collection from recent years are needed. Nevertheless, we recommend to reinforce TB surveillance programs in remote rural settings and border regions in Ecuador.

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.

A precision overview of genomic resistance screening in Ecuadorian isolates of Mycobacterium tuberculosis using web-based bioinformatics tools

INTRODUCTION

Tuberculosis (TB) is among the deadliest diseases worldwide, and its impact is mainly due to the continuous emergence of resistant isolates during treatment due to the laborious process of resistance diagnosis, nonadherence to treatment and circulation of previously resistant isolates of Mycobacterium tuberculosis. In this study, we evaluated the performance and functionalities of web-based tools, including Mykrobe, TB-profiler, PhyResSE, KvarQ, and SAM-TB, for detecting resistance in 88 Ecuadorian isolates of Mycobacterium tuberculosis drug susceptibility tested previously. Statistical analysis was used to determine the correlation between genomic and phenotypic analysis. Our results showed that with the exception of KvarQ, all tools had the highest correlation with the conventional drug susceptibility test (DST) for global resistance detection (98% agreement and 0.941 Cohen's kappa), while SAM-TB, PhyResSE, TB-profiler and Mykrobe had better correlations with DST for first-line drug analysis individually. We also identified that in our study, only 50% of mutations characterized by the web-based tools in the rpoB, katG, embB, pncA, gyrA and rrs regions were canonical and included in the World Health Organization (WHO) catalogue. Our findings suggest that SAM-TB, PhyResSE, TB-profiler and Mykrobe were efficient in determining canonical resistance-related mutations, but more analysis is needed to improve second-line detection. Improving surveillance programs using whole-genome sequencing tools for first-line drugs, MDR-TB and XDR-TB is essential to understand the molecular epidemiology of TB in Ecuador.

IMPORTANCE

Tuberculosis, an infectious disease caused by Mycobacterium tuberculosis, most commonly affects the lungs and is often spread through the air when infected people cough, sneeze, or spit. However, despite the existence of effective drug treatment, patient adherence, long duration of treatment, and late diagnosis have reduced the effectiveness of therapy and increased drug resistance. The increase in resistant cases, added to the impact of the COVID-19 pandemic, has highlighted the importance of implementing efficient and timely diagnostic methodologies worldwide. The significance of our research is in evaluating and identifying a more efficient and user-friendly web-based tool to characterize resistance in Mycobacterium tuberculosis by whole-genome sequencing, which will allow more routine application to improve TB strain surveillance programs locally.

Large genomic deletions delineate Mycobacterium tuberculosis L4 sublineages in South American countries

Mycobacterium tuberculosis (Mtb) is still one of the primary pathogens of humans causing tuberculosis (TB) disease. Mtb embraces nine well-defined phylogenetic lineages with biological and geographical disparities. The lineage L4 is the most globally widespread of all lineages and was introduced to America with European colonization. Taking advantage of many genome projects available in public repositories, we undertake an evolutionary and comparative genomic analysis of 522 L4 Latin American Mtb genomes. Initially, we performed careful quality control of public read datasets and applied several thresholds to filter out low-quality data. Using a genome de novo assembly strategy and phylogenomic methods, we spotted novel south American clades that have not been revealed yet. Additionally, we describe genomic deletion profiles of these strains from an evolutionary perspective and report Mycobacterium tuberculosis L4 sublineages signature-like gene deletions, some of the novel. One is a specific deletion of 6.5 kbp that is only present in sublineage 4.1.2.1. This deletion affects a complex group of 10 genes with putative products annotated, among others, as a lipoprotein, transmembrane protein, and toxin/antitoxin system proteins. The second novel deletion spans for 4.9 kbp and specific of a particular clade of the 4.8 sublineage and affects 7 genes. The last novel deletion affects 4 genes, extends for 4.8 kbp., and is specific to some strains within the 4.1.2.1 sublineage that are present in Colombia, Peru and Brasil.

Genomic evidence for adaptation to tuberculosis in the Andes before European contact

Most studies focusing on human high-altitude adaptation in the Andean highlands have thus far been focused on Peruvian populations. We present high-coverage whole genomes from Indigenous people living in the Ecuadorian highlands and perform multi-method scans to detect positive natural selection. We identified regions of the genome that show signals of strong selection to both cardiovascular and hypoxia pathways, which are distinct from those uncovered in Peruvian populations. However, the strongest signals of selection were related to regions of the genome that are involved in immune function related to tuberculosis. Given our estimated timing of this selection event, the Indigenous people of Ecuador may have adapted to Mycobacterium tuberculosis thousands of years before the arrival of Europeans. Furthermore, we detect a population collapse that coincides with the arrival of Europeans, which is more severe than other regions of the Andes, suggesting differing effects of contact across high-altitude populations.

Genomic Diversity of the Rarely Observed Genotype of the Mycobacterium tuberculosis Central Asian (CAS) Lineage 3 from North Brazil

Mycobacterium tuberculosis (Mtb) Central Asian Strain (CAS) Lineage 3 (L3) genotype is predominantly found in East-Africa, Central-Asia, Western-Asia, and South-Asia; however, a new spoligotyping CAS/SIT2545 was found in northern regions of Brazil. We aimed to characterize and describe the genetic diversity and perform a phylogenetic assessment of this novel genotype. We performed 24-MIRU-VNTR loci and Whole-genome sequencing (WGS) of six Brazilian isolates previously spoligotyped. The libraries were prepared using a Nextera-XT kit and sequenced in a NextSeq 550 Illumina instrument. We performed lineage assignment and genomic characterization. From publicly available genomes of Mtb L3 and other lineages, we created a robust dataset to run the MTBSeq pipeline and perform a phylogenetic analysis. MIRU-VNTR and WGS confirmed CAS/SIT2545 belongs to L3. Out of 1691 genomes, 1350 (79.83%) passed in quality control (genomic coverage > 95%). Strain 431 differed in 52 single nucleotide variants (SNV), confirming it does not belong to the same transmission chain. The eight genomes from a global dataset clustered closer to Brazilian strains differed in >52 SNVs. We hypothesized L3 and L1 were introduced in Brazilian Northern in the same historical event; however, there is a need for additional studies exploring the genetic diversity of Mtb Brazilian Northern.

Evolutionary history and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage, Tunisia

BACKGROUND

To infer the origin and spread of the Mycobacterium tuberculosis Latin American and Mediterranean (L4.3/LAM) sublineage in a Mediterranean country, Tunisia, where it predominates.

METHODS

We combined Bayesian (STRUCTURE) and maximum likelihood (MIGRAINE) estimation approaches based on a global 24-loci mycobacterial interspersed repetitive units-variable numbers of tandem repeats (MIRU-VNTR24) genotyping dataset consisting of 1573 L4.3/LAM clinical strains from four continents, including 252 isolates originating from Tunisia.

RESULTS

Phylogenetic analyses coupled with Bayesian estimations suggested that the most predominant L4.3/LAM subpopulation in Tunisia (65.07%), which is dominated by a single clonal complex, TUN4.3_CC1 (94.51%), has evolved from an ancestral pool that is restricted to Europe and Africa, contrasting with the remaining L4.3/LAM subpopulations whose ancestry was traced all over the word. Maximum likelihood analysis revealed that TUN4.3_CC1 has been undergoing a demographic expansion since 131 years ago (CI95%: 90.7-205), thus explaining its preponderance relative to the second most predominant CC, TUN4.3_CC2, whose population was found under contraction.

CONCLUSIONS

The preponderance of L4.3/LAM in Tunisia stems from a 130-year expansion process of a locally evolved clone.

Comparative genomics of drug-resistant strains of Mycobacterium tuberculosis in Ecuador

BACKGROUND

Tuberculosis is a serious infectious disease affecting millions of people. In spite of efforts to reduce the disease, increasing antibiotic resistance has contributed to persist in the top 10 causes of death worldwide. In fact, the increased cases of multi (MDR) and extreme drug resistance (XDR) worldwide remains the main challenge for tuberculosis control. Whole genome sequencing is a powerful tool for predicting drug resistance-related variants, studying lineages, tracking transmission, and defining outbreaks. This study presents the identification and characterization of resistant clinical isolates of Mycobacterium tuberculosis including a phylogenetic and molecular resistance profile study by sequencing the complete genome of 24 strains from different provinces of Ecuador.

RESULTS

Genomic sequencing was used to identify the variants causing resistance. A total of 15/21 isolates were identified as MDR, 4/21 as pre-XDR and 2/21 as XDR, with three isolates discarded due to low quality; the main sub-lineage was LAM (61.9%) and Haarlem (19%) but clades X, T and S were identified. Of the six pre-XDR and XDR strains, it is noteworthy that five come from females; four come from the LAM sub-lineage and two correspond to the X-class sub-lineage. A core genome of 3,750 genes, distributed in 295 subsystems, was determined. Among these, 64 proteins related to virulence and implicated in the pathogenicity of M. tuberculosis and 66 possible pharmacological targets stand out. Most variants result in nonsynonymous amino acid changes and the most frequent genotypes were identified as conferring resistance to rifampicin, isoniazid, ethambutol, para-aminosalicylic acid and streptomycin. However, an increase in the resistance to fluoroquinolones was detected.

CONCLUSION

This work shows for the first time the variability of circulating resistant strains between men and women in Ecuador, highlighting the usefulness of genomic sequencing for the identification of emerging resistance. In this regard, we found an increase in fluoroquinolone resistance. Further sampling effort is needed to determine the total variability and associations with the metadata obtained to generate better health policies.

Fingerprinting of Mycobacterium tuberculosis isolates by MIRU-VNTR genotyping and detection of isoniazid resistance by real-time PCR

Introduction. Tuberculosis (TB) is a great public health problem in developing countries such as Egypt. Genotyping of Mycobacterium tuberculosis isolates has a prominent role in the field of TB prevention.

Aim. This study aimed to evaluate real-time PCR using Minor Groove Binder (MGB) probes and to identify circulating lineages/sub-lineages of M. tuberculosis and their transmission patterns.

Hypothesis. We hypothesize that MIRU-VNTR technique is efficient in identifying circulating M. tuberculosis lineages in Egypt.

Methodology. Fifty sputum specimens positive for acid-fast bacilli were included. Isoniazid (INH) resistance was detected using the 1 % proportion method. Real-time PCR using MGB-probes was used for simultaneous detection of TB infection and INH resistance. Partial sequencing of the katG gene was used to confirm INH resistance results. A standard 15 Mycobacterial Interspersed Repetitive Unit Variable Number Tandem Repeat (15-MIRU-VNTR) approach was used for genotyping through the MIRU-VNTRplus online platform.

Results. Only seven specimens showed phenotypic resistance to INH. M. tuberculosis was detected in all samples, while a mutation in the katG gene codon 315 was detected only in five samples, which were also phenotypically INH-resistant. Sequencing of the katG gene showed codon 315 mutation genotypically and phenotypically in the five INH-resistant isolates. Molecular genotyping of M. tuberculosis isolates revealed that the majority of isolates (26/50, 52 %) belonged to the S family of lineage_4. A low clustering rate (2 %) was observed among our isolates. According to the Hunter-Gaston Discriminatory Index (HGDI), 11 MIRU-VNTR loci were highly or moderately discriminative, while four loci were less polymorphic.

Conclusion. MIRU-VNTR genotyping revealed a low clustering rate with a low recent transmission rate of M. tuberculosis strains in Alexandria, Egypt.

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

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

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

Spatial clustering and genetic diversity of Mycobacterium tuberculosis isolate among pulmonary tuberculosis suspected patients, Arsi Zone, Ethiopia

BACKGROUND

Tuberculosis remains a serious public health concern globally. The enormous social, economic, and health impacts of the diseases are attributed to the lack of updated data on the prevalence, geospatial distribution, population structures, and genotypic variants of the circulating M. tuberculosis.

METHODS

Structured questionnaire, mycobacterial culture, and standard 24-Mycobacterial Interspersed Repeated Units-Variable Number Tandem Repeats (MIRU-VNTR) were employed to collect sociodemographic characters, residence linked information, and genotype the isolates. The retrospective discrete Bernoulli model was used to identify the hot spot districts of sputum smear positivity, and Web-based Miru-VNTRPlus were used for the identification of lineages and sublineages.

RESULTS

Out of 832 presumptive pulmonary tuberculosis (PTB) suspects, 119 (14.3%) were smear-positive. In the multivariate binary logistic model, PTB suspected patients in the age groups of 7-25 and 25-34 and those from rural residents were 4.53 (AOR = 4.53; 95% CI 2.25-9.13), 3.00 (AOR = 3.00; 95% CI 1.41-6.35) and 1.65 (AOR = 1.65; 95% CI 1.01-2.70) times at higher risk of turning smear-positive. Eleven (47.8%) districts of Arsi Zone were shown to have a high rate of clustering (RR = 2.27; 95% CI 1.62-3.2) of smear-positive PTB. Of 72 isolates queried for the lineage assignment, 59 (81.9%) were classified into the previously known lineages and 13 (18.1%) were not assigned to any known lineages. Overall, 42 (58.3%) belong to M. tuberculosis lineage 4 (Euro-American), 16 (22.2%) M. tuberculosis lineage 3 (Delhi/CAS), and 1 (1.4%) M. tuberculosis Lineage 1 (Indo-Oceanic/ East Africa Indian). Further classification to the sublineage indicates that the predominant lineage was Delhi/CAS comprising 16 (22.2%) isolates followed by 15 (20.8%) isolates belonging to Haarlem. The remaining isolates were distributed as 13 (18.1%) TUR, 6 (8.3%) LAM, 4 (5.5%) URAL, 4 (4.5%) NEW-1 and 1 (1.4%) EAI.

CONCLUSION

Our study showed higher smear-positive results among PTB suspected patients and remarkable spatial variation across districts of Arsi Zone in smear-positive PTB. This information together with the genotypic features could be used as input for the efforts of designing control strategies.

Whole genome analysis of extensively drug resistant Mycobacterium tuberculosis strains in Peru

Peru has the highest burden of multidrug-resistant tuberculosis in the Americas region. Since 1999, the annual number of extensively drug-resistant tuberculosis (XDR-TB) Peruvian cases has been increasing, becoming a public health challenge. The objective of this study was to perform genomic characterization of Mycobacterium tuberculosis strains obtained from Peruvian patients with XDR-TB diagnosed from 2011 to 2015 in Peru. Whole genome sequencing (WGS) was performed on 68 XDR-TB strains from different regions of Peru. 58 (85.3%) strains came from the most populated districts of Lima and Callao. Concerning the lineages, 62 (91.2%) strains belonged to the Euro-American Lineage, while the remaining 6 (8.8%) strains belonged to the East-Asian Lineage. Most strains (90%) had high-confidence resistance mutations according to pre-established WHO-confident grading system. Discordant results between microbiological and molecular methodologies were caused by mutations outside the hotspot regions analysed by commercial molecular assays (rpoB I491F and inhA S94A). Cluster analysis using a cut-off ≤ 10 SNPs revealed that only 23 (34%) strains evidenced recent transmission links. This study highlights the relevance and utility of WGS as a high-resolution approach to predict drug resistance, analyse transmission of strains between groups, and determine evolutionary patterns of circulating XDR-TB strains in the country.