Genetic diversity, population structure and drug resistance of Mycobacterium tuberculosis in Peru

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.

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.

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.

Genomic signatures of pre-resistance in Mycobacterium tuberculosis

Recent advances in bacterial whole-genome sequencing have resulted in a comprehensive catalog of antibiotic resistance genomic signatures in Mycobacterium tuberculosis. With a view to pre-empt the emergence of resistance, we hypothesized that pre-existing polymorphisms in susceptible genotypes (pre-resistance mutations) could increase the risk of becoming resistant in the future. We sequenced whole genomes from 3135 isolates sampled over a 17-year period. After reconstructing ancestral genomes on time-calibrated phylogenetic trees, we developed and applied a genome-wide survival analysis to determine the hazard of resistance acquisition. We demonstrate that M. tuberculosis lineage 2 has a higher risk of acquiring resistance than lineage 4, and estimate a higher hazard of rifampicin resistance evolution following isoniazid mono-resistance. Furthermore, we describe loci and genomic polymorphisms associated with a higher risk of resistance acquisition. Identifying markers of future antibiotic resistance could enable targeted therapy to prevent resistance emergence in M. tuberculosis and other pathogens.

Dysglycemia is associated with Mycobacterium tuberculosis lineages in tuberculosis patients of North Lima-Peru

This study was performed to investigate the role of dysglycemia on the genetic diversity of Mycobacterium tuberculosis (MTB) among pulmonary tuberculosis (TB) patients to build scientific evidence about the possible mechanisms of TB transmission. MTB isolates obtained of patients affected by pulmonary tuberculosis from health care facilities of North Lima-Peru, were analyzed using whole genome sequencing and 24-locus mycobacterial interspersed repetitive-unit -variable-number tandem repeats (MIRU-VNTR). Subsequently, clinical and epidemiological characteristics were associated with clustering, lineages and comorbid conditions. The analysis carried out 112 pulmonary TB patients from various health centers in North Lima, 17 (15%) had diabetes mellitus (DM) and 33 (29%) had pre-diabetes (PDM). Latin American-Mediterranean, Haarlem and Beijing were the most frequent MTB lineages found in those patients. Previous TB (adjusted odds ratio [aOR] = 3.65; 95%CI: 1.32-17.81), age (aOR = 1.12; 95%CI: 1.03-1.45) and Beijing lineage (aOR = 3.53; 95%CI: 1.08-13.2) were associated with TB-DM comorbidity. Alcoholism (aOR = 2.92; 95%CI: 1.10-8.28), age (aOR = 1.05; 95%CI: 1.03-1.12) and Haarlem lineage (aOR = 2.54; 95%CI: 1.04-6.51) were associated with TB-PDM comorbidity. Beijing and Haarlem lineages were independently associated with TB-DM and TB-PDM comorbidities, respectively. Although these findings may be surprising, we must be cautious to suggest that dysglycemia could be associated with a highly clustering and predisposition of MTB lineages related to a serious impact on the severity of TB disease, which requires further research.

First approach to the population structure of Mycobacterium tuberculosis complex in the indigenous population in Puerto Nariño-Amazonas, Colombia

Introduction

Tuberculosis affects vulnerable groups to a greater degree, indigenous population among them.

Objective

To determine molecular epidemiology of clinical isolates of Mycobacterium tuberculosis circulating in an indigenous population through Spoligotyping and 24-loci MIRU-VNTR.

Methodology

A descriptive cross-sectional study was conducted in 23 indigenous communities of Puerto Nariño-Amazonas, Colombia. Recovered clinical isolates were genotyped. For genotyping analyzes global SITVIT2 database and the MIRU-VNTRplus web portal were used.

Results

74 clinical isolates were recovered. Genotyping of clinical isolates by spoligotyping determined 5 different genotypes, all of them belonged to Euro-American lineage. By MIRU-VNTR typing, a total of 14 different genotypes were recorded. Furthermore, polyclonal infection was found in two patients from the same community. The combination of the two methodologies determined the presence of 19 genotypes, 8 formed clusters with 63 clinical isolates in total. Based on epidemiological information, it was possible to establish a potential chain of active transmission in 10/63 (15.9%) patients.

Conclusions

High genomic homogeneity was determined in the indigenous population suggesting possible chains of active transmission. The results obtained showed that specific genotypes circulating among the indigenous population of Colombia are significantly different from those found in the general population.

Population structure and genetic diversity of Mycobacterium tuberculosis in Ecuador

Tuberculosis (TB) is a significant public health problem in Ecuador with an incidence of 43 per 100,000 inhabitants and an estimated multidrug-resistant-TB prevalence in all TB cases of 9%. Genotyping of Mycobacterium tuberculosis (MTBC) is important to understand regional transmission dynamics. This study aims to describe the main MTBC lineages and sublineages circulating in the country. A representative sample of 373 MTBC strains from 22 provinces of Ecuador, with data comprising geographic origin and drug susceptibility, were genotyped using 24 loci-MIRU-VNTR. For strains with an ambiguous sublineage designation, the lineage was confirmed by Regions of Difference analysis or by Whole Genome Sequencing. We show that lineage 4 is predominant in Ecuador (98.3% of the strains). Only 4 strains belong to lineages 2-sublineage Beijing and two strains to lineage 3-sublineage Delhi. Lineage 4 strains included sublineages LAM (45.7%), Haarlem (31.8%), S (13.1%), X (4.6%), Ghana (0.6%) and NEW (0.3%). The LAM sublineage showed the strongest association with antibiotic resistance. The X and S sublineages were found predominantly in the Coastal and the Andean regions respectively and the reason for the high prevalence of these strains in Ecuador should be addressed in future studies. Our database constitutes a tool for MIRU-VNTR pattern comparison of M. tuberculosis isolates for national and international epidemiologic studies and phylogenetic purposes.

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

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

Country-wide rapid screening for the Mycobacterium tuberculosis Beijing sublineage in Ecuador using a single-nucleotide polymorphism-polymerase chain reaction method

Background

Strains of the Beijing sublineage of Mycobacterium tuberculosis have caused large outbreaks of tuberculosis, often involving multidrug resistance strains and this genetically highly conserved family of strains predominates in some geographic areas. For most of the countries of Latin America, no country-wide studies about the prevalence of the Beijing lineage are available.

Methods

In this study, we determine the prevalence of the Beijing sublineage in Ecuador, using a large nation-wide sample of 991 isolates from the years 2014-2016 and with the strains, in case-related-proportional representation, emerging from most of the provinces of the country. The isolates were genotyped with asinglenucleotidespecific polymorphism (SNP) polymerase chain reaction for the Beijing sublineage. SNPpositive strains were confirmed as belonging to this lineage with 24 mycobacterial interspersed repetitive unitvariable number of tandem repeat and DNA sequencing.

Results

We identified only four Beijing isolates in this collection of 991 strains and calculated a prevalence rate of 0.43%.

Conclusions

Our study shows a limited dissemination of the Beijing strains in the Ecuadorian population. This in contrast with the neighbor countries of Peru and Colombia were locally a prevalence of up to 16% has been reported.

Exploring the "Latin American Mediterranean" family and the RDRio lineage in Mycobacterium tuberculosis isolates from Paraguay, Argentina and Venezuela

Background

The Latin American & Mediterranean (LAM) spoligotype family is one of the most successful genotype of Mycobacterium tuberculosis worldwide and particularly prevalent in South-America. Within this family, a sublineage named Region of Difference Rio (RD^Rio) was reported initially in Brazil and is characterized by a genomic deletion of about 26.3 kb. This lineage seems to show a specific adaptation to the Euro-Latin American population. In this context, we sought to evaluate the LAM family and the presence of the RD^Rio genotype in samples from three Latin American countries including Paraguay, Venezuela and Argentina. To detect LAM strains reliably we applied a typing scheme using spoligotyping, 12 loci MIRU-VNTR, the Ag85C¹⁰³ SNP and the regions of difference RD^Rio and RD174. IS6110-RFLP results were also used when available.

Results

Genotyping of 413 M. tuberculosis isolates from three Latin-American countries detected LAM (46%) and the ill-defined T clade (16%) as the most frequent families. The highest clustering rate was detected in the sample population from the city of Caracas in Venezuela. We observed considerable differences in the presence of the RD^Rio lineage, with high frequency in Caracas-Venezuela (55%) and low frequency in Buenos Aires-Argentina (11%) and Paraguay (10%). The molecular markers (RD174, Ag85C¹⁰³, MIRU02-MIRU40 signature) of the RD^Rio lineage were essentially confirmed. For the LAM family, the most polymorphic loci were MIRU40, MIRU31, MIRU10, MIRU26, MIRU16 and the least polymorphic MIRU24, MIRU20, MIRU04, MIRU23.

Conclusions

Our results suggest a differential adaptation of LAM-sublineages in neighboring populations and that RD^Rio strains spread regionally with different rates of distribution. The Ag85C SNP and RDs (RD174, RD^Rio) tested in this study can in fact facilitate molecular epidemiological studies of LAM strains in endemic settings and low-income countries.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1479-6) contains supplementary material, which is available to authorized users.

Tuberculosis Genetic Epidemiology: A Latin American Perspective

There are an estimated 10 million new cases of tuberculosis worldwide annually, with 282,000 new or relapsed cases each year reported from the Americas. With improvements in genome sequencing technology, it is now possible to study the genetic diversity of tuberculosis with much greater resolution. Although tuberculosis bacteria do not engage in horizontal gene transfer, the genome is far more variable than previously thought. The study of genome-wide variation in tuberculosis has improved our understanding of the evolutionary origins of tuberculosis, the arrival of tuberculosis in Latin America, the genetic determinants of drug resistance, and lineage-specific associations with important clinical phenotypes. This article reviews what is known about the arrival of tuberculosis in Latin America, the genetic diversity of tuberculosis in Latin America, and the genotypic determinants of clinical phenotypes.

Genetic diversity of drug and multidrug-resistant Mycobacterium tuberculosis circulating in Veracruz, Mexico

Background

Mexico is one of the most important contributors of drug and multidrug-resistant tuberculosis in Latin America; however, knowledge of the genetic diversity of drug-resistant tuberculosis isolates is limited.

Methods

In this study, the genetic structure of 112 Mycobacterium tuberculosis strains from the southeastern Mexico was determined by spoligotyping and 24-loci MIRU-VNTRs.

Findings

The results show eight major lineages, the most of which was T1 (24%), followed by LAM (16%) and H (15%). A total of 29 (25%) isolates were identified as orphan. The most abundant SITs were SIT53/T1 and SIT42/LAM9 with 10 isolates each and SIT50/H3 with eight isolates. Fifty-two spoligotype patterns, twenty-seven clusters and ten clonal complexes were observed, demonstrating an important genetic diversity of drug and multidrug-resistant tuberculosis isolates in circulation and transmission level of these aggravated forms of tuberculosis. Being defined as orphan or as part of an orphan cluster, was a risk factor for multidrug resistant-tuberculosis (OR 2.5, IC 1.05–5.86 and OR 3.3, IC 1–11.03, respectively). Multiple correspondence analyses showed association of some clusters and SITs with specific geographical locations.

Conclusions

Our study provides one of the most detailed description of the genetic structure of drug and multidrug-resistant tuberculosis strains in southeast Mexico, establishing for the first time a baseline of the genotypes observed in resistant isolates circulating, however further studies are required to better elucidate the genetic structure of tuberculosis in region and the factors that could be participating in their dispersion.

Spoligotype defined lineages of Mycobacterium tuberculosis and drug resistance: Merely a casual correlation?

Drug-resistant tuberculosis (TB) is a major challenge to TB control strategy worldwide. Analysis of genetic polymorphism among drug resistant Mycobacterium tuberculosis (MTB) strains may help provide some insight into the transmission dynamics of these strains. Spoligotyping is a widely used technique to identify genetic polymorphism, based on 43 known spacers interspersed between direct repeat regions. Considerable work has been done in various parts of the world using this technique to identify and analyse the polymorphic nature of MTB. Many studies have been carried out to determine the association of drug resistance with spoligotype defined lineages, and much data has been produced over the years. New information continues to be generated. This review aims to put together the findings of relevant studies in an attempt to understand the correlation of drug resistance with spoligotype defined lineages of MTB. This would help provide a perspective of the available data that can be used as a starting point to understand the molecular epidemiology of drug resistant TB.

[Characterization of clinical isolates of Mycobacterium tuberculosis from HIV positive individuals in Colombia, 2012]

INTRODUCTION

One third of the increase in tuberculosis cases is attributed to the spread of HIV. In 2012, 1,397 HIV-associated tuberculosis cases were reported in Colombia, i.e., 11.8% of the total cases. Molecular epidemiology tools help to understand the transmission of tuberculosis.

OBJECTIVE

To characterize clinical isolates of Mycobacterium tuberculosis derived from HIV-infected individuals, received at the Laboratorio Nacional de Referencia in the Instituto Nacional de Salud.

MATERIALS AND METHODS

This was a descriptive observational study. We analyzed 63 isolates of M. tuberculosis from HIV-infected individuals. Identification, drug susceptibility and genotyping assays were performed.

RESULTS

Of the new cases evaluated, three (5.0%) were resistant to isoniazid combined with streptomycin; two (3.3%) to rifampicin, and one (1.6%) to isoniazid. Previously treated cases were sensitive. No multidrug resistance was evident. Among the predominant genotypes, 20 isolates were (31.7%) LAM9, eight (12.7%), H1, and seven (11.1%), T1. Nineteen isolates corresponded to orphan patterns. One single grouping was observed among tested isolates. We found no statistically significantdifference between the proportions of the antituberculous drug resistance and genotypes.

CONCLUSION

We found resistant isolates to the most powerful drugs, rifampicin and isoniazid, among new cases, showing the transmission of resistant strains. Genetic families of M. tuberculosis LAM9, T1 and H1 correspond to those described in the general population. We detected no active transmission among studied isolates. More comprehensive studies are needed to assess the real situation of HIV associated tuberculosis in the country regarding sensitivity and transmission.

Analysis of Mycobacterium tuberculosis Genotypic Lineage Distribution in Chile and Neighboring Countries

Tuberculosis (TB), caused by the pathogen Mycobacterium tuberculosis (MTB), remains a disease of high importance to global public health. Studies into the population structure of MTB have become vital to monitoring possible outbreaks and also to develop strategies regarding disease control. Although Chile has a low incidence of MTB, the current rates of migration have the potential to change this scenario. We collected and analyzed a total of 458 M. tuberculosis isolates (1 isolate per patient) originating from all 15 regions of Chile. The isolates were genotyped using the spoligotyping method and the data obtained were analyzed and compared with the SITVIT2 database. A total of 169 different patterns were identified, of which, 119 patterns (408 strains) corresponded to Spoligotype International Types (SITs) and 50 patterns corresponded to orphan strains. The most abundantly represented SITs/lineages were: SIT53/T1 (11.57%), SIT33/LAM3 (9.6%), SIT42/LAM9 (9.39%), SIT50/H3 (5.9%), SIT37/T3 (5%); analysis of the spoligotyping minimum spanning tree as well as spoligoforest were suggestive of a recent expansion of SIT42, SIT50 and SIT37; all of which potentially evolved from SIT53. The most abundantly represented lineages were LAM (40.6%), T (34.1%) and Haarlem (13.5%). LAM was more prevalent in the Santiago (43.6%) and Concepción (44.1%) isolates, rather than the Iquique (29.4%) strains. The proportion of X lineage was appreciably higher in Iquique and Concepción (11.7% in both) as compared to Santiago (1.6%). Global analysis of MTB lineage distribution in Chile versus neighboring countries showed that evolutionary recent lineages (LAM, T and Haarlem) accounted together for 88.2% of isolates in Chile, a pattern which mirrored MTB lineage distribution in neighboring countries (n = 7378 isolates recorded in SITVIT2 database for Peru, Brazil, Paraguay, and Argentina; and published studies), highlighting epidemiological advantage of Euro-American lineages in this region. Finally, we also observed exclusive emergence of patterns SIT4014/X1 and SIT4015 (unknown lineage signature) that have hitherto been found exclusively in Chile, indicating that conditions specific to Chile, along with the unique genetic makeup of the Chilean population, might have allowed for a possible co-evolution leading to the success of these emerging genotypes.

Latin-American-Mediterranean lineage of Mycobacterium tuberculosis: Human traces across pathogen's phylogeography

Currently, Mycobacterium tuberculosis isolates of Latin-American Mediterranean (LAM) family may be detected far beyond the geographic areas that coined its name 15years ago. Here, we established the framework phylogeny of this geographically intriguing and pathobiologically important mycobacterial lineage and hypothesized how human demographics and migration influenced its phylogeography. Phylogenetic analysis of LAM isolates from all continents based on 24 variable number of tandem repeats (VNTR) loci and other markers identified three global sublineages with certain geographic affinities and defined by large deletions RD115, RD174, and by spoligotype SIT33. One minor sublineage (spoligotype SIT388) appears endemic in Japan. One-locus VNTR signatures were established for sublineages and served for their search in published literature and geographic mapping. We suggest that the LAM family originated in the Western Mediterranean region. The most widespread RD115 sublineage seems the most ancient and encompasses genetically and geographically distant branches, including extremely drug resistant KZN in South Africa and LAM-RUS recently widespread across Northern Eurasia. The RD174 sublineage likely started its active spread in Brazil; its earlier branch is relatively dominated by isolates from South America and the derived one is dominated by Portuguese and South/Southeastern African isolates. The relatively most recent SIT33-sublineage is marked with enigmatic gaps and peaks across the Americas and includes South African clade F11/RD761, which likely emerged within the SIT33 subpopulation after its arrival to Africa. In addition to SIT388-sublineage, other deeply rooted, endemic LAM sublineages may exist that remain to be discovered. As a general conclusion, human mass migration appears to be the major factor that shaped the M. tuberculosis phylogeography over large time-spans.

The Association between Mycobacterium Tuberculosis Genotype and Drug Resistance in Peru

Background

The comparison of Mycobacterium tuberculosis bacterial genotypes with phenotypic, demographic, geospatial and clinical data improves our understanding of how strain lineage influences the development of drug-resistance and the spread of tuberculosis.

Methods

To investigate the association of Mycobacterium tuberculosis bacterial genotype with drug-resistance. Drug susceptibility testing together with genotyping using both 15-loci MIRU-typing and spoligotyping, was performed on 2,139 culture positive isolates, each from a different patient in Lima, Peru. Demographic, geospatial and socio-economic data were collected using questionnaires, global positioning equipment and the latest national census.

Results

The Latin American Mediterranean (LAM) clade (OR 2.4, p<0.001) was significantly associated with drug-resistance and alone accounted for more than half of all drug resistance in the region. Previously treated patients, prisoners and genetically clustered cases were also significantly associated with drug-resistance (OR's 2.5, 2.4 and 1.8, p<0.001, p<0.05, p<0.001 respectively).

Conclusions

Tuberculosis disease caused by the LAM clade was more likely to be drug resistant independent of important clinical, genetic and socio-economic confounding factors. Explanations for this include; the preferential co-evolution of LAM strains in a Latin American population, a LAM strain bacterial genetic background that favors drug-resistance or the "founder effect" from pre-existing LAM strains disproportionately exposed to drugs.

Predominant Mycobacterium tuberculosis Families and High Rates of Recent Transmission among New Cases Are Not Associated with Primary Multidrug Resistance in Lima, Peru

Sputum samples from new tuberculosis (TB) cases were collected over 2 years as part of a prospective study in the northeastern part of Lima, Peru. To measure the contribution of recent transmission to the high rates of multidrug resistance (MDR) in this area, Mycobacterium tuberculosis complex (MTBc) isolates were tested for drug susceptibility to first-line drugs and were genotyped by spoligotyping and 15-locus mycobacterial interspersed repetitive-unit (MIRU-15)-variable-number tandem repeat (VNTR) analysis. MDR was found in 6.8% of 844 isolates, of which 593 (70.3%) were identified as belonging to a known MTBc lineage, whereas 198 isolates (23.5%) could not be assigned to these lineages and 12 (1.4%) represented mixed infections. Lineage 4 accounted for 54.9% (n = 463) of the isolates, most of which belonged to the Haarlem family (n = 279). MIRU-15 analysis grouped 551/791 isolates (69.7%) in 102 clusters, with sizes ranging from 2 to 46 strains. The overall high clustering rate suggests a high level of recent transmission in this population, especially among younger patients (odds ratio [OR], 1.6; P = 0.01). Haarlem strains were more prone to cluster, compared to the other families taken together (OR, 2.0; P < 0.0001), while Beijing (OR, 0.6; P = 0.006) and LAM (OR, 0.7; P = 0.07) strains clustered less. Whereas streptomycin-resistant strains were more commonly found in clusters (OR, 1.8; P = 0.03), clustering rates did not differ between MDR and non-MDR strains (OR, 1.8; P = 0.1). Furthermore, only 16/51 MDR strains clustered with other MDR strains, suggesting that patients with primary MDR infections acquired the infections mostly from index cases outside the study population, such as retreated cases.

Mycobacterium tuberculosis phylogeography in the context of human migration and pathogen's pathobiology: Insights from Beijing and Ural families

Here, I review the population structure and phylogeography of the two contrasting families of Mycobacterium tuberculosis, Beijing and Ural, in the context of strain pathobiology and human history and migration. Proprietary database (12-loci MIRU-VNTR profiles of 3067 Beijing genotype isolates) was subjected to phylogenetic and statistical analysis. The highest rate (90%) and diversity (HGI 0.80-0.95) of the Beijing genotype in North China suggest it to be its area of origin. Under VNTR-based MDS analysis the interpopulation genetic distances correlated with geography over uninterrupted landmasses. In contrast, large water distances together with long time generated remarkable outliers. Weak and less expected affinities of the distant M. tuberculosis populations may reflect hidden epidemiological links due to unknown migration. Association with drug-resistance or increased virulence/transmissibility along with particular human migration flows shape global dissemination of some Beijing clones. The paucity of data on the Ural genotype prevents from high-resolution analysis that was mainly based on the available spoligotyping data. The North/East Pontic area marked with the highest prevalence of the Ural family may have been the area of its origin and primary dispersal in Eurasia. Ural strains are not marked by increased pathogenic capacities, increased transmissibility and association with drug resistance (but most recent reports describe an alarming increase of MDR Ural strains in some parts of eastern Europe and northwestern Russia). Large-scale SNP or WGS population-based studies targeting strains from indigenous populations and, eventually, analysis of ancient DNA will better test these hypotheses. Host genetics factors likely play the most prominent role in differential dissemination of particular M. tuberculosis genotypes.

A first insight on the population structure of Mycobacterium tuberculosis complex as studied by spoligotyping and MIRU-VNTRs in Santiago, Chile

Tuberculosis (TB) remains a significant public health problem worldwide, but the ecology of the prevalent mycobacterial strains, and their transmission, can vary depending on country and region. Chile is a country with low incidence of TB, that has a geographically isolated location in relation to the rest of South American countries due to the Andes Mountains, but recent migration from neighboring countries has changed this situation. We aimed to assess the genotypic diversity of Mycobacterium tuberculosis complex (MTBC) strains in Santiago, Chile, and compare with reports from other Latin-American countries. We analyzed MTBC isolates from pulmonary tuberculosis cases collected between years 2008 and 2013 in Central Santiago, using two genotyping methods: spoligotyping and 12-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTRs). Data obtained were analyzed and compared to the SITVIT2 database. Mean age of the patients was 47.5 years and 61% were male; 11.6% were migrants. Of 103 strains (1 isolate/patient) included, there were 56 distinct spoligotype patterns. Of these, 16 strains (15.5%) corresponded to orphan strains in the SITVIT2 database, not previously reported. Latin American and Mediterranean (LAM) (34%) and T (33%) lineages were the most prevalent strains, followed by Haarlem lineage (16.5%). Beijing family was scarcely represented with only two cases (1.9%), one of them isolated from a Peruvian migrant. The most frequent clustered spoligotypes were SIT33/LAM3 (10.7%), SIT53/T1 (8.7%), SIT50/H3 (7.8%), and SIT37/T3 (6.8%). We conclude that LAM and T genotypes are the most prevalent genotypes of MTBC in Santiago, Chile, and together correspond to almost two thirds of analyzed strains, which is similar to strain distribution reported from other countries of Latin America. Nevertheless, the high proportion of SIT37/T3, which was rarely found in other Latin American countries, may underline a specific history or demographics of Chile related to probable human migrations and evolutions.

Characterization of the genetic diversity of extensively-drug resistant Mycobacterium tuberculosis clinical isolates from pulmonary tuberculosis patients in Peru

Background

Peru holds the fourth highest burden of tuberculosis in the Americas. Despite an apparently well-functioning DOTS control program, the prevalence of multidrug resistant tuberculosis (MDR-TB) continues to increase. To worsen this situation, cases of extensively drug resistance tuberculosis (XDR-TB) have been detected. Little information exists about the genetic diversity of drug-susceptible vs. MDR-TB and XDR-TB.

Methods

Cryopreserved samples of XDR strains from 2007 to 2009 (second semester), were identified and collected. Starting from 227 frozen samples, a total of 142 XDR-TB strains of Mycobacterium tuberculosis complex (MTBC; 1 isolate per patient) were retained for this study. Each strain DNA was analyzed by spoligotyping and the 15-loci Mycobacterial Interspersed Repetitive Unit (MIRU-15).

Results

Among the 142 isolates analyzed, only 2 samples (1.41%) could not be matched to any lineage. The most prevalent sublineage was Haarlem (43.66%), followed by T (27.46%), LAM (16.2%), Beijing (9.15%), and X clade (1.41%). Spoligotype analysis identified clustering for 128/142 (90.1%) isolates vs. 49/142 (34.5%) with MIRUs. Of the samples, 90.85% belonged to retreated patients. The drug resistant profile demonstrated that 62.67% showed resistance to injectable drugs capreomycin (CAP) and kanamycin (KAN) vs. 15.5% to CAP alone and 21.8% to KAN alone. The SIT219/T1 and SIT50/H3 were the most prevalent patterns in our study. The spoligoforest analysis showed that SIT53/T1 was at the origin of many of the T lineage strains as well as a big proportion of Haarlem lineage strains (SIT50/H3, followed by SIT47/H1, SIT49/H3, and SIT2375/H1), as opposed to the SIT1/Beijing strains that did not appear to evolve into minor Beijing sublineages among the XDR-TB strains.

Conclusion

In contrast with other Latin-American countries where LAM sublineage is the most predominant, we found the Haarlem to be the most common followed by T sublineage among the XDR-TB strains.

Functional Analysis Using Whole-Genome Sequencing of a Drug-Sensitive Mycobacterium tuberculosis Strain from Peru

We report the whole-genome sequence of a Latin American-Mediterranean (LAM) lineage drug-sensitive Mycobacterium tuberculosis strain from Peru, INS-SEN. The functional analysis revealed more mutations in secondary metabolite biosynthesis, transport, and catabolism (clusters of orthologous groups [COG] category Q) than for other LAM-sensitive strains. This study contributes to the understanding of the genomic diversity of drug-sensitive M. tuberculosis.

Genetic variability of Mycobacterium tuberculosis complex in patients with no known risk factors for MDR-TB in the North-Eastern part of Lima, Peru

Background

The aim of this study was to investigate the genetic diversity among Mycobacterium tuberculosis complex circulating in patients with no known risk factors for multi-drug resistant (MDR) tuberculosis (TB) living in a high MDR burden area and analyze the relationship between genotypes, primary drug resistance and age.

Methods

Samples were collected during January-July 2009. Isolates were tested for drug susceptibility to first-line drugs and were genotyped by spoligotyping and the 15-loci Mycobacterial Interspersed Repetitive Unit (MIRU15).

Results

Among the 199 isolates analyzed, 169 (84.9%) were identified in the SpolDB4.0 and 30 (15.1%) could not be matched to any lineage. The most prevalent lineage was Haarlem (29.6%), followed by T (15.6%), Beijing (14.1%), Latin American Mediterranean (12.6%) and U (8.5%). A few isolates belonged to the X and S clades (4.5%). Spoligotype analysis identified clustering among 148 of 169 isolates, whereas with MIRU15 all isolates were unique. Out of 197 strains; 31.5% were resistant to at least one drug, 7.5% were MDR and 22.3% showed any resistance to isoniazid.

Conclusion

In contrast with other Latin-American countries where LAM lineage is the most predominant, we found the spoligotype 50 from the Haarlem lineage as the most common. None of the prevailing lineages showed a significant association with age or resistance to isoniazid and/or rifampicin.

Multidrug-resistant tuberculosis in panama is driven by clonal expansion of a multidrug-resistant Mycobacterium tuberculosis strain related to the KZN extensively drug-resistant M. tuberculosis strain from South Africa

Multidrug-resistant tuberculosis (MDR-TB) is a significant health problem in Panama. The extent to which such cases are the result of primary or acquired resistance and the strain families involved are unknown. We performed whole-genome sequencing of a collection of 66 clinical MDR isolates, along with 31 drug-susceptible isolates, that were isolated in Panama between 2001 and 2010; 78% of the MDR isolates belong to the Latin American-Mediterranean (LAM) family. Drug resistance mutations correlated well with drug susceptibility profiles. To determine the relationships among these strains and to better understand the acquisition of resistance mutations, a phylogenetic tree was constructed based on a genome-wide single-nucleotide polymorphism analysis. The phylogenetic tree shows that the isolates are highly clustered, with a single strain (LAM9-c1) accounting for nearly one-half of the MDR isolates (29/66 isolates). The LAM9-c1 strain was most prevalent among male patients of working age and was associated with high mortality rates. Members of this cluster all share identical mutations conferring resistance to isoniazid (KatG S315T mutation), rifampin (RpoB S531L mutation), and streptomycin (rrs C517T mutation). This evidence of primary resistance supports a model in which MDR-TB in Panama is driven by clonal expansion and ongoing transmission of several strains in the LAM family, including the highly successful MDR strain LAM9-c1. The phylogenetic analysis also shows that the LAM9-c1 strain is closely related to the KwaZulu-Natal (KZN) extensively drug-resistant TB strain identified in KwaZulu-Natal, South Africa. The LAM9-c1 and KZN strains likely arose from a recent common ancestor that was transmitted between Panama and South Africa and had the capacity to tolerate an accumulation of multiple resistance mutations.

Genetic diversity of Mycobacterium tuberculosis in Peru and exploration of phylogenetic associations with drug resistance

Background

There is limited available data on the strain diversity of M tuberculosis in Peru, though there may be interesting lessons to learn from a setting where multidrug resistant TB has emerged as a major problem despite an apparently well-functioning DOTS control programme.

Methods

Spoligotyping was undertaken on 794 strains of M tuberculosis collected between 1999 and 2005 from 553 community-based patients and 241 hospital-based HIV co-infected patients with pulmonary tuberculosis in Lima, Peru. Phylogenetic and epidemiologic analyses permitted identification of clusters and exploration of spoligotype associations with drug resistance.

Results

Mean patient age was 31.9 years, 63% were male and 30.4% were known to be HIV+. Rifampicin mono-resistance, isoniazid mono-resistance and multidrug resistance (MDR) were identified in 4.7%, 8.7% and 17.3% of strains respectively. Of 794 strains from 794 patients there were 149 different spoligotypes. Of these there were 27 strains (3.4%) with novel, unique orphan spoligotypes. 498 strains (62.7%) were clustered in the nine most common spoligotypes: 16.4% SIT 50 (clade H3), 12.3% SIT 53 (clade T1), 8.3% SIT 33 (LAM3), 7.4% SIT 42 (LAM9), 5.5% SIT 1 (Beijing), 3.9% SIT 47 (H1), 3.0% SIT 222 (clade unknown), 3.0% SIT1355 (LAM), and 2.8% SIT 92 (X3). Amongst HIV-negative community-based TB patients no associations were seen between drug resistance and specific spoligotypes; in contrast HIV-associated MDRTB, but not isoniazid or rifampicin mono-resistance, was associated with SIT42 and SIT53 strains.

Conclusion

