Previous Article | Next Article ![]()
Journal of Clinical Microbiology, September 2005, p. 4558-4566, Vol. 43, No. 9
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.9.4558-4566.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Veterinary and Agrochemical Research Centre, Groeselenberg 99, 1180 Brussels, Belgium,1 Faculty of Veterinary Medicine, Ghent University, Salisburylaan 133, 9820 Merelbeke, Belgium,2 Department of Clinical Chemistry, Microbiology, and Immunology, Ghent University Hospital, De Pintelaan 185, 9000 Ghent, Belgium,3 Department of Microbiology, University of Brussels (VUB) Hospital, Laarbeeklaan 101, 1090 Brussels, Belgium4
Received 15 December 2004/ Returned for modification 14 February 2005/ Accepted 2 June 2005
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
To correctly differentiate all these species, a universal and fast identification technique would be extremely useful. Some promising methods have already been described (see, e.g., reference 25), but they do not yield digitized data, making exchange between laboratories difficult. An optimized tRNA gene PCR technique, originally described by Welsh and McClelland (27, 48), has been shown to be useful for correct and reproducible identification of very diverse bacterial species when combined with high-resolution electrophoresis (3-5, 10, 11, 23, 45). The technique is based on the amplification of spacer regions between tRNA genes using consensus tRNA gene primers. The amplified products are separated by electrophoresis for exact sizing, and the resulting species-specific peak profiles are subsequently archived in a database. Profiles obtained from an unknown sample can be compared with this data set while not-yet-included and/or newly described species can be added to expand the database further. We investigated the potential of this technique to correctly identify a large number of Acholeplasma, Mycoplasma, and Ureaplasma species.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Identification of strains. Most strains (i.e., MYP10 to MYP58 and MYP65 to MYP74) had been identified previously by means of phenotypical characteristics and the growth precipitation test using absorbed rabbit antisera (13, 36). After identification, the strains had been lyophilized in the presence of 20% sterile milk and stored at 4°C. U. urealyticum, U. parvum, and some of the M. hominis isolates were previously identified by their characteristic growth on A7 differential agar medium and by their ability to hydrolyze urea and arginine, respectively. Due to the numerous discrepancies with the results obtained in this study by means of tRNA gene PCR, most of the strains were reidentified. The identity of M. hyopneumoniae, M. hyorhinis, and M. flocculare was confirmed by a specific PCR (40). Also, the M. genitalium samples were identified by two specific PCR tests as described previously by Jensen et al. (19, 20). Most other isolates were reidentified using amplified rRNA gene restriction analysis (ARDRA) (unpublished data) and/or sequence analysis of the 16S rRNA genes (12, 46) (Table 1).
tRNA gene PCR and cluster analysis. tRNA gene PCR was performed as described previously using primers T5A (5'-AGTCCGGTGCTCTAACCAACTGAG) and T3B (5'-AGGTCGCGGGTTCGAATCC) (4, 27). Cluster analysis of the obtained tRNA gene PCR fingerprints was carried out by calculating a distance matrix using the differential base pair (DBP) algorithm (3) with a tolerance of 1.2 bp and including all peaks (i.e., no noise subtraction) from 50 to 500 base pairs in length. The DBP clustering algorithm was used to calculate similarity by taking the average of the two results that are obtained by dividing the number of tRNA gene-intergenic spacers in common between two strains by the total number of spacers of one strain, respectively, of the other strain. A similarity tree was constructed using the unweighted-pair group method using average linkages (PHYLIP, V3.6; J. Felsenstein, Department of Genome Sciences, University of Washington, Seattle, Wash.) and visualized using Treeview V1.6.6 (29).
Construction of a digital library. A digital library composed of consensus library entries was constructed. Each consensus library entry contained only peaks (amplified tRNA gene-intergenic spacers) present in all sample files of a particular species. Identification was carried out by comparing the fingerprint of a strain with all entries of the constructed library using the DBP algorithm. For identification, the DBP algorithm takes only peaks that are mutually present in the sample file and in the library entry into account, discarding the peaks only present in the sample file. For example, comparison of an unknown strain with 15 peaks, of which 10 are identical to all 10 peaks of a library entry, gives a similarity of 100%, as would have been the case for an unknown strain with 10 peaks, all identical to the 10 peaks of the library entry. An unknown strain with 8 peaks, of which all 8 are identical to 8 of the 10 peaks of the library entry, gives a value of only 80%. An unknown with 15 peaks, of which 8 are identical to an entry with 10 peaks, gives an identity of 80% as well. This method is (by experience) better suited for identification of unknown patterns, since intraspecific variability of the peaks with low intensity is better compensated for. To differentiate between resembling, but distinct, tRNA gene PCR patterns, the algorithm allows us to take absent peaks into account by adding a minus in front of the absent peak. The software counts the absence of the peak as a positive match, increasing the similarity score and increasing the reliability of the identification result.
Nucleotide sequence accession number. The 16S rRNA gene sequence for Acholeplasma sp. strain MYP14 was deposited in GenBank under accession number AY785356.
| RESULTS |
|---|
|
|
|---|
|
|
The three strains MYP31 to MYP33, received as M. columbinum, had identical patterns, clustering separately. One strain received as M. gallinaceum (MYP38) and two of the six strains received as M. gallisepticum (MYP41 and MYP42) had identical, unique tRNA gene PCR patterns and were identified by ARDRA as M. gallinarum.
Six strains (MYP13, MYP39 to MYP42, and MYP71) were received as M. gallisepticum. As explained, MYP41 and MYP42 were shown to be M. gallinarum. No amplification products were obtained from strains MYP13, MYP40, and MYP71, identified as M. gallisepticum according to ARDRA or sequence analysis. The tRNA gene PCR pattern of strain MYP39 showed a high similarity to members of the M. bovis-M. agalactiae cluster, although the tRNA gene pattern profile clearly had many additional peaks. Sequence analysis showed a mixed profile, confirming the contamination of this sample, which was therefore left out for further analysis.
The one M. glycophilum strain, MYP43, of which the identification was confirmed by ARDRA, was indistinguishable from strain MYP52, received as M. pullorum, which was identified by ARDRA as M. glycophilum.
The one M. lipofaciens strain (MYP49), confirmed as such by ARDRA, had a very specific tRNA gene PCR pattern. Two strains (MYP52 and MYP53) were received as M. pullorum. One (MYP52) was identified as M. glycophilum by ARDRA as discussed above. The other one (MYP53) was identified by sequence analysis as M. columborale, a species that was not included initially, which is consistent with the fact that this strain had a unique tRNA gene PCR pattern.
The one strain received as M. synoviae (MYP58) had a tRNA gene PCR pattern that strongly resembled that of the sole M. neurolyticum strain (MYP50), an identification that was confirmed by ARDRA.
Bovine, caprine, and ovine mollicutes. Four strains (MYP16 to MYP19) were received as M. agalactiae. The DNA extract obtained from M. agalactiae strain MYP16 was found to yield a separate tRNA gene PCR fingerprint (Fig. 1). This strain was later identified with ARDRA as M. arginini. The two genuine M. agalactiae strains (MYP17 and MYP19) had a very similar tRNA gene profile and clustered together with strain MYP18, which was later identified as M. bovis by means of ARDRA.
M. bovigenitalium strain MYP20 had a unique tRNA gene PCR pattern, and its identification could be confirmed by 16S rRNA gene sequence analysis.
Strains MYP22 to MYP23 and MYP67 to MYP68 had previously been identified as M. bovis. MYP68 clustered together with MYP66, which was received as M. bovirhinis and was indeed shown by ARDRA to be M. bovirhinis. The remaining three M. bovis strains clustered together with M. agalactiae strains MYP17 and MYP19.
Four strains were received as M. dispar (MYP34 to MYP37). MYP34 and MYP37 were clustered by tRNA gene PCR in the M. bovis-M. agalactiae group, a finding substantiated by ARDRA, which identified both strains as M. bovis. This identification was confirmed for MYP37 by sequencing. The remaining M. dispar strains, MYP35 and MYP36, clustered separately (Fig. 1) and were shown to be genuine M. dispar by ARDRA.
The so-called M. mycoides cluster comprises six species or subspecies of closely related mycoplasmas. The type strains of these ruminant mycoplasmas were included for analysis (MYP75 to MYP80). M. capricolum subsp. capricolum (MYP80) and in particular M. capricolum subsp. capripneumoniae (MYP76) clustered most closely to a Mycoplasma sp. bovine group 7 strain (MYP77). These profiles showed up to 30 peaks, with over 20 peaks in common. Likewise, the tRNA gene profiles of the three M. mycoides isolates showed several small peaks and clustered together, with M. mycoides subsp. capri (MYP78) and M. mycoides subsp. mycoides LC (large colony) (MYP79) most closely related to each other.
The one strain received as M. ovipneumoniae (MYP51) clustered within the M. bovis-M. agalactiae group (Fig. 1) and was identified as M. bovis by 16S rRNA gene sequence analysis.
The four M. putrefaciens strains (MYP54 to MYP57), confirmed as such by ARDRA, had very identical and characteristic tRNA gene PCR patterns.
Human mollicutes. Five samples were received as M. genitalium (MYP106 to MYP110). Strain MYP107 did not yield a fingerprint. Strikingly, there were two different groups observed within M. genitalium (Table 2) which did cluster separately but close to each other.
Three M. hominis cultures were included in this study. Two pure M. hominis isolates (MYP111 and MYP112), identified by means of 16S rRNA gene sequencing, had identical and specific tRNA gene PCR patterns composed of only two peaks (of 151.9 and 226.6 bp). One culture (MYP081) that was positive for M. hominis on A7 agar plates was clearly contaminated, since additional peaks of 56.0, 144.2, and 280.7 bp, shown in this study to be characteristic for U. urealyticum, were present. Only the tRNA gene PCR fingerprints of the pure isolates MYP111 and MYP112 were included in the cluster analysis.
One M. orale strain (MYP115) was obtained from a culture collection (National Collection of Type Cultures, United Kingdom) and had a very specific tRNA gene PCR profile, with over 25 characteristic peaks.
The tRNA gene PCR patterns of the three strains received as M. pneumoniae (MYP72 to MYP74), confirmed as such by ARDRA, were almost identical and were characterized by very short spacers (usually no longer than 77 bp).
The two M. salivarium strains (MYP113 and MYP114) had been isolated during studies of the complex microflora of tonsils and teeth, and their identity had been established by 16S rRNA sequencing. Their nearly identical and highly characteristic tRNA gene patterns clustered together.
A total of 24 Ureaplasma strains, received as U. urealyticum or U. parvum (MYP82 to MYP105), were included. For eight strains, no amplification could be obtained. Of the remaining samples, four (MYP82, MYP83, MYP104, and MYP105) had been received as contaminated with M. hominis, which was also apparent (as explained above when the results for M. hominis were presented) from the mixed tRNA gene PCR pattern that was obtained and which contained peaks characteristic to both species. These tRNA gene PCR fingerprints were not included in the cluster analysis. All other strains of both Ureaplasma species gave similar tRNA gene PCR patterns and clustered together.
Murine mollicutes. The one strain (MYP50) that was received as M. neurolyticum and confirmed as such by ARDRA had a tRNA gene PCR pattern that clustered together with that of strain MYP58, received as M. synoviae, but was shown to be M. neurolyticum, as described above.
Porcine mollicutes. The one strain that was received as Acholeplasma granularum (MYP15) was confirmed as such by ARDRA and could be identified easily by tRNA gene PCR, since its pattern was highly characteristic.
The three M. flocculare strains (MYP1 to MYP3) had a very similar pattern that made it possible to differentiate this species from all other species.
Five strains (MYP7 to MYP9, MYP44, and MYP45) were received as M. hyopneumoniae. A specific PCR identified MYP44 as M. hyorhinis, while MYP45 was identified as A. laidlawii as described below. MYP7, MYP8, and MYP9 had identical and characteristic tRNA gene PCR fingerprints, and their identity as M. hyopneumoniae was confirmed by ARDRA and by a specific PCR.
Strains MYP4 to MYP6, received as M. hyorhinis, were very much alike and had a typical pattern. The genuine M. hyorhinis strains clustered together with MYP44, which had been shown to be M. hyorhinis as well (see above).
Six strains (MYP59 to MYP64), received as M. hyosynoviae and confirmed as such by ARDRA, were very much alike with regard to their characteristic tRNA gene PCR patterns.
Other mollicutes. Five strains (MYP10 to MYP12, MYP14, and MYP65) were received as A. laidlawii. MYP14 had a unique tRNA gene PCR pattern and clustered separately. Sequence analysis showed significant differences with other known Acholeplasma spp., and the 16S rRNA gene sequence was submitted to GenBank. The other four strains clustered together with strain MYP45, which was received as M. hyopneumoniae but was also identified as A. laidlawii by means of ARDRA.
tRNA gene PCR-based identification. A digital library was constructed as described above. For A. laidlawii and M. genitalium, different tRNA gene PCR profiles were apparent, and for these species, two different consensus patterns were included in the library. All individual fingerprints (sample files including all peaks) were compared with this library using the similarity calculation designated DBP. The two M. agalactiae strains were indistinguishable from the eight M. bovis strains, having tRNA gene spacers with lengths of 57.4, 61.5, 67.3, 70.0, 78.6, 131.4, 144.2, 151.8, 159.5, and 257.4 bp in common. All other samples were identified correctly. Although the U. urealyticum and U. parvum strains grouped together during cluster analysis, they could clearly be distinguished on the basis of their specific tRNA gene PCR pattern. The strains MYP85 (serovar 1), MYP90 (serovar 6), and MYP98 (serovar 14) belong to U. parvum. These three stains, together with MYP103, for which serovar determination had not been carried out, had a peak of 279.5 bp (standard deviation, 0.1 bp) in common, whereas the U. urealyticum strains MYP88 (serovar 4), MYP89 (serovar 5), MYP91 (serovar 7), MYP93 (serovar 9), MYP94 (serovar 10), MYP95 (serovar 11), MYP96 (serovar 12), and MYP100 (serovar not determined) had a peak of 280.7 bp (standard deviation, 0.03 bp) in common.
| DISCUSSION |
|---|
|
|
|---|
In general, over 10 different peaks were visible in the tRNA gene PCR fingerprints of different Mollicutes species, which is a high number compared to those of most other bacteria. This is a somewhat unexpected finding, since mollicutes have only a limited number of tRNA genes. In view of the fact that the tRNA gene PCR technique applied with the ABI310 genetic analyzer (Applied Biosystems, Calif.) only takes into account small PCR fragments of less than 500 bp, the close proximity of tRNA genes or the possibly high rate of tRNA gene-like sequences (14) may partly explain these results. The presence of this high number of peaks makes the technique well applicable for the identification of a complex and diverse class of bacteria like the Mollicutes.
In addition, the tRNA gene PCR technique has been shown to be very reproducible (4). This is also apparent from our results, since nearly identical fingerprints were obtained even for strains that were received from different laboratories and isolated on different dates. Thanks to interlaboratory reproducibility (4) and digitized output data, the tRNA gene PCR fingerprints of more species and subspecies can be collected from different laboratories and published in a shared online database.
No amplified PCR fragments were observed in 12 cases. The reason for failure is unclear, but in the case of the one M. genitalium sample and the eight Ureaplasma samples, this is probably due to poor DNA quality or the possible presence of PCR inhibitors (1), since other strains of the same species were amplified without problems. For the three M. gallisepticum strains, the reason is less clear, since none of the strains yielded a tRNA gene PCR pattern despite the high quality of the DNA samples used, as can be concluded from the efficient amplification of the 16S rRNA gene. Alignment and comparison of all tRNA gene sequences of the fully sequenced M. gallisepticum strain R with those of other known Mycoplasma tRNA gene sequences did not reveal any exceptional differences (30). Also, the arrangement of tRNA gene clusters in the M. gallisepticum genome was very much like that of other mollicute species (38, 43).
For all other cases, correct identification was obtained, except for the indistinguishable tRNA gene patterns of M. bovis and M. agalactiae. The high similarity of the 16S rRNA gene sequences of both species has already been reported (24), and the close relatedness of both species is also reflected by the fact that M. bovis was first considered as a subspecies of M. agalactiae. Later studies involving DNA homology and serology led to the proposal that M. bovis should not be regarded as a subspecies of M. agalactiae but should be regarded as a distinct species (2). In contrast with the common peaks observed for both species, several minor differences were noted between individual tRNA gene PCR patterns. Importantly, genetic variability between M. bovis isolates was already demonstrated by McAuliffe et al. using several other molecular DNA techniques (26). Those authors observed two distinct groups of M. bovis isolates, and an earlier report also showed the existence of two distinct groups of M. agalactiae isolates based on antigenic profiles (39). Whether the presence of the minor peaks coincides with these subgroups is not yet known.
Still, tRNA gene PCR enabled differentiation between subgroups of a number of other species, as was the case for the strains of A. laidlawii. The two distinct tRNA gene PCR profiles may indicate the existence of two different genomic groups for this species. This has also been indicated by earlier pulsed-field gel electrophoresis results showing different genome sizes for two A. laidlawii strains (32) and especially by nucleic acid hybridization studies that demonstrated extensive genomic variation between different strains (42).
Based on tRNA gene PCR, we also observed genotypic diversity within M. genitalium. A recent report demonstrated the presence of a number of different M. genitalium genotypes (22), but possible correspondence with the different tRNA gene PCR groups established in this study remains to be studied.
In the dendrogram obtained (Fig. 1), some groups, like the M. hominis taxon, clustered close together, similar to phylogenetic data based on 16S rRNA gene sequence analyses, while other groups, like the M. neurolyticum taxon, were scattered throughout the dendrogram (31, 47). Therefore, although the use of tRNA gene PCR for phylogenetic studies of divergent species is limited, it can be a helpful tool in resolving taxonomic ambiguities for very related species with almost identical tRNA gene PCR patterns. For example, the patterns of the strains belonging to the species of the mycoides cluster closely resembled each other and were clearly different from all other species. The finding that Mycoplasma sp. bovine group 7 strain PG50 (MYP77) showed close resemblance to the M. capricolum strains is in accordance with suggestions to place these species together in an M. capricolum taxon (8, 9). Our results are also in agreement with the advice of the subcommittee on the taxonomy of Mollicutes, which favors the combination of M. mycoides subsp. capri and M. mycoides subsp. mycoides LC strains into one taxon (8) and a separate position for M. mycoides subsp. mycoides SC (small colony). Accordingly, tRNA gene PCR also supports the recent differentiation of Ureaplasma urealyticum into two distinct species (33). The clustering results indicate the close relationship between these species, but one of the intergenic tRNA spacers differs by 1 bp in length between both species.
The power of the technique was further demonstrated by the fact that tRNA gene PCR enabled detection of mixed cultures. For the mixed cultures of U. urealyticum and M. hominis, the presence of both species could be recognized because their specific patterns were included in the constructed consensus library.
In conclusion, tRNA gene PCR proved very useful for the identification of mollicute species. The unexpected high number of peaks provides the technique with a high discriminatory power. Although tRNA gene PCR is especially suited for the identification of unknown isolates, the technique can be a helpful tool to confirm current and future phylogenetic insights concerning the subdivision or merging of closely related species. Finally, we show that tRNA gene PCR can also be used to resolve mixed samples and to point to the existence of additional species.
| ACKNOWLEDGMENTS |
|---|
We kindly thank Branko Kokotovic (Danish Veterinary Institute, Copenhagen, Denmark), Lúcia Manso-Sillivan (CIRAD, Montpellier, France), and Jozef Bogaerts (Federal Agency of Health, Food Chain Security, and Environment, Brussels, Belgium) for supplying DNA of M. hyosynoviae isolates, the strains of the mycoides cluster, and the M. genitalium samples, respectively.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| Antimicrob. Agents Chemother. | Clin. Microbiol. Rev. |
|---|---|
| Clin. Vaccine Immunol. | ALL ASM JOURNALS |
|---|