Use of cluster-graphs from spoligotyping data to study genotype similarities and a comparison of three indices to quantify recent tuberculosis transmission among culture positive cases in French Guiana during a eight year period

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 Structuration, Demography and Evolutionary History of Mycobacterium tuberculosis LAM9 Sublineage in the Americas as Two Distinct Subpopulations Revealed by Bayesian Analyses

Tuberculosis (TB) remains broadly present in the Americas despite intense global efforts for its control and elimination. Starting from a large dataset comprising spoligotyping (n = 21183 isolates) and 12-loci MIRU-VNTRs data (n = 4022 isolates) from a total of 31 countries of the Americas (data extracted from the SITVIT2 database), this study aimed to get an overview of lineages circulating in the Americas. A total of 17119 (80.8%) strains belonged to the Euro-American lineage 4, among which the most predominant genotypic family belonged to the Latin American and Mediterranean (LAM) lineage (n = 6386, 30.1% of strains). By combining classical phylogenetic analyses and Bayesian approaches, this study revealed for the first time a clear genetic structuration of LAM9 sublineage into two subpopulations named LAM9C1 and LAM9C2, with distinct genetic characteristics. LAM9C1 was predominant in Chile, Colombia and USA, while LAM9C2 was predominant in Brazil, Dominican Republic, Guadeloupe and French Guiana. Globally, LAM9C2 was characterized by higher allelic richness as compared to LAM9C1 isolates. Moreover, LAM9C2 sublineage appeared to expand close to twenty times more than LAM9C1 and showed older traces of expansion. Interestingly, a significant proportion of LAM9C2 isolates presented typical signature of ancestral LAM-RD^Rio MIRU-VNTR type (224226153321). Further studies based on Whole Genome Sequencing of LAM strains will provide the needed resolution to decipher the biogeographical structure and evolutionary history of this successful family.

Mycobacterium tuberculosis Genotypes Determined by Spoligotyping to Be Circulating in Colombia between 1999 and 2012 and Their Possible Associations with Transmission and Susceptibility to First-Line Drugs

INTRODUCTION

Tuberculosis (TB) remains a primary public health problem worldwide. The number of multidrug-resistant tuberculosis (MDR TB) cases has increased in recent years in Colombia. Knowledge of M. tuberculosis genotypes defined by spoligotyping can help determine the circulation of genotypes that must be controlled to prevent the spread of TB.

OBJECTIVE

To describe the genotypes of M. tuberculosis using spoligotyping in resistant and drug-sensitive isolates and their possible associations with susceptibility to first-line drugs.

METHODS

An analytical observational study was conducted that included 741 isolates of M. tuberculosis from patients. The isolates originated from 31 departments and were obtained by systematic surveillance between 1999 and 2012.

RESULTS

In total 61.94% of the isolates were resistant to 1 or more drugs, and 147 isolates were MDR. In total, 170 genotypes were found in the population structure of Colombian M. tuberculosis isolates. The isolates were mainly represented by four families: LAM (39.9%), Haarlem (19%), Orphan (17%) and T (9%). The SIT42 (LAM 9) was the most common genotype and contained 24.7% of the isolates, followed by the genotypes SIT62 (Haarlem1), SIT53 (T1), and SIT50 (H3). A high clustering of isolates was evident with 79.8% of the isolates classified into 32 groups. The Beijing family was associated with resistant isolates, whereas the Haarlem and T families were associated with sensitive isolates. The Haarlem family was also associated with grouped isolates (p = 0.031).

CONCLUSIONS

A high proportion (approximately 80%) of isolates was found in clusters; these clusters were not associated with resistance to first-line drugs. The Beijing family was associated with drug resistance, whereas the T and Haarlem families were associated with susceptibility in the Colombian isolates studied.

Estimation of Recent Transmission of Mycobacterium Tuberculosis Strains among Iranian and Afghan Immigrants: A Cluster-Based Study

BACKGROUND

Iran has extended borders with high-TB burden countries (Afghanistan and Pakistan) and immigrations of these populations influences TB distribution in the region and threatens the control strategies. The aim of this study was to evaluate the extent of recent TB transmission among Iranian and Afghan cases.

MATERIALS AND METHODS

Spoligotyping and 15-locus variable number tandem repeat (VNTR) typing were applied to genotype 102 MTB isolates (2009 to 2010). Phylogenetic relationships were analysed by two methods: a cluster-graph method and a minimum spanning tree (MST) method. Furthermore, evaluation of recent TB transmission was assessed with three indices including, RTIn, RTIn-1 and TMI.

