Previous Article | Next Article 
Infection and Immunity, July 2002, p. 3978-3981, Vol. 70, No. 7
0019-9567/02/$04.00+0 DOI: 10.1128/IAI.70.7.3978-3981.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Genome Sequence Survey Identifies Unique Sequences and Key Virulence Genes with Unusual Rates of Amino Acid Substitution in Bovine Staphylococcus aureus
Lisa L. Herron,1 Rajit Chakravarty,1,
Christopher Dwan,2 J. Ross Fitzgerald,3 James M. Musser,3 Ernest Retzel,2 and Vivek Kapur1*
Department of Veterinary Pathobiology and Biomedical Genomics Center, University of Minnesota, St. Paul, Minnesota 55108,1
Computational Biology Center, University of Minnesota, Minneapolis, Minnesota 55455,2
Laboratory of Human Bacterial Pathogenesis, Rocky Mountain Laboratories, National Institute of Allergy and Infectious Disease, National Institutes of Health, Hamilton, Montana 598403
Received 28 March 2002/
Accepted 3 April 2002

ABSTRACT
Staphylococcus aureus is a major cause of mastitis in bovine
and other ruminant species. We here present the results of a
comparative genomic analysis between a bovine mastitis-associated
clone, RF122, and the recently sequenced human-associated clones,
Mu50 and N315, of
Staphylococcus aureus. A shotgun sequence
survey of

10% of the RF122 genome identified numerous unique
sequences and those with elevated rates of nonsynonymous substitution.
Taken together, these analyses show that there are notable differences
in the genomes of bovine mastitis-associated and human clones
of
S. aureus and provide a framework for the identification
of specific factors associated with host specificity in this
major human and animal pathogen.

TEXT
Staphylococcus aureus is a widely distributed animal and human
pathogen. Phenotypic analyses of
S. aureus strains isolated
from human and bovine hosts have revealed important differences
between the two populations of isolates, but the underlying
genetic basis for the variation among
S. aureus strains remains
unknown (
10,
14,
15). Multilocus enzyme electrophoretic analysis
of 2,077
S. aureus isolates from bovine, ovine, and human hosts
showed that only 6 of the 33 clusters included isolates from
more than one host species (
11). Population genetic analysis
of 357
S. aureus isolates recovered from cases of bovine mastitis
from global sources suggested that only eight clones represented
90% of the isolates, and these clones rarely overlapped with
those responsible for human disease (
5). To begin identifying
the nature and scope of the genetic differences between the
bovine and human staphylococcal genomes, we present here a comparative
analysis of the genome content of the common bovine
S. aureus clone RF122 (
4) and the recently sequenced human pathogenic
clones Mu50 and N315 (
6).
A random shotgun genomic library was prepared as described previously (7), with the exception that RF122 genomic DNA was mechanically sheared with a Hydroshear device (GeneMachines, San Carlos, Calif.). The first 850 random shotgun clones of RF122, encoding 285,077 unique bases and representing approximately 10% of the total genome, were sequenced by using dideoxy terminator chemistry on ABI 3700 DNA sequencers (ABI, Foster City, Calif.) and assembled into 600 contigs consisting of one to six sequences and ranging in size from 82 to 1,828 bases using phredPhrap (P. Green, University of Washington, Seattle) (details of individual genes, gene names, and categories are available at http://www.cbc.umn.edu/ResearchProjects/AGAC/Staph/staphhome.html). The average GC content of the initial sequences was 33.95%, closely approximating the GC content of the completed S. aureus genomes (6). The range of GC contents of individual contigs was 21 to 73%. Fourteen contigs had GC contents greater than 45%. Of these, two had no significant hits (P < e-15) to any GenBank submission based on tBLASTn searches, five had significant homology to S. aureus genomic DNA, and one was homologous to an S. aureus plasmid (1). The remaining five contigs had significant (P < e-15) homologies to sequences of organisms other than staphylococci and may reflect recent horizontal transfer events.
Identification of unique sequences in RF122.
Twelve sequences without homology to any existing sequence in the comprehensive microbial databases were identified in RF122. An additional seven contigs with no homology to S. aureus but with close matches to genes from other organisms were identified, including homologs to a Shigella transposon, a pathogenicity island element from Escherichia coli, a Rhodopseudomonas palustris ferredoxin transporter, and a hypothetical protein from Mesorhizobium loti. Other orthologs matched the gram-positive organisms Streptococcus agalactiae (sodA), Streptococcus pyogenes (putative phage terminase), Enterococcus faecalis (putative invertase), and Bacillus pumilus (Tn 10 transposon). Since the present analysis represents sequences from only
10% of the RF122 genome, these data suggest that there are likely to be 100 or more sequences that differentiate the bovine S. aureus clone from its human counterparts.
Recent analysis of RF122 led to the identification of a novel putative pathogenicity island, SaPIbov (4). The 15,981-bp island encodes several virulence factors and toxins (4). As expected, our RF122 contigs included several SaPIbov components. Two of the homologous contigs matched no sequences in either of the finished human S. aureus genomes, indicating possible bovine-specific components of SaPIbov. The genes in these contigs include homologs to sel, whose hypothetical protein product is similar to staphylococcal enterotoxin I, and to a phage terminase homolog.
We were able to assign putative functions to many of the gene products based on motifs identified in the homologs from the finished Mu50 and N315 S. aureus genomes and the genome of Bacillus subtilis (Fig. 1). In total, 28% of the identified homologs encoded proteins of unknown function. Other genes of interest include bap, a biofilm-related gene thus far identified only in bovine isolates, drp35, an antibiotic resistance determinant, and rot, a transcriptional repressor regulating toxin expression previously characterized only in human strains (2, 8, 9).
Patterns of nucleotide substitution.
To understand the level and nature of nucleotide variation between
RF122 and the recently sequenced human clones of
S. aureus,
amino acid substitution rates for the aligned sequences were
calculated according to the algorithm of Nei and Gojobori (
12).
For genes with integral roles in bacterial survival and fitness,
the rate of nonsynonymous substitution (dN) is often lower than
the rate of synonymous substitution (dS) because these genes
are under functional constraint and the majority of phenotypic
changes will be disadvantageous. The majority of dN:dS ratios
of the RF122
S. aureus sequences characterized thus far in comparison
with Mu50 and N315 fall in the range of 1:10 to 1:4, corresponding
to the most highly conserved sequences (Fig.
2A). The low average
ratio of dN:dS for RF122 compared to Mu50 and N315 indicates
that the majority of
S. aureus genes are under functional constraint
(Fig.
2B). Several of the highly conserved genes encoded virulence
factors, including clumping factor B, as well as metabolic genes
from energy, lipid, nucleic acid, and amino acid pathways (Table
1). The analysis also identified 19 orthologs with dN:dS ratios
greater than 1, indicating that mutations that result in a change
in the primary structure of the proteins have been conserved
in these sequences. While the majority of the protein products
from these genes were hypothetical, a urease subunit, an ammonium
transporter, the cell division protein FtsH, the surface protein
and antibiotic resistance determinant Drp35, and a leukotoxin
subunit were also included in this category. The leukotoxin
subunit is of particular interest in the context of host specificity,
as leukotoxins interact directly with host erythrocytes and
phagocytes to cause potentially catastrophic inflammation. It
has been suggested that leukotoxins may play a role in causing
clinical mastitis, and leukotoxin genes are present in all bovine
strains tested, but a causal relationship between leukotoxin
production and mastitis remains to be established (
3,
13). Should
leukotoxins be a significant factor in causing mastitis, positive
selection of phenotype-altering mutations in the F subunit of
S. aureus leukotoxin may be related to differences in the bovine
and human immunological environments.
Concluding comments.
In sum, this preliminary analysis of the genome of the common
clone of bovine
S. aureus reveals the presence of numerous genes
and sequences that differentiate this isolate from previously
characterized strains of the species. Our investigation also
provides strong evidence that the whole genome sequencing approach
will build a solid foundation for future investigation of the
mechanisms of virulence and host specificity of this important
human and animal pathogen.

ACKNOWLEDGMENTS
We thank Barbara May, Michael Paustian, and Ling-Ling Li for
their assistance with the preparation and analysis of the library
and sequences. We also appreciate the support of Megan Lillehei
and Todd Markovitz at the Advanced Genetic Analysis Center at
the University of Minnesota for assistance with automated DNA
sequencing.
This research is funded by a competitive award from the USDA-NRI Sustaining Animal Health and Well-Being Program (to V.K.). L.L.H. is supported by NIH NIGMS Biological Process Technology Institute fellowship GM08347.

FOOTNOTES
* Corresponding author. Mailing address: Department of Veterinary Pathobiology, University of Minnesota, 205 Veterinary Science, 1971 Commonwealth Ave., St. Paul, MN 55108. Phone: (612) 625-7712. Fax: (612) 625-5203. E-mail:
vkapur{at}umn.edu.

Editor: V. J. DiRita
Present address: Biomedical Engineering Center, University of Wisconsin, Madison, WI 53706. 

REFERENCES
1
- Altschul, S. F., T. L. Madden, A. A. Schaffer, J. Zhang, Z. Zhang, W. Miller, and D. J. Lipman. 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25:3389-3402.[Abstract/Free Full Text]
2
- Cucarella, C., C. Solano, J. Valle, B. Amorena, I. Lasa, and J. R. Penades. 2001. Bap, a Staphylococcus aureus surface protein involved in biofilm formation. J. Bacteriol. 183:2888-2896.[Abstract/Free Full Text]
3
- Fitzgerald, J. R., W. J. Meaney, P. J. Hartigan, and C. J. Smyth. 2000. Molecular population and virulence factor analysis of Staphylococcus aureus from bovine intramammary infection. J. Appl. Microbiol. 88:1028-1037.[CrossRef][Medline]
4
- Fitzgerald, J. R., S. R. Monday, T. J. Foster, G. A. Bohach, P. J. Hartigan, W. J. Meaney, and C. J. Smyth. 2001. Characterization of a putative pathogenicity island from bovine Staphylococcus aureus encoding multiple superantigens. J. Bacteriol. 183:63-70.[Abstract/Free Full Text]
5
- Kapur, V., W. M. Sischo, R. S. Greer, T. S. Whittam, and J. M. Musser. 1995. Molecular population genetic analysis of Staphylococcus aureus recovered from cows. J. Clin. Microbiol. 33:376-380.[Abstract]
6
- Kuroda, M., T. Ohta, I. Uchiyama, T. Baba, H. Yuzawa, I. Kobayashi, L. Cui, A. Oguchi, K. Aoki, Y. Nagai, J. Lian, T. Ito, M. Kanamori, H. Matsumaru, A. Maruyama, H. Murakami, A. Hosoyama, Y. Mizutani-Ui, N. K. Takahashi, T. Sawano, R. Inoue, C. Kaito, K. Sekimizu, H. Hirakawa, S. Kuhara, S. Goto, J. Yabuzaki, M. Kanehisa, A. Yamashita, K. Oshima, K. Furuya, C. Yoshino, T. Shiba, M. Hattori, N. Ogasawara, H. Hayashi, and K. Hiramatsu. 2001. Whole genome sequencing of meticillin-resistant Staphylococcus aureus. Lancet 357:1225-1240.[CrossRef][Medline]
7
- May, B. J., Q. Zhang, L-L. Li, M. L. Paustian, T. S. Whittam, and V. Kapur. 2000. Complete genomic sequence of Pasteurella multocida, Pm70. Proc. Natl. Acad. Sci. USA 98:3460-3465.[Abstract/Free Full Text]
8
- McNamara, P. J., K. C. Milligan-Monroe, S. Khalili, and R. A. Proctor. 2000. Identification, cloning, and initial characterization of rot, a locus encoding a regulator of virulence factor expression in Staphylococcus aureus. J. Bacteriol. 182:3197-3203.[Abstract/Free Full Text]
9
- Murakami, H., H. Matsumaru, M. Kanamori, H. Hayashi, and T. Ohta. 1999. Cell wall affecting antibiotics induce expression of a novel gene, drp35, in Staphylococcus aureus. Biochem. Biophys. Res. Commun. 264:348-351.[CrossRef][Medline]
10
- Musser, J. M., P. M. Schlievert, A. W. Chow, P. Ewan, B. N. Kreiswirth, V. T. Rosdahl, A. S. Naidu, W. Witte, and R. K. Selander. 1990. A single clone of Staphylococcus aureus causes the majority of cases of toxic shock syndrome. Proc. Natl. Acad. Sci. USA 87:225-229.[Abstract/Free Full Text]
11
- Musser, J. M., and R. K. Selander. 1990. Genetic analysis of natural populations of Staphylococcus aureus, p 59-67. In R. Novick and R. A. Skurray (ed.), Molecular biology of the staphylococci. VCH Publishers, New York, N.Y.
12
- Nei, M., and T. Gojobori. 1986. Simple methods for estimating the numbers of synonymous and nonsynonymous nucleotide substitutions. Mol. Biol. Evol. 3:418-426.[Abstract]
13
- Schuberth, H., C. Krueger, H. Zerbe, E. Bleckmann, and W. Leibold. 2001. Characterization of leukocytotoxic and superantigen-like factors produced by Staphylococcus aureus isolates from milk of cows with mastitis. Vet. Microbiol. 82:187-199.[CrossRef][Medline]
14
- Sompolinsky, D., Z. Samra, W. W. Karakawa, W. F. Vann, R. Schneerson, and Z. Malik. 1985. Encapsulation and capsular types in isolates of Staphylococcus aureus from different sources and relationship to phage types. J. Clin. Microbiol. 22:828-834.[Abstract/Free Full Text]
15
- Sordelli, D. O., F. R. Buzzola, M. I. Gomez, L. Steele-Moore, D. Berg, E. Gentilini, M. Catalano, A. J. Reitz, T. Tollersrud, G. Denamiel, P. Jeric, and J. C. Lee. 2000. Capsule expression by bovine isolates of Staphylococcus aureus from Argentina: genetic and epidemiologic analyses. J. Clin. Microbiol. 38:846-850.[Abstract/Free Full Text]
Infection and Immunity, July 2002, p. 3978-3981, Vol. 70, No. 7
0019-9567/02/$04.00+0 DOI: 10.1128/IAI.70.7.3978-3981.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
This article has been cited by other articles:
-
Lee, D.-S., Burd, H., Liu, J., Almaas, E., Wiest, O., Barabasi, A.-L., Oltvai, Z. N., Kapatral, V.
(2009). Comparative Genome-Scale Metabolic Reconstruction and Flux Balance Analysis of Multiple Staphylococcus aureus Genomes Identify Novel Antimicrobial Drug Targets. J. Bacteriol.
191: 4015-4024
[Abstract]
[Full Text]
-
Baba, T., Bae, T., Schneewind, O., Takeuchi, F., Hiramatsu, K.
(2008). Genome Sequence of Staphylococcus aureus Strain Newman and Comparative Analysis of Staphylococcal Genomes: Polymorphism and Evolution of Two Major Pathogenicity Islands. J. Bacteriol.
190: 300-310
[Abstract]
[Full Text]
-
Tuchscherr, L. P. N., Gomez, M. I., Buzzola, F. R., Calvinho, L. F., Lee, J. C., Sordelli, D. O.
(2007). Characterization of a New Variant of IS257 That Has Displaced the Capsule Genes within Bovine Isolates of Staphylococcus aureus. Infect. Immun.
75: 5483-5488
[Abstract]
[Full Text]
-
Smith, E. M., Green, L. E., Medley, G. F., Bird, H. E., Dowson, C. G.
(2005). Multilocus Sequence Typing of Staphylococcus aureus Isolated from High-Somatic-Cell-Count Cows and the Environment of an Organic Dairy Farm in the United Kingdom. J. Clin. Microbiol.
43: 4731-4736
[Abstract]
[Full Text]
-
Smith, E. M., Green, L. E., Medley, G. F., Bird, H. E., Fox, L. K., Schukken, Y. H., Kruze, J. V., Bradley, A. J., Zadoks, R. N., Dowson, C. G.
(2005). Multilocus Sequence Typing of Intercontinental Bovine Staphylococcus aureus Isolates. J. Clin. Microbiol.
43: 4737-4743
[Abstract]
[Full Text]
-
van Leeuwen, W. B., Melles, D. C., Alaidan, A., Al-Ahdal, M., Boelens, H. A. M., Snijders, S. V., Wertheim, H., van Duijkeren, E., Peeters, J. K., van der Spek, P. J., Gorkink, R., Simons, G., Verbrugh, H. A., van Belkum, A.
(2005). Host- and Tissue-Specific Pathogenic Traits of Staphylococcus aureus. J. Bacteriol.
187: 4584-4591
[Abstract]
[Full Text]
-
Tormo, M. A., Knecht, E., Gotz, F., Lasa, I., Penades, J. R.
(2005). Bap-dependent biofilm formation by pathogenic species of Staphylococcus: evidence of horizontal gene transfer?. Microbiology
151: 2465-2475
[Abstract]
[Full Text]
-
Cucarella, C., Tormo, M. A., Ubeda, C., Trotonda, M. P., Monzon, M., Peris, C., Amorena, B., Lasa, I., Penades, J. R.
(2004). Role of Biofilm-Associated Protein Bap in the Pathogenesis of Bovine Staphylococcus aureus. Infect. Immun.
72: 2177-2185
[Abstract]
[Full Text]
-
Gilot, P., van Leeuwen, W.
(2004). Comparative Analysis of agr Locus Diversification and Overall Genetic Variability among Bovine and Human Staphylococcus aureus Isolates. J. Clin. Microbiol.
42: 1265-1269
[Abstract]
[Full Text]