Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2002, p. 5281-5287, Vol. 68, No. 11
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.11.5281-5287.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Israel Oceanographic and Limnological Research Ltd., National Center for Mariculture, Eilat,1 Department of Clinical Microbiology, The Hebrew University-Hadassah Medical School, Jerusalem 91120, Israel2
Received 24 June 2002/ Accepted 7 August 2002
|
|
|---|
|
|
|---|
Traditional diagnosis of this pathogen requires that it be recovered on culture medium and identified by means of a battery of differential biochemical tests (16). This classical method, however, often fails to identify M. marinum conclusively. Not only may the morphology and biosynthetic capabilities of mycobacteria vary depending on culture conditions (20), but also not uncommon are strain varieties that do not quite fit the typical biochemical profile for the species. In fact, the Red Sea M. marinum isolates, by being scotochromogenic and unable to hydrolyze Tween 80 and deamidate pyrazinamide, were atypical (7) and could not be positively identified to the species level until a molecular approach was adopted (17).
Over the last two decades, the 16S rRNA gene has emerged as a good standard for determining phylogenetic relationships of bacteria (35, 36). Knibb et al. (17) identified M. marinum directly from infected fish by using PCR amplification and direct sequencing of 16S rRNA products, allowing at the same time proper taxonomic assignment and opening the way to molecular epidemiologic analysis. This gene is still considered a key standard for bacterial identification due to the wealth of data available from public sequence databases (14, 34). However, as more sequence information accumulated over time, it became evident that the resolution power of 16S rRNA sequences alone is often insufficient when closely related organisms are compared (19, 27). According to Palys et al. (19), protein-encoding genes may be more discriminative than those encoding rRNA, while the analysis of two or more unlinked loci would prevent bacterial misclassification due to possible homologous recombination with other taxa.
The hsp65 gene has been highly conserved during evolution (18), is present in all known Mycobacterium species, contains species-specific variations (21) but, on average, is more variable within the genus than the gene encoding 16S rRNA (15).
In the present work, a molecular characterization of 16S rRNA and hsp65 genes was carried out with M. marinum fish isolates from various geographically distant sources to determine whether combined analysis of the two genes increases the detection of intraspecific variations. We also hoped to shed some light on the sources and epidemiology of the Israeli isolates.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Mycobacterium isolates and their hosts
|
DNA extraction.
M. marinum cells from fresh or lyophilized cultures were used. The cells from the cultures were first washed in 10 mM Tris (pH 8.0)-1 mM EDTA (pH 8.0) buffer. The resulting pellet was ground with 300 µl of grinding buffer (100 mM Tris-HCl [pH 9], 100 mM EDTA, 1% sodium dodecyl sulfate). The homogenate was incubated for 30 min at 70°C. Forty-two microliters of 8 M potassium acetate was added, and the mixture was placed on ice for 30 min and then centrifuged for 15 min at 4°C and 12 000 x g. In order to avoid traces of the pellet, the supernatant was transferred to a fresh tube and centrifuged again for 5 min. DNA was precipitated with 1 volume of isopropanol and allowed to stand for 15 min at room temperature. Pelleted DNA was washed in 70% ethanol twice, and the air-dried pellet was dissolved in 50 µl of double-distilled H2O. DNA quality was assessed by electrophoresis in 0.7% agarose and ethidium bromide staining. DNA quantity and purity were estimated with an RNA/DNA calculator (Gene Quant pro; Amersham, Cambridge, England).
PCRs.
PCRs were performed with programmable thermal controllers (PTC-100; MJ Research, San Francisco, Calif.) and a final volume of 50 µl containing 1 U of Taq DNA polymerase (Promega, Madison, Wis.), PCR buffer (Promega), 1.5 mM MgCl2, deoxynucleotide triphosphates each at a final concentration of 0.2 mM, 12.5 pmol of each primer, and 10 to 100 ng of template DNA. Typical cycling parameters were 1 min of denaturation at 94°C, 1 min of annealing at 50°C, and 1.5 min of extension at 72°C for 30 cycles. The initial denaturation step was extended to 4 min, and the final extension step was extended to 10 min.
The primers used in this study are listed in Table 2. The primers for the 16S rRNA gene were designed by comparing the 16S rRNA gene sequences of Mycobacterium species with the E. coli 16S rRNA gene sequence (GenBank accession no. J01859). Primers 246, 414, 266, and 1522 were designed by Böddinghaus et al. (2), while primers My1 and 614 were designed by us. The primers for the hsp65 gene were designed on the basis of the sequence of M. tuberculosis published by Shinnick (26). Primers TB11 and TB12 were designed by Telenti et al. (31), primer M2 was designed by Plikaytis et al. (21), and primers HSF1, HS12F, HS4F, and HSR1were designed by us.
|
View this table: [in a new window] |
TABLE 2. PCR primers
|
![]() View larger version (45K): [in a new window] |
FIG. 3. Enzyme restriction mapping. (A and B) A 614-bp PCR product from the 16S rRNA gene cleaved with restriction enzyme MslI. (A) Gel separation. (B) Restriction map. (C and D) A 796-bp product from the hsp65 gene cleaved with restriction enzyme DdeI. (C) Gel separation. (D) Restriction map. Lanes: M, 100-bp molecular size marker; 1, DL240490; 2, SR030597; 3, CC240299; 4, BB170200; 5, ATCC 927; 6, DLDK.
|
Twenty-one Mycobacterium fish isolates (and the isolate from the sea turtle) were selected for 16S rRNA (1,522 bp) sequencing, and 17 fish isolates were selected for hsp65 (1,648 bp) sequencing.
Phylogenetic analysis.
The sequences of all the isolates used in this study were compared with those of several closely related isolates available in the GenBank database (Table 3).
|
View this table: [in a new window] |
TABLE 3. Mycobacterium sp. sequences available in the GenBank database compared to the sequences of our isolates
|
![]() View larger version (32K): [in a new window] |
FIG. 1. Phylogenetic relationships among Mycobacterium species derived from the 16S rRNA gene sequences, as constructed by maximum-parsimony analysis. The tree was rooted by using N. asteroides as the outgroup. Numbers on the branches indicate bootstrap proportions (100 replicates). Bar, 1.0 nucleotide replacement. Available GenBank accession numbers are shown.
|
![]() View larger version (28K): [in a new window] |
FIG. 2. Phylogenetic relationships among Mycobacterium species derived from the hsp65 gene sequences (A) and the deduced sequences of the amino acids encoded by the Mycobacterium hsp65 gene (B), as constructed by maximum-parsimony analysis. The tree was rooted by using N. asteroides as the outgroup. Numbers on the branches indicate bootstrap proportions (100 replicates). Bar, 1.0 nucleotide or amino acid replacement. Available GenBank accession numbers are shown.
|
|
|
|---|
M. marinum marine isolates from wild Red Sea fish as well as fish farmed along both the Red Sea and the Mediterranean Sea coasts of Israel all shared an identical nucleotide sequence (group A) different from that of reference strain ATCC 927 (group C) by five bases (99.7% homology). The three Greek isolates from cultured Mediterranean Sea fish were different from each other: the first isolate (DL041200) shared an identical nucleotide sequence with the Israeli marine type, while the second isolate (DL045) was found to differ from it by only one base out of 1,522 bp (99.93% homology). Both clustered in group A. The third Greek isolate (DL049) had only 96.8% homology with the other two isolates and belonged to group F.
The two Israeli isolates from farmed freshwater fish species (koi [Cyprinus carpio] and silver perch [Bidyanus bidyanus]) shared an identical nucleotide sequence (group B), which differed by five bases from the sequence of the Israeli marine isolates in group A (99.7% homology) and by four bases from the sequence of ATCC 927 in group C (99.7% homology).
A sea bass (Dicentrarchus labrax) isolate from Denmark (DLDK), two snakehead (Channa striatus) isolates from Thailand (S4 and S267), one rainbow fish (Melanotaenia praecox) isolate from Taiwan (GenBank accession no. AF251565), one M. marinum isolate from an unspecified source in Germany (GenBank accession no. X52920), and an M. ulcerans isolate from a human infection (GenBank accession no. X88926) all shared high sequence homology to complete identity (99.8 to 100%) with ATCC 927 and clustered in group C.
Group G includes isolates from golden sturgeon (Acipenser baeri) from Holland (isolates 185408 and 185409), which exhibit 100% identity with Mycobacterium sp. strain Fuerth 1999 (GenBank accession no. AF152558), a clinical isolate from a human infection from an unspecified location.
hsp65 phylogenetic analysis.
The tree constructed from the hsp65 gene sequences (Fig. 2A) revealed 10 distinct groups. The marine and freshwater Israeli type isolates (groups A and B) showed the same clustering pattern as in the 16S rRNA tree and differed from each other by four bases (99.7% homology). They differed from strain ATCC 927 by nine bases (99.4% homology) and seven bases (99.6% homology), respectively.
The three Greek isolates showed the same clustering pattern as well: isolates DL041200 and DL045 shared an identical nucleotide sequence with the marine Israeli type isolates and clustered in group A, while the third isolate (DL049) belonged to group F.
In contrast to what was seen with the tree constructed from the 16S rRNA gene sequences, group C was further divided into three subgroups, which included M. ulcerans (subgroup CI) (99.4% homology with ATCC 927), the isolate from sea bass from Denmark (isolate DLDK) and ATCC 927 (subgroup CII) (100% homology), and the isolates from snakehead from Thailand (subgroup CIII) (99.6% homology with ATCC 927). The isolates from golden sturgeon from Holland (isolates 185408 and 185409) clustered in group G.
A comparison of the deduced sequences of the amino acids encoded by the Mycobacterium hsp65 gene is shown in Fig. 2B. All Israeli isolates, marine and freshwater, and two of the Greek isolates clustered in one group (group A). M. ulcerans belonged to group B, while the isolate from sea bass from Denmark (isolate DLDK), snakehead from Thailand (isolates S4 and S267), and ATCC 927 clustered together in group C. The isolates from golden sturgeon from Holland (isolates 185408 and 185409) clustered in group F, while the Greek isolate from sea bass (isolate DL049) belonged to group G.
The sequences of the Israeli type isolates (marine and freshwater) differed from the sequences of all published Mycobacterium spp. in both genes. Consequently, a new name, M. marinum strain Eilaticum DL240490, was proposed for the marine isolates (GenBank accession no. AF456238 for the 16S rRNA gene sequence and AF456468 for the hsp65 gene sequence). The Israeli freshwater isolates were named M. marinum strain Cyprinum CC240299 (GenBank accession no. AF456239 for the 16S rRNA gene sequence and AF456469 for the hsp65 gene sequence).
Restriction enzyme mapping.
Based on the sequence information obtained from both genes, a two-step restriction map was drawn (Fig. 3B and D). MslI digestion of a 614-bp PCR product from the 16S rRNA gene yielded cleavage products of 80 and 534 bp, which separated the Israeli type isolates from the other selected isolates, as they have a unique site for the enzyme (Fig. 3A, lanes 1 to 4). This group was further characterized by DdeI digestion of a 796-bp PCR product from the hsp65 gene, which yielded cleavage products of 182 and 614 bp. Having a unique site for this enzyme, the freshwater isolates could be differentiated from the marine isolates (Fig. 3C).
|
|
|---|
Direct sequencing of the 16S rRNA and hsp65 genes confirmed the identities of all the Israeli isolates that were characterized in this study as M. marinum. Phylogenetic analysis based on the 16S rRNA and hsp65 gene sequences clustered the Israeli type Mycobacterium isolates, discriminating this type from all other Mycobacterium isolates found in the literature. Site-specific nucleotide differences revealed that in the last decade, all M. marinum isolates in Israel have belonged to two strains, one that has been infecting marine fish (both cultured and wild) and one that has been infecting freshwater fish (cultured).
Bacterial strains have been considered to be within the same species if they have fewer than 5 to 15 base differences in the 16S rRNA gene (10). For mycobacteria, a difference of less than 5 nucleotides within the complete 16S rRNA gene, along with clear phenotypic differences, is supposed to indicate a genetically unique and distinct taxon (34). Accordingly, M. marinum and members of the M. tuberculosis complex lie on close branches (25, 33). However, M. marinum shares over 99.8% homology with M. ulcerans, with only two nucleotide differences (22, 33). In fact, Tønjum et al. (33) pointed out that the 16S rRNA gene is not useful for discriminating among M. marinum, M. ulcerans, and M. haemophilum. However, the diseases caused by M. marinum and M. ulcerans in humans differ greatly in their clinical, histopathologic, and epidemiologic aspects (3, 29).
The sequences of the hsp65 gene discriminated among M. ulcerans, M. marinum ATCC 927 (reference strain), and M. marinum from the Far East (Thailand; group C; GenBank accession no. AF456474), which could not be distinguished with the 16S rRNA gene sequences. Whole-genome techniques, such as IS2404 restriction fragment length polymorphism analysis and amplified fragment length polymorphism analysis, clustered M. ulcerans and M. marinum into separate groups as well (5).
In this study, the restriction enzyme map based on the combination of the two gene (16S rRNA and hsp65) sequences enabled us to differentiate not only the Israeli type isolates from all other M. marinum isolates but also the Israeli marine isolates from the Israeli freshwater isolates. PCR-restriction fragment length polymorphism analysis of the hsp65 gene (31), although an established discriminating tool used with mycobacteria by many laboratories (4, 8, 24), is not applicable to the M. marinum complex, as different strains have the same restriction sites. This problem was thus overcome with our method.
Our sequence of the 16S rRNA gene of M. marinum ATCC 927 (reference strain) (GenBank accession no. AF456240) is identical to that published by Rogall et al. (25) (GenBank accession no. X52920) and is in contrast (93 to 98% homology) to the partial sequence published by Talaat et al. (30) (GenBank accession no. U92088). Also, the 16S rRNA gene sequence of the isolates from golden sturgeon (A. baeri) from Holland (isolates 185408 and 185409) exhibited 100% identity with that of Mycobacterium sp. strain Fuerth 1999 (GenBank accession no. AF152558), whose sequence was in turn found by Turenne et al. (34) to be identical to that of M. chelonae ATCC 19237.
Because there were no records of M. marinum in the Red Sea before the report by Colorni (7) for sea bass (D. labrax) and because of the active fish trade between Israel and other Mediterranean countries in recent years (in particular, for fingerlings hatched in one region to be reared in another), it had been assumed that the pathogen was introduced into Israel via mariculture operations and had then spread to other farmed and native species in the Gulf of Eilat (A. Diamant and A. Colorni, Abstr. Proc. 7th Int. Conf. Eur. Assoc. Fish Pathol., p. 66, 1995). This hypothesis was supported by the facts that acid-fast bacteria, some of which were specifically identified as M. marinum, had already been detected in D. labrax in Europe (12; W. Verdonck, L. Lambrechts, and F. Ollevier, Abstr. Proc. 2nd Int. Colloq. Pathol. Mar. Aquacult., p. 149, 1986; S. Mellergaard and I. Dalsgaard, Abstr. 6th Int. Conf. Eur. Assoc. Fish Pathol., p. 60, 1993) and that disease in Israel was first detected in offspring of D. labrax which originated in France. While our results show that no linkage exists between the Israeli type and the northern European M. marinum isolates, the hypothesis that Israeli marine M. marinum was nevertheless imported from Europe received support by the identity of one of the Greek isolates to the Israeli marine isolates. On the other hand, the fact that in the last decade, the same strain of M. marinum, different from the freshwater isolates, infected all the Israeli marine fish species as well as a captive hawksbill sea turtle in the Red Sea may indicate that a local, Red Sea strain of M. marinum exists and found extremely susceptible hosts in the exotic species (mainly sea bass) imported for aquaculture. In fact, our results suggest that every isolate of M. marinum examined is endemic, with a unique genotype specific and circumscribed to its geographic region. Interestingly, Portaels et al. (23) and Chemlal et al. (6) have similarly shown for M. ulcerans that the polymorphisms in the 3' end of the 16S rRNA gene sequence and in the IS2404 restriction fragment length polymorphism pattern, respectively, are related to their geographic origins (Australia, Africa, and Central America). Further characterization of M. marinum isolates from additional geographic regions may clarify whether this is indeed the case.
Present address: Bribie Island Aquaculture Research Centre, Bribie Island, Queensland 4507, Australia. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»