Next Article in Journal
Data Augmentation Enhances Plant-Genomic-Enabled Predictions
Previous Article in Journal
Transcriptomic Profiling at the Maternal-to-Zygotic Transition in Leech, Helobdella austinensis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 15
Issue 3
10.3390/genes15030284
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Whole Genome Sequencing and Comparative Genomics of Six Staphylococcus schleiferi and Staphylococcus coagulans Isolates

by
Mohamed A. Abouelkhair
and
Stephen A. Kania
*
Department of Biomedical and Diagnostic Sciences, University of Tennessee College of Veterinary Medicine, Knoxville, TN 37996, USA
*
Author to whom correspondence should be addressed.
Genes 2024, 15(3), 284; https://doi.org/10.3390/genes15030284
Submission received: 16 January 2024 / Revised: 6 February 2024 / Accepted: 22 February 2024 / Published: 24 February 2024
(This article belongs to the Section Microbial Genetics and Genomics)
Browse Figures
Versions Notes

Abstract

:
Staphylococcus schleiferi and Staphylococcus coagulans, closely related bacterial species within the Staphylococcus genus, present a challenge in classification and diagnosis due to their close genetic proximity and overlapping phenotypic features. Moreover, our understanding of the virulence mechanisms in staphylococcal species, beyond the extensively studied Staphylococcus aureus, remains limited, underscoring the importance of using comparative data to enhance our insights into virulence within these bacterial species. This study employed a comprehensive approach, utilizing comparative genomics, to identify genomic distinctions between S. schleiferi and S. coagulans, aiming to address the challenges in the accurate classification and diagnosis of these organisms and identify unique features. Whole genome sequencing was performed on six clinical isolates, and their genomes were compared to identify variations in gene content and virulence factors. De novo assembly and annotation revealed two samples as S. coagulans and four samples as S. schleiferi. Analysis of the core genomes revealed conserved regions crucial for defining species identity, while accessory genomic elements contained unique genes, possibly impacting the pathogenicity of the species.

1. Introduction

Staphylococcus schleiferi is a Gram-positive bacterium within the genus Staphylococcus. It has the ability to cause infections in both veterinary and human medicine. S. coagulans has been identified as a commensal and opportunistic pathogen affecting primarily dogs, where it is frequently detected on the skin [1] and external ear canal [2,3]. Furthermore, it has been linked to external ear otitis [1,2] and pyoderma [2,3]. Although rare, there have been documented cases of S. coagulans producing opportunistic infections in immunocompromised humans [3,4]. The genetic similarity and common phenotypic traits of these two bacterial species have provided considerable challenges to correct classification and diagnosis [5,6]. Their clinical importance justifies the need for a better understanding of their genetic and pathogenic properties as well as the development of new methods of diagnosis, prophylaxis and treatment.
The former subspecies of S. schleiferi, S. schleiferi subsp. coagulans, and S. schleiferi subsp. schleiferi, have been reclassified as different species—S. schleiferi and S. coagulans, respectively [7]. This reclassification was based on genome-to-genome distance, average nucleotide identities, and partial 16S rRNA genes comparing S. schleiferi subsp. coagulans GA211T with S. schleiferi subsp. schleiferi ATCC 43808T. A more extensive study highlighted the importance of genetic-based diagnostics in light of the failure of current, commonly used laboratory methods, including matrix-assisted laser desorption ionization–time of flight (MALDI-TOF), to distinguish them [6]. Even the coagulase test, which has historically been used to distinguish the (former) subspecies, has been shown to be unreliable for this purpose because of false reactions caused by pseudo-coagulase activity. Moreover, 16S rRNA sequencing is often used to distinguish bacterial species; however, it has been shown that these sequences in public databases contain a large number of nucleotide errors and misclassifications [8,9], and even 16S rRNA genes, though thought to be highly conserved, can have variability within a single species [10]. To address this problem, we employ comparative genomics techniques that can be used to elucidate the complexities of gene content, pangenomes, and dynamic genetic components, identifying distinguishing traits that can serve as robust classifiers. The identification of these genetic markers provides the basis for the development of cost-effective methods to both distinguish species and also to overcome the difficulties in accurate classifications of clinical isolates.
In addition, this study helps fill a gap in our knowledge of staphylococcal pathogenicity pathways. Whereas numerous studies have been conducted on S. aureus [11,12,13,14,15,16,17] and S. pseudintermedius [18,19,20,21], the virulence mechanisms of closely related species such as S. schleiferi and S. coagulans remain largely unexplored. This study adds to our understanding of the pathogenic potential and mechanisms of the Staphylococcus genus through genomic analysis. Methicillin-resistance in S. coagulans and S. schleiferi, as with other staphylococci, is associated with an increasing trend toward multi-drug resistance and poses a significant treatment challenge [22,23,24]. Thus, it is essential that alternative treatments and prophylaxis be developed. Strategies for this process require the identification of virulence proteins, resistance mechanisms, and knowledge of strain conservation.
This study explores genomic variability in S. schleiferi and S. coagulans isolates. It addresses challenges in classification and diagnosis, examines virulence mechanisms within two important members of the Staphylococcus genus, and lays the groundwork for improved microbial identification and infectious disease management. These areas of research not only contribute to scientific understanding but also hold promise for advancements in clinical microbiology and patient care.

2. Materials and Methods

2.1. Bacterial Strains, Media, and Growth Conditions

Bacteria propagated in this study included three S. schleiferi isolates (S. schleiferi 192, S. schleiferi 205, and S. schleiferi 214), two S. coagulans (S. coagulans 191 and S. coagulans 196) submitted to diagnostic laboratories from human cases in the USA, and the S. schleiferi type strain ATCC43808 isolated from a human patient in France and obtained from a culture collection. The bacteria were streaked on blood agar plates, and a single colony from each isolate was inoculated into 5 mL of sterile trypticase soy broth (TSB) (BD Biosciences, San Jose, CA, USA; Cat. no. RS1-011-21) and grown overnight with shaking at 225 rpm in a 37 °C incubator.

2.2. DNA Extraction, Library Preparation, and Whole Genome Sequencing

DNA was extracted using the MasterPure DNA purification kit (Epicentre, Madison, WI, USA, cat. no. MCD85201) following the manufacturer’s protocol. The Nextera DNA sample prep kit (Illumina, Inc., San Diego, CA, USA) was used to prepare libraries for sequencing. Genomes were sequenced using the MiSeq platform (Illumina, Inc.) with two runs (2 × 75 bps) at the University of Tennessee Genomics Core facility. Sequences were trimmed using BBDuk and de novo assembled using Geneious Prime® 2023 [25]. The quality of the assembled genomes was determined using the quality assessment tool for genomic assemblies (QUAST). [26]. The sequences were annotated with the NCBI Prokaryotic Genome Annotation Pipeline version 4.6 (https://www.ncbi.nlm.nih.gov/genome/annotation_prok, accessed on 17 October 2023) using the best-placed reference protein set with GenMarkS+.

2.3. Comparative Genomics Analysis

For whole-genome-based taxonomic analysis, the genome sequence data were uploaded to the Type (Strain) Genome Server (TYGS). This bioinformatics platform is available at https://tygs.dsmz.de, accessed on 17 October 2023 [27]. The analysis included methodological features for genome-based classification and nomenclature of prokaryotes described by Meier-Kolthoff et al. [28]. Information on synonymy, nomenclature, and associated taxonomic data was obtained using the LPSN (the List of Prokaryotic names with Standing in Nomenclature) which is available at https://lpsn.dsmz.de, accessed on 17 October 2023 [28]. The results were obtained from TYGS on 18 December 2023. The closest type strain genomes were determined using a sequential two-step process. First, using the MASH algorithm, a comparison was made with our genomes and all type strain genomes contained in the TYGS database [29]. Second, 10 closely related type strains were identified based on their 16S rRNA gene sequences. RNAmmer was used to extract them from the uploaded genomes [30], and each sequence was then subjected to BLAST analysis [31] against the 16S rRNA gene sequence of each of the type strains contained within the TYGS database. This approach, based on their bitscore, was used to find the 50 most similar type strains for each query genome. Subsequently, the Genome BLAST Distance Phylogeny (GBDP) approach with the algorithm coverage and distance formula d5 was used to determine precise distances [32]. The 10 type strain genomes closest to each of our genomes was thus determined. To detect phylogenomic inference, pairwise comparisons of genomes were made using GBDP, and intergenomic distances were inferred using the trimming algorithm and distance formula d5 [32]. For each analysis, a total of 100 replicates were determined, and the protein sequences of entire proteomes were used to conduct an additional GBDP phylogenomic analysis, allowing a robust resolution of the phylogeny of remotely related strains. Digital confidence intervals and DNA: DNA hybridization (DDH) values were determined using GGDC 4.0 with the recommended settings [28,32]. The intergenomic distances obtained from this analysis were used to produce a balanced minimum evolution tree with branch support inferred from 100 pseudo-bootstrap replicates using FASTME 2.1.6.1 with SPR postprocessing [10]. The trees were rooted at the midpoint [33] and visualized with PhyD3 [34]. The species clustering was performed using a 70% digital DNA: DNA hybridization (dDDH) radius around the 20 type strains as described previously [27]. Subspecies clustering was determined using a 79% dDDH threshold as described by Meier-Kolthoff et al. [35].
The genomes assembled from our isolates were compared against a set of complete genomes of S. schleiferi and S. coagulans downloaded from the NCBI database (Table 1) using Anvi’o 8 [36].
Average nucleotide identity (ANI) values were calculated using Anvi’o 8 [36] via pyANI [37,38], and OrthoANI was used for comparing the genomic similarity between the coding regions of the genomes [39]. For the pangenomic analysis, we downloaded S. schleiferi (n = 50) and S. coagulans (n = 182) assemblies from NCBI (https://www.ncbi.nlm.nih.gov/, accessed on 20 December 2023). The assemblies were subjected to Prokka annotation [40] and the Roary pangenome pipeline [41].
Phylogenetic relationships between S. schleiferi isolates from this study and previously sequenced, assembled, and annotated S. schleiferi isolates were determined using CSI Phylogeny 1.4 [42] (Call SNPs and Infer Phylogeny) using the S. schleiferi type strain MGYG-HGUT-01437 (GCA_902374935.1) as the reference genome using the default settings [minimum depth at single-nucleotide polymorphism (SNP) positions: 10×; minimum relative depth at SNP positions: 10 bp; minimum SNP quality: 10%, minimum distance between SNPs (prune): 30; minimum read mapping quality: 25 and minimum Z-score: 1.96].

3. Results

3.1. Genomic Features of Human Staphylococcus schleiferi and Staphylococcus coagulans

The genome size, GC content, predicted coding sequences, and predicted RNAs of the three S. schleiferi isolates (S. schleiferi 192, S. schleiferi 205, and S. schleiferi 214), two S. coagulans (S. coagulans 191 and S. coagulans 196) submitted to diagnostic laboratories from human cases in the USA, and the S. schleiferi type strain ATCC43808 isolated from a human patient in France are listed in (Table 2).

3.2. Type-Based Species Clustering

The Type (Strain) Genome Server was used for genome-based analysis because of the high quality of its comprehensive, curated database of type strains [28]. This analysis yielded 18 species clusters, and S. coagulans 191 and S. coagulans 196 were classified as S. coagulans, while S. schleiferi 192, S. schleiferi 205, S. schleiferi 214, and S. schleiferi type strain ATCC43808 were classified as S. schleiferi (Figure 1).
Pangenome analysis was generated with Anvi’o-8. The layers represent individual genomes organized by their phylogenomic relationships based on single-copy core genes (Figure 2). In the layers, dark colors indicate the presence of a gene cluster and light colors its absence. The average nucleotide identity (ANI) heatmap illustrates the genetic relatedness of S. schleiferi and S. coagulans, with nucleotide identities ranging between 95% and 100% (Figure 2). Specifically, S. coagulans 191 and S. coagulans 196 form a cluster with other closely related strains, including S. coagulans 2317-03, S. coagulans OT1-1, S. coagulans 5909-02, S. coagulans TSCC54, S. coagulans strain 1031336 reference genome GCA_022557135.1, and S. coagulans 2142-05 (Figure 2). On the other hand, S. schleiferi 192, S. schleiferi 205, S. schleiferi 214, and the S. schleiferi type strain ATCC43808 are grouped together and share a high genetic similarity with S. schleiferi MGYG-HGUT-01437 (Figure 2).
The clustering dendrogram shows the relationship of contigs based on their sequences and distribution across samples. ANI heatmap results for S. schleiferi and S. coagulans vary between 95 and 100 nucleotide identities %. The phylogenomic tree was reconstructed using single-copy genes. In the layers, dark colors indicate the presence of a gene cluster and light colors its absence. R stands for reference genomes that were used in the analysis.
S. schleiferi 192, S. schleiferi 205, S. schleiferi 214, and S. schleiferi ATCC 43808T were compared with the 46 publicly available S. schleiferi genomes using an SNP-based phylogenetic analysis (Figure 3). The percentage of the reference genome contained within all isolates was found to be 93.19%, and 2,310,197 single nucleotide polymorphisms (SNPs) positions were found in the analyzed genomes (Figure 3).
A genome-wide association study (GWAS) conducted using Roary and Scoary (Figure 4 and Figure 5) delineated species-specific genetic markers associated with S. schleiferi and S. coagulans. The comprehensive analysis revealed a distinct repertoire of genes unique to each species, with particular emphasis on the sialidase B (nanB) gene, Zn-dependent alcohol/formaldehyde dehydrogenase (FrmA), metal-dependent amidase/aminoacylase/carboxypeptidase (AbgB), membrane protease (YdiL), DNA glycosylase (YcaQ), and hypothetical protein 24 exclusively linked to S. schleiferi isolates and the chromate transport protein (chrA) gene identified as a unique feature of S. coagulans isolates.

3.3. Nucleotide Sequence Accession Numbers

The whole genome sequences of humans S. schleiferi, S. coagulans 191, S. schleiferi 192, S. coagulans 196, S. schleiferi 205, S. schleiferi 214, and S. schleiferi ATCC 43808T, have been deposited at DDBJ/ENA/GenBank under the accession numbers PNRJ00000000, POVG00000000, POVH00000000, POVI00000000, POVJ00000000 and POVK00000000, respectively.
The S. coagulans 191, S. schleiferi 192, S. coagulans 196, S. schleiferi 205, Staphylococcus schleiferi 214, and S. schleiferi ATCC 43808T accessions described in this paper are PNRJ01000000, POVG01000000, POVH01000000, POVI01000000, POVJ01000000 and POVK01000000, respectively.

4. Discussion

In this study, S. schleiferi and S. coagulans were found to be distinct from each other while also phylogenetically separated from other species on their own branch. This finding, along with difficulties distinguishing them from each other and from other pathogens, in clinical laboratories, justifies our in-depth analysis and comparison of their genomic sequences.
Most staphylococcal studies have focused on the pathogenesis, the development of vaccines, and epidemiological typing tools for S. aureus and S. pseudintermedius. Whereas other less studied staphylococcus species, such as S. schleiferi and S. coagulans, cause fewer infections, they are nonetheless important pathogens. S. coagulans may be underrecognized and underdiagnosed because, as Gram positive cocci that are coagulase positive, they may be misidentified as S. aureus or other staphylococci. Furthermore, S. schleiferi has a high incidence of methicillin resistance, with some reports indicating it may be as high as 57% in canine clinical isolates. Although S. schleiferi and S. coagulans are primarily associated with disease in dogs and occasional zoonotic infections, they have a global distribution and are also present on mink, birds, racoons, foxes, horses, goats, seals, and squirrels.
The goal of this research was to improve our understanding of the S. schleiferi and S. coagulans genomic landscapes, lay the groundwork for the development of an effective S. schleiferi and S. coagulans vaccine, and advance our overall capabilities to develop staphylococcal vaccines in naturally infected hosts. The comparative genomic study identified species-specific indicators, including the chromate transport protein (chrA) gene, which is only present in S. coagulans isolates. This finding is consistent with previous studies and provides a new foundation for discriminating between the two species. Genes found in S. schleiferi but not S. coagulans isolates include the sialidase B (neuraminidase B) (nanB) gene, Zn-dependent alcohol/formaldehyde dehydrogenase (FrmA), metal-dependent amidase/aminoacylase/carboxypeptidase (AbgB), membrane protease (YdiL), DNA glycosylase (YcaQ), and the hypothetical protein 24 gene. In addition to their use in the development of molecular tests, such as PCR, to discriminate between the species, these putative findings may represent some important virulence factors that can be targeted with vaccines or methods of gene regulation to neutralize their activity. Neuraminidase B plays an important role in the evasion of mucus binding in Streptococcus pneumoniae [43]. It controls expression of the neuraminidase A gene (nanA) through the release of sialic acid that upregulates nanA expression. Hammond et al. postulated that NanA and NanB act together to allow S. pneumoniae to evade mucosal defenses.
Formaldehyde is produced by S. aureus via the alternative heme degradation pathway to acquire iron [44]. FrmA is involved in formaldehyde detoxification; however, its role in S. schleiferi infections and iron acquisition has not been studied, although iron uptake is associated with virulence in Staphylococcus [45].
AbgB encodes a carboxypeptidase that is a member of the M20 family of metallopeptidases [46]. This family includes a methicillin-resistant S. aureus antibiotic resistance factor, p-aminobenzoyl-glutamate hydrolase.
YdiL encodes a membrane protease. Proteases serve multiple functions in protecting Staphylococcus from their hosts’ immune defenses. They inhibit phagocyte-mediated killing, inactivate components of the immune system, and interfere with epithelial defenses [47].
Mitomycins, such as mitomycin C, are bifunctional alkylators that are able to kill bacteria by crosslinking complementary strands of DNA and preventing DNA replication. They accomplish this indirectly through the reduction in quinone. YcaQ protects bacteria from crosslinking compounds by unhooking the crosslinked DNA [48].
The ChrA protein, whose gene was identified in S. coagulans, serves as an efflux pump. It transports chromate from bacteria but also has a number of homologs, many of which are not well characterized [49]. Given the importance of efflux pumps in antimicrobial resistance, especially in biofilms, it is possible that ChrA or a homolog play this role in S. schleiferi [50].
These genes and their associated products represent potentially important virulence factors and mechanisms of resistance. Additional research will be required to determine what role, if any, these proteins play in staphylococcal infections and effective vaccines. Additionally, our study has resulted in the identification of numerous hypothetical proteins, highlighting the complexity of the genetic landscape in both S. schleiferi and S. coagulans. These uncharacterized proteins constitute a viable area for future research. Further characterization and investigation of these putative proteins (Table S1) has the potential to reveal novel processes, virulence factors, or adaptive properties that may contribute to the distinct characteristics of these Staphylococcus species. This research not only improves our understanding of the genetic differences between S. schleiferi and S. coagulans, but it also underscores the importance of continuing to explore the undiscovered parts of their genomes.
The identification of these unexplored proteins may provide a better understanding of the biological complexities that define the pathogenic potential and adaptive strategies of S. schleiferi and S. coagulans. Despite the useful insights garnered from this study, it is important to recognize its limitations. One significant limitation is the small number of S. schleiferi isolates included in the study. The majority of these isolates came from similar geographical areas, primarily the Netherlands and Australia. This geographical bias may have an impact on the overall diversity captured in the study, thereby limiting the generalizability of the findings to a broader global setting. Furthermore, the majority of the S. coagulans isolates were obtained from Scotland and the Netherlands. This geographic concentration poses challenges regarding the representation of the global genetic diversity within S. coagulans populations. Bacterial populations can vary significantly depending on factors such as host populations, the environment, and local selective pressures. As a result, the observed species-specific genes may not represent the whole range of genetic adaptations that may exist in S. schleiferi and S. coagulans populations from different geographical areas. The identified species-specific genes should be functionally characterized in the future, with an emphasis on their roles in pathogenesis and host interactions.
Comparative transcriptomics and proteomics could provide insights into gene expression patterns and protein functions, further elucidating the adaptive strategies employed by these Staphylococcal species.

5. Conclusions

This study adds to the field of bacterial genomics by identifying species-specific genes in Staphylococcus coagulans and Staphylococcus Schleifer. The identified genes act as molecular signatures, potentially providing the basis for improved diagnostic techniques. The study also emphasizes the need to consider the pangenome rather than a single reference genome in order to gain a more complete knowledge of genetic variation within and between species. While this work provides useful insights into the genomes of S. schleiferi and S. coagulans, the limitations due to geographic bias highlight the necessity of ongoing research efforts that include a more diverse and internationally representative collection of isolates.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/genes15030284/s1, Table S1: Gene clusters identified in S. schleiferi and S. coagulans isolates.

Author Contributions

Conceptualization, S.A.K. and M.A.A.; methodology, M.A.A.; formal analysis, M.A.A.; investigation, M.A.A.; resources, S.A.K.; data curation, M.A.A. and S.A.K.; writing—original draft preparation, M.A.A.; writing—review and editing, S.A.K.; visualization, M.A.A.; project administration, S.A.K.; funding acquisition, S.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the University of Tennessee Center of Excellence in Livestock Diseases and Human Health.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available in GenBank as noted.

Acknowledgments

We acknowledge the guidance of David Bemis in conducting this study.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Yamashita, K.; Shimizu, A.; Kawano, J.; Uchida, E.; Haruna, A.; Igimi, S. Isolation and characterization of staphylococci from external auditory meatus of dogs with or without otitis externa with special reference to Staphylococcus schleiferi subsp. coagulans isolates. J. Vet. Med. Sci. 2005, 67, 263–268. [Google Scholar] [CrossRef]
  2. May, E.R.; Hnilica, K.A.; Frank, L.A.; Jones, R.D.; Bemis, D.A. Isolation of Staphylococcus schleiferi from healthy dogs and dogs with otitis, pyoderma, or both. J. Am. Vet. Med. Assoc. 2005, 227, 928–931. [Google Scholar] [CrossRef]
  3. Paterson, G.K. Genomic epidemiology of the opportunistic pathogen Staphylococcus coagulans from companion dogs. J. Med. Microbiol. 2021, 70, 001407. [Google Scholar] [CrossRef]
  4. Kumar, D.; Cawley, J.J.; Irizarry-Alvarado, J.M.; Alvarez, A.; Alvarez, S. Case of Staphylococcus schleiferi subspecies coagulans endocarditis and metastatic infection in an immune compromised host. Transpl. Infect Dis. 2007, 9, 336–338. [Google Scholar] [CrossRef]
  5. Igimi, S.; Takahashi, E.; Mitsuoka, T. Staphylococcus schleiferi subsp. coagulans subsp. nov., isolated from the external auditory meatus of dogs with external ear otitis. Int. J. Syst. Bacteriol. 1990, 40, 409–411. [Google Scholar] [CrossRef]
  6. Naing, S.Y.; Duim, B.; Broens, E.M.; Schweitzer, V.; Zomer, A.; van der Graaf-van Bloois, L.; van der Meer, C.; Stellingwerff, L.; Fluit, A.C.; Wagenaar, J.A. Molecular Characterization and Clinical Relevance of Taxonomic Reassignment of Staphylococcus schleiferi Subspecies into Two Separate Species, Staphylococcus schleiferi and Staphylococcus coagulans. Microbiol. Spectr. 2023, 11, e0467022. [Google Scholar] [CrossRef]
  7. Madhaiyan, M.; Wirth, J.S.; Saravanan, V.S. Phylogenomic analyses of the Staphylococcaceae family suggest the reclassification of five species within the genus Staphylococcus as heterotypic synonyms, the promotion of five subspecies to novel species, the taxonomic reassignment of five Staphylococcus species to Mammaliicoccus gen. nov., and the formal assignment of Nosocomiicoccus to the family Staphylococcaceae. Int. J. Syst. Evol. Microbiol. 2020, 70, 5926–5936. [Google Scholar]
  8. Ashelford, K.E.; Chuzhanova, N.A.; Fry, J.C.; Jones, A.J.; Weightman, A.J. At least 1 in 20 16S rRNA sequence records currently held in public repositories is estimated to contain substantial anomalies. Appl. Environ. Microbiol. 2005, 71, 7724–7736. [Google Scholar] [CrossRef]
  9. Vandenesch, F.; Lebeau, C.; Bes, M.; Lina, G.; Lina, B.; Greenland, T.; Benito, Y.; Brun, Y.; Fleurette, J.; Etienne, J. Clotting activity in Staphylococcus schleiferi subspecies from human patients. J. Clin. Microbiol. 1994, 32, 388–392. [Google Scholar] [CrossRef]
  10. Lefort, V.; Desper, R.; Gascuel, O. FastME 2.0: A Comprehensive, Accurate, and Fast Distance-Based Phylogeny Inference Program. Mol. Biol. Evol. 2015, 32, 2798–2800. [Google Scholar] [CrossRef]
  11. Abouelkhair, M.A.; Rifkin, R.E.; Grzeskowiak, R.M.; Biris, A.S.; Anderson, D.E.; Bemis, D.A.; Kania, S.A. Complete Genome Sequences of Four Staphylococcus aureus Sequence Type 398 Isolates from Four Goats with Osteomyelitis. Microbiol. Resour. Announc. 2018, 6, e01599-17. [Google Scholar] [CrossRef]
  12. Billings, C.; Rifkin, R.; Abouelkhair, M.; Jones, R.D.; Bow, A.; Kolape, J.; Rajeev, S.; Kania, S.; Anderson, D.E. In vitro and in vivo assessment of caprine origin Staphylococcus aureus ST398 strain UTCVM1 as an osteomyelitis pathogen. Front. Cell Infect Microbiol. 2022, 12, 1015655. [Google Scholar] [CrossRef]
  13. Abebe, A.A.; Birhanu, A.G. Methicillin Resistant Staphylococcus aureus: Molecular Mechanisms Underlying Drug Resistance Development and Novel Strategies to Combat. Infect Drug Resist 2023, 16, 7641–7662. [Google Scholar] [CrossRef]
  14. Aung, M.S.; Urushibara, N.; Kawaguchiya, M.; Ohashi, N.; Hirose, M.; Kimura, Y.; Kudo, K.; Ito, M.; Kobayashi, N. Molecular Epidemiological Characterization of Methicillin-Resistant Staphylococcus aureus from Bloodstream Infections in Northern Japan: Increasing Trend of CC1 and Identification of ST8-SCCmec IVa USA300-Like Isolate Lacking Arginine Catabolic Mobile Element. Microb. Drug Resist 2023, 30, 63–72. [Google Scholar]
  15. Kim, Y.H.; Park, J.; Chung, H.S. Genetic characterization of tetracycline-resistant Staphylococcus aureus with reduced vancomycin susceptibility using whole-genome sequencing. Arch. Microbiol. 2023, 206, 24. [Google Scholar] [CrossRef]
  16. Ko, D.S.; Kim, N.H.; Kim, E.K.; Ha, E.J.; Ro, Y.H.; Kim, D.; Choi, K.S.; Kwon, H.J. Comparative genomics of bovine mastitis-origin Staphylococcus aureus strains classified into prevalent human genotypes. Res. Vet. Sci. 2021, 139, 67–77. [Google Scholar] [CrossRef]
  17. Makarova, O.; Johnston, P.; Rodriguez-Rojas, A.; El Shazely, B.; Morales, J.M.; Rolff, J. Genomics of experimental adaptation of Staphylococcus aureus to a natural combination of insect antimicrobial peptides. Sci. Rep. 2018, 8, 15359. [Google Scholar] [CrossRef]
  18. Abouelkhair, M.A.; Bemis, D.A.; Giannone, R.J.; Frank, L.A.; Kania, S.A. Identification, cloning and characterization of SpEX exotoxin produced by Staphylococcus pseudintermedius. PLoS ONE 2019, 14, e0220301. [Google Scholar] [CrossRef]
  19. Abouelkhair, M.A.; Frank, L.A.; Bemis, D.A.; Giannone, R.J.; Kania, S.A. Staphylococcus pseudintermedius 5’-nucleotidase suppresses canine phagocytic activity. Vet. Microbiol. 2020, 246, 108720. [Google Scholar] [CrossRef]
  20. Abouelkhair, M.A.; Thompson, R.; Riley, M.C.; Bemis, D.A.; Kania, S.A. Complete Genome Sequences of Three Staphylococcus pseudintermedius Strains Isolated from Botswana. Genome Announc. 2018, 6, 10-128. [Google Scholar] [CrossRef]
  21. Phophi, L.; Abouelkhair, M.; Jones, R.; Henton, M.; Qekwana, D.N.; Kania, S.A. The molecular epidemiology and antimicrobial resistance of Staphylococcus pseudintermedius canine clinical isolates submitted to a veterinary diagnostic laboratory in South Africa. PLoS ONE 2023, 18, e0290645. [Google Scholar] [CrossRef]
  22. Griffeth, G.C.; Morris, D.O.; Abraham, J.L.; Shofer, F.S.; Rankin, S.C. Screening for skin carriage of methicillin-resistant coagulase-positive staphylococci and Staphylococcus schleiferi in dogs with healthy and inflamed skin. Vet. Dermatol. 2008, 19, 142–149. [Google Scholar] [CrossRef]
  23. Huse, H.K.; Miller, S.A.; Chandrasekaran, S.; Hindler, J.A.; Lawhon, S.D.; Bemis, D.A.; Westblade, L.F.; Humphries, R.M. Evaluation of Oxacillin and Cefoxitin Disk Diffusion and MIC Breakpoints Established by the Clinical and Laboratory Standards Institute for Detection of mecA-Mediated Oxacillin Resistance in Staphylococcus schleiferi. J. Clin. Microbiol. 2018, 56, e01653-17. [Google Scholar] [CrossRef]
  24. Roberts, S.; O’Shea, K.; Morris, D.; Robb, A.; Morrison, D.; Rankin, S. A real-time PCR assay to detect the Panton Valentine Leukocidin toxin in staphylococci: Screening Staphylococcus schleiferi subspecies coagulans strains from companion animals. Vet. Microbiol. 2005, 107, 139–144. [Google Scholar] [CrossRef]
  25. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef]
  26. Gurevich, A.; Saveliev, V.; Vyahhi, N.; Tesler, G. QUAST: Quality assessment tool for genome assemblies. Bioinformatics 2013, 29, 1072–1075. [Google Scholar] [CrossRef]
  27. Meier-Kolthoff, J.P.; Goker, M. TYGS is an automated high-throughput platform for state-of-the-art genome-based taxonomy. Nat. Commun. 2019, 10, 2182. [Google Scholar] [CrossRef]
  28. Meier-Kolthoff, J.P.; Carbasse, J.S.; Peinado-Olarte, R.L.; Goker, M. TYGS and LPSN: A database tandem for fast and reliable genome-based classification and nomenclature of prokaryotes. Nucleic. Acids Res. 2022, 50, D801–D807. [Google Scholar] [CrossRef]
  29. Ondov, B.D.; Treangen, T.J.; Melsted, P.; Mallonee, A.B.; Bergman, N.H.; Koren, S.; Phillippy, A.M. Mash: Fast genome and metagenome distance estimation using MinHash. Genome Biol. 2016, 17, 132. [Google Scholar] [CrossRef]
  30. Lagesen, K.; Hallin, P.; Rodland, E.A.; Staerfeldt, H.H.; Rognes, T.; Ussery, D.W. RNAmmer: Consistent and rapid annotation of ribosomal RNA genes. Nucleic Acids Res. 2007, 35, 3100–3108. [Google Scholar] [CrossRef]
  31. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST+: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef]
  32. Meier-Kolthoff, J.P.; Auch, A.F.; Klenk, H.P.; Goker, M. Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinform. 2013, 14, 60. [Google Scholar] [CrossRef]
  33. Farris, J.S. Estimating Phylogenetic Trees from Distance Matrices. Am. Nat. 1972, 106, 645–668. [Google Scholar] [CrossRef]
  34. Kreft, L.; Botzki, A.; Coppens, F.; Vandepoele, K.; Van Bel, M. PhyD3: A phylogenetic tree viewer with extended phyloXML support for functional genomics data visualization. Bioinformatics 2017, 33, 2946–2947. [Google Scholar] [CrossRef]
  35. Meier-Kolthoff, J.P.; Hahnke, R.L.; Petersen, J.; Scheuner, C.; Michael, V.; Fiebig, A.; Rohde, C.; Rohde, M.; Fartmann, B.; Goodwin, L.A.; et al. Complete genome sequence of DSM 30083(T), the type strain (U5/41(T)) of Escherichia coli, and a proposal for delineating subspecies in microbial taxonomy. Stand. Genomic Sci. 2014, 9, 2. [Google Scholar] [CrossRef]
  36. Eren, A.M.; Kiefl, E.; Shaiber, A.; Veseli, I.; Miller, S.E.; Schechter, M.S.; Fink, I.; Pan, J.N.; Yousef, M.; Fogarty, E.C.; et al. Community-led, integrated, reproducible multi-omics with anvi’o. Nat. Microbiol. 2021, 6, 3–6. [Google Scholar] [CrossRef]
  37. Pritchard, L.; Glover, R.H.; Humphris, S.; Elphinstone, J.G.; Toth, I.K. Genomics and taxonomy in diagnostics for food security: Soft-rotting enterobacterial plant pathogens. Anal. Methods 2016, 8, 12–24. [Google Scholar] [CrossRef]
  38. Pritchard, L.; Harrington, B.; Cock, P.; Davey, R.; Waters, N.; Sharma, A.; Balázs Brankovics, Y.T.; Esen, Ö. widdowquinn/pyani: Application and Python Module. 2022. Available online: https://zenodo.org/record/6927794 (accessed on 20 December 2023).
  39. Lee, I.; Ouk Kim, Y.; Park, S.-C.; Chun, J. OrthoANI: An improved algorithm and software for calculating average nucleotide identity. Int. J. Syst. Evol. Microbiol. 2016, 66, 1100–1103. [Google Scholar] [CrossRef]
  40. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef]
  41. Page, A.J.; Cummins, C.A.; Hunt, M.; Wong, V.K.; Reuter, S.; Holden, M.T.G.; Fookes, M.; Falush, D.; Keane, J.A.; Parkhill, J. Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 2015, 31, 3691–3693. [Google Scholar] [CrossRef]
  42. Kaas, R.S.; Leekitcharoenphon, P.; Aarestrup, F.M.; Lund, O. Solving the problem of comparing whole bacterial genomes across different sequencing platforms. PLoS ONE 2014, 9, e104984. [Google Scholar] [CrossRef]
  43. Hammond, A.J.; Binsker, U.; Aggarwal, S.D.; Ortigoza, M.B.; Loomis, C.; Weiser, J.N. Neuraminidase B controls neuraminidase A-dependent mucus production and evasion. PLoS Pathog. 2021, 17, e1009158. [Google Scholar] [CrossRef]
  44. Osman, D.; Piergentili, C.; Chen, J.; Sayer, L.N.; Uson, I.; Huggins, T.G.; Robinson, N.J.; Pohl, E. The Effectors and Sensory Sites of Formaldehyde-responsive Regulator FrmR and Metal-sensing Variant. J. Biol. Chem. 2016, 291, 19502–19516. [Google Scholar] [CrossRef]
  45. Heilbronner, S.; Holden, M.T.; van Tonder, A.; Geoghegan, J.A.; Foster, T.J.; Parkhill, J.; Bentley, S.D. Genome sequence of Staphylococcus lugdunensis N920143 allows identification of putative colonization and virulence factors. FEMS Microbiol. Lett. 2011, 322, 60–67. [Google Scholar] [CrossRef]
  46. Jamdar, S.N.; Are, V.N.; Navamani, M.; Kumar, S.; Nagar, V.; Makde, R.D. The members of M20D peptidase subfamily from Burkholderia cepacia, Deinococcus radiodurans and Staphylococcus aureus (HmrA) are carboxydipeptidases, primarily specific for Met-X dipeptides. Arch Biochem. Biophys. 2015, 587, 18–30. [Google Scholar] [CrossRef]
  47. Singh, V.; Phukan, U.J. Interaction of host and Staphylococcus aureus protease-system regulates virulence and pathogenicity. Med. Microbiol. Immunol. 2019, 208, 585–607. [Google Scholar] [CrossRef]
  48. Bradley, N.P.; Washburn, L.A.; Christov, P.P.; Watanabe, C.M.H.; Eichman, B.F. Escherichia coli YcaQ is a DNA glycosylase that unhooks DNA interstrand crosslinks. Nucleic Acids Res. 2020, 48, 7005–7017. [Google Scholar]
  49. Jimenez-Mejia, R.; Campos-Garcia, J.; Cervantes, C. Membrane topology of the chromate transporter ChrA of Pseudomonas aeruginosa. FEMS Microbiol. Lett. 2006, 262, 178–184. [Google Scholar] [CrossRef]
  50. Soto, S.M. Role of efflux pumps in the antibiotic resistance of bacteria embedded in a biofilm. Virulence 2013, 4, 223–229. [Google Scholar] [CrossRef]
