<title>Background</title> <p>Enterococci have emerged over the past few decades as the second to third most common cause of nosocomial infections, including urinary tract and soft tissue infections, bacteremia, and endocarditis [##REF##2404568##1##, ####REF##17609597##2##, ##REF##17889954##3####17889954##3##]. They are well equipped to thrive in environments with heavy antibiotic usage due to both their intrinsic resistance to antibiotics and their talent for swapping genetic information, which allows them to gain and share resistance determinants. Entecococcal infections are predominantly caused by <italic>E. faecalis </italic>and <italic>E. faecium</italic>. However, many, if not most, strains of these species are harmless commensals, with some enterococci being marketed in Europe to alleviate symptoms of irritable bowel syndrome and recurrent chronic sinusitis or bronchitis (Cylactin<sup>® </sup>and Fargo688<sup>® </sup>(<italic>E. faecium</italic>) and Symbioflor 1 (<italic>E. faecalis</italic>)). To differentiate the two faces of this organism, genome-wide comparisons are necessary. Although hundreds of microbial genomes have been sequenced, only two <italic>E. faecalis </italic>genomes have been reported (V583 as a clinical isolate [##REF##12663927##4##] and Symbioflor 1 as a commensal isolate [##REF##17466591##5##]), but only the V583 genome has been made publicly available. In this strain, more than one-quarter of the genome is mobile DNA, more than any other sequenced bacterial genome [##REF##12663927##4##]. The occurrence of multiple antibiotic resistance determinants in V583 [##REF##2554802##6##] makes it difficult to manipulate genetically. Moreover, the vancomycin resistance phenotype makes this strain more of a risk to handle in the laboratory. To avoid these issues, most laboratories use strain OG1 or its close derivatives. OG1 is a human isolate subsequently shown to cause dental caries in rats [##REF##807189##7##]. OG1RF is a <underline>r</underline>ifampicin and <underline>f</underline>usidic acid resistant derivative of OG1 [##REF##98769##8##,##REF##122512##9##]. By pulsed-field gel electrophoresis, Murray <italic>et al</italic>. [##REF##8349561##10##] estimated the size of the OG1RF genome as 2,825 kb and created a restriction map of the chromosome. Multilocus sequence typing (MLST) showed that OG1RF is clonally distinct from V583 (differs in six out of seven alleles of housekeeping genes) [##REF##16077117##11##] and characterization of regions flanking transposon insertions in OG1RF suggested that approximately 10% of their sequences differed [##REF##15489440##12##].</p> <p>OG1 and its derivatives have been successfully used over the past 20 years in various animal models, starting with the demonstration that it can cause caries in germ-free rats [##REF##807189##7##], and later to characterize factors important for <italic>E. faecalis </italic>virulence in a mouse model of peritonitis [##REF##9753005##13##], a rabbit model of endophthalmitis [##REF##12117982##14##], a rat model of endocarditis [##REF##16041002##15##] and in a mouse urinary tract infection model [##REF##17471437##16##]. OG1RF was also shown to be as virulent as V583 in the model host <italic>Caenorhabditis elegans </italic>[##REF##11535834##17##]. In addition to its virulence, the main reasons for the extensive use of OG1RF as a laboratory strain are that it does not carry plasmids, is readily transformable by electroporation, and is not resistant to commonly used antibiotics, other than rifampicin and fusidic acid. These resistances were serially selected in OG1 to provide strain markers [##REF##122512##9##]. The lack of resistance to common antibiotics facilitates the selection of plasmids, transposons, and allelic replacement markers introduced into the strain.</p> <p>Numerous factors important for virulence have been characterized in OG1RF. A recently described example are the Ebp pili, whose subunits are encoded by the <italic>ebp </italic>operon [##REF##17016560##18##] and whose genes are regulated by EbpR [##REF##17586623##19##]. A non-piliated mutant produces less biofilm than the parent strain and is attenuated in a rat model of endocarditis [##REF##17016560##18##] and in a murine urinary tract infection model [##REF##17471437##16##]. Also present is Ace, a member of the MSCRAMM (microbial surface component recognizing adhesive matrix molecules) family. The <italic>ace </italic>gene, like the <italic>ebp </italic>locus, is ubiquitous in <italic>E. faecalis </italic>and it occurs in at least four different forms that vary in the number of repeats of the B domain [##REF##10948146##20##]. Ace mediates conditional (that is, after growth at 46°C or in the presence of serum or collagen) adherence of <italic>E. faecalis </italic>to collagen type IV and to laminin [##REF##16926389##21##] and, in unpublished data, influences the ability of OG1RF to cause experimental endocarditis (KV Singh and BE Murray, unpublished observation). Finally, the Fsr system, a major positive and negative transcriptional regulator in OG1RF [##REF##16585749##22##], affects expression of several virulence factor genes, including <italic>gelE</italic>, which encodes gelatinase [##REF##10768947##23##], and contributes to infection in various animal models [##REF##16041002##15##,##REF##12228293##24##].</p> <p>The distinct MLST profile and the wide range of phenotypic and genotypic analyses of OG1RF, including many molecular genetic studies and experiments in various animal models, suggested that genomic analysis of this strain would prove insightful and would be useful to future studies. Thus, we analyzed the sequence of <italic>E. faecalis </italic>OG1RF. This revealed approximately 232 kb encoding 227 open reading frames (ORFs) that are unique to this important strain compared to V583. The unique regions were then characterized further.</p>
<title>Materials and methods</title> <title>Strains</title> <p><italic>E. faecalis </italic>OG1 is a strain of human origin (formerly designated 2SaR [##REF##807189##7##]) and was subsequently selected on rifampicin and fusidic acid [##REF##98769##8##,##REF##122512##9##] to generate OG1RF (deposited at the American Type Culture Collection (ATCC) under ATCC accession number 47077). V583 is a vancomycin resistant <italic>E. faecalis </italic>strain [##REF##2554802##6##], recovered from a blood culture of a patient hospitalized at the Barnes Hospital, St Louis, MO, USA in February 1987 (ATCC accession number 700802, NCBI complete genome accession number NC_004668). Bacteria were grown routinely at 37°C in BHI broth (Difco Laboratories, Detroit, MI, USA) or BHI agar unless otherwise indicated. Comparisons of OG1RF and V583 grown in broth (BHI, tryptic soy broth with glucose (TSBG), or BHI with 40% serum) did not reveal any obvious differences.</p> <title>DNA sequencing and genome assembly</title> <p>Genomic DNA was purified from cesium chloride (CsCl) gradients of whole cell lysates [##REF##8349561##10##]. DNA sequencing was performed by a combined approach using 454 Life Sciences pyrosequencing strategies [##REF##16056220##59##] and the Solexa approach [##REF##15165179##60##]. Read-pair information was used to create higher order scaffolds. Sanger sequencing was used for OG1RF whole gun sequencing and finishing. The coverage was 28× by the 454, 104× by Solexa, and 4.5× by Sanger sequencing technique. The assembly was done using Atlas [##REF##15060016##61##].</p> <title>Gene identification and annotation</title> <p>Gene prediction and manual annotation were performed as previously described [##REF##17895969##62##]. Glimmer [##REF##10556321##63##] and GeneMark [##REF##9461475##64##] were used independently to predict ORFs. Visualization of gene predictions was performed using the Genboree system [##UREF##0##65##] and the CONAN database [##REF##15317790##66##]. OG1RF-unique ORFs were analyzed with BLASTN and BLASTX. Protein sequences were analyzed by BLASTP versus the nr database at NCBI [##REF##15608222##67##]. When appropriate, other predictive tools were used as described previously [##REF##17895969##62##]. This whole genome shotgun project has been deposited at DDBJ/EMBL/GenBank under the project accession ABPI00000000. The version described in this paper is the first version, ABPI01000000. This project includes also the annotation of the ORFs unique to OG1RF. The OG1RF-unique ORFs are listed in Additional data file 1.</p> <title>Transposon mutations in OG1RF-unique sequences</title> <p>Following the creation of an <italic>E. faecalis </italic>Tn<italic>917 </italic>library [##REF##15489440##12##], 6,237 sequences representing the flanking regions of the transposon insertion sites were obtained and compared to the V583 genome by BLASTN. A total of 196 sequences were unique to OG1RF. Thirty-seven of the unique genes contained a transposon insertion. The locations of the transposon insertions are listed in Additional data file 1.</p> <title>Carbohydrate fermentation tests</title> <p>Forty-eight <italic>E. faecalis </italic>isolates, including OG1RF and V583, having different MLST profiles, pulsed field gel electrophoresis types or from various geographical origins, were streaked onto BHI agar and incubated overnight at 37°C. Five to ten colonies of each strain were resuspended in 100 μl of 0.9% saline in a microtiter plate and tested for fermentation using BBL™ Phenol Red Broth Base (Diagnostic Systems, Sparks, MD, USA) supplemented with agar and either 10 mM glucose (positive control), 10 mM dulcitol (negative control), or 10 mM myo-inositol (Sigma, St Louis, MO, USA). Plates were read after incubation at 37°C for 24 h; a yellow halo around the colony was considered positive for fermentation. <italic>iolB </italic>and <italic>iolG2 </italic>transposon mutants [##REF##15489440##12##] were also tested.</p> <title>PCR</title> <p>PCR was performed using genomic DNA purified using Bactozol™ (Molecular Research Center, Inc., Cincinnati, OH, USA), as recommended by the manufacturer. Specific PCR primer pairs (Additional data file 4) were used to assess the presence of the OG1RF-unique sequences and for confirmation of flanking DNA regions in common with V583.</p> <title>Competence assays</title> <p>To test strains for competence, overnight cultures, grown at 37°C in Todd-Hewitt broth, were diluted in Todd-Hewitt broth to an OD<sub>600 nm </sub>of 0.05 and then further diluted 10,000-fold in Todd-Hewitt broth to a final volume of 100 ml. After 30 minutes at 37°C, with shaking at 150 rpm, and every hour for 10 h, 0.5 ml samples were removed and 2.5 μg of plasmid DNA were added. The plasmids tested were pAM401 [##REF##3005240##48##] and pMSP3535VA [##REF##10964628##49##]. The samples were incubated for 2 h before plating on BHI or BHI plus antibiotic (chloramphenical 10 μg/ml for pAM401 or kanamycin 2 mg/ml for pMSP3535VA). Following overnight incubation at 37°C, the total numbers of CFU/ml recovered on selective agar for the plasmid were compared to the total number of CFU/ml (plated on BHI agar) for each time point.</p> <title>Biofilm assay and statistical analysis</title> <p>The biofilm assay was performed as described by Mohamed <italic>et al</italic>. [##REF##15155680##68##]. Each assay was performed using 16 wells on three occasions. The median was calculated using the 48 OD<sub>570 nm </sub>readings on data pooled from all experiments and statistical analysis was performed using a non-parametric <italic>t</italic>-test.</p> <title>Mouse peritonitis model</title> <p><italic>E. faecalis </italic>strains OG1RF and V583 [##REF##2554802##6##] were tested using a previously published method [##REF##9753005##13##]. Mice were injected intraperitoneally with appropriate dilutions of premixed bacteria/sterile rat fecal extract and were observed for five days. Two-fold dilutions of test bacteria in the range 10<sup>7</sup>-10<sup>9 </sup>CFU were used as the inocula for LD<sub>50 </sub>determination using 6-9 mice per inoculum group. Inocula CFU geometric mean values were obtained and used for LD<sub>50 </sub>calculation by the method of Reed and Muench [##UREF##1##69##].</p> <title>UTI model for competition assay and ID50 determination</title> <p><italic>E. faecalis </italic>strains OG1RF and V583 were tested in the UTI model as previously described [##REF##17471437##16##]. For the mixed inoculum experiments, an approximately 1:1 ratio of <italic>E. faecalis </italic>OG1RF:V583 at approximately 10<sup>3 </sup>CFU each was used. Two independent experiments, using 10 and 11 mice, respectively, were performed and the results were combined. The log<sub>10</sub>(CFU) of OG1RF and V583 per gram of tissue of each animal (kidney or bladder) from mixed infection were analyzed for significance by the paired <italic>t</italic>-test. For mono-infection, approximately 10<sup>3 </sup>CFU organisms grown in BHI + 40% horse serum were used for each strain independently and CFU obtained from kidney pairs (nine mice per strain) were analyzed for significance by the unpaired <italic>t</italic>-test. The minimum detectable limits of recovered bacteria were 10<sup>1 </sup>and 10<sup>2 </sup>CFU/gm of kidney pairs and urinary bladder homogenates, respectively.</p>
<title>Results and discussion</title> <title>General genome features</title> <p>The complete circular chromosome of OG1RF was found to be 2,739,625 bp with an average G+C content of 37.8%. The complete OG1RF sequence was obtained using three independent techniques (Solexa, the 454, and Sanger sequencing technique) with a higher than classic coverage (more than 100 times), diminishing the likelihood of sequencing-related frameshifts, base errors and/or misassembly. A comparison of our assembly of the closed OG1RF genome with the restriction map of OG1RF published by Murray <italic>et al</italic>. [##REF##8349561##10##] showed only minor variations (primarily an overestimation of 30 kb for the <italic>Sfi </italic>I fragment E, 540 kb versus 509 kb predicted from the sequence; Figure ##FIG##0##1##).</p> <p>We found 232 kb of OG1RF unique sequences distributed in 48 regions ranging from 101 bp to approximately 49 kb in length (Figure ##FIG##0##1##; Additional data file 1). Using the published DNA sequence of V583 as reference (NC_004668), OG1RF shares 2,474 ORFs as well as the 12 rRNA genes and 58 of 68 tRNA genes (Table ##TAB##0##1##). The 10 missing tRNA are localized in a region in V583 that has been replaced in OG1RF by a 49 kb region (see below). Surprisingly, the genomes align syntenically, as shown in Figure ##FIG##1##2##, despite the fact that 25% of the V583 genome is composed of mobile elements. Similarly, the presence of OG1RF-unique sequences has not affected the overall chromosomal arrangement. Some of the major insertions/deletions in the two genomes are shown in Figure ##FIG##1##2##, such as the absence of the pathogenicity island (PAI) in OG1RF and the presence of an approximately 49 kb fragment unique to OG1RF. However, most of the differences are small and cannot be visualized in this figure. Overall, we found 64 areas of divergence between the genomes that can be divided into 3 classes: an additional sequence present in OG1RF when compared with V583; a sequence replacement where a sequence in OG1RF differs from the sequence in V583; and the absence of a sequence from OG1RF when compared with V583.</p> <title>CRISPR loci</title> <p>The CRISPR (comprised of regularly interspaced short palindromic repeats) loci encoded by some bacterial strains is a recently described system that protects cells from infection with bacteriophage [##REF##17379808##25##, ####REF##18065539##26##, ##REF##18065545##27####18065545##27##]. The specificity of the phage resistance conferred by the CRISPR elements and CRISPR-associated genes (<italic>cas </italic>genes) is determined by spacer-phage sequence similarity. OG1RF carries two CRISPR elements: CRISPR1 (between the OG1RF homologue of EF0672 and EF0673) and CRISPR2 (between the OG1RF homologue of EF2062 and EF2063); CRISPR1 is linked to <italic>cas</italic>-like genes while CRISPR2 is not (Figure ##FIG##2##3##). Both OG1RF CRISPR elements are composed of 7 repeats of a 37 bp palindromic sequence with a 29 bp spacer. None of the 29 bp spacers (14 total) have homology to any sequences in GenBank. The CRISPR1-associated proteins belong to the Nmeni subtype [##REF##16292354##28##]. Species bearing this CRISPR/<italic>cas </italic>subtype have so far been found exclusively in bacteria that are vertebrate pathogens or commensals. The Nmeni subtype is characterized by the presence of four specific <italic>cas </italic>genes and a single copy of the repeat that is upstream of the first gene in the locus. The four <italic>cas </italic>genes encode Cas_csn1 (possible endonuclease), Cas1 (novel nuclease), Cas2 (conserved hypothetical protein), and Cas_csn2 (conserved hypothetical protein). The repeat upstream of <italic>cas_csn1 </italic>appears to have degenerated since it shares only 23 bp with the 37 bp repeat cluster downstream of the last gene. A unique feature of the OG1RF CRISPR1 locus is the presence of a gene downstream of the element, which encodes a hypothetical 119 amino acid transmembrane protein.</p> <p>The presence of the CRISPR loci among <italic>E. faecalis </italic>strains may be a powerful tool to avoid the load of prophage replication. To determine the distribution of the CRISPR1 locus in <italic>E. faecalis </italic>strains, 16 isolates of various MLST types were tested for the presence (PCR with primers specific for <italic>csn1 </italic>and <italic>cas1</italic>) or absence (PCR with primers overlapping the junction between EF0672 and EF0673) of the CRISPR1 locus (Table ##TAB##1##2##). Seven strains were <italic>cas </italic>positive, but negative for the junction and the remaining nine were positive only for the junction. This indicates that the location of the CRISPR1 locus appears to be conserved (between EF0672 and EF0673 when compared with the V583 genome). Interestingly, the two vancomycin resistant strains tested were both <italic>cas </italic>negative. It is appealing to postulate that the presence of the CRISPR locus in OG1RF may be the reason for the absence of prophage in this strain.</p> <title>A 14.8 kb region inserted in the 23.9 kb region containing <italic>fsrA </italic>and <italic>fsrB</italic></title> <p>Nakayama <italic>et al</italic>. [##REF##12039782##29##] described a conserved 23.9 kb chromosomal deletion when comparing <italic>fsrA-</italic>lacking/<italic>fsrC</italic><sup>+</sup>/<italic>gelE</italic><sup>+ </sup>strains (by PCR) from various origins with V583; the deleted sequences start in the middle of EF1841, include the <italic>fsrAB </italic>genes and end in the middle of the <italic>fsrC </italic>gene (EF1820). Loss of the <italic>fsr </italic>regulatory components results in a gelatinase-negative phenotype under routine test conditions despite the fact that these strains still carry the <italic>gelE </italic>gene [##REF##10768947##23##,##REF##12039782##29##]. The absence/presence of the 23.9 kb region, from EF1820/<italic>fsrC </italic>to EF1841, did not appear to correlate with the clinical origin of the isolates [##REF##15131223##30##]. In a more recent analysis of relationships between various <italic>E. faecalis </italic>strains, the 23.9 kb region was not detected in 86% of the strains of the clonal complex (CC)2, 58% of the CC9 strains, nor in any of the CC8 strains [##REF##17611618##31##]. The Symbioflor 1 strain, used as a probiotic, is one representative of the 7.4% of <italic>E. faecalis </italic>isolates that are missing the <italic>gelE </italic>gene in addition to the 23.9 kb region [##REF##17466591##5##,##REF##15131223##30##]. Our analysis of this area in OG1RF revealed the presence of an additional 14.8 kb fragment inserted between the corresponding EF1826 and EF1827 of OG1RF (confirmed by PCR; results not shown). In OG1RF, this 14.8 kb region contains two loci, a WxL locus (described below) and a seven-gene locus that may encode a possible ABC transporter with similarity to one annotated in <italic>Pediococcus pentosaceus</italic>.</p> <title>Components of the cell surface</title> <p>It has been shown in <italic>E. faecalis </italic>that at least one cell surface protein (Ace) is subject to domain variation [##REF##10948146##20##] and it has been postulated that domain variation may help bacteria escape the immune system. We found more polymorphisms in two families of <italic>E. faecalis </italic>proteins present on the cell surface: the MSCRAMMs and the WxL domain surface proteins. The MSCRAMMs are composed of two large regions, namely, the non-repeat A region (which is usually the ligand binding region for extracellular matrix molecules such as collagen or fibrinogen) and the B region (which typically contains repeated sub-domains). The B region of Ace contains five repeats in OG1RF, while it contains only four in V583 [##REF##10948146##20##]. We found two other MSCRAMM proteins that show polymorphisms in the number of their B repeats. OG1RF_0186 (corresponding to EF2505 of V583) has four repeats compared to seven in V583, and OG1RF_0165 (corresponding to EF2224 of V583) has eight repeats compared to five in V583. It has been proposed that the repeats are used as a stalk that projects the A region across the peptidoglycan and away from the cell surface [##REF##9799504##32##]. A hypothesis that the number of repeats may be proportional to the depth of the peptidoglycan has been proposed [##REF##9799504##32##]. However, OG1RF_0186 carries fewer repeats than EF2505 while Ace and OG1RF_0165 carry more repeats than their counterparts in V583, suggesting that our observation does not fit this hypothesis or that the peptidoglycan depth is not uniform. Apart from these three MSCRAMMs with B-repeat polymorphisms, we identified two unique MSCRAMM proteins in OG1RF: a homologue of EF0089 (OG1RF_0063, which shares 48% similarity) and a homologue of EF1896 (OG1RF_0039, which shares 75% similarity); both are located in the approximately 49 kb region unique to OG1RF described below (Figure ##FIG##0##1##; Additional data file 1).</p> <p>Another family of <italic>E. faecalis </italic>surface proteins includes the newly described WxL domain surface proteins. Siezen <italic>et al</italic>. [##REF##16723015##33##] reported a novel gene cluster encoding exclusively cell-surface proteins that is conserved in a subgroup of Gram-positive bacteria. Each gene cluster has at least one member of three gene families: a gene encoding a small LPxTG protein (approximately 120 amino acids); a gene encoding a member of the DUF916 transmembrane protein family; and a gene encoding a WxL domain surface protein. In addition, members of these gene families were found as singletons or associated with genes encoding other proteins (Additional data file 2). Recently, it was shown that the WxL domain attaches to the peptidoglycan on the cell surface [##REF##16963569##34##] and one member of this WxL domain family, the homologue of EF2686 in OG1RF (a probable internalin protein), was shown to be important for virulence in a mouse peritonitis model and is required for dissemination to the spleen and liver [##REF##17620355##35##]. OG1RF shares five complete WxL loci with V583 (EF0750-7, EF2682-6, EF2970-68, EF3181-8, and EF3248-53). OG1RF does not contain homologues of EF2248-54 (carrying instead the <italic>iol </italic>operon), though it has a novel WxL locus within the 14.8 kb unique region upstream of the <italic>fsr </italic>locus (Additional data file 2). In addition to the variation in the number of WxL loci, we also observed polymorphisms in six of the WxL domain surface proteins. For example, OG1RF_0213 shares 88% similarity with EF3188, while OG1RF_0224, OG1RF_0225, and OG1RF_0227 share 64-68% similarity with their V583 counterparts, EF3248, EF3250, and EF3252, respectively. Also, in place of EF3153, EF3154, and EF3155 (which share 70% similarity among themselves), were found non-distantly related homologues, OG1RF_0209 and OG1RF_0210, which share 60-80% similarity with EF3153, EF3154, and EF3155. It is interesting to note that while several of these WxL loci, including the EF0750 and EF3184 loci, were present by hybridization in all the strains (clinical or food isolates) tested by Lepage <italic>et al</italic>. [##REF##16980489##36##], other loci, including the EF3153 and EF3248 loci, were not detected in the majority of these strains. In addition, it appears that the EF3248 locus diverges in the Symbioflor 1 strain. When compared to V583, the sequence identity in this area between the two strains appears to be as low as 75% (depicted in Figure ##FIG##1##2## from reference [##REF##17466591##5##]).</p> <p>However, because the Symbioflor 1 genome sequence is not currently available, it was not possible to compare their respective sequences in more detail. Since these proteins are located at the surface of the cell, the low level of homology shared between them may be the result of antigenic variation. More analyses are required for a better understanding of the number, frequency and function of these WxL domain proteins and their possible relationship with the diversity of <italic>E. faecalis</italic>.</p> <p>Finally, as previously found using PCR, the <italic>cpsCDEFGHIJK </italic>operon capsule polysaccharide genes [##REF##15184433##37##] were confirmed here as missing, although OG1RF carries the <italic>cpsA </italic>and <italic>cpsB </italic>genes, which were proposed to be essential for <italic>E. faecalis </italic>since all strains tested by Hufnagel <italic>et al</italic>. [##REF##15184433##37##] carry these two genes. In OG1RF, the region that would encode the <italic>cps </italic>operon is only 59 bp in length and has no homology with V583. Thus, while V583 and OG1RF share much similarity between their surface components, there are unique differences that could potentially be important in affecting the behavior of the strains and might be useful for strain typing.</p> <title>Two-component regulatory systems</title> <p>OG1RF lacks four two-component systems found in V583. These are histidine kinase-response regulator (HK-RR)08, HK-RR12 located in the PAI, HK-RR16 and the <italic>vanB </italic>regulatory system HK-RR11 [##REF##12374813##38##]. However, an OG1RF-unique two-component system with high homology with the <italic>van</italic><sub>G </sub>locus was found at the location corresponding to the region between EF2860 and EF2861 in V583 (Table ##TAB##2##3##). OG1RF_0193 shares 82% similarity with VanR<sub>G </sub>and 81% similarity with VanR<sub>G2</sub>. Similarly, OG1RF_0192 shares 68% similarity with VanS<sub>G </sub>and VanS<sub>G2</sub>. A gene (OG1RF_0191) encoding an M15 family muramoyl pentapeptide carboxypeptidase is located downstream of these two-component regulatory genes (Figure ##FIG##3##4a##). The predicted carboxypeptidase (OG1RF_0191) shares 69% similarity over 179 amino acids with EF2297, a membrane-associated D, D-carboxypeptidase encoded by the <italic>vanB </italic>operon in V583. However, OG1RF_0191 lacks an identifiable transmembrane domain that is important to the VanY function and it is likely, therefore, that this protein may be a soluble D, D-carboxypeptidase/transpeptidase as seen in <italic>Streptomyces </italic>[##REF##3997832##39##] and <italic>Actinomadura </italic>[##REF##15987687##40##], and thus may not be involved in peptidoglycan metabolism. Consequently, it seems unlikely that this operon is a remnant of a vancomycin resistance operon in OG1RF, but rather part of a still unknown regulatory pathway.</p> <title>The <italic>iol </italic>operon</title> <p>OG1RF carries an <italic>iol </italic>operon while V583 does not. This operon encodes the factors necessary for the degradation of myo-inositol into glyceraldehyde-3P. Many soil and plant micro-organisms, including <italic>Bacillus subtilis </italic>[##REF##9226270##41##] (first <italic>iol </italic>operon identified), <italic>Klebsiella </italic>spp. [##REF##977784##42##], and cryptococci [##REF##394818##43##], have been reported to use myo-inositol as a sole carbon source. Myo-inositol, one of the nine isomers of the inositol group, belongs to the cyclitol group and is abundant in nature, particularly in the soil. The OG1RF <italic>iol </italic>operon appears to be closely related to ones described in <italic>Clostridium perfringens </italic>[##REF##15183876##44##] and <italic>Lactobacillus casei </italic>[##REF##17449687##45##]. In <italic>L. casei</italic>, the myo-inositol operon consists of ten genes with an upstream divergent regulator gene, <italic>iolR</italic>. In OG1RF, the operon appears to include ten genes, beginning with a probable transcriptional regulator (helix-turn-helix domain protein). Also, the OG1RF operon carries two copies of an <italic>iolG</italic>-like gene, which encodes inositol 2-dehydrogenase, the first enzyme of the myo-inositol degradation pathway (Figure ##FIG##4##5##). However, the order of the genes is not the same between <italic>E. faecalis </italic>and <italic>L. casei</italic>. In addition, <italic>iolH</italic>,<italic>iolJ </italic>and <italic>iolK</italic>, present in <italic>L. casei</italic>, are not present in OG1RF, nor are <italic>iolH </italic>and <italic>iolK </italic>present in the <italic>C. perfringens iol </italic>operon.</p> <p>Yebra <italic>et al</italic>. reported that <italic>L. casei </italic>was the sole member of the Lactobacillales to carry a functional <italic>iol </italic>operon [##REF##17449687##45##]. To survey <italic>E. faecalis</italic>, also a member of this order, for the presence of the <italic>iol </italic>operon, 48 isolates with different MLST and/or from various origins (including OG1RF and V583) were tested for myo-inositol fermentation; 23 of 48 isolates were positive. In addition, PCR verified the presence of <italic>iolE </italic>and <italic>iolR </italic>in these strains and in one negative for myo-inositol fermentation, indicating that the <italic>iol </italic>operon is not unique to OG1RF. To verify that the <italic>iol </italic>genes are responsible for the fermentation of myo-inositol in OG1RF, transposon insertion mutants [##REF##122512##9##] in the <italic>iolB </italic>and <italic>iolG2 </italic>genes of OG1RF were tested. Both mutants failed to ferment myo-inositol (data not shown), demonstrating that these genes are essential for myo-inositol fermentation. To investigate whether the <italic>iol </italic>operon was 'inserted into' or 'removed from' a putative ancestral strain, the sequences surrounding the <italic>iol </italic>genes were examined. In OG1RF, the <italic>iol </italic>operon is located between the equivalent of EF2239 and EF2352 when compared with V583. In V583, this region encodes probable prophage proteins and carries the <italic>vanB </italic>transposon, which confers vancomycin resistance. Since we did not identify any remnants of the <italic>iol </italic>operon in V583, it would appear that at least two independent events at the same location differentiate OG1RF and V583, suggesting that it is a hot region for rearrangement. This region between EF2239 and EF2352 (111 Kb) is also missing in the Symbioflor 1 strain (referred to as gap 2) [##REF##17466591##5##]. The possible junction and presence of unique sequence in this region, if investigated, was not mentioned in the publication. Nonetheless, preliminary analysis of other strains' genotypes in this area seemed to confirm the hypothesis of a hot region for rearrangement (data not shown).</p> <title>A homologue of Tn<italic>916 </italic>in OG1RF</title> <p>An analysis of the G+C content of OG1RF unique regions revealed several loci with a lower G+C content than the 37.8% average content of OG1RF. One of these is an approximately 49 kb fragment with a G+C content of 32.1% located between an rRNA operon and the homologue in OG1RF of EF1053, replacing 10 tRNA genes present in V583 (Figure ##FIG##0##1##). This fragment appears to be a patchwork composed of hypothetical genes, homologues of Tn<italic>916</italic>-associated genes and homologues of genes from other Gram-positive organisms, including <italic>Listeria</italic>, <italic>E. faecium</italic>, staphylococci, or lactococci (Additional data file 1). It is interesting to note that this region contains: a putative adhesin protein gene (OG1RF_0039) at one end of the fragment; homologues of 14 Tn<italic>916</italic>-associated genes (Tn<italic>916</italic>_2 to Tn<italic>916</italic>_12, Tn<italic>916</italic>_18 and Tn<italic>916</italic>_19, with an average of 70% similarity); and a gene encoding a putative integrase (OG1RF_0088) at the other end - these three features are also present in Tn<italic>5386 </italic>in <italic>E. faecium </italic>D344R [##REF##17408741##46##]. However, the approximately 49 kb fragment lacks an excisase gene and the probable lantibiotic ABC transporter genes present in Tn<italic>5386</italic>.</p> <title>An uninterrupted competence operon in OG1RF</title> <p>OG1RF contains what appears to be an intact competence operon while that of V583 appears to be non-functional. This operon in OG1RF is similar to a nine-gene operon described in <italic>Streptococcus mutans </italic>[##REF##15632435##47##], as shown in Figure ##FIG##5##6##. For example, the homologue in OG1RF of EF2046 shares 61% similarity with ComYA and the OG1RF homologue of EF2045 is 55% similar to ComYB. In <italic>S. mutans</italic>, only the first seven genes of the operon are essential for competence [##REF##15632435##47##]. In V583, the fourth gene of this operon (corresponding to OG1RF_0148) is interrupted by phage 4 (EF1896-EF2043); in addition, EF1984 contains a premature stop codon not found in the corresponding gene in OG1RF (OG1RF_0228).</p> <p>Natural competence has not been reported for <italic>E. faecalis</italic>. To assess the functionality of this operon in OG1RF, we evaluated the competence of cells in different phases of growth (early log growth to stationary phase) using pAM401 [##REF##3005240##48##] and pMSP3535VA [##REF##10964628##49##]. We were not able to show natural competence under the conditions tested. We have also noted that V583 is less transformable by electroporation than OG1RF. To investigate the possibility that directly or indirectly the <italic>com </italic>operon might be responsible for this phenotype, we also evaluated transformability by electroporation. When compared with OG1RF, transposon mutants [##REF##15489440##12##] in the OG1RF equivalent of EF2045 (encoding the <italic>comGB </italic>homologue) and in the OG1RF equivalent of EF1986 (encoding the <italic>comGF </italic>homologue) showed similar levels of transformability by electroporation (data not shown), implying that the difference in electroporation efficiency observed between OG1RF and V583 is not related to this locus.</p> <p>In <italic>Streptococcus pneumoniae </italic>[##REF##9701804##50##], the competence operon is tightly regulated by a quorum sensing two-component system (ComDE) and a competence-stimulating peptide (CSP; encoded by <italic>comC</italic>). We did not find any homologues of CSP in OG1RF. Two homologous ComDE sensor histidine kinase/response regulators were found in OG1RF, one of which is FsrC/FsrA. Based on our previous microarray data, the Fsr system does not regulate the <italic>comY </italic>operon, at least under our previously used conditions (mid-log phase growth to early stationary phase in brain heart infusion (BHI)) [##REF##16585749##22##]. The other ComDE homology is that with a two-component system unique to OG1RF (OG1RF_0199 and OG1RF_0198, respectively) that lies on a 4,706 bp unique fragment that maps between EF3114 and EF3115 in V583. This fragment also carries two genes (OG1RF_0200 and OG1RF_0201) encoding homologues of the YhaQ and YhaP sodium efflux ATP-binding cassette efflux/transporter proteins (Figure ##FIG##3##4b##). Although they are potential elements of a secretion apparatus, these two proteins do not share any homology at the protein level with the competence secretion apparatus ComAB of <italic>S. pneumoniae </italic>[##REF##11115120##51##] nor CslAB from <italic>S. mutans </italic>[##REF##11154426##52##]. Searching for a possible CSP in the vicinity of these genes, we identified a small ORF encoding 50 amino acids between <italic>yhaP </italic>and OG1RF199 and another encoding 20 amino acids downstream of OG1RF198. More analysis will be required to determine if there are conditions in which the OG1RF <italic>com </italic>operon is expressed and to determine whether or not this two-component system is involved in competence.</p> <title>Limited presence of mobile elements</title> <p>By probing a microarray of the V583 genome and plasmids with OG1RF genomic DNA, we previously estimated that only 75% of V583 ORFs were also present in OG1RF [##REF##16585749##22##]. Later, Aakra <italic>et al</italic>. [##REF##17220255##53##] compared nine strains, including OG1RF to V583, using comparative genomic hybridization. In these results, OG1RF appears to carry a few genes included in the PAI, and a few prophage genes. Using the complete genome sequence, we have now found that OG1RF lacks 639 genes and the three pTET plasmids described in V583. All but 45 of the missing genes are associated with putative mobile elements, such as the entire PAI, the recently described phages 1, 3, 4, 5, 6, and 7 [##REF##17611618##31##], and the approximately 111 kb area between genes EF2240 and EF2351 (including the <italic>vanB </italic>transposon) present in V583. The absence of these elements appears also to be a characteristic of the commensal strain Symbioflor 1, although because the genome was not completely finished, the possibility remains that some of these regions were not sequenced. In conclusion, other than the approximately 49 kb fragment containing a Tn<italic>916 </italic>homologue, it appears that OG1RF has only one additional possible mobile element derivative, namely the phage 2 proposed to be part of the core genome [##REF##17611618##31##].</p> <title>Fusidic acid and rifampicin resistance</title> <p>OG1RF was sequentially selected from OG1 for resistance to fusidic acid and rifampicin [##REF##122512##9##]. The mutation leading to rifampicin resistance was identified in the <italic>rpoB </italic>gene by Ozawa <italic>et al</italic>. [##REF##15907539##54##] and is caused by the A1467G mutation, which results in substitution of arginine for histidine at amino acid 489. The mutation also affected the clumping phenotype of <italic>traA </italic>mutants and this effect appears to be specific for the pAM373 system [##REF##15907539##54##]. All of the other 22 differences in <italic>rpoB </italic>between OG1RF and V583 are synonymous. Fusidic acid resistance is associated with mutation(s) in the <italic>fusA </italic>gene, which encodes elongation factor G. We compared <italic>fusA </italic>from OG1RF with that in V583 and identified two differences (C1368A and T1992C). The mutation T1992C is synonymous, while the mutation C1368A leads to the presence of glutamine (histidine in V583) at position 404 in OG1RF. Mutations in this region have been associated with fusidic acid resistance in <italic>Staphylococcus aureus </italic>[##REF##11036042##55##,##REF##11309125##56##], and thus the C1368A mutation is likely the cause of the fusidic acid resistance phenotype in OG1RF.</p> <title>Virulence and biofilm comparisons of OG1RF with V583</title> <p>When compared in the mouse peritonitis model, the LD<sub>50 </sub>values of V583 in different determinations were lower (4.8 × 10<sup>7 </sup>to 1.1 × 10<sup>8 </sup>colony forming units (CFU)/ml) than the LD<sub>50 </sub>values of OG1RF (1.2 × 10<sup>8 </sup>to 4.8 × 10<sup>8 </sup>CFU/ml). However, at comparable inoculum, OG1RF (4 × 10<sup>8 </sup>CFU/ml) showed more rapid mortality versus V583 (5 × 10<sup>8 </sup>CFU/ml) in the first 48 hours (<italic>P </italic>= 0.0034; Additional data file 3). In a urinary tract infection model administering mixed equal inocula of V583 and OG1RF, OG1RF significantly outnumbered V583 in kidney with geometric means of 1.3 × 10<sup>4 </sup>CFU/gm for OG1RF versus 1.9 × 10<sup>2 </sup>CFU/gm for V583 (<italic>P </italic>= 0.0005); in urinary bladder homogenates, the geometric mean CFU/gm was 1.7 × 10<sup>3 </sup>for OG1RF versus 6.6 × 10<sup>1 </sup>for V583 (<italic>P </italic>= 0.003; Figure ##FIG##6##7a##). Similarly, in mono-infection, the geometric mean CFU/gm of OG1RF in kidneys was 9.4 × 10<sup>3 </sup>versus 4 × 10<sup>1 </sup>for V583 (<italic>P </italic>= 0.0035; Figure ##FIG##6##7b##). We also found that OG1RF produced 20% more biofilm (<italic>P </italic>&lt; 0.01) than V583 at 24 hours (results not shown). These results, together with the previous results in <italic>C. elegans </italic>[##REF##11535834##17##], demonstrate that OG1RF, although lacking what was thought to be important for virulence (PAI, plasmids, prophages), is as pathogenic as V583 in at least three assays.</p>
<title>Results and discussion</title> <title>General genome features</title> <p>The complete circular chromosome of OG1RF was found to be 2,739,625 bp with an average G+C content of 37.8%. The complete OG1RF sequence was obtained using three independent techniques (Solexa, the 454, and Sanger sequencing technique) with a higher than classic coverage (more than 100 times), diminishing the likelihood of sequencing-related frameshifts, base errors and/or misassembly. A comparison of our assembly of the closed OG1RF genome with the restriction map of OG1RF published by Murray <italic>et al</italic>. [##REF##8349561##10##] showed only minor variations (primarily an overestimation of 30 kb for the <italic>Sfi </italic>I fragment E, 540 kb versus 509 kb predicted from the sequence; Figure ##FIG##0##1##).</p> <p>We found 232 kb of OG1RF unique sequences distributed in 48 regions ranging from 101 bp to approximately 49 kb in length (Figure ##FIG##0##1##; Additional data file 1). Using the published DNA sequence of V583 as reference (NC_004668), OG1RF shares 2,474 ORFs as well as the 12 rRNA genes and 58 of 68 tRNA genes (Table ##TAB##0##1##). The 10 missing tRNA are localized in a region in V583 that has been replaced in OG1RF by a 49 kb region (see below). Surprisingly, the genomes align syntenically, as shown in Figure ##FIG##1##2##, despite the fact that 25% of the V583 genome is composed of mobile elements. Similarly, the presence of OG1RF-unique sequences has not affected the overall chromosomal arrangement. Some of the major insertions/deletions in the two genomes are shown in Figure ##FIG##1##2##, such as the absence of the pathogenicity island (PAI) in OG1RF and the presence of an approximately 49 kb fragment unique to OG1RF. However, most of the differences are small and cannot be visualized in this figure. Overall, we found 64 areas of divergence between the genomes that can be divided into 3 classes: an additional sequence present in OG1RF when compared with V583; a sequence replacement where a sequence in OG1RF differs from the sequence in V583; and the absence of a sequence from OG1RF when compared with V583.</p> <title>CRISPR loci</title> <p>The CRISPR (comprised of regularly interspaced short palindromic repeats) loci encoded by some bacterial strains is a recently described system that protects cells from infection with bacteriophage [##REF##17379808##25##, ####REF##18065539##26##, ##REF##18065545##27####18065545##27##]. The specificity of the phage resistance conferred by the CRISPR elements and CRISPR-associated genes (<italic>cas </italic>genes) is determined by spacer-phage sequence similarity. OG1RF carries two CRISPR elements: CRISPR1 (between the OG1RF homologue of EF0672 and EF0673) and CRISPR2 (between the OG1RF homologue of EF2062 and EF2063); CRISPR1 is linked to <italic>cas</italic>-like genes while CRISPR2 is not (Figure ##FIG##2##3##). Both OG1RF CRISPR elements are composed of 7 repeats of a 37 bp palindromic sequence with a 29 bp spacer. None of the 29 bp spacers (14 total) have homology to any sequences in GenBank. The CRISPR1-associated proteins belong to the Nmeni subtype [##REF##16292354##28##]. Species bearing this CRISPR/<italic>cas </italic>subtype have so far been found exclusively in bacteria that are vertebrate pathogens or commensals. The Nmeni subtype is characterized by the presence of four specific <italic>cas </italic>genes and a single copy of the repeat that is upstream of the first gene in the locus. The four <italic>cas </italic>genes encode Cas_csn1 (possible endonuclease), Cas1 (novel nuclease), Cas2 (conserved hypothetical protein), and Cas_csn2 (conserved hypothetical protein). The repeat upstream of <italic>cas_csn1 </italic>appears to have degenerated since it shares only 23 bp with the 37 bp repeat cluster downstream of the last gene. A unique feature of the OG1RF CRISPR1 locus is the presence of a gene downstream of the element, which encodes a hypothetical 119 amino acid transmembrane protein.</p> <p>The presence of the CRISPR loci among <italic>E. faecalis </italic>strains may be a powerful tool to avoid the load of prophage replication. To determine the distribution of the CRISPR1 locus in <italic>E. faecalis </italic>strains, 16 isolates of various MLST types were tested for the presence (PCR with primers specific for <italic>csn1 </italic>and <italic>cas1</italic>) or absence (PCR with primers overlapping the junction between EF0672 and EF0673) of the CRISPR1 locus (Table ##TAB##1##2##). Seven strains were <italic>cas </italic>positive, but negative for the junction and the remaining nine were positive only for the junction. This indicates that the location of the CRISPR1 locus appears to be conserved (between EF0672 and EF0673 when compared with the V583 genome). Interestingly, the two vancomycin resistant strains tested were both <italic>cas </italic>negative. It is appealing to postulate that the presence of the CRISPR locus in OG1RF may be the reason for the absence of prophage in this strain.</p> <title>A 14.8 kb region inserted in the 23.9 kb region containing <italic>fsrA </italic>and <italic>fsrB</italic></title> <p>Nakayama <italic>et al</italic>. [##REF##12039782##29##] described a conserved 23.9 kb chromosomal deletion when comparing <italic>fsrA-</italic>lacking/<italic>fsrC</italic><sup>+</sup>/<italic>gelE</italic><sup>+ </sup>strains (by PCR) from various origins with V583; the deleted sequences start in the middle of EF1841, include the <italic>fsrAB </italic>genes and end in the middle of the <italic>fsrC </italic>gene (EF1820). Loss of the <italic>fsr </italic>regulatory components results in a gelatinase-negative phenotype under routine test conditions despite the fact that these strains still carry the <italic>gelE </italic>gene [##REF##10768947##23##,##REF##12039782##29##]. The absence/presence of the 23.9 kb region, from EF1820/<italic>fsrC </italic>to EF1841, did not appear to correlate with the clinical origin of the isolates [##REF##15131223##30##]. In a more recent analysis of relationships between various <italic>E. faecalis </italic>strains, the 23.9 kb region was not detected in 86% of the strains of the clonal complex (CC)2, 58% of the CC9 strains, nor in any of the CC8 strains [##REF##17611618##31##]. The Symbioflor 1 strain, used as a probiotic, is one representative of the 7.4% of <italic>E. faecalis </italic>isolates that are missing the <italic>gelE </italic>gene in addition to the 23.9 kb region [##REF##17466591##5##,##REF##15131223##30##]. Our analysis of this area in OG1RF revealed the presence of an additional 14.8 kb fragment inserted between the corresponding EF1826 and EF1827 of OG1RF (confirmed by PCR; results not shown). In OG1RF, this 14.8 kb region contains two loci, a WxL locus (described below) and a seven-gene locus that may encode a possible ABC transporter with similarity to one annotated in <italic>Pediococcus pentosaceus</italic>.</p> <title>Components of the cell surface</title> <p>It has been shown in <italic>E. faecalis </italic>that at least one cell surface protein (Ace) is subject to domain variation [##REF##10948146##20##] and it has been postulated that domain variation may help bacteria escape the immune system. We found more polymorphisms in two families of <italic>E. faecalis </italic>proteins present on the cell surface: the MSCRAMMs and the WxL domain surface proteins. The MSCRAMMs are composed of two large regions, namely, the non-repeat A region (which is usually the ligand binding region for extracellular matrix molecules such as collagen or fibrinogen) and the B region (which typically contains repeated sub-domains). The B region of Ace contains five repeats in OG1RF, while it contains only four in V583 [##REF##10948146##20##]. We found two other MSCRAMM proteins that show polymorphisms in the number of their B repeats. OG1RF_0186 (corresponding to EF2505 of V583) has four repeats compared to seven in V583, and OG1RF_0165 (corresponding to EF2224 of V583) has eight repeats compared to five in V583. It has been proposed that the repeats are used as a stalk that projects the A region across the peptidoglycan and away from the cell surface [##REF##9799504##32##]. A hypothesis that the number of repeats may be proportional to the depth of the peptidoglycan has been proposed [##REF##9799504##32##]. However, OG1RF_0186 carries fewer repeats than EF2505 while Ace and OG1RF_0165 carry more repeats than their counterparts in V583, suggesting that our observation does not fit this hypothesis or that the peptidoglycan depth is not uniform. Apart from these three MSCRAMMs with B-repeat polymorphisms, we identified two unique MSCRAMM proteins in OG1RF: a homologue of EF0089 (OG1RF_0063, which shares 48% similarity) and a homologue of EF1896 (OG1RF_0039, which shares 75% similarity); both are located in the approximately 49 kb region unique to OG1RF described below (Figure ##FIG##0##1##; Additional data file 1).</p> <p>Another family of <italic>E. faecalis </italic>surface proteins includes the newly described WxL domain surface proteins. Siezen <italic>et al</italic>. [##REF##16723015##33##] reported a novel gene cluster encoding exclusively cell-surface proteins that is conserved in a subgroup of Gram-positive bacteria. Each gene cluster has at least one member of three gene families: a gene encoding a small LPxTG protein (approximately 120 amino acids); a gene encoding a member of the DUF916 transmembrane protein family; and a gene encoding a WxL domain surface protein. In addition, members of these gene families were found as singletons or associated with genes encoding other proteins (Additional data file 2). Recently, it was shown that the WxL domain attaches to the peptidoglycan on the cell surface [##REF##16963569##34##] and one member of this WxL domain family, the homologue of EF2686 in OG1RF (a probable internalin protein), was shown to be important for virulence in a mouse peritonitis model and is required for dissemination to the spleen and liver [##REF##17620355##35##]. OG1RF shares five complete WxL loci with V583 (EF0750-7, EF2682-6, EF2970-68, EF3181-8, and EF3248-53). OG1RF does not contain homologues of EF2248-54 (carrying instead the <italic>iol </italic>operon), though it has a novel WxL locus within the 14.8 kb unique region upstream of the <italic>fsr </italic>locus (Additional data file 2). In addition to the variation in the number of WxL loci, we also observed polymorphisms in six of the WxL domain surface proteins. For example, OG1RF_0213 shares 88% similarity with EF3188, while OG1RF_0224, OG1RF_0225, and OG1RF_0227 share 64-68% similarity with their V583 counterparts, EF3248, EF3250, and EF3252, respectively. Also, in place of EF3153, EF3154, and EF3155 (which share 70% similarity among themselves), were found non-distantly related homologues, OG1RF_0209 and OG1RF_0210, which share 60-80% similarity with EF3153, EF3154, and EF3155. It is interesting to note that while several of these WxL loci, including the EF0750 and EF3184 loci, were present by hybridization in all the strains (clinical or food isolates) tested by Lepage <italic>et al</italic>. [##REF##16980489##36##], other loci, including the EF3153 and EF3248 loci, were not detected in the majority of these strains. In addition, it appears that the EF3248 locus diverges in the Symbioflor 1 strain. When compared to V583, the sequence identity in this area between the two strains appears to be as low as 75% (depicted in Figure ##FIG##1##2## from reference [##REF##17466591##5##]).</p> <p>However, because the Symbioflor 1 genome sequence is not currently available, it was not possible to compare their respective sequences in more detail. Since these proteins are located at the surface of the cell, the low level of homology shared between them may be the result of antigenic variation. More analyses are required for a better understanding of the number, frequency and function of these WxL domain proteins and their possible relationship with the diversity of <italic>E. faecalis</italic>.</p> <p>Finally, as previously found using PCR, the <italic>cpsCDEFGHIJK </italic>operon capsule polysaccharide genes [##REF##15184433##37##] were confirmed here as missing, although OG1RF carries the <italic>cpsA </italic>and <italic>cpsB </italic>genes, which were proposed to be essential for <italic>E. faecalis </italic>since all strains tested by Hufnagel <italic>et al</italic>. [##REF##15184433##37##] carry these two genes. In OG1RF, the region that would encode the <italic>cps </italic>operon is only 59 bp in length and has no homology with V583. Thus, while V583 and OG1RF share much similarity between their surface components, there are unique differences that could potentially be important in affecting the behavior of the strains and might be useful for strain typing.</p> <title>Two-component regulatory systems</title> <p>OG1RF lacks four two-component systems found in V583. These are histidine kinase-response regulator (HK-RR)08, HK-RR12 located in the PAI, HK-RR16 and the <italic>vanB </italic>regulatory system HK-RR11 [##REF##12374813##38##]. However, an OG1RF-unique two-component system with high homology with the <italic>van</italic><sub>G </sub>locus was found at the location corresponding to the region between EF2860 and EF2861 in V583 (Table ##TAB##2##3##). OG1RF_0193 shares 82% similarity with VanR<sub>G </sub>and 81% similarity with VanR<sub>G2</sub>. Similarly, OG1RF_0192 shares 68% similarity with VanS<sub>G </sub>and VanS<sub>G2</sub>. A gene (OG1RF_0191) encoding an M15 family muramoyl pentapeptide carboxypeptidase is located downstream of these two-component regulatory genes (Figure ##FIG##3##4a##). The predicted carboxypeptidase (OG1RF_0191) shares 69% similarity over 179 amino acids with EF2297, a membrane-associated D, D-carboxypeptidase encoded by the <italic>vanB </italic>operon in V583. However, OG1RF_0191 lacks an identifiable transmembrane domain that is important to the VanY function and it is likely, therefore, that this protein may be a soluble D, D-carboxypeptidase/transpeptidase as seen in <italic>Streptomyces </italic>[##REF##3997832##39##] and <italic>Actinomadura </italic>[##REF##15987687##40##], and thus may not be involved in peptidoglycan metabolism. Consequently, it seems unlikely that this operon is a remnant of a vancomycin resistance operon in OG1RF, but rather part of a still unknown regulatory pathway.</p> <title>The <italic>iol </italic>operon</title> <p>OG1RF carries an <italic>iol </italic>operon while V583 does not. This operon encodes the factors necessary for the degradation of myo-inositol into glyceraldehyde-3P. Many soil and plant micro-organisms, including <italic>Bacillus subtilis </italic>[##REF##9226270##41##] (first <italic>iol </italic>operon identified), <italic>Klebsiella </italic>spp. [##REF##977784##42##], and cryptococci [##REF##394818##43##], have been reported to use myo-inositol as a sole carbon source. Myo-inositol, one of the nine isomers of the inositol group, belongs to the cyclitol group and is abundant in nature, particularly in the soil. The OG1RF <italic>iol </italic>operon appears to be closely related to ones described in <italic>Clostridium perfringens </italic>[##REF##15183876##44##] and <italic>Lactobacillus casei </italic>[##REF##17449687##45##]. In <italic>L. casei</italic>, the myo-inositol operon consists of ten genes with an upstream divergent regulator gene, <italic>iolR</italic>. In OG1RF, the operon appears to include ten genes, beginning with a probable transcriptional regulator (helix-turn-helix domain protein). Also, the OG1RF operon carries two copies of an <italic>iolG</italic>-like gene, which encodes inositol 2-dehydrogenase, the first enzyme of the myo-inositol degradation pathway (Figure ##FIG##4##5##). However, the order of the genes is not the same between <italic>E. faecalis </italic>and <italic>L. casei</italic>. In addition, <italic>iolH</italic>,<italic>iolJ </italic>and <italic>iolK</italic>, present in <italic>L. casei</italic>, are not present in OG1RF, nor are <italic>iolH </italic>and <italic>iolK </italic>present in the <italic>C. perfringens iol </italic>operon.</p> <p>Yebra <italic>et al</italic>. reported that <italic>L. casei </italic>was the sole member of the Lactobacillales to carry a functional <italic>iol </italic>operon [##REF##17449687##45##]. To survey <italic>E. faecalis</italic>, also a member of this order, for the presence of the <italic>iol </italic>operon, 48 isolates with different MLST and/or from various origins (including OG1RF and V583) were tested for myo-inositol fermentation; 23 of 48 isolates were positive. In addition, PCR verified the presence of <italic>iolE </italic>and <italic>iolR </italic>in these strains and in one negative for myo-inositol fermentation, indicating that the <italic>iol </italic>operon is not unique to OG1RF. To verify that the <italic>iol </italic>genes are responsible for the fermentation of myo-inositol in OG1RF, transposon insertion mutants [##REF##122512##9##] in the <italic>iolB </italic>and <italic>iolG2 </italic>genes of OG1RF were tested. Both mutants failed to ferment myo-inositol (data not shown), demonstrating that these genes are essential for myo-inositol fermentation. To investigate whether the <italic>iol </italic>operon was 'inserted into' or 'removed from' a putative ancestral strain, the sequences surrounding the <italic>iol </italic>genes were examined. In OG1RF, the <italic>iol </italic>operon is located between the equivalent of EF2239 and EF2352 when compared with V583. In V583, this region encodes probable prophage proteins and carries the <italic>vanB </italic>transposon, which confers vancomycin resistance. Since we did not identify any remnants of the <italic>iol </italic>operon in V583, it would appear that at least two independent events at the same location differentiate OG1RF and V583, suggesting that it is a hot region for rearrangement. This region between EF2239 and EF2352 (111 Kb) is also missing in the Symbioflor 1 strain (referred to as gap 2) [##REF##17466591##5##]. The possible junction and presence of unique sequence in this region, if investigated, was not mentioned in the publication. Nonetheless, preliminary analysis of other strains' genotypes in this area seemed to confirm the hypothesis of a hot region for rearrangement (data not shown).</p> <title>A homologue of Tn<italic>916 </italic>in OG1RF</title> <p>An analysis of the G+C content of OG1RF unique regions revealed several loci with a lower G+C content than the 37.8% average content of OG1RF. One of these is an approximately 49 kb fragment with a G+C content of 32.1% located between an rRNA operon and the homologue in OG1RF of EF1053, replacing 10 tRNA genes present in V583 (Figure ##FIG##0##1##). This fragment appears to be a patchwork composed of hypothetical genes, homologues of Tn<italic>916</italic>-associated genes and homologues of genes from other Gram-positive organisms, including <italic>Listeria</italic>, <italic>E. faecium</italic>, staphylococci, or lactococci (Additional data file 1). It is interesting to note that this region contains: a putative adhesin protein gene (OG1RF_0039) at one end of the fragment; homologues of 14 Tn<italic>916</italic>-associated genes (Tn<italic>916</italic>_2 to Tn<italic>916</italic>_12, Tn<italic>916</italic>_18 and Tn<italic>916</italic>_19, with an average of 70% similarity); and a gene encoding a putative integrase (OG1RF_0088) at the other end - these three features are also present in Tn<italic>5386 </italic>in <italic>E. faecium </italic>D344R [##REF##17408741##46##]. However, the approximately 49 kb fragment lacks an excisase gene and the probable lantibiotic ABC transporter genes present in Tn<italic>5386</italic>.</p> <title>An uninterrupted competence operon in OG1RF</title> <p>OG1RF contains what appears to be an intact competence operon while that of V583 appears to be non-functional. This operon in OG1RF is similar to a nine-gene operon described in <italic>Streptococcus mutans </italic>[##REF##15632435##47##], as shown in Figure ##FIG##5##6##. For example, the homologue in OG1RF of EF2046 shares 61% similarity with ComYA and the OG1RF homologue of EF2045 is 55% similar to ComYB. In <italic>S. mutans</italic>, only the first seven genes of the operon are essential for competence [##REF##15632435##47##]. In V583, the fourth gene of this operon (corresponding to OG1RF_0148) is interrupted by phage 4 (EF1896-EF2043); in addition, EF1984 contains a premature stop codon not found in the corresponding gene in OG1RF (OG1RF_0228).</p> <p>Natural competence has not been reported for <italic>E. faecalis</italic>. To assess the functionality of this operon in OG1RF, we evaluated the competence of cells in different phases of growth (early log growth to stationary phase) using pAM401 [##REF##3005240##48##] and pMSP3535VA [##REF##10964628##49##]. We were not able to show natural competence under the conditions tested. We have also noted that V583 is less transformable by electroporation than OG1RF. To investigate the possibility that directly or indirectly the <italic>com </italic>operon might be responsible for this phenotype, we also evaluated transformability by electroporation. When compared with OG1RF, transposon mutants [##REF##15489440##12##] in the OG1RF equivalent of EF2045 (encoding the <italic>comGB </italic>homologue) and in the OG1RF equivalent of EF1986 (encoding the <italic>comGF </italic>homologue) showed similar levels of transformability by electroporation (data not shown), implying that the difference in electroporation efficiency observed between OG1RF and V583 is not related to this locus.</p> <p>In <italic>Streptococcus pneumoniae </italic>[##REF##9701804##50##], the competence operon is tightly regulated by a quorum sensing two-component system (ComDE) and a competence-stimulating peptide (CSP; encoded by <italic>comC</italic>). We did not find any homologues of CSP in OG1RF. Two homologous ComDE sensor histidine kinase/response regulators were found in OG1RF, one of which is FsrC/FsrA. Based on our previous microarray data, the Fsr system does not regulate the <italic>comY </italic>operon, at least under our previously used conditions (mid-log phase growth to early stationary phase in brain heart infusion (BHI)) [##REF##16585749##22##]. The other ComDE homology is that with a two-component system unique to OG1RF (OG1RF_0199 and OG1RF_0198, respectively) that lies on a 4,706 bp unique fragment that maps between EF3114 and EF3115 in V583. This fragment also carries two genes (OG1RF_0200 and OG1RF_0201) encoding homologues of the YhaQ and YhaP sodium efflux ATP-binding cassette efflux/transporter proteins (Figure ##FIG##3##4b##). Although they are potential elements of a secretion apparatus, these two proteins do not share any homology at the protein level with the competence secretion apparatus ComAB of <italic>S. pneumoniae </italic>[##REF##11115120##51##] nor CslAB from <italic>S. mutans </italic>[##REF##11154426##52##]. Searching for a possible CSP in the vicinity of these genes, we identified a small ORF encoding 50 amino acids between <italic>yhaP </italic>and OG1RF199 and another encoding 20 amino acids downstream of OG1RF198. More analysis will be required to determine if there are conditions in which the OG1RF <italic>com </italic>operon is expressed and to determine whether or not this two-component system is involved in competence.</p> <title>Limited presence of mobile elements</title> <p>By probing a microarray of the V583 genome and plasmids with OG1RF genomic DNA, we previously estimated that only 75% of V583 ORFs were also present in OG1RF [##REF##16585749##22##]. Later, Aakra <italic>et al</italic>. [##REF##17220255##53##] compared nine strains, including OG1RF to V583, using comparative genomic hybridization. In these results, OG1RF appears to carry a few genes included in the PAI, and a few prophage genes. Using the complete genome sequence, we have now found that OG1RF lacks 639 genes and the three pTET plasmids described in V583. All but 45 of the missing genes are associated with putative mobile elements, such as the entire PAI, the recently described phages 1, 3, 4, 5, 6, and 7 [##REF##17611618##31##], and the approximately 111 kb area between genes EF2240 and EF2351 (including the <italic>vanB </italic>transposon) present in V583. The absence of these elements appears also to be a characteristic of the commensal strain Symbioflor 1, although because the genome was not completely finished, the possibility remains that some of these regions were not sequenced. In conclusion, other than the approximately 49 kb fragment containing a Tn<italic>916 </italic>homologue, it appears that OG1RF has only one additional possible mobile element derivative, namely the phage 2 proposed to be part of the core genome [##REF##17611618##31##].</p> <title>Fusidic acid and rifampicin resistance</title> <p>OG1RF was sequentially selected from OG1 for resistance to fusidic acid and rifampicin [##REF##122512##9##]. The mutation leading to rifampicin resistance was identified in the <italic>rpoB </italic>gene by Ozawa <italic>et al</italic>. [##REF##15907539##54##] and is caused by the A1467G mutation, which results in substitution of arginine for histidine at amino acid 489. The mutation also affected the clumping phenotype of <italic>traA </italic>mutants and this effect appears to be specific for the pAM373 system [##REF##15907539##54##]. All of the other 22 differences in <italic>rpoB </italic>between OG1RF and V583 are synonymous. Fusidic acid resistance is associated with mutation(s) in the <italic>fusA </italic>gene, which encodes elongation factor G. We compared <italic>fusA </italic>from OG1RF with that in V583 and identified two differences (C1368A and T1992C). The mutation T1992C is synonymous, while the mutation C1368A leads to the presence of glutamine (histidine in V583) at position 404 in OG1RF. Mutations in this region have been associated with fusidic acid resistance in <italic>Staphylococcus aureus </italic>[##REF##11036042##55##,##REF##11309125##56##], and thus the C1368A mutation is likely the cause of the fusidic acid resistance phenotype in OG1RF.</p> <title>Virulence and biofilm comparisons of OG1RF with V583</title> <p>When compared in the mouse peritonitis model, the LD<sub>50 </sub>values of V583 in different determinations were lower (4.8 × 10<sup>7 </sup>to 1.1 × 10<sup>8 </sup>colony forming units (CFU)/ml) than the LD<sub>50 </sub>values of OG1RF (1.2 × 10<sup>8 </sup>to 4.8 × 10<sup>8 </sup>CFU/ml). However, at comparable inoculum, OG1RF (4 × 10<sup>8 </sup>CFU/ml) showed more rapid mortality versus V583 (5 × 10<sup>8 </sup>CFU/ml) in the first 48 hours (<italic>P </italic>= 0.0034; Additional data file 3). In a urinary tract infection model administering mixed equal inocula of V583 and OG1RF, OG1RF significantly outnumbered V583 in kidney with geometric means of 1.3 × 10<sup>4 </sup>CFU/gm for OG1RF versus 1.9 × 10<sup>2 </sup>CFU/gm for V583 (<italic>P </italic>= 0.0005); in urinary bladder homogenates, the geometric mean CFU/gm was 1.7 × 10<sup>3 </sup>for OG1RF versus 6.6 × 10<sup>1 </sup>for V583 (<italic>P </italic>= 0.003; Figure ##FIG##6##7a##). Similarly, in mono-infection, the geometric mean CFU/gm of OG1RF in kidneys was 9.4 × 10<sup>3 </sup>versus 4 × 10<sup>1 </sup>for V583 (<italic>P </italic>= 0.0035; Figure ##FIG##6##7b##). We also found that OG1RF produced 20% more biofilm (<italic>P </italic>&lt; 0.01) than V583 at 24 hours (results not shown). These results, together with the previous results in <italic>C. elegans </italic>[##REF##11535834##17##], demonstrate that OG1RF, although lacking what was thought to be important for virulence (PAI, plasmids, prophages), is as pathogenic as V583 in at least three assays.</p>
<title>Conclusion</title> <p><italic>E. faecalis </italic>OG1RF carries a number of unique loci compared to V583. Those of particular interest include new surface proteins (MSCRAMMs and WxL domain proteins), an operon encoding myo-inositol utilization, an intact competence operon, and two CRISPR elements. The CRISPR elements may be of particular significance when one considers that most of what is missing from OG1RF compared to V583 consists of mobile genetic elements (MGEs), including 6 phages or remnants thereof. The presence of the CRISPR elements in OG1RF provides a tantalizing, but as yet unproven, explanation for the discordance in the number of mobile elements between these two strains.</p> <p>The acquisition of MGE is believed to be an important mechanism by which the species <italic>E. faecalis </italic>had been able to generate genetic diversity and, thereby, highly variable phenotypes [##REF##12663927##4##]. It has been proposed that the ability of <italic>E. faecalis </italic>to cause healthcare related infections is associated with these MGEs [##REF##12663927##4##,##REF##17466591##5##]. This hypothesis was supported by several studies that have highlighted the importance of virulence determinants carried by these mobile elements, such as cytolysin [##REF##15638770##57##] by the PAI. However, more recent results from Aakra <italic>et al</italic>. [##REF##17220255##53##] and Lepage <italic>et al</italic>. [##REF##16980489##36##] demonstrate that these factors may be present in harmless strains while absent in clinical isolates, indicating that <italic>E. faecalis </italic>virulence is not dependent on any single virulence factor. Indeed, few studies have compared the virulence pattern of strains from various origins. The increased ability of OG1RF to infect kidneys and to produce biofilm, despite the absence of MGEs and their associated virulence factors, was surprising. Different possibilities can be proposed relating to the factors important for these differences in enterococcal infections. One of these is that virulence in the assays used may be linked to the shared core genome of these two strains, with the differences arising from the unique genes. On the other hand, virulence could be associated primarily with the genes unique to each strain, but with each set being able to complement the absence of the other. It seems most likely that virulence, and/or some combination of virulence and fitness, is caused by the expression of a mixture of both the core and unique genes. It is also important to remember that <italic>E. faecalis </italic>is a well adapted commensal, carrying the genes necessary to survive and to colonize the gut, and that a subset, particularly MLST CC2 and CC9 [##REF##16880002##58##], predominates among hospital acquired infections. It may be that these clonal complexes are not more virulent <italic>per se</italic>, as defined in the assays described here, but rather are better able to survive and/or colonize hospitalized patients, taking advantage of factors that predispose to nosocomial infections such as urinary or venous catheters, or mucositis, among others.</p> <p>Sequencing of more <italic>E. faecalis </italic>strains may facilitate our understanding of the path from commensalism to pathogenicity, a crucial prerequisite for designing therapeutic interventions directed to control an organism that is already resistant to a large spectrum of antibiotics.</p>
<p>This is an open access article distributed under the terms of the Creative Commons Attribution License (<ext-link ext-link-type="uri" xlink:href="http://creativecommons.org/licenses/by/2.0"/>), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.</p> <p>A comparison of two strains of the hospital pathogen <italic>Enterococcus faecalis</italic> suggests that mediators of virulence differ between strains and that virulence does not depend on mobile gene elements</p> <title>Background</title> <p><italic>Enterococcus faecalis </italic>has emerged as a major hospital pathogen. To explore its diversity, we sequenced <italic>E. faecalis </italic>strain OG1RF, which is commonly used for molecular manipulation and virulence studies.</p> <title>Results</title> <p>The 2,739,625 base pair chromosome of OG1RF was found to contain approximately 232 kilobases unique to this strain compared to V583, the only publicly available sequenced strain. Almost no mobile genetic elements were found in OG1RF. The 64 areas of divergence were classified into three categories. First, OG1RF carries 39 unique regions, including 2 CRISPR loci and a new WxL locus. Second, we found nine replacements where a sequence specific to V583 was substituted by a sequence specific to OG1RF. For example, the <italic>iol </italic>operon of OG1RF replaces a possible prophage and the <italic>vanB </italic>transposon in V583. Finally, we found 16 regions that were present in V583 but missing from OG1RF, including the proposed pathogenicity island, several probable prophages, and the <italic>cpsCDEFGHIJK </italic>capsular polysaccharide operon. OG1RF was more rapidly but less frequently lethal than V583 in the mouse peritonitis model and considerably outcompeted V583 in a murine model of urinary tract infections.</p> <title>Conclusion</title> <p><italic>E. faecalis </italic>OG1RF carries a number of unique loci compared to V583, but the almost complete lack of mobile genetic elements demonstrates that this is not a defining feature of the species. Additionally, OG1RF's effects in experimental models suggest that mediators of virulence may be diverse between different <italic>E. faecalis </italic>strains and that virulence is not dependent on the presence of mobile genetic elements.</p>
<title>Abbreviations</title> <p>ATCC, American type culture collection; BHI, brain heart infusion; CC, clonal complex; CFU, colony forming units; CRISPR, comprised of regularly interspaced short palindromic repeats; CSP, competence-stimulating peptide; HK-RR, histidine kinase-response regulator; MGE, mobile genetic element; MSCRAMM, microbial surface component recognizing adhesive matrix molecules; MLST, multilocus sequence typing; ORF, open reading frame; PAI, pathogenicity island.</p> <title>Authors' contributions</title> <p>GMW, DAG, and BEM designed the study. AB performed much of the post-annotation analysis and non-animal experiments, and wrote the draft of the manuscript. KVS performed the animal experiments. AB, DAG, XQ, JS, SY, AM, KAF, JG, CAA, YS, SRN, MZ, VPP, SC, and SKH annotated the genome. XQ and HJ contributed bioinformatics support. YD, SD-R, CB, HS, GC, GW, DM, LC, and RAG composed the sequencing and finishing team. DAG, BEM, SKH, and GMW assisted in critical review of the manuscript. All authors read and approved the final manuscript.</p> <title>Additional data files</title> <p>The following additional data files are available with the online version of this paper. Additional data file ##SUPPL##0##1## is a list of the ORFs unique to OG1RF compared to V583 with their OG1RF locus tag, location in the genome, and definition. Additional data file ##SUPPL##1##2## is a list of genes encoding proteins with a WxL domain in OG1RF and/or V583. Additional data file ##SUPPL##2##3## shows the results of the mouse peritonitis model using OG1RF and V583, with the statistical analysis. Additional data file ##SUPPL##3##4## is a list of the significant primers used in this study.</p> <title>Supplementary Material</title>
<title>Acknowledgements</title> <p>We would like to express our sincere thanks to J Hernandez, S Wang, Z Li, D Ngo, and L Hemphill for their help during the sequencing process. We also would like to thank JM Urbach, Massachusetts General Hospital, Boston, MA for helping localize the transposon insertions in the OG1RF unique sequences. This research was supported by grant R21 AI064470 from the National Institutes of Health to GMW and by NIH grant R37 AI47923 from the Division of Microbiology and Infectious Diseases, NIAID, to BEM.</p>
<fig position="float" id="F1"><label>Figure 1</label><caption><p>Map of the OG1RF chromosome. The following features are displayed (from the inside out): restriction maps using <italic>Sfi</italic>I, <italic>Asc</italic>I, and <italic>Not</italic>I (black) from Murray <italic>et al</italic>. [##REF##8349561##10##] overlaid with the digestion profile predicted from the sequence (red); G+C content in percentage in green; the total OG1RF-unique genes are shown in purple with those in (+) orientation labeled in blue, and those in (-) orientation labeled in red.</p></caption></fig> <fig position="float" id="F2"><label>Figure 2</label><caption><p>Dot plot of OG1RF versus V583 generated by BLASTN. The dot plot was generated by aligning the OR1RF genome against the V583 genome using BLASTN (e-value 1e-10). The alignment pairs were plotted according to their genome coordinates. The visible areas of divergences are labeled using 'Δ ' to indicate a sequence absent in OG1RF and '∇ ' to indicate a sequence unique to OG1RF (locus tag OG1_xxxx) when compared with V583 (locus tag EFxxxx). Phages 1, 3, 4, 5, 6, 7 of V583 (φ1 to 7; see [##REF##17611618##31##]) and the PAI locations, all of which are missing from OG1RF, are also indicated.</p></caption></fig> <fig position="float" id="F3"><label>Figure 3</label><caption><p>The two CRISPR loci of OG1RF. <bold>(a)</bold> The CRISPR1 locus. The CRISPR1 element is represented with a hatched box while the CRISPR1 associated genes are represented in orange; the white arrows indicate ORFs present in both OG1RF and V583. The black diamonds represent the 37 bp repeat sequences, while the open boxes with a number indicate the 29 bp unique sequences. <bold>(b)</bold> The CRISPR2 locus containing only a CRISPR element. <bold>(c)</bold> CRISPR consensus and unique sequences. The underlined bases indicate mismatches at these locations. The sequences numbered 1 to 14 represent the unique sequences located in the CRISPR1 and CRISPR2 elements.</p></caption></fig> <fig position="float" id="F4"><label>Figure 4</label><caption><p>Two-component systems unique to OG1RF. <bold>(a)</bold> Two-component system with homology to the Van<sub>G </sub>system. <bold>(b) </bold>Two-component system with homology to the <italic>comCD </italic>genes of <italic>S. pneumoniae</italic>. The two-component system (OG1RF_0198 and OG1RF_0199) is indicated in light blue; the two ORFs encoding potential transporter proteins (OG1RF_0200 and OG1RF_0201) are represented in pink. In green are indicated two small ORFs encoding polypeptides of less than 51 amino acids. The white arrows indicate ORFs also present in V583.</p></caption></fig> <fig position="float" id="F5"><label>Figure 5</label><caption><p>The <italic>iol </italic>operon. The <italic>iol </italic>genes are labeled based on the homology/conserved motif of their encoded proteins with known enzymes necessary for myo-inositol degradation. For all strains, the described or probable regulator is represented in blue. <italic>E. faecalis </italic>OG1RF: the <italic>iol </italic>operon is represented in yellow, OG1RF_0166 (green arrow) located downstream of the <italic>iol </italic>operon encodes a probable PTS IIC component, while the white arrows indicate ORFs also present in V583. For <italic>B. subtilis </italic>168, <italic>C. perfringens </italic>strain 13, and <italic>L. casei </italic>BL23, the <italic>iol </italic>genes are represented in green, orange and purple, respectively. <italic>C. perfringens iol </italic>mRNA transcript includes five other genes encoding proteins whose functions do not appear to be related to myo-inositol degradation; these genes are represented in gray.</p></caption></fig> <fig position="float" id="F6"><label>Figure 6</label><caption><p>The OG1RF competence operon and its similarity with the competence operon of <italic>S. mutans</italic>. The ORFs essential for natural competence in <italic>S. mutans </italic>are shown in green as well as their homologues in OG1RF and V583. The ORF corresponding to the homologue of ComYD was not described in V583 [##REF##12663927##4##], due to the presence of a probable prophage (EF1986-EF2043). The premature stop codon in EF1984 in V583 is indicated with an asterisk. <italic>ackA/</italic>EF1983 is represented in orange. The proteins encoded by the ORFs represented in white do not share any features of the known competence proteins or homology between <italic>S. mutans </italic>and <italic>E. faecalis</italic>; in <italic>S. mutans</italic>, <italic>ackA </italic>and <italic>ytxK </italic>are co-transcribed with the <italic>comY </italic>genes [##REF##15632435##47##].</p></caption></fig> <fig position="float" id="F7"><label>Figure 7</label><caption><p>Comparison of OG1RF and V583 in a mouse urinary tract infection model. <bold>(a)</bold> Mixed infection by wild-type <italic>E. faecalis </italic>strains OG1RF and V583 in the kidneys and urinary bladders of mice (n = 21; competition assay). Data are expressed as the log<sub>10</sub>(CFU)/gm for OG1RF or V583; the log<sub>10</sub>(CFU)/gm for both kidneys were combined and averaged from two independent experiments. Black solid diamonds and triangles represent <italic>E. faecalis </italic>strains OG1RF and V583, respectively, for kidney homogenates, and empty diamonds and triangles represent OG1RF and V583, respectively, for urinary bladder homogenates. Horizontal bars represent geometric means. Log<sub>10</sub>(CFU) were compared for statistical significance by the paired <italic>t</italic>-test. The minimum detection limit in these experiments was 10<sup>1 </sup>and 10<sup>2 </sup>CFU/gm of kidney and urinary bladder homogenates, respectively. <bold>(b) </bold>Mono-infection using <italic>E. faecalis </italic>strains OG1RF or V583 in the kidneys of mice (10<sup>3 </sup>CFU per mice, n = 9). Data are expressed as log<sub>10</sub>(CFU)/gm for OG1RF recovered from kidney homogenates 48 h after infection; the log<sub>10</sub>(CFU)/gm for a kidney pair were combined and averaged. Black and white triangles represent OG1RF and V583, respectively. Horizontal bars represent geometric means. The CFU of V583 recovered from kidneys was significantly reduced compared to the CFU of OG1RF, as determined by the unpaired <italic>t</italic>-test.</p></caption></fig>
<table-wrap position="float" id="T1"><label>Table 1</label><caption><p>General features of OG1RF compared to V583</p></caption><table frame="hsides" rules="groups"><thead><tr><td/><td align="left">V583</td><td/><td align="left">OG1RF</td></tr></thead><tbody><tr><td align="left">General features</td><td/><td/><td/></tr><tr><td align="left">Size (base pairs)</td><td align="left">3,218,031</td><td/><td align="left">2,739,633</td></tr><tr><td align="left">G+C content (%)</td><td align="left">37.5</td><td/><td align="left">37.8</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left">rRNA genes</td><td align="left">12</td><td/><td align="left">12</td></tr><tr><td align="left">tRNA genes</td><td align="left">68</td><td/><td align="left">58</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left">Genes common to both strains</td><td/><td align="left">2,474*</td><td/></tr><tr><td/><td/><td/><td/></tr><tr><td align="left">Genes unique to OG1RF</td><td/><td/><td/></tr><tr><td align="left">Similar to known proteins</td><td/><td/><td align="left">114<sup>†</sup></td></tr><tr><td align="left">Conserved hypotheticals</td><td/><td/><td align="left">50</td></tr><tr><td align="left">No database match</td><td/><td/><td align="left">63</td></tr><tr><td align="left">Total</td><td/><td/><td align="left">227</td></tr><tr><td/><td/><td/><td/></tr><tr><td align="left">Total number of ORFs</td><td align="left">3,113</td><td/><td align="left">2,701<sup>‡</sup></td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T2"><label>Table 2</label><caption><p>Frequency of the CRISPR locus among <italic>E. faecalis</italic></p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">Name</td><td align="left">Other</td><td align="left">Origin</td><td align="left">Source/reference</td><td align="left">MLST</td><td align="left">ErmR*</td><td align="left">VanR<sup>†</sup></td><td align="left">cas<sup>‡</sup></td><td align="left">EF0672-3<sup>§</sup></td></tr></thead><tbody><tr><td align="left">TX4002</td><td align="left">OG1RF</td><td align="left">Human</td><td align="left">[##REF##98769##8##,##REF##122512##9##]</td><td align="left">1</td><td align="left">-</td><td align="left">-</td><td align="left">+</td><td align="left">-</td></tr><tr><td align="left">TX2708</td><td align="left">V583</td><td align="left">Clinical isolate</td><td align="left">[##REF##2554802##6##]</td><td align="left">6<sup>¶</sup></td><td align="left">+</td><td align="left">+</td><td align="left">-</td><td align="left">+</td></tr><tr><td align="left">TX2144</td><td align="left">E1840</td><td align="left">Clinical isolate</td><td align="left">Ruiz-Garbajosa P.<sup>#</sup></td><td align="left">40</td><td align="left">+</td><td align="left">-</td><td align="left">+</td><td align="left">-</td></tr><tr><td align="left">TX2135</td><td align="left">E1795</td><td align="left">Hospital survey</td><td align="left">Ruiz-Garbajosa P.</td><td align="left">44</td><td align="left">-</td><td align="left">-</td><td align="left">-</td><td align="left">+</td></tr><tr><td align="left">TX2137</td><td align="left">E1798</td><td align="left">Hospital survey</td><td align="left">Ruiz-Garbajosa P.</td><td align="left">16</td><td align="left">+</td><td align="left">-</td><td align="left">+</td><td align="left">-</td></tr><tr><td align="left">TX2141</td><td align="left">E1825</td><td align="left">Clinical isolate</td><td align="left">Ruiz-Garbajosa P.</td><td align="left">25</td><td align="left">-</td><td align="left">-</td><td align="left">-</td><td align="left">+</td></tr><tr><td align="left">TX2140</td><td align="left">E1803</td><td align="left">Hospital survey</td><td align="left">Ruiz-Garbajosa P.</td><td align="left">38</td><td align="left">-</td><td align="left">-</td><td align="left">-</td><td align="left">+</td></tr><tr><td align="left">TX2138</td><td align="left">E1801</td><td align="left">Hospital survey</td><td align="left">Ruiz-Garbajosa P.</td><td align="left">48</td><td align="left">-</td><td align="left">-</td><td align="left">-</td><td align="left">+</td></tr><tr><td align="left">TX2146</td><td align="left">E1844</td><td align="left">Clinical isolate</td><td align="left">Ruiz-Garbajosa P.</td><td align="left">61</td><td align="left">-</td><td align="left">-</td><td align="left">-</td><td align="left">+</td></tr><tr><td align="left">TX2139</td><td align="left">E1802</td><td align="left">Hospital survey</td><td align="left">Ruiz-Garbajosa P.</td><td align="left">35</td><td align="left">+</td><td align="left">-</td><td align="left">+</td><td align="left">-</td></tr><tr><td align="left">TX4240</td><td align="left">A0826</td><td align="left">Pig</td><td align="left">Jensen L.</td><td align="left">98</td><td align="left">+</td><td align="left">-</td><td align="left">+</td><td align="left">-</td></tr><tr><td align="left">TX4247</td><td align="left">E1876</td><td align="left">Pig</td><td align="left">Gaastra W.</td><td align="left">20</td><td align="left">+</td><td align="left">-</td><td align="left">+</td><td align="left">-</td></tr><tr><td align="left">TX4245</td><td align="left">E1872</td><td align="left">Dog</td><td align="left">Gaastra W.</td><td align="left">16</td><td align="left">+</td><td align="left">-</td><td align="left">+</td><td align="left">-</td></tr><tr><td align="left">TX4243</td><td align="left">E0252</td><td align="left">Calf</td><td align="left">Mevius D.</td><td align="left">23</td><td align="left">+</td><td align="left">+</td><td align="left">-</td><td align="left">+</td></tr><tr><td align="left">TX4255</td><td align="left">A0808</td><td align="left">Clinical isolate</td><td align="left">Kawalec M.</td><td align="left">88</td><td align="left">-</td><td align="left">-</td><td align="left">-</td><td align="left">+</td></tr><tr><td align="left">TX4259</td><td align="left">A1006</td><td align="left">Clinical isolate</td><td align="left">Kawalec M.</td><td align="left">135</td><td align="left">-</td><td align="left">-</td><td align="left">-</td><td align="left">+</td></tr></tbody></table></table-wrap> <table-wrap position="float" id="T3"><label>Table 3</label><caption><p>OG1RF-unique regulators</p></caption><table frame="hsides" rules="groups"><thead><tr><td align="left">OG1RF</td><td align="left">Description</td><td align="left">Best hit</td><td align="left">Size*</td><td align="left">Comments</td></tr></thead><tbody><tr><td align="left">OG1RF_0070</td><td align="left">Transcriptional regulator</td><td align="left">116512576</td><td align="left">102</td><td align="left">-</td></tr><tr><td align="left">OG1RF_0073</td><td align="left">LytR family response regulator</td><td align="left">81428169</td><td align="left">151</td><td align="left">-</td></tr><tr><td align="left">OG1RF_0120</td><td align="left">BglG family transcriptional antiterminator</td><td align="left">47095712</td><td align="left">494</td><td align="left">Probable regulator of the downstream PTS system</td></tr><tr><td align="left">OG1RF_0138</td><td align="left">Transcriptional regulator</td><td align="left">116493423</td><td align="left">219</td><td align="left">Probable transcriptional regulator of the downstream ABC superfamily transporter</td></tr><tr><td align="left">OG1RF_0143</td><td align="left">GntR family transcriptional regulator</td><td align="left">82745913</td><td align="left">236</td><td align="left">Probable regulator of the downstream PTS system</td></tr><tr><td align="left">OG1RF_0175</td><td align="left">DNA binding protein</td><td align="left">15890504</td><td align="left">293</td><td align="left">Probable regulator of the <italic>iol </italic>operon</td></tr><tr><td align="left">OG1RF_0192</td><td align="left">Sensor histidine kinase VanS<sub>G</sub></td><td align="left">119635646</td><td align="left">371</td><td align="left">Best homology with Van<sub>G </sub>and</td></tr><tr><td align="left">OG1RF_0193</td><td align="left">Response regulator VanR<sub>G</sub></td><td align="left">119635645</td><td align="left">235</td><td align="left">Van<sub>G2 </sub>two-component systems.</td></tr><tr><td/><td/><td/><td/><td align="left">OG1RF_0192 and OG1RF_0193 appear cotranscribed with a gene encoding a M15 family muramoylpentapeptide carboxypeptidase</td></tr><tr><td align="left">OG1RF_0198</td><td align="left">Response regulator</td><td align="left">47567135</td><td align="left">240</td><td align="left">Best homology with AgrA from <italic>Bacillus cereus </italic>G9241. However, no presence of AgrB or AgrD homologues in the vicinity. Also similar to ComE of <italic>S. pneumoniae </italic>(52% similarity)</td></tr><tr><td align="left">OG1RF_0199</td><td align="left">Sensor histidine kinase</td><td align="left">47567134</td><td align="left">443</td><td align="left">Best homology with AgrC from <italic>Bacillus cereus </italic>G9241. Also similar to ComD of <italic>Streptococcus pneumoniae </italic>(48% similarity)</td></tr><tr><td align="left">OG1RF_0220</td><td align="left">Probable endoribonuclease MazF</td><td align="left">69244828</td><td align="left">114</td><td align="left">Toxin-antitoxin described in <italic>E.</italic></td></tr><tr><td align="left">OG1EF_0221</td><td align="left">Probable antitoxin MazE</td><td align="left">69244829</td><td align="left">77</td><td align="left"><italic>coli </italic>and recently on an <italic>E. faecium </italic>plasmid</td></tr></tbody></table></table-wrap>
<supplementary-material content-type="local-data" id="S1"><caption><title>Additional data file 1</title><p>OG1RF locus tag, location in the genome, and definition are given.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S2"><caption><title>Additional data file 2</title><p>Genes encoding proteins with a WxL domain in OG1RF and/or V583.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S3"><caption><title>Additional data file 3</title><p>Results of the mouse peritonitis model using OG1RF and V583, with the statistical analysis.</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="S4"><caption><title>Additional data file 4</title><p>The significant primers used in this study.</p></caption></supplementary-material>
<table-wrap-foot><p>*The assessment of the genes common to both strains is based on the homology at the DNA level with the ORFs described for V583 (source TIGR [##UREF##2##70##]). The BLASTN cutoff e-value was 1e-5. <sup>†</sup>This number includes the proteins with domain polymorphism (see text for details). <sup>‡</sup>Estimated number of ORFs calculated by adding the OG1RF-unique ORFs to the number of ORFs shared with V583.</p></table-wrap-foot> <table-wrap-foot><p>*Erythromycin resistance was tested at 5 μg/ml. <sup>†</sup>Vancomycin resistance was tested at 10 μg/ml. <sup>‡</sup>Two sets of primers were used to detect the <italic>cas </italic>genes (<italic>cas1 </italic>and <italic>csn1</italic>). <sup>§</sup>This set of primers amplifies the junction between EF0672 and EF0673 where the CRISPR1 locus is inserted in OG1RF. <sup>¶</sup>CC2. <sup># </sup>Ruiz-Garbajosa P. (Spain), Jensen L. (Denmark), Gaastra W. and Mevius D. (Netherland), and Kawalec M. (Poland).</p></table-wrap-foot> <table-wrap-foot><p>*Amino acids</p></table-wrap-foot>
<graphic xlink:href="gb-2008-9-7-r110-1"/> <graphic xlink:href="gb-2008-9-7-r110-2"/> <graphic xlink:href="gb-2008-9-7-r110-3"/> <graphic xlink:href="gb-2008-9-7-r110-4"/> <graphic xlink:href="gb-2008-9-7-r110-5"/> <graphic xlink:href="gb-2008-9-7-r110-6"/> <graphic xlink:href="gb-2008-9-7-r110-7"/>
<media xlink:href="gb-2008-9-7-r110-S1.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-7-r110-S2.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-7-r110-S3.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media> <media xlink:href="gb-2008-9-7-r110-S4.pdf" mimetype="application" mime-subtype="pdf"><caption><p>Click here for file</p></caption></media>
[{"article-title": ["Genboree"]}, {"surname": ["Reed", "Muench"], "given-names": ["L", "H"], "article-title": ["A simple method of estimating fifty per cent end points."], "source": ["Am J Hygiene"], "year": ["1938"], "volume": ["27"], "fpage": ["493"], "lpage": ["497"]}, {"article-title": ["TIGR"]}]
Genome Biol. 2008 Jul 8; 9(7):R110