Previous Article | Next Article ![]()
Journal of Clinical Microbiology, May 2003, p. 2015-2026, Vol. 41, No. 5
0095-1137/03/$08.00+0 DOI: 10.1128/JCM.41.5.2015-2026.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
Food Animal Health Research Program, Ohio Agricultural Research and Development Center, The Ohio State University, Wooster,1 Animal Disease Diagnostic Laboratory, Ohio Department of Agriculture, Reynoldsburg,4 Department of Veterinary Preventive Medicine, The Ohio State University, Columbus, Ohio,2 National Animal Disease Center, Agricultural Research Service, U.S. Department of Agriculture, Ames, Iowa,6 Department of Molecular Microbiology, University of Central Florida, Orlando, Florida,5 Biomedical Genomics Center and Departments of Microbiology and Veterinary Pathobiology, University of Minnesota, St. Paul, Minnesota3
Received 8 November 2002/ Returned for modification 8 January 2003/ Accepted 31 January 2003
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Comprehensive analysis of the molecular diversity and comparative molecular pathology of M. avium subsp. paratuberculosis will help establish the degree of heterogeneity in strains isolated from a variety of host species. The extent of strain sharing across a variety of hosts will reflect the degree of interspecies transmission. It will also identify any strain similarities between JD and CD. Until recently, an IS900-based restriction fragment length polymorphism (RFLP) fingerprinting has been applied to study molecular epidemiology of M. avium subsp. paratuberculosis isolates. Extensive analyses of the IS900-RFLP patterns have identified that JD in cattle and other species such as goats and rabbits (3, 20, 33) is caused by indistinguishable strains. Two M. avium subsp. paratuberculosis isolates from a CD patient were also shown to carry a bovine-like IS900-RFLP fingerprint (47). JD in sheep appears to be caused by a different M. avium subsp. paratuberculosis strain (3, 47). However, the occurrence of a "sheep strain" in cattle with JD has been reported, indicating that interspecies transmission cannot be ruled out (49). Use of other independent fingerprinting and molecular diversity analysis tools to verify strain sharing is essential to provide substantial evidence of interspecies transmission.
The control and eradication of JD is severely hindered by prolonged incubation time, the presence of undetected subclinical cases, the absence of M. avium subsp. paratuberculosis-specific diagnostic tools, and the lack of knowledge of strain diversity. Current diagnostic methods include isolation from fecal and tissue specimens, enzyme-linked immunosorbent assay (ELISA), and IS900-based PCR. Although culturing the organisms is considered a "gold standard," this method is fraught with difficulties. For example, the organism takes at least 12 to 16 weeks to grow to detectable levels, and even the most sensitive culture methods have only 50% sensitivity (43). Serologic tests such as agar gel immunodiffusion (18), complement fixation (25), and ELISA are limited in their use because of low specificity and sensitivity since antibodies may not be detectable either due to anergy or until their late appearance in the pathogenesis of JD (32, 42).
The conventional wisdom is that M. avium subsp. paratuberculosis differs from other members of the M. avium complex (MAC) in having 14 to 18 copies of IS900 inserted into conserved loci in its genome. Insertion element IS900 is a 1,451-bp segment that lacks inverted terminal repeats and has features characteristic of the IS900 family such as homologous transposases and particular insertion sites (19). IS900-based PCR identification techniques have been routinely used for M. avium subsp. paratuberculosis diagnosis (22, 50), and a number of studies have used DNA probes based on IS900 for M. avium subsp. paratuberculosis detection in tissues (29, 50), fecal specimens (17, 43), and milk samples (35). There is growing concern about the zoonotic potential of M. avium subsp. paratuberculosis, and several studies have attempted to show a causal relationship between M. avium subsp. paratuberculosis and CD by demonstration of IS900 in tissue (22, 37) and milk (9, 31, 35) samples of CD patients. Additionally, based on the presence of IS900, several studies have implied that bovine milk may serve as a mode of transmission of M. avium subsp. paratuberculosis (9, 35). However, IS900-like elements have been found in MAC strains (30) and in some atypical mycobacteria (10, 14, 28). This implies that diagnosis of M. avium subsp. paratuberculosis based on detection of IS900 alone needs to be evaluated with caution. Hence, there is a need for identification and large-scale analysis of additional M. avium subsp. paratuberculosis-specific molecular targets other than IS900 to confirm diagnosis. We report here a molecular test that exploits the fact that IS900 is integrated into conserved loci in the genome and that these integration sites are unique to M. avium subsp. paratuberculosis. Microtiter plate hybridization to integration site-specific probes was used to confirm specificity of PCR amplification (40). A recently identified M. avium subsp. paratuberculosis unique sequence 251 (2) was also evaluated as a target sequence for use in molecular diagnosis. The presence of hsp65 polymorphism was used to confirm the identification of M. avium subsp. paratuberculosis (15, 26). Two different molecular typing methods, multiplex PCR of 28 IS900 integration loci (MPIL) (5) and amplified fragment length polymorphism (AFLP) analysis were used to fingerprint M. avium subsp. paratuberculosis isolates obtained from a variety of geographical localities and host species. An IS1311 PCR and restriction endonuclease analysis (PCR-REA) was used to differentiate between bovine and ovine M. avium subsp. paratuberculosis isolates (46) within MPIL and AFLP clusters.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
|
Genetic fingerprinting by MPIL. Genetic fingerprints of a subset of the M. avium subsp. paratuberculosis and MAC isolates were generated by using a method described previously (5) with several modifications. Briefly, six sets of multiplex PCRs referred as 9R1, 5R2, 4L1, 4L2, 3L3, and 3L4 were carried out as described earlier. The modifications to the published technique included the use of primer Loc14L in set 5R2 (instead of primer Loc14R) and primer Loc14R (6) in set 3L3 (instead of primer Loc14L). All primers were used at a concentration of 0.2 µM except the common primer, which was used at a concentration of 1 µM in a reaction volume of 50 µl that contained 5 U of Taq polymerase, 1x PCR buffer containing 1.5 mM MgCl2 (Qiagen), a 200 µM concentration of deoxyribonucleotide triphosphates (Amersham Pharmacia Biotech, Inc., Piscataway, N.J.), 10% dimethyl sulfoxide (DMSO), and 1 µl (ca. 50 ng/µl) of genomic DNA. The thermal cycler parameters were as described previously (5). The PCR products were electrophoresed at 100 V for 4 h in a 2% agarose gel prestained with ethidium bromide and visualized on a UV transilluminator (Alpha Innotech Corporation, San Leandro, Calif.).
Analysis of MPIL fingerprints. The resulting fingerprints were converted into binary data where "1" indicated the presence and "2" indicated the absence of a specific band. Cluster analysis was done by using Molecular Evolutionary Genetics Analysis (MEGA; version 2.1) by the unweighted pair group method with averages (UPGMA) (39). The distance matrix for input into the MEGA program was created from the binary data by using ETDIV and ETMEGA (http://foodsafe.msu.edu/whittam/programs/).
AFLP of M. avium subsp. paratuberculosis. AFLP was performed as described in the AFLP Microbial Fingerprinting Kit protocol (Perkin-Elmer Applied Biosystems Division, Foster City, Calif.). In brief, the M. avium subsp. paratuberculosis DNA was restricted with endonucleases EcoRI and MseI. The restricted fragments were then ligated to EcoRI and MseI adapters. The restriction-ligation reaction mixture included 100 ng of genomic DNA, 5 U of EcoRI (New England Biolabs, Inc. [NEB], Beverly, Mass.), 5 U of MseI (NEB), 1 U of DNA ligase (NEB), and a 0.4 µM concentration of adapter pairs in T4 DNA ligase buffer (50 mM Tris-HCl [pH 7.8], 10 mM dithiothreitol, 1 mM ATP, 25 µg of bovine serum albumin). Preselective amplification was carried out in a final reaction volume of 20 µl that contained 4-µl portions of the products of the restriction-ligation reaction, 125 pmol of AFLP preselective primer pairs (EcoRI preselective primer [5'-GACTGCGTACCAATTC] and MseI preselective primer [5'-GATGAGTCCTGAGTAA]), and 15 µl of AFLP core mix (PE-ABD). PCRs were performed in a thermal cycler (GeneAmp PCR system 9700; Perkin-Elmer) with an initial denaturation at 72°C for 2 min, followed by 20 PCR cycles of denaturation at 94°C for 1 s, annealing at 56°C for 30 s, an extension step at 72°C for 2 min, and a final extension step at 60°C for 30 min. Selective amplification was performed with the EcoRI selective primer containing a fluorescent label on the 5' end with additional base C (carboxyfluorescein dye, FAM-EcoRI-C [FAM-5'-GACTGCGTACCAATTCC]) and MseI selective primer with no additional base (MseI-O [5'-GATGAGTCCTGAGTAA]) in a final reaction volume of 10 µl containing 1.5 µl of the preselective amplification products, 125 pmol of MseI-O primer, 25 pmol of FAM-EcoRI-C primer, and 7.5 µl of AFLP core mix (PE-ABD) in a thermal cycler (GeneAmp PCR system 9700). Aliquots (1 µl) of each amplification product were mixed with 2 µl of sterile water and subjected to electrophoresis in 5% Long Ranger gels (LMC Bioproducts, Rockland, Maine) by using the ABI model 377 automated DNA fragment analyzer (Perkin-Elmer). An internal size standard (GS-500, ROX size standard) was included in each lane. The AFLP fingerprint profiles were then automatically analyzed with GeneScan (Perkin-Elmer). Sizes of amplified products were then tabulated and exported for cluster analysis by using Molecular Analyst software (Bio-Rad, Hercules, Calif.). To analyze AFLP fingerprints, DNA fragments of from 50 to 500 bp were included for cluster analysis.
Computer assisted analysis of AFLP fingerprints. Molecular Analyst software (Bio-Rad) was used to compare the AFLP fingerprint profiles among M. avium subsp. paratuberculosis isolates. The program automatically computed the similarity for each pair of fingerprints on the basis of band positions by using the Dice coefficient (SD). Pairs of isolates with a similarity coefficient (SD) of 0.95 were considered similar and not distinguishable. Cluster analysis among isolates was performed by the UPGMA linkage method.
Concordance analysis. Kappa statistics (27) were calculated for each concordant pair of the isolates (n = 100) that were analyzed by both MPIL and AFLP techniques to assess the degree of agreement between the two fingerprinting techniques.
PCR amplification and hybridization for IS900 integration loci L1 and L9. Two IS900 integration sites designated locus L1 (left-side integration site of IS900 into an unknown open reading frame [ORF]) and locus L9 (left-side integration of IS900 into alkA) were amplified by using PCR. The primers used were L1 (5'-CCCGTGACAAGGCCGAAGA-3'), L9 (5'-CGGCCCTGGCGTTCCTATG-3'), and the biotinylated common forward primer 900R (5'-/5Bio/ACGCTGTCTACCACCCCGCA-3') (5). Primers were used at a final concentration of 0.2 µM in a 50-µl PCR master mix also including 1.25 U of Taq polymerase, 1x PCR buffer containing 1.5 mM MgCl2 (Qiagen), and a 200 µM concentration of deoxyribonucleotide triphosphates (Amersham Pharmacia Biotech). A total of 5 µl (ca. 50 ng/µl) of genomic DNA was used in each reaction. The thermal cycler parameters used for amplification included an initial 95°C incubation step for 15 min for Taq polymerase activation, followed by 35 cycles of denaturation at 94°C for 15 s, annealing at 58°C for 20 s, and extension at 72°C for 20 s, with a final extension step at 72°C for 7 min. A PCR blank was included for each batch and used as negative control for molecular analyses. The amplification products were detected in two independent ways. First, they were visualized on a UV transilluminator (Alpha Innotech) after electrophoresis at 125 V for 45 min in 1.5% agarose gels prestained with 5% ethidium bromide. Second, the biotinylated amplicons were detected in a reverse hybridization reaction on a microtiter plate as described previously (40). The probes were designed to target the unique integration regions (spanning IS900 into the genes of insertion) located within the amplicons. The probes used were PrL1 (spanning the IS900-unknown ORF region; 5'-TACACACAGCCGCCATACACTTCGCTTCATGCCCTTACG-3') and PrL9 (spanning the IS900-alkA region; 5'-AGCCATCGAAGGGGATTCTCAATGACGTTGTCAA-3'). The two integration loci were detected simultaneously in microtiter wells (Corning Incorporated, Corning, N.Y.) coated with a mixture of both probes as described previously (40). DNA extraction blanks, PCR amplification negatives, and hybridization negatives were used in all analyses as negative controls. The cutoff optical density was set at 0.2 on the basis of analysis of several known negatives and M. avium subsp. paratuberculosis positives.
PCR amplification of M. avium subsp. paratuberculosis unique sequence 251. The primers used to amplify the entire 417-bp ORF of M. avium subsp. paratuberculosis unique sequence 251 have been published elsewhere (2). PCRs were performed in 25-µl volumes, which included 0.625 U of Taq polymerase, 1x PCR buffer, 2.5 mM MgCl2, a 200 µM concentration of deoxyribonucleotide triphosphates (Amersham Pharmacia Biotech), a 0.2 µM concentration of each primer, 5% DMSO, and 2 µl (ca. 50 ng/µl) of genomic DNA. The thermal cycler parameters used for amplification included an initial 95°C incubation step for 15 min for Taq polymerase activation, followed by 35 cycles of denaturation at 94°C for 20 s, annealing at 53°C for 30 s, and extension at 72°C for 30 s, with a final extension step at 72°C for 7 min. Amplicons were visualized after electrophoresis at 125 V for 45 min in a 1.5% agarose gel prestained with ethidium bromide.
PCR-REA for polymorphisms in hsp65 and IS1311 genes. New primers were designed to target previously described (15, 46) polymorphic regions of hsp65 and IS1311. A 248-bp fragment of hsp65 was amplified with the primers hsp3 (5'-GCCGCTGCTGATCATCGCCGA--3') and hsp4 (5'-CCTTGGTGACGACGACCTT--3'), whereas a 749-bp region of IS1311 was amplified with the primers IS1311F (5'-GCGTGAGGCTCTGTGGTGAA--3') and IS1311R (5'-TCAGAGATCACCAGCTGCAC--3'). PCR was carried out in a reaction volume of 50 µl. The master mix included 1.25 U of Taq polymerase, a 0.2 µM concentration of each primer, 1x PCR buffer containing 1.5 mM MgCl2 (Qiagen), a 200 µM concentration of deoxyribonucleotide triphosphates (Amersham Pharmacia Biotech), 5% DMSO, and 2 µl (ca. 50 ng/µl) of genomic DNA. The thermal cycler parameters used were as described for integration loci PCR. A PCR blank was included in each batch and used as negative control for molecular analyses. After the amplification, 10 µl of each PCR product was digested with 4 U of restriction enzyme at 37°C for 1 h according to the supplier's instructions. PstI (New England Biolabs) was used to detect polymorphism in nucleotide 861 of the hsp65 ORF (GenBank accession no. X74518, ORF from 125 to 1750), and HinfI (New England Biolabs) was used to detect the IS1311 polymorphism in nucleotide 167 of the IS1311 ORF (GenBank accession no. AJ223974, ORF from 35 to 1219). Both digested and undigested amplicons were electrophoresed in adjacent wells of 1.5% agarose gels and visualized on a UV transilluminator (Alpha Innotech).
DNA sequencing. PCR products were sequenced by using standard dye terminator chemistry and analyzed on an automated DNA sequencer (Perkin-Elmer ABI 377). Sequences were aligned against sequences in GenBank database and analyzed by using Editseq and MegAlign programs (DNASTAR, Inc., Madison, Wis.).
| RESULTS |
|---|
|
|
|---|
|
|
|
= 0.2) and A18-Z2 (
= 0.25).
PCR-REA for polymorphisms in IS1311 gene.
A C
T polymorphism in nucleotide 167 in ORF (GenBank accession no. AJ223974, ORF from 35 to 1219) of insertion element IS1311 is unique to bovine M. avium subsp. paratuberculosis isolates, whereas ovine M. avium subsp. paratuberculosis strains are identical to M. avium subsp. avium in this locus (46). This results in the gain of an HinfI restriction site in M. avium subsp. paratuberculosis of bovine origin. A panel of 190 isolates (Table 1) was analyzed for polymorphism in the IS1311 gene. Undigested amplicons from all MAC isolates and an atypical Mycobacterium (M. xenopi) showed the presence of a single 749-bp band after PCR. After digestion with HinfI all cattle, goat, mouse, and deer M. avium subsp. paratuberculosis isolates and five MAC isolates showed presence of the predicted four bands (67, 85, 218, and 379 bp), while all ovine M. avium subsp. paratuberculosis isolates and 50 of 55 MAC isolates showed the presence of three bands (85, 285, and 379 bp), indicating the absence of the restriction site (Fig. 5). The other five atypical mycobacterium and nonmycobacterium isolates did not amplify the IS1311 segment. Amplification specificity was confirmed by sequencing a bovine, an ovine, and a MAC isolate. All three isolates showed homology to the published sequences (GenBank accession numbers AJ308375, AJ223975, and U16276), as well as the expected polymorphism. One bovine isolate, which was L1-L9 and 251 positive, showed a different pattern when the restriction enzyme products were electrophoresed (Fig. 5). The undigested product was of the expected band length (749 bp). The nucleotide sequence of this product revealed that this isolate had numerous mutations resulting in a total of six restriction sites, only one of which was in the expected position. Restriction digestion resulted in seven fragments (31, 63, 85, 95, 130, 148, and 195 bp) instead of the expected three or four.
|
Presence of M. avium subsp. paratuberculosis-specific target 251 in known M. avium subsp. paratuberculosis isolates. A subset of 328 isolates from the culture collection was analyzed for the presence of M. avium subsp. paratuberculosis unique sequence 251 (Table 1). Agarose gel electrophoresis showed the presence of a 417-bp band for all 113 known M. avium subsp. paratuberculosis isolates and 5 MAC isolates that were also positive by L1-L9 PCR. The locus was absent in all other mycobacterial (n = 27), atypical mycobacterial (n = 6), and nonmycobacterial bacterial (n = 121) isolates analyzed. The PCR products from nine bovine (Ohio), three human (CD), and three ovine (Ohio) isolates were sequenced. The nucleotide sequences had 100% identity with the published ORF of locus 251 (GenBank accession no. AF445445, ORF from 123 to 539), indicating the stability of this target and its usefulness in JD diagnostics.
PCR-REA for polymorphisms in hsp65 gene.
A subset of 272 isolates (Table 1) that were analyzed for the presence of L1 and L9 integration sites were reflexed for identity confirmation by PCR-REA analysis of this locus. A G
T polymorphism in nucleotide 861 in the ORF (GenBank accession no. X74518) of the hsp65 gene is unique to M. avium subsp. paratuberculosis and leads to the loss of a PstI restriction site in the species. Prior to digestion, all mycobacterial isolates analyzed showed the presence of a single 248-bp band when visualized by agarose gel electrophoresis. After digestion with PstI, all M. avium subsp. paratuberculosis isolates retained the full 248-bp band, while all MAC and other bacterial isolates showed the presence of two bands (135 and 113 bp), indicating the presence of the restriction site. All six atypical mycobacterium and 10 M. tuberculosis complex isolates showed a MAC-like pattern. The nonmycobacterial isolates included in the analysis did not amplify the region of interest.
| DISCUSSION |
|---|
|
|
|---|
Molecular diversity in M. avium subsp. paratuberculosis. Mutations in structural genes (24) or antigenic variations (16) can serve as useful molecular markers to study epidemiology and natural history of an organism. However, in well-conserved genomes such as that of mycobacteria (41), these markers provide limited information in strain typing. In mycobacterial genomes, the degree of diversity in the number and site of integration of transposable genetic elements or single nucleotide polymorphisms in genes associated with drug resistance or host immunological pressure may be more useful than those that rely on genome-wide single nucleotide polymorphisms as indicators of strain variation. Several molecular epidemiological studies have demonstrated the usefulness of insertion sequences in determining strain distribution of mycobacteria such as IS1245 in M. avium (45) and IS6110 in M. tuberculosis (44) and IS711 in Brucella sp. (4). Other repetitive elements have also been used in strain typing (12). Seminal studies that identified IS900 in M. avium subsp. paratuberculosis strains suggested that it was exclusively present in M. avium subsp. paratuberculosis (19). Hence, IS900 has been the marker of choice for most fingerprinting studies reported (11, 34, 47). This exploits the fact that the IS900 elements show a high degree of target sequence specificity resulting in similar fingerprints in epidemiologically related isolates (19). IS900-RFLP analyses have broadly divided the M. avium subsp. paratuberculosis strains into either sheep or cattle and other ruminant subtypes, designated S and C, respectively (3, 8, 47). In contrast, isolation of an S strain from cattle has been reported (49). Most investigations have identified limited strain sharing between cattle and sheep hosts, whereas little variation in fingerprints is evident between cattle and other animals, including rabbits and goats, and M. avium subsp. paratuberculosis isolates from a human with CD. Evidence of strain sharing between bovine and human M. avium subsp. paratuberculosis isolates is of special interest since this implies the existence of a potential animal reservoir for CD. Several studies have implicated an inability to compare fingerprints from sheep isolates due to difficulties encountered in culturing them. Thus, an indication of limited strain sharing between bovine and ovine isolates does not rule out interspecies transmission since there may be a bias in the analysis due to under-representation of the slow-growing ovine strains. In summary, there has been inadequate data to support the idea of extensive strain sharing and interspecies transmission. Besides being elaborate and expensive, fingerprinting by RFLP requires large amounts of unsheared genomic DNA extracted from relatively large cultures. There are multiple reports (10, 14, 28, 30) of the presence of IS900-like elements in mycobacteria other than M. avium subsp. paratuberculosis. Sequence analysis of a few MAC isolates with the commonly used IS900 primers (3) showed the presence of IS900 elements (data not shown). Nucleotide sequences indicated significant homology between these regions and the IS900 element in M. avium subsp. paratuberculosis. This supports the notion that the use of IS900 probe as a tool in molecular typing could result in false positives corresponding to IS900-related elements in mycobacteria other than M. avium subsp. paratuberculosis (8, 38).
Both techniques reported in the present study require very small quantities of genomic DNA and are applicable to degraded DNA samples not suitable for RFLP analysis. AFLP analysis results with 104 isolates from distinct geographic regions and hosts showed that 72% of the M. avium subsp. paratuberculosis, including bovine, caprine, and murine strains, fell into either one of the two major clusters reiterating the low degree of heterogeneity at the genomic level. This suggests that only a few strains may be responsible for disease in bovine hosts and that interspecies transmission could have occurred since these clusters also included M. avium subsp. paratuberculosis strains isolated from other hosts. In agreement with previous studies, the two major branches included only one of the four ovine isolates, while the other two ovine strains clustered into distinct branches in close proximity to a bovine isolate. In contrast to earlier reports (33, 47) and our MPIL data, AFLP fingerprints of the human M. avium subsp. paratuberculosis isolates were unique and did not cluster with either the bovine or ovine strains, suggesting that they were relatively more heterogeneous at the genome level.
In MPIL analysis 78% of the U.S. bovine isolates, including those from Argentina and Nova Scotia, fell within one major cluster, indicating a significant degree of uniformity within isolates infecting bovine hosts. M. avium subsp. paratuberculosis strains from other ruminants, such as goat and deer, and one mouse strain also fell within this cluster. Similar to AFLP, this major cluster included only 2 of the 16 ovine isolates analyzed. The analysis identified clear genetic diversity between ovine isolates, whereas there were limited differences within strains from bovine hosts across several geographic localities. M. avium subsp. paratuberculosis isolates from the human CD patients were more diverse with phylogenetic proximities at various levels to strains isolated from both bovine and ovine host species, an indication of a greater level of diversity. Although this suggests a close association of the human and animal strains, it does not provide direct evidence for a causal role of M. avium subsp. paratuberculosis in CD. Comprehensive comparative genomics with larger numbers of human isolates and reproduction of JD in animals by the human strains will be required to demonstrate, albeit indirectly, a causal association between M. avium subsp. paratuberculosis and CD.
The possibility of nonspecific amplification of Loc5R and Loc5L primers in M. avium subsp. avium strains due to homology in desA1 and desA2 genes have been described (5). Similar nonspecific amplifications may have been the case in the low-copy strains of our collection. Conversion of this gel-based analysis into a microarray by using integration site probes would provide a more accurate picture of the fingerprinting by MPIL.
Concordance analysis showed that AFLP fingerprinting was able to discriminate between isolates clustered together by MPIL analysis. Similarly, MPIL fingerprinting discriminated between isolates in one cluster by AFLP analysis. Kappa scores for the major concordant pairs show that fingerprints generated by MPIL and AFLP analysis are independent. This indicates the need to use more than one fingerprinting technique to distinguish between the M. avium subsp. paratuberculosis strains. Restricted allelic variation within this mycobacterial subspecies may limit the application of any single fingerprinting method in epidemiological studies of JD. Alternately, the use of other combinations of restriction enzymes (such as ApaI and TaqI) or application of a multiplex-AFLP protocol (23) may provide a better indication of diversity.
Analysis of the isolates for presence of polymorphism in IS1311 by PCR-REA facilitated the classification of the M. avium subsp. paratuberculosis isolates into bovine or ovine strains. The gain of the HinfI restriction site due to a single nucleotide polymorphism, which is unique to the bovine strain, is indicative of limited strain sharing between the two hosts. One bovine isolate from an Ohio farm showed a difference in the density of the restricted bands compared to other isolates in our analysis. This is consistent with the speculation (46) that M. avium subsp. paratuberculosis isolates may vary in the number of IS1311 insertion elements carrying the polymorphism. The well-conserved nature of the mycobacterial genome indicates that the probability of several mutations within one gene is very low. Hence, the bovine M. avium subsp. paratuberculosis isolate with multiple HinfI restriction sites within IS1311 is of special interest. Although the absence of the C
T polymorphism indicates that it is an ovine strain, it clustered with other bovine isolates by MPIL analysis.
Integration loci as M. avium subsp. paratuberculosis unique targets. IS900 is a 1,451-bp segment that lacks inverted repeats and has features characteristic of the IS900 family such as homologous transposases and particular insertion sites (19). Detection of a portion of IS900 sequence is routinely used to detect the presence of M. avium subsp. paratuberculosis in field and clinical samples (37, 50). However, the presence of IS900-like sequences has been demonstrated in non-M. avium subsp. paratuberculosis mycobacteria, including M. avium subsp. avium (30, 38), M. cookii (14), M. marinum, and M. scrofulaceum (10, 28). Because of the importance of differentiating these closely related mycobacteria in clinical specimens there is a need for molecular targets based on genes other than IS900 to confirm the presence of M. avium subsp. paratuberculosis.
The present study reports the development of a multiplex PCR to amplify two distinct integration sites of the insertion element IS900. This exploits the fact that the IS900 elements integrate into conserved loci in the M. avium subsp. paratuberculosis genome (19). Although geographically restricted exceptions cannot be ruled out, this large-scale analysis indicates that the integration sites of IS900 at loci L1 and L9 were unique to M. avium subsp. paratuberculosis and are consistently present in M. avium subsp. paratuberculosis from diverse locations and hosts. These two integration sequences were absent in non-M. avium subsp. paratuberculosis mycobacteria, nonspecific bacterial DNA, and in 80 of 85 MAC strains analyzed, 33 of which had low numbers (one to six) of IS900 integration sites as determined by multiplex PCR. The five MAC isolates that showed presence of both the integration sites were reclassified as M. avium subsp. paratuberculosis on the basis of M. avium subsp. paratuberculosis-like pattern for six distinct molecular targets: L1-L9, 251, hsp65 PCR-REA, IS1311 PCR-REA, MPIL, and AFLP (two isolates) cluster analysis. The isolates were confirmed to be mycobactin J dependent. This implies that these sites could serve as valuable tools in the molecular diagnostics of M. avium subsp. paratuberculosis and can significantly reduce the time to diagnosis of M. avium subsp. paratuberculosis from 16 to 23 weeks to 2 to 10 days in subclinical infections. Routine use of the integration sites to detect M. avium subsp. paratuberculosis will prevent the reporting of false-positive IS900 PCR results. A recent study in our laboratory of feces from field animals (A. Ozbek, F. C. Michel, M. Strother, A. S. Motiwala, B. Byrum, W. Shulaw, C. G. Thornton, and S. Sreevatsan, unpublished data) established the feasibility of application of this method in the development of a sensitive assay for primary clinical samples.
To confirm the stability of the integration loci as diagnostic targets and to evaluate the possibility of reflexing ambiguous isolates to hsp65 PCR-REA, we analyzed a subset of bacterial isolates for polymorphism in hsp65 gene. This subset included M. avium subsp. paratuberculosis strains obtained from a variety of host species. Although this technique proved to be very valuable for pure cultures, subsequent analysis of M. avium subsp. paratuberculosis from broth cultures showed that it is not as efficient for mixed cultures from feces or broth due to the highly conserved nature of this gene across many species. Hence, hsp65 analysis was restricted to pure cultures, whereas locus 251 was used as an additional marker for broth cultures.
Recently identified unique M. avium subsp. paratuberculosis sequences (2) show promise of serving as target sequences in molecular diagnosis of M. avium subsp. paratuberculosis. One of these sequences, locus 251, was consistently detected with the L1-L9 integration sites in 118 of 121 known positive isolates. Sequence analysis from bovine, ovine, and human isolates show 100% homology with the published sequence, indicating the possibility of applying this locus in a multiplex PCR-hybridization format, along with LI-L9 integration sites to aid in the unambiguous diagnosis of JD. Although the specificity in amplification of the locus from M. avium subsp. paratuberculosis isolates in our sample was absolute, we note that the forward primer (251R) used for the amplification of locus 251 revealed homology (15 of 20 nucleotides at the 3' end) to plasmid sequences from the enteric bacterium Salmonella spp. on BLASTN searches with the public databases (GenBank accession numbers AF250878 and AL513383). Hence, positive results of amplification, as seen on a gel, must be interpreted with caution, especially when fecal samples and primer pairs are used as described in the present study. Postamplification detection on a solid phase (such as a microtiter plate) through hybridization to an internal sequence is likely to improve the specificity of the results.
In conclusion, the present study identified a high degree of genetic similarity within the bovine isolates, indicating that only a few very closely related clones of M. avium subsp. paratuberculosis may be responsible for widespread infection in cattle, other ruminants, and possibly wildlife. This study established strain sharing among isolates from a variety of hosts and geographic locations. In addition, we identified multiple M. avium subsp. paratuberculosis unique diagnostic targets and protocols that may provide a facile approach to unambiguously detect M. avium subsp. paratuberculosis. Future comparative genomics and pathogenesis studies with the bovine and ovine M. avium subsp. paratuberculosis strains and with human M. avium subsp. paratuberculosis and MAC isolates will be required to further elucidate the natural history of this organism.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles: