|
|
||||||||



* Department of Veterinary and Biomedical Sciences,
Department of Veterinary Population Medicine, and
Department of Animal Science, University of Minnesota, St. Paul 55108
1 Corresponding author: beyxx001{at}umn.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: mastitis Klebsiella diversity
| INTRODUCTION |
|---|
|
|
|---|
Coliform bacteria can cause clinical or subclinical mastitis. Cows with clinical mastitis can have abnormal milk, swollen quarters, anorexia, and fever. In subclinical mastitis there are no visible signs of disease. Inflammation can only be monitored by in vitro diagnostic tests. Previously, coliform bacteria have been shown to be responsible for approximately half of the acute mastitis episodes (Erskine et al., 1991). The control and prevention of coliform mastitis remains a challenge to dairy producers. Control measures that are applied for contagious mastitis pathogens are often ineffective in controlling mastitis caused by coliform bacteria. Presently, the most efficient strategy for controlling coliform mastitis is prevention through appropriate management and prophylactic immunization (Bradley and Green, 2000).
The number of mastitis cases caused by coliform bacteria has increased worldwide (Schukken et al., 1989). Specifically, K. pneumoniae has been found to affect many dairy herds (Silva and Costa, 2001; Sampimon et al., 2006). Braman et al. (1973) found 33 capsular types of K. pneumoniae in 12 herds. Within each herd as many as 13 capsular types were found, indicating there were different variants of the same organism present. In another study, Kikuchi et al. (1995) reported diversity in plasmid profiles among 47 K. pneumoniae isolates from mastitic milk. However, plasmid profile typing is not as discriminatory as chromosomal DNA genotyping because plasmids may be lost, gained spontaneously, or exchanged between strains, making it difficult to determine the potential number of variants present.
Mastitis caused by K. pneumoniae can be particularly severe compared with mastitis caused by Escherichia coli, due to its poor response to antibiotic therapy, rapid evolution to toxic shock, and death (Silva and Costa, 2001). Infection of the gland with Klebsiella species has been reported to result in higher losses than infections due E. coli (Grohn et al., 2004). Thus, infection with Klebsiella spp. profoundly affects the profitability of dairy farmers due to discarded milk, costs of antibiotic treatment, and extra labor, death, or culling of infected animals and decreased fertility (Hansen et al., 2004). Genotyping of K. pneumoniae isolated from clinical mastitis is expected to increase our understanding of the epidemiology of this organism and perhaps assist in the development of appropriate strategies for the prevention and control of udder infections in dairy cows. The goal of this study was to measure the genetic diversity of K. pneumoniae from cases of clinical mastitis within a dairy herd, using 3 genotyping methods [repetitive DNA sequence genotyping-PCR (rep-PCR), pulsed-field gel electrophoresis (PFGE), and multilocus sequence typing (MLST)] to fingerprint and genotype these isolates.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Repetitive DNA Sequence Genotyping
Bacterial genomic DNA was extracted using Prep-Man Ultra Reagent (Applied Biosystems, Foster City, CA). Briefly, frozen cultures were streaked onto MacConkey agar plates and incubated. Two individual colonies were suspended in 300 µL of PBS and pelleted by centrifugation for 3 min at 7,500 x g. The supernatant was discarded, 200 µL of PrepMan was added, and the tubes placed in a heat block for 10 min at 100°C. After incubation, the solution was fractionated by centrifugation for 3 min at 7,500 x g, and the supernatant containing DNA and an equal amount of sterile DNase-free water was added. The rep-PCR fingerprints were obtained using a boxA1R primer (5'CTACGGCAAGGCGACGCTGACG 3'). The 25-µL PCR mixture contained 25 mM MgCl2, 10x Buffer II, 25 pmol of boxA1R primer, 100 mM dNTP mix, AmpliTaq Gold DNA polymerase, and 2 µL of DNA template. Thermocycler PCR conditions were 95°C for 7 min, followed by 30 cycles consisting of 94°C for 1 min, 66°C for 8 min, 71°C for 1 min, and an extension of 71°C for 15 min followed by storage at 4°C. These conditions were previously described by Goldberg et al. (2006). Samples were separated by electrophoresis using a 2% agarose gel at 4°C for 5 to 6 h at 62 V (6 V/cm) and stained in ethidium bromide/1x Tris-acetate-EDTA buffer. Gel images were electronically captured using Labworks 4.0 Image Acquisition and Analysis Software (UVP Inc., Upland, CA).
Pulsed-Field Gel Electrophoresis
The PFGE protocol was based on procedures previously described by Cho et al. (2006) in addition to those described in the CHEF Genomic DNA Plug Kit (BioRad, Hercules, CA). Klebsiella pneumoniae isolates were grown overnight in 10 mL of brain heart infusion broth at 37°C. Cells were pelleted by centrifugation and resuspended in 10 mL of cell suspension buffer (100 mM Tris:100 mM EDTA, pH 8.0) followed by centrifugation. The cells were washed a second time and resuspended in 5 mL of cell suspension buffer. Agarose plugs were prepared by adding 250 µL of 2% CleanCut agarose (BioRad) to 250 µL of the cell suspension with gentle mixing at 56°C. One hundred microliters of this mixture was transferred to plug molds and allowed to solidify. The solidified agarose plugs were incubated in 300 µL of freshly prepared lysozyme solution (0.01 mg/µL final concentration) for 2 h at 37°C. After this, plugs were incubated overnight at 56°C in proteinase K solution (0.24 U/µL final concentration). Following incubation, the plugs were washed with sterile water with Tris-EDTA buffer pH 8.0 at 50°C. After washing, half of the plugs were digested with XbaI (New England BioLabs Inc., Beverly, MA) for 1.5 to 2 h at 37°C. After digestion, the Xba buffer was removed, and about one-quarter of each plug was loaded into each well. A lambda DNA ladder (New England BioLabs) was used as a molecular size standard. Electrophoresis was carried out in 0.5x Tris-borate-EDTA buffer at 14°C with an initial switch time of 1 s and a final switch time of 30 s for 23 h at 6V/cm. The gel was stained with ethidium bromide and destained with distilled water. Gel images were viewed and captured with Labworks 4.0 Image Acquisition and Analysis Software (UVP Inc.).
Rep-PCR and PFGE Fingerprint Analysis
Gel images were analyzed using BioNumerics software version 4.0 (Applied Maths BVBA, Sint-Martens-Latem, Belgium). For rep-PCR, the band positions in each gel were normalized using the 1 kb DNA ladder (Hi-Lo DNA Marker, Minnesota Molecular Inc., Minneapolis, MN), and for PFGE the DNA band normalization was accomplished by using the 1,000 kb lambda DNA ladder (New England BioLabs Inc., Beverly, MA). Band matching and isolate similarity using rep-PCR and PFGE was accomplished using Dice band-based coefficient of similarity, which provides the most accurate similarity results when compared with visual inspection of the fingerprint patterns (Carriço et al., 2005).
The Dice similarity coefficient formula used in the analysis was
![]() |
where Sij = similarity; ni = number of unique bands occurring in pattern i; nj = number of unique bands occurring in pattern j; and nij = number of bands shared between the 2 patterns.
Tolerance and optimization settings were 1.5 and 1.0% for rep-PCR and 3.0 and 1.5% for PFGE. These settings were determined by evaluating different ranges (0 to 6%) and selecting the ones with greatest accuracy when compared with visual inspection of fingerprints in the gel picture (Carriço et al., 2005). A dendrogram was constructed using the unweighted pair group method with arithmetic means, which employs a sequential clustering algorithm in which the relationships are identified in order of similarity and the dendrogram is built in a stepwise manner (Durbin et al., 1998). The PCR reproducibility was verified from 16 distinct isolates by freshly extracting DNA 4 separate times. A dendrogram containing 4 fingerprints of the same isolate was constructed for K. pneumoniae using the BioNumerics software and the settings mentioned above. The percent similarity among fingerprints for each dendrogram was recorded and the average calculated across dendrograms. Reproducibility showed an average of 93.6% fingerprint similarity for the same isolate and was used, along with cluster pattern, to conservatively establish the 90% similarity cut-off for genotype definition used in this study. No other mathematical analysis was used in establishing the 90% cut-off.
Multilocus Sequence Typing
Amplification of 7 housekeeping genes, rpoB (beta-subunit of RNA polymerase), gapA (glyceraldehyde 3-phosphate dehydrogenase), mdh (malate dehydrogenase), pgi (phosphoglucose isomerase), phoE (phosphorine E), infB (translation initiation factor 2), and tonB (periplasmic energy transducer), was accomplished using PCR. Two K. pneumoniae reference strains, MGH 78578 (ATCC 700721 and in MLST database) and strain IA565 (which is avirulent by intratracheal inoculation in mice; kindly provided by Steve Clegg, Iowa City, IA), were included as positive controls. Primer sets for gene amplification and sequencing were described previously (Diancourt et al., 2005). The PCR amplification was performed using AmpliTaq Gold DNA polymerase (Applied Biosystems) and standard thermocycler parameters. The annealing temperatures of each primer were described at the K. pneumoniae MLST database (http://pubmlst.org/kpneumoniae). Amplified products were purified by ultrafiltration with Montage PCR centrifugal filter devices (Millipore, Bedford, MA). Nucleotide sequencing was performed at the Advanced Genetic Analysis Center (AGAC) University of Minnesota, St. Paul. A minimum of 2 chromatograms was performed for each amplified gene. The software MEGA 3.1 (Kumar et al., 2004) was used for editing the nucleotide sequences and to evaluate the significance of branching within trees by bootstrap analysis (500 replicates) using the maximum parsimony method. Each distinct sequence within a locus was assigned a unique allele number. Each nucleotide sequence was analyzed by BLASTn algorithm against the K. pneumoniae MLST database (http://pubmlst.org/kpneumoniae) in search of a matching allele number. Nucleotide sequences having no matching allele in the MLST database were considered new and had their respective allele numbers and sequence type (ST) assigned by the curator of the K. pneumoniae MLST database Web site. Once all of the allele numbers and sequence types were assigned, cluster analysis was performed by BioNumerics software version 4.0 (Applied Maths BVBA, Martens-Latem, Belgium) using categorical coefficient. In this method, alleles containing identical nucleotide sequence were considered a match and alleles containing any nucleotide difference (irrespective of the number of repeats) were considered a no match. Additionally, the numbers assigned to each allele type does not imply any relationship between ST.
Statistical Analysis
Spearman rank coefficient was used to calculate the correlation between a matrix of fingerprints similarity and a matrix of temporal distances (number of days between dates of sample collection) using the Mac-Anova statistical analysis program (School of Statistics, University of Minnesota, Minneapolis, MN). A P-value of
0.05 was considered significant for the temporal correlation. Simpsons diversity index was used for measurement of genetic diversity and discriminatory power of the typing techniques. Calculation of Simpsons diversity index was performed using the BioNumerics software and it was based on the following formula:
![]() |
where N is the total number of isolates in the sample population, S the total number of genotypes described, and nj is the number of isolates belonging to the jth genotype. The genetic diversity is based on the probability that 2 K. pneumoniae isolates chosen at random will be of a different genotype. Higher Simpsons diversity index indicates a high discriminatory power of the technique and genetic diversity within the population (Hunter and Gaston, 1988).
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
Using a 90% cut-off for similarity in the rep-PCR banding patterns, 21 distinct fingerprint patterns were identified among the 84 isolates and were assigned a genotype. The Simpsons diversity index was 93.5%, indicating a wide genetic diversity and a high degree of discrimination capability of rep-PCR. Each genotype was composed of between 1 and 12 isolates. Analysis of 3 distinct isolates originating from the same cow showed identical rep-PCR patterns 46.4% of the time. Two of the colonies had the same genotype pattern 35.7% of the time. The remaining 17.9% showed different genotype patterns for all isolates indicating that more than one K. pneumoniae variant can be isolated from the mastitic bovine mammary gland. The possibility of sample contamination exists, but to reduce this possibility only those samples from which K. pneumoniae was the only organism isolated were included in this study and the individual responsible for sample collection had been extensively trained in how to aseptically collect milk samples. It is possible that errors inherent in the PCR technique may have partially accounted for the different genotypes encountered in K. pneumoniae isolated from the same mastitis case sample.
Two dairy cows participating in this study died with acute mastitis. Isolates originating from cows that died (716 and 1481) possessed rep-PCR genotypes distinct from each other. Other cows infected with the same genotypes did not die and remained in the herd. It is possible that in addition to Klebsiellas virulence traits, factors such as stage of lactation, immune and nutritional status, age, and genetics may have contributed to the severity of these mastitis cases.
Spearman rank coefficient was used to calculate the correlation between a matrix of genetic similarity and a matrix of temporal distances. The association of the genetic similarity matrix and the temporal distance matrix could provide more information on the dynamics of K. pneumoniae mastitis. The results showed a significant negative association between the 2 matrices (R = –0.1525; P = 0.0246), indicating that isolates that were collected closer in time had higher similarity than those collected farther apart in time. The correlation was not very strong but was significantly different from zero. The great genetic diversity may explain the low correlation observed. The results suggest that temporal distance seems to play a role in the K. pneumoniae genotypes isolated from clinical cases of mastitis.
Pulsed-Field Gel Electrophoresis Genotyping
To verify results obtained by rep-PCR, a subset of 43 isolates was selected to be typed by PFGE. The enzyme XbaI was used for PFGE typing because it is recommended by the Centers for Disease Control and Prevention. The selection of isolates for PFGE typing was based on similarities in rep-PCR pattern. Out of 21 rep-PCR genotypes, 15 were represented in PFGE. We also included isolates that originated from the same clinical case: 3 per case (n = 12), 2 per case (n = 16). The remaining isolates were singly represented (n = 15), meaning that only 1 isolate per clinical case was analyzed. Using a 90% cut-off, 23 distinct PFGE genotypes were detected in the 43 isolates and each genotype contained 1 to 6 isolates. Klebsiella pneumoniae isolates typed by PFGE showed greater genetic diversity than observed with rep-PCR. It was interesting to note that, with both techniques, some of the isolates from the same clinical mastitis cases exhibited more than 1 genotype pattern (29.0 and 17.9% by PFGE and rep-PCR, respectively), supporting the result that more than one K. pneumoniae variant can be isolated from the infected mammary gland and confirming the high discriminatory power of PFGE. Isolates that originated from the same mastitis cases (45.2%) clustered similarly by rep-PCR and by PFGE. The Simpsons diversity index for PFGE was 96.1%, indicating the wide diversity of K. pneumoniae within the dairy herd. Once again, isolates from the 2 cows that died of acute mastitis showed distinct genotype patterns by PFGE. The isolates clustering with the acute mastitis cases by PFGE were different from those identified by rep-PCR. Pulsed-field gel electrophoresis showed higher discriminatory power (Simpsons diversity index = 96.1%) when compared with rep-PCR (Simpsons diversity index = 93.5%).
Multilocus Sequence Typing Genotyping
Due to the great genetic diversity observed and the possible subjective interpretation of the results by DNA band-based techniques, genetic diversity among the K. pneumoniae isolates was analyzed using MLST, which is characterized as having high reproducibility and discriminatory power. A subset of 29 isolates was analyzed by MLST. Thirteen PFGE genotypes were typed by MLST, which included isolates from the same mastitis cases, 3 per case (n = 9) and 2 per case (n = 8), and single isolates from different mastitis cases (n = 12). Twenty-nine K. pneumoniae mastitis isolates and 2 type strains used as positive controls (MGH78578 and IA565) were analyzed by MLST. This analysis resulted in a total of 23 ST measured as a composite of 7 genetic loci studied. Excluding ST formed by the 2 reference strains, MLST analysis had a Simpsons diversity index of 97.0%. This high level of diversity is supported by bootstrap analysis, which also had high values. Clinical mastitis isolates composed 21 ST, whereas 2 distinct ST composed the 2 reference strains. As in rep-PCR and PFGE, isolates originating from the same clinical mastitis case had distinct ST. Five cases of mastitis with 2 isolates per case exhibited distinct ST for each isolate. Three isolates originating from the same clinical mastitis sample presented the same ST. Comparable variability was shown by rep-PCR for the same set of samples in which 3 of them had the same rep-PCR genotype and 2 had distinct rep-PCR genotype. Two samples having distinct ST by MLST showed the same genotype by rep-PCR. One of the isolates (716A), which caused acute mastitis (evolution to toxemic shock and death), shared the same ST with isolate 2261A. Both isolates clustered together by MLST and rep-PCR.
Nucleotide variation was observed in all 7 genes with each gene having 4 to 13 distinct allelic types. Five ST (98, 104, 110, 112, and 116) and 5 of the allelic types found were new, novel, and did not present patterns found in the MLST database for K. pneumoniae (http://pubmlst.org/kpneumoniae). The genes rpoB, gapA, and infB contained the least amount of allelic variation with each having 4 to 5 different allelic types. The mdh and pgi genes contained 6 allelic types. Genes phoE and tonB contained the largest amount of allele variation with 9 and 13 allelic types, respectively. The novel ST and allelic types are shown in Table 1
. Novel allelic types were found for all genes except rpoB, pgi, and phoE. Two new allelic types were found for tonB and gapA; mdh and infB had one new allelic type each. Considering that all isolates originated from a single herd, a significant amount of nucleotide diversity was observed in the housekeeping genes. Our results support those previously reported by Diancourt et al. (2005) that tonB exhibited the greatest amount of nucleotide variation. Three of the clinical mastitis isolates from cows 716, 2,261, and 3,825 shared ST (14 and 45) with human clinical isolates and one fecal isolate listed in the MLST database for K. pneumoniae. It is possible that isolates from clinical cases of mastitis have virulence attributes capable of causing human disease.
|
|
We were unable to identify 1 or 2 unique clones within a K. pneumoniae population associated with clinical mastitis. We have shown, with the use of 3 distinct methods (rep-PCR, PFGE, and MLST) that clinical mastitis can be caused by a genetically diverse group of K. pneumoniae and that more than 1 genotype can be isolated from cow(s) with clinical mastitis.
Struve and Krogfelt (2004), using a mouse infection model, found that K. pneumoniae isolates of environmental origin were as virulent as those isolated from a human infection. Klebsiella pneumoniae had the highest growth in bedding materials when compared with other environmental bacterial species (Zehner et al., 1986). Studies have shown that Klebsiella species had extended survival rates in the environment and were more resistant to solar radiation than Escherichia coli (McCambridge and McMeekin, 1981; Lopez-Torres et al., 1987). It is possible that the hearty and ubiquitous nature of Klebsiella associated with the particular conditions within a dairy environment contributed to the large genetic diversity observed, which is consistent with previous reports (Braman et al., 1973; Kikuchi et al., 1995). Other environmental mastitis causing bacteria, such as Streptococcus uberis, have also shown a wide genetic variation (Wieliczko et al., 2002). Similar genetic diversity has been found in isolates of E. coli (Yang et al., 2004). The results of the present investigation help explain, in part, the difficulty in controlling and preventing environmental mastitis. The presence of such a genetically complex and diverse microbial population makes it difficult to target specific subpopulations of bacteria for control within the dairy environment. Klebsiella pneumoniae isolated from the dairy herd in this study is genetically diverse and may represent a snapshot of the diversity in the population. Still, it is possible that there are many potential sources of K. pneumoniae within a dairy herd. A larger study including isolates from several herds is necessary to confirm such findings.
Cleanliness in a dairy facility can be difficult to maintain due to the constant production of organic matter such as feces, milk, and urine and the influence of external factors such as the weather. The presence of such a genetically complex microbial population associated with the challenges of maintaining good sanitation makes it difficult to control environmental mastitis pathogens. Klebsiella pneumoniae in the environment is subject to diverse conditions (temperature, pH, nutrients, moisture, and exposure to other microbes), which may lead to the horizontal transfer of DNA, mutations, or both, which could account for some of the genetic variation observed (Rocha, 2004). Additionally, K. pneumoniae could have variable genome sizes that account for some of the differences observed in rep-PCR and PFGE. This hypothesis is based on E. coli genome size variations (Bergthorsson and Ochman, 1998). The combination of all factors mentioned above may be responsible for the wide genetic variation encountered in this K. pneumoniae population.
The majority of molecular epidemiologic studies of K. pneumoniae have been performed with isolates from humans, especially outbreaks in neonatal and adult intensive care units. Contrary to our findings, K. pneumoniae isolated from hospitals follow a clonal pattern, and a single clone was usually associated with an outbreak. Some epidemiologic studies have reported up to 4 distinct K. pneumoniae clones circulating in hospital outbreaks (Dashti et al., 2006). Unlike a hospital, a dairy facility is subjected to factors that are difficult to control. Some of these factors include the weather, animal crowding, bedding management, moisture, cleanliness of installations, lack of consistency in pre-milking cow preparation, introduction of new animals into the herd, birds, and rodents. It is possible that a very diverse population of K. pneumoniae would have access to dairy facilities resulting in increased genetic diversity.
| CONCLUSION |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication November 20, 2006. Accepted for publication April 18, 2007.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. A. Munoz, G. J. Bennett, C. Ahlstrom, H. M. Griffiths, Y. H. Schukken, and R. N. Zadoks Cleanliness Scores as Indicator of Klebsiella Exposure in Dairy Cows J Dairy Sci, October 1, 2008; 91(10): 3908 - 3916. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. G. Paulin-Curlee, S. Sreevatsan, R. S. Singer, R. Isaacson, J. Reneau, R. Bey, and D. Foster Molecular Subtyping of Mastitis-Associated Klebsiella pneumoniae Isolates Shows High Levels of Diversity Within and Between Dairy Herds J Dairy Sci, February 1, 2008; 91(2): 554 - 563. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. A. Munoz, F. L. Welcome, Y. H. Schukken, and R. N. Zadoks Molecular Epidemiology of Two Klebsiella pneumoniae Mastitis Outbreaks on a Dairy Farm in New York State J. Clin. Microbiol., December 1, 2007; 45(12): 3964 - 3971. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |