Analysis of rate of change of IS6110 RFLP patterns of Mycobacterium tuberculosis based on serial patient isolates

Characterization of a trinucleotide repeat sequence (CGG)5 and potential use in restriction fragment length polymorphism typing of Mycobacterium tuberculosis

The genomes of 28 bacterial strains, including mycobacterial species Mycobacterium tuberculosis and Mycobacterium bovis, were analyzed for the presence of a special class of microsatellite, that of trinucleotide repeat sequences (TRS). Results of a search of all 10 possible TRS motifs (i.e., CCT, CGG, CTG, GAA, GAT, GTA, GTC, GTG, GTT, and TAT) with five or more repeating units showed that (CGG)(5) was highly represented within the genomic DNA of M. tuberculosis and M. bovis. Most of the (CGG)(5) repeats in the genome were within the open reading frames of two large gene families encoding PE_PGRS and PPE proteins that have the motifs Pro-Glu (PE) and Pro-Pro-Glu (PPE). (CGG)(5)-probed Southern hybridization showed that some mycobacterial species, such as Mycobacterium marinum, Mycobacterium kansasii, and Mycobacterium szulgai, possess many copies of (CGG)(5) in their genomes. Analysis of clinical isolates obtained from Tokyo and Warsaw with both IS6110 and (CGG)(5) probes showed that there is an association between the fingerprinting patterns and the geographic origin of the isolates and that (CGG)(5) fingerprinting patterns were relatively more stable than IS6110 patterns. The (CGG)(5) repeat is a unique sequence for some mycobacterial species, and (CGG)(5) fingerprinting can be used as an epidemiologic method for these species as well as IS6110 fingerprinting can. If these two fingerprinting methods are used together, the precise analysis of M. tuberculosis isolates will be accomplished. (CGG)(5)-based fingerprinting is particularly useful for M. tuberculosis isolates with few or no insertion elements and for the identification of other mycobacterial species when informative probes are lacking.

Evidence for positive selection on Mycobacterium tuberculosis within patients

BACKGROUND

While the pathogenesis and epidemiology of tuberculosis are well studied, relatively little is known about the evolution of the infectious agent Mycobacterium tuberculosis, especially at the within-host level. The insertion sequence IS6110 is a genetic marker that is widely used to track the transmission of tuberculosis between individuals. This and other markers may also facilitate our understanding of the disease within patients.

RESULTS

This article presents three lines of evidence supporting the action of positive selection on M. tuberculosis within patients. The arguments are based on a comparison between empirical findings from molecular epidemiology, and population genetic models of evolution. Under the hypothesis of neutrality of genotypes, 1) the mutation rate of the marker IS6110 is unusually high, 2) the time it takes for substitutions to occur within patients is too short, and 3) the amount of polymorphism within patients is too low.

CONCLUSIONS

Empirical observations are explained by the action of positive selection during infection, or alternatively by very low effective population sizes. I discuss the possible roles of antibiotic treatment, the host immune system and extrapulmonary dissemination in creating opportunities for positive selection.

Molecular fingerprinting of Mycobacterium bovis subsp. caprae isolates from central Europe

To study the dissemination of Mycobacterium bovis subsp. caprae, 79 European isolates from cattle, humans, and other hosts were examined by spoligotyping and IS6110 restriction fragment length polymorphism (RFLP) analysis. Among a total of 11 different spoligotypes identified, type C1 proved to be predominant (n = 62). Five of the spoligotypes are described for the first time. A total of 43 different RFLP types were identified, thus allowing further differentiation for epidemiological tracking. Isolates from a series of outbreaks in one village proved to be of the same spoligotype and of identical or closely related RFLP types.

Definition of the Beijing/W lineage of Mycobacterium tuberculosis on the basis of genetic markers

Mycobacterium tuberculosis Beijing genotype strains are highly prevalent in Asian countries and in the territory of the former Soviet Union. They are increasingly reported in other areas of the world and are frequently associated with tuberculosis outbreaks and drug resistance. Beijing genotype strains, including W strains, have been characterized by their highly similar multicopy IS6110 restriction fragment length polymorphism (RFLP) patterns, deletion of spacers 1 to 34 in the direct repeat region (Beijing spoligotype), and insertion of IS6110 in the genomic dnaA-dnaN locus. In this study the suitability and comparability of these three genetic markers to identify members of the Beijing lineage were evaluated. In a well-characterized collection of 1,020 M. tuberculosis isolates representative of the IS6110 RFLP genotypes found in The Netherlands, strains of two clades had spoligotypes characteristic of the Beijing lineage. A set of 19 Beijing reference RFLP patterns was selected to retrieve all Beijing strains from the Dutch database. These reference patterns gave a sensitivity of 98.1% and a specificity of 99.7% for identifying Beijing strains (defined by spoligotyping) in an international database of 1,084 strains. The usefulness of the reference patterns was also assessed with large DNA fingerprint databases in two other European countries and for identification strains from the W lineage found in the United States. A standardized definition for the identification of M. tuberculosis strains belonging to the Beijing/W lineage, as described in this work, will facilitate further studies on the spread and characterization of this widespread genotype family of M. tuberculosis strains.

Molecular epidemiology of tuberculosis and recent developments in understanding the epidemiology of tuberculosis

The development of molecular typing techniques during the last decade has led to remarkable improvements in epidemiological studies of tuberculosis (TB) and other infectious disease. In this paper, the burden of TB disease, current issues about disease control in developed countries including Australia and the contribution of molecular typing to improved understanding of TB epidemiology are reviewed.

Transposition rates of Mycobacterium tuberculosis IS6110 restriction fragment length polymorphism patterns

To determine the rate at which IS6110 restriction fragment length polymorphism (RFLP) patterns in Mycobacterium tuberculosis change over time, we applied a smooth nonparametric survival model to several data sets, including data from previous publications on the rate of change. The results strongly suggest a simple parametric model, with an instantaneous change at time zero and essentially a zero rate of change thereafter. Our interpretation of the results is that at the time of collection of the first isolate, more than one strain is present. We speculate that the selection of mutant strains is most likely during rapid growth, revival of the dormant bacteria, and/or adaptation to a new host. The parameter most accurately describing changing RFLP patterns is the proportion of isolates with band changes, rather than the half-life or the rate of change.

Inhibitor-associated transposition events in Corynebacterium glutamicum

In up to 100% of all bacteria grown in the presence of initially inhibitory concentrations of five diverse inhibitors, an extra copy of the resident insertion element IS 31831 was found in specific chromosomal regions, the sites of which apparently depended on the inhibitor used. Thus, in nine out of nine independently isolated cyanide-associated transpositions, the acquired copy was located within an ORF encoding a protein related to the hypothetical but conserved protein YeiH of Escherichia coli. A putative Sox box upstream of the yeiH gene implicates superoxide as a potential regulator of the gene, a possibility further supported by the finding that superoxide dismutase (SodA) is overexpressed in cells cultured in cyanide-containing medium. Neither the cyanide-associated nor any of the other transposition mutations appeared to confer any discernible phenotypic advantage upon cells grown in the presence or absence of the inhibitors, as revealed most stringently by mixed-cell experiments. An alternative, albeit heterodox, explanation for the emergence of the mutants postulates a very high rate of transpositional activity in the presence of inhibitors. The initial emergence of the mutants was found to depend crucially upon the cell density. Thus, when growth medium was supplemented with 50 mM fluoropyruvate and inoculated to a density of 2 x 10(7) cfu/ml, single colonies with heterogeneous restriction fragment length polymorphisms (RFLPs) were routinely isolated at a frequency of 6 to 16% after 1-2 days of incubation. After 3 days, 10-36% of the colonies showed RFLPs, but the type was now dominated by the fluoropyruvate-specific RFLP, which, at higher resolution, invariably proved to be heterogeneous. This heterogeneity proved that these specific mutants were of multiple origin, indicating that clonal enrichment was irrelevant to their emergence. It is suggested that the presence of the inhibitor induces the development of hyper-transpositional activity, which is regulated by a soluble bacterial product.

Stability of polymorphic GC-rich repeat sequence-containing regions of Mycobacterium tuberculosis

Mycobacterium tuberculosis cultures were subjected to DNA fingerprinting with IS6110- and polymorphic GC-rich sequence (PGRS)-containing probes. The PGRS banding patterns remained highly stable during multiple cultures of specimens from one disease episode (0.5% changed) and during transmission in patients with close contact (1.9% changed). Characteristic PGRS-restriction fragment length polymorphism motifs for different strain groupings may indicate distant evolutionary events leading to the differentiation of M. tuberculosis strain lineages.

Stable association between strains of Mycobacterium tuberculosis and their human host populations

Mycobacterium tuberculosis is an important human pathogen in virtually every part of the world. Here we investigate whether distinct strains of M. tuberculosis infect different human populations and whether associations between host and pathogen populations are stable despite global traffic and the convergence of diverse strains of the pathogen in cosmopolitan urban centers. The recent global movement and transmission history of 100 M. tuberculosis isolates was inferred from a molecular epidemiologic study of tuberculosis that spans 12 years. Genetic relationships among these isolates were deduced from the distribution of large genomic deletions, which were identified by DNA microarray and confirmed by PCR and sequence analysis. Phylogenetic analysis of these deletions indicates that they are unique event polymorphisms and that horizontal gene transfer is extremely rare in M. tuberculosis. In conjunction with the epidemiological data, phylogenies reveal three large phylogeographic regions. A host's region of origin is predictive of the strain of tuberculosis he or she carries, and this association remains strong even when transmission takes place in a cosmopolitan urban center outside of the region of origin. Approximate dating of the time since divergence of East Asian and Philippine clades of M. tuberculosis suggests that these lineages diverged centuries ago. Thus, associations between host and pathogen populations appear to be highly stable.

Analysis of genetic polymorphisms affecting the four phospholipase C (plc) genes in Mycobacterium tuberculosis complex clinical isolates

The Mycobacterium tuberculosis genome contains four highly related genes which present significant similarity to Pseudomonas aeruginosa genes encoding phospholipase C enzymes. Three of these genes, plcA, plcB and plcC, are organized in tandem (locus plcABC). The fourth gene, plcD, is located in a different region. This study investigates variations in plcABC and plcD genes in clinical isolates of M. tuberculosis, Mycobacterium africanum and ‘Mycobacterium canettii’. Genetic polymorphisms were examined by PCR, Southern blot hybridization, sequence analysis and RT-PCR. Seven M. tuberculosis isolates contain insertions of IS6110 elements within plcA, plcC or plcD. In 19 of 25 M. tuberculosis isolates examined, genomic deletions were identified, resulting in loss of parts of genes or complete genes from the plcABC and/or plcD loci. Partial plcD deletion was observed in one M. africanum isolate. In each case, deletions were associated with the presence of a copy of the IS6110 element and in all occurrences IS6110 was transposed in the same orientation. A mechanism of deletion resulting from homologous recombination of two copies of IS6110 was recognized in a group of genetically related M. tuberculosis isolates. Five M. tuberculosis isolates presented major polymorphisms in the plcABC and plcD regions, along with loss of expression competence that affected all four plc genes. Phospholipase C is a well-known bacterial virulence factor. The precise role of phospholipase C in the pathogenicity of M. tuberculosis is unknown, but considering the potential importance that the plc genes may have in the virulence of the tubercle bacillus, the study of isolates cultured from patients with active tuberculosis bearing genetic variations affecting these genes may provide insights into the significance of phospholipase C enzymes for tuberculosis pathogenicity.

Use of genetic distance as a measure of ongoing transmission of Mycobacterium tuberculosis

The stability of the genotypic marker IS6110, used to define the epidemiology of Mycobacterium tuberculosis, is one of the most important factors influencing the interpretation of DNA fingerprint data. We propose that evolved strains should be considered together with clustered strains to represent chains of ongoing transmission. For the present study we used a large set of fingerprint data for strains collected between 1992 and 1998 from residents of a community with a high incidence of tuberculosis in Cape Town, South Africa. Interstrain genetic distances were calculated by counting the banding pattern mismatches in the IS6110 DNA fingerprints of different isolates. These data demonstrate that the propensity to change by one or two bands is independent of the IS6110 copy number. Hence, the genetic distance between pairs of isolates can be simply expressed as the number of differences in the banding patterns. From this foundation, a data set which identifies newly evolved strains has been generated. Inclusion of these evolved strains into various molecular epidemiological calculations significantly increased the estimate of ongoing transmission in this study setting. The indication is that nearly all cases of tuberculosis in this community are due to ongoing transmission. This has important implications for tuberculosis control, as it indicates that the control measures used at present are unable to reduce the level of transmission. This technique may also be applicable to the study of low-incidence tuberculosis outbreaks as well as the analysis of epidemiological data from other disease epidemics.

The transmission of tuberculosis in the light of new molecular biological approaches

This review briefly summarises the recent achievements in tuberculosis epidemiology associated with the introduction of molecular methods, and considers the implications of these methods for the understanding of occupational tuberculosis transmission. Special attention is paid to the relative contribution of recently transmitted tuberculosis; risk factors for recent transmission; and the occurrence and frequency of exogenous reinfection. There is a need for occupational epidemiological studies, which should combine the methods of "classical" epidemiology with those of molecular epidemiology.

Molecular differentiation of Mycobacterium bovis isolates. Review of main techniques and applications

Until recently, none of the Mycobacterium bovis typing techniques permitted a satisfactory differentiation of isolates. During the last 10 years, the genome of pathogenic mycobacteria has been extensively studied, and phylogenetic analyses have shown that all (except Mycobacterium avium) belong to a single genetic species: the Mycobacterium tuberculosis complex. This increase in knowledge about the genome of these bacteria has lead to the discovery of molecular markers that allow us to differentiate isolates. Because of the phylogenetic proximity of the strains, even if most of these markers have been discovered in M. tuberculosis, they could be successfully adapted to the other bacteria of the M. tuberculosis complex, especially M. bovis. The most common markers in use today are the IS6110 insertion sequence, the direct repeat (DR) region, the poly(GC) rich (PGRS) sequences and the variable number tandem repeats (VNTR) sequences. The corresponding typing techniques are briefly described, and current knowledge of polymorphism and marker stability is detailed. If molecular markers are to offer wide perspectives for field studies, these two characteristics (polymorphism and stability) must be taken into account when choosing the marker(s) used in a study. In this context, examples of the application of molecular typing techniques for M. bovis are reviewed, on the one hand with epidemiological studies for which the major problem is the comparison between isolates and, on the other, with more general studies about the population genetics of M. bovis in a given country, and about its history and its phylogeny.

Recent transmission of tuberculosis in Madrid: application of capture–recapture analysis to conventional and molecular epidemiology

BACKGROUND

Population-based studies using a combination of molecular techniques and conventional epidemiological methods have been used to study the dynamics of tuberculosis (TB) transmission but the relative utility of each technique has not yet been established.

METHODS

A prospective population-based molecular and epidemiological study of patients diagnosed with TB was conducted in three urban districts of Madrid (Spain) during 1997-1999. Analysis was performed using the capture-recapture method including covariates in which conventional epidemiological data and the information on clustered cases obtained by DNA fingerprinting were regarded as independent and complementary procedures.

RESULTS

The estimate obtained by molecular analysis alone, that 31.6% of TB cases were due to recent transmission, was revised to 44.8% (95% CI: 31.4-58.2) using the capture-recapture method. The estimated completeness of the combined databases for identification of recent transmission was 59.2%. Underestimation of the true prevalence of recent transmission was higher with conventional epidemiology than molecular analysis, particularly for patients <35 years old and those with a history of imprisonment.

CONCLUSIONS

In this study, use of the capture-recapture technique allowed us to combine epidemiological information obtained by conventional and molecular methods to quantify the number of cases of recently transmitted TB in the community and identify specific populations at high risk of disease. This information is clearly important because such groups are a prime target for improved TB control measures. In the long term, this combination of techniques may contribute significantly to control the spread of TB.

Widespread pyrazinamide-resistant Mycobacterium tuberculosis family in a low-incidence setting

An unusually high prevalence of pyrazinamide (PZA) monoresistance in Mycobacterium tuberculosis has been observed in Quebec. In the absence of a recognized outbreak, we hypothesized that these isolates most likely represented reactivation of an old endemic strain in this low-incidence area. A case-control study of 77 PZA-resistant isolates with a specific Quebec mutation and 253 PZA-susceptible control M. tuberculosis isolates was undertaken. By molecular analysis, all 77 case isolates shared a unique mutation profile in the pncA gene which was not present in control isolates. While control isolates manifested diverse IS6110 restriction fragment length polymorphism (RFLP) patterns, spoligotypes, and major genetic groups, case isolates had similar but nonidentical IS6110 RFLP patterns, had common spoligotypes, and were confined to one major genetic group, suggesting a common clonal ancestor. By epidemiologic and geographic analyses, however, there were no significant differences between the cases and the controls. We conclude that a clonally related family of PZA-monoresistant M. tuberculosis isolates in Quebec represents historic rather than recent transmission.

Locating transposable element polymorphisms in bacterial genomes

Although whole-genome sequencing is greatly extending our knowledge of the genetic capacity of those bacterial species, it is only directly informative for the particular strain sequenced. Many bacterial species exhibit more or less genetic polymorphism within their populations and characterising this variety is an extremely important way of elucidating the biology of these species. Often genomic polymorphisms are associated with multicopy elements, particularly transposable elements. We describe a novel method that efficiently characterises the sequences of such polymorphisms. We have optimised heminested inverse PCR (hINVPCR) to assess the diversity of insertional polymorphisms of a transposable element (IS6110) in clinical isolates of Mycobacterium tuberculosis. To increase the yield of information, genomic DNA was digested with different endonucleases (Bsp1286I, HaeII or PvuI), and primers based on both the 5' and 3' ends of IS6110 were used to amplify and determine the genomic sequence upstream (or downstream) of the transposable element. We found that both the choice of restriction enzyme and the use of primers at both ends of the transposable element significantly increased the diversity of the insertion sites identified. Band stabbing was incorporated into the method as an alternative to cloning in order to screen large number of isolates at a sequence level in a rapid and labour-efficient fashion. We describe some of the purposes to which such data can be put.

Estimating change rates of genetic markers using serial samples: applications to the transposon IS6110 in Mycobacterium tuberculosis

In infectious disease epidemiology, it is useful to know how quickly genetic markers of pathogenic agents evolve while inside hosts. We propose a modular framework with which these genotype change rates can be estimated. The estimation scheme requires a model of the underlying process of genetic change, a detection scheme that filters this process into observable quantities, and a monitoring scheme that describes the timing of observations. We study a linear "birth-shift-death" model for change in transposable element genotypes, obtaining maximum-likelihood estimators for various detection and monitoring schemes. The method is applied to serial genotypes of the transposon IS6110 in Mycobacterium tuberculosis. The estimated birth rate of 0.0161 (events per copy of the transposon per year) and death rate of 0.0108 are both significantly larger than the estimated shift rate of 0.0018. The sum of these estimates, which corresponds to a "half-life" of 2.4 years for a typical strain that has 10 copies of the element, substantially exceeds a previous estimate of 0.0135 total changes per copy per year. We consider experimental design issues that enable the precision of estimates to be improved. We also discuss extensions to other markers and implications for molecular epidemiology.

High-throughput epidemiologic typing in clinical microbiology

Mapping, and ultimately preventing, the dissemination of infectious agents is an important topic in public health. Newly developed molecular-microbiological methods have contributed significantly to recent advances in the efficient tracking of the nosocomial and environmental spread of microbial pathogens. Not only has the application of novel technologies led to improved understanding of microbial epidemiology, but the concepts of population structure and dynamics of many of the medically significant microorganisms have advanced significantly also. Currently, genetic identification of microbes is also within the reach of clinical microbiology laboratory professionals including those without specialized technology research interests. This review summarizes the possibilities for high-throughput molecular-microbiological typing in adequately equipped medical microbiology laboratories from both clinical and fundamental research perspectives. First, the development and application of methods for large-scale comparative typing of serially isolated microbial strains are discussed. The outcome of studies employing these methods allows for long-term epidemiologic surveillance of infectious diseases. Second, recent methods enable an almost nucleotide-by-nucleotide genetic comparison of smaller numbers of strains, thereby facilitating the identification of the genetic basis of, for instance, medically relevant microbiological traits. Whereas the first approach provides insights into the dynamic spread of infectious agents, the second provides insights into intragenomic dynamics and genetic functionality. The current state of technology is summarized, and future perspectives are sketched.

Annual Mycobacterium tuberculosis infection risk and interpretation of clustering statistics

Several recent studies have used proportions of tuberculosis cases sharing identical DNA fingerprint patterns (i.e., isolate clustering) to estimate the extent of disease attributable to recent transmission. Using a model of introduction and transmission of strains with different DNA fingerprint patterns, we show that the properties and interpretation of clustering statistics may differ substantially between settings. For some unindustrialized countries, where the annual risk for infection has changed little over time, 70% to 80% of all age groups may be clustered during a 3-year period, which underestimates the proportion of disease attributable to recent transmission. In contrast, for a typical industrialized setting (the Netherlands), clustering declines with increasing age (from 75% to 15% among young and old patients, respectively) and underestimates the extent of recent transmission only for young patients. We conclude that, in some settings, clustering is an unreliable indicator of the extent of recent transmission.

Stability of variable-number tandem repeats of mycobacterial interspersed repetitive units from 12 loci in serial isolates of Mycobacterium tuberculosis

Variable number tandem repeats (VNTRs) of elements named mycobacterial interspersed repetitive units (MIRUs) have previously been identified in 12 minisatellite loci of the Mycobacterium tuberculosis genome. These markers allow reliable high-throughput genotyping of M. tuberculosis and represent a portable approach to global molecular epidemiology of M. tuberculosis. To assess their temporal stability, we genotyped 123 serial isolates, separated by up to 6 years and belonging to a variety of distinct IS6110 restriction fragment length polymorphism (RFLP) families, from 56 patients who had positive sputum cultures. All 12 MIRU VNTR loci were completely identical within the groups of serial isolates in 55 out of 56 groups (98.2%), although 11 pairs of isolates from the same patients with conserved MIRU VNTRs displayed slightly different IS6110 RFLP profiles. In a single case, serial isolates with an unchanged IS6110 RFLP profile showed a change in 1 out of 12 MIRU VNTR loci. These results indicate that MIRU VNTRs are stable over time and therefore are suitable for reliable follow-up of patients chronically infected with tuberculosis over long periods. Moreover, they support MIRU VNTR genotyping as a powerful first-line method followed by subtyping by IS6110 RFLP to define ongoing transmission clusters.

Impact of genotyping of Mycobacterium tuberculosis on public health practice in Massachusetts

Massachusetts was one of seven sentinel surveillance sites in the National Tuberculosis Genotyping and Surveillance Network. From 1996 through 2000, isolates from new patients with tuberculosis (TB) underwent genotyping. We describe the impact that genotyping had on public health practice in Massachusetts and some limitations of the technique. Through genotyping, we explored the dynamics of TB outbreaks, investigated laboratory cross-contamination, and identified Mycobacterium tuberculosis strains, transmission sites, and accurate epidemiologic links. Genotyping should be used with epidemiologic follow-up to identify how resources can best be allocated to investigate genotypic findings.

Molecular epidemiology of tuberculosis in a sentinel surveillance population

We conducted a population-based study to assess demographic and risk-factor correlates for the most frequently occurring Mycobacterium tuberculosis genotypes from tuberculosis (TB) patients. The study included all incident, culture-positive TB patients from seven sentinel surveillance sites in the United States from 1996 to 2000. M. tuberculosis isolates were genotyped by IS6110-based restriction fragment length polymorphism and spoligotyping. Genotyping was available for 90% of 11923 TB patients. Overall, 48% of cases had isolates that matched those from another patient, including 64% of U.S.-born and 35% of foreign-born patients. By logistic regression analysis, risk factors for clustering of genotypes were being male, U.S.-born, black, homeless, and infected with HIV; having pulmonary disease with cavitations on chest radiograph and a sputum smear with acid-fast bacilli; and excessive drug or alcohol use. Molecular characterization of TB isolates permitted risk correlates for clusters and specific genotypes to be described and provided information regarding cluster dynamics over time.

Mycobacterium microti infection (vole tuberculosis) in wild rodent populations