Two spoligotypes were associated with MDR particularly amongst patients with HIV. The MDR-HIV association was significantly reduced after controlling for SIT42 and SIT53 status; residual confounding may explain the remaining apparent association. These data are suggestive of a prolonged, clonal, hospital-based outbreak of MDR disease amongst HIV patients but do not support a hypothesis of strain-specific propensity for the acquisition of resistance-conferring mutations.

Mycobacterium tuberculosis mutation rate estimates from different lineages predict substantial differences in the emergence of drug-resistant tuberculosis

A critical question in tuberculosis control is why some strains of Mycobacterium tuberculosis are preferentially associated with multiple drug resistances. We demonstrate that M. tuberculosis strains from Lineage 2 (East Asian lineage and Beijing sublineage) acquire drug resistances in vitro more rapidly than M. tuberculosis strains from Lineage 4 (Euro-American lineage) and that this higher rate can be attributed to a higher mutation rate. Moreover, the in vitro mutation rate correlates well with the bacterial mutation rate in humans as determined by whole genome sequencing of clinical isolates. Finally, using a stochastic mathematical model, we demonstrate that the observed differences in mutation rate predict a substantially higher probability that patients infected with a drug susceptible Lineage 2 strain will harbor multidrug resistant bacteria at the time of diagnosis. These data suggest that interventions to prevent the emergence of drug resistant tuberculosis should target bacterial as well as treatment-related risk factors.

Genotypes of Mycobacterium tuberculosis in patients at risk of drug resistance in Bolivia

Bolivia ranks among the 10 Latin American countries with the highest rates of tuberculosis (TB) and multidrug resistant (MDR) TB. In view of this, and of the lacking information on the population structure of Mycobacterium tuberculosis in the country, we explored genotype associations with drug resistance and clustering by analyzing isolates collected in 2010 from 100 consecutive TB patients at risk of drug resistance in seven of the nine departments in which Bolivia is divided. Fourteen isolates were MDR, 29 had other drug resistance profiles, and 57 were pansusceptible. Spoligotype family distribution was: Haarlem 39.4%, LAM 26.3%, T 22.2%, S 2.0%, X 1.0%, orphan 9.1%, with very low intra-family diversity and absence of Beijing genotypes. We found 66 different MIRU-VNTR patterns; the most frequent corresponded to Multiple Locus Variable Analysis (MLVA) MtbC15 patterns 860, 372 and 873. Twelve clusters, each with identical MIRU-VNTR and spoligotypes, gathered 35 patients. We found no association of genotype with drug resistant or MDR-TB. Clustering associated with SIT 50 and the H3 subfamily to which it belongs (p<0.0001). The largest cluster involved isolates from three departments and displayed a genotype (SIT 50/MLVA 860) previously identified in Bolivian migrants into Spain and Argentina suggesting that this genotype is widespread among Bolivian patients. Our study presents a first overview of M. tuberculosis genotypes at risk of drug resistance circulating in Bolivia. However, results should be taken cautiously because the sample is small and includes a particular subset of M. tuberculosis population.

Genetic diversity and transmission characteristics of Beijing family strains of Mycobacterium tuberculosis in Peru

Beijing family strains of Mycobacterium tuberculosis have attracted worldwide attention because of their wide geographical distribution and global emergence. Peru, which has a historical relationship with East Asia, is considered to be a hotspot for Beijing family strains in South America. We aimed to unveil the genetic diversity and transmission characteristics of the Beijing strains in Peru. A total of 200 Beijing family strains were identified from 2140 M. tuberculosis isolates obtained in Lima, Peru, between December 2008 and January 2010. Of them, 198 strains were classified into sublineages, on the basis of 10 sets of single nucleotide polymorphisms (SNPs). They were also subjected to variable number tandem-repeat (VNTR) typing using an international standard set of 15 loci (15-MIRU-VNTR) plus 9 additional loci optimized for Beijing strains. An additional 70 Beijing family strains, isolated between 1999 and 2006 in Lima, were also analyzed in order to make a longitudinal comparison. The Beijing family was the third largest spoligotyping clade in Peru. Its population structure, by SNP typing, was characterized by a high frequency of Sequence Type 10 (ST10), which belongs to a modern subfamily of Beijing strains (178/198, 89.9%). Twelve strains belonged to the ancient subfamily (ST3 [n=3], ST25 [n=1], ST19 [n=8]). Overall, the polymorphic information content for each of the 24 loci values was low. The 24 loci VNTR showed a high clustering rate (80.3%) and a high recent transmission index (RTI(n-1)=0.707). These strongly suggest the active and on-going transmission of Beijing family strains in the survey area. Notably, 1 VNTR genotype was found to account for 43.9% of the strains. Comparisons with data from East Asia suggested the genotype emerged as a uniquely endemic clone in Peru. A longitudinal comparison revealed the genotype was present in Lima by 1999.