RESULTS

Using molecular typing, 35 different spoligotypes were detected among the studied isolates. Seventy seven cases (75.4%) were distributed into 10 clusters and the remaining 25 (24.5%) isolates had a unique pattern. The cluster sizes also ranged from 2 to 21 isolates. The most frequent spoligotype in our populations belong to Haarlem (n=30, 29.4%) followed by CAS (n= 29, 28.4%) and Beijing (n=16, 15.6%) lineages. The used indices give the following values: RTIn = 0.75, RTIn-1 = 0.65 and TMI = 0.24.

CONCLUSION

The low rate of TB transmission in our findings (24%) showed that the mode of TB transmission in Iran is mostly associated with reactivation of a previous TB infection and that recently a transmitted disease has a minor role. However, the increasing incidence of the intra-community transmission in recent years highlights the need for establishing new strategies for control of TB.

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.

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.

Spoligotyping of Mycobacterium tuberculosis isolates among pulmonary tuberculosis patients in Amhara Region, Ethiopia

The aim of this study was to characterize Mycobacterium tuberculosis isolates circulating in the Amhara Region of Ethiopia. Sputum samples were collected from new pulmonary tuberculosis (TB) patients in the Region. Genotyping of mycobacterial DNA was performed by spoligotyping and isolates were assigned to families using the SpolDB4 and the model-based program 'Spotclust'. A high level of diversity was found among the 237 isolates. Sixty-five different spoligopatterns were obtained. The T (30.8%), Central Asian (CAS; 21.1%) and U (17.7%) families were the predominant isolates comprising 69.6% of the total strains. Eighty-five per cent of the U lineage belonged to Spoligo-International-Type (SIT) 910 and SIT 1729. Only a few of these strains are included in SpolDB4 to date. Of the total strains, 41 (17.3%) were unique and have not been described in SpolDB4 to date. This study indicated that the TB epidemic in Amhara Region, Ethiopia, is characterized by the circulation of numerous M. tuberculosis strain families. The high proportion of SIT 910 and SIT 1729 strains may indicate an increase in the importance of these lineages in the transmission of TB in the study region.

Molecular diversity of Mycobacterium tuberculosis isolates from patients with tuberculosis in Honduras

Background

Tuberculosis persists as a public health problem in Honduras. A better knowledge of the molecular characteristics of Mycobacterium tuberculosis strains will contribute to understand the transmission dynamics of the disease within the country. The aim of this study was to provide an insight of the genetic biodiversity of M. tuberculosis clinical isolates collected in Honduras between 1994 and 2002. Genotyping was performed using spoligotyping and RFLP. The spoligotypes obtained were compared with the SITVIT2 proprietary database of the Pasteur Institute of Guadeloupe.

Results

Spoligotyping grouped 84% of the isolates into 27 clusters (2 to 43 strains per cluster). Of the 44 shared international types (SITs) identified among the Honduran stains, 8 SITs were newly identified either within the present study or after match with an orphan type previously identified in the SITVIT2 database. In addition, 16 patterns corresponded to orphan, previously unreported isolates.

The Latin American Mediterranean (LAM) lineage was the most common in this study; 55% of the strains belonged to this family. Other genotypes found were Haarlem (16%), T (16%), X-clade (6%), Unknown signature (5%) and S (1%). Only one Beijing strain was identified (0.5%).

We observed a high degree of diversity after characterizing the 43 isolates belonging to the main spoligotyping cluster (SIT 33, LAM3) with IS6110-RFLP. A total of 35 different RFLP-fingerprints were detected, of which 6 patterns corresponded to the same number of clusters comprising 14 strains.

Conclusions

The findings obtained in this study show that tuberculosis transmission in Honduras is due to modern M. tuberculosis lineages with high level of biodiversity.

Tuberculosis and leprosy in perspective

Two of humankind's most socially and psychologically devastating diseases, tuberculosis and leprosy, have been the subject of intensive paleopathological research due to their antiquity, a presumed association with human settlement and subsistence patterns, and their propensity to leave characteristic lesions on skeletal and mummified remains. Despite a long history of medical research and the development of effective chemotherapy, these diseases remain global health threats even in the 21st century, and as such, their causative agents Mycobacterium tuberculosis and M. leprae, respectively, have recently been the subject of molecular genetics research. The new genome-level data for several mycobacterial species have informed extensive phylogenetic analyses that call into question previously accepted theories concerning the origins and antiquity of these diseases. Of special note is the fact that all new models are in broad agreement that human TB predated that in other animals, including cattle and other domesticates, and that this disease originated at least 35,000 years ago and probably closer to 2.6 million years ago. In this work, we review current phylogenetic and biogeographic models derived from molecular biology and explore their implications for the global development of TB and leprosy, past and present. In so doing, we also briefly review the skeletal evidence for TB and leprosy, explore the current status of these pathogens, critically consider current methods for identifying ancient mycobacterial DNA, and evaluate coevolutionary models.

Mycobacterium tuberculosis ecology in Venezuela: epidemiologic correlates of common spoligotypes and a large clonal cluster defined by MIRU-VNTR-24

BACKGROUND

Tuberculosis remains an endemic public health problem, but the ecology of the TB strains prevalent, and their transmission, can vary by country and by region. We sought to investigate the prevalence of Mycobacterium tuberculosis strains in different regions of Venezuela. A previous study identified the most prevalent strains in Venezuela but did not show geographical distribution nor identify clonal genotypes. To better understand local strain ecology, we used spoligotyping to analyze 1298 M. tuberculosis strains isolated in Venezuela from 1997 to 2006, predominantly from two large urban centers and two geographically distinct indigenous areas, and then studied a subgroup with MIRU-VNTR 24 loci.

RESULTS

The distribution of spoligotype families is similar to that previously reported for Venezuela and other South American countries: LAM 53%, T 10%, Haarlem 5%, S 1.9%, X 1.2%, Beijing 0.4%, and EAI 0.2%. The six most common shared types (SIT's 17, 93, 605, 42, 53, 20) accounted for 49% of the isolates and were the most common in almost all regions, but only a minority were clustered by MIRU-VNTR 24. One exception was the third most frequent overall, SIT 605, which is the most common spoligotype in the state of Carabobo but infrequent in other regions. MIRU-VNTR homogeneity suggests it is a clonal group of strains and was named the "Carabobo" genotype. Epidemiologic comparisons showed that patients with SIT 17 were younger and more likely to have had specimens positive for Acid Fast Bacilli on microscopy, and patients with SIT 53 were older and more commonly smear negative. Female TB patients tended to be younger than male patients. Patients from the high incidence, indigenous population in Delta Amacuro state were younger and had a nearly equal male:female distribution.

CONCLUSION

Six SIT's cause nearly half of the cases of tuberculosis in Venezuela and dominate in nearly all regions. Strains with SIT 17, the most common pattern overall may be more actively transmitted and SIT 53 strains may be less virulent and associated with reactivation of past infections in older patients. In contrast to other common spoligotypes, strains with SIT 605 form a clonal group centered in the state of Carabobo.

Models of deletion for visualizing bacterial variation: an application to tuberculosis spoligotypes

BACKGROUND

Molecular typing methods are commonly used to study genetic relationships among bacterial isolates. Many of these methods have become standardized and produce portable data. A popular approach for analyzing such data is to construct graphs, including phylogenies. Inferences from graph representations of data assist in understanding the patterns of transmission of bacterial pathogens, and basing these graph constructs on biological models of evolution of the molecular marker helps make these inferences. Spoligotyping is a widely used method for genotyping isolates of Mycobacterium tuberculosis that exploits polymorphism in the direct repeat region. Our goal was to examine a range of models describing the evolution of spoligotypes in order to develop a visualization method to represent likely relationships among M. tuberculosis isolates.

RESULTS

We found that inferred mutations of spoligotypes frequently involve the loss of a single or very few adjacent spacers. Using a second-order variant of Akaike's Information Criterion, we selected the Zipf model as the basis for resolving ambiguities in the ancestry of spoligotypes. We developed a method to construct graphs of spoligotypes (which we call spoligoforests). To demonstrate this method, we applied it to a tuberculosis data set from Cuba and compared the method to some existing methods.

CONCLUSION

We propose a new approach in analyzing relationships of M. tuberculosis isolates using spoligotypes. The spoligoforest recovers a plausible history of transmission and mutation events based on the selected deletion model. The method may be suitable to study markers based on loci of similar structure from other bacteria. The groupings and relationships in the spoligoforest can be analyzed along with the clinical features of strains to provide an understanding of the evolution of spoligotypes.