| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2008, p. 114-124, Vol. 74, No. 1
0099-2240/08/$08.00+0 doi:10.1128/AEM.01373-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.

Department of Biotechnology and Environmental Biology, RMIT University, Bundoora West Campus, Bundoora, Melbourne, Victoria 3083, Australia,1 University of Melbourne, Rural Veterinary Unit, Maffra, Victoria 3860, Australia,2 Institute for Animal Health, Compton Laboratory, Compton RG20 7NN, United Kingdom,3 School of Veterinary Medicine and Science, University of Nottingham, Sutton Bonington, Leicestershire LE12 5RD, United Kingdom4
Received 21 June 2007/ Accepted 1 November 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
Although the introduction of mastitis control schemes has been successful in reducing the incidence of contagious mastitis caused by Streptococcus agalactiae and Staphylococcus aureus in dairy herds (32), these measures have had little impact on environmental pathogens, notably Streptococcus uberis and Escherichia coli (3, 4, 26, 29). In an earlier study, the incidence of clinical mastitis in Australian herds caused by S. uberis was estimated at 22.7% of cases (58). The ability to control these infections depends on a detailed knowledge of the epidemiology of the organism and the management of the herd and its environment.
Several molecular typing methods have been used to investigate the epidemiology of S. uberis. These methods include restriction fragment length polymorphism (23, 27), pulsed-field gel electrophoresis (PFGE) (9, 40, 57), randomly amplified polymorphic DNA fingerprinting (61), and multilocus sequence typing (MLST) (6, 42, 62). These studies suggested that transmission of S. uberis occurs from environmental sources (40, 61, 63) and that feces could be a reservoir (61, 63). An earlier MLST study identified specific sequence types (STs) belonging to major clonal complexes (CCs) from milk samples and the cow's environment and suggested that these STs have the ability to survive in the environment and to establish intramammary infections (42).
Several investigators have used PFGE to demonstrate that S. uberis is a highly diverse species (9, 36, 40, 57) with many different PFGE types present on a single farm, suggesting that the species is behaving as an opportunistic pathogen. However, chronically infected cows often harbor the same PFGE type (36, 40, 41, 57), and there has been a single report demonstrating (61) and others suggesting cow-to-cow transmission of a single clone (1, 40). These observations suggest the possibility that some strains are either hypervirulent, hypertransmissible between cows, or able to survive in a host, for example, by evading the host immune response. It has not been conclusively proven, however, whether certain clones with enhanced virulence are responsible for mastitis. In addition, despite numerous epidemiological studies using PFGE, no clear evidence with regard to the relatedness of S. uberis isolates from different herds has been reported in the literature. This lack of evidence can be attributed partly to the inherent difficulties associated with the interlaboratory comparison of PFGE-based epidemiological studies.
MLST differs fundamentally from PFGE and most other molecular typing methods, being based on nucleotide sequence data from approximately 500 bp of housekeeping genes that have been shown to accumulate sequence variation slowly. Accordingly, MLST, in contrast to PFGE, which uses whole genomic DNA, is much less affected by recent rearrangement of the genome by recombination. Consequently, MLST has greater utility for determining the recent ancestral lineage and the relatedness of individual strains. In addition, MLST utilizes multiple genes of different sizes for analysis, which provides substantial discriminatory power for subtyping. In some species, inclusion of one or more virulence-associated genes (multi-virulence-locus sequence typing) can be used to further increase the discriminatory power (33, 34, 54, 64). Both MLST and multi-virulence-locus sequence typing profiles are unambiguous and can be represented by a number of digits corresponding to the allelic number of each of the loci used. This greatly facilitates interlaboratory comparisons and the study of global epidemiology (12).
Two MLST methods have been developed for the typing of S. uberis. The initial method (62) used the pauA gene, which encodes a virulence factor (plasminogen activator A), and the gapC gene, which encodes a vaccine target (glyceraldehyde-3-phosphate dehydrogenase) (18). The pauA gene has been shown to exhibit a high ratio of nonsynonymous substitutions to synonymous substitutions (dN/dS) (>1.0) (62), indicating positive selection, and therefore may not provide accurate phylogenetic information (6). In addition, the initial MLST scheme uses only six loci (62), and as such, it may not provide sufficient resolution for a large data set compared with the seven loci used in the more recent MLST scheme (6). For these reasons, we chose to use the MLST scheme developed by Coffey et al. in 2006 (6).
In the present study, we performed MLST analyses on a collection of Australian S. uberis isolates from cases of clinical and subclinical mastitis and from cows with low somatic cell counts in their milk. The main aim of this study was to investigate the molecular epidemiology of S. uberis mastitis in Australia to ascertain whether there was any association of a particular MLST with mastitis. In addition, we investigated the evolutionary lineage of Australian isolates in relation to the S. uberis global population as a whole by comparing the MLST profiles of Australian isolates with those currently available in the MLST database from different countries.
| MATERIALS AND METHODS |
|---|
|
|
|---|
The geographic distribution and date of collection of the samples were as follows: 29 isolates were collected from 24 cows (28 quarters) from a dairy farm in Newry (eastern Victoria, Australia). Of these, 16 isolates were collected from a cow with clinical mastitis (n = 1), cows with subclinical mastitis (n = 5), and cows with low somatic cell counts in their milk (n = 10) in October 2004, and 13 isolates from cows with clinical mastitis were collected between November and December 2005. Twelve isolates from 12 quarters of seven cows in Timboon (western Victoria, Australia) were collected September 2004. Five isolates were obtained from a centralized veterinary diagnostic laboratory (Gippsland Vetnostics, Traralgon, Victoria, Australia) between April and May 2006. Consequently, the latter samples were originally isolated from quarter milk samples (n = 5) from dairy cattle (n = 5) from various parts of Victoria, Australia.
Culture of milk sample and identification of S. uberis.
Quarter milk samples were collected according to Australian guidelines (5). Samples were kept at 4 °C during transportation and cultured immediately on arrival at the laboratory. Milk was cultured by spreading 50 µl over the surface of one horse blood agar plate and one Edward's agar plate (Oxoid Ltd., Basingstoke, United Kingdom). Plates were incubated at 37°C for 48 h.
All isolates were initially identified as S. uberis on the basis of the following: Gram stain, catalase reaction, hydrolysis of esculin and hippurate, or API 20 Strep (bioMérieux Vitek, Hazelwood, MO). In addition, S. uberis species-specific PCRs (20) for both the 16S rRNA and 23S rRNA genes were performed for primary species identification. The PCR conditions and the primers used are listed in Table 1. To further confirm species identification, 16S rRNA gene sequencing was also performed as described previously (42). The 16S rRNA gene sequences were compared with the reference sequence of S. uberis, strain HN1 (AB023576.1) using BLASTN at the NCBI (http://www.ncbi.nlm.nih.gov/BLAST/), and the individual isolates were designated as S. uberis if the 16S rRNA sequence was 100% identical to the reference sequence.
|
The GelCompar II software (Applied Maths, Kortrijk, Belgium) was used to perform cluster analysis and construct a dendrogram using the unweighted pair group method using arithmetic averages (UPGMA) based on Dice coefficients using 0.75% optimization and 1% tolerance parameters. For analysis, PFGE patterns showing <80% similarity were assigned to different types, whereas patterns with similarity values of >80% and <100% were assigned to subtypes.
MLST analysis.
MLST was performed using the primers described by Coffey et al. in 2006 (6) and the modified cycling conditions as detailed in Table 1. The S. uberis reference strain, ATCC 700407, was used as a positive control. The allelic profiles and sequence types (STs) of the individual isolates were determined using the query functions implemented in the S. uberis MLST database (http://pubmlst.org/suberis). The clonal complex for each ST in the MLST database was assigned by the curator and in this study referred to as global clonal complex (GCC).
Data analysis.
Simpson's index of diversity (D) with 95% confidence intervals (95% CI) of the MLST and the PFGE techniques used in this study were determined as described previously (19, 24).
The dN/dS ratios, the G+C ratios, and the number and proportion of variable nucleotide sites for each of the seven loci were analyzed using Sequence Type Analysis and Recombinational Tests, version 2 (START2) (28).
Lineage analyses of Australian isolates were performed by constructing a UPGMA-based dendrogram using START2 (28). The eBURST version 3 program (http://eburst.mlst.net/) was used to investigate the population structure of the Australian S. uberis isolates and the relatedness of the Australian isolates in the global population of S. uberis. For eBURST analyses, an individual CC was defined as a group of isolates that shared alleles at six of the seven loci (12, 14, 49).
The standardized index of association (ISA) (21) was calculated using START2. ISA measures the degree of linkage equilibrium between allelic profiles and estimates the rate of recombination (30, 48, 50). The rates of recombination to mutation and the per site ratio of recombination to mutation parameter were determined as described previously (16). Statistical analyses were performed using chi-square tests using SPSS 13.0 software (SPSS Inc., Chicago, IL).
Investigation of yqiL-negative Australian isolates.
In an attempt to amplify the yqiL gene of the four Australian isolates which were PCR negative using the original primers (6), two additional primer pairs for the yqiL gene were designed using Primer3 (44). The primer specifications and amplification conditions used are detailed in Table 1.
To confirm the absence of the yqiL gene, Southern blot hybridization was performed. For Southern blot hybridization, 1 µg of genomic DNA was digested with HindIII (Promega, Madison, WI) and separated by agarose gel electrophoresis at 80 V for 90 min in a 1.5% (wt/vol) agarose gel. The DNA fragments were subsequently transferred by Southern blotting to a positively charged nylon membrane (Immobilon-Ny+; Amersham, Little Chalfont, United Kingdom) and hybridized with a yqiL gene probe using standard procedures. The yqiL gene probe was prepared using the digoxigenin DNA labeling kit (Roche Diagnostics, Basel, Switzerland) in accordance with the manufacturer's instructions with a yqiL PCR product generated from a reference yqiL-positive S. uberis isolate used as a positive control.
The flanking region of the yqiL gene was further investigated to determine whether a possible virulence factor gene or surface protein gene was present. Any such association would place the locus under high selective pressure making it unsuitable for MLST. The yqiL flanking region of the S. uberis genome sequence obtained from the S. uberis genome project (http://www.sanger.ac.uk./Projects/S_uberis/) was used for this analysis.
The S. uberis yqiL gene was queried against the S. uberis genome sequence using BLASTN and BLASTX to ascertain whether possible yqiL homologues that could possibly complement for yqiL-negative isolates were present.
Detection of hasA and hasC genes by PCR and Southern blot hybridization.
Amplification of the capsular genes, hasA and hasC, was performed using a modification of a previously described method (17). The PCR conditions and the primers used are listed in Table 1. The isolates that were negative for hasA by PCR and the strains that were positive for hasA by PCR (control strains) were further examined by Southern blot hybridization using a digoxigenin-labeled hasA gene probe generated from a reference S. uberis strain using the same procedure as that used for the yqiL probe described earlier.
Nucleotide sequence accession numbers.
The nucleotide sequences of the MLST alleles from the Australian S. uberis isolates have been deposited in GenBank as a data set with the accession numbers EF672733 to EF672747.
| RESULTS |
|---|
|
|
|---|
Statistical analysis based on the geographic origin of the samples against types of global clonal complex demonstrated that while more than 50% of our samples were isolated from farms in Newry, Australia, there was no significant difference (P = 0.57) in the distribution of strains by geographic location of sampling against types of GCC.
PFGE analysis.
Analysis of the Australian S. uberis isolates by PFGE clearly demonstrated that they exhibited considerable diversity. In total, 42 PFGE types (including subtypes) were identified from the 46 isolates (Table 2). There were four examples of identical PFGE types being isolated from different cows on the same farm (three cases) or from different quarters of the same cow (one case). In no instance were identical PFGE types observed for cows from different farms.
|
The most common GCC was GCC ST143 (26%), followed by GCC ST5 (13%) and GCC ST86 (9%); 24 of the isolates (52%) were not assigned to a specific GCC (Table 2). Thirty-one out of 33 STs, 4 yqiL-negative isolates, and 15 of the alleles had not been previously recorded in the MLST database and were identified in Australia for the first time.
The STs with the largest number of isolates were ST60 and ST155, with both consisting of four isolates (8.7%) (Table 2). All four ST60 isolates had different PFGE types (including subtypes), whereas three ST155 isolates had an identical PFGE type (PFGE type 8a), and the fourth was a different subtype (type 8b) (Table 2). The ST of four Australian ST60 isolates was identical to the ST of isolates originally identified in the United Kingdom (41). One ST184 isolate had a MLST profile identical to that of an isolate from New Zealand (42). In addition, two ST194 isolates recovered from different farms located 100 km apart in Victoria, Australia, were shown by PFGE analysis not to be clones.
(ii) Analysis of the loci of Australian isolates.
Analysis of the STs (n = 33) identified in the Australian data set (excluding the four yqiL-negative isolates, as START2 cannot perform locus analysis using partial allelic profiles) showed the dN/dS ratios for each locus were substantially <1.0, ranging from 0.000 (gki) to 0.139 (tpi). The G+C ratios for each locus ranged from 35.7% (tdk) to 45.6% (recP). The average number of alleles at each locus was 7.1, and the proportion of variable nucleotide sites in each locus ranged from 1.1% (recP) to 4.4% (tdk). The dN/dS ratios, the G+C ratios, and the proportions of variable nucleotide sites were comparable to those in a previous study (6).
(iii) Lineage analysis using UPGMA and eBURST.
The majority of the Australian isolates (85%) were resolved (<0.3 linkage distance) into three clusters (Fig. 1). The level of 0.3 linkage distance was arbitrarily set. This is not an absolute limit, but it distinguished two main clusters and an association between GCC and disease status. Cluster I (n = 31) consisted of two-thirds of the Australian isolates, and 90% of the isolates in this cluster were from cows with either clinical or subclinical mastitis. There were two possible clonal complexes in cluster I identified by eBURST analysis, and these were termed local clonal complex ST60 (LCC ST60) and ST277 (LCC ST277) (Fig. 2). ST60 was the founder of LCC ST60, which consisted of seven isolates, four single-locus variants (SLVs) and two double-locus variants, while LCC ST277 consisted of three isolates and two SLVs (Fig. 2). The majority of STs (80%) belonging to the two LCCs belonged to GCC ST5 and GCC ST143 (Fig. 1).
|
|
All six isolates in the GCC ST5 and approximately 90% of the isolates of GCC ST143 were associated with clinical and subclinical mastitis (Table 3). In contrast, none of the isolates of GCC ST86 were associated with clinical mastitis, and 75% of these isolates were associated with cows with low somatic cell counts. Chi-square analysis of the association of disease status and GCC found that it was statistically significant (P = 0.006) (Table 3).
|
|
(ii) Ratio of recombinational to mutational events and the per site ratio of recombination to mutation parameter.
All of the MLST STs available as of February 2007 in the S. uberis MLST database (n = 593) were used for this analysis. A total of 9 possible CCs and 60 SLVs were identified following BURST analysis (Fig. 4). While 59 of the SLVs arose by recombination, surprisingly, only one variant yqiL allele arose by point mutation (rate of recombinational to mutational events of 59:1) (Fig. 4; Table 4) (16). While 21 variant alleles differed at a single nucleotide position from corresponding alleles in founders of the CCs, 20 of these were assigned as arising by recombinational events, since the variant alleles were present in unrelated STs, which differed at three or more loci, in the S. uberis MLST database (Fig. 4; Table 4) (16).
|
|
Investigation of yqiL-negative isolates.
All of the 46 Australian isolates yielded a PCR product for six of the seven loci following PCR (6). Four isolates (686-5, 686-10, 3064-1, and 3217-3) failed to amplify the yqiL gene with the original primers (6) and with the new primers designed for this study. The new yqiL primers successfully amplified yqiL from the control S. uberis. Southern blot analysis confirmed that these isolates lacked the yqiL gene (data not shown). BLASTN and BLASTX analyses of the yqiL gene against the S. uberis genome sequence in the Sanger Institute database identified a single yqiL homologue that exhibited 38% amino acid (57/148) and 57% nucleotide sequence identity (171/299) with the S. uberis yqiL allele 1. This open reading frame (ORF) was located in the region from 1795067 bp to 1796317 bp of the genome of S. uberis strain 0140J, 400 kb distant from yqiL. On further analysis using BLASTP, the yqiL ORF homologue exhibited extensive identity with both the thiolase B gene from Oceanobacillus iheyensis (46% amino acid identity [186/404]) (52) and the acetyl coenzyme A (acetyl-CoA) acetyltransferase gene from Clostridium difficile (46% amino acid identity [186/402]) (46). A search of the Conserved Domain Database at NCBI (35) showed that the S. uberis yqiL ORF homologue and both the thiolase B and the acetyl-CoA acetyltransferase genes contained several conserved domains characteristic of enzymes involved in the synthesis of acetoacyl-CoA (cd00751) (25, 60).
Using ORF Finder (NCBI, http://www.ncbi.nlm.nih.gov/gorf/gorf.html) and BLASTX (NCBI, http://www.ncbi.nlm.nih.gov/BLAST), a penicillin binding protein was found 15 kb upstream of the yqiL gene.
hasA and hasC PCR.
While all of the 46 Australian S. uberis isolates possessed the hasC gene, only 39 (85%) of the isolates possessed the hasA gene (Table 2). Three out of 10 hasA PCR-negative isolates and the three hasA PCR-positive control isolates were positive by Southern blotting using the hasA probe. Chi-square analysis of the association between disease status and the presence of hasA was statistically significant (P = 0.014) (Table 5). In addition, chi-square analysis of the association between GCCs and the presence of hasA was statistically significant (P = 0.005) (Table 6).
|
|
| DISCUSSION |
|---|
|
|
|---|
The data presented in this paper strongly suggest that the isolates in cluster I, in particular those belonging to GCC ST5 and GCC ST143, were highly associated with clinical and subclinical mastitis. The dendrogram derived following UPGMA for lineage analysis (Fig. 1) clearly distinguished a cluster of isolates associated with clinical and subclinical mastitis (cluster I) from a cluster of isolates from cows with low somatic cell counts in their milk (cluster III). Ninety percent of isolates belonging to cluster I were from cows with clinical and subclinical mastitis. The vast majority (94%) of the Australian isolates belonging to the GCC ST5 and GCC ST143 were distributed in cluster I, and there was a statistically significant association between GCC type and disease status (P = 0.006) (Table 3). In addition, eBURST analysis (Fig. 2) identified two possible LCCs in cluster I, suggesting the presence of genetically closely related strains in the cluster, supporting a previous study that demonstrated low levels of heterogeneity in isolates from clinical mastitis (27). The isolates belonging to GCC ST5 and GCC ST143 may possess virulence factors promoting invasion of host tissue, survival in the host environment or evasion of the host immune response, and internalization in the mammary gland cells (53). This hypothesis was further supported by hasA analysis. In this study, 95% of the isolates in GCC ST5 and GCC ST143 possessed hasA, whereas only 25% of the isolates in GCC ST86 possessed hasA gene with chi-square analysis demonstrating a significant link between the possession of hasA and the GCC type (P = 0.005) (Table 6). These results are in agreement with previous studies (6, 42), which suggested a link between both GCC ST5s and GCC ST143s with hasA. Although hasA is not thought to be directly involved in the pathogenicity of S. uberis (17), it nevertheless appears to be a virulence marker gene.
In contrast, cluster III, which was genetically distinct from cluster I (Fig. 1), consisted of isolates mainly from cows with low somatic cell counts and comprised all of the isolates belonging to GCC ST86. This finding was further corroborated following eBURST analysis (Fig. 3). These observations suggest that isolates in GCC ST86, although forming a clonal complex in the global S. uberis population, were not associated with mastitis. Accordingly, it may be possible that these isolates may have the ability to colonize in a host environment and evade the host immune system but lack pathogenicity, behaving like a commensal species. This hypothesis is further supported by the results of a recent study which demonstrated that most of the isolates in GCC ST86 were less associated with clinical mastitis (42).
There was an inconsistency in cluster III where one isolate (ST264) belonging to GCC ST143 grouped with GCC ST86s (Fig. 1). This could be due to the assigning system used by the MLST database, which differs from UPGMA which was used in Fig. 1 (6). ST264 shares four out of seven allelic profiles with both ST143 and ST86 and therefore could be assigned either to GCC ST143 or to GCC ST86 in the database; however, ST264 also shares five of seven profiles with its nearest neighbor (ST262) and therefore should possibly be assigned to GCC ST86 by the MLST database (Fig. 1).
Australian LCCs were distributed between GCC ST5 and GCC ST143 in the global population of S. uberis (Fig. 3), suggesting that these isolates were closely associated with both GCC ST5 and ST143. This is in contrast to a previous study which demonstrated that STs identified in United Kingdom and New Zealand were clearly separated into two GCC types, GCC ST5 and ST143 (42).
ST60, the founder of LCC ST60, was initially identified in the United Kingdom and has not been found in any other country. In a similar manner, ST184 was originally identified in New Zealand, and neither of these STs has been demonstrated outside their country of origin. Accordingly, the detection of both ST60 and ST184 in our Australian data set represents the first report of the occurrence of identical STs in different countries (New Zealand) and continents (Europe). These findings support the hypothesis that ancestors of Australian S. uberis are related to isolates in the United Kingdom and New Zealand. Specifically, the eBURST results (Fig. 3) suggest ST60 or ST184 could represent a recent ancestor of Australian isolates highly associated with mastitis and may represent key pathogenic lineages in the Australian S. uberis population.
It should be noted that not all STs belonging to the three major GCCs (Table 2) are represented in Fig. 3. Five STs (ST153, ST217, ST259, ST264, and ST272) do not appear in Fig. 3. These apparent inconsistencies were due to the assigning system used by the MLST database which differs from the eBURST analysis (6, 42). To avoid any possible confusion, the classification of GCC for each ST by the MLST database was used for all analyses.
The ISA analyses for both of the Australian data set and all isolates in the MLST database demonstrated linkage disequilibrium, suggesting a clonal population structure, but the low ISA also implied that substantial recombination had occurred in the S. uberis population. A comparison of the ISA between other species demonstrated that the ISA of S. uberis was lower than that observed in Neisseria meningitidis (ISA = 0.14) and Campylobacter jejuni (ISA = 0.256) which are known for their high rate of inter- and intraspecies recombination (13, 15, 30, 50). These results are in agreement with earlier studies using maximum likelihood tree analysis which demonstrated substantial recombination in the evolutionary history (6) and significant reticulate evolution detected between the loci used for MLST (62).
In addition, empirical recombination analysis demonstrated a substantial recombination rate among the entire MLST database including our Australian data set. The ratio of recombinational to mutational events (59:1) was at least five times higher than the rate reported for N. meningitidis and Streptococcus pneumoniae (16). This finding further supports the occurrence of substantial recombination in the S. uberis population.
In a previous study, S. uberis was shown to possess phages (22), and although transduction could be responsible for recombination, the exact mechanisms of recombination in S. uberis are still unknown. However, given the substantial recombination observed in both the Australian and global S. uberis populations, it is possible that frequent recombination of genomic DNA between isolates could lead to the genetic diversity as demonstrated by the MLST and PFGE results in this and previous studies (9, 40, 57). Consequently, frequent recombination would render phylogenetic analysis using PFGE data of very little significance unless the samples were very closely related. In addition, it has been demonstrated that frequent recombination within a bacterial population decreased the accuracy of eBURST analysis for identification of links between ancestors and descendants in a population, suggesting that additional analyses should be used in conjunction with eBURST to confirm the validity of the inferred relationships (56).
The yqiL locus is commonly used in MLST assays for various bacteria including Staphylococcus aureus (10), Staphylococcus epidermidis (55), Streptococcus pyogenes (11), Enterococcus faecalis (45), and Streptococcus uberis (6). Of the four isolates in our Australian collection shown to lack yqiL, two had identical 16S rRNA gene sequences to the reference sequence, whereas the other two yqiL-negative isolates differed only at 1 or 2 nt. Phylogenetic analysis demonstrated that they were only distantly related to Streptococcus parauberis and Streptococcus iniae with S. uberis being the closest species (data not shown). It is therefore highly unlikely that the Australian yqiL-negative isolates represent different species. This was further substantiated insofar as all of the other species identification tests, such as PCRs (20), biochemical tests, and API Strep 20 performed on the yqiL-negative isolates were consistent with S. uberis.
The investigation of the flanking region of the S. uberis yqiL gene revealed a penicillin binding protein 15 kb upstream of the yqiL locus. Since this possible surface protein was more than 6 kb distant from the yqiL gene, it was clear of the standard of locus selection (7) and would not be expected to affect the selection of the yqiL locus. While the mechanisms of yqiL deletions are unknown, it is possible that given the high rate of recombination observed in the S. uberis population in both this and previous studies (6), the deletion of the yqiL gene may have occurred by recombination. This is further substantiated by the occurrence of nonstandard yqiL alleles (yqiL alleles 4, 11, and 23 which possess a 1-nt deletion within the yqiL locus) in the S. uberis MLST database. Furthermore, strains of S. pyogenes lacking yqiL have been reported (38, 43) and recorded in the S. pyogenes MLST database (http://spyogenes.mlst.net/misc/info.asp) as nonstandard alleles.
The occurrence of yqiL-negative isolates remains enigmatic, and it is still unknown how the four yqiL-negative Australian strains could survive without what is considered to be a housekeeping gene. The identification of a yqiL homologue with homology to similar genes also involved in the synthesis of acetyl-CoA from diverse species of bacteria (25, 60) suggests that our four yqiL-negative isolates might utilize such a homologue to replace yqiL function. Since yqiL is absent from some strains of S. uberis, its status as a housekeeping gene is now in doubt.
In summary, MLST analysis of Australian isolates of S. uberis has clearly identified a cluster of specific STs highly associated with clinical and subclinical mastitis and a cluster of specific STs associated with cows with low somatic cell counts. Specifically, particular groups of STs, GCC ST5 and GCC ST143, were highly associated with clinical and subclinical mastitis. It is therefore postulated that such STs may represent a lineage of virulent S. uberis. In contrast, all of the isolates belonging to GCC ST86 were associated with low-cell-count cows and may represent a less virulent ST lineage.
The use of MLST has greatly facilitated interlaboratory comparison and has identified specific STs that were initially characterized in different countries and continents but are currently present in the Australian S. uberis population. Although four of the Australian isolates were shown to lack yqiL, the current MLST assay has been of significant utility and has identified possible highly virulent STs among the genetically diverse population of S. uberis both within Australia and in the global population as a whole.
The use of MLST will greatly facilitate the identification and characterization of pathogenic S. uberis isolates and specific STs that can be exploited for further research and vaccine development.
| ACKNOWLEDGMENTS |
|---|
We thank the staff at the Maffra Veterinary Centre; the staff at Gippsland Vetnostics, Traralgon, Victoria, Australia; veterinarian W. F. Morgan; the dairy farmers in Timboon and Newry in Australia for the provision of milk samples; and Geoff Hogg and Angelo Zaia, Microbiological Diagnostic Unit Public Health Laboratory, University of Melbourne, for access to sequencing facilities. We acknowledge the technical assistance of Wei Gao and Elizabeth Grabsch, Microbiology Department, Austin Hospital, and we thank Ben Fry, Flinders University, for consultation.
| FOOTNOTES |
|---|
Published ahead of print on 16 November 2007. ![]()
| REFERENCES |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||