Mycobacterium microti (vole tuberculosis) infections in small wild mammals were first described more than 60 years ago in several populations in Great Britain. Few studies of vole tuberculosis have been undertaken since then, and little is known about the relationship between M. microti isolates originating from different populations or at different times or of the prevalence of this infection in wild rodent populations, despite human cases of M. microti infections being increasingly reported. In this study, field voles (Microtus agrestis), bank voles (Clethrionomys glareolus), and wood mice (Apodemus sylvaticus) were found to be infected, with up to 8% having external tuberculous signs, in wild populations in Northumberland and Cheshire, England. Spoligotyping applied directly to the clinical material simultaneously detected and typed M. microti bacteria in skin lesions, lymph glands, and internal abcesses. IS6110 restriction fragment length polymorphism typing of cultured bacteria was used to compare these isolates with previously isolated strains from both animals and humans. This demonstrated that although the current rodent isolates were distinct from those isolated from voles in the 1930s in Great Britain, they had a high degree of similarity to these strains and were distinct from the M. microti isolates from humans, a pig, and a ferret from The Netherlands. Thus, M. microti infection seems to be widespread in wild rodent populations, but more studies are needed to understand how M. microti might be transmitted from animals to humans and to determine better the zoonotic risk posed.

Mycobacterium africanum subtype II is associated with two distinct genotypes and is a major cause of human tuberculosis in Kampala, Uganda

The population structure of 234 Mycobacterium tuberculosis complex strains obtained during 1995 and 1997 from tuberculosis patients living in Kampala, Uganda (East Africa), was analyzed by routine laboratory procedures, spoligotyping, and IS6110 restriction fragment length polymorphism (RFLP) typing. According to biochemical test results, 157 isolates (67%) were classified as M. africanum subtype II (resistant to thiophen-2-carboxylic acid hydrazide), 76 isolates (32%) were classified as M. tuberculosis, and 1 isolate was classified as classical M. bovis. Spoligotyping did not lead to clear differentiation of M. tuberculosis and M. africanum, but all M. africanum subtype II isolates lacked spacers 33 to 36, differentiating them from M. africanum subtype I. Moreover, spoligotyping was not sufficient for differentiation of isolates on the strain level, since 193 (82%) were grouped into clusters. In contrast, in the IS6110-based dendrogram, M. africanum strains were clustered into two closely related strain families (Uganda I and II) and clearly separated from the M. tuberculosis isolates. A further characteristic of both M. africanum subtype II families was the absence of spoligotype spacer 40. All strains of family I also lacked spacer 43. The clustering rate obtained by the combination of spoligotyping and RFLP IS6110 analysis was similar for M. africanum and M. tuberculosis, as 46% and 49% of the respective isolates were grouped into clusters. The results presented demonstrate that M. africanum subtype II isolates from Kampala, Uganda, belong to two closely related genotypes, which may represent unique phylogenetic branches within the M. tuberculosis complex. We conclude that M. africanum subtype II is the main cause of human tuberculosis in Kampala, Uganda.

Use of molecular techniques to distinguish between treatment failure and exogenous reinfection with Mycobacterium tuberculosis

We investigated the means by which drug resistance emerges among drug-susceptible Mycobacterium tuberculosis strains during antituberculosis therapy. Patients who experienced failure of treatment for active pulmonary tuberculosis, who initially received diagnoses of infection with drug-susceptible M. tuberculosis, and who had had at least 3 isolates tested for drug susceptibility were selected from a 6-year period in the Estonian National Reference Laboratory archive. Eleven patients from whom 35 sequential isolates of M. tuberculosis had been obtained were recruited into the study. Their clinical data and treatment charts were analyzed and correlated with drug-susceptibility patterns and IS6110 restriction fragment-length polymorphism (RFLP) profiles. Six patients excreted isogenic drug-susceptible M. tuberculosis strains, whereas, in the other 5 patients, the isolated strain shifted from a susceptible to a resistant phenotype. In all cases, this shift correlated to a shift in RFLP pattern, which showed reinfection with a new strain. Exogenous reinfection with drug-resistant M. tuberculosis may be misinterpreted as the emergence of drug resistance if molecular testing techniques are not used.

Calculation of the stability of the IS6110 banding pattern in patients with persistent Mycobacterium tuberculosis disease

The interpretation of molecular epidemiologic data of Mycobacterium tuberculosis infection is dependent on the understanding of the stability and evolutionary characteristics of the DNA fingerprinting marker used to classify clinical isolates. This study investigated the stability of the IS6110 banding pattern in serial tuberculosis isolates collected from patients resident in an area with a high incidence of tuberculosis. Evolutionary changes were observed in 4% of the strains, and a half-life (t(1/2)) of 8.74 years was calculated, assuming a constant rate of change over time. This rate may be composed of a high rate of change seen during the early disease phase (t(1/2) = 0.57 years) and a low rate of change seen in the late disease phase (t(1/2) = 10.69 years). The early rate probably reflects change occurring during active growth prior to therapy, while the low late rate may reflect change occurring during or after treatment. We demonstrate that the calculation of these rates is strongly influenced by the time interval between onset of disease and sputum sampling. These calculations are further complicated by partial replacement of the original strain population, resulting in the sporadic appearance of clonal variants in sputum specimens. Therefore, the true extent of genetic diversity may be underestimated within each host, thereby influencing molecular epidemiological data used to establish transmission chains.

Evolution of the IS6110-based restriction fragment length polymorphism pattern during the transmission of Mycobacterium tuberculosis

Interpretation of the molecular epidemiological data of Mycobacterium tuberculosis is dependent on the validity of the assumptions that have been made. It is assumed that the IS6110 banding pattern is sufficiently stable to define epidemiological events representing ongoing transmission. However, molecular epidemiological data also support the observation that the IS6110 banding pattern may change over time. Factors affecting this rate may include the nature and duration of disease in a host and the opportunity to experience different host environments during the transmission cycle. To estimate the rate of IS6110 change occurring during the process of transmission, M. tuberculosis isolates from epidemiologically linked patients were genotypically characterized by restriction fragment length polymorphism (RFLP) analysis. The identification of IS6110 banding pattern changes during ongoing transmission suggested that a rate could be estimated. IS6110 change was significantly associated with strains with >5 IS6110 elements (P = 0.013) and was not observed in low-copy-number isolates. The minimum rate of appearance of variant strains was calculated to be 0.14 variant cases per source-case per year. This data suggest that clustering of isolates based on identical RFLP patterns is expected to underestimate transmission in patients infected with high-copy-number isolates. A model based on the rate of appearance of both variant and invariant strains demonstrates that the genotypically defined population structure may change by 18.6% during the study period of approximately 6.5 years. The implications for the use of RFLP data for epidemiologic study are discussed.

Molecular evolution of Mycobacterium tuberculosis: phylogenetic reconstruction of clonal expansion

SETTING

M. tuberculosis isolates were collected from patients attending health clinics in a high incidence urban community and in a low incidence rural setting in South Africa.

OBJECTIVE

To reconstruct the evolutionary history of a group of closely related M. tuberculosis isolates using IS6110, DRr and MTB484(1) restriction fragment length polymorphism (RFLP) data.

DESIGN

Mycobacterium tuberculosis isolates containing an average of ten IS6110 elements, with a similarity index of > or = 65% were genotypically classified by DNA fingerprinting using the IS6110 derived probes IS-3' and IS-5', as well as the DRr and MTB484(1) probes, in combination with PvuII or Hinfl endonuclease digestion. These RFLP data were subjected to phylogenetic analysis using both genetic distance and parsimony algorithms.

RESULTS

Phylogenetic analysis predicted the existence of two independently evolving lineages, possibly evolving from a common ancestral strain. The topology of the phylogenetic tree was supported by comprehensive bootstrapping and the specific partitioning of DNA methylation phenotypes. The observed difference in the branch lengths of the two lineages may suggest differential evolutionary rates. Isolates collected from different geographical regions demonstrate independent evolution, suggesting that it is highly unlikely that strains have been recently transmitted between the two regions. The number of evolutionary events identified in this strain family differs significantly from that of previously characterized strain families, implying that evolutionary rate may be strain family dependent.

CONCLUSION

Based on this analysis we propose that the algorithm used to calculate recent epidemiological events should be revised to incorporate the evolutionary characteristics of individual strain families, thereby enhancing the accuracy of molecular epidemiological calculations.

Spoligotype diversity of Mycobacterium bovis strains isolated in France from 1979 to 2000

The molecular fingerprints of 1,349 isolates of Mycobacterium bovis received between 1979 and August 2000 at Agence Française de Sécurité Sanitaire des Aliments (Afssa) have been obtained by spoligotyping. The majority of the isolates (1,266) were obtained from cattle living in France. An apparently high level of heterogeneity was observed between isolates. One hundred sixty-one spoligotypes were observed in total, of which 153 were from French isolates. The two predominant spoligotypes, designated BCG-like and GB54, accounted for 26 and 12% of the isolates, respectively. In addition, 84% of the spoligotypes were found fewer than 10 times. Analysis of the results by clustering and parsimony-based algorithms revealed that the majority of the spoligotypes were closely related. The predominant spoligotype was identical to that of the vaccine strain Mycobacterium bovis BCG, which was isolated in France at the end of the 19th century. Some spoligotypes were closely associated with restricted geographical areas. Interestingly, some spoligotypes, which were frequently observed in France, were also observed in neighboring countries. Conversely, few spoligotypes were common to France and England, and those that were shared were observed at very different frequencies. This last point illustrates the potential role for an international data bank, which could help trace the spread of M. bovis across national borders.

Epidemiological evidence of the spread of a Mycobacterium tuberculosis strain of the Beijing genotype on Gran Canaria Island

Molecular epidemiological studies suggest that particular Mycobacterium tuberculosis strains have an enhanced capacity to spread within a community. One strain, the Beijing genotype, has been associated with outbreaks in a number of communities throughout the world. IS6110 restriction fragment length polymorphism (RFLP) analysis was performed on M. tuberculosis isolates from 566 of the 721 patients (78.5%) diagnosed with tuberculosis (TB) on Gran Canaria Island from 1993 to 1996, as well as 35% of isolates from 1991-1992 (85 strains). RFLP identification of the family of strains of the Beijing genotype was confirmed by spoligotyping. Medical records of all patients were reviewed and epidemiological links were identified. Of 566 M. tuberculosis isolates from 1993 to 1996 with RFLP available, 72% belonged to clusters. The largest contained 75 cases and was caused by a strain of the Beijing genotype that was introduced to the island in 1993. It was found in 10 patients in 1993 (5.5%), 12 in 1994 (8.1%), 18 in 1995 (16.4%), and 35 in 1996 (27.1%). Epidemiological linkage was confirmed for 68% of cases. This study has demonstrated rapid dissemination of this strain of the Beijing genotype. This genotype might play an important role in the future of the worldwide tuberculosis epidemic.

Demonstration of reinfection and reactivation in HIV-negative autopsied cases of secondary tuberculosis: multilesional genotyping of Mycobacterium tuberculosis utilizing IS 6110 and other repetitive element-based DNA fingerprinting

SETTING

Secondary tuberculosis may follow reinfection or endogenous reactivation. The design of effective preventive and treatment protocols requires knowledge about the relative importance of these mechanisms in specific communities. Molecular typing of Mycobacterium tuberculosis has permitted linkage of cases and demonstration of patterns of inter- or intrapatient strain diversity correlating with reinfection and reactivation.

OBJECTIVE

The use of DNA fingerprinting to examine intrapatient strain diversity in autopsied, HIV-negative individuals resident in a high incidence community.

DESIGN

Autopsy (12 cases) and pneumonectomy (one case) permitted multilesional sampling for bacterial culture and comparative genotyping. Mycobacterium tuberculosis strains were typed using probes directed against IS 6110, the direct repeat sequence (DR) and MTB484(1).

RESULTS

In two patients, the demonstration of pulmonary infection by two distinct strains suggested dual infection of these individuals. In one other case, the strain isolated from a Ghon-focus was identical to that obtained from all secondary cavitating lesions, which suggested reactivation of the primary infection. In the remaining cases, all isolates were identical, but primary lesions could not be identified with certainty or were culture negative. One of these cases showed evidence of strain evolution.

CONCLUSIONS

Multilesional strain genotyping suggested both reinfection and reactivation in a series of HIV-negative autopsied patients with secondary tuberculosis. In most (11/13) cases, the same strain (or clonal variant) was present in all lesions, pulmonary and extrapulmonary.

Optimal estimation of transposition rates of insertion sequences for molecular epidemiology

Outbreaks of infectious disease can be confirmed by identifying clusters of DNA fingerprints among bacterial isolates from infected individuals. This procedure makes assumptions about the underlying properties of the genetic marker used for fingerprinting. In particular, it requires that each fingerprint changes sufficiently slowly within an individual that isolates from separate individuals infected by the same strain will exhibit similar or identical fingerprints. We propose a model for the probability that an individual's fingerprint will change over a given period of time. We use this model together with published data in order to estimate the fingerprint change rate for IS6110 in human tuberculosis, obtaining a value of 0.0139 changes per copy per year. Although we focus on insertion sequences (IS), our method applies to other fingerprinting techniques such as pulsed-field gel electrophoresis (PFGE). We suggest sampling intervals that produce the least error in estimates of the fingerprint change rate, as well as sample sizes that achieve specified levels of error in the estimate.

Molecular epidemiology of tuberculosis: recent developments and applications

The standard method for the typing of Mycobacterium tuberculosis is still IS6110 restriction fragment length polymorphism (RFLP). This method has been widely used and has provided information on the variety and distribution of tuberculosis strain types across the globe. Recently, IS6110 RFLP has been used to investigate the question of reinfection versus reactivation, examine the existence of multiple infection, and track the spread of multidrug-resistant tuberculosis. There have also been efforts to increase our understanding of the biologic characteristics of IS6110. These studies have resulted in a clearer understanding of fingerprinting data and increased our understanding of the evolution and pathogenicity of this organism.

Comparison of variable number tandem repeat and IS6110-restriction fragment length polymorphism analyses for discrimination of high- and low-copy-number IS6110 Mycobacterium tuberculosis isolates

The present study was designed to evaluate the use of variable number tandem repeat (VNTR) and IS6110-restriction fragment length polymorphism (RFLP) analyses in combination as a two-step strategy for discrimination (as measured by the Hunter-Gaston Discrimination Index [HGDI]) of both high- and low-copy-number IS6110 Mycobacterium tuberculosis isolates compared to IS6110-RFLP alone with an unselected collection of isolates. Individually, IS6110-RFLP fingerprinting produced six clusters that accounted for 69% of the low-copy-number IS6110 isolates (five clusters) and 5% of the high-copy-number IS6110 isolates (one cluster). A total of 39% of all the isolates were clustered (HGDI = 0.97). VNTR analysis generated a total of 35 different VNTR allele profile sets from 93 isolates (HGDI = 0.938). Combining IS6110-RFLP analysis with VNTR analysis reduced the overall percentage of clustered isolates to 29% (HGDI = 0.988) and discriminated a further 27% of low-copy-number isolates that would have been clustered by IS6110-RFLP alone. The use of VNTR analysis as an initial typing strategy facilitates further analysis by IS6110-RFLP, and more importantly, VNTR analysis subdivides some IS6110-RFLP-defined clusters containing low- and single-copy IS6110 isolates.

Molecular epidemiology of methicillin-resistant Staphylococcus aureus strains: state of affairs and tomorrow' s possibilities

Methicillin-resistant strains of Staphylococcus aureus (MRSA) have posed a clinical threat for nearly 40 years. During these years, an array of additional technologies suited for identification of MRSA below the species level has become available. The technologies, whether they assess phenotype or genotype, provide data that can be used for elucidation of the routes of dissemination of individual MRSA types. This review summarizes the current state of affairs with respect to the quality of the various laboratory techniques and includes descriptions of novel strategies such as binary typing and multilocus sequence typing (MLST). Drawbacks of procedures will be compared, and the value of molecular typing in the elucidation of complex biological phenomena, such as epidemicity, carriage, and reduced vancomycin susceptibility, will be indicated. Means for integrated assessment of bacterial biology, epidemiology, and population structure will be discussed.

High-resolution minisatellite-based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology

The worldwide threat of tuberculosis to human health emphasizes the need to develop novel approaches to a global epidemiological surveillance. The current standard for Mycobacterium tuberculosis typing based on IS6110 restriction fragment length polymorphism (RFLP) suffers from the difficulty of comparing data between independent laboratories. Here, we propose a high-resolution typing method based on variable number tandem repeats (VNTRs) of genetic elements named mycobacterial interspersed repetitive units (MIRUs) in 12 human minisatellite-like regions of the M. tuberculosis genome. MIRU-VNTR profiles of 72 different M. tuberculosis isolates were established by PCR analysis of all 12 loci. From 2 to 8 MIRU-VNTR alleles were identified in the 12 regions in these strains, which corresponds to a potential of over 16 million different combinations, yielding a resolution power close to that of IS6110-RFLP. All epidemiologically related isolates tested were perfectly clustered by MIRU-VNTR typing, indicating that the stability of these MIRU-VNTRs is adequate to track outbreak episodes. The correlation between genetic relationships inferred from MIRU-VNTR and IS6110-RFLP typing was highly significant. Compared with IS6110-RFLP, high-resolution MIRU-VNTR typing has the considerable advantages of being fast, appropriate for all M. tuberculosis isolates, including strains that have a few IS6110 copies, and permitting easy and rapid comparison of results from independent laboratories. This typing method opens the way to the construction of digital global databases for molecular epidemiology studies of M. tuberculosis.

Molecular epidemiology of tuberculosis and other mycobacterial infections: main methodologies and achievements

In the last decade, DNA fingerprint techniques have become available to study the interperson transmission of tuberculosis and other mycobacterial infections. These methods have facilitated epidemiological studies at a population level. In addition, the species identification of rarely encountered mycobacteria has improved significantly. This article describes the state of the art of the main molecular typing methods for Mycobacterium tuberculosis complex and non-M. tuberculosis complex (atypical) mycobacteria. Important new insights that have been gained through molecular techniques into epidemiological aspects and diagnosis of mycobacterial diseases are highlighted.

Genetic heterogeneity in Mycobacterium tuberculosis isolates reflected in IS6110 restriction fragment length polymorphism patterns as low-intensity bands

Mycobacterium tuberculosis isolates with identical IS6110 restriction fragment length polymorphism (RFLP) patterns are considered to originate from the same ancestral strain and thus to reflect ongoing transmission. In this study, we investigated 1,277 IS6110 RFLP patterns for the presence of multiple low-intensity bands (LIBs), which may indicate infections with multiple M. tuberculosis strains. We did not find any multiple LIBs, suggesting that multiple infections are rare in the Netherlands. However, we did observe a few LIBs in 94 patterns (7.4%) and examined the nature of this phenomenon. With single-colony cultures it was found that LIBs mostly represent mixed bacterial populations with slightly different RFLP patterns. Mixtures were expressed in RFLP patterns as LIBs when 10 to 30% of the DNA analyzed originated from a bacterial population with another RFLP pattern. Presumably, a part of the LIBs did not represent mixed bacterial populations, as in some clusters all strains exhibited LIBs in their RFLP patterns. The occurrence of LIBs was associated with increased age in patients. This may reflect either a gradual change of the bacterial population in the human body over time or IS6110-mediated genetic adaptation of M. tuberculosis to changes in the environmental conditions during the dormant state or reactivation thereafter.

Tuberculosis treatment failure and drug resistance--same strain or reinfection?

Tuberculosis patients may have Mycobacterium tuberculosis in their sputum at the end of treatment, and may show new drug resistance, due to either inadequate treatment of the original episode or reinfection with a new strain during therapy. In a cohort study of mineworkers with tuberculosis in South Africa, 57 of 438 patients had positive sputum cultures 6 months after recruitment in 1995. Of the 31 patients who initially had fully sensitive strains, 3 developed multidrug resistance (MDR) and 3 single-drug resistance (SDR). Of the 6 who started with SDR, 3 became MDR. HIV infection was not associated with drug resistance at enrollment or 6 months later. We compared pairs of DNA fingerprints from isolates of M. tuberculosis at recruitment and 6 months later in the 48 patients for whom we had both available. In 45, the pairs were identical. In 1 patient, although both isolates were fully sensitive, the later fingerprint had 1 less band (transposition). In 2 pairs, the fingerprint patterns were completely different: one seemed to be the result of laboratory error and the other was a true reinfection with an MDR strain. Despite a high risk of infection, with a moderate proportion of background drug-resistant strains (11% SDR, 6% MDR), reinfection is not a common cause of treatment failure or drug resistance at 6 months.

Molecular surveillance of European quinolone-resistant clinical isolates of Pseudomonas aeruginosa and Acinetobacter spp. using automated ribotyping

Nosocomial isolates of Pseudomonas aeruginosa and Acinetobacter spp. exhibit high rates of resistance to antibiotics and are often multidrug resistant. In a previous study (D. Milatovic, A. Fluit, S. Brisse, J. Verhoef, and F. J. Schmitz, Antimicrob. Agents Chemother. 44:1102-1107, 2000), isolates of these species that were resistant to sitafloxacin, a new advanced-generation fluoroquinolone with a high potency and a broad spectrum of antimicrobial activity, were found in high proportion in 23 European hospitals. Here, we investigate the clonal diversity of the 155 P. aeruginosa and 145 Acinetobacter spp. sitafloxacin-resistant isolates from that study by automated ribotyping. Numerous ribogroups (sets of isolates with indistinguishable ribotypes) were found among isolates of P. aeruginosa (n = 34) and Acinetobacter spp. (n = 16), but the majority of the isolates belonged to a limited number of major ribogroups. Sitafloxacin-resistant isolates (MICs > 2 mg/liter, used as a provisional breakpoint) showed increased concomitant resistance to piperacillin, piperacillin-tazobactam, ceftriaxone, ceftazidime, amikacin, gentamicin, and imipenem. The major ribogroups were repeatedly found in isolates from several European hospitals; these isolates showed higher levels of resistance to gentamicin and imipenem, and some of them appeared to correspond to previously described multidrug-resistant international clones of P. aeruginosa (serotype O:12) and Acinetobacter baumannii (clones I and II). Automated ribotyping, when used in combination with more discriminatory typing methods, may be a convenient library typing system for monitoring future epidemiological dynamics of geographically widespread multidrug-resistant bacterial clones.

Stability of IS6110 restriction fragment length polymorphism patterns of Mycobacterium tuberculosis strains in actual chains of transmission

The stability of IS6110 restriction fragment length polymorphism (RFLP) patterns of Mycobacterium tuberculosis strains in actual transmission chains has been assessed by analyzing the variability of IS6110 RFLP patterns of strains in fingerprint clusters that have been confirmed by classical epidemiological data. Forty susceptible and 35 drug-resistant (including 17 multidrug-resistant) M. tuberculosis strains obtained from 75 patients living in Germany have been analyzed. The epidemiological relationship among strains within the fingerprint clusters has been verified by family contacts (14 clusters) or by contact tracing of the public health offices (7 clusters). The time spans between the first and the last isolate of one cluster ranged from less than 1 to 29 months. Of the 75 strains only 1 showed a one-band variation when compared to the other nine isolates grouped in the same cluster, corresponding with a rate of change of approximately 1.9% per possible transmission (one index patient per cluster was subtracted from the total number of isolates). These results confirm a high degree of stability of IS6110 RFLP patterns of transmitted M. tuberculosis strains. Furthermore, the data presented indicate that isolates with identical IS6110 DNA fingerprint patterns are a good indicator for the rate of recent transmission in a study population.

Detection of a previously unamplified spacer within the DR locus of Mycobacterium tuberculosis: epidemiological implications

Spoligotyping, a method based on the variability of distribution of the 43 inter-direct repeat (DR) spacers of Mycobacterium tuberculosis and Mycobacterium bovis BCG, is useful to study the molecular epidemiology of bovine and human tuberculosis. Recently, a major family of M. tuberculosis clinical isolates named the Haarlem family, which did not contain spacers 31 and 33 to 36, was reported in a multicenter study. Independently, a data bank containing all the published spoligotypes showed that the two most prevalent spoligotypes in the world differed only by the presence or absence of spacer 31. A careful analysis of the DR locus sequence led us to hypothesize that spacer 31 may not have been amplified in some isolates with the primer sets DRa and DRb currently used for spoligotyping. Consequently, a modified spoligotyping method based on different combinations of the 36-bp DR and IS6110 primers was devised that was able to discriminate between the left and the right parts of the DR locus and demonstrated the presence of the previously unamplified spacer 31 for some of the clinical isolates. By analogy, we suggest that a single-spacer difference in some epidemiologically linked cases of tuberculosis may simply arise due to the insertion of an extra copy of IS6110 within the DR locus, leading to its asymmetrical disruption and subsequent lack of the DRa or DRb targets. The influence of the IS6110 preferential insertion sites within the DR locus on spoligotyping results should be further investigated.