|
|
||||||||
1 Department of Veterinary Pathobiology, The Royal Veterinary and Agricultural University, Stigbøjlen 4, DK-1870 Frederiksberg C, Denmark
2 Center for Biological Sequence Analysis, BioCentrum-DTU, Technical University of Denmark, Building 208, DK-2800 Lyngby, Denmark
3 Department of Wildlife, The National Veterinary Institute, SE-751 89 Uppsala, Sweden
6 Corresponding author (email: miki{at}life.ku.dk)
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Based on phenotypic traits, ruminant strains classified as trehalose-negative M. [P.] haemolytica could be divided into 12 biogroups or taxa (Bisgaard and Mutters, 1986; Bisgaard et al., 1986). Interestingly, the 16S rRNA phylogeny indicates that different biogroups are present among two or more species, suggesting that members of the same biogroup are not necessarily closely related (Angen et al., 1999; Larsen et al., 2007). Previous works have shown that the ß-hemolytic phenotype correlates with the pathogenic potential of strains belonging to genus Mannheimia. The results of Angen et al. (1999) indicated that most strains isolated from lesions were ß-hemolytic on bovine blood agar plates (except strains belonging to M. granulomatis), whereas most strains belonging to the nonpathogenic species M. ruminalis were nonhemolytic. Larsen et al. (2007) used Southern blot analysis, sheep blood agar plates, and a representative strain collection to show that the lktA genotype and the corresponding ß-hemolytic phenotype were present in all Mannheimia species, including all strains belonging to M. haemolytica and M. glucosida, M. granulomatis, and M. varigena, but only in a fraction of strains belonging to M. ruminalis. Murphy et al. (1995) showed that the ß-hemolytic phenotype of M. haemolytica correlates with the presence of an active leukotoxin (lkt) operon on the chromosome. The structural leukotoxin (LktA) protein of M. haemolytica, encoded by lktA (Lo et al., 1987), belongs to the E. coli HlyA-like subfamily of cytotoxic RTX (repeats in toxin) proteins and is essential in both evasion and exploitation of the adaptive immune system during pulmonary infection (Petras et al., 1995; Tatum et al., 1998; Highlander et al., 2000). These results support the predicted correlation between pathogenic potential and the ß-hemolytic phenotype.
In this study, we characterize eight atypical Mannheimia isolates obtained from lesions in roe deer (Capreolus capreolus) by using a polyphasic strategy that combines phenotypic and genotypic methods. The implications of our results for proper identification and classification of Mannheimia strains are discussed.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The strains used in this study and their geographic origin, time of isolation, and corresponding pathological lesions are given in Table 1
. A total of eight atypical Mannheimia strains were isolated from lesions in roe deer submitted to the Department of Wildlife, The National Veterinary Institute, Uppsala, Sweden for post-mortem examination.
|
All isolates were phenotypically characterized as previously described by Bisgaard et al. (1991). A total of 82 different phenotypic tests have been used. To determine presence/inactivation of the lkt operons, we screened the strains for the ß-hemolytic phenotype on sheep blood agar plates as previously described (Murphy et al., 1995).
Reconstruction of phylogenetic trees
The relationships of the strains in this study were inferred with 16S rRNA; these sequences have been used successfully for systematic studies in this group (Angen et al., 1999; Larsen et al., 2007). Sequences for 41 strains were obtained from GenBank (Table 2
). The primers previously described by Angen et al. (1999) were used to amplify and directly sequence 16S rRNA from two roe deer isolates (B 234/94 and B 1829/84). The 16S rRNA sequences were aligned using both Dialign 2 (Morgenstern et al., 1999) and the NAST server (DeSantis et al., 2006), which uses a prealigned set of 16S rRNA sequences as the basis for aligning the users sequences. The two methods yielded identical alignments. This was expected, because the average pairwise sequence similarity in the full set of 43 sequences was 97.5% and even the least similar pair of sequences was 92.5% identical, supporting that alignment errors were unlikely. Phylogenetic trees were reconstructed with the use of Bayesian techniques as implemented in the program MrBayes version 3.1.1 (Huelsenbeck and Ronquist, 2001; Ronquist and Huelsenbeck, 2003). The best-fitting model was GTR+ invgamma based on the Akaike Information Criterion (AIC) computed by the program MrModeltest (version 2.2, Department of Systematic Zoology, Uppsala University, Uppsala, Sweden; http://www.ebc.uu.se/systzoo/staff/nylander.html). Markov chain Monte Carlo (MCMC) was run for 10,000,000 generations with four chains while sampling once every 100 generations. Convergence was confirmed by comparing the results of two independent runs. The program Tracer (version 1.3; Department of Zoology, Oxford University, Oxford, UK; http://evolve.zoo.ox.ac.uk/software.html?id=tracer) was used to determine burn-in and also for further confirmation of proper mixing and adequate run length. A burn-in of 1,000,000 generations (10,000 samples) was used in both cases. The distribution of 16S rRNA trees resulting from MCMC was summarized in the form of a consensus tree with all compatible bipartitions included by using the authors own software in a manner that is essentially identical to what is obtained when using MrBayes sumt command with the setting contype = allcompat, except that branch lengths were averaged over all trees; setting the branch length to zero for those trees that did not contain the corresponding bipartition, as suggested by Felsenstein (2003). The tree was rooted based on a maximum-likelihood analysis reported by Angen et al. (1999).
|
We analyzed the distribution of the lktA gene among roe deer isolates by using Southern blot as previously described by Larsen et al. (2007). The primers manpop_UP and manpop_ DOWN previously described by Larsen et al. (2007) were used to amplify the +845/+1,302 region of lktA. The reaction conditions were 2.5 U Taq polymerase, 16 mM (NH4)2SO4, 67 mM Tris-HCl, 0.01% Tween-20, 2.5 mM Mg2SO4, each primer at 0.5 mM, and each nucleotide at 0.1 mM. The cycling conditions were initial denaturation at 94 C followed by 25 cycles of 94 C for 30 s, 52 C for 30 sec, and 72 C for 30 sec, finishing with extension at 72 C for 10 min.
Sequence data
Sequences have been deposited in the GenBank database (accession numbers (AY425296 and AY425297).
| RESULTS |
|---|
|
|
|---|
The strains used in this study and their geographic origin, time of isolation, and corresponding lesions are given in Table 1
. It should be noted that the strains were isolated from six different locations within a 300-km radius from Norrkö ping (58°36'N, 16°11'E) in southeastern Sweden. The isolates were obtained during a period from 1980 to 1995. Six of the isolates were obtained from animals <1 yr old, whereas a single isolate originated from a 7-yr-old female. Lesions varied from acute purulent bronchopneumonia (n = 2), purulent pneumonia and pleuritis (n = 2), glossitis (n = 1), purulent stomatitis and glossitis (n = 1), and phlegmon and abscess formation (n = 1).
Phenotypic characterization
All isolates were gram-negative, non-motile, catalase and oxidase positive, and decomposed glucose fermentatively in Hugh and Leifson medium. All isolates were porphyrin, alanine aminopeptidase and phosphatase positive, and reduced nitrate. In addition, acid was produced from glycerol, D-ribose, D-xylose, meso-inositol, D-mannitol, D-fructose, D-galactose, D-glucose without gas formation, lactose, sucrose, and raffinose. Positive reactions were also obtained for all isolates in ONPG (ß-galactosidase) and PGUA (ß-glucuronidase). The ß-hemolytic phenotype on bovine and ovine blood agar plates was present in six isolates, whereas isolate B 1828/84 was nonhemolytic. Late acid production from L-fucose was observed for two isolates, whereas a single isolate showed a positive reaction in the ONPX test (ß-xylosidase).
A negative reaction was observed for all isolates in the following tests: symbiotic growth, Simmons citrate, malonate, H2S/TSI, KCN growth, MR 37 C, VP 37 C, formation of gas from nitrate, urease, arginine dihydrolase, lysine and ornithine decarboxylase, phenylalanine deaminase, indole, geletinase, hydrolysis of Tween 20 and Tween 80, McConkey agar growth, pigment formation, and production of acid from meso-erythritiol, adonitol, D-arabitol, xylitol, L-arabinose, D-arabinose, L-xylose, dulcitol, D-sorbitol, D-fucose, D-mannose, L-rhamnose, L-sorbose, cellobiose, maltose, D-melibiose, trehalose, D-melezitose, dextrin, D-glycogen, inulin, aesculin, amygdalin, arbutin, gentiobiose, salicin, D-turanose, and ß-N-CH3-glucosamid. Negative reactions were also observed in NPG (ß-glucosidase), PNPG (
-glucosidase), ONPF (
-fucosidase),
-galactosidase, and
-mannosidase. Phenotypic characters separating the isolates from biogroups belonging to M. granulomatis and the type strain of M. ruminalis (HPA92T) are given in Table 3
.
|
BLAST searches of the NCBI non-redundant (nr) database with the 16S rRNA sequences of isolates B 234/94 and B 1829/84 as queries suggested that they are most closely related to M. granulomatis. Because the closest BLAST hits are not necessarily each others closest relative (Koski and Golding, 2001), we reconstructed the phylogeny with the use of these sequences, along with sequences of 41 Mannheimia strains. The 16S rRNA data set included 1,257 bp after removal of ambiguous bases. Both runs converged and the consensus tree with all compatible bipartitions included and with branch lengths averaged over all trees, setting the branch length to zero for those trees that did not contain the corresponding bipartition, is shown in Figure 1
, along with posterior probabilities (PP). This phylogram is largely topologically concordant with the Bayesian analysis reported by Larsen et al. (2007) and suggest that the two roe deer isolates form a cluster within M. granulomatis (PP = 0.90%). There was strong PP support for monophyly of the roe deer isolates (PP = 100%), supporting that they belong to a distinct lineage of M. granulomatis (PP = 100%).
|
The lktA genotype and the corresponding ß-hemolytic phenotype are given in Table 1
. It should be noted that isolate P 4737/80 died during storage and could not be included in these analyses. Southern blot analysis indicated that the lktA gene was present in all seven isolates. The +845/+1,302 region of the lktA gene was amplified by PCR in the six ß-hemolytic isolates. We could not amplify lktA from the nonhemolytic isolate B 1829/84, presumably because substitutions have resulted in imperfect matches of the primers.
| DISCUSSION |
|---|
|
|
|---|
The bovine M. granulomatis strains W4672/1 (biogroup 3J) and P1135/26T (M. [P.] granulomatis) were their closest neighbors on the 16S rRNA tree, whereas the leporine strains Ph13 (Bisgaard taxon 20 biovar 1) and BJ1680.3 (Bisgaard taxon 20 biovar 2) constituted a more distantly related group (Fig. 1
). The results presented here suggest the divergence of at least three distinct M. granulomatis lineages that may have adapted to cervine, bovine, and leporine hosts, respectively.
Prior to this work, nonhemolytic Mannheimia strains had never been isolated outside the bovine and ovine rumen and they had never been associated with disease conditions, supporting the predicted correlation between pathogenic potential and the ß-hemolytic phenotype. However, our analyses showed that isolate B 1829/84 was nonhemolytic, although the lktA gene was present on the chromosome. The presence of an inactivated lkt operon has been observed in a single M. ruminalis strain (HPA88), whereas chromosomal deletions of ~7,400 bp, corresponding to the entire lkt operon, have occurred in the remaining nonhemolytic M. ruminalis strains (Larsen et al., in press). If one assumes that the evolving M. granulomatis populations adapted to different ecologic niches (host species and anatomic niches) offering different conditions for growth, then this diversity could imply that different pools of unused genes have been lost in those populations. This hypothesis predicts the existence of variable morphologic and physiologic traits between strains belonging to the same species and constitutes the most parsimonious explanation for the observed phenotypic variation within M. granulomatis and other Mannheimia species (Angen et al., 1997a, b). These results support the use of a polyphasic strategy, combining phenotypic and genotypic methods, for proper identification and classification of Mannheimia strains.
| ACKNOWLEDGMENTS |
|---|
| FOOTNOTES |
|---|
| LITERATURE CITED |
|---|
|
|
|---|
———, B. AALBAEK, E. FALSEN, J. E. OLSEN, AND M. BISGAARD. 1997b. Relationships among strains classified with the ruminant Pasteurella haemolytica-complex using quantitative evaluation of phenotypic data. Zentralblatt für Bakteriologie 285: 459–479.[Medline]
———, D. A. CAUGANT, J. E. OLSEN, AND M. BISGAARD. 1997c. Genotypic relationships among strains classified under the (Pasteurella) haemolytica-complex as indicated by ribotyping and multilocus enzyme electrophoresis. Zentralblatt für Bakteriologie 286: 333–354.[Medline]
———, R. MUTTERS, D. A. CAUGANT, J. E. OLSEN, AND M. BISGAARD. 1999. Taxonomic relationships of the [Pasteurella] haemolytica complex as evaluated by DNA-DNA hybridizations and 16S rRNA sequencing with proposal of Mannheimia haemolytica gen. nov., comb. nov., Mannheimia granulomatis comb. nov., Mannheimia glucosida sp. nov., Mannheimia ruminalis sp. nov. and Mannheimia varigena sp. nov. International Journal of Systematic Bacteriology 49 Pt 1: 67–86.
BISGAARD, M., AND R. MUTTERS. 1986. Re-investigations of selected bovine and ovine strains previously classified as Pasteurella haemolytica and description of some new taxa within the Pasteurella haemolytica-complex. Acta Pathologica, Microbiologica et Immunologica Scandinavica B 94: 185–193.
———, J. E. PHILLIPS, AND W. MANNHEIM. 1986. Characterization and identification of bovine and ovine Pasteurellaceae isolated from the oral cavity and rumen of apparently normal cattle and sheep. Acta Pathologica, Microbiologica et Immunologica Scandinavica B 94: 9–17.
———, S. B. HOUGHTON, R. MUTTERS, AND A. STENZEL. 1991. Reclassification of German, British and Dutch isolates of so-called Pasteurella multocida obtained from pneumonic calf lungs. Veterinary Microbiology 26: 115–124.[Medline]
DAVIES, R. L., S. ARKINSAW, AND R. K. SELANDER. 1997. Evolutionary genetics of Pasteurella haemolytica isolates recovered from cattle and sheep. Infection and Immunity 65: 3585–3593.
DESANTIS, T. Z. JR., P. HUGENHOLTZ, K. KELLER, E. L. BRODIE, N. LARSEN, Y. M. PICENO, R. PHAN, AND G. L. ANDERSEN. 2006. NAST: A multiple sequence alignment server for comparative analysis of 16S rRNA genes. Nucleic Acids Research 34: W394–399.
DEVRIESE, L. A., M. BISGAARD, J. HOMMEZ, E. UYTTEBROEK, R. DUCATELLE, AND F. HAESEBROUCK. 1991. Taxon 20 (Fam Pasteurellaceae) infections in European brown hares (Lepus europaeus). Journal of Wildlife Disease 27: 685–687.[Abstract]
FELSENSTEIN, J. 2003. Inferring phylogenies. Sinauer Associates, Sunderland, Massachusetts, 664 pp.
HIGHLANDER, S. K., N. D. FEDOROVA, D. M. DUSEK, R. PANCIERA, L. E. ALVAREZ, AND C. RINEHART. 2000. Inactivation of Pasteurella (Mannheimia) haemolytica leukotoxin causes partial attenuation of virulence in a calf challenge model. Infection and Immunity 68: 3916–3922.
HUELSENBECK, J. P., AND F. RONQUIST. 2001. MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics 17: 754–755.
KOSKI, L. B., AND G. B. GOLDING. 2001. The closest BLAST hit is often not the nearest neighbor. Journal of Molecular Evolution 54: 540–542.
LARSEN, J., A. G. PEDERSEN, H. CHRISTENSEN, M. BISGAARD, Ø. ANGEN, P. AHRENS, AND J. E. OLSEN. 2007. Evidence for vertical inheritance and loss of the leukotoxin operon in genus Mannheimia. Journal of Molecular Evolution. http://dx.doi.org/10.1007/s00239-006-0065.3.
LO, R. Y., C. A. STRATHDEE, AND P. E. SHEWEN. 1987. Nucleotide sequence of the leukotoxin genes of Pasteurella haemolytica A1. Infection and Immunity 55: 1987–1996.
MORGENSTERN, B. 1999. DIALIGN 2: Improvement of the segment-to-segment approach to multiple sequence alignment. Bioinformatics 15: 211–218.
MURPHY, G. L., L. C. WHITWORTH, K. D. CLINKENBEARD, AND P. A. CLINKENBEARD. 1995. Haemolytic activity of the Pasteurella haemolytica leukotoxin. Infection and Immunity 63: 3209–3212.
PETRAS, S. F., M. CHIDAMBARAM, E. F. ILLYES, S. FROSHAUER, G. M. WEINSTOCK, AND C. P. REESE. 1995. Antigenic and virulence properties of Pasteurella haemolytica leukotoxin mutants. Infection and Immunity 63: 1033–1039.
RIET-CORREA, F., S. L. LADEIRA, G. B. ANDRADE, AND G. R. CARTER. 2000. Lechiguana (focal proliferative fibrogranulomatous panniculitis) in cattle. Veterinary Research Communications 24: 557–572.[Medline]
RONQUIST, F., AND J. P. HUELSENBECK. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics 19: 1572–1574.
TATUM, F. M., R. E. BRIGGS, S. S. SREEVATSAN, E. S. ZEHR, H. S. LING, L. O. WHITELEY, T. R. AMES, AND S. K. MAHESWARAN. 1998. Construction of an isogenic leukotoxin deletion mutant of Pasteurella haemolytica serotype 1: Characterization and virulence. Microbial Pathogenesis 24: 37–46.[Medline]
Received for publication 6 July 2006.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |