Identifying Mycobacterium tuberculosis complex strain families using spoligotypes

Molecular typing of Mycobacterium tuberculosis: a review of current methods, databases, softwares, and analytical tools

Studies on the epidemiology and clinical relevance of Mycobacterium tuberculosis complex (MTBC) have immensely benefited from molecular typing methods, associated software applications, and bioinformatics tools. Over the last two decades, the Pasteur Institute of Guadeloupe has developed a range of bioinformatic resources, including databases and software, to advance understanding of TB epidemiology. Traditional methods, such as IS6110-RFLP, MIRU-VNTR typing, and spoligotyping, have been instrumental but are increasingly supplanted by more precise and high-throughput techniques. These typing methods offer relatively good discrimination and reproducibility, making them popular choices for epidemiological studies. However, the advent of whole-genome sequencing (WGS) has revolutionized Mycobacterium tuberculosis complex (MTBC) typing, providing unparalleled resolution and data analysis depth. WGS enables the identification of single nucleotide polymorphisms and other genetic variations, facilitating robust phylogenetic reconstructions, and detailed outbreak investigations. This review summarizes current molecular typing methods, as well as databases and software tools used for MTBC data analysis. A comprehensive comparison of available tools and databases is provided to guide future research on the epidemiology of TB and pathogen-associated variables (drug resistance or virulence) and public health initiatives.

Clinical Strains of Mycobacterium tuberculosis Representing Different Genotype Families Exhibit Distinct Propensities to Adopt the Differentially Culturable State

Growing evidence points to the presence of differentially culturable tubercle bacteria (DCTB) in clinical specimens from individuals with active tuberculosis (TB) disease. These bacteria are unable to grow on solid media but can resuscitate in liquid media. Given the epidemiological success of certain clinical genotype families of Mycobacterium tuberculosis, we hypothesize that different strains may have distinct mechanisms of adaptation and tolerance. We used an in vitro carbon starvation model to determine the propensity of strains from lineages 2 and 4 that included the Beijing and LAM families respectively, to generate DCTB. Beijing strains were associated with a greater propensity to produce DCTB compared to LAM strains. Furthermore, LAM strains required culture filtrate (CF) for resuscitation whilst starved Beijing strains were not dependent on CF. Moreover, Beijing strains showed improved resuscitation with cognate CF, suggesting the presence of unique growth stimulatory molecules in this family. Analysis of starved Beijing and LAM strains showed longer cells, which with resuscitation were restored to a shorter length. Cell wall staining with fluorescent D-amino acids identified strain-specific incorporation patterns, indicating that cell surface remodeling during resuscitation was distinct between clinical strains. Collectively, our data demonstrate that M. tuberculosis clinical strains from different genotype lineages have differential propensities to generate DCTB, which may have implications for TB treatment success.

Strain structure analysis of Mycobacterium tuberculosis circulating among HIV negative, positive and drug resistant TB patients attending chest clinics in Western Kenya

BACKGROUND

Despite global tuberculosis (TB) interventions, the disease remains one of the major public health concerns. Kenya is ranked 15th among 22 high burden TB countries globally.

METHODS

A cross-sectional study was conducted in Western Kenya, which comprises 10 counties. A multistage sampling method was used where a single sub-county was randomly selected followed by sampling two high volume health facility from each sub-county. Identification of spoligotype profiles and their family distribution and lineage level were achieved by comparison with SITVIT database.

RESULTS

Lineage distribution pattern revealed that the most predominant lineage was CAS 220 (39.8%) followed by Beijing 128 (23.1%). The other lineages identified were T, LAM, H, X, S and MANU which were quantified as 87 (15.7%), 67 (12.1%), 16 (2.8%), 10 (1.8%), 8 (1.4%) and 5 (0.9%) respectively. CAS and Beijing strains were the most predominant lineage in both HIV negative and positive TB patients. The Beijing lineage was also the most predominant in resistant M. tuberculosis strains as compared to wild type. A total of 12 (2.0%) were orphaned M. tuberculosis strains which were spread across all the 10 counties of the study site. In multivariate logistic regression adjusting for potential cofounders three potential risk factors were significant. HIV status (OR = 1.52, CI = 0.29-3.68 and P value of 0.001), Alcohol use (OR = 0.59, CI = 0.43-3.12 and P-value =0.001) and cross border travel (OR = 0.61, CI = 0.49-3.87 and P value = 0.026). Most M. tuberculosis clinical isolates showed genetic clustering with multivariate logistic regression indicating three potential risk factors to clustering. HIV status (OR = 1.52, CI = 0.29-3.68 and P value of 0.001), Alcohol use (OR = 0.59, CI = 0.43-3.12 and P-value =0.001) and cross border travel (OR = 0.61, CI = 0.49-3.87 and P value = 0.026).

CONCLUSION

There exist diverse strains of M. tuberculosis across the 10 counties of Western Kenya. Predominant distribution of clustered genotype points to the fact that most TB cases in this region are as a result of resent transmission other than activation of latent TB.

Latent tuberculosis and computational biology: A less-talked affair

Tuberculosis (TB) is a pervasive and devastating air-borne disease caused by the organisms belonging to the Mycobacterium tuberculosis (Mtb) complex. Currently, it is the global leader in infectious disease-related death in adults. The proclivity of TB to enter the latent state has become a significant impediment to the global effort to eradicate TB. Despite decades of research, latent tuberculosis (LTB) mechanisms remain poorly understood, making it difficult to develop efficient treatment methods. In this review, we seek to shed light on the current understanding of the mechanism of LTB, with an accentuation on the insights gained through computational biology. We have outlined various well-established computational biology components, such as omics, network-based techniques, mathematical modelling, artificial intelligence, and molecular docking, to disclose the crucial facets of LTB. Additionally, we highlighted important tools and software that may be used to conduct a variety of systems biology assessments. Finally, we conclude the article by addressing the possible future directions in this field, which might help a better understanding of LTB progression.

Resistance to Second-Line Anti-TB Drugs in Cambodia: A Phenotypic and Genetic Study

Background

Due to the emergence of Mycobacterium tuberculosis (M.tb) clinical isolates resistant to most potent first-line drugs (FLD), second-line drugs (SLD) are being prescribed more frequently. We explore the genetic characteristics and molecular mechanisms of M.tb isolates phenotypically resistant to SLD, including pre-extensively drug-resistant (pre-XDR) and extensively drug-resistant (XDR) isolates.

Methods

Drug-resistant (DR) M.tb isolates collected from 2012 to 2017 were tested using sequencing and phenotypic drug susceptibility testing. Genotypes were determined to explore their links with SLD resistance patterns.

Results

Of the 272 DR M.tb isolates, 6 non-multidrug resistant (non-MDR) isolates were fluoroquinolones (FQ)-resistant, 3 were XDR and 16 were pre-XDR (14 resistant to FQ and 2 to second-line injectable drugs). The most frequent mutations in FQ-resistant and second-line injectable drugs resistant isolates were gyrA D94G (15/23) and rrs a1401g (3/5), respectively. Seventy-five percent of pre-XDR isolates and 100% of XDR isolates harbored mutations conferring resistance to pyrazinamide. All XDR isolates belonged to the Beijing genotype, of which one, named XDR+, was resistant to all drugs tested. One cluster including pre-XDR and XDR isolates was observed.

Conclusion

This is the first description of SLD resistance in Cambodia. The data suggest that the proportion of XDR and pre-XDR isolates remains low but is on the rise compared to previous reports. The characterization of the XDR+ isolate in a patient who refused treatment underlines the risk of transmission in the population. In addition, genotypic results show, as expected, that the Beijing family is the main involved in pre-XDR and XDR isolates and that the spread of the Beijing pre-XDR strain is capable of evolving into XDR strain. This study strongly indicates the need for rapid interventions in terms of diagnostic and treatment to prevent the spread of the pre-XDR and XDR strains and the emergence of more resistant ones.

Genetic Diversity and Molecular Epidemiology of Mycobacterium tuberculosis in Roraima State, Brazil

National border areas are special places for the spread of Mycobacterium tuberculosis (MTB). These regions concentrate vulnerable populations and constant population movements. Understanding the dynamics of the transmission of MTB is fundamental to propose control measures and to monitor drug resistance. We conducted a population-based prospective study of tuberculosis (TB) to evaluate molecular characteristics of MTB isolates circulating in Roraima, a state on the border of Venezuela and Guyana. Eighty isolates were genotyped by IS6110-RFLP (restriction fragment length polymorphism), spoligotyping, and 24-locus mycobacterial interspersed repetitive unit-variable number of repeats tandem (MIRU-VNTR). Drug susceptibility tests were performed by using the proportion method and GeneXpert^® MTB/RIF (Cepheid, Sunnyvale, CA). Isolates showing a phenotypic resistance profile were submitted to polymerase chain reaction (PCR) and sequencing. Spoligotyping showed 40 distinct patterns with a high prevalence of Latin-American and Mediterranean (LAM), Haarlem (H), and the "ill-defined" T clades. Mycobacterial interspersed repetitive unit -VNTR and IS6110-RFLP showed clustering rates of 21.3% and 30%, respectively. Drug resistance was detected in 11 (15.1%) isolates, and all were found to have primary resistance; among these, six (8.2%) isolates were streptomycin mono-resistant, four (5.4%) isoniazid mono-resistant, and one (1.3%) multidrug resistant. This is the first study on the molecular epidemiology and drug resistance profile of MTB from Roraima. Herein, we describe high diversity of genetic profiles circulating in this region that may be driven by the introduction of new strain types because of large population flow in this region. In summary, our results showed that analyses of these circulating strains can contribute to a better understanding of TB epidemiology in the northern Brazilian border and be useful to establish public health policies on TB prevention.

Molecular epidemiological study of multidrug-resistant tuberculosis isolated from sputum samples in Eastern Cape, South Africa

Drug-resistant tuberculosis prevalence is still a global challenge. Making it imperative to examine the molecular epidemiology of drug resistant tuberculosis. Molecular epidemiology methods can evaluate transmission patterns and risk factors, ascertain transmission cases of multidrug-resistant tuberculosis (MDR-TB) and furthermore determine transmission patterns in a human populace. This work focuses on MDR-TB isolates in distinguishing them into several species and genotyping the MDR-TB isolates, mainly for epidemiological studies using the genomic regions of difference and the spoligotyping techniques. A total of 184 deoxyribonucleic acid isolated from sputum samples that showed resistance against the two major first-line anti-tuberculosis drugs (Rifampicin and Isoniazid) were examined. The deoxyribonucleic acid samples were amplified with primers specific for each flanking region of the genomic regions of difference for the identification of different MTBC species. Isolates were further characterized into different lineages using the spoligotyping commercial kit. The M. tuberculosis species was detected in 83.7% (154/184) of the deoxyribonucleic acid isolates, followed by the M. caprae in 8.7% (16/184) and the least detected species was the M. africanum in 2.2% (4/184). Nineteen spoligotype international types (SITs) were identified in this study. The pre-existing shared types were from 94.6% (174/184) isolates with 1.1% (2/184) isolates recognized as orphans and 4.3% (8/184) isolates were not found in the SITVIT database. The predominant family (spoligotype) was the Beijing with 67.4% (124/184) strains. This study gives a general overview of drug resistant strains and the circulating strains in the Eastern Cape, South Africa and it shows that the common Mycobacteria in the province is the Beijing strain.

Purifying Selective Pressure Suggests the Functionality of a Vitamin B12 Biosynthesis Pathway in a Global Population of Mycobacterium tuberculosis

Mycobacterium tuberculosis is one of the deadliest and most challenging pathogens to study in current microbiological research. One of the issues that remains to be resolved is the importance of cobalamin in the metabolism of M. tuberculosis. The functionality of a vitamin B12 biosynthesis pathway in M. tuberculosis is under dispute, and the ability of this pathogen to scavenge vitamin B12 from the host is unknown. Here, we quantified the ratios of nonsynonymous and synonymous nucleotide substitution rates (dN/dS) in the genes involved in vitamin B12 biosynthesis and transport and in genes encoding cobalamin-dependent enzymes in nearly four thousand strains of M. tuberculosis. We showed that purifying selection is the dominant force acting on cobalamin-related genes at the levels of individual codons, genes and groups of genes. We conclude that cobalamin-related genes may not be essential but are adaptive for M. tuberculosis in clinical settings. Furthermore, the cobalamin biosynthesis pathway is likely to be functional in this species.

First insights into the genetic characteristics and drug resistance of Mycobacterium tuberculosis population collected during the first national tuberculosis prevalence survey of Lao PDR (2010-2011)

Background

In Lao People’s Democratic Republic (PDR), tuberculosis (TB) prevalence was estimated at 540/100,000 in 2011. Nevertheless, little is known about the genetic characteristics and anti-TB drug resistance of the Mycobacterium tuberculosis population. The main objective of this work was to study the genetic characteristics and drug resistance of M. tuberculosis population collected during the first National TB Prevalence Survey (TBPS) of Lao PDR (2010–2011).

Methods

Two hundred and twenty two isolates collected during TBPS (2010–2011) were analyzed with the GenoType MTBDRplus test for M. tuberculosis identification and drug resistance detection. Then, 206 of the 222 isolates were characterized by spoligotyping and MIRU-VNTR typing.

Results

Among the 222 M. tuberculosis isolates, 11 were mono-resistant to isoniazid and 2 were resistant to isoniazid and rifampicin (MDR-TB), using the GenoType MTBDRplus test. Among the 202 genetically characterized isolates, the East African-Indian (EAI) family was predominant (76.7%) followed by the Beijing (14.4%) and T (5.5%) families. EAI isolates came from all the country provinces, whereas Beijing isolates were found mainly in the northern and central provinces. A higher proportion of Beijing isolates was observed in people younger than 35 years compared to EAI. Moreover, the percentage of drug resistance was higher among Beijing (17.2%) than EAI (5.2%) isolates, and the two MDR-TB isolates belonged to the Beijing family. Combined analysis of the MIRU-VNTR and spoligotyping results (n = 202 isolates) revealed an estimated clustering rate of 11% and the occurrence of mini-outbreaks of drug-resistant TB caused by Beijing genotypes.

Conclusions

The EAI family, the ancient and endemic family in Asia, is predominant in Lao PDR whereas the prevalence of Beijing, the most harmful M. tuberculosis family for humans, is still low, differently from neighboring countries. However, its association with drug resistance, its presence in young patients and its potential association with recent transmission suggest that the Beijing family could change TB epidemiological pattern in Lao PDR. Therefore, efficient TB control and surveillance systems must be maintained and reinforced to prevent the emergence of highly transmissible and drug-resistant strains in Lao PDR, as observed in neighboring countries.

Electronic supplementary material

The online version of this article (10.1186/s12879-019-4435-z) contains supplementary material, which is available to authorized users.

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.

Extensive global movement of multidrug-resistantM. tuberculosis strains revealed by whole-genome analysis

Background

While the international spread of multidrug-resistant (MDR) Mycobacterium tuberculosis strains is an acknowledged public health threat, a broad and more comprehensive examination of the global spread of MDR-tuberculosis (TB) using whole-genome sequencing has not yet been performed.

Methods

In a global dataset of 5310 M. tuberculosis whole-genome sequences isolated from five continents, we performed a phylogenetic analysis to identify and characterise clades of MDR-TB with respect to geographic dispersion.

Results

Extensive international dissemination of MDR-TB was observed, with identification of 32 migrant MDR-TB clades with descendants isolated in 17 unique countries. Relatively recent movement of strains from both Beijing and non-Beijing lineages indicated successful global spread of varied genetic backgrounds. Migrant MDR-TB clade members shared relatively recent common ancestry, with a median estimate of divergence of 13–27 years. Migrant extensively drug-resistant (XDR)-TB clades were not observed, although development of XDR-TB within migratory MDR-TB clades was common.

Conclusions

Application of genomic techniques to investigate global MDR migration patterns revealed extensive global spread of MDR clades between countries of varying TB burden. Further expansion of genomic studies to incorporate isolates from diverse global settings into a single analysis, as well as data sharing platforms that facilitate genomic data sharing across country lines, may allow for future epidemiological analyses to monitor for international transmission of MDR-TB. In addition, efforts to perform routine whole-genome sequencing on all newly identified M. tuberculosis, like in England, will serve to better our understanding of the transmission dynamics of MDR-TB globally.

Evaluation of the 15 and 24-loci MIRU-VNTR genotyping tools with spoligotyping in the identification of Mycobacterium tuberculosis strains and their genetic diversity in molecular epidemiology studies

BACKGROUND

The transmission dynamics of Mycobacterium tuberculosis (Mtb) using various genotyping tools has been studied globally and a particular tool for genotyping Mtb is the mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR). Tuberculosis (TB) remains an important public health problem worldwide and Sri Lanka being a country of tourist destination; because of major development projects undergoing, it has a high proportion of tourists and immigrants from Asia and Europe that are characterized with highest TB incidences and drug-resistant clinical isolates. Hence, in order to address the question of Mtb genetic diversity, we investigated the discriminatory power of both MIRU-VNTR typing of 15 and 24 loci with spoligotyping to differentiate Mtb isolates.

METHOD

Acid-fast bacilli positive sputum samples (n = 150) from first visit patients were collected. Decontamination of sputum and extraction of genomic DNA were carried out using standard techniques. The isolates were characterized by MIRU-VNTR for both the 15 and 24 loci and spoligotyping.

RESULTS

In our study population, MIRU-VNTR 15 and 24 loci did not show a significant difference among the identified M. tuberculosis strains. However, MIRU 24 loci yielded an additional strain LAM, which is of T1 origin. 15 loci strain grouping had more clusters of strains grouped together while 24 loci differentiated the same cluster of strains into distinct strain types.

CONCLUSION

We conclude that the use of 15-locus MIRU-VNTR typing is sufficient for a first-line epidemiological study to genotype M. tuberculosis, but the additional discriminatory power of 24 loci MIRU-VNTR has been able to differentiate samples within highly homologous groups.

Genetic diversity of Mycobacterium tuberculosis in south coastal Karnataka, India, using spoligotyping

Background & objectives:

Despite high occurrence of tuberculosis in India very little information is available about the genetic diversity of Mycobacterium tuberculosis isolates prevailing in coastal Karnataka, India. Thus, the present study was undertaken to explore the genetic biodiversity of M. tuberculosis isolates prevailing in south coastal region of Karnataka (Udupi District), India.

Methods:

A total of 111 Mycobacterial isolates were cultured in Lowenstein Jensen (LJ) medium and after obtaining growth, DNA was extracted and spoligotyping was performed. SITVIT WEB database was used to locate families of spoligotypes.

Results:

On analyzing the hybridization results of all 111 isolates on SITVIT WEB database 57 (51.35%) isolates were clustered into 11 Spoligotype International Types (SIT). The largest cluster of 14 (12.61%) isolates was SIT-48 (EAI1-SOM), followed by SIT-1942 (CAS1-Delhi) with 11 isolates (9.9%) and SIT-11 with seven (6.30%). Moreover, 23 isolates (20.72%) had unique spoligotypes and 31 (27.92%) were orphans. Spotclust analysis revealed that majority (67%) of orphan isolates were variants of CAS (37%) and EAI-5 (34%).

Interpretation & conclusions:

The present study revealed high biodiversity among the circulating isolates of M. tuberculosis in this region with the presence of mixed genotypes earlier reported from north and south India along with certain new genotypes with unique SITs. The study highlights the need for further longitudinal studies to explore the genetic diversity and to understand the transmission dynamics of prevailing isolates.

Characterization of Mutations Conferring Resistance to Rifampin in Mycobacterium tuberculosis Clinical Strains

Resistance of Mycobacterium tuberculosis to rifampin (RMP), mediated by mutations in the rpoB gene coding for the beta-subunit of RNA polymerase, poses a serious threat to the efficacy of clinical management and, thus, control programs for tuberculosis (TB). The contribution of many individual rpoB mutations to the development and level of RMP resistance remains elusive. In this study, the incidence of mutations throughout the rpoB gene among 115 Mycobacterium tuberculosis clinical isolates, both resistant and susceptible to RMP, was determined. Of the newly discovered rpoB mutations, the role of three substitutions in the causation of RMP resistance was empirically tested. The results from in vitro mutagenesis experiments were combined with the assessment of the prevalence of rpoB mutations, and their reciprocal co-occurrences, across global M. tuberculosis populations. Twenty-two different types of mutations in the rpoB gene were identified and distributed among 58 (89.2%) RMP-resistant strains. The MICs of RMP were within the range of 40 to 800 mg/liter, with MIC₅₀ and MIC₉₀ values of 400 and 800 mg/liter, respectively. None of the mutations (Gln429His, Met434Ile, and Arg827Cys) inspected for their role in the development of RMP resistance produced an RMP-resistant phenotype in isogenic M. tuberculosis H37Rv strain-derived mutants. These mutations are supposed to compensate for fitness impairment incurred by other mutations directly associated with drug resistance.

RDRioMycobacterium tuberculosis lineage in the Brazil/Paraguay/Argentina triple border

The high tuberculosis (TB) incidence rates, the closeness of the cities and the high migration flux on the Brazil/Paraguay/Argentina border deserves an in-depth study, using Mycobacterial Interspersed Repetitive Unit (MIRU) and Spoligotyping genetic markers to explore the impact of the Mycobacterium tuberculosis RD^Rio lineage on disease transmission and resistance to anti-TB drugs in this setting. Although without the totality of M. tuberculosis isolates causing TB in this studied setting, a number of 97 isolates obtained from sputa samples culture of patients with confirmed TB, from 2013 to 2015, were submitted to 24 loci MIRU, Spoligotyping, detection of RD^Rio lineage and detection of mutation related to isoniazid and rifampicin resistance by MTBDRplus/DNA STRIP. In this sample, it was observed high clonal variability of circulating M. tuberculosis isolates causing TB in Brazilian cities bordering Paraguay and Argentina. The percentage of RD^Rio lineage causing TB in this setting was 15.46%, and lower than the detected in different areas of Brazil. According to 24 loci MIRU, the major MIRU International Type (MIT) related with RD^Rio lineage were MIT 26, MIT 738, MIT 601 with four, two and one isolates, respectively. Eight isolates with RD^Rio marker were classified as orphans. The mainly Spoligofamily related with RD^Rio lineage was LAM1 and LAM9 and no relationship between RD^Rio lineage and resistance in M. tuberculosis isolates circulating in this setting could be established. This work is pioneer in studying the dynamics of RD^Rio lineage transmission on the Brazil/Paraguay/Argentina border and deserves further studies to analyze the real contribution of the RD^Rio lineage in outbreaks and the risk of significant development of MDR-TB in the setting studied.

Risk factors associated with cluster size of Mycobacterium tuberculosis (Mtb) of different RFLP lineages in Brazil

BACKGROUND

Tuberculosis (TB) transmission is influenced by patient-related risk, environment and bacteriological factors. We determined the risk factors associated with cluster size of IS6110 RFLP based genotypes of Mycobacterium tuberculosis (Mtb) isolates from Vitoria, Espirito Santo, Brazil.

METHODS

Cross-sectional study of new TB cases identified in the metropolitan area of Vitoria, Brazil between 2000 and 2010. Mtb isolates were genotyped by the IS6110 RFLP, spoligotyping and RDRio. The isolates were classified according to genotype cluster sizes by three genotyping methods and associated patient epidemiologic characteristics. Regression Model was performed to identify factors associated with cluster size.

RESULTS

Among 959 Mtb isolates, 461 (48%) cases had an isolate that belonged to an RFLP cluster, and six clusters with ten or more isolates were identified. Of the isolates spoligotyped, 448 (52%) were classified as LAM and 412 (48%) as non-LAM. Our regression model found that 6-9 isolates/RFLP cluster were more likely belong to the LAM family, having the RDRio genotype and to be smear-positive (adjusted OR = 1.17, 95% CI 1.08-1.26; adjusted OR = 1.25, 95% CI 1.14-1.37; crude OR = 2.68, 95% IC 1.13-6.34; respectively) and living in a Serra city neighborhood decrease the risk of being in the 6-9 isolates/RFLP cluster (adjusted OR = 0.29, 95% CI, 0.10-0.84), than in the others groups. Individuals aged 21 to 30, 31 to 40 and > 50 years were less likely of belonging the 2-5 isolates/RFLP cluster than unique patterns compared to individuals < 20 years of age (adjusted OR = 0.49, 95% CI 0.28-0.85, OR = 0.43 95% CI 0.24-0.77and OR = 0. 49, 95% CI 0.26-0.91), respectively. The extrapulmonary disease was less likely to occur in those infected with strains in the 2-5 isolates/cluster group (adjustment OR = 0.45, 95% CI 0.24-0.85) than unique patterns.

CONCLUSIONS

We found that a large proportion of new TB infections in Vitoria is caused by prevalent Mtb genotypes belonging to the LAM family and RDRio genotypes. Such information demonstrates that some genotypes are more likely to cause recent transmission. Targeting interventions such as screening in specific areas and social risk groups, should be a priority for reducing transmission.

The TB Portals: an Open-Access, Web-Based Platform for Global Drug-Resistant-Tuberculosis Data Sharing and Analysis

The TB Portals program is an international consortium of physicians, radiologists, and microbiologists from countries with a heavy burden of drug-resistant tuberculosis working with data scientists and information technology professionals. Together, we have built the TB Portals, a repository of socioeconomic/geographic, clinical, laboratory, radiological, and genomic data from patient cases of drug-resistant tuberculosis backed by shareable, physical samples. Currently, there are 1,299 total cases from five country sites (Azerbaijan, Belarus, Moldova, Georgia, and Romania), 976 (75.1%) of which are multidrug or extensively drug resistant and 38.2%, 51.9%, and 36.3% of which contain X-ray, computed tomography (CT) scan, and genomic data, respectively. The top Mycobacterium tuberculosis lineages represented among collected samples are Beijing, T1, and H3, and single nucleotide polymorphisms (SNPs) that confer resistance to isoniazid, rifampin, ofloxacin, and moxifloxacin occur the most frequently. These data and samples have promoted drug discovery efforts and research into genomics and quantitative image analysis to improve diagnostics while also serving as a valuable resource for researchers and clinical providers. The TB Portals database and associated projects are continually growing, and we invite new partners and collaborations to our initiative. The TB Portals data and their associated analytical and statistical tools are freely available at https://tbportals.niaid.nih.gov/.

Application of whole genome sequence analysis to the study of Mycobacterium tuberculosis in Nunavut, Canada

Canada has one of the lowest rates of tuberculosis (TB) in the world, however, among certain sub-populations, disease incidence rates approach those observed in sub-Saharan Africa, and other high incidence regions. In this study, we applied mycobacterial interspersed repetitive unit (MIRU) variable number of tandem repeat (VNTR) and whole genome sequencing (WGS) to the analysis of Mycobacterium tuberculosis isolates obtained from Northern communities in the territory of Nunavut. WGS was carried out using the Illumina MiSeq, with identified variants used to infer phylogenetic relationships and annotated to infer functional implications. Additionally, the sequencing data from these isolates were augmented with publically available WGS to evaluate data from the Nunavut outbreak in the broader Canadian context. In this study, isolates could be classified into four major clusters by MIRU-VNTR analysis. These could be further resolved into sub-clusters using WGS. No evidence for antimicrobial resistance, either genetic or phenotypic, was observed in this cohort. Among most subjects with multiple samples, reactivation/incomplete treatment likely contributed to recurrence. However, isolates from two subjects appeared more likely to have occurred via reinfection, based on the large number of genomic single nucleotide variants detected. Finally, although quite distinct from previously reported Canadian MTB strains, isolates obtained from Nunavut clustered most closely with a cohort of samples originating in the Nunavik region of Northern Quebec. This study demonstrates the benefit of using WGS for discriminatory analysis of MTB in Canada, especially in high incidence regions. It further emphasizes the importance of focusing epidemiological intervention efforts on interrupting transmission chains of endemic TB throughout Northern communities, rather than relying on strategies applied in regions where the majority of TB cases result from importation of foreign strains.

Proteomics of Culture Filtrate of Prevalent Mycobacterium tuberculosis Strains: 2D-PAGE Map and MALDI-TOF/MS Analysis

Although diverse efforts have been done to identify biomarkers for control of tuberculosis using laboratory strain Mycobacterium tuberculosis H37Rv, the disease still poses a threat to mankind. There are many emerging M. tuberculosis strains, and proteomic profiling of these strains might be important to find out potential targets for diagnosis and/or prevention of tuberculosis. We evaluated the comparative proteomic profiling of culture filtrate (CF) proteins from prevalent M. tuberculosis strains (Central Asian or Delhi type; CAS1_Del, East African-Indian; EAI-3 and Beijing family) by 2D polyacrylamide gel electrophoresis and matrix-assisted laser desorption ionization-time-of-flight mass spectrometry. As a result, we could identify 12 CF proteins (Rv0066c, Rv1310, Rv3375, Rv1415, Rv0567, Rv1886c, Rv3803c, Rv3804c, Rv2031c, Rv1038c, Rv2809, and Rv1911c), which were consistently increased in all prevalent M. tuberculosis strains, and interestingly, two CF proteins (Rv2809, Rv1911c) were identified with unknown functions. Consistent increased intensity of these proteins suggests their critical role for survival of prevalent M. tuberculosis isolates, and some of these proteins may also have potential as diagnostic and vaccine candidates for tuberculosis, which needs to be further explored by immunological analysis.

Genotypic characterization directly applied to sputum improves the detection of Mycobacterium africanum West African 1, under-represented in positive cultures

BACKGROUND

This study aimed to compare the prevalence of Mycobacterium tuberculosis complex (MTBc) lineages between direct genotyping (on sputum) and indirect genotyping (on culture), to characterize potential culture bias against difficult growers.

METHODOLOGY/PRINCIPAL FINDINGS

Smear-positive sputa from consecutive new tuberculosis patients diagnosed in Cotonou, (Benin) were included, before patients had started treatment. An aliquot of decontaminated sputum was used for direct spoligotyping, and another aliquot was cultured on Löwenstein Jensen (LJ) medium (90 days), for indirect spoligotyping. After DNA extraction, spoligotyping was done according to the standard method for all specimens, and patterns obtained from sputa were compared versus those from the derived culture isolates. From 199 patient's sputa, 146 (73.4%) yielded a positive culture. In total, direct spoligotyping yielded a pattern in 98.5% (196/199) of the specimens, versus 73.4% (146/199) for indirect spoligotyping on cultures. There was good agreement between sputum- and isolate derived patterns: 94.4% (135/143) at spoligotype level and 96.5% (138/143) at (sub)lineage level. Two of the 8 pairs with discrepant pattern were suggestive of mixed infection in sputum. Ancestral lineages (Lineage 1, and M. africanum Lineages 5 and 6) were less likely to grow in culture (OR = 0.30, 95%CI (0.14 to 0.64), p = 0.0016); especially Lineage 5 (OR = 0.37 95%CI (0.17 to 0.79), p = 0.010). Among modern lineages, Lineage 4 was over-represented in positive-culture specimens (OR = 3.01, 95%CI (1.4 to 6.51), p = 0.005).

CONCLUSIONS/ SIGNIFICANCE

Ancestral lineages, especially M. africanum West African 1 (Lineage 5), are less likely to grow in culture relative to modern lineages, especially M. tuberculosis Euro-American (Lineage 4). Direct spoligotyping on smear positive sputum is effective and efficient compared to indirect spoligotyping of cultures. It allows for a more accurate unbiased determination of the population structure of the M. tuberculosis complex.

TRIAL REGISTRATION

ClinicalTrials.gov NCT02744469.

First Insight into a Nationwide Genotypic Diversity of Mycobacterium tuberculosis among Previously Treated Pulmonary Tuberculosis Cases in Benin, West Africa

BACKGROUND

Molecular studies on tuberculosis (TB) are rare in low-resource countries like Benin, where data on molecular study on previously treated TB cases is unavailable.

MATERIALS AND METHODS

From January to December 2014, all smear- and culture-positive previously treated pulmonary TB patients from all TB clinics were systematically recruited. Drug susceptibility testing and spoligotyping were performed on all isolates.

RESULTS

Of the 100 patients recruited, 71 (71.0%) were relapse cases and 24 (24.0%) were failure cases, while 5 (5.0%) were default cases. Resistance rate to any first-line drug was 40.0%, while 12.0% of strains were multidrug-resistant (MDR) and no strain was extensively drug-resistant (XDR). A total of 40 distinct spoligotypes were found to be corresponding to a genotypic diversity of 40.0%. ST61 was the most predominant spoligotype with prevalence of 33.0%. In all, 31 single spoligotypes and nine clusters were observed with 2 to 33 strains per cluster giving a clustering rate of 69.0%. Euro-American (Lineage 4) was the most prevalent lineage (74.0%) and Lineage 2 was associated with resistance to streptomycin.

CONCLUSION

This first insight into genetic diversity of previously treated pulmonary TB patients in Benin showed a relatively high genetic diversity of Mycobacterium tuberculosis.

Genetic diversity of Mycobacterium tuberculosis isolates from central India

Background & objectives:

There is a paucity of data available on genetic biodiversity of Mycobacterium tuberculosis isolates from central India. The present study was carried out on isolates of M. tuberculosis cultured from diagnostic clinical samples of patients from Bhopal, central India, using spoligotyping as a method of molecular typing.

Methods:

DNA was extracted from 340 isolates of M. tuberculosis from culture, confirmed as M. tuberculosis by molecular and biochemical methods and subjected to spoligotyping. The results were compared with the international SITVIT2 database.

Results:

Sixty five different spoligo international type (SIT) patterns were observed. A total of 239 (70.3%) isolates could be clustered into 25 SITs. The Central Asian (CAS) and East African Indian (EAI) families were found to be the two major circulating families in this region. SIT26/CAS1_DEL was identified as the most predominant type, followed by SIT11/EAI3_IND and SIT288/CAS2. Forty (11.8%) unique (non-clustered) and 61 (17.9%) orphan isolates were identified in the study. There was no significant association of clustering with clinical and demographic characteristics of patients.

Interpretation & conclusions:

Well established SITs were found to be predominant in our study. SIT26/CAS1_DEL was the most predominant type. However, the occurrence of a substantial number of orphan isolates may indicate the presence of active spatial and temporal evolutionary dynamics within the isolates of M. tuberculosis.

Genetic diversity, transmission dynamics and drug resistance of Mycobacterium tuberculosis in Angola

Tuberculosis (TB) poses a serious public health problem in Angola. No surveillance data on drug resistance is available and nothing is known regarding the genetic diversity and population structure of circulating Mycobacterium tuberculosis strains. Here, we have genotyped and evaluated drug susceptibility of 89 Mycobacterium tuberculosis clinical isolates from Luanda. Thirty-three different spoligotype profiles corresponding to 24 different Shared International Types (SIT) and 9 orphan profiles were detected. SIT 20 (LAM1) was the most prevalent (n = 16, 18.2%) followed by SIT 42 (LAM9; n = 15, 17.1%). Overall, the M. tuberculosis population structure in this sample was dominated by LAM (64.8%) and T (33.0%) strains. Twenty-four-loci MIRU-VNTR analysis revealed that a total of 13 isolates were grouped in 5 distinct clusters. Drug susceptibility data showed that 22 (24.7%) of the 89 clinical isolates were resistant to one or more antibacillary drugs of which 4 (4.5%) were multidrug resistant. In conclusion, this study demonstrates a high predominance of LAM strains circulating in the Luanda setting and the presence of recent transmission events. The rate and the emergence dynamics of drug resistant TB found in this sample are significant and highlight the need of further studies specifically focused on MDR-TB transmission.

Mycobacterium tuberculosis lineages and anti-tuberculosis drug resistance in reference hospitals across Viet Nam

BACKGROUND

Mycobacterium tuberculosis, the tuberculosis (TB) pathogen, despite a low level of genetic diversity, has revealed a high variety of biological and epidemiological characteristics linked to their lineages, such as transmissibility, fitness and propensity to acquire drug resistance. This has important implications for the epidemiology of TB. We conducted this first countrywide cross-sectional study to identify the prevalent M. tuberculosis lineages and to assess their epidemiological associations and their relation to drug resistance. The study was conducted among isolates acquired in reference hospitals across Vietnam. Isolates with drug susceptibility testing profiles were identified for their lineages by spoligotyping. Logistic regression was used to investigate the association of M. tuberculosis lineages with location, age and sex of the patients and drug resistance levels.

RESULTS

Results showed that the most prevalent lineage was Beijing (55.4 %), followed by EAI (27.5 %), T (6.4 %), LAM (1.3 %), Haarlem (1 %) and Zero type (0.3 %). The proportion of Beijing isolates in the North (70.4 %) and the South (68 %) was higher than in the Centre (28 %) (OR = 1.7 [95 % CI: 1.4-2.0], p < 0.0001), whereas the proportion of EAI isolates in the North (7.1 %) and the South (17 %) was much lower compared with the Centre (59 %) (OR = 0.5 [95 % CI: 0.4-0.6], p < 0.0001). Overall, Beijing isolates were the most likely to be drug-resistant and EAI isolates were the least likely to be drug-resistant, except in the South of Vietnam where EAI is also highly drug-resistant. The proportion of Beijing isolates was significantly higher (p < 0.01), and the proportion of EAI isolates was significantly lower (p < 0.05) in younger patients. The proportion of drug-resistance was higher in isolates collected from male patients and from patients in the middle age groups.

CONCLUSIONS

The findings suggest ongoing replacement of EAI lineage, which is mainly more drug-susceptible with highly drug-resistant Beijing lineage in all studied regions of Vietnam. Male patients of working ages should be the focus for better control to prevent the emergence of drug-resistant TB.

A Mycobacterial Perspective on Tuberculosis in West Africa: Significant Geographical Variation of M. africanum and Other M. tuberculosis Complex Lineages

BACKGROUND

Phylogenetically distinct Mycobacterium tuberculosis lineages differ in their phenotypes and pathogenicity. Consequently, understanding mycobacterial population structures phylogeographically is essential for design, interpretation and generalizability of clinical trials. Comprehensive efforts are lacking to date to establish the West African mycobacterial population structure on a sub-continental scale, which has diagnostic implications and can inform the design of clinical TB trials.

METHODOLOGY/PRINCIPAL FINDINGS

We collated novel and published genotyping (spoligotyping) data and classified spoligotypes into mycobacterial lineages/families using TBLineage and Spotclust, followed by phylogeographic analyses using statistics (logistic regression) and lineage axis plot analysis in GenGIS, in which a phylogenetic tree constructed in MIRU-VNTRplus was analysed. Combining spoligotyping data from 16 previously published studies with novel data from The Gambia, we obtained a total of 3580 isolates from 12 countries and identified 6 lineages comprising 32 families. By using stringent analytical tools we demonstrate for the first time a significant phylogeographic separation between western and eastern West Africa not only of the two M. africanum (West Africa 1 and 2) but also of several major M. tuberculosis sensu stricto families, such as LAM10 and Haarlem 3. Moreover, in a longitudinal logistic regression analysis for grouped data we showed that M. africanum West Africa 2 remains a persistent health concern.

CONCLUSIONS/SIGNIFICANCE

Because of the geographical divide of the mycobacterial populations in West Africa, individual research findings from one country cannot be generalized across the whole region. The unequal geographical family distribution should be considered in placement and design of future clinical trials in West Africa.

Determinants of the Sympatric Host-Pathogen Relationship in Tuberculosis

Major contributions from pathogen genome analysis and host genetics have equated the possibility of Mycobacterium tuberculosis co-evolution with its human host leading to more stable sympatric host–pathogen relationships. However, the attribution to either sympatric or allopatric categories depends on the resolution or grain of genotypic characterization. We explored the influence on the sympatric host-pathogen relationship of clinical (HIV infection and multidrug-resistant tuberculosis [MDRTB]) and demographic (gender and age) factors in regards to the genotypic grain by using spacer oligonucleotide typing (spoligotyping) for classification of M. tuberculosis strains within the Euro-American lineage. We analyzed a total of 547 tuberculosis (TB) cases, from six year consecutive sampling in a setting with high TB-HIV coinfection (32.0%). Of these, 62.0% were caused by major circulating pathogen genotypes. The sympatric relationship was defined according to spoligotype in comparison to the international spoligotype database SpolDB4. While no significant association with Euro-American lineage was observed with any of the factors analyzed, increasing the resolution with spoligotyping evidenced a significant association of MDRTB with sympatric strains, regardless of the HIV status. Furthermore, distribution curves of the prevalence of sympatric and allopatric TB in relation to patients’ age showed an accentuation of the relevance of the age of onset in the allopatric relationship, as reflected in the trimodal distribution. On the contrary, sympatric TB was characterized by the tendency towards a typical (standard) distribution curve. Our results suggest that within the Euro-American lineage a greater degree of genotyping fine-tuning is necessary in modeling the biological processes behind the host-pathogen interplay. Furthermore, prevalence distribution of sympatric TB to age was suggestive of host genetic determinisms driven by more common variants.

Genomics and Machine Learning for Taxonomy Consensus: The Mycobacterium tuberculosis Complex Paradigm

Infra-species taxonomy is a prerequisite to compare features such as virulence in different pathogen lineages. Mycobacterium tuberculosis complex taxonomy has rapidly evolved in the last 20 years through intensive clinical isolation, advances in sequencing and in the description of fast-evolving loci (CRISPR and MIRU-VNTR). On-line tools to describe new isolates have been set up based on known diversity either on CRISPRs (also known as spoligotypes) or on MIRU-VNTR profiles. The underlying taxonomies are largely concordant but use different names and offer different depths. The objectives of this study were 1) to explicit the consensus that exists between the alternative taxonomies, and 2) to provide an on-line tool to ease classification of new isolates. Genotyping (24-VNTR, 43-spacers spoligotypes, IS6110-RFLP) was undertaken for 3,454 clinical isolates from the Netherlands (2004-2008). The resulting database was enlarged with African isolates to include most human tuberculosis diversity. Assignations were obtained using TB-Lineage, MIRU-VNTRPlus, SITVITWEB and an algorithm from Borile et al. By identifying the recurrent concordances between the alternative taxonomies, we proposed a consensus including 22 sublineages. Original and consensus assignations of the all isolates from the database were subsequently implemented into an ensemble learning approach based on Machine Learning tool Weka to derive a classification scheme. All assignations were reproduced with very good sensibilities and specificities. When applied to independent datasets, it was able to suggest new sublineages such as pseudo-Beijing. This Lineage Prediction tool, efficient on 15-MIRU, 24-VNTR and spoligotype data is available on the web interface “TBminer.” Another section of this website helps summarizing key molecular epidemiological data, easing tuberculosis surveillance. Altogether, we successfully used Machine Learning on a large dataset to set up and make available the first consensual taxonomy for human Mycobacterium tuberculosis complex. Additional developments using SNPs will help stabilizing it.

Clonal diversity and drug resistance in Mycobacterium tuberculosis isolated from extra-pulmonary samples in central India--a pilot study

In India, extrapulmonary tuberculosis (EPTB) accounts for 10 - 15% of all types of tuberculosis. To identify and compare predominant spoligotypes and drug-resistance patterns in strains of Mycobacterium tuberculosis isolated from extrapulmonary and pulmonary specimens in central India, drug susceptibility testing and spoligotyping were carried out. Spoligotyping data was analyzed using SITVIT2 database. ST11/EAI3_Ind with 33% isolates among extrapulmonary specimens and ST26/CAS1_DEL with 28% isolates among pulmonary specimens were the most predominant lineages. Multidrug resistance was found in 5.5% of the strains isolated from extrapulmonary specimens in contrast to 17% isolated from pulmonary specimens.

Prevalence of Tuberculosis, Drug Susceptibility Testing, and Genotyping of Mycobacterial Isolates from Pulmonary Tuberculosis Patients in Dessie, Ethiopia

Due to their initially seemingly high cost, timely diagnosis and effective treatment of tuberculosis (TB) are usually hampered by lack or shortage of resources in many high TB burden countries. However, the benefits of effective treatment can eventually outweigh those of empirical treatment. Here, a cross-sectional study was conducted on samples from smear-positive new and retreatment TB patients. Data on sociodemographic and HIV status were collected. Samples were cultured for identification, conventional drug sensitivity testing, and molecular typing by deletion typing and spoligotyping. The results showed the youth were disproportionately affected. New cases were being treated following general treatment guidelines only. Monoresistance or multiple drug resistance was found in 16.5% of new patients. Spoligotyping showed that there were 44 patterns with families H3 and T1 (lineage 4) and CAS-Delhi (lineage 3) being dominant. Some rare patterns from lineage 7 were also found. Spoligotype pattern, HIV positivity, and previous treatment were not associated with drug resistance. That the vast majority of the patients were new cases and young and the large number of these patients with mono- or multiple drug resistance indicate that most TB cases are due to recent transmissions and that urgent actions are needed to curb the transmissions.

Genetic diversity of Mycobacterium tuberculosis isolated from tuberculosis patients in the Serengeti ecosystem in Tanzania

This study was part of a larger cross-sectional survey that was evaluating tuberculosis (TB) infection in humans, livestock and wildlife in the Serengeti ecosystem in Tanzania. The study aimed at evaluating the genetic diversity of Mycobacterium tuberculosis isolates from TB patients attending health facilities in the Serengeti ecosystem. DNA was extracted from 214 sputum cultures obtained from consecutively enrolled newly diagnosed untreated TB patients aged ≥18 years. Spacer oligonucleotide typing (spoligotyping) and Mycobacterium Interspersed Repetitive Units and Variable Number Tandem Repeat (MIRU-VNTR) were used to genotype M. tuberculosis to establish the circulating lineages. Of the214 M. tuberculosis isolates genotyped, 55 (25.7%) belonged to the Central Asian (CAS) family, 52 (24.3%) were T family (an ill-defined family), 38 (17.8%) belonged to the Latin American Mediterranean (LAM) family, 25 (11.7%) to the East-African Indian (EAI) family, 25 (11.7%) comprised of different unassigned (‘Serengeti’) strain families, while 8 (3.7%) belonged to the Beijing family. A minority group that included Haarlem, X, U and S altogether accounted for 11 (5.2%) of all genotypes. MIRU-VNTR typing produced diverse patterns within and between families indicative of unlinked transmission chains. We conclude that, in the Serengeti ecosystem only a few successful families predominate namely CAS, T, LAM and EAI families. Other types found in lower prevalence are Beijing, Haarlem, X, S and MANU. The Haarlem, EAI_Somalia, LAM3 and S/convergent and X2 subfamilies found in this study were not reported in previous studies in Tanzania.

Whole Genome Sequence of Polyresistant Mycobacterium tuberculosis CWCFVRF PRTB 19 Sputum Isolate from Chennai, India, Closely Clustering with East African Indian 5 Genogroup

We announce the draft genome sequence of a polyresistant Mycobacterium tuberculosis strain (CWCFVRF PRTB 19) isolated from the sputum of a clinically suspected tuberculosis patient, and it closely clusters to the East African Indian 5 (EAI5) lineage.

Screening for streptomycin resistance-conferring mutations in Mycobacterium tuberculosis clinical isolates from Poland

Currently, mutations in three genes, namely rrs, rpsL, and gidB, encoding 16S rRNA, ribosomal protein S12, and 16S rRNA-specific methyltransferase, respectively, are considered to be involved in conferring resistance to streptomycin (STR) in Mycobacterium tuberculosis. The aim of this study was to investigate the spectrum and frequency of these mutations in M. tuberculosis clinical isolates, both resistant and susceptible to STR. Sixty-four M. tuberculosis isolates recovered from as many TB patients from Poland in 2004 were included in the study. Within the sample were 50 multidrug-resistant (32 STR-resistant and 18 STR-susceptible) and 14 pan-susceptible isolates. Preliminary testing for STR resistance was performed with the 1% proportion method. The MICs of STR were determined by the Etest method. Mutation profiling was carried out by amplifying and sequencing the entire rrs, rpsL, and gidB genes. Non-synonymous mutations in either rrs or rpsL gene were detected in 23 (71.9%) of the STR-resistant and none of the STR-susceptible isolates. Mutations in the gidB gene were distributed among 12 (37.5%) STR-resistant and 13 (40.6%) STR-susceptible isolates. Four (12.5%) STR-resistant isolates were wild-type at all three loci examined. None of the rrs, rpsL or gidB mutations could be linked to low, intermediate or high level of STR resistance. In accordance with previous findings, the gidB 47T→G (L16R) mutation was associated with the Latin American-Mediterranean genotype family, whereas 276A→C (E92D) and 615A→G (A205A) mutations of the gidB gene were associated with the Beijing lineage. The study underlines the usefulness of rrs and rpsL mutations as molecular markers for STR resistance yet not indicative of its level. The gidB polymorphisms can serve as phylogenetic markers.

Increased level of acute phase reactants in patients infected with modern Mycobacterium tuberculosis genotypes in Mwanza, Tanzania

BACKGROUND

There is increasing evidence to suggest that different Mycobacterium tuberculosis lineages cause variations in the clinical presentation of tuberculosis (TB). Certain M. tuberculosis genotypes/lineages have been shown to be more likely to cause active TB in human populations from a distinct genetic ancestry. This study describes the genetic biodiversity of M. tuberculosis genotypes in Mwanza city, Tanzania and the clinical presentation of the disease caused by isolates of different lineages.

METHODS

Two-hundred-fifty-two isolates from pulmonary TB patients in Mwanza, Tanzania were characterized by spoligotyping, and 45 isolates were further characterized by mycobacterium interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR). The patients' level of the acute phase reactants AGP, CRP and neutrophil counts, in addition to BMI, were measured and compared to the M. tuberculosis lineage of the infectious agent for each patient.

RESULTS

The most frequent genotype was ST59 (48 out of 248 [19.4%]), belonging to the Euro-American lineage LAM11_ZWE, followed by ST21 (CAS_KILI lineage [44 out of 248 [17.7%]). A low degree of diversity (15.7% [39 different ST's out of 248 isolates]) of genotypes, in addition to a high level of mixed M. tuberculosis sub-populations among isolates with an unreported spoligotype pattern (10 out of 20 isolates [50.0%]) and isolates belonging to the ST53 lineage (13 out of 25 [52%]) was observed. Isolates of the 'modern' (TbD1-) Euro-American lineage induced higher levels of α1-acid glycoprotein (β = 0.4, P = 0.02; 95% CI [0.06-0.66]) and neutrophil counts (β = 0.9, P = 0.02; 95% CI [0.12-1.64]) and had lower BMI score (β = -1.0, P = 0.04; 95% CI[-1.89 - (-0.03)]). LAM11_ZWE ('modern') isolates induced higher levels of CRP (β = 24.4, P = 0.05; 95% CI[0.24-48.63]) and neutrophil counts (β = 0.9, P = 0.03; 95% CI[0.09-1.70]).

CONCLUSION

The low diversity of genotypes may be explained by an evolutionary advantage of the most common lineages over other lineages combined with optimal conditions for transmission, such as overcrowding and inadequate ventilation. The induction of higher levels of acute phase reactants in patients infected by 'modern' lineage isolates compared to 'ancient' lineages may suggest increased virulence among 'modern' lineage isolates.

Predicting Mycobacterium tuberculosis complex clades using knowledge-based Bayesian networks

We develop a novel approach for incorporating expert rules into Bayesian networks for classification of Mycobacterium tuberculosis complex (MTBC) clades. The proposed knowledge-based Bayesian network (KBBN) treats sets of expert rules as prior distributions on the classes. Unlike prior knowledge-based support vector machine approaches which require rules expressed as polyhedral sets, KBBN directly incorporates the rules without any modification. KBBN uses data to refine rule-based classifiers when the rule set is incomplete or ambiguous. We develop a predictive KBBN model for 69 MTBC clades found in the SITVIT international collection. We validate the approach using two testbeds that model knowledge of the MTBC obtained from two different experts and large DNA fingerprint databases to predict MTBC genetic clades and sublineages. These models represent strains of MTBC using high-throughput biomarkers called spacer oligonucleotide types (spoligotypes), since these are routinely gathered from MTBC isolates of tuberculosis (TB) patients. Results show that incorporating rules into problems can drastically increase classification accuracy if data alone are insufficient. The SITVIT KBBN is publicly available for use on the World Wide Web.

Molecular assessment, drug-resistant profile, and spacer oligonucleotide typing (spoligotyping) of Mycobacterium tuberculosis strains from Tamaulipas, México

BACKGROUND

Tuberculosis remains a serious global health problem involving one-third of the world population. A wide diversity of Mycobacterium tuberculosis strains cause about 1.5 million deaths/year worldwide, but in developing countries, the genetic diversity of M. tuberculosis strains remains largely unknown. We conducted a first insight into the population diversity of M. tuberculosis strains from Tamaulipas, Mexico.

METHODS

Seventy-two M. tuberculosis strains were identified and genetic diversity determined by spoligotyping. Drug sensibility testing and punctual mutations in inhA, ahpC, rpoB, and katG genes were assessed.

RESULTS

Spoligotyping analysis showed a higher prevalence of LAM9 > T1 > Haarlem3 subfamilies among 53 spoligotype patterns. Unexpectedly, five Beijing strains conforming four unique spoligopatterns were recovered. The more frequently isolated strains (LAM9 and T1), but none of the Beijing strains, were found resistant to INH or RIF. Also, no drug resistance was found among Haarlem3 isolates. The katG(315) gene mutation was found in 83% of INH-resistant strains, whereas rpoB(526) were associated in only 43% of RIF M. tuberculosis drug-resistant strains.

CONCLUSIONS

This and other studies report a high rate of orphan spoligotypes, which highlights the need for genotyping implementation as a routine technique for better understanding of M. tuberculosis strains in developing countries such as Mexico.

Draft Genome Sequence of a Multidrug-Resistant Clinical Isolate of Mycobacterium tuberculosis Belonging to a Novel Spoligotype

We describe the genome sequencing and analysis of a multidrug-resistant (MDR) clinical isolate of Mycobacterium tuberculosis, strain OSDD105 from India, belonging to a novel spoligotype.

Genotype diversity of Mycobacterium isolates from children in Jimma, Ethiopia

BACKGROUND

Paediatric tuberculosis (TB) is poorly addressed in Ethiopia and information about its magnitude and the genotype distribution of the causative Mycobacterium tuberculosis strains responsible for its spread are scanty.

METHODS

Gastric lavage or sputum samples were collected from consecutively enrolled TB suspect children visiting Jimma University Hospital in 2011 and cultured on Middlebrook 7H11 and Löwenstein-Jensen media. Acid fast bacterial (AFB) isolates were subjected to molecular typing targeting regions of difference (RDs), 16S rDNA gene and the direct repeat (DR) region using multiplex polymerase chain reaction (mPCR), gene sequencing and spoligotyping, respectively. Molecular drug susceptibility testing of M. tuberculosis isolates was performed by Genotype®MTBDRplus line probe assay (LPA) (Hain Life Sciences, Germany).

RESULTS

Gastric lavage (n = 43) or sputum (n = 58) samples were collected from 101 children and 31.7% (32/101) of the samples were positive for AFB by microscopy, culture and/or PCR. Out of 25 AFB isolates, 60% (15/25) were identified as M. tuberculosis by PCR, and 40% isolates (10/25) were confirmed to be non-tuberculous mycobacteria (NTM) by genus typing and 16S rDNA gene sequencing. Lineage classification assigned the M. tuberculosis strains into Euro-American (EUA, 66.7%; 10/15), East-African-Indian (EAI; 2/15), East-Asian (EA; 1/15) and Indio-Oceanic (IO; 1/15) lineages. Seven M. tuberculosis strains were new to the SpolDB4 database. All of the M. tuberculosis isolates were susceptible to isoniazid (INH) and rifampicin (RIF), except for one strain (of spoligotype SIT-149 or T3_ETH family) which had a mutation at the inhA locus which often confers resistance to INH (low level) and ethionamide.

CONCLUSIONS

Analysis of the genetic population structure of paediatric M. tuberculosis strains suggested similarity with that of adults, indicating an on-going and active transmission of M. tuberculosis from adults to children in Ethiopia. There were no multidrug-resistant TB (MDR-TB) strains among the isolates.

SNP typing reveals similarity in Mycobacterium tuberculosis genetic diversity between Portugal and Northeast Brazil

Human tuberculosis is an infectious disease caused by bacteria from the Mycobacterium tuberculosis complex (MTBC). Although spoligotyping and MIRU-VNTR are standard methodologies in MTBC genetic epidemiology, recent studies suggest that Single Nucleotide Polymorphisms (SNP) are advantageous in phylogenetics and strain group/lineages identification. In this work we use a set of 79 SNPs to characterize 1987 MTBC isolates from Portugal and 141 from Northeast Brazil. All Brazilian samples were further characterized using spolygotyping. Phylogenetic analysis against a reference set revealed that about 95% of the isolates in both populations are singly attributed to bacterial lineage 4. Within this lineage, the most frequent strain groups in both Portugal and Brazil are LAM, followed by Haarlem and X. Contrary to these groups, strain group T showed a very different prevalence between Portugal (10%) and Brazil (1.5%). Spoligotype identification shows about 10% of mis-matches compared to the use of SNPs and a little more than 1% of strains unidentifiability. The mis-matches are observed in the most represented groups of our sample set (i.e., LAM and Haarlem) in almost the same proportion. Besides being more accurate in identifying strain groups/lineages, SNP-typing can also provide phylogenetic relationships between strain groups/lineages and, thus, indicate cases showing phylogenetic incongruence. Overall, the use of SNP-typing revealed striking similarities between MTBC populations from Portugal and Brazil.

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.

Mycobacterium tuberculosis DNA fingerprint clusters and its relationship with RD(Rio) genotype in Brazil

Mycobacterium tuberculosis (Mtb) strains designated as RD(Rio) are responsible for a large cluster of new cases of tuberculosis (TB) in Rio de Janeiro. They were previously shown to be associated with severe manifestations of TB. Here, we used three genotyping methods (IS6110 RFLP, spoligotyping, and multiplex PCR) to characterize RD(Rio) and non-RD(Rio) strains from the metropolitan area of Vitória, State of Espirito Santo in southeast Brazil to determine strain diversity and transmission patterns. Strains with identical IS6110 RFLP patterns were considered to belong to a cluster indicative of recent transmission. Between 2000 and 2010, we identified 5470 new TB patients and genotyped 981 Mtb strains. Of these, 376 (38%) were RD(Rio). By RFLP, 180 (48%) of 376 RD(Rio) strains and 235 (40%) of 593 non-RD(Rio) strains belonged to RFLP cluster pattern groups (p = 0.023). Simpson's diversity index based on RFLP patterns was 0.96 for RD(Rio) and 0.98 for non-RD(Rio) strains. Thus, although RD(Rio) strains appear to be comprised of a fewer number of RFLP genotypes, they represent a heterogeneous group. While TB cases caused by RD(Rio) appear more likely to be due to recent transmission than cases caused by non-RD(Rio) strains, the difference is small. These observations suggest that factors other than inherent biological characteristic of RD(Rio) lineages are more important in determining recent transmission, and that public health measures to interrupt new transmissions need to be emphasized for TB control in Vitória.

Molecular detection of nontuberculous mycobacteria: advantages and limits of a broad-range sequencing approach

The isolation of nontuberculous mycobacteria (NTM) from clinical specimens has become very common in recent years. Such organisms are typically environmental and occasionally pathogenic for humans and animals. Standard diagnosis of mycobacterial infections relies on direct examination and culture. However, molecular tools are now available which allow quicker and more accurate diagnosis. Detection of NTM can be performed directly from clinical samples, although identification is mostly carried out after isolation. Sequencing of genomic targets (such as 16S rRNA, ITS, rpoB or hsp65) allows accurate and rapid identification, but has some technical limitations. A brief summary of the molecular methods available for NTM identification and a discussion of the problems associated with the use of sequencing analysis together with a description of available algorithms for NTM identification are the major objectives of this review.

High prevalence of Beijing and EAI4-VNM genotypes among M. tuberculosis isolates in northern Vietnam: sampling effect, rural and urban disparities

A total of 221 isolates of M. tuberculosis were sampled from hospitals and the general population in the northern plain of Vietnam, one of the most populated region of the country. Genotypic composition and diversity were characterized, and we investigated how they are affected by sampling (hospital vs. general population), correcting for potential confounding effects (location, age and gender of the patients). Spoligotyping and 12 MIRU-VNTR typing were used as first line. Then 15 MIRU-VNTR standard set was used, making 21 MIRU-VNTR typing for the clustered isolates. Result showed that 8 lineages and 13 sub-lineages were circulating in the region. The most predominant lineages were Beijing (38.5%) and EAI (38.5%). Others appeared with small proportions H (1.4%), LAM (1.8%), T (8.1%), X (0.9%), MANU (2.3%), and Zero (0.4%). Higher clustering rate was found in the hospital samples (17.9% in urban and 19.2% in rural areas) compared to the population ones (0%). The typical Vietnamese EAI4-VNM sub-lineage of EAI lineage accounted for 67% of EAI strains and was associated with older ages. Beijing genotypes were associated with younger, urban population and were characterized by high clustering rates. These characteristics strongly suggest that Beijing strains are invading the population, replacing the local EAI-VNM4, thus predicting a more serious tuberculosis situation in the future in the absence of more effective control strategies.

High prevalence of shared international type 53 among Mycobacterium tuberculosis complex strains in retreated patients from Côte d'Ivoire

BACKGROUND

Genotyping methods are useful tools to provide information on tuberculosis epidemic. They can allow a better response from health authorities and the implementation of measures for tuberculosis control. This study aimed to identify the main lineages and clades of Mycobacterium tuberculosis complex strains circulating in Côte d'Ivoire.

METHODS/MAIN FINDINGS

Strains isolated from sputum samples of patients ongoing retreatment from all the country were characterized by spoligotyping and by MIRU-VNTR. Profiles obtained by spoligotyping were first compared to the SITVIT/SpolDB4 database for family assignment. Of 194 strains analysed, 146 (75.3%) belonged to the T lineage. The most predominant spoligotype was the shared international type 53 with 135 strains (69.6%). In contrast with neighbouring countries, LAM (11 strains, 5.7%) and H (9 strains 4.6%) lineages were slightly represented. Only 3 Beijing strains (1.5%) and 4 strains of Mycobacterium africanum (2%) were found. Analysis of the results obtained with MIRU-VNTR revealed also a high level of clustering.

CONCLUSION/SIGNIFICANCE

The population of Mycobacterium tuberculosis complex strains among retreatment cases in Côte d'Ivoire exhibits a low diversity, allowing to assume recent transmission and locally based infection.

Clonal expansion of both modern and ancient genotypes of Mycobacterium tuberculosis in southern Taiwan

We present the first comprehensive analysis of Mycobacterium tuberculosis isolates circulating in the Kaohsiung region of southern Taiwan. The major spoligotypes found in the 224 isolates studied were Beijing lineages (n = 97; 43.3%), EAI lineages (n = 72; 32.1%) and Haarlem lineages (n = 18; 8.0%). By 24 MIRU-VNTR typing, 174 patterns were identified, including 24 clusters of 74 isolates and 150 unique patterns. The combination of spoligotyping and 12-MIRU-VNTR revealed that 129 (57.6%) of the 224 isolates were clustered in 18 genotypes. Moreover, 63.6% (7/11) of infected persons younger than 30 years had a Beijing strain, which could suggest recent spread among younger persons by this family of TB strains in Kaohsiung. Among the 94 Beijing family (SIT1, SIT250 and SIT1674) isolates further analyzed for SNPs by mass spectrometry, the most frequent strain found was ST10 (n = 49; 52%), followed by ST22 (n = 17; 18%) and ST19 (n = 11; 12%). Among the EAI-Manila family isolates analyzed by region deletion-based subtyping, the most frequent strain found was RD type 1 (n = 63; 87.5%), followed by RD type 2 (n = 9; 12.5%). In our previous study, the proportion of modern Beijing strains (52.5%) in northern Taiwan was significantly higher than the proportion of EAI strains (11%). In contrast, in the present study, EAI strains comprised up to 32% of Beijing strains in southern Taiwan. In conclusion, both ‘modern’ (Beijing) and ‘ancient’ (EAI) M. tuberculosis strains are prevalent in the Kaohsiung region, perhaps suggesting that both strains are somehow more adapted to southern Taiwan. It will be interesting to investigate the dynamics of the lineage composition by different selection pressures.

TB-Lineage: an online tool for classification and analysis of strains of Mycobacterium tuberculosis complex

This paper formulates a set of rules to classify genotypes of the Mycobacterium tuberculosis complex (MTBC) into major lineages using spoligotypes and MIRU-VNTR results. The rules synthesize prior literature that characterizes lineages by spacer deletions and variations in the number of repeats seen at locus MIRU24 (alias VNTR2687). A tool that efficiently and accurately implements this rule base is now freely available at http://tbinsight.cs.rpi.edu/run_tb_lineage.html. When MIRU24 data is not available, the system utilizes predictions made by a Naïve Bayes classifier based on spoligotype data. This website also provides a tool to generate spoligoforests in order to visualize the genetic diversity and relatedness of genotypes and their associated lineages. A detailed analysis of the application of these tools on a dataset collected by the CDC consisting of 3198 distinct spoligotypes and 5430 distinct MIRU-VNTR types from 37,066 clinical isolates is presented. The tools were also tested on four other independent datasets. The accuracy of automated classification using both spoligotypes and MIRU24 is >99%, and using spoligotypes alone is >95%. This online rule-based classification technique in conjunction with genotype visualization provides a practical tool that supports surveillance of TB transmission trends and molecular epidemiological studies.

Microevolution of extensively drug-resistant tuberculosis in Russia

Extensively drug-resistant (XDR) tuberculosis (TB), which is resistant to both first- and second-line antibiotics, is an escalating problem, particularly in the Russian Federation. Molecular fingerprinting of 2348 Mycobacterium tuberculosis isolates collected in Samara Oblast, Russia, revealed that 72% belonged to the Beijing lineage, a genotype associated with enhanced acquisition of drug resistance and increased virulence. Whole-genome sequencing of 34 Samaran isolates, plus 25 isolates representing global M. tuberculosis complex diversity, revealed that Beijing isolates originating in Eastern Europe formed a monophyletic group. Homoplasic polymorphisms within this clade were almost invariably associated with antibiotic resistance, indicating that the evolution of this population is primarily driven by drug therapy. Resistance genotypes showed a strong correlation with drug susceptibility phenotypes. A novel homoplasic mutation in rpoC, found only in isolates carrying a common rpoB rifampicin-resistance mutation, may play a role in fitness compensation. Most multidrug-resistant (MDR) isolates also had mutations in the promoter of a virulence gene, eis, which increase its expression and confer kanamycin resistance. Kanamycin therapy may thus select for mutants with increased virulence, helping preserve bacterial fitness and promoting transmission of drug-resistant TB strains. The East European clade was dominated by two MDR clusters, each disseminated across Samara. Polymorphisms conferring fluoroquinolone resistance were independently acquired multiple times within each cluster, indicating that XDR TB is currently not widely transmitted.