Previous Article | Next Article ![]()
Applied and Environmental Microbiology, November 2007, p. 7199-7209, Vol. 73, No. 22
0099-2240/07/$08.00+0 doi:10.1128/AEM.01629-07
Copyright © 2007, American Society for Microbiology. All Rights Reserved.

INRA, Centre de Clermont-Ferrand Theix, Unité Microbiologie, Qualité et Sécurité des Aliments, 63122 Saint-Genès Champanelle, France
Received 17 July 2007/ Accepted 13 September 2007
|
|
|---|
|
|
|---|
Genomes, and consequently phenotypes, can be highly heterogeneous in bacterial strains, even if they belong to the same species. Genomic comparative studies between all Staphylococcus sequenced genomes showed an average 78% identity (35). The 22% genetic variability observed within the Staphylococcus genus is also observed within species. Indeed, genetic diversity between 36 S. aureus strains using whole-genome DNA microarrays showed that 78% of genes were common to all strains and 22% of genes were strain specific and might play a role in adaptation to specialized niches (7). In a previous study with genomic DNA restriction and pulsed-field gel electrophoresis (PFGE) analysis of seven S. xylosus strains of various origins, we showed genomic diversity within this species, with variations of chromosome size and ribosomal operon (rrn) number (5). S. xylosus isolated from opportunistic infections seemed to have chromosomes up to 11% larger than strains isolated from the food environment (5). Genomic differences were also found between S. xylosus strains isolated from meat products by PFGE and randomly amplified polymorphic DNA analyses (5), but much remains to be learned about the genetic content differences. Strain-specific genes of S. xylosus are of major interest in explaining the versatility of this species. We expanded comparisons of S. xylosus strains by using suppressive and subtractive hybridization (SSH). This method provides comprehensive surveys of genome differences among closely related strains by facilitating the identification of DNA segments present in one genome (tester) but absent in another (driver). It has been widely used for the analysis of bacterial pathogen genomes to discover new epidemiological markers, virulence factors, and host specificity determinants (38). SSH was also successfully used to study genome plasticity or intraspecies genomic diversity (23, 30).
In this work, SSH was carried out between three S. xylosus tester strains: S. xylosus S04002, used as a starter culture, S. xylosus S04009, isolated from cow mastitis, and S. xylosus 00-1747, responsible for mouse dermatitis (39). S. xylosus C2a, a commensal of human skin, was used as a driver. The chromosomal organization of the S. xylosus C2a strain has been reported (5) and has constituted the starting point for studying the genetic diversity within this species. Specific DNA fragments of 29.7, 25.4, and 22.8 kb were isolated from the three SSH libraries obtained from S. xylosus S04002, S04009, and 00-1747, respectively. Their distributions among 20 S. xylosus strains were analyzed, and their locations in the respective chromosomes of the S. xylosus tester strains were determined.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. S. xylosus strains used in this study and their origins
|
Chromosomal DNA and plasmid preparation.
Chromosomal DNA from S. xylosus was isolated by the method of Marmur (18). To isolate plasmid DNA from S. xylosus, cells were lysed with 16 µg/ml lysostaphin, and the QIAprep Miniprep kit (Qiagen, Courtaboeuf, France) was subsequently applied.
Suppressive-subtractive hybridizations and differential screening.
Three subtracted DNA libraries were generated following the recommendations for the Clontech PCR-Select bacterial genome subtraction kit (BD Biosciences Clontech, Erembodegem, Belgium). Genomic DNAs of S. xylosus strains S04002, S04009, and 00-1747 were used as the testers and that of S. xylosus C2a was used as a driver. To increase the proportion of chromosomal subtracted DNA, the endogenous plasmids of the tester strains were also used as drivers. Genomic and plasmid DNAs were digested by the appropriate RsaI enzyme. Suppression-subtractive hybridizations were carried out as recommended by the supplier except that hybridization temperatures were set at 61°C instead 63°C due to the low G+C content of Staphylococcus DNA. The pool of subtracted fragments was amplified by PCR using the Advantage cDNA polymerase (BD Biosciences Clontech). PCR products generated in each library were cloned into the pGEM-T Easy vector (Promega) according to the supplier's recommendations and transformed into Escherichia coli TG1 competent cells, which were then cultured on Luria-Bertani agar plates containing 100 µg/ml ampicillin, 50 µg/ml 5-bromo-4-chloro-3-indolyl-β-D-galactopyranoside, and 0.5 mM isopropyl thiogalactoside. Resulting transformants were individually grown in Luria-Bertani medium containing 100 µg/ml ampicillin in 96-well microtiter plates and stored in 20% glycerol at –70°C. Subtracted library clones were randomly selected as templates for PCR amplification using nested primers 1 and 2R, provided in the genome subtraction kit. PCR products were analyzed by agarose gel electrophoresis to verify the size of the inserts and the presence of one distinct amplified product. PCR products of each subtracted library were screened for tester-specific fragments using dot blot hybridization. They were arrayed in duplicate onto positively charged nylon membranes (Hybond N+; Amersham Biosciences, Orsay, France). One µg of RsaI-restricted genomic DNA from the driver and each tester DNA was randomly labeled with the Dig-High Prime system (Roche Applied Science, Neuilly sur Seine, France) and then used as probes to separately hybridize the blots. The digoxigenin-labeled probes were detected by the Dig color detection kit (Roche).
The distribution of tester-specific fragments was also evaluated among S. xylosus strains. Genomic DNAs from 20 strains (Table 1) were used as probes, as described above. The distribution of tester-specific fragments was analyzed by a complete linkage hierarchical cluster analysis using Microsoft Excel software with the XLStat add-in. A distance matrix was calculated by using Sokal and Sneath's (1) similarity coefficient to take into account positive (occurrence) and negative (nonoccurrence) matches.
DNA sequencing and analysis.
Plasmid DNA containing tester-specific inserts was isolated from individual clones using the NucleoSpin Plasmid QuickPure kit (Macherey-Nagel, Hoerdt, France). Sequencing reactions were performed using the BigDye terminator as recommended by the supplier (Applied Biosystems, Courtaboeuf, France) and analyzed by using an ABI 310 automated DNA sequencer. Bioinformatic analyses were performed from Web-based servers. The sequences were analyzed by BLASTX and/or BLASTN searches against the NCBI GenBank database. Based on the analysis, redundant sequences were identified and sequence assembly was done with the CAP3 program (http://bioinformatics.iastate.edu/aat/sas.html). ORF Finder and BLASTX analyses were performed on contigs. When no similarity was found, a PSI-BLAST analysis was performed.
Determination of chromosomal locations of subtracted sequences.
For Southern blotting, digested chromosomal DNAs of S. xylosus S04002, S04009, and 00-1747 strains separated by PFGE were passively transferred onto positively charged nylon membranes (Hybond-N+; Amersham Biosciences). Equal amounts of PCR amplicons of each tester-specific insert were pooled and randomly labeled using the Dig-High Prime system (Roche Applied Science). After hybridizations, detection was performed using the Dig chemiluminescence detection kit following the supplier's instructions (Roche Applied Science). To confirm the localization of tester-specific inserts, chromosomal DNA restriction fragments of interest were excised from the PFGE gel, rinsed three times for 30 min with 10 mM Tris-HCl, pH 8.5, and digested with RsaI (Promega) according to the supplier's recommendations. DNA fragments were thus extracted from agarose and purified before labeling with the AlkPhos direct labeling kit (Amersham Biosciences) and hybridized to the dot blot containing the subtracted unique sequence as described for the differential screening.
Nucleotide sequence accession numbers.
Insert sequences of the subtracted libraries have been deposited in GenBank and given accession numbers ER935832 to ER935953.
|
|
|---|
![]() View larger version (42K): [in a new window] |
FIG. 1. Unweighted pair group method analysis dendrogram with Dice coefficient of PFGE profiles of 24 SmaI macro-restricted S. xylosus chromosomes.
|
|
View this table: [in a new window] |
TABLE 2. S. xylosus S04002-, S04009-, and 00-1747-specific sequences with significant protein matches in the GenBank database
|
In the S. xylosus S04009 strain genome, fragments corresponding to five genes involved in myo-inositol metabolism as well as the Na+:myo-inositol cotransporter were subtracted (Table 2). Myo-inositol is a hexitol which is a constituent of virtually every living cell; it is synthesized by plants and is essential for humans and animals. Different bacterial species (Klebsiella spp., Salmonella spp., Pseudomonas spp., Bacillus subtilis, etc.) are able to use myo-inositol as carbohydrate source, but genes implicated in this catabolism have been described in few species (16). In Bacillus subtilis, two iol operons were described, including an operon of 10 genes, iolABCDEFGHIJ, and an operon of two genes, iolRS (40). Until now, these genes had never been described in a Staphylococcus genus.
Five ORFs of the S. xylosus 00-1747 strain were related to lysine and glycine metabolism genes, none of which corresponded to staphylococcal genes (Table 2). Sarcosine oxidase and betaine aldehyde dehydrogenase were shown to be implicated in metabolism of glycine betaine, an effective osmoprotectant (1, 4). Only betaine aldehyde dehydrogenase and dihydropicolinate synthase were identified in actual staphylococcal sequenced genomes, but these genes showed no significant match at the DNA level with these ORFs. Sarcosine oxidase and saccharopine dehydrogenase have never been observed in staphylococcal genomes. The fact that the S. xylosus 00-1747 strain contained numerous nonstaphylococcal homology ORFs might be explained by its origin. This strain was isolated from mouse dermatitis in Korea. The particular origin of this strain suggests that ecological and/or geographical divergences in S. xylosus species may exist.
From the S. xylosus meat starter S04002 subtraction, fragments corresponding to fermentative metabolism were found: an acetolactate synthase which transforms pyruvate in acetolactate, its regulator AlsR, and an acetoin dehydrogenase which transforms acetoin into diacetyl, as well as a lactate permease also implicated in anaerobic metabolism (Table 2). Acetoin and diacetyl are aroma compounds. Their production in fermented meat products has been characterized with a dairy product odor (20). Some S. xylosus strains produce these aroma compounds more than other staphylococcus species (34). This property is variable between S. xylosus strains, and a correlation between strain persistence in ripened meat and the capacity of strains to produce acetoin has been shown (9). In lactic acid bacteria isolated from fermented sausages, the production of acetoin is stimulated by a low pH and a low sugar concentration, conditions that could be found at the end of the ripening process (9). Recently, Fuchs et al. (8) showed that in S. aureus the production of acetoin and lactate was stimulated by anaerobic growth conditions.
The large number of differential genes predicted to be involved in metabolism strongly suggested a wide adaptation of S. xylosus strains to their environment and to the availability of the different substrates.
Horizontal transfer in S. xylosus S04009.
The S. xylosus S04009 strain isolated from cow mastitis showed 12 specific subtracted fragments corresponding to phage-related sequences (Table 2). The presence of at least three different bacteriophages, 77, 2638A, and EW, was deduced, and an IS605/IS200 transposase was also detected. Three phage-related proteins were found in S. xylosus S04002 and one in 00-1747 subtracted fragments. Phages in S. aureus are well studied: 37 phage genomes are reported in the GenBank database (3, 17). To our knowledge, there are few studies of phages from coagulase-negative staphylococci, and only two phages from Staphylococcus epidermidis have been completely sequenced and molecularly characterized (3). Other putative evidence of horizontal transfer in S. xylosus S04009 is the presence of specific sequences corresponding to antibiotic resistance as methicillin, fosfomycin, and macrolide resistance proteins (Table 2).
Particular S. xylosus surface structures.
Specific surface structures were found in the S. xylosus tester strains. The meat starter strain S. xylosus S04002 possessed a specific fragment corresponding to a 67-kDa myosin-cross-reactive streptococcal antigen-like protein. This protein family is thought to have structural features in common with the beta chain of the class II antigens. In Streptococcus pyogenes it reacts with antimyosin antibodies and plays a role in the pathogenesis of streptococcal infections (13). This gene is found in other Staphylococcus sequenced genomes, such as S. aureus, S. epidermidis, or Staphylococcus haemolyticus, but its function in those bacteria is still unknown. The S. xylosus 00-1747 strain implicated in mouse dermatitis contained a specific fragment resembling a Toll-interleukin receptor (TIR). Such a eukaryotic signal-mimicking molecule has been described in Salmonella enterica serovar Enteriditis (24). It was demonstrated that it is an important virulence factor in vivo by modulating host defense mechanisms involved in regulation of NF-
B and caspase activation. Experiments in mice infected with a Salmonella strain deficient in TIR showed reduced lethality. Predicted TIR-like domains are widespread and show an unusual distribution in bacterial species. A lateral transfer of this gene has been supposed (24).
Specific fragments shared by S. xylosus S04002, S04009, and 00-1747.
The presence of common fragments in the three S. xylosus tester strain chromosomes was detected by dot blot hybridization using genomic DNA of the three strains as a probe. The three strains shared only six specific fragments. Those six common fragments corresponded to arsenic resistance genes with the arsenate reductase and the arsenical efflux pump membrane protein, a riboflavin kinase, a hypothetical phage infection protein, and two fragments which were similar to an autolysin, a mannosyl-glycoprotein endo-β-N-acetylglucosamidase (Table 2).
Among the 122 specific fragments, 66 hybridized with S. xylosus S04009 genomic DNA. This mastitis isolate contained the highest number of specific fragments in its genome.
The S. xylosus strains S04009 and 00-1747 isolated from mastitis and dermatitis, respectively, shared the highest number of specific fragments, with 29 in common. They had specific sequences similar to a hexitol-specific PTS system, probably sorbitol specific. Indeed, other fragments related to sorbitol assimilation were isolated, such as a sorbitol dehydrogenase and the putative regulator of the sorbitol operon. Those two commensal strains also shared numerous specific DNA sequences implicated in resistance to antibiotics and other toxic components, such as tellurium, but also antibiotic synthesis-related sequences, such as a monooxygenase and a lipopeptide, micrococcin P1. These two strains also had in common the feoB gene, implicated in ferrous ion transport. This gene was identified in S. aureus and S. epidermidis genomes but had not been characterized (10). The ability to acquire iron (Fe) in a host is known to be an essential attribute of most pathogenic bacteria (2).
The meat starter strain S04002 and the strain S04009 isolated from mastitis shared 20 fragments. They have in common the iolABCDE myo-inositol genes but not the iolG gene, which was present in S04002 but absent in S04009. IolG corresponds to a myo-inositol 2-dehydrogenase, which is the first enzyme of myo-inositol catabolism, degrading the incoming myo-inositol into 2-keto-myo-inositol. The iolFHIJ and iolRS genes described in B. subtilis were not detected in the subtracted libraries (40). The gene iolR corresponds to a negative regulator of the iol operon, iolF is related to a myo-inositol transporter, and iolE codes for a 2-keto-myo-inositol dehydratase, while the functions of the other iol genes are still unclear in B. subtilis (40).
The meat starter S. xylosus S04002 had only nine fragments in common with the mouse dermatitis strain 00-1747. The biotin synthesis-related sequences obtained after genomic subtraction of these two strains seemed to be drastically different at the DNA level, and none of these fragments hybridized with the two genomic DNAs.
Conservation of tester-specific fragments within S. xylosus strains.
The extent of polymorphism among S. xylosus strains for the specific fragments was assessed. The RsaI-digested genomic DNAs of 20 different S. xylosus strains, 10 used as starters in the food industry and 10 isolated from cow or goat mastitis (Table 1), were hybridized with the 122 specific fragment sequences of the three tester strains. Four specific fragments (D39, D129, D137, and D158) from the S. xylosus 00-1747 library and four (Mb33, Mb66, M85, and M130) from the S. xylosus S04009 library were absent in the genomic DNA of the 20 S. xylosus strains. No specific fragment was present in all the S. xylosus genomes. The subtracted fragment most commonly found, which hybridized with 21 of the 23 S. xylosus strains (including the three tester strains), was fragment D294, corresponding to a hypothetical membrane protein and a CHAP domain corresponding to an amidase domain frequently found in autolysin precursor. A hierarchic ascendant classification gathered the strains into two clusters named I and II (Fig. 2).
![]() View larger version (11K): [in a new window] |
FIG. 2. Distribution of tester-specific fragments within S. xylosus strains.
|
Cluster II comprised 18 strains that shared a few fragments with cluster I, and they shared the E154 fragment, which was absent in cluster I except for one strain. This fragment corresponded to the accessory gene regulator C (agrC), which is a staphylococcal conserved regulator. It has been observed that there is a high genetic variability of the agr locus in Staphylococcus species and particularly in the agrC sequence (26). This variability was observed not only between Staphylococcus species but also between subpopulations of a single species (6). In S. aureus, agrD and agrC polymorphism permitted definition of a group broadly correlated with genotype which was correlated with pathotype (26). The existence of at least two different genotypic clusters (I and II) within the S. xylosus species was confirmed by the argC polymorphism distribution observed.
Cluster II could be divided into two subgroups named IIa and IIb. Subgroup IIa comprised the tester strain S04002 and the S01001 and S01008 strains of meat origin and six other strains isolated from mastitis (Fig. 2; Table 1). Within this heterogeneous subgroup, only two fragments were found to be in common: D294, described previously to be the most shared fragment between the 23 S. xylosus strains, and the fragment E154, already mentioned above. Subgroup IIb comprised eight S. xylosus strains used as starters (Fig. 2; Table 1) and one S. xylosus isolated from cow mastitis S04016. Twenty-one fragments were common between all these strains. They corresponded to the six myo-inositol metabolism genes (iolABCDEG), all the fermentative metabolism genes (als, acetoin reductase, and lactate permease genes), and the lipase gene (E179). The presence of iol genes in the majority of starter S. xylosus strains and their absence in the actual staphylococcal sequenced genomes indicate a specific and fortuitous acquisition of these genes by some S. xylosus strains. The presence of fermentative metabolism may be an adaptation to anaerobic growth conditions. In S. aureus, the overexpression of the genes implicated in butanediol fermentation and lactate production were recently observed under fermentative conditions (8). Adaptation to anaerobic conditions and lipase activity could procure a technological advantage for S. xylosus used in food processing.
Location of tester-specific fragments on their respective chromosomes.
The location of each specific fragment from the three DNA subtraction libraries was analyzed by Southern hybridization. Genomic DNAs of the S. xylosus S04002, S04009, and 00-1747 strains were digested with I-CeuI or SmaI, resolved by PFGE (Fig. 3), and blotted for Southern hybridization with a labeled pool of the specific fragments of each library. The specific fragments were located on the
1,100-kb I-CeuI fragment of each tester strain chromosome and on a few SmaI fragments. The specific fragments from the S. xylosus S04002 library were located on the 375-kb and 125-kb SmaI fragments, those of the S. xylosus S04009 were mainly on the 300-kb but also on the 150-kb, 100-kb, and 50-kb SmaI fragments, and the fragments of S. xylosus 00-1747 were located on the 400-kb and 100-kb SmaI fragments (Fig. 3). To confirm the probable colocalization of the SSH library-specific fragments, a second hybridization was performed using PFGE macro-restricted fragments as probes. The distribution of all the specific fragments from the three libraries is shown in Fig. 3. For each of the S. xylosus S04002, S04009, and 00-1747 subtracted libraries, more than 97.5% of inserts were located in the
1,100 kb I-CeuI fragment. Most of these fragments were located on the 300-/400-kb SmaI fragments on each S. xylosus tester strain, confirming that numerous specific fragments were colocalized on the S. xylosus S04002, S04009, and 00-1747 chromosomes. In a previous study we demonstrated by positioning the dnaA and gyrA genes, which are often found near the chromosome origin of replication of bacteria (oriC), that in S. xylosus C2a the origin of replication oriC was in the 1,138-kb I-CeuI and 365-kb SmaI fragments (5). To check if the colocalization region corresponded to the oriC of the three tester strain chromosomes, we hybridized dnaA and gyrA probes on I-CeuI or SmaI restriction fragments of S. xylosus S04002, S04009, and 00-1747 strains (data not shown). The genes dnaA and gyrA were located on the
1,100-kb I-CeuI fragment and on 300-/400-kb SmaI fragments of the three tester strains. Takeuchi et al. (35) observed in a comparative analysis of the S. haemolyticus genome with the seven S. aureus and the two S. epidermidis sequenced genomes that a region named the "oriC environ" emerged as a chromosomal region unique to each staphylococcal species. It dictates many biological features that characterize each staphylococcal species. For example, in the oriC environ of S. aureus, important virulence determinants such as protein A (spa), coagulase (coa), and the capsule locus [cap5(8)A to -P] were present. This oriC environ region was defined in the downstream region 30 kb from the oriC. It contained the staphylococcal cassette chromosomes, which might serve as a vehicle for introduction of exogenous genes into this region as well as insertion sequence elements and other recombinases serving to excise unused genes. A species-specific upstream region was also observed in some Staphylococcus strains. The authors' theory was that certain useful foreign genes were preserved by being transferred upstream of the oriC, protecting them from insertion sequence-mediated deletion occurring in the oriC environ. It was observed that the oriC environ was similar within a species but diverse among the different species, and the seven S. aureus sequenced strains had the downstream part of the oriC environ in common (35). In this study, we showed that the polymorphism of S. xylosus strains was principally centered on this oriC environ region, a species-specific region and also a strain-specific region in S. xylosus. Analysis of the differential genetic content of different S. xylosus demonstrated that this region is implicated in staphylococcal speciation but also actively participates in species evolution by the acquisition of new metabolism functions, like assimilation of new carbohydrate sources, antibiotic resistance, or factors possibly implicated in virulence. These strain-specific genes may be used as markers in epidemiological or evolutionary studies. The genomic evolution capacity is essential for ubiquitous species such as S. xylosus. Sequencing of this strain-specific "oriC environ" region in numerous Staphylococcus species and numerous strains within a species should provide a better understanding of staphylococcal genome plasticity and the genetic fitness of this genus in response to specialized niches.
![]() View larger version (58K): [in a new window] |
FIG. 3. PFGE patterns of the S. xylosus 00-1747, S04009, and S04002 tester strains. The percentage of specific fragments that hybridized with restriction fragments is detailed for each restriction pattern.
|
Emilie Dordet-Frisoni is the recipient of a fellowship from the French Ministry of "Education Nationale et Recherche." This work was financially supported by the ANR project "Genoferment" ANR-05-PNRA-020.
Published ahead of print on 21 September 2007. ![]()
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»