Mycobacterium tuberculosis lineage 1 genetic diversity in Pará, Brazil, suggests common ancestry with east-African isolates potentially linked to historical slave trade.

Lineage 1 (L1) is one of seven Mycobacterium tuberculosis complex (MTBC) lineages. The objective of this study was to improve the complex taxonomy of L1 using phylogenetic SNPs, and to look for the origin of the main L1 sublineage prevalent in Para, Brazil. We developed a high-throughput SNPs-typing assay based on 12-L1-specific SNPs. This assay allowed us to experimentally retrieve SNP patterns on nine of these twelve SNPs in 277 isolates previously tentatively assigned to L1 spoligotyping-based sublineages. Three collections were used: Pará-Brazil (71); RIVM, the Netherlands (102), Madagascar (104). One-hundred more results were generated in Silico using the PolyTB database. Based on the final SNPs combination, the samples were classified into 11 clusters (C1-C11). Most isolates within a SNP-based cluster shared a mutual spoligotyping-defined lineage. However, L1/EAI1-SOM (SIT48) and L1/EAI6-BGD1 (SIT591) showed a poor correlation with SNP data and are not monophyletic. L1/EAI8-MDG and L1/EAI3-IND belonged to C5; this result suggests that they share a common ancestor. L1.1.3/SIT129, a spoligotype pattern found in SNPs-cluster C6, was found to be shared between Pará/Brazil and Malawi. SIT129 was independently found to be highly prevalent in Mozambique, which suggests a migration history from East-Africa to Brazil during the 16th-18th slave trade period to Northern Brazil.

Evaluation of spoligotyping, SNPs and customised MIRU-VNTR combination for genotyping Mycobacterium tuberculosis clinical isolates in Madagascar.

BACKGROUND: Combining different molecular typing methods for Mycobacterium tuberculosis complex (MTBC) can be a powerful tool for molecular epidemiology-based investigation of TB. However, the current standard method that provides high discriminatory power for such a combination, mycobacterial interspersed repetitive units-variable numbers of tandem repeats typing (MIRU-VNTR), is laborious, time-consuming and often too costly for many resource-limited laboratories. We aimed to evaluate a reduced set of loci for MIRU-VNTR typing in combination with spoligotyping and SNP-typing for routine molecular epidemiology of TB. METHOD: Spoligotyping and SNP-typing, in combination with the 15 loci MIRU-VNTR typing, were first used to type clinical MTBC isolates (n = 158) from Madagascar. A step by step reduction of MIRU-VNTR loci number was then performed according to the Hunter and Gaston Discriminatory Index (HGDI) and to the Principal component analysis (PCA) correlation with the spoligotype profiles to evaluate the discrimination power inside the generated spoligotype clusters. The 15 MIRU-VNTR was used as reference and SNP-typing was used to determine the main MTBC lineages. RESULTS: Of the 158 clinical isolates studied, the SNP-typing classified 23 into Lineage 1 (14.6%), 31 into Lineage 2 (19.6%), 23 into Lineage 3 (14.6%) and 81 into Lineage 4 strains (51.3%). 37 different spoligotypes profiles were obtained, 15 of which were unique and 20 in clusters. 15-loci MIRU-VNTR typing revealed 144 different genotypes: 132 isolates had a unique MIRU-VNTR profile and 27 isolates were grouped into 12 clusters. After a stepwise reduction of the MIRU-VNTR loci number within each main spoligotype families, three different sets composed of 5 customised MIRU-VNTR loci had a similar discrimination level to the reference 15 loci MIRU-VNTR in lineage 1, lineage 2 and lineage 3. For lineage 4, a set of 4 and 3 MIRU-VNTR loci were proposed to subtype the Harleem and LAM spoligotype families, respectively. For the T spoligotype family, a set of 5 MIRU-VNTR loci was proposed. CONCLUSION: According to the lineages and the spoligotype families, the number of MIRU-VNTR loci can be reduced to get an optimal classification of MTBC. These customized sets of MIRU-VNTR loci reduce workload and save resources while maintaining optimal discriminatory power.

Assessment of tuberculosis spatial hotspot areas in Antananarivo, Madagascar, by combining spatial analysis and genotyping.

BACKGROUND: Tuberculosis (TB) remains a public health problem in Madagascar. A crucial element of TB control is the development of an easy and rapid method for the orientation of TB control strategies in the country. Our main objective was to develop a TB spatial hotspot identification method by combining spatial analysis and TB genotyping method in Antananarivo. METHODS: Sputa of new pulmonary TB cases from 20 TB diagnosis and treatment centers (DTCs) in Antananarivo were collected from August 2013 to May 2014 for culture. Mycobacterium tuberculosis complex (MTBC) clinical isolates were typed by spoligotyping on a Luminex® 200 platform. All TB patients were respectively localized according to their neighborhood residence and the spatial distribution of all pulmonary TB patients and patients with genotypic clustered isolates were scanned respectively by the Kulldorff spatial scanning method for identification of significant spatial clustering. Areas exhibiting spatial clustering of patients with genotypic clustered isolates were considered as hotspot TB areas for transmission. RESULTS: Overall, 467 new cases were included in the study, and 394 spoligotypes were obtained (84.4%). New TB cases were distributed in 133 of the 192 Fokontany (administrative neighborhoods) of Antananarivo (1 to 15 clinical patients per Fokontany) and patients with genotypic clustered isolates were distributed in 127 of the 192 Fokontany (1 to 13 per Fokontany). A single spatial focal point of epidemics was detected when ignoring genotypic data (p = 0.039). One Fokontany of this focal point and three additional ones were detected to be spatially clustered when taking genotypes into account (p < 0.05). These four areas were declared potential TB transmission hotspots in Antananarivo and will be considered as priority targets for surveillance in the future. CONCLUSION: This method, combining spatial analysis and TB genotyping will now be used for further focused clinical and epidemiological studies in Madagascar and will allow better TB control strategies by public health authorities.

Biomarkers for risk of developing active tuberculosis in contacts of TB patients: a prospective cohort study.

Identifying those Mycobacterium tuberculosis latent-infected individuals most at risk of developing active tuberculosis (TB) using routine clinical and laboratory tests remains a huge challenge in TB control efforts. We conducted a prospective longitudinal study of clinical and laboratory markers associated with the risk of developing active TB in contacts with latent M. tuberculosis infection.HIV-negative household contacts (n=296) of pulmonary TB patients underwent monitoring of clinical features, full blood cell counts, tuberculin skin text (TST) and chest radiography performed regularly during 18 months of follow-up. Paired statistical tests, a Kaplan-Meier analysis and Cox proportional hazard modelling were performed on variables between contacts progressing or not progressing to active TB.The appearance of TB disease symptoms in contacts was significantly associated with an elevated peripheral percentage of blood monocytes (adjusted hazard ratio (aHR) 6.25, 95% CI 1.63-23.95; p<0.01), a ≥14 mm TST response (aHR 5.72, 95% CI 1.22-26.80; p=0.03) and an increased monocyte:lymphocyte ratio (aHR 4.97, 95% CI 1.3-18.99; p=0.03). Among contacts having TST ≥14 mm, a strong association with risk of progression to TB was found with an elevated blood monocyte percentage (aHR 8.46, 95% CI 1.74-41.22; p<0.01).Elevated percentage of peripheral blood monocytes plus an elevated TST response are potential biomarkers for identifying contacts of TB patients at highest risk of developing active TB.

Study of the BCG Vaccine-Induced Cellular Immune Response in Schoolchildren in Antananarivo, Madagascar.

OBJECTIVE: Although the Bacillus Calmette-Guérin vaccine (BCG) protects young children against serious forms of TB, protection against pulmonary TB is variable. We assessed BCG vaccine-induced cellular immune responses and determined for how long they could be detected during childhood in Antananarivo, Madagascar. METHODS: We assessed BCG vaccine-induced cellular immune responses by TST and IGRA (in-house ELISPOT assay) using BCG and PPD as stimulation antigen, and compared results between vaccinated and non-vaccinated schoolchildren of two age groups, 6-7 and 13-14 years old. RESULTS: Three hundred and sixty-three healthy schoolchildren were enrolled. TST was performed on 351 children and IGRA on 142. A high proportion (66%; 229/343) of the children had no TST reactivity (induration size 0 mm). TST-positive responses (≥15 mm) were more prevalent among 13-14 year-old (31.7%) than 6-7 year old (16.5%) children, both in the non-vaccinated (43% vs. 9%, p<0.001) and vaccinated (29% vs. 13%, p=0.002) subgroups. There were no significant differences in TST responses between vaccinated and non-vaccinated children in either of the age groups. The IGRA response to BCG and to PPD stimulation was not significantly different according to BCG vaccination record or to age group. A high rate (15.5%; 22/142) of indeterminate IGRA responses was observed. There was very poor agreement between TST and IGRA-PPD findings (k= 0.08) and between TST and IGRA-BCG findings (k= 0.02). CONCLUSION: Analysis of TST and IGRA response to stimulation with BCG and PPD revealed no difference in immune response between BCG-vaccinated and non-vaccinated children; also no decrease of the BCG vaccine-induced cellular immune response over time was observed. We conclude that TST and IGRA have limitations in assessing a role of BCG or tuberculosis-related immunity.

