Genomic analysis of prevalent Staphylococcus epidermidis multidrug-resistant strains isolated during eight years in a single children hospital in México City.

Laboratorio de Patogenicidad Bacteriana, Departamento de Salud Pública, Facultad de Medicina, Universidad Nacional Autónoma de México, Ciudad de México, México
Centro de Ciencias Genómicas, Universidad Nacional Autónoma de México, Cuernavaca, Morelos, México
Instituto Nacional de Medicina Genómica, Mexico City, México
DOI
10.7287/peerj.preprints.27693v2
Subject Areas
Genomics, Microbiology, Infectious Diseases
Keywords
Staphylococcus epidermidis, genomes, pangenome, prophages, CRISPR, Recombination, antibiotic resistance, insertion sequences, clonal structure
Copyright
© 2019 Cabrera-Contreras et al.
Licence
This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, reproduction and adaptation in any medium and for any purpose provided that it is properly attributed. For attribution, the original author(s), title, publication source (PeerJ Preprints) and either DOI or URL of the article must be cited.
Cite this article
Cabrera-Contreras R, Santamaría RI, Bustos P, Martínez-Flores I, Meléndez E, Morelos R, Barbosa-Amezcua M, González-Covarrubias V, Silva-Herzog E, Soberón X, González V. 2019. Genomic analysis of prevalent Staphylococcus epidermidis multidrug-resistant strains isolated during eight years in a single children hospital in México City. PeerJ Preprints 7:e27693v2

Abstract

Staphylococcus epidermidis is a human commensal and pathogen worldwide distributed. In this work, we surveyed for multi-resistant S. epidermidis strains in eight years at a children health-care unit in México City. Multidrug-resistant S. epidermidis were present in all years of the study. Resistance to methicillin, beta-lactams, fluoroquinolones, and macrolides were included. To understand the genetic basis of antibiotic resistance and its association with virulence and gene exchange, we sequenced the genomes of 17 S. epidermidis isolates. Whole-genome nucleotide identities between all the pairs of S. epidermidis strains were about 97% to 99%. We inferred a clonal structure and eight Multilocus Sequence Types (MLST´s) in the S. epidermidis sequenced collection. The profile of virulence includes genes involved in biofilm formation and phenol-soluble modulins (PSMs). However, half of the S. epidermidis analyzed lacked the icaoperon for biofilm formation. Likely, they are commensal S. epidermidis strains but multi-antibiotic resistant. Uneven distribution of insertion sequences, phages, and CRISPR-Cas immunity phage systems suggest frequent horizontal gene transfer. Rates of recombination between S. epidermidis strains were more prevalent than the mutation rate and affected the whole genome. Therefore, the multidrug-resistance, independently of the pathogenic traits, might explain the persistence of specific highly adapted S. epidermidis clonal lineages in nosocomial settings.

Author Comment

In the previous version, we try to answer a general question about what are the S. epidermidis flowing in the children hospital during 8 years. The reviewers main concern was that we use a limited number of isolates and how these isolates were selected. We also lack of a clear explanation of why we performed the experiments with multidrug-resistant strains. Then, the question was not adequately responded. These and other details were addressed in the new reviewed version.

The new version focus on the study of the frequency of the multi resistant isolates during eight years, and if these isolates correspond to one or few clones, and if they are pathogenic or commensal. We added data on a survey of staphylococci along eight years. The results showed that S. epidermidis multi resistant to more than 8 antibiotics were the most common found in these years. Our results indicate that multiresistance is a property shared by pathogenic and commensal isolates. Then, it could the factor that explain the permanency of clones along thee years in the same hospital setting.

Supplemental Information

Completeness of SE INPer genome collection

DOI: 10.7287/peerj.preprints.27693v2/supp-2

Genomes of SEfrom the GenBank used in this work

DOI: 10.7287/peerj.preprints.27693v2/supp-3

Prophage sequences and their homologs found in SEINPer genomes

DOI: 10.7287/peerj.preprints.27693v2/supp-4

CRISPR-Cas elements present in SE INPer strains

DOI: 10.7287/peerj.preprints.27693v2/supp-5

Insertion sequences present in SEINPer strains

DOI: 10.7287/peerj.preprints.27693v2/supp-6

Survey of staphylococci at Instituto Nacional de Perinatología, México City Alone eight years

A. Species classification and proportion. B. Origin of the isolates and proportion. C. Isolation sites proportion.

DOI: 10.7287/peerj.preprints.27693v2/supp-7

Frequency cumulation of antibiotic resistant strains in Staphylococcus

A. S. epidermidis. B. S. aureus. The absolute number of antibiotic resistances (x-axis) by the number of strains (y-axis) was counted from 2006 to 2013.

DOI: 10.7287/peerj.preprints.27693v2/supp-8

Genome identity (ANIm) between pairs of SE INPer strains

Pairwise whole genome alignments were done with Mummer within the JSspecies program (Richter et al., 2016). The percent of average nucleotide alignment (ANI) was illustrated by a heat-map constructed with ggPlot2 in R (see methods). ANI > 99% are in red color. S. epidermidis ATCC 12228 was included for comparison. S10 strain had ANI = 97% respect all the other SE strains.

DOI: 10.7287/peerj.preprints.27693v2/supp-9

Pangenome model of SE strains

The pangenome model of 29 SE strains was performe with GET_HOMOLOGUES (Contreras-Moreira and Vinuesa, 2013) as described in methods. A. Pangenome size (number of gene family clusters, Y axis) as a function of the number of SE genomes (X axis). B. Core genome according to the Tettelin equations.

DOI: 10.7287/peerj.preprints.27693v2/supp-10

Clonal relationships of the STs detected in INPer strains respect to the ST database

Alleles for the seven proteins used in the S. epidermidisMLST scheme (Thomas et al., 2007) were looked at the Staphylococcus epidermidis MLST database ( https://pubmlst.org/sepidermidis/ ; Table 1; see methods) (Feil et al., 2004). The clonal relationships among STs were determined by eBURST( http://eburst.mlst.net ). Six out of 8 ST complexes assigned to the SEINPer strains are denoted by numbers in violet color.

DOI: 10.7287/peerj.preprints.27693v2/supp-11

Phylogenetic tree of the SCCmec recombinases

A. ccrA. B. ccrB. C. ccrB

DOI: 10.7287/peerj.preprints.27693v2/supp-12

Proposed structure of the combined cassette SCCmec IV-VIII

The segment corresponding to SCCmec cassettes contained within the contig 12 of the S07 strain is annotated according to the best blastN matches against the nr database of the Genebank

DOI: 10.7287/peerj.preprints.27693v2/supp-13