- Extended genome report
- Open Access
Complete genome sequence of Staphylococcus aureus, strain ILRI_Eymole1/1, isolated from a Kenyan dromedary camel
Standards in Genomic Sciencesvolume 10, Article number: 109 (2015)
We report the genome of a Staphylococcus aureus strain (ILRI_Eymole1/1) isolated from a nasal swab of a dromedary camel (Camelus dromedarius) in North Kenya. The complete genome sequence of this strain consists of a circular chromosome of 2,874,302 bp with a GC-content of 32.88 %. In silico annotation predicted 2755 protein-encoding genes and 76 non-coding genes. This isolate belongs to MLST sequence type 30 (ST30). Phylogenetic analysis based on a subset of 283 core genes revealed that it falls within the human clonal complex 30 (CC30) S. aureus isolate cluster but is genetically distinct. About 79 % of the protein encoding genes are part of the CC30 core genome (genes common to all CC30 S. aureus isolates), ~18 % were within the variable genome (shared among multiple but not all isolates) and ~ 3 % were found only in the genome of the camel isolate. Among the 85 isolate-specific genes, 79 were located within putative phages and pathogenicity islands. Protein encoding genes associated with bacterial adhesion, and secretory proteins that are essential components of the type VII secretion system were also identified. The complete genome sequence of S. aureus strain ILRI_Eymole1/1 has been deposited in the European Nucleotide Archive under the accession no LN626917.1.
S. aureus is a bacterial species which has been isolated from diverse hosts including humans, other mammals and birds [1, 2]. In humans, it is persistently present in the nares of approximately 20 % of all individuals and intermittently carried by nearly 30 % individuals . S. aureus has been reported to be a common cause of wound infections, pneumonia and bacteraemia in humans in Kenya [4, 5]. In small and large ruminants and pseudo ruminants such as dromedary camels (Camelus dromedaries), S. aureus causes mastitis and therefore negatively impacts the productivity of the dairy industry worldwide [6, 7].
Zoonotic transmission of S. aureus has also been reported [8, 9]. In arid and semi arid regions of the Greater Horn of Africa, camels represent an important and valuable livestock species that provides a significant percentage of the population with animal protein, particularly from milk . Moreover, camel milk is often consumed raw without proper heat-treatment, which increases the risk of acquiring infections with zoonotic pathogens [11, 12].
Currently our knowledge of bacterial pathogens in camels is rather limited . S. aureus has been reported to cause infections of the skin, udder, eyes and joints [14–17] in camels. In North Kenya between 1999 and 2004, the prevalence of S. aureus in camels has been reported as 54 % in closed skin abscesses, 36 % in open skin abscesses, 39 % in skin necrosis and 31 % in lymph node abscesses . A recent survey reports the prevalence of intramammary infections (IMI) associated with S. aureus as 11 % in lactating camels in Kenya . A study has also reported genotype data and identified ‘candidate’ virulence factors of S. aureus strains in Middle Eastern camels . Here we present the complete genome sequence, annotation and comparative analysis of the S. aureus ST30 strain ILRI_Eymole1/1 isolated from a nasal swab of a dromedary camel in Kenya.
Classification and features
The S. aureus strain ILRI_Eymole1/1 was isolated in Kenya in 2004 from a nasal swab of a camel. It was identified as a member of the Staphylococcus aureus species on the basis of standard microbiological procedures  combined with a species-specific PCR . S. aureus is a Gram-positive, coccus shaped, non-motile, nonspore forming and facultative anaerobic bacterium. S. aureus were grown on agar. Agar pieces were cut out and fixed in 150 mM HEPES, pH 7.35, containing 1.5 % formaldehyde and 1.5 % glutaraldehyde for 30 min at room temperature and at 4 °C over night. After dehydration in acetone and critical point drying, cells were gold sputtered and observed in a Philips SEM 505. Images were acquired using 10 kV at 10.000×/20 nm spot size or 40.00×/10 nm spot size. The bacterial cells are 0.5 to 1.0 mm in diameter, and occurs either singly or in the form of pairs or clusters (Fig. 1). The culture produces smooth, circular, glistening colonies of diameter > 5 mm. It produces a grey pigment. The general features of S. aureus strain Eymole1/1 are presented in Table 1 and Additional file 1: Table S1. The optimal growth temperature range is 37–42 °C. Tolerance to NaCl was tested in liquid medium, LB with NaCl concentrations between 0 and 4 M NaCl. Cells were grown overnight at 37 °C.
Carbohydrate utilization was tested using ID 32 STAPH, a standardized system for the identification of the genera Staphylococcus , Micrococcus etc. (bioMérieux, Inc, Box 15969,Durham, NC 27704-0969 / USA). These tests showed positive results for glucose, fructose, mannose, maltose, lactose, trehalose, sucrose and turanose.
The sequence type of the S. aureus isolate was determined using a previously described MLST dataset . ILRI_Eymole1/1 belongs to ST30 MLST group. A BLASTn search  of all five copies of 16S rRNA sequence of ILRI_Eymole1/1 using default parameters revealed 99–100 % identity (with 98–100 % coverage) with all available S. aureus genomes in the database. The phylogenetic relationship was established using the 16S rRNA sequences of the type strains defining the genus Staphylococcus (accession numbers are provided in Additional file 2: Table S2). In addition, the 16S rRNA sequences of 9 S. aureus isolates (NC_021554, NC_017333, NC_017349, NC_022113, NC_002952, NC_017342, NC_002758, NC_002745, NC_020529) were extracted from the genome sequences, and a neighbor joining phylogenetic tree was constructed with MEGA v.6.06 (Fig. 2). The tree illustrates the close relationship of S. aureus ILRI_Eymole1/1 with S. aureus isolates from ST 30, 36, 5, 45 and the S. aureus type strain L36472 (Fig. 2). The position relative to other species within the genus Staphylococcus is also illustrated. Bacillus subtilis type strain DSM10 was used as an outgroup for the genus Staphylococcus .
Genome sequencing information
Genome project history
Twenty three strains of S. aureus have been isolated from healthy and diseased camels in East Africa using standard methods (1). The strains were isolated using primary cultivation on Columbia blood agar plates (Oxoid, UK) and were sub-cultured on mannitol-salt agar plates (Oxoid, UK). Afterwards the strains were subjected to multi locus sequence typing (2). Four strains belonged to sequence type 30, previously characterized in humans. The other isolates had novel sequence types that were likely to be camel specific. We selected the strain ILRI-Eymole1/1 for subsequent analysis since we wanted to elucidate the relationship of S. aureus isolates from camels and humans in order to project a zoonotic potential. S. aureus ILRI_Eymole1/1 was isolated from a nasal swab taken (transport Amies swab w/o charcoal, Copan, Italy) from a dromedary camel calf with rhinitis in Kenya in 2004.
Growth conditions and genomic DNA preparation
The strain was grown in 10 ml liquid Brain heart medium (Carl Roth, Germany) at 37 °C and 200 rpm overnight. The strain was grown in 10 ml liquid Brain heart medium (Carl Roth, Germany) at 37 °C and 200 rpm overnight. The bacterial cells were pelleted using centrifugation at 5000 × g for 20 min. The supernatant was discarded and cells were subjected to genomic DNA isolation using the PureLink™ Genomic DNA Mini Kit (Invitrogen, USA) according to vendor’s instructions. The DNA was quantified using the Qubit® 3.0 Fluorometer (Thermo Scientific, Kenya) and the Qubit® dsDNA BR Assay Kit (Thermo Scientific, Kenya). The DNA concentration was 84.6 ng/ul, the 260/280 and 260/230 ratios were 1.49 and 0.56, respectively. To remove impurities, the DNA was further cleaned using a ratio of 1.6 AMPure beads (ref).
Genome sequencing and assembly
Genome sequencing of S. aureus ILRI_Eymole1/1 was performed using the Illumina Genome Analyzer GAIIx platform. A 300 bp paired-end library with an average insert size of 550 bp was sequenced. The software MIRA v 4.0  was used to assemble the S. aureus ILRI_Eymole1/1 genome, using as the input 1,154,246 Illumina paired-end reads. The de novo genome assembly generated a total of 118 contigs with average coverage of 109 × and average quality of 83 (Table 2). The whole genome alignment tools Mauve  and MUMmer v 3.2.2  were used to order contigs of length greater than 1000 bp (69 contigs) against a reference genome sequence MRSA252/NC_002952 . A complete genome sequence was obtained by joining the ordered contigs on the basis of their overlaps. The assembly output ACE file was viewed and analyzed in Tablet viewer version 1.13.05.17 .
The complete genome sequence of S. aureus ILRI_Eymole1/1 was annotated using RAST . Ribosomal RNA genes were identified using RNAmmer server v 1.2 , and the tRNA genes were predicted using tRNAscan-SE v 1.21 . The COG genes and associated functional categories information were downloaded from the COG database . The COG categories were assigned to the ILRI_Eymole1/1 genome annotation using blastp v 2.2.28  against the COG genes collection/ myva. Genes with more than 40 % amino acid identity and with e-values of less than 0.00005 were classified as putative homologues of genes within the COG database, and functional categories were assigned. Signal peptides were predicted using the SignalP v 4.1 server , and the transmembrane helices/ membrane spanning domains were identified using TopPred2 . Phage-like sequences were predicted using PHAST .
The S. aureus ILRI_Eymole1/1 genome is a circular chromosome of 2,874,302 bp with a GC-content of 32.88 %. A total of 2831 genes were predicted comprising 2755 protein encoding genes, 60 tRNA genes and 16 rRNA genes (Table 3, Fig. 3). Five copies of both 16S and 23S rRNA genes and six copies of 5S rRNA genes were identified. Among the predicted protein encoding genes, 652 (23.66 %) were hypothetical proteins. A total of 162 genes (5.88 %) were predicted to encode proteins with secretory signal peptides (potentially targeted to the secretory pathway) and 1040 (37.75 %) were genes encoding proteins with transmembrane helices or membrane spanning proteins. A total of 2054 (74.56 %) predicted genes were assigned to COG functional categories, while 701 (25.44 %) were not present within the COG collection (Table 4).
Insights from the genome sequence
We performed a comparative analysis of the camel S. aureus ILRI Eymole1/1 isolate of sequence type 30 with 16 previously sequenced ST30 S. aureus isolates, two ST36 methicillin resistance Staphylococcus aureus isolates MRSA252, EMRSA16 and one ST431 S. aureus isolate M809, which together belong to the clonal complex 30 (CC30). ST36 and ST431 are single locus MLST variants of ST30. Previously sequenced S. aureus complete genome sequences were downloaded from the NCBI FTP site  (accession numbers are provided in Additional file 3: Table S3), and CC30 isolates were selected by analyzing their house keeping genes using the S. aureus MLST database . The collection of draft S. aureus CC30 genomes was derived from previous studies [36, 37].
Core genome analysis and COG classification
All 20 S. aureus CC30 genomes were annotated using the RAST server, and amino acid sequences of protein encoding genes from all CC30 genomes were used for the core genome analysis. Blastp searching of protein sequences of all CC30 isolates was carried out using local blastp v 2.2.28 . The genes matching in all CC30 genomes with > =80 % identity, e-value < 0.00005, and alignment length > = 50 % were classified as core genes using custom scripts (Additional file 4: Supplementary material S4). The core genes were further analyzed for their COG functional classification using a matching criterion of > = 40 % identity and an e-value < 0.00005. Among 2163 core genes, 1810 (83.68 %) were present in COG database, whereas 353 (16.32 %) were not present in COG database. The functional classification of these genes is shown in Fig. 4.
S. aureus ILRI_Eymole1/1 variable genes (shared with some of the CC30 S. aureus genomes), and isolate-specific genes were also identified. We identified 2163 core genes (78.51 % of the total protein encoding genes), 507 (18.40 %) variable protein-encoding genes and 85 (3.09 %) isolate-specific genes.
The colonization and adhesion of S. aureus to the nasal epithelial cells is thought to be mediated by surface proteins ClfB, IsdA and the serine-aspartic acid repeat proteins SdrC and SdrE. A published study demonstrated that a mutant lacking these four proteins did not exhibit the adherence phenotype . S. aureus ILRI_Eymole1/1 possesses genes encoding fibrinogen-binding protein ClfB (CEH27447), adhesin proteins SdrC (CEH25318) and SdrE (CEH25319), in common with a subset of the CC30 S. aureus genomes. A gene encoding Heme/ Iron regulated surface protein IsdA (CEH26009) was present among the core protein repertoire of ILRI_Eymole1/1 isolate, and is known to be important for S. aureus infection of human skin, through mediating resistance to skin innate defense mechanisms .
S. aureus ILRI_Eymole1/1 possessed genes encoding many fibrinogen-binding proteins, including clumping factor/ fibrinogen binding protein ClfA, (CEH26520: variable gene), fibronectin/ fibrinogen binding protein FnBP (CEH25930: core gene), extracellular fibrinogen-binding protein Efb (CEH25981: core gene), fibronectin binding protein FnbB (CEH27312: variable gene), fibronectin binding protein FnbA (CEH27314: variable gene), and clumping factor ClfB, fibrinogen binding protein (CEH27447: variable gene). These FnBPs bind the host fibronectin receptor β1-integrins to promote S. aureus invasion of various mammalian cells including epithelial cells, endothelial cells and fibroblasts. These cells do not require specific co-receptors for S. aureus .
A gene encoding an extracellular adherence protein of broad specificity Eap/Map (CEH26760: core gene) was also identified in the camel S. aureus isolate. This protein has been reported to be involved in S. aureus internalization into the host cells. Eap is known to be responsible for agglutination of bacterial cells by rebinding to the surface of S. aureus . It shows dual affinity for the cell surface plasma proteins as well as the bacterial surface. Eap plays a complementary role together with FnBP, in the internalization and long time persistence of S. aureus within eukaryotic cells. It was found to be a key component of the novel internalization pathway that works either in parallel with, or in addition to, the FnBP dependent internalization pathway .
Sec-independent Ess secretion pathway/ Type VII secretion system, T7SS
Many Gram-positive bacterial species, including S. aureus , secrete exotoxins or virulence factors across the membrane, through signal peptides or the Sec translocon. A Sec-independent translocation of these factors has also been reported in Gram-positive bacteria. Human S. aureus has been shown to secret the ESAT-6-like secretory proteins EsxA and EsxB. The genes encoding these proteins cluster in the genome as an operon together with several additional genes to form a secretion system, known as the type 7 secretion system (T7SS) that is involved in bacterial pathogenicity [42, 43]. The S. aureus isolate ILRI_Eymole1/1 possessed genes encoding proteins related to T7SS; which were also present in all other CC30 S. aureus genomes. These encoded the secretory antigen precursor SsaA (CEH25002), ESAT-6/Esx family secreted protein EsxA/YukE (CEH25003), putative secretion accessory protein EsaA/YueB (CEH25004), putative secretion system component EssA (CEH25005), putative secretion accessory protein EsaB/YukD (CEH25006), putative secretion system component EssB/YukC (CEH25007), and a FtsK/SpoIIIE family protein, together with putative secretion system component EssC/YukA (CEH25008).
Isolate specific genes
Out of the total 85 isolate specific genes encoded by S. aureus ILRI_Eymole1/1, 79 genes (92.94 %) were clustered into six large insertions. Four insertions were putative bacteriophages comprising four complete phages with sizes of 52.5 kb, 30 kb, 60.3 kb and 58.8 kb, respectively. All phage sequences possessed attL and attR integration sequences at the forward and reverse ends. Superantigen pathogenicity islands (SaPI) are mobile genetic elements in Gram-positive bacteria including S. aureus that carry genes associated with superantigens, virulence, resistance and metabolic functions; also named as S. aureus pathogenicity islands ‘SaPI’. These are known for their strong association with temperate phages and result in high transfer frequencies . Two insertions constituted complete SaPI islands (SaPIcam1 and SaPIcam2) at positions 426,323-443,273 and 758,187-774,130. These were confirmed by the identification of forward and reverse sequences at the 5′–3′ ends of previously characterized SaPIs, namely SaPIbov and SaPImw2 . SaPIcam1 and SaPIcam2 both had integrase and terminase encoding genes at their termini. SaPIcam1 also possessed an HTH-type transcriptional regulator LrpC adjacent to 3′ end. The SaPIcam2 contained a candidate superantigen tst gene (toxic shock syndrome toxin 1 TSST-1, as part of variable gene content) located adjacent to the 3′ end. The ‘SaPI2’ island in the CC30 isolates encodes a TSST-1 gene, and these have a clonal association with CC30 nasal infective and bacteremia causing isolates [36, 46]. Among S. aureus CC30 isolates analyzed in this study, only ILRI_Eymole1/1, EMRSA16, A017934_97, Btn1260 and MN8 genomes contained ‘SaPI2’, encoding the tst gene. All other CC30 genomes possessed a ‘SaPI1’ island, and therefore encoded Ear, a secretory protein, at the 3′ end. The presence of the high level of isolate specific genes (92.94 %) in these phage insertions and the SaPI islands strongly suggests the acquisition of these genes through lateral gene transfer from either phages or heterologous bacterial species harboring these insertions.
Phylogeny using polymorphic set of core genes
We determined the phylogenetic relationship among the isolates using a stringently defined set of 283 core genes that were shared among the CC30 isolates. Two S. aureus genomes from ST1 and ST5 were also included in this analysis as outgroups. The core genes were defined among these 22 S. aureus genomes using the criteria of (identity > = 95 % and < 100 %), (e-value < 0.00005) and (alignment length > = 90 %). Duplicate copies of genes were filtered out, resulting in the final total of 283 core genes. Multiple sequence alignment of the concatenated sequences of these genes was performed using the Mugsy aligner , generating an alignment comprising 316,359 nucleotides from each isolate. We estimated a maximum likelihood phylogeny using PhyML v. 3.0 . The General Time-Reversible (GTR) model was used, where the base frequencies and the relative substitution rates between them were estimated by maximizing the likelihood of the phylogeny. For estimating the tree topology both nearest neighbor interchange and subtree pruning and regrafting methods were used. One hundred bootstrap replicates were run (Fig. 5a and b).
In both rooted and unrooted trees (Fig. 5a and b), the human CC30 isolates group in three clusters, in agreement with a published study . The camel S. aureus isolate ILRI_Eymole1/1 clusters in the CC30 (Fig. 5a), but is genetically distant from human CC30 S. aureus isolates (Fig. 5b).
Here we report the first genome of a S. aureus isolated from Camelus dromedarius . Our analysis shows that a high proportion of isolate-specific genes were located in putative phage insertions and SaPI islands in the camel isolate clearly distinguished it from human isolates. The analysis based on a polymorphic set of core genes clearly shows that the camel S .aureus isolate belongs to ST30 but this isolate has greater genetic difference when compared to human isolates. Therefore, we consider the likelihood of exchange between camel and human populations low. However, this is the complete genome of a single S. aureus from a camel. The analysis of additional S. aureus isolates from camels and humans living in the same area, followed by a detailed comparative and phylogenetic analysis will underpin improved understanding of host adaptation and zoonotic potential.
Clonal Complex 30
Methicillin Resistance Staphylococcus aureus
Methicillin Susceptible/Sensitive Staphylococcus aureus
Sung JM-L, Lloyd DH, Lindsay JA. Staphylococcus aureus host specificity: comparative genomics of human versus animal isolates by multi-strain microarray. Microbiology. 2008;154:1949–59.
Smyth DS, Feil EJ, Meaney WJ, Hartigan PJ, Tollersrud T, Fitzgerald JR, et al. Molecular genetic typing reveals further insights into the diversity of animal-associated Staphylococcus aureus. J Med Microbiol. 2009;58:1343–53.
Wertheim HFL, Melles DC, Vos MC, Van Leeuwen W, Van Belkum A, Verbrugh HA, et al. Subscription Information: Review The role of nasal carriage in Staphylococcus aureus infections. Lancet Infect Dis. 2005;5:751–62.
Maina EK, Kiiyukia C, Wamae CN, Waiyaki PG, Kariuki S. Characterization of methicillin-resistant Staphylococcus aureus from skin and soft tissue infections in patients in Nairobi, Kenya. Int J Infect Dis. 2013;17:e115–9.
Ladhani S. Bacteraemia due to Staphylococcus aureus. Arch Dis Child. 2004;89:568–71.
Fitzgerald JR. Livestock-associated Staphylococcus aureus: origin, evolution and public health threat. Trends Microbiol. 2012;20:192–8.
Guinane CM, Ben Zakour NL, Tormo-Mas MA, Weinert LA, Lowder BV, Cartwright RA, et al. Evolutionary genomics of Staphylococcus aureus reveals insights into the origin and molecular basis of ruminant host adaptation. Genome Biol Evol. 2010;2:454–66.
Christou L. The global burden of bacterial and viral zoonotic infections. Clin Microbiol Infect. 2011;17:326–30.
Petersen A, Stegger M, Heltberg O, Christensen J, Zeuthen A, Knudsen L, et al. Epidemiology of methicillin-resistant Staphylococcus aureus carrying the novel mecC gene in Denmark corroborates a zoonotic reservoir with transmission to humans. Clin Microbiol Infect. 2013;19:E16–22.
Abdurahman OAS. Udder health and milk quality among camels in the Errer valley of eastern Ethiopia. In: Livestock Research for Rural Development, vol. 18. 2006. p. 8.
Gautret P, Benkouiten S, Gaillard C, Parola P, Brouqui P. Camel milk-associated infection risk perception and knowledge in French Hajj pilgrims. Vector Borne Zoonotic Dis. 2013;13:425–7.
Sprague LD, Al-Dahouk S, Neubauer H. A review on camel brucellosis: a zoonosis sustained by ignorance and indifference. Pathog Glob Health. 2012;106:144–9.
Fischer A, Liljander A, Kaspar H, Muriuki C, Fuxelius H-H, Bongcam-Rudloff E, et al. Camel Streptococcus agalactiae populations are associated with specific disease complexes and acquired the tetracycline resistance gene tetM via a Tn916-like element. Vet Res. 2013;44:86.
Monecke S, Ehricht R, Slickers P, Wernery R, Johnson B, Jose S, et al. Microarray-based genotyping of Staphylococcus aureus isolates from camels. Vet Microbiol. 2011;150:309–14.
Younan M, Bornstein S. Papers & Articles Lancefield group B and C streptococci in East African camels (Camelus dromedarius). Vet Rec. 2007;160:330–5.
Younan M, Ali Z, Bornstein S, Müller W. Application of the California mastitis test in intramammary Streptococcus agalactiae and Staphylococcus aureus infections of camels (Camelus dromedarius) in Kenya. Prev Vet Med. 2001;51:307–16.
Wernery U, Kinne J, Schuster RK. Camelid Infectious Disorders. Paris, France: World Organisation for Animal Health OIE; 2014. p. 500.
Carter GRCJ. Diagnostic procedures in veterinary bacteriology and mycology. 5th ed. San Diego, Calif: Academic; 1990. p. 620.
Hussain M, Von Eiff C, Sinha B, Joost I, Herrmann M, Peters G, et al. eap Gene as Novel Target for Specific Identification of Staphylococcus aureus. J Clin Microbiol. 2008;46:470–6.
Enright MC, Day NPJ, Davies CE, Peacock SJ. Multilocus sequence typing for characterization of methicillin- resistant and methicillin-susceptible clones of Staphylococcus aureus. J Clin Microbiol. 2000;38:1008–15.
Altschul SF, Gish W, Miller W, Myers EW, Lipmanl DJ. Basic Local Alignment Search Tool. J Mol Biol. 1990;215:403–10.
Chevreux B, Wetter T, Suhai S. Genome sequence assembly using trace signals and additional sequence information. In: Computer science and biology: proceedings of the German Conference on Bioinformatics, vol. 99. Göttingen, Germany: GCB; 1999. p. 45–56.
Darling AE, Mau B, Perna NT. Progressive Mauve: multiple genome alignment with gene gain, loss and rearrangement. PLoS One. 2010;5:e11147.
Kurtz S, Phillippy A, Delcher AL, Smoot M, Shumway M, Antonescu C, et al. Versatile and open software for comparing large genomes. Genome Biol. 2004;5:R12.
Holden MTG, Feil EJ, Lindsay JA, Peacock SJ, Day NPJ, Enright MC, et al. Complete genomes of two clinical Staphylococcus aureus strains : Evidence for the rapid evolution of virulence and drug resistance. Proc Natl Acad Sci U S A. 2004;101:9786–91.
Milne I, Stephen G, Bayer M, Cock PJA, Pritchard L, Cardle L, et al. Using Tablet for visual exploration of second-generation sequencing data. Brief Bioinform. 2013;14:193–202.
Aziz RK, Bartels D, Best AA, DeJongh M, Disz T, Edwards RA, et al. The RAST Server: rapid annotations using subsystems technology. BMC Genomics. 2008;9:75–89.
Lagesen K, Hallin P, Rødland EA, Staerfeldt H-H, Rognes T, Ussery DW. RNAmmer: consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007;35:3100–8.
Lowe TM, Eddy SR. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucl Acids Res. 1997;25:955–64.
Tatusov RL, Galperin MY, Natale DA, Koonin EV. The COG database: a tool for genome-scale analysis of protein functions and evolution. Nucleic Acids Res. 2000;28:33–6.
BLAST+ executables. ftp://ftp.ncbi.nlm.nih.gov/blast/executables/blast+/.
Petersen TN, Brunak S, von Heijne G, Nielsen H. SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat Methods. 2011;8:785–6.
Von Heijne G. Membrane hydrophobicity protein structure prediction analysis and the positive-inside. J Mol Biol. 1992;225:487–94.
Zhou Y, Liang Y, Lynch KH, Dennis JJ, Wishart DS. PHAST: A Fast Phage Search Tool. Nucl Acids Res. 2011;39:W347–52.
NCBI FTP site. https://www.ncbi.nlm.nih.gov/Ftp/.
McGavin MJ, Arsic B, Nickerson NN. Evolutionary blueprint for host- and niche-adaptation in Staphylococcus aureus clonal complex CC30. Front Cell Infect Microbiol. 2012;2:1–13.
Harro JM, Daugherty S, Bruno VM, Jabra-rizk MA, Rasko DA, Shirtliff E. Draft Genome Sequence of the Methicillin-Resistant Staphylococcus aureus Isolate MRSA-M2. Genome Announc. 2013;1:e00037–12.
Corrigan RM, Miajlovic H, Foster TJ. Surface proteins that promote adherence of Staphylococcus aureus to human desquamated nasal epithelial cells. BMC Microbiol. 2009;9:22.
Clarke SR, Mohamed R, Bian L, Routh AF, Kokai-Kun JF, Mond JJ, et al. The Staphylococcus aureus surface protein IsdA mediates resistance to innate defenses of human skin. Cell Host Microbe. 2007;1:199–212.
Sinha B, Francois P, Que Y-A, Hussain M, Heilmann C, Moreillon P, et al. Heterologously expressed Staphylococcus aureus fibronectin-binding proteins are sufficient for invasion of host cells. Infect Immun. 2000;68:6871–8.
Haggar A, Hussain M, Lo H, Herrmann M, Norrby-teglund A, Flock J. Extracellular Adherence Protein from Staphylococcus aureus Enhances Internalization into Eukaryotic Cells. Infect Immun. 2003;71:2310–7.
Burts ML, Williams WA, DeBord K, Missiakas DM. EsxA and EsxB are secreted by an ESAT-6-like system that is required for the pathogenesis of Staphylococcus aureus infections. Proc Natl Acad Sci U S A. 2005;102:1169–74.
Anderson M, Chen Y-H, Butler EK, Missiakas DM. EsaD, a secretion factor for the Ess pathway in Staphylococcus aureus. J Bacteriol. 2011;193:1583–9.
Ubeda C, Maiques E, Barry P, Matthews A, Tormo MA, Lasa I, et al. SaPI mutations affecting replication and transfer and enabling autonomous replication in the absence of helper phage. Mol Microbiol. 2008;67:493–503.
Sato’o Y, Omoe K, Ono HK, Nakane A, Hu D-L. A novel comprehensive analysis method for Staphylococcus aureus pathogenicity islands. Microbiol Immunol. 2013;57:91–9.
Holtfreter S, Grumann D, Schmudde M, Nguyen HTT, Eichler P, Strommenger B, et al. Clonal distribution of superantigen genes in clinical Staphylococcus aureus isolates. J Clin Microbiol. 2007;45:2669–80.
Angiuoli SV, Salzberg SL. Mugsy: fast multiple alignment of closely related whole genomes. Bioinformatics. 2011;27:334–42.
Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.
Woese CR, Kandler O, Wheelis ML. Towards a natural system of organisms: Proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci U S A. 1990;87:4576–9.
Gibbons NE, Murray RGE. Proposals Concerning the Higher Taxa of Bacteria. Int J Syst Bacteriol 1978;1:1–6
Schleifer K-H. Phylum XIII. Firmicutes. In: De Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, Schleifer KH WW, editors. Bergey’s Manual of Systematic Bacteriology, Second Edition, Volume 3. New York: Springer; 2009. p. 19.
Ludwig W, Schleifer K-H, Whitman WB. Class I. Bacilli class nov. In: De Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, Schleifer KH WW, editors. Bergey’s Manual of Systematic Bacteriology, Second Edition, Volume 3. Springer-Verlag, New York; 2009. p. 19–20
Anon. List of new names and new combinations previously effectively, but not validly, published. List no. 132. Int J Syst Evol Microbiol. 2010;60:469–72.
De Vos P. Order I. Bacillales Prévot 1953, 60AL. In: De Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, Schleifer KH WW, editors. Bergey’s Manual of Systematic Bacteriology, Second Edition, Volume 3. Springer-Verlag, New York; 2009. p. 20.
Yakoubou S, Xu D, Côté J. Phylogeny of the Order Bacillales inferred from 3′ 16S rDNA and 5′ 16S-23S ITS nucleotide sequences. Nat Sci. 2010;2:990–7.
Skerman VBD, McGowan V, Sneath PHA. Approved lists of bacterial names. Int J Syst Bacteriol. 1980;30:225–420.
Prévot AR. Dictionnaire des Bactéries Pathogènes. 2nd ed. In: Hauderoy P, Ehringer G, Guillot G, J. M, Prévot AR, Rosset D, Urbain A, editors. Paris: Masson et Cie; 1953. p. 1–692.
Schleifer K-H, Bell JA. Family VIII. Staphylococcaceae fam. nov. In: De Vos P, Garrity G, Jones D, Krieg NR, Ludwig W, Rainey FA, Schleifer KH WW, editors. Bergey’s Manual of Systematic Bacteriology, Second Edition, Volume 3. Springer-Verlag, New York; 2009. p. 392.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene Ontology: tool for the unification of biologybiology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Carver T, Thomson N, Bleasby A, Berriman M, Parkhill J. DNAPlotter: circular and linear interactive genome visualization. Bioinformatics. 2009;25:119–20.
Gene Ontology Consortium. Gene ontology consortium: going forward. Nucleic acids research 43, no. D1 (2015): D1049-D1056. http://www.geneontology.org
Tamura, Koichiro, Glen Stecher, Daniel Peterson, Alan Filipski, Sudhir Kumar. "MEGA6: molecular evolutionary genetics analysis version 6.0."Molecular biology and evolution 30, no. 12 (2013): 2725-729.
Carver, Tim, Nick Thomson, Alan Bleasby, Matthew Berriman, and Julian Parkhill. "DNAPlotter: circular and linear interactive genome visualization."Bioinformatics 25, no. 1 (2009): 119-20.
This work was supported by the Consultative Group of International Agricultural Research (CGIAR) program on Agriculture for Nutrition and Health. Erik Bongcam-Rudloff was partially supported by the Project ALLBIO grant number EU FP7; 289452; KBBE.2011.3.6-02. Higher Education Commission (HEC) of Pakistan supported Saima Zubair. The Centrum of International Migration (CIM) supported Anne Fischer and Anne Liljander. Support by BILS (Bioinformatics Infrastructure for Life Sciences) and UPPMAX/UPPNEX are gratefully acknowledged. We thank Gerhard Preiss for excellent maintenance and help with electron microscopes and Andrea Kofink-Germershausen and Sabine Fiedler for excellent technical assistance.
The authors declare that they have no competing interests.
IR and MY collected the data; AL, JM, JH performed the laboratory work; SZ, AF, HG, MM, MIM performed the analysis; EB-R and JJ designed the study; and SZ, AF, RBP, JJ wrote the manuscript. All authors read and approved the final manuscript.
Associated MIGS record. (DOC 73 kb)
The 16S rRNA sequences of the type strains of genus Staphylococcus used in phylogenetic tree (Fig. 1). (DOC 47 kb)
General genomic features of twenty CC30 S. aureus genomes. (DOC 53 kb)
Supplementary data. (DOCX 143 kb)