The use of microbead-based spoligotyping for Mycobacterium tuberculosis complex to evaluate the quality of the conventional method: providing guidelines for Quality Assurance when working on membranes.

BACKGROUND: The classical spoligotyping technique, relying on membrane reverse line-blot hybridization of the spacers of the Mycobacterium tuberculosis CRISPR locus, is used world-wide (598 references in Pubmed on April 8th, 2011). However, until now no inter-laboratory quality control study had been undertaken to validate this technique. We analyzed the quality of membrane-based spoligotyping by comparing it to the recently introduced and highly robust microbead-based spoligotyping. Nine hundred and twenty-seven isolates were analyzed totaling 39,861 data points. Samples were received from 11 international laboratories with a worldwide distribution. METHODS: The high-throughput microbead-based Spoligotyping was performed on CTAB and thermolyzate DNA extracted from isolated Mycobacterium tuberculosis complex (MTC) strains coming from the genotyping participating centers. Information regarding how the classical Spoligotyping method was performed by center was available. Genotype discriminatory analyses were carried out by comparing the spoligotypes obtained by both methods. The non parametric U-Mann Whitney homogeneity test and the Spearman rank correlation test were performed to validate the observed results. RESULTS: Seven out of the 11 laboratories (63%), perfectly typed more than 90% of isolates, 3 scored between 80-90% and a single center was under 80% reaching 51% concordance only. However, this was mainly due to discordance in a single spacer, likely having a non-functional probe on the membrane used. The centers using thermolyzate DNA performed as well as centers using the more extended CTAB extraction procedure. Few centers shared the same problematic spacers and these problematic spacers were scattered over the whole CRISPR locus (Mostly spacers 15, 14, 18, 37, 39, 40). CONCLUSIONS: We confirm that classical spoligotyping is a robust method with generally a high reliability in most centers. The applied DNA extraction procedure (CTAB or thermolyzate) did not affect the results in this study. However performance was center-dependent, suggesting that training is a key component in quality assurance of spoligotyping. Overall, no particular spacer yielded a higher degree of deviating results, suggesting that errors occur randomly either in the process of re-using membranes, or during the reading of the results and transferring of data from the film to a digital file. Last, the performance of the microbead-based method was excellent as previously shown by Cowan et al. (J. Clin. Microbiol. 2004) and Zhang et al. (J. Med. Microbiol. 2009) and demonstrated the proper detection of spacer 15 that is known to occasionally give weak signals in the classical spoligotyping.

Variation in gamma interferon responses to different infecting strains of Mycobacterium tuberculosis in acid-fast bacillus smear-positive patients and household contacts in Antananarivo, Madagascar.

The majority of healthy individuals exposed to Mycobacterium tuberculosis will not develop tuberculosis (TB), though many may become latently infected. More precise measurement of the human immune response to M. tuberculosis infection may help us understand this difference and potentially identify those subjects most at risk of developing active disease. Gamma interferon (IFN-gamma) production has been widely used as a proxy marker to study infection and to examine the human immune response to specific M. tuberculosis antigens. It has been suggested that genetically distinct M. tuberculosis strains may invoke different immune responses, although how these differences influence the immune responses and clinical outcome in human tuberculosis is still poorly understood. We therefore evaluated the antigen-specific IFN-gamma production responses in peripheral blood mononuclear cells from two cohorts of subjects recruited in Antananarivo, Madagascar, from 2004 to 2006 and examined the influence of the infecting M. tuberculosis strains on this response. The cohorts were sputum-positive index cases and their household contacts. Clinical strains isolated from the TB patients were typed by spoligotyping. Comparison of the IFN-gamma responses with the spoligotype of the infecting clinical strains showed that "modern" M. tuberculosis strains, like Beijing and Central Asian (CAS) strains, tended to induce lower IFN-gamma responses than "ancient" strains, like East African-Indian (EAI) strains, in index cases and their household contacts. These results suggest that new strains may have evolved to induce a host response different from that of ancient strains. These findings could have important implications in the development of therapeutic and diagnostic strategies.

Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology.

BACKGROUND: The Direct Repeat locus of the Mycobacterium tuberculosis complex (MTC) is a member of the CRISPR (Clustered regularly interspaced short palindromic repeats) sequences family. Spoligotyping is the widely used PCR-based reverse-hybridization blotting technique that assays the genetic diversity of this locus and is useful both for clinical laboratory, molecular epidemiology, evolutionary and population genetics. It is easy, robust, cheap, and produces highly diverse portable numerical results, as the result of the combination of (1) Unique Events Polymorphism (UEP) (2) Insertion-Sequence-mediated genetic recombination. Genetic convergence, although rare, was also previously demonstrated. Three previous international spoligotype databases had partly revealed the global and local geographical structures of MTC bacilli populations, however, there was a need for the release of a new, more representative and extended, international spoligotyping database. RESULTS: The fourth international spoligotyping database, SpolDB4, describes 1939 shared-types (STs) representative of a total of 39,295 strains from 122 countries, which are tentatively classified into 62 clades/lineages using a mixed expert-based and bioinformatical approach. The SpolDB4 update adds 26 new potentially phylogeographically-specific MTC genotype families. It provides a clearer picture of the current MTC genomes diversity as well as on the relationships between the genetic attributes investigated (spoligotypes) and the infra-species classification and evolutionary history of the species. Indeed, an independent Naïve-Bayes mixture-model analysis has validated main of the previous supervised SpolDB3 classification results, confirming the usefulness of both supervised and unsupervised models as an approach to understand MTC population structure. Updated results on the epidemiological status of spoligotypes, as well as genetic prevalence maps on six main lineages are also shown. Our results suggests the existence of fine geographical genetic clines within MTC populations, that could mirror the passed and present Homo sapiens sapiens demographical and mycobacterial co-evolutionary history whose structure could be further reconstructed and modelled, thereby providing a large-scale conceptual framework of the global TB Epidemiologic Network. CONCLUSION: Our results broaden the knowledge of the global phylogeography of the MTC complex. SpolDB4 should be a very useful tool to better define the identity of a given MTC clinical isolate, and to better analyze the links between its current spreading and previous evolutionary history. The building and mining of extended MTC polymorphic genetic databases is in progress.

Usefulness of restriction fragment length polymorphism and spoligotyping for epidemiological studies of Mycobacterium bovis in Madagascar: description of new genotypes.

Tuberculosis is highly prevalent in cattle in Madagascar. An epidemiological study based on genotyping of Mycobacterium bovis and its transmission to humans was carried out. The restriction fragment length polymorphism (IS6110 and DR markers) and spoligotyping were used to assess the genetic diversity of strains from different regions of Madagascar. One of these strains was isolated from goat, the other strains were isolated from zebu cattle. Nine IS6110 profiles, 20 DR profiles and 12 spoligotypes were obtained. About 90% of all isolates gave a single IS6110 band at about 1.8 kb. Most strains had the same spoligotype. M. bovis strains commonly lack spacers 39-43, and all Malagasy strains also lacked spacers 3-5, 8-10 and 16. This pattern has not been reported elsewhere. DR was the most discriminatory of the three markers. The patterns obtained with the three markers were combined to identify 34 different genotypes, one of which was found in 35% of the strains. No region-specific M. bovis genotype was identified, but the genotyping of 18 M. bovis strains isolated from patients showed that the human and bovine strains were identical, suggesting possible human contamination from zebu cattle.

A study of spoligotyping-defined Mycobacterium tuberculosis clades in relation to the origin of peopling and the demographic history in Madagascar.

Despite well-developed tuberculosis (TB) control policies in Madagascar, the incidence of TB remains high and is estimated at about 100 new cases per 100000 inhabitants. This paper describes genetic characteristics of TB bacilli in Madagascar. Using an international spoligotyping database, SpolDB4, we also attempted to identify the origin of strains circulating in Madagascar. DNA polymorphism of 333 Mycobacterium tuberculosis complex isolates was assessed. A total of 301 isolates belonging to 60 spoligotyping-defined clusters were found, whereas 32 isolates harbored orphan patterns. By comparison with the international database, we identified a new genetic group of closely genetically related M. tuberculosis strains which we suggested to be specific from Madagascar. Most of them belonging to the East-African-Indian (EAI) superfamily of strains that are responsible for 14% of total TB cases (shared types ST1514-1525). These strains are closely related to the most prevalent shared type ST109, whose distribution is mainly confined to Madagascar. The observed distribution of genotypes shows that principal genetic group 1 strains (EAI, Beijing, CAS, Afri, "Manu") is high (35.4%) suggesting an ancient evolutionary history of tuberculosis in Madagascar, in relation to the origin of peopling and the demographic history.

Snapshot of moving and expanding clones of Mycobacterium tuberculosis and their global distribution assessed by spoligotyping in an international study.

The present update on the global distribution of Mycobacterium tuberculosis complex spoligotypes provides both the octal and binary descriptions of the spoligotypes for M. tuberculosis complex, including Mycobacterium bovis, from >90 countries (13,008 patterns grouped into 813 shared types containing 11,708 isolates and 1,300 orphan patterns). A number of potential indices were developed to summarize the information on the biogeographical specificity of a given shared type, as well as its geographical spreading (matching code and spreading index, respectively). To facilitate the analysis of hundreds of spoligotypes each made up of a binary succession of 43 bits of information, a number of major and minor visual rules were also defined. A total of six major rules (A to F) with the precise description of the extra missing spacers (minor rules) were used to define 36 major clades (or families) of M. tuberculosis. Some major clades identified were the East African-Indian (EAI) clade, the Beijing clade, the Haarlem clade, the Latin American and Mediterranean (LAM) clade, the Central Asian (CAS) clade, a European clade of IS6110 low banders (X; highly prevalent in the United States and United Kingdom), and a widespread yet poorly defined clade (T). When the visual rules defined above were used for an automated labeling of the 813 shared types to define nine superfamilies of strains (Mycobacterium africanum, Beijing, M. bovis, EAI, CAS, T, Haarlem, X, and LAM), 96.9% of the shared types received a label, showing the potential for automated labeling of M. tuberculosis families in well-defined phylogeographical families. Intercontinental matches of shared types among eight continents and subcontinents (Africa, North America, Central America, South America, Europe, the Middle East and Central Asia, and the Far East) are analyzed and discussed.

Extrapulmonary and pulmonary tuberculosis in antananarivo (madagascar): high clustering rate in female patients.

Antananarivo, the capital city of Madagascar, has an endemic focus of tuberculosis (TB). We specifically studied patients with extrapulmonary TB (EPTB) and grouped patients according to infected body site. The strains were characterized by IS6110 fingerprinting and compared with those isolated from patients with pulmonary TB (PTB) during the same period in order to determine the possible association between the genotype and the clinical expression of TB. A total of 316 TB patients were included in this study: 151 individuals with EPTB, 10 with both PTB and EPTB, and 155 with PTB alone. Pleural TB was the major EPTB localization (77%) and was found more often in older patients, while PTB or EPTB in which the localization was other than pleural (other EPTB) was found in younger patients. The male-to-female ratio was slightly higher in pleural TB patients (3.06:1) than in patients with other EPTB (1.35:1). There was no significant difference in the BCG status among patients with PTB, pleural TB, and other EPTB. Analysis of IS6110 patterns showed that 167 patients (52.8%) were assigned to 37 clusters of 2 to 34 patients. Analysis of the IS6110 clusters and the IS6110 families did not show any association with a particular clinical expression of the disease. Patients with PTB or other EPTB were more likely to have strains with one IS6110 copy than patients with pleural TB. The clustering rate was found to be significantly higher in female patients (62%) than in male patients (48%) (P = 0.029), suggesting that Malagasy women were more likely to progress to disease after infection than men.

Global distribution of Mycobacterium tuberculosis spoligotypes.

We present a short summary of recent observations on the global distribution of the major clades of the Mycobacterium tuberculosis complex, the causative agent of tuberculosis. This global distribution was defined by data-mining of an international spoligotyping database, SpolDB3. This database contains 11708 patterns from as many clinical isolates originating from more than 90 countries. The 11708 spoligotypes were clustered into 813 shared types. A total of 1300 orphan patterns (clinical isolates showing a unique spoligotype) were also detected.