Previous Article | Next Article ![]()
Journal of Clinical Microbiology, May 2002, p. 1705-1708, Vol. 40, No. 5
0095-1137/02/$04.00+0 DOI: 10.1128/JCM.40.5.1705-1708.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
MRC Centre for Molecular and Cellular Biology, Department of Medical Biochemistry,1 Department of Pediatrics and Child Health, Faculty of Health Sciences, University of Stellenbosch, Tygerberg 7505, South Africa,2 Royal Netherlands Tuberculosis Association, 2510 CC The Hague, The Netherlands,3 Department of Medicine, McGill University Health Centre, Montreal H3G 1A4, Canada4
Received 3 December 2001/ Returned for modification 4 January 2002/ Accepted 30 January 2002
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The molecular epidemiology of Mycobacterium tuberculosis is most frequently based on the restriction fragment length polymorphism (RFLP) banding patterns generated by Southern hybridization with the probe IS6110 (10, 12). In community settings isolates with identical IS6110 genotypes are grouped together into clusters and are thought to represent recent epidemiological events (5), while isolates with IS6110 banding patterns unrelated to any others in the database (isolates with unique banding patterns) are thought to reflect reactivation of latent infection (1, 9). This methodology has been used in numerous settings to quantify the relative proportion of recent epidemiological events and thereby estimate the extent of recently transmitted disease (1, 9, 14).
In order to estimate the stability of the IS6110 RFLP patterns, studies have examined serial isolates collected from patients with persistent disease (3, 6, 15) and have demonstrated that the IS6110 banding pattern may change over time in a subset of these patients. When survival analysis was applied to the RFLP data collected from patients with persistent disease in The Netherlands, De Boer et al. (3) calculated that the half-life of the IS6110 banding pattern was on the order of 3.2 years. A similar rate was calculated for isolates originating from San Francisco, Calif., although this was restricted to isolates obtained at least 90 days apart, to exclude the possibility of variable patterns representing coinfection rather than change over time (3). In both data sets, the estimated rate of IS6110 pattern change was felt to be of the right magnitude to allow the use of IS6110 as a marker for molecular epidemiological studies (3).
However, the calculated rate of change may not be uniform in two groups of patients: those with multiple isolates at the time of diagnosis and those in whom a new pattern develops during or after treatment. In this study we have investigated the stability of the IS6110 banding pattern in a group of patients for whom more than one isolate of M. tuberculosis was collected during the course of disease or during retreatment. The results are discussed in the context of the applicability of calculating a rate of change in these patient groups and how such change may influence molecular epidemiological calculations.
| MATERIALS AND METHODS |
|---|
|
|
|---|
RFLP generation and GelCompar analysis. The entire culture representing all the possible clonal variants present in each isolate was used for the DNA isolation, and this DNA was genotypically classified according to the internationally standardized DNA fingerprinting protocol, using the IS-3' probe (10, 12). The Southern blot autoradiographs were normalized and the IS-3' bands were assigned using GelCompar software (version 4.1). The assignments were visually checked by two independent persons, and only bands with an intensity of >20% of the other bands were scored as representing IS6110-mediated evolutionary events (4). Replicative transposition was identified when the evolved strain showed an additional IS6110 hybridizing band, while deletion of an IS6110 element was identified when an IS6110 hybridizing band was absent from the evolved strain. Variation in the electrophoretic mobility of an IS6110 hybridizing band was classified as a band shift representing a mutational event in the chromosomal domain flanking the IS6110 element. Banding pattern changes suspected of being the result of partial digestion of methylated restriction sites (13) were excluded.
Study population. Patients with persistent disease from whom DNA fingerprints from more than one isolate were available were selected from the complete data set. This excluded isolates for which laboratory error (11), contamination, or loss of viability had been noted. Patients who had been assigned as being reinfected with a genotypically unrelated strain were treated as having two separate cases of disease if serial isolates were available for both infections. On average, three isolates were collected from each case of disease, and this was independent of the tuberculosis control program at the time of sample collection. Serial isolates collected from patients who presented with a genotypically identical or related strain after therapy were treated as a having a single case of disease. Selected patients were divided into two categories: (i) the group of patients infected with a strain which remained genotypically identical during the course of sampling was termed the "invariant group," and (ii) the group of patients infected with a strain which changed during the course of sampling was termed the "variant group." Change in the IS6110 banding pattern was identified by the appearance, disappearance, or electrophoretic shift of IS-3' hybridizing elements. To avoid a possible bias in the variant patient group, variant strains with more than four changes in the IS6110 banding pattern were considered to reflect reinfection rather than evolution and were excluded from all calculations. In addition, if any variant strain (independent of the number of changes in the IS6110 banding pattern) was found to be present in the community prior to appearance in the patient, the patient was excluded from the study. In such a case, the patient may have been reinfected from a community contact. This stringent inclusion criteria excludes the possibility that the IS6110 banding pattern may revert to a previous evolutionary state or that identical banding patterns could have evolved convergently. It is envisaged that the above exclusion criteria will lead to an overestimate of the stability of the IS6110 banding pattern.
The sampling interval (in days) for patients infected with invariant strains was calculated from the date when the first isolate was collected to the date when the last isolate was collected in each patient. The sampling interval (in days) for patients infected with variant strains was calculated from the date when the first isolate was collected to the date when the first variant isolate was identified (3). The overall rate of change was calculated using survival analysis according to the method described by de Boer et al. (3). Similarly, the early rate of change was calculated using cases with interisolate intervals of
90 days, while the late rate of change was calculated using cases with interisolate intervals of >90 days.
Statistical analysis. Fisher's exact test was used to test whether changes in the IS6110 banding pattern were associated with the number of IS6110 insertions present in the evolving strain. This test was also used to identify associations between previous therapy and the appearance of a variant strain. Chi-square analysis was used to identify differences between the study sample set and the sample set from The Netherlands (3).
| RESULTS |
|---|
|
|
|---|
Patients with serial isolates of M. tuberculosis. A total of 346 patients had two or more M. tuberculosis isolates, representing serial isolates from 351 active cases of tuberculosis. Analysis of the serial isolates by IS-3' DNA fingerprinting identified 335 cases (95.4%) in which the serial M. tuberculosis isolates remained genotypically constant over the sampling period (range, 0 to 2,203 days), termed invariant strains (Table 1). Serial isolates from 16 cases showed an IS6110 banding pattern change over the sampling period, termed variant strains. Two of the patients infected with a variant strain were excluded, as the variant strain was found to have been present in the community prior to appearing in the patient, thereby suggesting possible exogenous reinfection from a community source. The remaining 14 variant cases were sampled over a period of 0 to 1,141 days (Table 1). Comparing patients with invariant strains to those with variant strains did not reveal any differences in the demographic variables. However, a previous history of disease was weakly associated with the appearance of variant strains (Fisher's exact test, P = 0.045) (Table 1), which may suggest that the evolution of clonal variants occurred either during the latent disease phase after the previous infection or during the reactivation process.
|
Analysis of the RFLP data of the variant isolates showed that in 10 of 14 cases (71.4%) the IS6110 banding patterns changed by replicative transposition, while in 2 of 14 cases (14.3%) a hybridizing band was deleted. In 2 of 14 cases (14.3%), an electrophoretic shift was observed, suggestive of a chromosomal mutation. For four patients the variant strain disappeared from subsequent serial isolates, demonstrating incomplete clonal replacement. This suggests that two clonal variant populations may be present in the patient and that a single sputum specimen may not reflect the true genetic diversity of the bacterial population in the host.
Survival analysis (3) of these data estimated the overall half-life of the IS6110 banding pattern in our serial isolates to be 8.74 years (95% CI, 7.51 to 10.45 years) (Fig. 1). However, the survival plot suggests the presence of two distinct rates of change (Fig. 1). In 8 of the 14 patients (57%) both the initial and variant isolate were sampled within the first 20 days (median of 10 days), suggesting that the banding pattern change occurred early in the disease process, during the latent phase of a previous infection, or during reactivation. This was followed by only partial clonal replacement at the time of diagnosis or shortly thereafter. In the remaining six patients (43%), variant isolates were only sampled after >250 days, suggesting that these isolates evolved late in the disease phase. Since all strains sampled between 20 days and 250 days were invariant, this provides further evidence for the argument that these late variant strains evolve over the course of time, rather than being due to under-sampling at an earlier point. Taking into account these two populations of observations, survival analysis estimated a half-life of 0.57 years (95% CI, 0.45 to 0.77 years) for the early group and 10.69 years (95% CI, 8.22 to 15.31 years) for the late group (Fig. 1).
|
| DISCUSSION |
|---|
|
|
|---|
In this study, the proportion of IS6110 banding pattern changes identified in serial isolates collected from patients with persistent disease was comparable to those reported in The Netherlands (3) but lower than those reported in Germany (6) and San Francisco (15). The comparable proportion of variant isolates seen in this study and the study from The Netherlands suggests that the rate of change is independent of the incidence of disease. However, in keeping with the observations of other groups, the ability to observe change appears to depend on the number of IS6110 elements present in the precursor strain, as strains with fewer than five IS6110 elements were genotypically invariant.
On further analysis, the emergence of these variant strains does not appear to be uniform, suggesting the presence of two distinct populations. The first type of change would occur early in the disease process, prior to antituberculosis treatment. This change might then represent evolution influenced by active growth or adaptation to the new host environment. It is also possible that change was induced in certain previously treated patients either during the latent growth phase or during the subsequent reactivation process. The second type of change would occur in the late disease phase or after treatment. This lower rate of change may be a consequence of a lower growth rate of the bacilli. Similar results were seen in data from The Netherlands (3). Although we have calculated an early rate, it is highly probable that this rate will be an underestimate of the stability of the banding pattern, as most of the IS6110 band changes will have occurred prior to the sampling of the first isolate. This notion is supported by the low in vivo growth rate of M. tuberculosis, requiring extended periods before strain population replacement would be observed. Furthermore, if there are long delays in the period between onset of early active disease and the first collection of sputa, then it is possible that newly evolved strains may be missed and thereby excluded from the rate calculations.
Due to the uncertainty of when the early evolutionary events occurred, we attempted to calculate a rate of change for variant isolates appearing after 90 days (15). The rate estimate of 10.69 years for variant isolates appearing after 90 days was significantly lower than the 2 years calculated from the San Francisco data (3). Review of the different rate calculation methods (3, 8) shows that these methods assume that the rate of appearance of variant strains is constant over time. If this assumption is incorrect, as is suggested by the present study, then the rate calculated by these methods (3, 8) will be sensitive to the proportion of (invariant) cases with extended interisolate time intervals arising due to either noncompliance, treatment failure, or treatment interruption. Comparison of our data with those from The Netherlands shows that there is a significant difference in the proportion of patients with invariant isolates after 210 days (this study, 21.2%; study from The Netherlands, 6.6%; P < 0.0001). This implies that if treatment efficacy is higher, the proportion of serial isolates with long intervals will be smaller, which may result in a higher rate of change estimate. Consequently, the accuracy of rate calculations is limited by the local epidemiology, and therefore rate estimates will be proportional to the effectiveness of the tuberculosis control program.
The epidemiological importance of calculating an evolutionary rate for isolates collected from patients with persistent disease is unclear, given that change is observed only in a very small and unique subset of the patients analyzed (3, 6, 15). Therefore, these data cannot be easily extrapolated to the overall study population. Furthermore, the true extent of genetic diversity within a study population may be masked by incomplete clonal replacement at the time of sputum collection, as extremely low numbers of clonal variants in the sputum specimens will not be detected due to the limited sensitivity of the DNA fingerprinting method. Therefore, strain variants will not be identified for patients from whom inadequate numbers of sputum specimens have been collected. Accordingly, it may be difficult to identify chains of transmission based only on genotypic identity (clustering of identical IS6110 banding patterns).
Although this study has highlighted the limitations of this sample set, we conclude that this patient set may provide information on how growth rate and/or antituberculosis therapy could influence the IS6110-mediated evolutionary rate. Furthermore, we do not exclude the possibility that change in the IS6110 banding pattern may be more frequent in other patient groups. To gain a greater insight into the rate of IS6110 change it will be essential to determine how transmission influences the evolutionary process (7).
| ACKNOWLEDGMENTS |
|---|
E. Engelke, S. Carlini, and M. De Kock are thanked for their technical assistance.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Antimicrob. Agents Chemother. | Clin. Microbiol. Rev. |
|---|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|