|
|
||||||||
1 School of Chemistry, Faculty of Engineering and Physical Sciences, The University of Manchester, PO Box 88, Sackville Street, Manchester M60 1QD, UK
2 TB Research Group, Veterinary Laboratories Agency (VLA), Weybridge, New Haw, Addlestone, Surrey KT15 3NB, UK
Correspondence
Catherine L. Winder
C.Winder{at}manchester.ac.uk
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
The ideal microbial typing technique would be easy to use, automated, with intuitive nomenclature and have good inter-laboratory reproducibility (Roring et al., 2004
). Restriction fragment length polymorphisms (RFLP) (Goyal et al., 1997
) of the IS6110 element are generally used in discriminating members of the M. tuberculosis complex; however, this technique has limited use in M. bovis due to the low copy number of IS6110. Therefore, alternative genetic typing methods have been used in the analysis of M. bovis, with spacer oligonucleotide typing, spoligotyping (Kamerbeek et al., 1997
), the current method of choice. Spoligotyping focuses on identifying repeat variation in the direct repeat (DR) locus, which is composed of multiple, virtually identical 36 bp repeats interspersed with unique DNA spacer sequences (direct variant repeats; DVRs) of similar size. Different spoligotype patterns are identified by the presence or absence of DVRs. In the UK ten major spoligotypes account for
90 % of all isolates, with the two dominant spoligotypes 9 and 17 representing 65 % of all isolates (Inwald et al., 2002
). The dominance of two main types and the tendency of the spoligotype pattern to alter slowly with time (Smith et al., 2003
) reduce the usefulness of spoligotyping as an epidemiological tool.
Alternative genotypic typing techniques identifying sets of tandem repeats, such as exact tandem repeats (ETRs), mycobacterial interspersed repetitive units (MIRUs) and variable number tandem repeats (VNTRs) have been described as independent tools for the discrimination of members of the M. tuberculosis complex (Roring et al., 2004
). VNTRs have been successfully used in the discrimination of M. bovis isolates, and the combination of spoligotyping with VNTR analysis vastly improves the levels of discrimination (Gibson et al., 2004
). Indeed, our previous work has shown that VNTR analysis can successfully subdivide spoligotype 9 into clusters of related strains; however, the VNTR technique requires optimization and the level of discrimination depends upon the loci applied and test panel chosen for the study. Different degrees of discrimination may be appropriate for different studies, with some loci useful for population-based studies whereas more discriminating loci may be required to monitor outbreak analysis (Roring et al., 2004
).
In this study we investigated the ability of the vibrational spectroscopic technique Fourier-transform infrared spectroscopy (FT-IR) to produce robust biochemical fingerprints' of M. bovis isolates. Numerous advantages of using this method of discrimination over traditional microbiological and molecular techniques are documented; these include speed, minimal sample preparation, automation and relatively low expense. The use of multivariate statistical techniques to analyse the infrared spectra has proved successful in the discrimination of the bacteria to subspecies levels (Maquelin et al., 2003
; Naumann et al., 1991
; Timmins et al., 1998
; Winder et al., 2004
). To date, the discrimination of members of the Mycobacterium genus by FT-IR has not been reported in the literature. Moreover, we show that the FT-IR method successfully clustered strains of the same molecular type, suggesting that molecular types share distinct phenotypic characteristics that were discriminated by FT-IR. The implications of this result are discussed.
| METHODS |
|---|
|
|
|---|
|
Reflectance-based FT-IR.
A 10x10 cm aluminium plate was rinsed with acetone and dried at 50 °C for 10 min prior to use. FT-IR analysis was initially performed on the plate without samples to provide a reference reading for each well. The plate was loaded onto the motorized stage of a reflectance TLC accessory (Winder et al., 2004
) attached to a Bruker IFS28 FT-IR spectrometer (Bruker Spectrospin). This was equipped with a mercury-cadmium-telluride (MCT) detector, which was cooled with liquid N2. A 5 µl aliquot of each sample was evenly applied in triplicate onto the 400-well plate (the so-called machine replicates). The plate was then dried at 50 °C for 30 min before FT-IR analysis was performed. Spectra were collected over the wavelength range of 4000600 cm1 under the control of an IBM-compatible computer programmed with Opus 2.1 running under IBM OS/2 Warp, which was provided by the manufacturers. Spectra were acquired with a resolution of 4 cm1 and 256 spectra were co-added and averaged. The collection time for the spectra was approximately 10 s per sample. The spectra are displayed in terms of absorbance, which was calculated from the reflectance-absorbance spectra using Opus software.
Transmission-based FT-IR.
Prior to use a 96-well zinc selenide plate was rinsed with 2-propanol and deionized water and allowed to dry at room temperature. Aliquots (15 µl) of each bacterial sample were evenly applied in triplicate onto the plate and dried at 50 °C for 30 min. The plate was loaded onto a motorized microplate module HTS-XT attached to an Equinox 55 module (Bruker Optics). The motorized module of this instrument introduces the plate into the airtight optics of the instrument, in which tubes of desiccant are contained to remove moisture (Harrigan et al., 2004
). A deuterated triglycine sulphate (DTGS) detector was employed for transmission measurements of the samples to be acquired. Spectra were collected over the wavelength range 4000900 cm1 (a reduced range is used in comparison to the reflectance measurements due to the zinc selenide focusing element) under the control of a computer programmed with Opus 4, operated under MS Windows 2000. Spectra were acquired at a resolution of 4 cm1 and 64 spectra were co-added and averaged to improve the signal-to-noise ratio. The collection time for each spectrum was approximately 1 min and the spectra were displayed in terms of absorbance.
Preprocessing.
The ASCII data were imported into Matlab version 6 (The MathWorks). To minimize problems arising from baseline shifts, Matlab was used to correct for CO2 vibrations [the CO2 peaks at 24032272 cm1 and 683656 cm1 (when present) were removed and filled with a trend] and windows of the spectra likely containing H2O vibrations were smoothed with a window of 55 cm1 to reduce noise (40003615 cm1 and 19981763 cm1). The spectra were normalized such that the smallest recorded absorbance was set to 0 and the highest was set to 1 for each spectrum and then the second derivatives (Savitzky & Golay, 1964
), with a window of 9), were used for cluster analysis.
Cluster analyses.
The unsupervised clustering method principal components analysis (PCA; Jolliffe, 1986
) was performed on the spectra to reduce the dimensionality of the multivariate data whilst preserving the variance, prior to discriminant function analysis (DFA). DFA is a supervised technique that discriminates between groups on the basis of the retained principal components with a priori knowledge of which spectra were replicates. As this process uses information on the machine replicates for each isolate it does not bias the results in any way (MacFie et al., 1978
; Windig et al., 1983
). Hence the cluster analysis would not bias the clustering of the isolates according to their respective spoligotypes because the a priori knowledge is based on the isolate replicates. DFA was programmed to minimize within group variance and maximize between group variance. Finally, hierarchical cluster analysis (HCA) was used to construct a dendrogram from the a priori group centres in the PC-DFA space, using scaled Euclidean distances as described by Goodacre et al. (1998)
, and dendrograms were produced using the mean linkage clustering algorithm (Manly, 1994
).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
PC-DFA and HCA were used to produce ordination (scatter) plots and dendrograms respectively. Whilst both the ordination plots and dendrograms were used to investigate the discriminatory ability of FT-IR, for brevity and clarity the dendrograms are predominantly shown for reflection-based FT-IR (Fig. 2
) and transmission-based measurements (see Fig. 4
).
|
|
Reflectance-based FT-IR
The resulting dendrogram for the reflectance FT-IR illustrates the clear discrimination of spoligotypes 20, 25 and 35 (Fig. 2
). The PC-DFA ordination plot of the first analysis of the composite dendrogram is shown in Fig. 3
(a), to illustrate the clustering of spoligotypes 20, 25 and 35. Fig. 3(b)
indicates the mean centres of the a priori groups and the circles indicate a 95 %
2 confidence region around the point. The overlapping circles are shown for spoligotypes 20, 25 and 35; the circles representing the other spoligotypes are not shown for clarity of the figure. The analysis shows the isolates representing the chosen spoligotypes cluster within the confidence regions. Removal of the three discriminated spoligotypes and the subsequent reanalysis of the data demonstrated the dispersal of spoligotype 9 throughout the HCA space (data not shown). The successive removal of spoligotype 9 and reanalysis of the data illustrates the recovery of spoligotypes 10 and 22. A fourth analysis was then performed on the remaining data following the removal of spoligotypes 10 and 22. In this final dendrogram, the representative samples of spoligotype 17 were recovered together; however, the isolates of spoligotypes 12 were recovered with M. bovis spoligotypes 11 and 13 isolates.
|
Comparison of the two FT-IR methods
In the initial analysis of both the reflectance and transmission data, three distinct branches were observed according to three spoligotypes. Spoligotypes 25 and 35 were recovered in both analyses but the third branch varied according to the different FT-IR techniques. The relationship between the spoligotypes analysed in this study illustrates spoligotypes 25 and 35 to be closer to the ancestor, in which all the spacers are present (Fig. 5
). Spoligotype 10 was resolved for the transmission data whereas spoligotype 20 was initially resolved in the analysis of the reflectance data. The isolates representing spoligotype 9 were intermixed with the samples representing the other spoligotypes in the analysis for both the reflectance and transmission data. The evolutionary history indicates that a number of spoligotypes (10, 11, 12, 13, 17 and 22) evolved after spoligotype 9 (Fig. 5
). The remaining spoligotypes (with the exceptions of 12 and 13) were resolved in the third analysis for the transmission data. A fourth analysis was necessary to recover the majority of the reflectance-collected samples according to spoligotype. In both analyses, the isolates representing spoligotypes 11, 12 and 13 were poorly resolved in the analysis. The representative samples of spoligotypes 12 and 13 (transmission data) were analysed separately; the isolates of spoligotype 12 were recovered in a cluster, but those from spoligotype 13 appeared diffuse in the analysis.
|
Analysis of UK dominant M. bovis isolates
As reported above, spoligotypes 9 and 17 are the two dominant spoligotypes in the UK, accounting for 65 % of all isolates from cattle (Inwald et al., 2002
). We therefore focused the final analysis on just type 9 and 17 isolates. The resulting PC-DFA scores plot is shown in Fig. 6
, illustrating the differentiation of five clusters; this was confirmed by the 95 %
2 confidence region. It can be clearly seen that a tight cluster of the isolates representing spoligotype 17 was observed, which is perhaps not surprising as these are a single VNTR group (7555*33.1), whilst four of the five isolates representing spoligotype 9 in the analysis were dispersed throughout the PC-DFA space [note that each of the spoligotype 9 isolates was of a different VNTR group (Table 1
)]. This clearly reflects what was observed in the dendrograms (Figs 2 and 4![]()
) discussed above. A similar pattern was observed in the analysis of the reflectance data (data not shown). The remaining representative of spoligotype 9 (9d, VNTR group 7554*33.1) clustered with the isolates of spoligotype 17 (VNTR group 7555*33.1). This may possibly result from the close phylogenetic relationship of these isolates according to the VNTR phylogeny as seen by us (to be reported elsewhere).
|
Hence, the analysis of the strains using FT-IR allowed us to generate clusters that showed excellent agreement with clustering performed using genetic markers; this suggests that strains of the same molecular type have conserved phenotypic traits. This is intriguing when considered in the light of our knowledge of the population structure of M. bovis in Great Britain, which found that the frequency of recovery of spoligotypes 9 and 17 isolated in the UK cannot be explained by random mutation and drift (Smith et al., 2003
). Instead these strains have undergone a clonal expansion, rising to a high frequency in the population either through a favourable mutation or by colonization of a new geographical location (Smith et al., 2003
). The results reported here suggest that the success of M. bovis clonal groupings is due to favourable mutations in these strains that may affect virulence, transmissibility, or survival outside the host. Our FT-IR data therefore provide us with a unique opportunity to examine the relationship between clonal groups and phenotype.
Conclusions
The data collected by transmission FT-IR illustrated better differentiation of the M. bovis isolates according to spoligotypes than that collected by reflectance FT-IR. The phenotypic analysis in this study indicates homogeneity within spoligotypes 10, 17, 20, 22, 25 and 35, while spoligotypes 11, 12 and 13 appeared diffuse and difficult to cluster in the analysis. Tight clustering of molecular types detected by metabolic fingerprinting suggests that these clonal groupings share distinct phenotypic traits that may have a significant role in their success as pathogens. Isolates of spoligotype 9 appear very heterogeneous relative to other spoligotypes, an observation that was supported by VNTR analysis and reflects the wide variation within this group. Therefore this study describes the first use of FT-IR for the discrimination of M. bovis strains. As well as being a robust and high-throughput method that produces strain clusters reflecting those generated by genetic methods, FT-IR has the added benefit of providing the first insights into phenotypegenotype links in M. bovis clones. This study therefore sets the baseline for a more in-depth phenotypic analysis of M. bovis molecular types.
| ACKNOWLEDGEMENTS |
|---|
| REFERENCES |
|---|
|
|
|---|
Clifton-Hadley, R. S., Wilesmith, J. W., Richards, M. S., Upton, A. & Johnston, S. (1995). The occurance of Mycobacterium bovis infection in cattle in and around an area subject to extensive badger (Meles meles) control. Epidemiol Infect 114, 179193.[Medline]
Delahay, R. J., Cheeseman, C. L. & Clifton-Hadley, R. S. (2001). Wildlife disease reservoirs: the epidemiology of Mycobacterium bovis infection in the European badger (Meles meles) and other British mammals. Tuberculosis 81, 4349.
Gibson, A. L., Hewinson, G., Goodchild, T., Watt, B., Story, A., Inwald, J. & Drobniewski, F. A. (2004). Molecular epidemiology of disease due to Mycobacterium bovis in humans in the United Kingdom. J Clin Microbiol 42, 431434.
Goodacre, R., Timmins, E. M., Burton, R., Kaderbhai, N., Woodward, A. M., Kell, D. B. & Rooney, P. J. (1998). Rapid identification of urinary tract infection bacteria using hyperspectral whole-organism fingerprinting and artificial neural networks. Microbiology 144, 11571170.[Abstract]
Goyal, M., Saunders, N. A., VanEmbden, J. D. A., Young, D. B. & Shaw, R. J. (1997). Differentiation of Mycobacterium tuberculosis isolates by spoligotyping and IS6110 restriction fragment length polymorphism. J Clin Microbiol 35, 647651.[Abstract]
Griffin, J. M., Williams, D. H., Kelly, G. E., Clegg, T. A., O'Boyle, I., Collins, J. D. & More, S. J. (2005). The impact of badger removal on the control of tuberculosis in cattle herds in Ireland. Prev Vet Med 67, 237266.[CrossRef][Medline]
Harrigan, G. G., LaPlante, R. H., Cosma, G. N., Cockerell, G., Goodacre, R., Maddox, J. F., Luyendyk, J. P., Ganey, P. E. & Roth, R. A. (2004). Applications of high-throughput Fourier-transform infrared spectroscopy in toxicology studies: contribution to a study on the development of an animal model for idiosyncratic toxicity. Toxicol Lett 146, 197205.[CrossRef][Medline]
Inwald, A., Hinds, J., Dale, J., Palmer, S., Butcher, P., Hewinson, R. G. & Gordon, S. V. (2002). Microarray-based comparative genomics: genome plasticity in Mycobacterium bovis. Comp Funct Genomics 3, 342344.[CrossRef]
Jarvis, R. M. & Goodacre, R. (2004). Discrimination of bacteria using surface-enhanced Raman spectroscopy. Anal Chem 76, 4047.[Medline]
Jolliffe, I. (1986). Principal Component Analysis. New York: Springer-Verlag.
Kamerbeek, J., Schouls, L., Kolk, A. & 8 other authors (1997). Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol 35, 907914.[Abstract]
Lopez-Diez, E. C. & Goodacre, R. (2004). Characterization of microorganisms using UV resonance Raman spectroscopy and chemometrics. Anal Chem 76, 585591.[Medline]
MacFie, H. J. H., Gutteridge, C. S. & Norris, J. R. (1978). Use of canonical variates in differentiation of bacteria by pyrolysis gas-liquid chromatography. J Gen Microbiol 104, 6774.[Medline]
Manly, B. F. J. (1994). Multivariate Statistical Methods. London: Chapman & Hall.
Maquelin, K., Kirschner, C., Choo-Smith, L. P. & 7 other authors (2003). Prospective study of the performance of vibrational spectroscopies for rapid identification of bacterial and fungal pathogens recovered from blood cultures. J Clin Microbiol 41, 324329.
Naumann, D., Helm, D. & Labischinski, H. (1991). Microbiological characterizations by FT-IR spectroscopy. Nature 351, 8182.[CrossRef][Medline]
Roring, S., Scott, A. N., Hewinson, R. G., Neill, S. D. & Skuce, R. A. (2004). Evaluation of variable number tandem repeat (VNTR) loci in molecular typing of Mycobacterium bovis isolates from Ireland. Vet Microbiol 101, 6573.[CrossRef][Medline]
Savitzky, A. & Golay, M. J. E. (1964). Smoothing and differentiation of data by simpli®ed least squares procedures. Anal Chem 36, 16271633.[CrossRef]
Smith, N. H., Dale, J., Inwald, J., Palmer, S., Gordon, S. V., Hewinson, R. G. & Smith, J. M. (2003). The population structure of Mycobacterium bovis in Great Britain: clonal expansion. Proc Natl Acad Sci U S A 100, 1527115275.
Timmins, E. M., Howell, S. A., Alsberg, B. K., Noble, W. C. & Goodacre, R. (1998). Rapid differentiation of closely related Candida species and strains by pyrolysis mass spectrometry and Fourier transform-infrared spectroscopy. J Clin Microbiol 36, 367374.
Winder, C. L., Carr, E., Goodacre, R. & Seviour, R. (2004). The rapid identification of Acinetobacter species using Fourier transform infrared spectroscopy. J Appl Microbiol 96, 328339.[CrossRef][Medline]
Windig, W., Haverkamp, J. & Kistemaker, P. G. (1983). Interpretation of sets of pyrolysis mass spectra by discriminant analysis and graphical rotation. Anal Chem 55, 8188.[CrossRef]
Received 10 March 2006;
revised 1 June 2006;
accepted 5 June 2006.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| INT J SYST EVOL MICROBIOL | MICROBIOLOGY | J GEN VIROL |
| J MED MICROBIOL | ALL SGM JOURNALS | |