Igor Mokrousov

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

BackgroundThe 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.ResultsThe 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.ConclusionOur 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.

SITVITWEB--a publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology.

Among various genotyping methods to study Mycobacterium tuberculosis complex (MTC) genotypic polymorphism, spoligotyping and mycobacterial interspersed repetitive units-variable number of DNA tandem repeats (MIRU-VNTRs) have recently gained international approval as robust, fast, and reproducible typing methods generating data in a portable format. Spoligotyping constituted the backbone of a publicly available database SpolDB4 released in 2006; nonetheless this method possesses a low discriminatory power when used alone and should be ideally used in conjunction with a second typing method such as MIRU-VNTRs for high-resolution epidemiological studies. We hereby describe a publicly available international database named SITVITWEB which incorporates such multimarker data allowing to have a global vision of MTC genetic diversity worldwide based on 62,582 clinical isolates corresponding to 153 countries of patient origin (105 countries of isolation). We report a total of 7105 spoligotype patterns (corresponding to 58,180 clinical isolates) - grouped into 2740 shared-types or spoligotype international types (SIT) containing 53,816 clinical isolates and 4364 orphan patterns. Interestingly, only 7% of the MTC isolates worldwide were orphans whereas more than half of SITed isolates (n=27,059) were restricted to only 24 most prevalent SITs. The database also contains a total of 2379 MIRU patterns (from 8161 clinical isolates) from 87 countries of patient origin (35 countries of isolation); these were grouped in 847 shared-types or MIRU international types (MIT) containing 6626 isolates and 1533 orphan patterns. Lastly, data on 5-locus exact tandem repeats (ETRs) were available on 4626 isolates from 59 countries of patient origin (22 countries of isolation); a total of 458 different VNTR patterns were observed - split into 245 shared-types or VNTR International Types (VIT) containing 4413 isolates) and 213 orphan patterns. Datamining of SITVITWEB further allowed to update rules defining MTC genotypic lineages as well to have a new insight into MTC population structure and worldwide distribution at country, sub-regional and continental levels. At evolutionary level, the data compiled may be useful to distinguish the occasional convergent evolution of genotypes versus specific evolution of sublineages essentially influenced by adaptation to the host. This database is publicly available at: http://www.pasteur-guadeloupe.fr:8081/SITVIT_ONLINE.

Snapshot of Moving and Expanding Clones of Mycobacterium tuberculosis and Their Global Distribution Assessed by Spoligotyping in an International Study

ABSTRACT 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.

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 11,708 patterns from as many clinical isolates originating from more than 90 countries. The 11,708 spoligotypes were clustered into 813 shared types. A total of 1,300 orphan patterns (clinical isolates showing a unique spoligotype) were also detected.

Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage

Mycobacterium tuberculosis strains of the Beijing lineage are globally distributed and are associated with the massive spread of multidrug-resistant (MDR) tuberculosis in Eurasia. Here we reconstructed the biogeographical structure and evolutionary history of this lineage by genetic analysis of 4,987 isolates from 99 countries and whole-genome sequencing of 110 representative isolates. We show that this lineage initially originated in the Far East, from where it radiated worldwide in several waves. We detected successive increases in population size for this pathogen over the last 200 years, practically coinciding with the Industrial Revolution, the First World War and HIV epidemics. Two MDR clones of this lineage started to spread throughout central Asia and Russia concomitantly with the collapse of the public health system in the former Soviet Union. Mutations identified in genes putatively under positive selection and associated with virulence might have favored the expansion of the most successful branches of the lineage.

High Prevalence of KatG Ser315Thr Substitution among Isoniazid-Resistant Mycobacterium tuberculosis Clinical Isolates from Northwestern Russia, 1996 to 2001

ABSTRACT A total of 204 isoniazid (INH)-resistant strains of Mycobacterium tuberculosis isolated from different patients in the northwestern region of Russia from 1996 to 2001 were screened by a PCR-restriction fragment length polymorphism (RFLP) assay. This assay uses HapII cleavage of an amplified fragment of the katG gene to detect the transversion 315AGC→ACC (Ser→Thr), which is associated with INH resistance. This analysis revealed a 93.6% prevalence of the katG S315T mutation in strains from patients with both newly and previously diagnosed cases of tuberculosis (TB). This mutation was not found in any of 57 INH-susceptible isolates included in the study. The specificity of the assay was 100%; all isolates that contained the S315T mutation were classified as resistant by a culture-based susceptibility testing method. The Beijing genotype, defined by IS6110-RFLP analysis and the spacer oligonucleotide typing (spoligotyping) method, was found in 60.3% of the INH-resistant strains studied. The katG S315T shift was more prevalent among Beijing genotype strains than among non-Beijing genotype strains: 97.8 versus 84.6%, respectively, for all isolates, including those from patients with new and previously diagnosed cases, isolated from 1999 to 2001 and 100.0 versus 86.5%, respectively, for isolates from patients with new cases isolated from 1996 to 2001. The design of this PCR-RFLP assay allows the rapid and unambiguous identification of the katG 315ACC mutant allele. The simplicity of the assay permits its implementation into routine practice in clinical microbiology laboratories in regions with a high incidence of TB where this mutation is predominant, including northwestern Russia.

Detection of embB306 Mutations in Ethambutol-Susceptible Clinical Isolates of Mycobacterium tuberculosis from Northwestern Russia: Implications for Genotypic Resistance Testing

ABSTRACT A total of 183 epidemiologically unlinked Mycobacterium tuberculosis isolates collected in the St. Petersburg area of Russia from 1996 to 2001 were screened for alterations in codon 306 of the embB gene; mutations in this codon are reported to confer resistance to ethambutol (EMB). The embB306 mutations were detected in 14 (48.3%) of 29 EMB-resistant strains and, quite surprisingly, in 48 (31.2%) of 154 EMB-susceptible strains. A discrepancy between the results of phenotypic and genotypic EMB resistance tests was restricted to the strains already resistant to other antitubercular (anti-TB) drugs. In particular, 40 (60%) of the 69 EMB-susceptible strains resistant to rifampin, isoniazid, and streptomycin but none of the 43 pansusceptible strains harbored an embB306 mutation. We hypothesize that the phenomenon observed could reflect the presence of a target other than EmbB for the drug in tubercle bacilli; this unknown target could be sensitized and affected, sensu lato, by EMB during treatment with other first-line anti-TB drugs. Comparison with DNA fingerprinting data showed that, irrespectively of the phenotypic susceptibility profiles, 46 (50.6%) of 91 Beijing family strains and 16 (17.4%) of 92 strains of other genotypes had a mutation in embB306.

Detection of Isoniazid-Resistant Mycobacterium tuberculosis Strains by a Multiplex Allele-Specific PCR Assay Targeting katG Codon 315 Variation

ABSTRACT We describe a simple multiplex allele-specific (MAS)-PCR assay to detect mutations in the second base of the katG gene codon 315, including AGC→ACC and ACA (Ser→Thr) substitutions that confer resistance to isoniazid (INH) in Mycobacterium tuberculosis clinical isolates. The 315 ACC allele is found in the majority of Inhr strains worldwide, especially in areas with a high incidence of tuberculosis. The 315 ACA allele is characteristic of the New York City multidrug-resistant (MDR) strain W and its progenies in the United States. The mutations in katG315 are revealed depending on the presence or absence of an indicative fragment amplified from the wild-type allele of this codon. Initially optimized on the purified DNA samples, the assay was then tested on crude cell lysates and auramine-stained sputum slide preparations with the same reproducibility and interpretability of profiles generated by agarose gel electrophoresis. The MAS-PCR assay can be used for the detection of resistance to INH in clinical laboratories in regions with a high prevalence of MDR M. tuberculosis strains.

Allele-Specific rpoB PCR Assays for Detection of Rifampin-Resistant Mycobacterium tuberculosis in Sputum Smears

ABSTRACT We describe an allele-specific PCR assay to detect mutations in three codons of the rpoB gene (516, 526, and 531) in Mycobacterium tuberculosis strains; mutations in these codons are reported to account for majority of M. tuberculosis clinical isolates resistant to rifampin (RIF), a marker of multidrug-resistant tuberculosis (MDR-TB). Three different allele-specific PCRs are carried out either directly with purified DNA (single-step multiplex allele-specific PCR), or with preamplified rpoB fragment (nested allele-specific PCR [NAS-PCR]). The method was optimized and validated following analysis of 36 strains with known rpoB sequence. A retrospective analysis of the 287 DNA preparations from epidemiologically unlinked RIF-resistant clinical strains from Russia, collected from 1996 to 2002, revealed that 247 (86.1%) of them harbored a mutation in one of the targeted rpoB codons. A prospective study of microscopy-positive consecutive sputum samples from new and chronic TB patients validated the method for direct analysis of DNA extracted from sputum smears. The potential of the NAS-PCR to control for false-negative results due to lack of amplification was proven especially useful in the study of these samples. The developed rpoB-PCR assay can be used in clinical laboratories to detect RIF-resistant and hence MDR M. tuberculosis in the regions with high burdens of the MDR-TB.

Molecular Characterization of Ofloxacin-Resistant Mycobacterium tuberculosis Strains from Russia

ABSTRACT In this work, we studied the variation in the gyrA and gyrB genes in ofloxacin- and multidrug-resistant Mycobacterium tuberculosis strains circulating in northwest Russia. Comparison with spoligotyping data suggested that similar to the spread of multidrug-resistant tuberculosis, the spread of fluoroquinolone-resistant tuberculosis in Russia may be due, at least partly, to the prevalence of the Beijing genotype in a local population of M. tuberculosis.