Resolving lineage assignation on Mycobacterium tuberculosis clinical isolates classified by spoligotyping with a new high-throughput 3R SNPs based method

Origin and dispersal of the Mycobacterium tuberculosis Haarlem genotype: Clues from its phylogeographic landscape and human migration

The Haarlem family belongs to the Euro-American phylogenetic lineage of Mycobacterium tuberculosis and is one of the globally spread genotypes of this important human pathogen. In spite of the sporadic observations on drug resistance and peculiar virulence profile, Haarlem remains in the shade of other M. tuberculosis genotypes. I analyzed genotyping data of the Haarlem genotype in light of its pathogenic properties and relevant human migration, to gain insight into its origin, evolutionary history, and current spread. Central Europe is marked with a very high prevalence of both major Haarlem subclades ancestral H3/SIT50 and derived H1, jointly making 33-41% in Czechia, Austria, and Hungary. There is a declining gradient of Haarlem beyond central Europe with 10-18% in Italy, France, Belgium, 10-13% in the Balkan countries and Turkey. Placing the available genetic diversity and ancient DNA data within the historical context, I hypothesize that M. tuberculosis Haarlem genotype likely originated in Central Europe and its primary long-term circulation occurred within the area of the former Austria/Austria-Hungary Empire in the 14th-19th centuries. The genotype is not highly transmissible and its spread was driven by long-term human migration. The European colonial expansion (when accompanied by a sufficient volume of migration) was a vehicle of its secondary dissemination. I conclude that human migration and its lack thereof (but not strain pathobiology) was a major driving force that shaped the population structure of this global lineage of M. tuberculosis. At the same time, Haarlem strains appear over-represented in some ethnic groups which warrants in-depth experimental research.

CRISPR-Cas systems: role in cellular processes beyond adaptive immunity

Clustered regularly interspaced short palindromic repeats and associated Cas proteins (CRISPR-Cas) are the only known adaptive immune system in prokaryotes. CRISPR-Cas system provides sequence-specific immunity against invasion by foreign genetic elements. It carries out its functions by incorporating a small part of the invading DNA sequence, termed as spacer into the CRISPR array. Although the CRISPR-Cas systems are mainly responsible for adaptive immune functions, their alternative role in the gene regulation, bacterial pathophysiology, virulence, and evolution has started to unravel. In several species, these systems are revealed to regulate the processes beyond adaptive immunity by employing various components of CRISPR-Cas machinery, independently or in combination. The molecular mechanisms entailing the regulatory processes are not clear in most of the instances. In this review, we have discussed some well-known and some recently established noncanonical functions of CRISPR-Cas system and its fast-extending applications in other biological processes.

RDRioMycobacterium tuberculosis strains associated with isoniazid resistance in Northern Mexico

BACKGROUND

It has been established that the genomic background of Mycobacterium tuberculosis may influence disease progression, in particular for the Beijing family and the Latin American and Mediterranean (LAM)/RD^Rio strains. The purpose of this study was to evaluate the prevalence of the LAM/RD^Rio genotype in cases of tuberculosis from Mexico and their drug susceptibility profile.

METHODS

Two hundred eighteen M. tuberculosis isolates were screened by 43-spacer spoligotyping. The LAM/RD^Rio genotype was confirmed by multiplex PCR, and the drug susceptibility testing was carried out in solid Löwenstein-Jensen media.

RESULTS

Among the LAM strains identified, 24 (63.1%) were confirmed as M. tuberculosis RD^Rio. All RD^Rio strains shared the RD174 deletion, that was associated with isoniazid resistance (p=0.0264).

CONCLUSIONS

We documented for the first time the isolation of the LAM/RD^Rio genotype in pulmonary cases of tuberculosis in Mexico, and we found resistance to the first-line anti-tuberculosis drug isoniazid in these strains.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSION

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

Ubiquitous and multifaceted: SIT53 spoligotype does not correlate with any particular family of Mycobacterium tuberculosis

Euro-American lineage (Lineage 4) of Mycobacterium tuberculosis comprises genetically and geographically diverse families that differ in their clinical and/or epidemiological capacities. Due to the characteristic structure of the CRISPR locus (presence of almost all 43 classical spacers except for deleted signals 33-36), spoligotype SIT53 takes the basal position in the evolution of this lineage. In the SITVIT database, it is assigned to the "ill-defined" T family and T1 subfamily. Here, I analyzed the phylogenetic diversity of SIT53 isolates and discussed interconnected terminological issues concerning M. tuberculosis population structure. The 24-MIRU-VNTR profiles of 266 SIT53 isolates from Europe, Asia, Africa, and South America were submitted to the phylogenetic analysis jointly with reference profiles of different families from MIRU-VNTRplus database. Under this analysis, SIT53 isolates were clustered within different and distant families such as Ghana, NEW-1 (L4.5), TUR (L4.2.2.1), etc whereas many remained unclassified within L4. This confirms the evolutionarily basal position of this spoligotype and in its turn, this demonstrates that SIT53 does not correspond to any particular family of M. tuberculosis. Instead, different SIT53 subpopulations with evolutionarily stable and unchanged CRISPR locus gave rise to different and distant families that in many instances evolved through long-term allopatric evolution.

Mycobacterium tuberculosis complex in wildlife: Review of current applications of antemortem and postmortem diagnosis

Tuberculosis (TB) is a chronic inflammatory and zoonotic disease caused by Mycobacterium tuberculosis complex (MTBC) members, which affects various domestic animals, wildlife, and humans. Some wild animals serve as reservoir hosts in the transmission and epidemiology of the disease. Therefore, the monitoring and surveillance of both wild and domestic hosts are critical for prevention and control strategies. For TB diagnosis, the single intradermal tuberculin test or the single comparative intradermal tuberculin test, and the gamma-interferon test, which is regarded as an ancillary test, are used. Postmortem examination can identify granulomatous lesions compatible with a diagnosis of TB. In contrast, smears of the lesions can be stained for acid-fast bacilli, and samples of the affected organs can be subjected to histopathological analyses. Culture is the gold standard test for isolating mycobacterial bacilli because it has high sensitivity and specificity compared with other methods. Serology for antibody detection allows the testing of many samples simply, rapidly, and inexpensively, and the protocol can be standardized in different laboratories. Molecular biological analyses are also applicable to trace the epidemiology of the disease. In conclusion, reviewing the various techniques used in MTBC diagnosis can help establish guidelines for researchers when choosing a particular diagnostic method depending on the situation at hand, be it disease outbreaks in wildlife or for epidemiological studies. This is because a good understanding of various diagnostic techniques will aid in monitoring and managing emerging pandemic threats of infectious diseases from wildlife and also preventing the potential spread of zoonotic TB to livestock and humans. This review aimed to provide up-to-date information on different techniques used for diagnosing TB at the interfaces between wildlife, livestock, and humans.

RDRioMycobacterium tuberculosis strains associated with isoniazid resistance in Northern Mexico

BACKGROUND

It has been established that the genomic background of Mycobacterium tuberculosis may influence disease progression, in particular for the Beijing family and the Latin American and Mediterranean (LAM)/RDRio strains. The purpose of this study was to evaluate the prevalence of the LAM/RDRio genotype in cases of tuberculosis from Mexico and their drug susceptibility profile.

METHODS

Two hundred eighteen M. tuberculosis isolates were screened by 43-spacer spoligotyping. The LAM/RDRio genotype was confirmed by multiplex PCR, and the drug susceptibility testing was carried out in solid Löwenstein-Jensen media.

RESULTS

Among the LAM strains identified, 24 (63.1%) were confirmed as M. tuberculosis RDRio. All RDRio strains shared the RD174 deletion, that was associated with isoniazid resistance (p=0.0264).

CONCLUSIONS

We documented for the first time the isolation of the LAM/RDRio genotype in pulmonary cases of tuberculosis in Mexico, and we found resistance to the first-line anti-tuberculosis drug isoniazid in these strains.

Most LAM Mycobacterium tuberculosis strains in Venezuela, but not SIT605, belong to the RDRio subfamily

Tuberculosis is a global public health problem that is resurgent in Venezuela, with 13 thousand estimated new cases in 2018. Strains of the Mycobacterium tuberculosis RD^Rio, subfamily belong to the Latín American Mediterranean (LAM) family and are a major cause of TB in Rio de Janeiro, Brazil. LAM strains predominate in Venezuela, where spoligotype SIT605 is common, but surprisingly not found elsewhere. We sought to assess the presence of RD^Rio strains in tuberculosis patients in different regions of Venezuela and determine whether SIT605 also belongs to the RD^Rio subfamily. Using spoligotyping and MIRU-VNTR 24 loci, we identified 86 clinical LAM and SIT605 isolates from the Venezuelan capital Caracas and several Venezuelan states. Region of difference deletion loci RD174 and RD^Rio, and also IS1561 were used to identify strains of the RD^Rio subfamily, while IS6110 at position 932,204 and the Ag85C¹⁰³ polymorphism were used to validate SIT 605 as a LAM family strain. We found that 69.8% of the isolates were RD^Río, including 94.3% of strains isolated in Caracas, 17.9% isolated in the state of Carabobo, the two strains analyzed from Delta Amacuro, and one each from Sucre, Apure and Aragua states. RD^Rio was in 100% of: SIT17 (LAM 2); SIT20 (LAM 1); SITs 93, 1694, 1696, 960, 1367 (LAM 5); and SITs 216 (LAM 9); but only 75% of SIT42 (LAM 9) strains. Thus, most of the LAM strains in Venezuela belong to the RD^Río subfamily. SIT 605 strains, although LAM, are not in the RD^Río subfamily.

Mycobacterium bovis prevalence affects the performance of a commercial serological assay for bovine tuberculosis in African buffaloes

The endemic presence of bovine tuberculosis (BTB) in African buffaloes in South Africa has severe consequences for BTB control in domestic cattle, buffalo ranching and wildlife conservation, and poses a potential risk to public health. This study determined the BTB prevalence in free-ranging buffaloes in two game reserves and assessed the influence of the prevalence of mycobacterial infections on the performance of a commercial cattle-specific serological assay for BTB (TB ELISA). Buffaloes (n = 997) were tested with the tuberculin skin test and TB ELISA; a subset (n = 119) was tested longitudinally. Culture, PCR and sequencing were used to confirm infection with M. bovis and/or non-tuberculous mycobacteria (NTM). Prevalence of BTB, but not NTM, influenced the TB ELISA performance. Multiple testing did not increase test confidence. The findings strongly illustrate the need for development of novel assays that can supplement existing assays for a more comprehensive testing scheme for BTB in African buffaloes.

CRISPR-cas system: biological function in microbes and its use to treat antimicrobial resistant pathogens

The development of antibiotic resistance in bacteria is a major public health threat. Infection rates of resistant pathogens continue to rise against nearly all antimicrobials, which has led to development of different strategies to combat the antimicrobial resistance. In this review, we discuss how the newly popular CRISPR-cas system has been applied to combat antibiotic resistance in both extracellular and intracellular pathogens. We also review a recently developed method in which nano-size CRISPR complex was used without any phage to target the mecA gene. However, there is still challenge to practice these methods in field against emerging antimicrobial resistant pathogens.

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.

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.

Genetic diversity of Mycobacterium tuberculosis complex strains isolated from livestock workers and cattle in Nigeria

Molecular typing techniques are useful in understanding tuberculosis epidemiology; yet, they have been under-utilised at the human-animal interface in Nigeria. Sixty-four Mycobacterium tuberculosis complex (MTBC) isolates including 42 M. tuberculosis, 13 M. bovis and nine M. africanum obtained from livestock workers (LW, n = 47) and their cattle (n = 17) in three geographical zones of Nigeria were genotyped to identify and evaluate the genetic diversity of the circulating MTBC using spoligotyping. Distribution into clades of M. tuberculosis revealed; 45.3% Uganda I- [SIT46- cattle: 1; LW: 28], 14.1% Latin American Mediterranean- [SIT61, cattle: 1; LW: 8], and 1.6% T- [SIT53—LW: 1]. The M. bovis strains were 6.3% SB0944 [cattle: 4] and 1.6% each of SB0300, SB1026, SB1027 and SB1439 [cattle: 4]. Seventeen MTBC isolates [cattle: 7; LW: 10] yielded 14 new spoligotype patterns including three M. tuberculosis strains (three isolates), five M. bovis strains (five isolates) and six M. africanum strains (nine isolates), two of which belonged to MAF1. Only few families namely, the not previously described Uganda I-, LAM and SB0944 are predominant among the LW and cattle, with other types in lower prevalences. The strain population structure indicates an intriguing diversity and possible zoonotic linkage with consequences for TB control in the country. The need to employ newer molecular techniques such as Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeats and whole genome sequence to decipher circulating MTBC strains in Nigeria is advocated.

Molecular drug susceptibility testing and strain typing of tuberculosis by DNA hybridization

In Mycobacterium tuberculosis (Mtb) the detection of single nucleotide polymorphisms (SNPs) is of high importance both for diagnostics, since drug resistance is primarily caused by the acquisition of SNPs in multiple drug targets, and for epidemiological studies in which strain typing is performed by SNP identification. To provide the necessary coverage of clinically relevant resistance profiles and strain types, nucleic acid-based measurement techniques must be able to detect a large number of potential SNPs. Since the Mtb problem is pressing in many resource-poor countries, requiring low-cost point-of-care biosensors, this is a non-trivial technological challenge. This paper presents a proof-of-concept in which we chose simple DNA-DNA hybridization as a sensing principle since this can be transferred to existing low-cost hardware platforms, and we pushed the multiplex boundaries of it. With a custom designed probe set and a physicochemical-driven data analysis it was possible to simultaneously detect the presence of SNPs associated with first- and second-line drug resistance and Mtb strain typing. We have demonstrated its use for the identification of drug resistance and strain type from a panel of phylogenetically diverse clinical strains. Furthermore, reliable detection of the presence of a minority population (<5%) of drug-resistant Mtb was possible.

Macro-geographical specificities of the prevailing tuberculosis epidemic as seen through SITVIT2, an updated version of the Mycobacterium tuberculosis genotyping database

In order to provide a global overview of genotypic, epidemiologic, demographic, phylogeographical, and drug resistance characteristics related to the prevailing tuberculosis (TB) epidemic, we hereby report an update of the 6th version of the international genotyping database SITVIT2. We also make all the available information accessible through a dedicated website (available at http://www.pasteur-guadeloupe.fr:8081/SITVIT2). Thanks to the public release of SITVIT2 which is currently the largest international multimarker genotyping database with a compilation of 111,635 clinical isolates from 169 countries of patient origin (131 countries of isolation, representing 1032 cities), our major aim is to highlight macro- and micro-geographical cleavages and phylogeographical specificities of circulating Mycobacterium tuberculosis complex (MTBC) clones worldwide. For this purpose, we retained strains typed by the most commonly used PCR-based methodology for TB genotyping, i.e., spoligotyping based on the polymorphism of the direct repeat (DR) locus, 5-loci Exact Tandem Repeats (ETRs), and MIRU-VNTR minisatellites used in 12-, 15-, or 24-loci formats. We describe the SITVIT2 database and integrated online applications that permit to interrogate the database using easy drop-down menus to draw maps, graphics and tables versus a long list of parameters and variables available for individual clinical isolates (year and place of isolation, origin, sex, and age of patient, drug-resistance, etc.). Available tools further allow to generate phylogenetical snapshot of circulating strains as Lineage-specific WebLogos, as well as minimum spanning trees of their genotypes in conjunction with their geographical distribution, drug-resistance, demographic, and epidemiologic characteristics instantaneously; whereas online statistical analyses let a user to pinpoint phylogeographical specificities of circulating MTBC lineages and conclude on actual demographic trends. Available associated information on gender (n = 18,944), age (n = 16,968), drug resistance (n = 19,606), and HIV serology (n = 2673), allowed to draw some important conclusions on TB geo-epidemiology; e.g. a positive correlation exists between certain Mycobacterium tuberculosis lineages (such as CAS and Beijing) and drug resistance (p-value<.001), while other lineages (such as LAM, X, and BOV) are more frequently associated with HIV-positive serology (p-value<.001). Besides, availability of information on the year of isolation of strains (range 1759-2012), also allowed to make tentative correlations between drug resistance information and lineages - portraying probable evolution trends over time and space. To conclude, the present approach of geographical mapping of predominant clinical isolates of tubercle bacilli causing the bulk of the disease both at country and regional level in conjunction with epidemiologic and demographic characteristics allows to shed new light on TB geo-epidemiology in relation with the continued waves of peopling and human migration.

Pathogenomics of Virulence Traits of Plesiomonas shigelloides That Were Deemed Inconclusive by Traditional Experimental Approaches

One of the major challenges of modern medicine includes the failure of conventional protocols to characterize the pathogenicity of emerging pathogens. This is particularly apparent in the case of Plesiomonas shigelloides. Although a number of infections have been linked to this microorganism, experimental evidence of its virulence factors (VFs), obtained by traditional approaches, is somewhat inconclusive. Hence, it remains unclear whether P. shigelloides is a true or opportunistic one. In the current study, four publicly available whole-genome sequences of P. shigelloides (GN7, NCTC10360, 302-73, and LS1) were profiled using bioinformatics platforms to determine the putative candidate VFs to characterize the bacterial pathogenicity. Overall, 134 unique open reading frames (ORFs) were identified that were homologous or orthologous to virulence genes identified in other pathogens. Of these, 52.24% (70/134) were jointly shared by the strains. The numbers of strain-specific virulence traits were 4 in LS1; 7 in NCTC10360; 10 in 302-73; and 15 in GN7. The pathogenicity islands (PAIs) common to all the strains accounted for 24.07% ORFs. The numbers of PAIs exclusive to each strain were 8 in 302-73; 11 in NCTC10360; 14 in GN7; and 18 in LS1. A PAI encoding Vibrio cholerae ToxR-activated gene d protein was specific to 302-73, GN7, and NCTC10360 strains. Out of 33 antibiotic multi-resistance genes identified, 16 (48.48%) genes were intrinsic to all strains. Further, 17 (22.08%) of 77 antibiotic resistance islands were found in all the strains. Out of 23 identified distinct insertion sequences, 13 were only harbored by strain LS1. The number of intact prophages identified in the strains was 1 in GN7; 2 in 302-73; and 2 in NCTC10360. Further, 1 CRISPR element was identified in LS1; 2 in NCTC10360; and 8 in 302-73. Fifteen (78.95%) of 19 secretion systems and secretion effector variants were identified in all the strains. In conclusion, certain P. shigelloides strains might possess VFs associated with gastroenteritis and extraintestinal infections. However, the role of host factors in the onset of infections should not be undermined.

Genotypic diversity of Mycobacterium tuberculosis in Buenos Aires, Argentina

Buenos Aires is an overpopulated port city historically inhabited by people of European descent. Together with its broader metropolitan area, the city exhibits medium tuberculosis rates, and receives migrants, mainly from tuberculosis highly endemic areas of Argentina and neighboring countries. This work was aimed to gain insight into the Mycobacterium tuberculosis population structure in two suburban districts of Buenos Aires which are illustrative of the overall situation of tuberculosis in Argentina. The Lineage 4 Euro-American accounted for >99% of the 816 isolates analyzed (one per patient). Frequencies of spoligotype families were T 35.9%, LAM 33.2%, Haarlem 19.5%, S 3.2%, X 1.5%, Ural 0.7%, BOV 0.2%, Beijing 0.2%, and Cameroon 0.2%. Unknown signatures accounted for 5.3% isolates. Of 55 spoligotypes not matching any extant shared international type (SIT) in SITVIT database, 22 fitted into 15 newly-issued SITs. Certain autochthonous South American genotypes were found to be actively evolving. LAM3, which is wild type for RD^rio, was the predominant LAM subfamily in both districts and the RD^rio signature was rare among autochthonous, newly created, SITs and orphan patterns. Two genotypes that are rarely observed in neighboring countries ̶ SIT2/H2 and SIT159/T1 Tuscany ̶ were conspicuously represented in Argentina. The infrequent Beijing patterns belonged to Peruvian patients. We conclude that the genotype diversity observed reflects the influence of the Hispanic colonization and more recent immigration waves from Mediterranean and neighboring countries. Unlike in Brazil, the RD^rio type does not play a major role in the tuberculosis epidemic in Buenos Aires.

Molecular Characterization of Mycobacterium tuberculosis Strains with TB-SPRINT

We evaluated Tuberculosis-Spoligo-Rifampicin-Isoniazid Typing (TB-SPRINT), a microbead-based method for spoligotyping and detection of rifampicin and isoniazid resistance in Mycobacterium tuberculosis. For that, 67 M. tuberculosis complex strains were retrospectively selected. Membrane-based spoligotyping, restriction fragment length polymorphism, DNA sequencing/pyrosequencing of rpoB, katG, and inhA promoter, TB-SPRINT, and SNP typing were performed. Concordance between spoligotyping methods was 99.6% (2,785/2,795 spoligotype data points). For most of the discordant cases, the same lineage was assigned with both methods. Concordance between phenotypic drug susceptibility testing and TB-SPRINT for detecting rifampicin and isoniazid resistance was 98.4% (63/64) and 93.8% (60/64), respectively. Concordance between DNA sequencing/pyrosequencing and TB-SPRINT for detecting mutations in rpoB, katG, and inhA were 98.4% (60/61), 100% (64/64), and 96.9% (62/64), respectively. In conclusion, TB-SPRINT is a rapid and easy-to-perform assay for genotyping and detecting drug resistance in a single tube; therefore, it may be a useful tool to improve epidemiological surveillance.

Mycobacterium bovis and Other Uncommon Members of the Mycobacterium tuberculosis Complex

Since its discovery by Theobald Smith, Mycobacterium bovis has been a human pathogen closely related to animal disease. At present, M. bovis tuberculosis is still a problem of importance in many countries and is considered the main cause of zoonotic tuberculosis throughout the world. Recent development of molecular epidemiological tools has helped us to improve our knowledge about transmission patterns of this organism, which causes a disease indistinguishable from that caused by Mycobacterium tuberculosis. Diagnosis and treatment of this mycobacterium are similar to those for conventional tuberculosis, with the important exceptions of constitutive resistance to pyrazinamide and the fact that multidrug-resistant and extremely drug-resistant M. bovis strains have been described. Among other members of this complex, Mycobacterium africanum is the cause of many cases of tuberculosis in West Africa and can be found in other areas mainly in association with immigration. M. bovis BCG is the currently available vaccine for tuberculosis, but it can cause disease in some patients. Other members of the M. tuberculosis complex are mainly animal pathogens with only exceptional cases of human disease, and there are even some strains, like "Mycobacterium canettii," which is a rare human pathogen that could have an important role in the knowledge of the evolution of tuberculosis in the history.

Molecular characterization of Mycobacterium tuberculosis isolates from Tehran, Iran by restriction fragment length polymorphism analysis and spoligotyping

INTRODUCTION Characterization of Mycobacterium tuberculosis (MTB) isolates by DNA fingerprinting has contributed to tuberculosis (TB) control. The aim of this study was to determine the genetic diversity of MTB isolates from Tehran province in Iran. METHODS MTB isolates from 60 Iranian and 10 Afghan TB patients were fingerprinted by standard IS6110-restriction fragment length polymorphism (RFLP) analysis and spoligotyping. RESULTS The copy number of IS6110 ranged from 10-24 per isolate. The isolates were classified into 22 clusters showing ≥ 80% similarity by RFLP analysis. Fourteen multidrug-resistant (MDR) isolates were grouped into 4 IS6110-RFLP clusters, with 10 isolates [71% (95% CI: 45-89%)] in 1 cluster, suggesting a possible epidemiological linkage. Eighteen Iranian isolates showed ≥ 80% similarity with Afghan isolates. There were no strains with identical fingerprints. Spoligotyping of 70 isolates produced 23 distinct patterns. Sixty (85.7%) isolates were grouped into 13 clusters, while the remaining 10 isolates (14.2%) were not clustered. Ural (formerly Haarlem4) (n = 22, 31.4%) was the most common family followed by Central Asian strain (CAS) (n = 18, 25.7%) and T (n = 9, 12.8%) families. Only 1strain was characterized as having the Beijing genotype. Among 60 Iranian and 10 Afghan MTB isolates, 25% (95% CI: 16-37) and 70% (95% CI: 39-89) were categorized as Ural lineage, respectively. CONCLUSIONS A higher prevalence of Ural family MTB isolates among Afghan patients than among Iranian patients suggests the possible transmission of this lineage following the immigration of Afghans to Iran.

Mycobacterium tuberculosis lineage 4 comprises globally distributed and geographically restricted sublineages

Generalist and specialist species differ in the breadth of their ecological niches. Little is known about the niche width of obligate human pathogens. Here we analyzed a global collection of Mycobacterium tuberculosis lineage 4 clinical isolates, the most geographically widespread cause of human tuberculosis. We show that lineage 4 comprises globally distributed and geographically restricted sublineages, suggesting a distinction between generalists and specialists. Population genomic analyses showed that, whereas the majority of human T cell epitopes were conserved in all sublineages, the proportion of variable epitopes was higher in generalists. Our data further support a European origin for the most common generalist sublineage. Hence, the global success of lineage 4 reflects distinct strategies adopted by different sublineages and the influence of human migration.

Methodological and Clinical Aspects of the Molecular Epidemiology of Mycobacterium tuberculosis and Other Mycobacteria

Molecular typing has revolutionized epidemiological studies of infectious diseases, including those of a mycobacterial etiology. With the advent of fingerprinting techniques, many traditional concepts regarding transmission, infectivity, or pathogenicity of mycobacterial bacilli have been revisited, and their conventional interpretations have been challenged. Since the mid-1990s, when the first typing methods were introduced, a plethora of other modalities have been proposed. So-called molecular epidemiology has become an essential subdiscipline of modern mycobacteriology. It serves as a resource for understanding the key issues in the epidemiology of tuberculosis and other mycobacterial diseases. Among these issues are disclosing sources of infection, quantifying recent transmission, identifying transmission links, discerning reinfection from relapse, tracking the geographic distribution and clonal expansion of specific strains, and exploring the genetic mechanisms underlying specific phenotypic traits, including virulence, organ tropism, transmissibility, or drug resistance. Since genotyping continues to unravel the biology of mycobacteria, it offers enormous promise in the fight against and prevention of the diseases caused by these pathogens. In this review, molecular typing methods for Mycobacterium tuberculosis and nontuberculous mycobacteria elaborated over the last 2 decades are summarized. The relevance of these methods to the epidemiological investigation, diagnosis, evolution, and control of mycobacterial diseases is discussed.

Turkish and Japanese Mycobacterium tuberculosis sublineages share a remote common ancestor

Two geographically distant M. tuberculosis sublineages, Tur from Turkey and T3-Osaka from Japan, exhibit partially identical genotypic signatures (identical 12-loci MIRU-VNTR profiles, distinct spoligotyping patterns). We investigated T3-Osaka and Tur sublineages characteristics and potential genetic relatedness, first using MIRU-VNTR locus analysis on 21 and 25 samples of each sublineage respectively, and second comparing Whole Genome Sequences of 8 new samples to public data from 45 samples uncovering human tuberculosis diversity. We then tried to date their Most Recent Common Ancestor (MRCA) using three calibrations of SNP accumulation rate (long-term=0.03SNP/genome/year, derived from a tuberculosis ancestor of around 70,000years old; intermediate=0.2SNP/genome/year derived from a Peruvian mummy; short-term=0.5SNP/genome/year). To disentangle between these scenarios, we confronted the corresponding divergence times with major human history events and knowledge on human genetic divergence. We identified relatively high intrasublineage diversity for both T3-Osaka and Tur. We definitively proved their monophyly; the corresponding super-sublineage (referred to as "T3-Osa-Tur") shares a common ancestor with T3-Ethiopia and Ural sublineages but is only remotely related to other Euro-American sublineages such as X, LAM, Haarlem and S. The evolutionary scenario based on long-term evolution rate being valid until T3-Osa-Tur MRCA was not supported by Japanese fossil data. The evolutionary scenario relying on short-term evolution rate since T3-Osa-Tur MRCA was contradicted by human history and potential traces of past epidemics. T3-Osaka and Tur sublineages were found likely to have diverged between 800y and 2000years ago, potentially at the time of Mongol Empire. Altogether, this study definitively proves a strong genetic link between Turkish and Japanese tuberculosis. It provides a first hypothesis for calibrating TB Euro-American lineage molecular clock; additional studies are needed to reliably date events corresponding to intermediate depths in tuberculosis phylogeny.

Combined Genotypic, Phylogenetic, and Epidemiologic Analyses of Mycobacterium tuberculosis Genetic Diversity in the Rhône Alpes Region, France

BACKGROUND

The present work relates to identification and a deep molecular characterization of circulating Mycobacterium tuberculosis complex (MTBC) strains in the Rhône-Alpes region, France from 2000 to 2010. It aimed to provide with a first snapshot of MTBC genetic diversity in conjunction with bacterial drug resistance, type of disease and available demographic and epidemiologic characteristics over an eleven-year period, in the south-east of France.

METHODS

Mycobacterium tuberculosis complex (MTBC) strains isolated in the Rhône-Alpes region, France (n = 2257, 1 isolate per patient) between 2000 and 2010 were analyzed by spoligotyping. MIRU-VNTR typing was applied on n = 1698 strains (with full results available for 974 strains). The data obtained were compared with the SITVIT2 database, followed by detailed genotyping, phylogenetic, and epidemiologic analyses in correlation with anonymized data on available demographic, and epidemiologic characteristics, and location of disease (pulmonary or extrapulmonary TB).

RESULTS

The most predominant spoligotyping clusters were SIT53/T1 (n = 346, 15.3%) > SIT50/H3 (n = 166, 7.35%) > SIT42/LAM9 (n = 125, 5.5%) > SIT1/Beijing (n = 72, 3.2%) > SIT47/H1 (n = 71, 3.1%). Evolutionary-recent strains belonging to the Principal Genetic Group (PGG) 2/3, or Euro-American lineages (T, LAM, Haarlem, X, S) were predominant and represented 1768 or 78.33% of all isolates. For strains having drug resistance information (n = 1119), any drug resistance accounted for 14.83% cases vs. 1.52% for multidrug resistance (MDR); and was significantly more associated with age group 21-40 years (p-value<0.001). Extra-pulmonary TB was more common among female patients while pulmonary TB predominated among men (p-value<0.001; OR = 2.16 95%CI [1.69; 2.77]). Also, BOV and CAS lineages were significantly well represented in patients affected by extra-pulmonary TB (p-value<0.001). The origin was known for 927/2257 patients: 376 (40.6%) being French-born vs. 551 (59.4%) Foreign-born. French patients were significantly older (mean age: 58.42 yrs 95%CI [56.04; 60.80]) than Foreign-born patients (mean age: 42.38 yrs. 95%CI [40.75; 44.0]).

CONCLUSION

The study underlined the importance of imported TB cases on the genetic diversity and epidemiologic characteristics of circulating MTBC strains in Rhône-Alpes region, France over a large time-period. It helps better understand intricate relationships between certain lineages and geographic origin of the patients, and pinpoints genotypic and phylogenetic specificities of prevailing MTBC strains. Lastly, it also demonstrated a slow decline in isolation of M. africanum lineage in this region between 2000 and 2010.

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

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

CRISPR-Based Typing and Next-Generation Tracking Technologies

Bacteria occur ubiquitously in nature and are broadly relevant throughout the food supply chain, with diverse and variable tolerance levels depending on their origin, biological role, and impact on the quality and safety of the product as well as on the health of the consumer. With increasing knowledge of and accessibility to the microbial composition of our environments, food supply, and host-associated microbiota, our understanding of and appreciation for the ratio of beneficial to undesirable bacteria are rapidly evolving. Therefore, there is a need for tools and technologies that allow definite, accurate, and high-resolution identification and typing of various groups of bacteria that include beneficial microbes such as starter cultures and probiotics, innocuous commensals, and undesirable pathogens and spoilage organisms. During the transition from the current molecular biology-based PFGE (pulsed-field gel electrophoresis) gold standard to the increasingly accessible omics-level whole-genome sequencing (WGS) N-gen standard, high-resolution technologies such as CRISPR-based genotyping constitute practical and powerful alternatives that provide valuable insights into genome microevolution and evolutionary trajectories. Indeed, several studies have shown potential for CRISPR-based typing of industrial starter cultures, health-promoting probiotic strains, animal commensal species, and problematic pathogens. Emerging CRISPR-based typing methods open new avenues for high-resolution typing of a broad range of bacteria and constitute a practical means for rapid tracking of a diversity of food-associated microbes.

Molecular Epidemiology of Mycobacterium bovis in Humans and Cattle

Bovine tuberculosis (bTB), caused by Mycobacterium bovis (M. bovis), is a serious re-emerging disease in both animals and humans. The evolution of the Multi- and Extensively drug-resistant M. bovis strains (MDR-TB and XDR-TB) represents a global threat to public health. Worldwide, the disease is responsible for great economic losses in the veterinary field, serious threat to the ecosystem, and about 3.1% of human TB cases, up to 16% in Tanzania. Only thorough investigation to understand the pathogen's epidemiology can help in controlling the disease and minimizing its threat. For this purpose, various tools have been developed for use in advanced molecular epidemiological studies of bTB, either alone or in combination with standard conventional epidemiological approaches. These techniques enable the analysis of the intra- and inter-species transmission dynamics of bTB. The delivered data can reveal detailed insights into the source of infection, correlations among human and bovine isolates, strain diversity and evolution, spread, geographical localization, host preference, tracing of certain virulence factors such as antibiotic resistance genes, and finally the risk factors for the maintenance and spread of M. bovis. They also allow for the determination of epidemic and endemic strains. This, in turn, has a significant diagnostic impact and helps in vaccine development for bTB eradication programs. The present review discusses many topics including the aetiology, epidemiology and importance of M. bovis, the prevalence of bTB in humans and animals in various countries, the molecular epidemiology of M. bovis, and finally applied molecular epidemiological techniques.

Genetic diversity and molecular epidemiology of multidrug-resistant Mycobacterium tuberculosis in Minas Gerais State, Brazil

Background

We aimed to characterize the genetic diversity of drug-resistant Mycobacterium tuberculosis (MTb) clinical isolates and investigate the molecular epidemiology of multidrug-resistant (MDR) tuberculosis from Minas Gerais State, Brazil.

Methods

One hundred and four MTb clinical isolates were assessed by IS6110-RFLP, 24-locus mycobacterial interspersed repetitive units variable-number tandem repeats (MIRU-VNTR), TB-SPRINT (simultaneous spoligotyping and rifampicin-isoniazid drug-resistance mutation analysis) and 3R-SNP-typing (analysis of single-nucleotide polymorphisms in the genes involved in replication, recombination and repair functions).

Results

Fifty-seven different IS6110-RFLP patterns were found, among which 50 had unique patterns and 17 were grouped into seven clusters. The discriminatory index (Hunter and Gaston, HGDI) for RFLP was 0.9937. Ninety-nine different MIRU-VNTR patterns were found, 95 of which had unique patterns and nine isolates were grouped into four clusters. The major allelic diversity index in the MIRU-VNTR loci ranged from 0.6568 to 0.7789. The global HGDI for MIRU-VNTR was 0.9991. Thirty-two different spoligotyping profiles were found: 16 unique patterns (n = 16) and 16 clustered profiles (n = 88). The HGDI for spoligotyping was 0.9009. The spoligotyped clinical isolates were phylogenetically classified into Latin-American Mediterranean (66.34 %), T (14.42 %), Haarlem (5.76 %), X (1.92 %), S (1.92 %) and U (unknown profile; 8.65 %). Among the U isolates, 77.8 % were classified further by 3R-SNP-typing as 44.5 % Haarlem and 33.3 % LAM, while the 22.2 % remaining were not classified. Among the 104 clinical isolates, 86 were identified by TB-SPRINT as MDR, 12 were resistant to rifampicin only, one was resistant to isoniazid only, three were susceptible to both drugs, and two were not successfully amplified by PCR. A total of 42, 28 and eight isolates had mutations in rpoB positions 531, 526 and 516, respectively. Correlating the cluster analysis with the patient data did not suggest recent transmission of MDR-TB.

Conclusions

Although our results do not suggest strong transmission of MDR-TB in Minas Gerais (using a classical 100 % MDR-TB identical isolates cluster definition), use of a smoother cluster definition (>85 % similarity) does not allow us to fully eliminate this possibility; hence, around 20–30 % of the isolates we analyzed might be MDR-TB transmission cases.

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.

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

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

High-Throughput CRISPR Typing of Mycobacterium tuberculosis Complex and Salmonella enterica Serotype Typhimurium

Spoligotyping was developed almost 18 years ago and still remains a popular first-lane genotyping technique to identify and subtype Mycobacterium tuberculosis complex (MTC) clinical isolates at a phylogeographic level. For other pathogens, such as Salmonella enterica, recent studies suggest that specifically designed spoligotyping techniques could be interesting for public health purposes. Spoligotyping was in its original format a reverse line-blot hybridization method using capture probes designed on "spacers" and attached to a membrane's surface and a PCR product obtained from clustered regularly interspaced short palindromic repeats (CRISPRs). Cowan et al. and Fabre et al. were the first to propose a high-throughput Spoligotyping method based on microbeads for MTC and S. enterica serotype Typhimurium, respectively. The main advantages of the high-throughput Spoligotyping techniques we describe here are their low cost, their robustness, and the existence (at least for MTC) of very large databases that allow comparisons between spoligotypes from anywhere.

A robust SNP barcode for typing Mycobacterium tuberculosis complex strains

Strain-specific genomic diversity in the Mycobacterium tuberculosis complex (MTBC) is an important factor in pathogenesis that may affect virulence, transmissibility, host response and emergence of drug resistance. Several systems have been proposed to classify MTBC strains into distinct lineages and families. Here, we investigate single-nucleotide polymorphisms (SNPs) as robust (stable) markers of genetic variation for phylogenetic analysis. We identify ~92k SNP across a global collection of 1,601 genomes. The SNP-based phylogeny is consistent with the gold-standard regions of difference (RD) classification system. Of the ~7k strain-specific SNPs identified, 62 markers are proposed to discriminate known circulating strains. This SNP-based barcode is the first to cover all main lineages, and classifies a greater number of sublineages than current alternatives. It may be used to classify clinical isolates to evaluate tools to control the disease, including therapeutics and vaccines whose effectiveness may vary by strain type.

Genetic variation in Mycobacterium tuberculosis complex (MTBC) bacteria is responsible for differences in factors such as virulence and transmissibility. Here, the authors analyse the genomes of 1,601 MTBC isolates from diverse geographic locations and identify 62 SNPs that may be used to resolve lineages and sublineages of these strains.

Optimizing multiplex SNP-based data analysis for genotyping of Mycobacterium tuberculosis isolates

Background

Multiplex ligation-dependent probe amplification (MLPA) is a powerful tool to identify genomic polymorphisms. We have previously developed a single nucleotide polymorphism (SNP) and large sequence polymorphisms (LSP)-based MLPA assay using a read out on a liquid bead array to screen for 47 genetic markers in the Mycobacterium tuberculosis genome. In our assay we obtain information regarding the Mycobacterium tuberculosis lineage and drug resistance simultaneously. Previously we called the presence or absence of a genotypic marker based on a threshold signal level. Here we present a more elaborate data analysis method to standardize and streamline the interpretation of data generated by MLPA. The new data analysis method also identifies intermediate signals in addition to classification of signals as positive and negative. Intermediate calls can be informative with respect to identifying the simultaneous presence of sensitive and resistant alleles or infection with multiple different Mycobacterium tuberculosis strains.

Results

To validate our analysis method 100 DNA isolates of Mycobacterium tuberculosis extracted from cultured patient material collected at the National TB Reference Laboratory of the National Center for Tuberculosis and Lung Diseases in Tbilisi, Republic of Georgia were tested by MLPA. The data generated were interpreted blindly and then compared to results obtained by reference methods. MLPA profiles containing intermediate calls are flagged for expert review whereas the majority of profiles, not containing intermediate calls, were called automatically. No intermediate signals were identified in 74/100 isolates and in the remaining 26 isolates at least one genetic marker produced an intermediate signal.

Conclusion

Based on excellent agreement with the reference methods we conclude that the new data analysis method performed well. The streamlined data processing and standardized data interpretation allows the comparison of the Mycobacterium tuberculosis MLPA results between different experiments. All together this will facilitate the implementation of the MLPA assay in different settings.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-572) contains supplementary material, which is available to authorized users.

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.

Insights into the population structure of Mycobacterium tuberculosis using spoligotyping and RDRio in a southeastern Brazilian prison unit

Tuberculosis (TB) is still a serious public health problem, continuing to be an important threat for confined populations. We used spoligotyping to estimate the genotypic clades of Mycobacterium tuberculosis isolates from inmates in two blocks in a southeastern Brazilian prison unit, with TB incidence rate of 8185/100.000. The Latin American Mediterranean (LAM) clade is well represented in the country, and the LAM specific molecular markers, RD(Rio) large sequence polymorphism and the SNP on the Rv3062 [ligB(1212)], were used to characterize spoligotype signatures from prison isolates. Typing of RD(Rio) and ligB increase LAM clade from 66.7% (n=72/108) to 69.4% (n=75). The LAM2 SIT17 (n=23) and SIT179 (n=12) signatures comprised one third of all isolates, followed by Haarlem (11.5%, n=12), T (8.7%, n=9) and X (5.7%, n=6) clades. Strains with unknown signatures represented 5.5% (n=6), and four (3.7%) did not match any lineage. We observed RD(Rio) among 64 (59.2%) isolates, and 54 (50%) were of the LAM clade. In particular, the LAM2/RD(Rio) sub-lineage was significantly associated with clustering (p=0.02) and its frequency was higher (32%) when compared to that of the previous general TB cases in RJ (4.29%). Overall cluster frequency defined by spoligotyping/IS6110-RFLP was 62%. The two evolutionary markers helped to evaluate some LAM signature misconceptions and demonstrate that LAM2/RD(Rio) was found with high frequency, hitherto being unnoticed. All these data, allied to high clustering, imply that public health measures to minimize the escalation of TB in prison is essential, and both spoligotyping as well as RD(Rio) would be useful tools to monitor the effects of the measures with respect to M. tuberculosis lineage variation.

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.

Bacterial strain typing

Over the course of the past several decades, rapid advancements in molecular technologies have revolutionized the practice of public health microbiology, and have fundamentally changed the nature, accuracy, and timeliness of laboratory data for outbreak investigation and response. Whole-genome sequencing, in particular, is becoming an increasingly feasible and cost-effective approach for near real-time high-resolution strain typing, genomic characterization, and comparative analyses. This review discusses the current state of the art in bacterial strain typing for outbreak investigation and infectious disease surveillance, and the impact of emerging genomic technologies on the field of public health microbiology.

The role of CRISPR-Cas systems in virulence of pathogenic bacteria

Clustered regularly interspaced short palindromic repeats (CRISPR) and CRISPR-associated (Cas) genes are present in many bacterial and archaeal genomes. Since the discovery of the typical CRISPR loci in the 1980s, well before their physiological role was revealed, their variable sequences have been used as a complementary typing tool in diagnostic, epidemiologic, and evolutionary analyses of prokaryotic strains. The discovery that CRISPR spacers are often identical to sequence fragments of mobile genetic elements was a major breakthrough that eventually led to the elucidation of CRISPR-Cas as an adaptive immunity system. Key elements of this unique prokaryotic defense system are small CRISPR RNAs that guide nucleases to complementary target nucleic acids of invading viruses and plasmids, generally followed by the degradation of the invader. In addition, several recent studies have pointed at direct links of CRISPR-Cas to regulation of a range of stress-related phenomena. An interesting example concerns a pathogenic bacterium that possesses a CRISPR-associated ribonucleoprotein complex that may play a dual role in defense and/or virulence. In this review, we describe recently reported cases of potential involvement of CRISPR-Cas systems in bacterial stress responses in general and bacterial virulence in particular.

Mycobacterium tuberculosis Latin American-Mediterranean family and its sublineages in the light of robust evolutionary markers

Mycobacterium tuberculosis has a clonal population structure, and the Latin American-Mediterranean (LAM) family is one of the largest and most widespread within this species, showing evidence for remarkable pathobiology and a confusing phylogeny. Here, we applied robust phylogenetic markers to study the evolution of the LAM family and its major sublineages circulating in Russia and neighboring countries. A total of 250 M. tuberculosis isolates were confirmed to belong to the LAM family based on the analysis of the LAM-specific single-nucleotide polymorphisms (SNPs) in the Rv3062 and Rv0129c genes. At this stage, the family status was rectified for 121 isolates misleadingly assigned by CRISPR spoligotyping to non-LAM families (T1- or T5-RUS1). Consequently, the reestimated LAM prevalence rate increased 2-fold in Russia and Kazakhstan and 4-fold in Belarus. The majority (91.8 to 98.7%) of the LAM isolates from all three countries belonged to the LAM-RUS sublineage. In contrast, the Ibero-American LAM RD-Rio sublineage was identified in only 7 Russian isolates. Taken together, our findings and further analyses suggest a monophyletic origin of LAM-RUS: at a historically distant time, in Russia, in a small founding bacterial/human population. Its dissemination pattern and high prevalence rate in Northern Eurasia may indicate a long-term coexistence of the LAM-RUS sublineage and local human populations hypothetically leading to coadaptation and reduced pathogenicity of the relatively more ancient clones, such as spoligotype international type 254 (SIT254), compared to the more recent SIT252 and SIT266 clones. In contrast, rare LAM RD-Rio isolates were likely brought to Russia through occasional human contact. The spread of RD-Rio strains is not as global as commonly claimed and is determined largely by human migration flows (rather than by pathobiological properties of these strains). Consequently, a host population factor appears to play a major role in shaping the in situ dissemination pattern of the imported strains in an autochthonous population.

Genotyping of ancient Mycobacterium tuberculosis strains reveals historic genetic diversity

The evolutionary history of the Mycobacterium tuberculosis complex (MTBC) has previously been studied by analysis of sequence diversity in extant strains, but not addressed by direct examination of strain genotypes in archaeological remains. Here, we use ancient DNA sequencing to type 11 single nucleotide polymorphisms and two large sequence polymorphisms in the MTBC strains present in 10 archaeological samples from skeletons from Britain and Europe dating to the second–nineteenth centuries AD. The results enable us to assign the strains to groupings and lineages recognized in the extant MTBC. We show that at least during the eighteenth–nineteenth centuries AD, strains of M. tuberculosis belonging to different genetic groups were present in Britain at the same time, possibly even at a single location, and we present evidence for a mixed infection in at least one individual. Our study shows that ancient DNA typing applied to multiple samples can provide sufficiently detailed information to contribute to both archaeological and evolutionary knowledge of the history of tuberculosis.

Overview and phylogeny of Mycobacterium tuberculosis complex organisms: implications for diagnostics and legislation of bovine tuberculosis

Members of the Mycobacterium tuberculosis complex (MTBC) cause a serious disease with similar pathology, tuberculosis; in this review, bovine tuberculosis will be considered as disease caused by any member of the MTBC in bovids. Bovine tuberculosis is responsible for significant economic loss due to costly eradication programs and trade limitations and poses a threat to both endangered and protected species as well as to public health. We here give an overview on all members of the MTBC, focusing on their isolation from different animal hosts. We also review the recent advances made in elucidating the evolutionary and phylogenetic relationships of members of the MTBC. Because the nomenclature of the MTBC is controversial, its members have been considered species, subspecies or ecotypes, this review discusses the possible implications for diagnostics and the legal consequences of naming of new species.