Genes 15 00284 g001
Figure 1. Tree inferred with FastME 2.1.6.1 [10] from GBDP distances calculated from genome sequences. The branch lengths are scaled in terms of GBDP distance. The tree was rooted at the midpoint [33].
Figure 1. Tree inferred with FastME 2.1.6.1 [10] from GBDP distances calculated from genome sequences. The branch lengths are scaled in terms of GBDP distance. The tree was rooted at the midpoint [33].
Genes 15 00284 g001
Genes 15 00284 g002
Figure 2. Static images from the Anvi’o-8 interactive display for the S. schleiferi and S. coagulans dataset.
Figure 2. Static images from the Anvi’o-8 interactive display for the S. schleiferi and S. coagulans dataset.
Genes 15 00284 g002
Genes 15 00284 g003
Figure 3. Phylogenetic tree based on the SNP analysis conducted in CSI Phylogeny 1.4 comparing S. schleiferi 192, S. schleiferi 205, S. schleiferi 214, and S. schleiferi ATCC 43808T genomes with 46 publicly available S. schleiferi genomes.
Figure 3. Phylogenetic tree based on the SNP analysis conducted in CSI Phylogeny 1.4 comparing S. schleiferi 192, S. schleiferi 205, S. schleiferi 214, and S. schleiferi ATCC 43808T genomes with 46 publicly available S. schleiferi genomes.
Genes 15 00284 g003
Genes 15 00284 g004
Figure 4. Pangenome analysis of 182 S. coagulans genomes using Roary. The left panel represents the Bayesian phylogenetic tree based on S. coagulans core gene single-nucleotide polymorphisms (SNPs) and Roary matrix. The right panel represents the matrix where the accessory and core genes were either present or absent. Isolates from different countries are color-coded as follows: Netherlands: #e4524a, USA: #5e4fa2, Scotland: #f2f9a9, Thailand: #f67d47, South Korea: #bce6a0, United Kingdom: #75c8a4, Japan: #1e90bf, Russia: #fdbf6c, and India: #fdeea1.
Figure 4. Pangenome analysis of 182 S. coagulans genomes using Roary. The left panel represents the Bayesian phylogenetic tree based on S. coagulans core gene single-nucleotide polymorphisms (SNPs) and Roary matrix. The right panel represents the matrix where the accessory and core genes were either present or absent. Isolates from different countries are color-coded as follows: Netherlands: #e4524a, USA: #5e4fa2, Scotland: #f2f9a9, Thailand: #f67d47, South Korea: #bce6a0, United Kingdom: #75c8a4, Japan: #1e90bf, Russia: #fdbf6c, and India: #fdeea1.
Genes 15 00284 g004
Genes 15 00284 g005
Figure 5. Pangenome analysis for 50 S. schleiferi genomes using Roary. The left side represents the Bayesian phylogenetic tree based on S. schleiferi core gene single-nucleotide polymorphisms (SNPs) and Roary matrix. The right side represents the matrix where the accessory and core genes were either present or absent. Isolates from different countries are color-coded as follows: Netherlands:#e4524a, Australia: #eef8a2, USA: #00a4bd, France: #a1d9a4, United Kingdom: #5e4fa2, Switzerland: #fde999 and Not provided: #fca55b. Isolates from different hosts are color-coded as follows: Homo sapiens: #5e4fa2, Camel: #bce6a0, and Not provided: #fdbf6c.
Figure 5. Pangenome analysis for 50 S. schleiferi genomes using Roary. The left side represents the Bayesian phylogenetic tree based on S. schleiferi core gene single-nucleotide polymorphisms (SNPs) and Roary matrix. The right side represents the matrix where the accessory and core genes were either present or absent. Isolates from different countries are color-coded as follows: Netherlands:#e4524a, Australia: #eef8a2, USA: #00a4bd, France: #a1d9a4, United Kingdom: #5e4fa2, Switzerland: #fde999 and Not provided: #fca55b. Isolates from different hosts are color-coded as follows: Homo sapiens: #5e4fa2, Camel: #bce6a0, and Not provided: #fdbf6c.
Genes 15 00284 g005
Table 1. Accession numbers for the genome sequences retrieved from GenBank, utilized in the Anvio-8 analysis of S. schleiferi and S. coagulans isolates.
Table 1. Accession numbers for the genome sequences retrieved from GenBank, utilized in the Anvio-8 analysis of S. schleiferi and S. coagulans isolates.
IsolateGenBank Accession Number
S. schleferi MGYG-HGUT-01437GCA_902374935.1
S. coagulans 2317-03GCA_001188915.1
S. coagulans OT1-1GCA_004026205.1
S. coagulans 5909-02GCA_001188875.1
S. coagulans TSCC54GCA_001548255.1
S. coagulans reference genomeGCA_022557135.1
S. coagulans 2142-05GCA_001188895.1
All genomes were downloaded from NCBI (ftp://ftp.ncbi.nlm.nih.gov, accessed on 20 December 2023).
Table 2. Genomic features of human Staphylococcus schleiferi.
Table 2. Genomic features of human Staphylococcus schleiferi.
StrainWGS Accession No *No of ContigsN50 (bp)Genome Length (bp)G+C Content (%)Predicted Coding SequencesPredicted RNAs
Staphylococcus coagulans 191PNRJ0000000051138,8932,508,13335.73229472
Staphylococcus schleiferi 192POVG0000000010259,7862,452,48735.87220374
Staphylococcus coagulans 196POVH0000000056110,2792,508,60435.74229974
Staphylococcus schleiferi 205POVI000000009256,9582,468,34235.92221876
Staphylococcus schleiferi 214POVJ000000009457,2472,469,69935.92221676
Staphylococcus schleiferi ATCC 43808TPOVK000000008856,9382,469,63835.92221873
* WGS—whole genome sequence.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Abouelkhair, M.A.; Kania, S.A. Whole Genome Sequencing and Comparative Genomics of Six Staphylococcus schleiferi and Staphylococcus coagulans Isolates. Genes 2024, 15, 284. https://doi.org/10.3390/genes15030284

AMA Style

Abouelkhair MA, Kania SA. Whole Genome Sequencing and Comparative Genomics of Six Staphylococcus schleiferi and Staphylococcus coagulans Isolates. Genes. 2024; 15(3):284. https://doi.org/10.3390/genes15030284

Chicago/Turabian Style

Abouelkhair, Mohamed A., and Stephen A. Kania. 2024. "Whole Genome Sequencing and Comparative Genomics of Six Staphylococcus schleiferi and Staphylococcus coagulans Isolates" Genes 15, no. 3: 284. https://doi.org/10.3390/genes15030284

APA Style

Abouelkhair, M. A., & Kania, S. A. (2024). Whole Genome Sequencing and Comparative Genomics of Six Staphylococcus schleiferi and Staphylococcus coagulans Isolates. Genes, 15(3), 284. https://doi.org/10.3390/genes15030284

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop