Front. Cell. Infect. Microbiol.Frontiers in Cellular and Infection MicrobiologyFront. Cell. Infect. Microbiol.2235-2988Frontiers Media S.A.10.3389/fcimb.2023.1194254Cellular and Infection MicrobiologyOriginal ResearchAtopic dermatitis-derived Staphylococcus aureus strains: what makes them special in the interplay with the hostConteAntonietta Lucia1†BrunettiFrancesca1†MarazzatoMassimiliano1*LonghiCatia1MauriziLinda1RaponiGiammarco1PalamaraAnna Teresa23GrassiSara4‡ConteMaria Pia1‡1Department of Public Health and Infectious Diseases, Sapienza University of Rome, Rome, Italy2Department of Infectious Diseases, Istituto Superiore di Sanità, Rome, Italy3Department of Public Health and Infectious Diseases, Sapienza University of Rome, Laboratory Affiliated to Institute Pasteur Italia- Cenci Bolognetti Foundation, Rome, Italy4Dermatology Clinic, Department of Clinical Internal, Anesthesiological and Cardiovascular Sciences, Sapienza University of Rome, Rome, Italy
Edited by: Marat R. Sadykov, University of Nebraska Medical Center, United States
Reviewed by: Roman J. Nowicki, Medical University of Gdansk, Poland; Vineeth Rajan, Central Institute of Fisheries Technology (ICAR), India
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Background
Atopic dermatitis (AD) is a chronic inflammatory skin condition whose pathogenesis involves genetic predisposition, epidermal barrier dysfunction, alterations in the immune responses and microbial dysbiosis. Clinical studies have shown a link between Staphylococcus aureus and the pathogenesis of AD, although the origins and genetic diversity of S. aureus colonizing patients with AD is poorly understood. The aim of the study was to investigate if specific clones might be associated with the disease.
Methods
WGS analyses were performed on 38 S. aureus strains, deriving from AD patients and healthy carriers. Genotypes (i.e. MLST, spa-, agr- and SCCmec-typing), genomic content (e.g. virulome and resistome), and the pan-genome structure of strains have been investigated. Phenotypic analyses were performed to determine the antibiotic susceptibility, the biofilm production and the invasiveness within the investigated S. aureus population.
Results
Strains isolated from AD patients revealed a high degree of genetic heterogeneity and a shared set of virulence factors and antimicrobial resistance genes, suggesting that no genotype and genomic content are uniquely associated with AD. The same strains were characterized by a lower variability in terms of gene content, indicating that the inflammatory conditions could exert a selective pressure leading to the optimization of the gene repertoire. Furthermore, genes related to specific mechanisms, like post-translational modification, protein turnover and chaperones as well as intracellular trafficking, secretion and vesicular transport, were significantly more enriched in AD strains. Phenotypic analysis revealed that all of our AD strains were strong or moderate biofilm producers, while less than half showed invasive capabilities.
Conclusions
We conclude that in AD skin, the functional role played by S. aureus may depend on differential gene expression patterns and/or on post-translational modification mechanisms rather than being associated with peculiar genetic features.
Atopic dermatitis (AD) is a chronic inflammatory skin disease, characterized by intense itching and recurrent eczematous lesions, that affect 15-20% of children and 2-5% of adults in industrialized countries (Sybilski et al., 2015; Hadi et al., 2021). It usually occurs during early childhood with onset in 45% of cases within the first six months of the child’s life and 60% within the first year of life (Weidinger et al., 2018). There are many contributors to the pathogenesis of AD, including genetic defects in the epidermal barrier that can be explained, at least in part, by inherited mutations in keratinocyte proteins such as filaggrin and loricrin (Kim et al., 2008; Howell et al., 2009; Callewaert et al., 2020), an imbalance towards a T-helper-2 (Th2) immune response in which the cytokines IL-4, IL-5 and IL-13 play a major role (Homey et al., 2006; Brandt and Sivaprasad, 2011), and microbial dysbiosis that increases the susceptibility to AD (Bjerre et al., 2017). Numerous studies have shown that over 90% of patients with AD have Staphylococcus aureus skin colonization (Bieber, 2008; Kong et al., 2012; Totté et al., 2016; Meylan et al., 2017; Altunbulakli et al., 2018).
S. aureus is an important human pathogen in both community and hospital settings but it can also establish a commensal relationship with humans. The anterior nares represent the main habitat for S. aureus, which persistently colonized 20–30% of the population (Peacock et al., 2001; van Belkum et al., 2009). Metagenomic analysis has shown that the disease activity is significantly associated with the dysbiosis of skin microbiota, which is characterized by a reduced microbial diversity during eczema flare, allowing the proliferation of S. aureus (Kong et al., 2012). The most colonized sites are the lesional skin (56–96.2%), the nose (46.1–64.1%), and the non-lesional skin (28–39%). Interestingly, 65–77.3% of AD patients were colonized both in the anterior nares and on the skin, whereas only 10.2% of healthy control subjects were colonized on the skin (Na et al., 2012; Totté et al., 2016; Alsterholm et al., 2017; Clausen et al., 2017). Moreover, some reports found that 55% of AD patients were persistent carriers of S. aureus (Alsterholm et al., 2017). However, its role in the pathogenesis of AD remains poorly understood. S. aureus employs multiple mechanisms to colonize the skin. They can produce adhesins, such as i) fibronectin-binding proteins (FnBP) A and B, expressed on the surface to promote the bacterial adherence to fibronectin (Cho et al., 2001), ii) the cell-wall-anchored (CWA) protein ClfB, that binds to the cornified envelope proteins loricrin and cytokeratin (Leonard et al., 2019), iii) the CWA protein IsdA, that renders the bacteria much more resistant to the inhibitory effects of the fatty acids, a property that may contribute to colonization of AD skin (Clarke et al., 2007), and iv) the S. aureus Protein A (spa), involved in pro-inflammatory responses (Gómez et al., 2004). Recently, it has been evidenced a predominance of S. aureus strains forming biofilms in AD lesions and this has been directly related to disease severity (Di Domenico et al., 2018). Biofilm represents a survival strategy protecting bacteria from antimicrobials and immune-mediated clearance (Brandwein et al., 2016). S. aureus expresses a plethora of secreted virulence factors (VF) that contribute to the epidermal barrier disruption, such as the phenol-soluble modulins (PSMs) peptides, that includes the PMSα and PMSβ families along with the δ-toxin. These factors are highly cytotoxic to a wide variety of cells (e.g. keratinocytes), contribute to the biofilm development (which is important for staphylococcal colonization and persistence) and trigger the skin inflammation upon the epidermal colonization (Peschel and Otto, 2013; Nakagawa et al., 2017). PSMs are regulated by the quorum sensing system, which is called Accessory Gene Regulator (agr) (Le and Otto, 2015). This regulation is essential for the timing of VFs expression during infection (Nakamura et al., 2020). Superantigens, such as the staphylococcal enterotoxin (SE) SEA, SEB, SEC, SED, and toxic shock syndrome toxin-1 (TSST-1) were the most common toxin genes carried by S. aureus isolated from AD patients (Kim et al., 2009). TSST-1, SEA, and SEB activate polyclonal T cells and subsequently cause T cell-mediated inflammation in AD lesions (Park et al., 2016). About 34% of the patients were colonized by α-toxin-producing S. aureus. The pore-forming α-toxin is likely to play an important role in disrupting the skin barrier by inducing keratinocyte cell death and Th2 cytokine production (Breuer et al., 2005). S. aureus also produces extracellular proteases that contribute to primary sensitization to allergens, abrogating the epidermal permeability barrier (Hirasawa et al., 2010). In recent years, it has become evident that S. aureus is a facultative intracellular pathogen, able to invade and survive in a range of cell types. Major invasion factors include the fibronectin-binding proteins (FnBPs), staphylococcal autolysin (Atl), and the lipoprotein-like (lpl) gene cluster (Jevon et al., 1999; Sinha and Fraunholz, 2010; Nguyen et al., 2015; Huitema et al., 2021). Intracellular localization, protecting S. aureus from the immune system and most antimicrobial treatments, is associated with chronic or relapsing staphylococcal infections (Horn et al., 2018). Moreover, among S. aureus strains colonizing AD patients, the percentage of methicillin-resistant S. aureus is 4–13 times higher than in a healthy population; however, methicillin-resistant S. aureus isolates seems to be much lower in adult patients as compared to the children (Kim et al., 2009; Lo et al., 2010; Pascolini et al., 2011; Rojo et al., 2014). Typing of S. aureus at the genetic level by multi-locus sequence typing (MLST) revealed that ST188, ST1, ST5, and ST513 were the most frequent STs identified in AD patients. Furthermore, clonal complexes (CCs) (based on the spa-type) showed that CC1 was associated with AD patients and was more prevalent in those carrying mutations in the gene encoding filaggrin (FLG) (Clausen et al., 2017), whereas CC30 was more common in healthy subjects (Clausen et al., 2019). However, the distribution of STs and staphylococcal CCs in AD patients points to significant heterogeneity, and no specific clone/clones prevailed in this group of patients. Conventional AD treatments involve emollients (medical moisturizers) and topical anti-inflammatory corticosteroids and calcineurin inhibitors (Carr, 2013). A plethora of pathogenesis-based treatments have been developed, including monoclonal antibodies that neutralize IgE or block specific interleukins or their receptors. The most advanced therapy is the IL-4 receptor blocker dupilumab which inhibits the function of IL-4 and IL-13, essential mediators of the Th2 pathway (Callewaert et al., 2020). In this study, we conducted a genotypic and phenotypic characterization of S. aureus isolated from AD patients during eczema flares and in the post-flare recovery phase, after treatment with dupixent® (dupilumab), to deepen our current knowledge about the features of the S. aureus strains associated with the disease.
Materials and methods
Fifteen (7 females and 8 males) patients, mean age of 30.2 ± 8.9 years old (range: 20-45), were enrolled at the Clinic of Dermatology, Department of Clinical Internal, Anesthesiologic and Cardiovascular Sciences, “Sapienza” University of Rome from September 2019 to February 2021.
The inclusion criteria to the study were the following: i) Caucasian Ethnicity; ii) age between 20 and 45 years old; iii) moderate-to-severe disease severity assessed through the Scoring Atopic Dermatitis “score index” (SCORAD); iv) eligibility for treatment with dupixent® (dupilumab).
The exclusion criteria were: i) local/systemic antibiotics, immunomodulatory or immunosuppressive therapies within 4 weeks prior to enrollment; ii) phototherapy treatments (nb-UVB, PUVA) within 4 weeks prior to enrollment; iii) participation to research trials with experimental treatments with other monoclonal antibodies; iv) past medical history of bone marrow transplant, intestinal diseases, diabetes mellitus or congenital immunodeficiencies.
Disease severity assessment and sample collection
At the time of enrollment, after signing the informed consent to the study, all the patients were subjected to clinical evaluation to assess the study eligibility and the degree of disease severity, according to the SCORAD index. The clinical diagnosis of AD was based on historical features, morphology and distribution of skin lesions, and associated clinical signs. The SCORAD index range lies between 0 and 103. Based on the SCORAD index results, AD patients have been classified into mild (≤25), moderate (26–50), and severe (≥ 51) forms.
Sampling was performed for each patient at two different stages of the disease: acute phase, characterized by an exacerbation of the disease on any skin region before starting treatment (t0); post-acute phase defined at least 4 weeks after reaching clinical remission with treatment (t1). All AD patients enrolled received dupixent® (dupilumab) as per protocol (600 mg subcutaneous injection once, then 300 mg every 2 weeks). Twenty-four hours before all sampling time points, patients had to avoid washing and the use of topical emollients at the sampling skin site. A sterile plastic device consisting of the open end of a cylinder with a diameter of 8 cm was applied to the skin area to be sampled. All samples were taken at the level of the popliteal fossa, which is divided into two parts (lesional and non-lesional) and each area was rubbed for 30 s with a sterile swab moistened with 1x PBS. At the same time, sampling at the level of the nasal choana using a sterile swab moistened with 1x PBS was carried out.
<italic>S. aureus</italic> isolation, identification and collection
Skin swabs were placed in tubes containing 3 ml of tryptic soy broth (TSB, Oxoid, Italy), immediately capped, and shaken for 30 s using a vortex mixer to suspend the bacterial cells. A volume of 100 µl of each swab suspension was plated onto selective Columbia agar with 5% of sheep blood (CAP, Oxoid Italy), Mannitol Salt Agar (MSA, Oxoid Italy), and Tryptic Soy Agar (TSA, Oxoid Italy) plates. Nasal swabs were immediately plated on an elective and selective growth medium. After incubation, the bacterial colonies were identified and counted.
From those samples where S. aureus was detected, a single colony was amplified and stored in a maintenance freeze medium (15% glycerol) (Oxoid) at −80 °C until use. The bacterial load from infection sites was determined by the conventional standard plate count method. The bacterial identification was performed by MALDI-TOF (Matrix-Assisted Laser Desorption Ionization Time-Of-Flight, VITEK® MS, BioMérieux) at the Laboratory of Microbiology Analysis of the “Policlinico Umberto I” Hospital of Rome. Seventeen S. aureus strains, previously isolated from the anterior nares of healthy carriers and part of the in-house culture collection, were used as control strains.
DNA extraction
S. aureus isolates from AD patients and healthy carriers were grown on TSA plates for 24 h at 37°C. Then, a single colony was inoculated in TSB and incubated for 18 h at 37°C. After incubation, 1.5 ml of each bacterial culture was centrifuged at 10000 rpm for 10 min.
DNA extraction on pelleted bacteria was conducted using “QIAamp DNA Mini Kit” (Qiagen, Germany), as the manufacturer’s instructions.
Random amplified polymorphic DNA-polymerase chain reaction
RAPD-PCR reactions were performed using the OLP 13 primer (5’-ACCGCCTGCT-3’) (Zare et al., 2019), and a 25μl reaction mix, containing 3mM MgCl2 (Bio Line, London, UK), 0.5 mM dNTPs (New England BioLabs, UK), 1 µM of primer, 1 U of Taq DNA polymerase (Bio Taq DNA polymerase, Meridian Bioscience, USA) in 1X PCR buffer (Meridian Bioscience, USA) and 3 μl of template dsDNA, was prepared. PCR was performed as follows: 1 cycle at 94°C for 5 min, followed by 40 cycles at 93°C for 1 min, 37°C for 1.5 min, and 72°C for 1 min, with a final extension at 72°C for 7 min (Mastercycler pro, Eppendorf, Germany). Amplification products were solved by electrophoresis on 1% agarose-TBE gel, stained with 5µl of “Midori Green Advance DNA Stain” (Nippon Genetics Europe GmbH, Germany), visualized under UV transilluminator (UVP Inc., Cambridge) and photographed with DigiDoc-It (UVP, Cambridge) system. To guarantee the reproducibility of RAPD-PCR profiles, the reactions were carried out in triplicate.
RAPD-PCR profiles analysis
The DNA banding patterns of RAPD-PCR assays were analyzed using TotalLab TL120 Trace version 2006 (Nonlinear Dynamics). The Dice coefficient of similarity was calculated on the base of presence/absence of bands while the unweighted pair group method with arithmetic averages (UPGMA) as agglomerative method. Cluster analysis has been carried out by using the statistical excel plugin XLstat 7.5 (Addinsoft, USA). To distinguish clonally different groups, the similarity percentage cut-off was set at 80%.
Whole genome sequencing analysis
Whole Genome Sequencing (WGS) analysis were performed on the S. aureus genomes, which were isolated from both AD patients and healthy carriers and sequenced on the MiSeq platform, using the Nextera XT library preparation kit v3 according to the manufacturer’s instructions (Illumina, San Diego, CA, USA).
The quality of the short reads was evaluated using FastQc v. 0.11.5 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/), combined with MultiQC v.1.11 (Ewels et al., 2016). Using the Shovill wrapper (https://github.com/tseemann/shovill), the reads were trimmed by Trimmomatic v. 0.39 (Bolger et al., 2014) and then de novo assembled by SPAdes assembler v.3.15.3 (Bankevich et al., 2012). The open reading frames were predicted and annotated by Prokka v.1.14.6 (Seemann, 2014), using the dedicated databases for the Staphylococcus genus associated with the software.
The average Nucleotide Identity (ANI) analysis and the Ribosomal Multilocus Sequence Typing (rMLST) were performed using fastANI v. 1.32 (Jain et al., 2018) and querying the PubMLST dedicated database (Jolley et al., 2012), respectively, in order to verify the taxonomy assignment of sequenced strains.
A kmer search was performed by AgrVate v.1.0.1 (Raghuram et al., 2022) on the agr-group database associated with the software to determine the agr locus type of the isolated genomes. The spa-type was established using DNAGear (AL-Tam et al., 2012). The MLST was performed using MLST software (https://github.com/tseemann/mlst) v. 2.19.0, scanning the assembled genomes against the PubMLST scheme for Staphylococcus aureus (Jolley et al., 2018). Only exact matches were kept. Alleles characterized by a sequence coverage of 100% and a sequence identity lower than 100% were further investigated mapping the paired-end short reads of the isolate on the most similar reference alleles, using Bowtie2 v.2.3.5.1 (Langmead and Salzberg, 2012). The Staphylococcal cassette chromosome mec (SCCmec) typing was performed by querying the extended database of SCCmecFinder (Kaya et al., 2018), to characterize the methicillin-susceptibility of the isolates.
To determine the presence of genes associated to antibiotic resistance and VFs, the sequences were aligned against the Comprehensive Antibiotic Resistance Database (CARD) (Alcock et al., 2019), ResFinder (Bortolaia et al., 2020) and Virulence Factor Database (VFDB) (Chen et al., 2016), using Abricate (https://github.com/tseemann/abricate) v. 0.9.9. The identified genes were then filtered, setting a threshold of 80% for both sequence identity and sequence coverage.
A core-genome based phylogenetic tree was constructed using the Bacterial Pan Genome Analysis (BPGA) software package based on unweighted pair group method with arithmetic mean (UPGMA) and setting 80% as sequence identity cut-off (Chaudhari et al., 2016). In addition to our isolates, S. aureus USA300 (WGS Project accession: JADMJC000000000) was also included in the analysis while the S. epidermidis strain ATCC12228 (NCBI accession: NC_004461) was considered as the outgroup.
To investigate the presence of the bsa (bacteriocin of S. aureus) operon and the dprA (DNA processing protein A) gene, the reference sequences of such loci were extracted from the publicly available genomes of the S. aureus Newman strain (Genbank accession: AP009351) and the S. aureus N315 (Genbank accession: BA000018) strain, respectively, and aligned against the newly assembled S. aureus genomes, using the stand-alone version of the Blast algorithm v.2.12.0 (https://blast.ncbi.nlm.nih.gov/Blast.cgi).
A k-mer search by RUCS (Thomsen et al., 2017) was performed to identify unique and core genomic sequences associated with a single group of assemblies.
R analysis
Data produced from the previous steps of the analysis were filtered, manipulated and tested for differences using R v.4.1.2 (https://www.R-project.org/) and its available libraries data.table (https://cran.r-project.org/package=data.table), reshape2 (Wickham, 2007) and dunn.test (https://cran.r-project.org/package=dunn.test). The results were then displayed using ggplot (https://ggplot2.tidyverse.org), ggtree (Yu et al., 2017), ggtext (https://cran.r-project.org/package=ggtext), ggpubr (https://cran.r-project.org/package=ggpubr) and gridExtra (https://cran.r-project.org/package=gridExtra).
Functions in the vegan package (https://cran.r-project.org/package=vegan) were applied on a binary matrix representing the gene clusters identified in each assembly, in order to compute a distance matrix and to perform a Principal Coordinate Analysis (PcoA); the Jaccard similarity index and the Lingoes correction for negative eigenvalues were applied to two steps of the analysis, respectively.
The Pan-genome composition was evaluated for each assembly and compared between the groups. A functional annotation was performed in order to assign a Cluster of Orthologous Groups (COG) category (https://www.ncbi.nlm.nih.gov/research/cog) to each gene cluster. The relative frequencies of the COG categories were computed for every component of the pan-genome.
Antibiotic susceptibility test
Antibiotic susceptibility profiles of S. aureus strains were defined by the guidelines of the European Committee on Antimicrobial Susceptibility Testing (EUCAST, v. 6.0) and performed by the automated microbial system Vitek-2 (BioMérieux), at the Laboratory of Microbiology Analysis of the Policlinico Umberto I, “Sapienza” University of Rome. The antibiotic susceptibility profile was defined by testing the following antibiotics: benzylpenicillin, oxacillin, gentamicin, levofloxacin, erythromycin, clindamycin, linezolid, daptomycin, teicoplanin, vancomycin, tetracycline, tigecycline, fusidic acid, mupirocin, rifampicin, trimethoprim/sulfamethoxazole.
Biofilm assay
S. aureus strains were tested for their ability to form biofilm as previously described (Stepanović et al., 2007). S. aureus ATCC 6538P was used as biofilm positive control.
Biofilms were quantified by measuring the absorbance at λ 570 nm (A570). The average optical density (OD) value was calculated and cut-off values (ODc) were established. ODc is defined as the average of OD of negative control + three standard deviations. According to their absorbance, isolates were defined as no biofilm producer OD ≤ ODc; weak biofilm producer ODc < OD ≤ (2 X ODc); moderate biofilm producer (2 X ODc) < OD ≤ (4 X ODc); strong biofilm producer (4ODc) < OD.
Cell line and culture conditions
The human immortalized keratinocyte cell line, HaCaT, was cultured in Dulbecco’s modified Eagle’s medium high glucose (DMEM, Corning, Italy) supplemented with 10% inactivated fetal bovine serum (FBS, Gibco, Italy), and 1% penicillin/streptomycin (Sigma-Aldrich, Italy) and maintained in 5% of CO2 atmosphere at 37°C.
Invasion assays
Invasion assays were performed in twenty-four- well culture plates seeded with HaCat cells (2x105 cells/ml) and incubated at 37°C in 5% CO2 for 48 h. Cell monolayers were infected with bacterial suspensions (multiplicity of infection [MOI] of 10) in 0.5 ml of cell culture medium devoid of antibiotics. Bacteria were centrifuged onto cellular monolayers two times, each at 1200 rpm for 2.5 min, and incubated for 1 h at 37°C. After infection time, cell monolayers were washed, and incubated with fresh medium containing 200 μg/ml of gentamicin (Sigma-Aldrich, Italy) for 1 h, to kill extracellular bacteria. Finally, cells were lysed with Triton X-100 (0.1% v/v) and plated on TSA plates to determine the number of intracellular vital bacteria (calculated as cfu/ml).
Strains were considered invasive when the percentage of the ratio between the number of intracellular bacteria and the initial inoculum was ≥0.1%. Data represent the mean of three independent experiments in duplicate.
Statistical tests
An initial descriptive analysis was performed with tables and plots corresponding to the type of qualitative or quantitative variables. Statistical analysis was carried out using both parametric and non-parametric tests according to the analyzed variable. Particularly, non-parametric univariate statistical tests, X2 and Fisher’s test, were used for the discrete variables. For continuous variables, data were tested for normality using the Shapiro-Wilk Normality test, followed by the pairwise comparisons using T-test or Mann–Whitney U test. Comparisons among multiple groups were tested using Kruskal–Wallis test followed by pairwise Dunn’s post hoc tests. Levene’s test was used to assess the equality of variances of the groups.
The permutational multiple analysis of variance (PERMANOVA) test was calculated on distance matrices to assess the presence of statistically significant partitions between groups, by performing 1000 permutations.
The statistical association between ordinal and/continuous variables was assessed by Kendall’s rank correlation analysis.
When necessary, p-values were corrected, by using the Benjamin–Hochberg false discovery rate (FDR) procedure to account for multiple hypothesis testing. In each case, a p-value ≤ 0.05 was considered statistically significant.
ResultsPatients’ characteristics and clinical response
Fifteen patients, were included in the study. At the time of the enrollment, 6 out of 15 patients were affected by a severe form of AD (SCORAD index≥75), while 9 out of 15 patients presented with a moderate form of AD (26 ≥ SCORAD index < 75). After the sampling, all the patients started a systemic treatment with dupilumab as per protocol (600 mg subcutaneous injection once, then 300 mg every 2 weeks). All the patients had a good response to treatment in an average time of 5.1 months with complete clinical remission of the disease in 5 out of 15 patients, while 9 out of 15 patients had a significant reduction of severity (SCORAD < 25). After the treatment, 3 of the patients with moderate AD at t0 showed a complete clinical remission, while 6 patients improved to mild AD. Concerning the severe cases of AD at t0, 2 patients showed a complete clinical remission while 4 improved to mild AD form. The only patient who didn’t show clinical remission was excluded from the study.
General colonization pattern
A total of 50 S. aureus isolates, of which 37 were from AD patients at t0 and 13 at t1, were obtained by microbiological analysis. S. aureus was detected at t0 i) both in the anterior nares and on the lesional skin in 80% (12/15) of the patients and ii) on all of the sampled anatomical sites in 40% (6/15), whereas iii) cutaneous colonization alone was found in 20% (3/15) of the patients.
At t1, 21.4% (3/14) of patients were nasally and skin colonized; 35.71% (5/14) showed S. aureus colonization only in the nose and 14.2% (2/14) only on the skin, while in four patients, S. aureus was completely eradicated.
To further examine the correlation between S. aureus skin colonization and AD severity, a quantitative detection of S. aureus was carried out. At t0, S. aureus was the predominant bacterial species on the lesional skin, whose bacterial load resulted in 104- 106 cfu/ml. Instead, a lower load of S. aureus (10 -102 cfu/ml) was detected at t1, on the skin of AD patients. A significant positive correlation was found between the S. aureus load and disease status (severe, moderate, or mild) (τb =0.837; p-value < 0.0001).
A single <italic>S. aureus</italic> strain dominates in each AD patient
A RAPD-PCR analysis was carried out to determine the genetic relatedness of 50 S. aureus isolates from the skin and nose of 15 non-related AD patients. By using 80% similarity, a total of 21 distinctive clonal groups (strains) were detected (data not shown), 17 from AD patients at t0 and 4 from AD patients at t1. Furthermore, in 33.3% (5/15) of the patients, the t0 S. aureus isolates from the nose, the lesional skin and non lesional skin of the same subject belonged to the same clonal group; in 33.3% (5/15) of the patients, S. aureus isolates from the nose and the lesional skin of the same subject belonged to the same clonal group; in 13.3% (2/15) of the patients, S. aureus isolates from the lesional and non-lesional skin of the same subject belonged to the same clonal group; in 6.6% (1/15), S. aureus was recovered only in the lesional skin; in 6.6% (1/15) of the patients, isolates from nose and lesional skin showed no genetic relatedness and in 6.6% (1/15), the nasal and the lesional skin isolates belonged to the same clonal group, which was different from the one of the non lesional skin isolate.
At t1 S. aureus was recovered only from the nose of 57.1% (5/14) patients and these isolates resulted clonal to those recovered at t0 from the same patient. In 14.2% (2/14) patients, S. aureus was recovered at t1 only from the skin and these isolates were non-clonal with those isolated at t0 from the same patient. Both nasal and skin colonization at t1 was found in 21.4% (3/14) patients and the nasal isolates resulted clonal to those recovered at t0. Only in one patient, the t1 skin isolate was clonal to those recovered at t0. Overall, the examination of the intra-host genetic heterogeneity of the colonizing S. aureus showed a clonal expansion and diffusion of a single strain that can persist in the post-flare.
Whole genome sequencing analysis
WGS analysis was performed on 38 S. aureus strains (21 from AD patients and 17 from healthy carriers). Genomic assemblies obtained from high-quality short reads and their features are summarized in Table S1. The typical %GC content and genome size of S. aureus organisms, as well as the taxonomy assessment analysis supported the taxonomy assignment of the strains.
<italic>In silico</italic> typing analysis
All S. aureus strains were typed according to MLST, spa-, agr- (Table S2), and SCCmec-typing.
<italic>spa</italic>-typing
The spa-typing analysis evidenced the presence of 27 different types within the population, of which 3 were never reported before. No predominant types were observed, while the most prevalent ones were t1451, t091, and t0148 which were found in 13.1% (5/38), 10.5% (4/38), and 7.9% (3/38) of the population, respectively.
Among the AD-associated strains, 13 different spa-types were identified, where t1451 and t091 were the two most prevalent types among the diseased subjects, each one characterizing 19% (4/21) of the strains. Likewise, both t359 and t008 were detected in 9.5% (2/21) of the AD-associated strains.
Concerning the healthy carrier strains, 14 different spa-types were identified and the most common types were t148 and t012, which were found in 17.6% (3/17) and 11.7% (2/17) of the strains, respectively. The latter spa-types (i.e. t4022, t153, t1994, t166, t331, t1078, t2777, t1509, t160, t449, t1451) were identified only once, composing up to 64.7% (11/17) of the group.
Three novel spa-types were identified. The AD strains 2pSAt0LS and 4pSAt0LS showed two novel spa-types, having similar repeat sequences to those of t449 and t148, respectively. A third novel spa-type was detected from the healthy carrier strain 29cSA, whose repeat sequence is similar to t2777. (Table S3)
No statistical differences were found in the distribution of spa-types between strains isolated from individuals with AD and those from healthy carriers.
MLST
Analogously to the spa-typing, the MLST evidenced a high ST diversity due to the identification of 24 distinct types, of which 4 were never reported before. The most frequent type was ST7, which was observed in 13.1% (5/38) of the population. Among the 21 AD S. aureus strains, 12 different STs were identified. The group-wise evaluation showed that ST7, ST398, ST3510, ST5, and ST8 composed 23.8% (5/21), 9.5% (2/21), 9.5% (2/21), 9.5% (2/21), and 9.5% (2/21) of the AD strains, respectively. In three cases, S. aureus strains deriving from the same individual showed different STs. Furthermore, two AD strains isolated from the same patient were characterized by a novel arcC allele, whose sequence differs from allele 3 for a single C-T transition at position 303.
Among the S. aureus strains isolated from healthy carriers, 12 different STs were identified. In particular, ST72, ST30, ST34, and ST45 were observed in 17.6% (3/17), 11.7% (2/17), 11.7% (2/17), and 11.7% (2/17) of the strains, respectively. Two new alleles of aroE and arcC genes, respectively, were identified within the healthy carrier group. In particular, the first allele differed from aroE allele 5 with respect to a C-A transversion at position 420, while the second one differed from arcC 1 for a T-C transition at position 402. No statistically significant differences were observed in the distribution of STs between the two groups of strains.
<italic>agr</italic>-typing
agr-I represented the most frequent allotype, which was found in 63.1% (24/38) of the strains. A lower prevalence was observed for the agr-II and agr-III, which composed 21% (8/38) and 13.1% (5/38) of the strains, respectively, while only one healthy carrier strain was classified as agr-IV. No statistically significant differences were found in the distribution of agr-types between strains from AD and healthy carriers.
SCC<italic>mec</italic>-typing
According to the SCCmec-typing, the population is exclusively characterized by methicillin-susceptible S. aureus strains.
Virulence factor genes
In the present study, the presence of VF genes was determined by aligning the genomes of our strains with the sequences present in the dedicated database VFDB.
Among the 100 VF genes identified in the whole population, 80 were found in ≥ 50% of the strains, as shown in Figures S1, S2. Virtually all VF genes associated with immune evasion, adhesion, and toxic activity (exoenzymes and toxins) were found in 100% of the population. By comparing the distribution of the VF coding sequences between the groups, no statistically significant differences have been evidenced, except for some enterotoxin genes (i.e. sec, sei, sell, sem, sen, seo, and seu) (Figure S3), that resulted significantly more prevalent in the healthy carriers group.
Antimicrobial resistance genes identification
All the antimicrobial resistance (AMR) genes in the reference databases were tested for their presence and we noticed that each strain was characterized by at least arlR, arlS, norA, mgrA, tet38, LmrS, mepA, mepR, cls, pgsA, mprF, walK, rpoB, fusA, gyrA, parC and parE. NorA, Tet38, and LmrS efflux pumps are members of the major facilitator superfamily (MFS). NorA expression is regulated by the two-component system ArlRS and MgrA, and mediates the resistance to fluoroquinolones (e.g. ciprofloxacin). MgrA also regulates the expression of tet38 gene (Truong-Bolduc et al., 2005), whose overexpression confers resistance to tetracycline (Chen and Hooper, 2018). Instead, the LmrS multidrug efflux pump is involved in the resistance mechanisms to aminoglycosides, macrolides, phenicols, diaminopyrimidine, and oxazolidinone (Floyd et al., 2010). Previous studies suggest that mutations on specific genes can mediate the resistance to tigecycline (mepA and mepR) (McAleese et al., 2005; Dabul et al., 2018; Fang et al., 2020), to daptomycin (cls, pgsA, mprF, and walK) (Peleg et al., 2012; Song et al., 2013), to rifampicin (rpoB) (Wichelhaus et al., 1999), to fusidic acid (fusA) (Lannergård et al., 2009), to trimethoprim/sulfamethoxazole (dfrA and folP) (Nurjadi et al., 2020) and to levofloxacin (gyrA, parC and parE)(Sanfilippo et al., 2011).
As shown in Figure S4, further genes were identified in some of our strains, although no statistically significant differences were detected between the groups. The aac(6’)-Ie-aph(2’’)-Ia and ant(4’)-Ib genes mediate the aminoglycoside resistance and they were detected among AD-associate strains in 4.76% (1/21) and in 9.52% (2/21), respectively, whereas only ant(4’)-Ib was identified among healthy carriers in 5.88% (1/17) of the strains.
qacA is another gene encoding for a MFS efflux pump and it confers resistance to several fluoroquinolones (Alam et al., 2003). This gene was found only in 5.88% (1/17) of AD-associated strains. The analysis revealed also the presence of AMR genes involved in the resistance to erythromycin (ermC, ermT, mphC and msrA), to clindamycin (ermC and ermT), to benzylpenicillin (blaZ), to fosfomycin (fosB) and to fusidic acid (fusB). Among the AD-associated strains, ermT was found in 19% (4/21), mphC in 9.52% (2/21), msrA in 9.5% (2/21), blaZ in 47.6% (10/21) and fosB in 42.8% (9/21). Instead, among the strains isolated from healthy carriers we detected ermC and ermT in 11.76% (2/17) and 5.88% (1/17), respectively. fusB was found in 5.88% (1/17), msrA in 5.88% (1/17), mupA in 5.88% (1/17), blaZ in 82.35% (14/17) and fosB in 58.8% (10/17).
A wide antibiotic susceptibility characterizes all <italic>S. aureus</italic> strains
All 38 S. aureus strains underwent automated susceptibility testing with the Vitek-2 system using EUCAST breakpoints. As shown in Table S4, most of the strains were susceptible to a wide range of antibiotics, regardless of isolation origin. In particular, none of the S. aureus strains was resistant to gentamicin, linezolid, daptomycin, tetracycline, tigecycline, and rifampicin. Furthermore, no methicillin-resistant S. aureus strains were detected in both groups. Such finding agrees with the results of the in silico SCCmec typing. Resistance to benzylpenicillin is the most common AMR in the population, characterizing 68.42% (26/38) of strains, and it was more diffused among healthy carrier strains than AD ones (82% vs 57%). Both groups were also characterized by erythromycin-resistant strains, accounting for 26.31% (10/38) of the population, as well as clindamycin-resistant strains, accounting for 21.05% (8/38). Resistances to oxacillin and to teicoplanin were observed only among AD strains, in 10% (2/21) and in 5% (1/21) of strains, respectively. No statistically significant differences were found between the groups for any of the antibiotics.
It is worth noticing that multi-drug resistance was observed exclusively in the healthy carrier strain 12cSA (resistant to erythromycin, benzylpenicillin, mupirocin, clindamycin, and fusidic acid) and in the AD strain 14pSAt1S (resistant to erythromycin, benzylpenicillin, oxacillin, and trimethoprim/sulfamethoxazole).
An overall agreement between the WGS-predicted resistance and the automated phenotypic susceptibility testing was noticed in the bacterial population (Table S5). All strains were susceptible to gentamicin, linezolid, and vancomycin, whose corresponding AMR genes were not detected in the population. Instead, both the mupirocin- and fusidic acid-resistance genes were exclusively found in the S. aureus 12cSA strain, which was the only resistant strain to these antibiotics.
Concordance between the presence of at least one AMR gene associated with a specific antibiotic resistance mechanism and the corresponding resistant phenotype was observed for erythromycin, clindamycin, benzylpenicillin, oxacillin, and for teicoplanin, with some exceptions. For instance, i) 7pSAt0LS is the only benzylpenicillin susceptible-strain carrying the blaZ gene. ii) No corresponding AMR genes were found in some strains showing an intermediate or resistant phenotype to clindamycin [30% (3/10)], benzylpenicillin [7,70% (2/26)], erythromycin [9.09% (1/11)], teicoplanin [100% (1/1)] and oxacillin [100% (2/2)). iii] The automated phenotypic susceptibility testing revealed that all the strains were tetracycline-susceptible, although tet38 and its regulator mgrA were ubiquitously found in the population.
In this study, we could not assess the concordance between genotypic and phenotypic profiles with respect to those antibiotics whose resistance is associated with a mutation on the specific genes (levofloxacin, rifampicin, tigecycline, daptomycin, and trimethoprim/sulfamethoxazole). Further studies will be conducted to investigate the presence of these genetic mutations.
Core-genome phylogeny
A phylogenetic reconstruction has been performed in order to investigate the genetic relationship among strains and to reveal a potential association with the presence of the disease. As shown in Figure 1, the population was split into 11 different clusters, each one characterized by strains sharing a similar spa, agr, and ST types as well as a similar set of VF and AMR genes. Furthermore, we observed that the strains deriving from the same patient tended to cluster together, although this was not the case for 4pSAt0LS and 4pSAt1S, 9pSAt0LS and 9pSAt0NLS, and 2pSAt0LS and 2pSAt0N. No statistically significant associations were observed between the identified clusters and the presence of the disease.
Core phylogenetic tree of the 38 S. aureus strains. The S. aureus USA300 strain was also included in the analysis, while the S. epidermidis strain ATCC12228, has been considered as the outgroup. Branch lengths (-log10 scale) expressed on the tree are proportional to the phylogenetic distances. Different colors were used to highlight the 11 clusters.
Pangenomic analysis
The analysis of the pan-genome revealed 5939 groups of orthologous gene clusters, of which 2086 (52.69%) were present in all strains within the analyzed population (core-genome), 1468 (37.08%) were found in more than one strain but not in all (accessory-genome) while only a more limited number of OGC, 405 (10.23%), were identified in a single strain (unique-genome).
The sequential addition of S. aureus strains to the analysis highlighted a rapid decrease in size of the core-genome, as well as a continuous growth of the pan-genome. These trends revealed an “open” nature of the S. aureus pan-genome in agreement with the power curve, computed by the BPGA pipeline (Figure 2).
Pan-core plot depicting the trend of the pan-genome and the core-genome by the sequential addition of S. aureus strains to the analysis.
As shown in Figure 3A, the core-genome constituted the largest portion within the genome of each strain (median, IQR; 83.71%, 82.82%-84.24%), followed by the accessory genome (median, IQR; 16.14%, 15.13%-16.67%) and then by the unique genome (median, IQR; 0.16%, 0.05%-0.39%). No statistically significant differences were detected in the relative abundances of each component of the pan-genome, between the groups of strains isolated from the lesional skin of AD patients and those obtained from healthy subjects (Figure 3B). However, the healthy carrier strains were characterized by a significantly higher variability in the relative abundance of both core and accessory genes with respect to the strains coming from AD subjects (p-valuecore = 0.013; p-valueaccessory = 0.014).
(A) Genomic composition of each S. aureus strain. (B) Relative abundance of the core, accessory, and unique genes. (C) Within-group diversity, measured by the Jaccard dissimilarity index. Values can range from 0 (i.e. the organisms share the same set of genes) to 1 (i.e. no genes in common). Comparisons between AD lesional strains and healthy carriers-related strains were carried out by Mann-Whitney U-test, considering a p-value ≤ 0.05 as statistically significant. (****p-value ≤ 0.0001). (D) Between-diversity analysis expressed by the Principal coordinate analysis, based on the Jaccard index and performed on AD lesional skin and healthy carriers-related strains. For each Principal coordinate, the variance explained is expressed by numbers within parentheses.
The within-group genomic diversity, assessed on the basis of the Jaccard dissimilarity index (thus on the number of genes shared by each pair of strains), evidenced a higher within-group variability in the population of S. aureus strains isolated from healthy carriers than that obtained from AD (p-value < 0.0001) (Figure 3C). No significant partitions have been determined between the groups as evidenced by results obtained from principal coordinate analysis (PCoA), followed by the PERMANOVA test (Figure 3D). Furthermore, no genomic sequences have been exclusively found in either the AD-related or the healthy carriers-related strains.
<italic>bsa</italic> operon and <italic>dprA</italic> gene search
The bsa operon encodes for a lantibiotic bacteriocin and it was detected in only two strains from each group. The aligned genes were characterized by sequence coverage >98% and sequence identity >85% with respect to reference sequences (Table S6).
On the other hand, virtually all strains (34/38) carried the dprA gene, known to affect the ability to assimilate foreign DNA (sequence coverage >99% and sequence identity >98%).
Functional enrichment analysis
A functional enrichment analysis was performed, allowing the classification of the identified gene families according to the 20 COG categories. As shown in Figure 4, the most prevalent categories (>33%) in the core-genome were related to the Transport and metabolism of amino acids (E), nucleotides (F), carbohydrates (G), coenzymes (H), lipids (I), and inorganic ions (P). Transcription (K) and Replication, recombination and repair (L) mechanisms were the most enriched COG categories in the accessory- (13.96% and 17.49%, respectively) and in the unique-genome (23.73% and 18.11%, respectively). (Figure 4) General function prediction only (R) represents the top enriched COG category in both AD- (12.45%; 0.18%) and healthy carriers (12.28%; 0.29%) groups (Figure S5).
Functional distribution of the core, accessory and unique genes of the population, according to COG categories.
Interestingly, genes related to post-translational modification, protein turnover and chaperones (O), and intracellular trafficking, secretion, and vesicular transport (U) resulted in statistically more abundant among the AD lesional strains (p-valueO = 0.0039 and p-valueU = 0.043) (Figure S5).
AD strains are strong biofilm producers
Biofilm is an important factor for S. aureus skin colonization, and most of the isolates from AD lesional skin resulted in strong biofilm producers (Geoghegan et al., 2018). Accordingly, our results revealed that 67% (14/21) of AD strains were strong biofilm producers, 33% (7/21) were moderate biofilm producers, and no weak biofilm producers were found. In contrast, healthy subjects did not carry strong biofilm-producing strains, while moderate and weak biofilm producers accounted for 47% (8/17) and 53% (9/17) of the group, respectively. Statistically significant differences in the percentage of strong and weak biofilm producers were found between AD and healthy carrier strains (p valuestrong ≤ 0.0001 and p valueweak≤ 0.001).
Invasion ability could be not essential for AD strains
S. aureus is able to invade and survive in different host cells, including keratinocytes. To evaluate the invasive behavior of S. aureus strains from AD patients and healthy carriers, gentamicin protection assays on HaCat cells were performed. A significantly higher percentage of invasive S. aureus strains was observed among healthy carriers (88%) with respect to that recovered in AD patients (43%) (p-value ≤ 0.01).
Discussion
AD is one of the most frequent chronic inflammatory skin diseases characterized by recurrent eczematous lesions, itch and dry skin (David Boothe et al., 2017). The complex and multifactorial pathophysiology of AD involves a strong genetic predisposition, epidermal dysfunction, T-cell driven inflammation and skin microbiome abnormalities. In particular, S. aureus represents a dominant colonizer and pathogen, especially of the AD lesional sites, and it contributes to AD pathogenesis in many ways (e.g. skin barrier disruption and direct proinflammatory effects) (Geoghegan et al., 2018). However, the exact roles played by this microorganism in vivo are not completely understood.
Herein, we performed a genotypic and phenotypic characterization of patient-derived strains, sampled at two different stages of the disease, namely at acute phase before starting dupixent® (dupilumab) therapy (t0) and at post-acute phase after reaching clinical remission with treatment (t1).
S. aureus was detected at t0 i) both in the anterior nares and on the lesional skin in 80% (12/15) of the patients and ii) on all of the sampled anatomical site in 40% (6/15), whereas iii) a cutaneous colonization alone was found in 20% (3/15) of the patients. In agreement with previous reports (Tauber et al., 2016; Totté et al., 2016), our results confirmed S. aureus as the predominant bacterial species on the inflamed skin of AD patients.
Following the dupilumab therapy, about 57% of patients still presented S. aureus colonization on the nose and 35.7% on the skin sites, while a complete eradication of S. aureus from all of the sampled body sites was observed in 28.57% of cases.
A significant positive correlation was found between the load of S. aureus on lesional skin and the severity of AD, defined by the objective SCORAD. S. aureus was more abundant in patients with moderate or severe disease (SCORAD ≥ 25). After the therapy, the significant improvement of the disease activity was accompanied by the simultaneous eradication of both nasal and skin colonizations, in 28% of AD patients. However, more than half of the AD patients were characterized by persisting S. aureus nasal colonization, confirming the anterior nares as the primary reservoir for S. aureus and as a major source for extra-nasal auto-transmission (Patel, 2001).
RAPD-PCR analysis on the AD-associated S. aureus isolates showed that in 80% of patients the nose, the lesional and the non-lesional skin, were colonized by the same strain during the flare, suggesting an intrapersonal spread of the same clone from one body site to another. The nasal isolates recovered at t1 resulted clonal to those sampled at t0 in 57% of patients, while the skin isolates were different from those sampled at t0 in 28% of patients, confirming the effectiveness of the treatment with dupilumab to reduce the S. aureus-induced inflammation.
Based on the intra-host genetic heterogeneity of S. aureus, we can also speculate about the presence of patient-specific strains that persist in the nose in the post-flare, constituting a potential reservoir for skin reinfections (Chiu et al., 2009). These findings should be carefully evaluated since they imply serious limitations in the long term efficacy of monoclonal antibodies-based therapeutic strategies
Specific therapies, targeting also the S. aureus nasal colonization, should be taken into consideration in this context.
Differently from the healthy controls S. aureus strains, the ones isolated from AD skin have been associated with higher levels of skin inflammation (Byrd et al., 2017). Such findings evidence that the contribution of S. aureus to the complexity of the disease may be strain-dependent. Herein, we performed an in silico genotyping analysis on the S. aureus strains, evaluating their spa-types, agr-types and STs, to establish whether specific clones of S. aureus prevail among AD patients. The high degree of genetic heterogeneity in terms of spa-types and STs of our S. aureus population evidenced the absence of a prevailing genotype linked to AD. Analogously, the large genetic variability of S. aureus strains isolated from AD patients has been reported by previous studies (Capoluongo et al., 2001; Chung et al., 2008; Kim et al., 2009). S. aureus possesses a vast arsenal of VFs that facilitate its invasion into and its survival within the human host. The expression of these VFs is controlled by the agr-system (Bernabè et al., 2021). The agr operon in S. aureus includes agrA, agrB, agrC, and agrD genes, and it could be divided into four types (namely agr-I to -IV), according to the sequences of agrC and agrD genes. Analogously to previous works (Zhao et al., 2016; Nakamura et al., 2020; Tayebi et al., 2020), our findings highlighted that agr-I is the most predominant type among S. aureus strains, suggesting that agr-type is not associated with AD. As Nakamura and colleagues point out (Nakamura et al., 2020), the expression of a functional agr system is required for epidermal colonization and the induction of AD-like inflammation, suggesting a pivotal role of the agr virulence in the development of AD in humans. Thus, the attention should be focused on the retention of a functional agr system rather than on the sequence type of agrC and agrD.
VF produced by S. aureus have a crucial role in the inflammation process and in skin barrier dysfunction in AD (Taskapan and Kumar, 2000; Tomi et al., 2005; Syed et al., 2015; Geoghegan et al., 2018). Sequence alignment between our strains and the reference database evidenced that most of the detected VF genes were widely spread in the population (prevalence ≥ 50%). In particular, all VFs involved in the immune evasion, adhesion, and toxic activity (exoenzymes and toxins) were found in the whole population, suggesting a common pathogenic potentiality that is independent from pathophysiological context. Aziz and colleagues suggest that specific enterotoxins may play a role in AD (Aziz et al., 2020). Instead, we found that genes encoding for enterotoxins were significantly more prevalent among healthy carrier strains. Taken together, these evidences may suggest that the role played by enterotoxins is not fundamental to S. aureus in AD. Overall, our results are consistent with previous studies, in which no CCs, agr-types, and VFs have been consistently linked to the pathogenesis of infection or colonization in AD (Benito et al., 2016; Fleury et al., 2017; Geoghegan et al., 2018).
Both genotypic and phenotypic analysis were performed in our study in order to assess the antibiotic resistance profile of all the 38 S. aureus strains. The WGS analysis evidenced the presence of a discrete collection of AMR genes, most of which were ubiquitously found in the bacterial population. Regardless of isolation origin, a widespread susceptibility to several antibiotics was revealed by the phenotypic analysis. Only methicillin-susceptible S. aureus strains were detected in the population, according to both of the analyses. Taken together, these results highlighted an overall agreement between the WGS-predicted resistance and phenotypic susceptibility tests, demonstrating a high concordance to 10 different agents (gentamicin, fusidic acid, mupirocin, linezolid, vancomycin, erythromycin, clindamycin, benzylpenicillin, oxacillin and teicoplanin), with some exceptions. According to the phenotypic test, the S. aureus strain 7pSAt0LS was benzylpenicillin-susceptible, although it carried the blaZ gene. As previously reported (Ivanovic et al., 2023), we could hypothesize that incompleteness of the bla operon could be responsible for the susceptible phenotype of this strain. No corresponding AMR genes were found in some strains showing an intermediate or resistant phenotype to clindamycin, benzylpenicillin, erythromycin, teicoplanin and oxacillin. However, we may reasonably assume that these intermediate or resistant phenotypes are likely mediated by AMR genes not yet determined, due to the current limited repertoire of sequences available to date. Phenotypic tests revealed that all the strains were tetracycline-susceptible, although tet38 and its regulator mgrA were ubiquitously found in the population. This apparent incongruence may be explained by the study conducted by Ding and colleagues (Ding et al., 2008), that revealed a distinct expression profile of the S. aureus tet38 efflux pump between the in vitro growth model and the in vivo infection model in mice, hypothesizing that environmental triggers, other than antimicrobial substances, can play an important role in the overexpression of tet38. Successively, the authors (Chen and Hooper, 2018) evaluated whether the selective overexpression of tet38 has an effect on response of S. aureus to tetracycline in the in vitro and in vivo models, observing a reduced efficacy of the antibiotic only in the abscess (in vivo model).
In our study, we could not assess the concordance between genotypic and phenotypic profiles to several antibiotics (i.e. levofloxacin, rifampicin, tigecycline, daptomycin and trimethoprim/sulfamethoxazole), whose resistance is due to mutations on specific genes. Further studies will be conducted to investigate the presence of these genetic mutations. Although genomic analysis represents a promising tool for the determination of the resistance profile of microbial strains, our findings evidence that the phenotypic characterization of strains still plays a primary role in the evaluation of the antimicrobial susceptibility.
A considerable degree of variety within the S. aureus population was demonstrated by the results of the reconstruction of phylogenetic relationships among strains. It is likely that the diversity pattern does not reflect the natural division of the species into subgroups and, in this context, belonging to a certain spa or MLST type appears to be the most appropriate trait explaining the genetic clustering of strains. The availability of a number of putative non-mutually exclusive mechanisms that might contribute to the patterns seen in the data could be used to explain these findings. First, there might be considerable biological obstacles preventing genetic exchange between lineages of the S. aureus population; these obstacles include restriction-modification (RM) systems, which operate as a major roadblock to DNA transfer (Corvaglia et al., 2010). The genetic exchange may be promoted or prevented by several systems that restrict phage and plasmid infection, such as CRISPR-cas systems (Marraffini and Sontheimer, 2008), phage-control systems (Ram et al., 2012; Ram et al., 2014), and variety in cell surface phage receptors (Winstel et al., 2014; Li et al., 2015). Other mechanisms, including plasmid incompatibility, functional restraints, or even functional redundancy of genes, can strengthen these genetic barriers and restrict the maintenance of DNA sequences acquired by clones. The existence of barriers across lineages also increases the chance that within each of these subgroups, recombination events may be responsible for the observed genetic similarity among strains. Another idea is that the existing lineages may be the isolated offspring of several extinct failed lineages, thus the consequence of the formation of a few extremely successful groups.
Locally, it is frequently noted that only a small number of CCs or STs dominate at any given time, but this pattern appears to be broken by the invasion of new dominant clones (Blanc et al., 2007; Chambers and Deleo, 2009), demonstrating that different clones coexist and compete for the same niche(s) (Sakwinska et al., 2009). In the S. aureus literature, examples of epidemic strain invasion and replacement are still accumulating (Uhlemann et al., 2014).
The Pangenome analysis of the S. aureus population revealed that, in each strain, approximately half of the pangenome has been represented by core genes. Although we expect that this evidence could be influenced by the limited sample size of the analyzed population, it is to note that the obtained results underline a very high degree of shared genetic sequence among the S. aureus population. A possible explanation to these findings is represented by the main role that the horizontal gene transfer plays among S. aureus strains, together with high rates of mobile genetic elements exchange during their co-colonization (McCarthy et al., 2014). According to previous observations, the core-genome of the S.aureus population is impacted by core genome transfer (CGT) supported by broad- and fine-scale trends in homoplasy attributable to macro-recombination, most likely driven by mobile elements (Everitt et al., 2014). The high degree with which the sequences are acquired and preserved in the genome of S. aureus leads to the formation of a gene reservoir which allows the single strains to colonize and adapt to a large variety of ecological niches, as demonstrated by the majority of cases, where the same clone predominates in inflamed and non-inflamed tissue. In this context, however, it should be highlighted that exposure to inflammatory conditions, such as those present in the skin lesions of AD subjects, is associated with a significant reduction in the variability in the gene content of the S. aureus pangenome, in line with the results obtained from the present study. On one hand, this reduction indicates the inflammatory conditions act to optimize the genetic content of the strains through selective pressure. On the other hand, it allows us to speculate about the presence of a certain degree of functional redundancy within the genome of the species.
Functional enrichment analysis showed that genes related to post-translational modifications (PTMs), protein turnover and chaperones, intracellular trafficking, secretion and vesicular transport resulted statistically more abundant among the AD strains. Although few studies have been conducted so far, evidence suggests that PTM mechanisms in S. aureus may be involved in the pathogenic lifestyle, as well as, in other cellular processes. For instance, glycosylation is involved in the cell wall teichoic acid (WTA) biosynthesis pathway; WTAs play an important role in the bacterial physiology, resistance to antimicrobial molecules, host interaction, virulence and biofilm formation (Mistretta et al., 2019). Likewise, phosphorylation is a common PTM mechanism, that is involved in the central metabolic processes and in the regulation of the efflux-pump NorA, by acting on the MgrA protein; however, the role of phosphorylation in staphylococcal pathogenesis, as well as the role of other PTMs, needs to be further investigated, especially regarding AD (Ohlsen and Donat, 2010). Moreover, bacteria are able to manipulate the host immune response through the secretion of membrane vesicles, that play an essential role in the delivery of bacterial effectors promoting virulence, biofilm formation, signal transduction, cytotoxicity, resistance to antimicrobial skin fatty acids and in skin inflammatory disorders (Gurung et al., 2011; Alpdundar Bulut et al., 2020; Kengmo Tchoupa and Peschel, 2020; Staudenmaier et al., 2022). Several studies have demonstrated the importance of staphylococcal biofilms in the pathogenesis of AD (Allen et al., 2014; Sonesson et al., 2017; Di Domenico et al., 2018; Haasbroek et al., 2022). Most of our AD strains were strong biofilm producers, confirming the importance of biofilm in i) the increased resistance to the host immune responses, in ii) the reduced susceptibility to antimicrobials and in iii) inflammation, thus, promoting skin colonization and persistence. Previous works have shown that S. aureus can invade host cells and persist intracellularly in cell culture models, including keratinocytes (Lowy, 2000; Sendi and Proctor, 2009; Sayedyahossein et al., 2015). Our results revealed that only a limited number of AD strains were invasive if compared to healthy carriers, suggesting that the ability to invade could be not essential for AD strains. Overall results evidence that genotypic analysis fails to identify specific features associable to AD strains and partially this may be due to the relatively low number of strains analyzed. However, the phenotypic characterization highlights substantial differences between the two groups regarding biofilm production and invasiveness of strains.
Conclusion
Although phenotypic differences were found between AD and healthy carrier S. aureus strains, our results evidenced no predominant genotypic traits in AD strains, suggesting the absence of S. aureus clones shared among AD patients.
We could speculate that the functional role played by S. aureus in AD likely depends on other factors, such as a differential gene expression pattern on the AD skin sites or PTM mechanisms that are involved in the S. aureus pathogenic lifestyle.
Data availability statement
The data analyzed in this study is subject to the following licenses/restrictions: Original data will be available on request. Requests to access these datasets should be directed to francesca.brunetti@uniroma1.it. The data presented in the study are deposited in the DDBJ/ENA/GenBank repositories, accession numbers:JARZEJ000000000, JARZEK000000000, JARZEL000000000, JARZEM000000000, JARZEN000000000, JARZEO000000000, JARZEP000000000, JARZEQ000000000, JARZER000000000, JARZES000000000, JARZET000000000, JARZEU000000000, JARZEV000000000, JARZEW000000000, JARZEX000000000, JARZEY000000000, JARZEZ000000000, JARZFA000000000, JARZFB000000000, JARZFC000000000, JARZFD000000000, JARZFE000000000, JARZFF000000000, JARZFG000000000, JARZFH000000000, JARZFI000000000, JARZFJ000000000, JARZFK000000000, JARZFL000000000, JARZFM000000000, JARZFN000000000, JARZFO000000000, JARZFP000000000, JARZFQ000000000, JARZFR000000000, JARZFS000000000, JARZFT000000000, JARZFU000000000.
Ethics statement
The studies involving human participants were reviewed and approved by Ethics Committee of Sapienza University of Rome. The patients/participants provided their written informed consent to participate in this study.
Author contributions
MPC, ATP and SG conceived and supervised the study. SG enrolled the patients and collected the samples. ALC, LM and CL analyzed the samples and performed the experiments. GR performed bacterial identification and the antibiotic susceptibility tests. FB and MM performed bioinformatic and statistical analysis. ALC and FB wrote the draft manuscript. ALC, FB, MPC and MM wrote and reviewed the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by Ateneo Grant 2018, RG1181642C2C4B17. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2023.1194254/full#supplementary-material
ReferencesAlamM. M.KobayashiN.UeharaN.WatanabeN. (2003). Analysis on distribution and genomic diversity of high-level antiseptic resistance genes qacA and qacB in human clinical isolates of Staphylococcus aureus. Microb. Drug Resist.9, 109–121. doi: 10.1089/107662903765826697AlcockB. P.RaphenyaA. R.LauT. T. Y.TsangK. K.BouchardM.EdalatmandA.. (2019). CARD 2020: antibiotic resistome surveillance with the comprehensive antibiotic resistance database. Nucleic Acids Res.48, D517–D525. doi: 10.1093/nar/gkz935AllenH. B.VazeN. D.ChoiC.HailuT.TulbertB. H.CusackC. A.. (2014). The presence and impact of biofilm-producing staphylococci in atopic dermatitis. JAMA Dermatol.150, 260–265. doi: 10.1001/jamadermatol.2013.8627Alpdundar BulutE.Bayyurt KocabasB.YazarV.AykutG.GulerU.SalihB.. (2020). Human gut commensal membrane vesicles modulate inflammation by generating M2-like macrophages and myeloid-derived suppressor cells. J. Immunol.205, 2707–2718. doi: 10.4049/jimmunol.2000731AlsterholmM.StrömbeckL.LjungA.KaramiN.WidjestamJ.GillstedtM.. (2017). Variation in Staphylococcus aureus colonization in relation to disease severity in adults with atopic dermatitis during a five-month follow-up. Acta Derm. Venereol.97, 802–807. doi: 10.2340/00015555-2667AL-TamF.BrunelA. S.BouzinbiN.CorneP.BañulsA. L.ShahbazkiaH. R. (2012). DNAGear–a free software for spa type identification in Staphylococcus aureus. BMC Res. Notes5, 642. doi: 10.1186/1756-0500-5-642AltunbulakliC.ReigerM.NeumannA. U.Garzorz-StarkN.FlemingM.HuelpueschC.. (2018). Relations between epidermal barrier dysregulation and staphylococcus species–dominated microbiome dysbiosis in patients with atopic dermatitis. J. Allergy Clin. Immunol.142, 1643–1647.e12. doi: 10.1016/j.jaci.2018.07.005AzizF.HisatsuneJ.YuL.KajimuraJ.Sato’oY.OnoH. K.. (2020). Staphylococcus aureus isolated from skin from atopic-dermatitis patients produces staphylococcal enterotoxin y, which predominantly induces T-cell receptor vα-specific expansion of T cells. Infect. Immun.88, e00360-19. doi: 10.1128/IAI.00360-19BankevichA.NurkS.AntipovD.GurevichA. A.DvorkinM.KulikovA. S.. (2012). SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing. J. Comput. Biol.19, 455–477. doi: 10.1089/cmb.2012.0021BenitoD.AspirozC.GilaberteY.SanmartínR.Hernández-MartinÁ.AlonsoM.. (2016). Genetic lineages and antimicrobial resistance genotypes in Staphylococcus aureus from children with atopic dermatitis: detection of clonal complexes CC1, CC97 and CC398. J. Chemother.28, 359–366. doi: 10.1179/1973947815Y.0000000044BernabèG.Dal PraM.RoncaV.PaulettoA.MarzaroG.SaluzzoF.. (2021). A novel aza-derivative inhibits agr quorum sensing signaling and synergizes methicillin-resistant Staphylococcus aureus to clindamycin. Front. Microbiol.12. doi: 10.3389/fmicb.2021.610859BieberT. (2008). Atopic dermatitis. N. Engl. J. Med.358, 1483–1494. doi: 10.1056/NEJMra074081BjerreR. D.BandierJ.SkovL.EngstrandL.JohansenJ. D. (2017). The role of the skin microbiome in atopic dermatitis: a systematic review. Br. J. Dermatol.177, 1272–1278. doi: 10.1111/bjd.15390BlancD. S.PetignatC.WengerA.KuhnG.ValletY.FracheboudD.. (2007). Changing molecular epidemiology of methicillin-resistant Staphylococcus aureus in a small geographic area over an eight-year period. J. Clin. Microbiol.45, 3729–3736. doi: 10.1128/JCM.00511-07BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a flexible trimmer for illumina sequence data. Bioinformatics30, 2114–2120. doi: 10.1093/bioinformatics/btu170BortolaiaV.KaasR. S.RuppeE.RobertsM. C.SchwarzS.CattoirV.. (2020). ResFinder 4.0 for predictions of phenotypes from genotypes. J. Antimicrob. Chemother.75, 3491–3500. doi: 10.1093/jac/dkaa345BrandtE. B.SivaprasadU. (2011). Th2 cytokines and atopic dermatitis. J. Clin. Cell. Immunol.2, 110. doi: 10.4172/2155-9899.1000110BrandweinM.SteinbergD.MeshnerS. (2016). Microbial biofilms and the human skin microbiome. NPJ Biofilms Microbiomes2, 3. doi: 10.1038/s41522-016-0004-zBreuerK.WittmannM.KempeK.KappA.MaiU.Dittrich-BreiholzO.. (2005). Alpha-toxin is produced by skin colonizing Staphylococcus aureus and induces a T helper type 1 response in atopic dermatitis. Clin. Exp. Allergy35, 1088–1095. doi: 10.1111/j.1365-2222.2005.02295.xByrdA. L.DemingC.CassadyS. K. B.HarrisonO. J.NgW. I.ConlanS.. (2017). Staphylococcus aureus and Staphylococcus epidermidis strain diversity underlying pediatric atopic dermatitis. Sci. Transl. Med.9, eaal4651. doi: 10.1126/scitranslmed.aal4651CallewaertC.NakatsujiT.KnightR.KosciolekT.VrbanacA.KotolP.. (2020). IL-4Rα blockade by dupilumab decreases Staphylococcus aureus colonization and increases microbial diversity in atopic dermatitis. J. Invest. Dermatol.140, 191–202.e7. doi: 10.1016/j.jid.2019.05.024CapoluongoE.GiglioA. A.LavieriM. M.Lesnoni-La ParolaI.FerraroC.CristaudoA.. (2001). Genotypic and phenotypic characterization of Staphylococcus aureus strains isolated in subjects with atopic dermatitis. higher prevalence of exfoliative b toxin production in lesional strains and correlation between the markers of disease intensity and colonization density. J. Dermatol. Sci.26, 145–155. doi: 10.1016/s0923-1811(00)00171-7CarrW. W. (2013). Topical calcineurin inhibitors for atopic dermatitis: review and treatment recommendations. Paediatr. Drugs15, 303–310. doi: 10.1007/s40272-013-0013-9ChambersH. F.DeleoF. R. (2009). Waves of resistance: Staphylococcus aureus in the antibiotic era. Nat. Rev. Microbiol.7, 629–641. doi: 10.1038/nrmicro2200ChaudhariN. M.GuptaV. K.DuttaC. (2016). BPGA- an ultra-fast pan-genome analysis pipeline. Sci. Rep.6, 24373. doi: 10.1038/srep24373ChenC.HooperD. C. (2018). Effect of Staphylococcus aureus Tet38 native efflux pump on in vivo response to tetracycline in a murine subcutaneous abscess model. J. Antimicrob. Chemother.73, 720–723. doi: 10.1093/jac/dkx432ChenL.ZhengD.LiuB.YangJ.JinQ. (2016). VFDB 2016: hierarchical and refined dataset for big data analysis–10 years on. Nucleic Acids Res.44, D694–D697. doi: 10.1093/nar/gkv1239ChiuL. S.HoM. S.HsuL. Y.TangM. B. (2009). Prevalence and molecular characteristics of Staphylococcus aureus isolates colonizing patients with atopic dermatitis and their close contacts in Singapore. Br. J. Dermatol.160, 965–971. doi: 10.1111/j.1365-2133.2009.09038.xChoS. H.StricklandI.BoguniewiczM.LeungD. Y. (2001). Fibronectin and fibrinogen contribute to the enhanced binding of Staphylococcus aureus to atopic skin. J. Allergy Clin. Immunol.108, 269–274. doi: 10.1067/mai.2001.117455ChungH.-J.JeonH. S.SungH.KimM. N.HongS. J. (2008). Epidemiological characteristics of methicillin-resistant Staphylococcus aureus isolates from children with eczematous atopic dermatitis lesions. J. Clin. Microbiol.46, 991–995. doi: 10.1128/JCM.00698-07ClarkeS. R.MohamedR.BianL.RouthA. F.Kokai-KunJ. F.MondJ. J.. (2007). The Staphylococcus aureus surface protein IsdA mediates resistance to innate defenses of human skin. Cell Host Microbe1, 199–212. doi: 10.1016/j.chom.2007.04.005ClausenM. L.EdslevS. M.AndersenP. S.ClemmensenK.KrogfeltK. A.AgnerT. (2017). Staphylococcus aureus colonization in atopic eczema and its association with filaggrin gene mutations. Br. J. Dermatol.177, 1394–1400. doi: 10.1111/bjd.15470ClausenM. L.EdslevS. M.NørresletL. B.SørensenJ. A.AndersenP. S.AgnerT. (2019). Temporal variation of Staphylococcus aureus clonal complexes in atopic dermatitis: a follow-up study. Br. J. Dermatol.180, 181–186. doi: 10.1111/bjd.17033CorvagliaA. R.FrançoisP.HernandezD.PerronK.LinderP.SchrenzelJ. (2010). A type III-like restriction endonuclease functions as a major barrier to horizontal gene transfer in clinical Staphylococcus aureus strains. Proc. Natl. Acad. Sci. U. S. A.107, 11954–11958. doi: 10.1073/pnas.1000489107DabulA. N. G.Avaca-CruscaJ. S.Van TyneD.GilmoreM. S.CamargoI. L. B. C. (2018). Resistance in In vitro selected tigecycline-resistant methicillin-resistant Staphylococcus aureus sequence type 5 is driven by mutations in mepR and mepA genes. Microb. Drug Resist.24, 519–526. doi: 10.1089/mdr.2017.0279David BootheW.TarboxJ. A.TarboxM. B. (2017). Atopic dermatitis: pathophysiology. Adv. Exp. Med. Biol.1027, 21–37. doi: 10.1007/978-3-319-64804-0_3Di DomenicoE. G.CavalloI.BordignonV.PrignanoG.SperdutiI.GurtnerA.. (2018). Inflammatory cytokines and biofilm production sustain Staphylococcus aureus outgrowth and persistence: a pivotal interplay in the pathogenesis of atopic dermatitis. Sci. Rep.8, 9573. doi: 10.1038/s41598-018-27421-1DingY.OnoderaY.LeeJ. C.HooperD. C. (2008). NorB, an efflux pump in Staphylococcus aureus strain MW2, contributes to bacterial fitness in abscesses. J. Bacteriol.190, 7123–7129. doi: 10.1128/jb.00655-08EverittR. G.DidelotX.BattyE. M.MillerR. R.KnoxK.YoungB. C.. (2014). Mobile elements drive recombination hotspots in the core genome of Staphylococcus aureus. Nat. Commun.5, 3956. doi: 10.1038/ncomms4956EwelsP.MagnussonM.LundinS.KällerM. (2016). MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics32, 3047–3048. doi: 10.1093/bioinformatics/btw354FangR.SunY.DaiW.ZhengX.TianX.ZhangX.. (2020). Mutations in the MepRAB efflux system contribute to the in vitro development of tigecycline resistance in Staphylococcus aureus. J. Glob. Antimicrob. Resist.22, 631–636. doi: 10.1016/j.jgar.2020.06.005FleuryO. M.McAleerM. A.FeuillieC.Formosa-DagueC.SansevereE.BennettD. E.. (2017). Clumping factor b promotes adherence of Staphylococcus aureus to corneocytes in atopic dermatitis. Infect. Immun.85, e00994-16. doi: 10.1128/IAI.00994-16FloydJ. L.SmithK. P.KumarS. H.FloydJ. T.VarelaM. F. (2010). LmrS is a multidrug efflux pump of the major facilitator superfamily from Staphylococcus aureus. Antimicrob. Agents Chemother.54, 5406–5412. doi: 10.1128/aac.00580-10GeogheganJ. A.IrvineA. D.FosterT. J. (2018). Staphylococcus aureus and atopic dermatitis: a complex and evolving relationship. Trends Microbiol.26, 484–497. doi: 10.1016/j.tim.2017.11.008GómezM. I.LeeA.ReddyB.MuirA.SoongG.PittA.. (2004). Staphylococcus aureus protein a induces airway epithelial inflammatory responses by activating TNFR1. Nat. Med.10, 842–848. doi: 10.1371/journal.pone.0027958GurungM.MoonD. C.ChoiC. W.LeeJ. H.BaeY. C.KimJ.. (2011). Staphylococcus aureus produces membrane-derived vesicles that induce host cell death. PLoS One, 6, p.e27958. doi: 10.1371/journal.pone.0027958HaasbroekK.YagiM.YoneiY. (2022). Staphylococcus aureus biofilm inhibiting activity of advanced glycation endproduct crosslink breaking and glycation inhibiting compounds. Antibiotics (Basel)11, 1412. doi: 10.3390/antibiotics11101412HadiH. A.TarmiziA. I.KhalidK. A.GajdácsM.AslamA.JamshedS. (2021). The epidemiology and global burden of atopic dermatitis: a narrative review. Life (Basel)11, 936. doi: 10.3390/life11090936HirasawaY.TakaiT.NakamuraT.MitsuishiK.GunawanH.SutoH.. (2010). Staphylococcus aureus extracellular protease causes epidermal barrier dysfunction. J. Invest. Dermatol.130, 614–617. doi: 10.1038/jid.2009.257HomeyB.SteinhoffM.RuzickaT.LeungD. Y. (2006). Cytokines and chemokines orchestrate atopic skin inflammation. J. Allergy Clin. Immunol.118, 178–189. doi: 10.1016/j.jaci.2006.03.047HornJ.StelznerK.RudelT.FraunholzM. (2018). Inside job: Staphylococcus aureus host-pathogen interactions. Int. J. Med. Microbiol.308, 607–624. doi: 10.1016/j.ijmm.2017.11.009HowellM. D.KimB. E.GaoP.GrantA. V.BoguniewiczM.DeBenedettoA.. (2009). Cytokine modulation of atopic dermatitis filaggrin skin expression. J. Allergy Clin. Immunol.124, R7–R12. doi: 10.1016/j.jaci.2009.07.012HuitemaL.PhillipsT.AlexeevV.Tomic-CanicM.PastarI.IgouchevaO. (2021). Intracellular escape strategies of Staphylococcus aureus in persistent cutaneous infections. Exp. Dermatol.30, 1428–1439. doi: 10.1111/exd.14235IvanovicI.BossR.RomanòA.GuédonE.Le-LoirY.LuiniM.. (2023). Penicillin resistance in bovine Staphylococcus aureus: genomic evaluation of the discrepancy between phenotypic and molecular test methods. J. Dairy Sci.106, 462–475. doi: 10.3168/jds.2022-22158JainC.Rodriguez-RL. M.PhillippyA. M.KonstantinidisK. T.AluruS. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun.9, 5114. doi: 10.1038/s41467-018-07641-9JevonM.GuoC.MaB.MordanN.NairS. P.HarrisM.. (1999). Mechanisms of internalization of Staphylococcus aureus by cultured human osteoblasts. Infect. Immun.67, 2677–2681. doi: 10.1128/IAI.67.5.2677-2681.1999JolleyK. A.BlissC. M.BennettJ. S.BratcherH. B.BrehonyC.CollesF. M.. (2012). Ribosomal multilocus sequence typing: universal characterization of bacteria from domain to strain. Microbiology158, 1005–1015. doi: 10.1099/mic.0.055459-0JolleyK. A.BrayJ. E.MaidenM. C. J. (2018). Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications. Wellcome Open Res.3, 124. doi: 10.12688/wellcomeopenres.14826.1KayaH.HasmanH.LarsenJ.SteggerM.JohannesenT. B.AllesøeR. L.. (2018). SCCmecFinder, a web-based tool for typing of staphylococcal cassette chromosome mec in Staphylococcus aureus using whole-genome sequence data. mSphere3, e00612-17. doi: 10.1128/mSphere.00612-17Kengmo TchoupaA.PeschelA. (2020). Staphylococcus aureus releases proinflammatory membrane vesicles to resist antimicrobial fatty acids. mSphere5, e00804-20. doi: 10.1128/mSphere.00804-20KimB. E.LeungD. Y.BoguniewiczM.HowellM. D. (2008). Loricrin and involucrin expression is down-regulated by Th2 cytokines through STAT-6. Clin. Immunol.126, 332–337. doi: 10.1016/j.clim.2007.11.006KimD. W.ParkJ. Y.ParkK. D.KimT. H.LeeW. J.LeeS. J.. (2009). Are there predominant strains and toxins of Staphylococcus aureus in atopic dermatitis patients? genotypic characterization and toxin determination of s. aureus isolated in adolescent and adult patients with atopic dermatitis. J. Dermatol.36, 75–81. doi: 10.1111/j.1346-8138.2009.00592.xKongH. H.OhJ.DemingC.ConlanS.GriceE. A.BeatsonM. A.. (2012). Temporal shifts in the skin microbiome associated with disease flares and treatment in children with atopic dermatitis. Genome Res.22, 850–859. doi: 10.1101/gr.131029.111LangmeadB.SalzbergS. L. (2012). Fast gapped-read alignment with bowtie 2. Nat. Methods9, 357–359. doi: 10.1038/nmeth.1923LannergårdJ.NorströmT.HughesD. (2009). Genetic determinants of resistance to fusidic acid among clinical bacteremia isolates of Staphylococcus aureus. Antimicrob. Agents Chemother.53, 2059–2065. doi: 10.1128/AAC.00871-08LeK. Y.OttoM. (2015). Quorum-sensing regulation in staphylococci-an overview. Front. Microbiol.6. doi: 10.3389/fmicb.2015.01174LeonardA. C.PetrieL. E.CoxG. (2019). Bacterial anti-adhesives: inhibition of Staphylococcus aureus nasal colonization. ACS Infect. Dis.5, 1668–1681. doi: 10.1021/acsinfecdis.9b00193LiX.GerlachD.DuX.LarsenJ.SteggerM.KühnerP.. (2015). An accessory wall teichoic acid glycosyltransferase protects Staphylococcus aureus from the lytic activity of podoviridae. Sci. Rep.5, 17219. doi: 10.1038/srep17219LoW. T.WangS. R.TsengM. H.HuangC. F.ChenS. J.WangC. C. (2010). Comparative molecular analysis of meticillin-resistant Staphylococcus aureus isolates from children with atopic dermatitis and healthy subjects in Taiwan. Br. J. Dermatol.162, 1110–1116. doi: 10.1111/j.1365-2133.2010.09679.xLowyF. D. (2000). Is Staphylococcus aureus an intracellular pathogen? Trends Microbiol.8, 341–343. doi: 10.1016/s0966-842x(00)01803-5MarraffiniL. A.SontheimerE. J. (2008). CRISPR interference limits horizontal gene transfer in staphylococci by targeting DNA. Science322, 1843–1845. doi: 10.1126/science.1165771McAleeseF.PetersenP.RuzinA.DunmanP. M.MurphyE.ProjanS. J.. (2005). A novel MATE family efflux pump contributes to the reduced susceptibility of laboratory-derived Staphylococcus aureus mutants to tigecycline. Antimicrob. Agents Chemother.49, 1865–1871. doi: 10.1128/AAC.49.5.1865-1871.2005McCarthyA. J.LoefflerA.WitneyA. A.GouldK. A.LloydD. H.LindsayJ. A. (2014). Extensive horizontal gene transfer during Staphylococcus aureus co-colonization in vivo. Genome Biol. Evol.6, 2697–2708. doi: 10.1093/gbe/evu214MeylanP.LangC.MermoudS.JohannsenA.NorrenbergS.HohlD.. (2017). Skin colonization by Staphylococcus aureus precedes the clinical diagnosis of atopic dermatitis in infancy. J. Invest. Dermatol.137, 2497–2504. doi: 10.1016/j.jid.2017.07.834MistrettaN.BrossaudM.TellesF.SanchezV.TalagaP.RokbiB. (2019). Glycosylation of Staphylococcus aureus cell wall teichoic acid is influenced by environmental conditions. Sci. Rep.9, 3212. doi: 10.1038/s41598-019-39929-1NaS. Y.RohJ. Y.KimJ. M.TamangM. D.LeeJ. R. (2012). Analysis of colonization and genotyping of the exotoxins of Staphylococcus aureus in patients with atopic dermatitis. Ann. Dermatol.24, 413–419. doi: 10.5021/ad.2012.24.4.413NakagawaS.MatsumotoM.KatayamaY.OgumaR.WakabayashiS.NygaardT.. (2017). Staphylococcus aureus virulent PSMα peptides induce keratinocyte alarmin release to orchestrate IL-17-dependent skin inflammation. Cell Host Microbe22, 667–677.e5. doi: 10.1016/j.chom.2017.10.008NakamuraY.TakahashiH.TakayaA.InoueY.KatayamaY.KusuyaY.. (2020). Staphylococcus agr virulence is critical for epidermal colonization and associates with atopic dermatitis development. Sci. Transl. Med.12, eaay4068. doi: 10.1126/scitranslmed.aay4068NguyenM. T.KraftB.YuW.DemicriogluD. D.HertleinT.BurianM.. (2015). The νSaα specific lipoprotein like cluster (lpl) of s. aureus USA300 contributes to immune stimulation and invasion in human cells. PloS Pathog.11 (6), e1004984. doi: 10.1371/journal.ppat.1004984NurjadiD.ZizmannE.ChanthalangsyQ.HeegK.BoutinS. (2020). Integrative analysis of whole genome sequencing and phenotypic resistance toward prediction of trimethoprim-sulfamethoxazole resistance in Staphylococcus aureus. Front. Microbiol.11. doi: 10.3389/fmicb.2020.607842OhlsenK.DonatS. (2010). The impact of serine/threonine phosphorylation in Staphylococcus aureus. Int. J. Med. Microbiol.300, 137–141. doi: 10.1016/j.ijmm.2009.08.016ParkK. D.PakS. C.ParkK. K. (2016). The pathogenetic effect of natural and bacterial toxins on atopic dermatitis. Toxins (Basel)9, 3. doi: 10.3390/toxins9010003PascoliniC.SinagraJ.PecettaS.BordignonV.De SantisA.CilliL.. (2011). Molecular and immunological characterization of Staphylococcus aureus in pediatric atopic dermatitis: implications for prophylaxis and clinical management. Clin. Dev. Immunol.2011, 718708. doi: 10.1155/2011/718708PatelR. (2001). Association between nasal methicillin-resistant Staphylococcus aureus carriage and infection in liver transplant recipients. Liver Transpl.7, 752–753. doi: 10.1053/jlts.2001.0070752PeacockS. J.de SilvaI.LowyF. D. (2001). What determines nasal carriage of Staphylococcus aureus? Trends Microbiol.9, 605–610. doi: 10.1016/s0966-842x(01)02254-5PelegA. Y.MiyakisS.WardD. V.EarlA. M.RubioA.CameronD. R.. (2012). Whole genome characterization of the mechanisms of daptomycin resistance in clinical and laboratory derived isolates of Staphylococcus aureus. PloS One7, e28316. doi: 10.1371/journal.pone.0028316PeschelA.OttoM. (2013). Phenol-soluble modulins and staphylococcal infection. Nat. Rev. Microbiol.11, 667–673. doi: 10.1038/nrmicro3110RaghuramV.AlexanderA. M.LooH. Q.PetitR. A.GoldbergJ. B.ReadT. D. (2022). Species-wide phylogenomics of the Staphylococcus aureus agr operon revealed convergent evolution of frameshift mutations. Microbiol. Spectr.10, e0133421. doi: 10.1128/spectrum.01334-21RamG.ChenJ.KumarK.RossH. F.UbedaC.DamleP. K.. (2012). Staphylococcal pathogenicity island interference with helper phage reproduction is a paradigm of molecular parasitism. Proc. Natl. Acad. Sci. U. S. A.109, 16300–16305. doi: 10.1073/pnas.1204615109RamG.ChenJ.RossH. F.NovickR. P. (2014). Precisely modulated pathogenicity island interference with late phage gene transcription. Proc. Natl. Acad. Sci. U. S. A.111, 14536–14541. doi: 10.1073/pnas.1406749111RojoA.AguinagaA.MoneckeS.YusteJ. R.GastaminzaG.EspañaA. (2014). Staphylococcus aureus genomic pattern and atopic dermatitis: may factors other than superantigens be involved? Eur. J. Clin. Microbiol. Infect. Dis.33, 651–658. doi: 10.1007/s10096-013-2000-zSakwinskaO.KuhnG.BalmelliC.FrancioliP.GiddeyM.PerretenV.. (2009). Genetic diversity and ecological success of Staphylococcus aureus strains colonizing humans. Appl. Environ. Microbiol.75, 175–183. doi: 10.1128/AEM.01860-08SanfilippoC. M.HesjeC. K.HaasW.MorrisT. W. (2011). Topoisomerase mutations that are associated with high-level resistance to earlier fluoroquinolones in Staphylococcus aureus have less effect on the antibacterial activity of besifloxacin. Chemotherapy57, 363–371. doi: 10.1159/000330858SayedyahosseinS.XuS. X.RudkouskayaA.McGavinM. J.McCormickJ. K.DagninoL. (2015). Staphylococcus aureus keratinocyte invasion is mediated by integrin-linked kinase and Rac1. FASEB J.29, 711–723. doi: 10.1096/fj.14-262774SeemannT. (2014). Prokka: rapid prokaryotic genome annotation. Bioinformatics30, 2068–2069. doi: 10.1093/bioinformatics/btu153SendiP.ProctorR. A. (2009). Staphylococcus aureus as an intracellular pathogen: the role of small colony variants. Trends Microbiol.17, 54–58. doi: 10.1016/j.tim.2008.11.004SinhaB.FraunholzM. (2010). Staphylococcus aureus host cell invasion and post-invasion events. Int. J. Med. Microbiol.300, 170–175. doi: 10.1016/j.ijmm.2009.08.019SonessonA.PrzybyszewskaK.ErikssonS.MörgelinM.KjellströmS.DaviesJ.. (2017). Identification of bacterial biofilm and the Staphylococcus aureus derived protease, staphopain, on the skin surface of patients with atopic dermatitis. Sci. Rep.7, 8689. doi: 10.1038/s41598-017-08046-2SongY.RubioA.JayaswalR. K.SilvermanJ. A.WilkinsonB. J. (2013). Additional routes to Staphylococcus aureus daptomycin resistance as revealed by comparative genome sequencing, transcriptional profiling, and phenotypic studies. PloS One8, e58469. doi: 10.1371/journal.pone.0058469StaudenmaierL.FockenJ.SchlattererK.KretschmerD.SchittekB. (2022). Bacterial membrane vesicles shape Staphylococcus aureus skin colonization and induction of innate immune responses. Exp. Dermatol.31, 349–361. doi: 10.1111/exd.14478StepanovićS.VukovićD.HolaV.Di BonaventuraG.DjukićS.CirkovićI.. (2007). Quantification of biofilm in microtiter plates: overview of testing conditions and practical recommendations for assessment of biofilm production by staphylococci. APMIS115, 891–899. doi: 10.1111/j.1600-0463.2007.apm_630.xSybilskiA. J.RaciborskiF.LipiecA.TomaszewskaA.LusawaA.Samel-KowalikP.. (2015). Epidemiology of atopic dermatitis in Poland according to the epidemiology of allergic disorders in Poland (ECAP) study. J. Dermatol.42, 140–147. doi: 10.1111/1346-8138.12731SyedA. K.ReedT. J.ClarkK. L.BolesB. R.KahlenbergJ. M. (2015). Staphlyococcus aureus phenol-soluble modulins stimulate the release of proinflammatory cytokines from keratinocytes and are required for induction of skin inflammation. Infect. Immun.83, 3428–3437. doi: 10.1128/IAI.00401-15TaskapanM. O.KumarP. (2000). Role of staphylococcal superantigens in atopic dermatitis: from colonization to inflammation. Ann. Allergy Asthma Immunol.84, 3–10. doi: 10.1016/S1081-1206(10)62731-7TauberM.BalicaS.HsuC. Y.Jean-DecosterC.LauzeC.RedoulesD.. (2016). Staphylococcus aureus density on lesional and nonlesional skin is strongly associated with disease severity in atopic dermatitis. J. Allergy Clin. Immunol.137, 1272–1274.e3. doi: 10.1016/j.jaci.2015.07.052TayebiZ.GoudarziH.DadashiM.GoudarziM. (2020). Genotype distribution of methicillin-susceptible Staphylococcus aureus clinical isolates in Iran: high multiresistant clonal complex 8. BMC Res. Notes13, 277. doi: 10.1186/s13104-020-05127-wThomsenM. C. F.HasmanH.WesthH.KayaH.LundO. (2017). RUCS: rapid identification of PCR primers for unique core sequences. Bioinformatics33, 3917–3921. doi: 10.1093/bioinformatics/btx526TomiN. S.KränkeB.AbererE. (2005). Staphylococcal toxins in patients with psoriasis, atopic dermatitis, and erythroderma, and in healthy control subjects. J. Am. Acad. Dermatol.53, 67–72. doi: 10.1016/j.jaad.2005.02.034TottéJ. E.van der FeltzW. T.HennekamM.van BelkumA.van ZuurenE. J.PasmansS. G. (2016). Prevalence and odds of Staphylococcus aureus carriage in atopic dermatitis: a systematic review and meta-analysis. Br. J. Dermatol.175, 687–695. doi: 10.1111/bjd.14566Truong-BolducQ. C.DunmanP. M.StrahilevitzJ.ProjanS. J.HooperD. C. (2005). MgrA is a multiple regulator of two new efflux pumps in Staphylococcus aureus. J. Bacteriol.187, 2395–2405. doi: 10.1128/JB.187.7.2395-2405.2005UhlemannA. C.OttoM.LowyF. D.DeLeoF. R. (2014). Evolution of community- and healthcare-associated methicillin-resistant Staphylococcus aureus. Infect. Genet. Evol.21, 563–574. doi: 10.1016/j.meegid.2013.04.030van BelkumA.VerkaikN. J.de VogelC. P.BoelensH. A.VerveerJ.NouwenJ. L.. (2009). Reclassification of Staphylococcus aureus nasal carriage types. J. Infect. Dis.199, 1820–1826. doi: 10.1086/599119WeidingerS.BeckL. A.BieberT.KabashimaK.IrvineA. D. (2018). Atopic dermatitis. Nat. Rev. Dis. Primers4, 1. doi: 10.1038/s41572-018-0001-zWichelhausT. A.SchäferV.BradeV.BöddinghausB. (1999). Molecular characterization of rpoB mutations conferring cross-resistance to rifamycins on methicillin-resistant Staphylococcus aureus. Antimicrob. Agents Chemother.43, 2813–2816. doi: 10.1128/aac.43.11.2813WickhamH. (2007). Reshaping data with the reshape package. J. Stat. Software21, 1–20. doi: 10.18637/jss.v021.i12WinstelV.XiaG.PeschelA. (2014). Pathways and roles of wall teichoic acid glycosylation in Staphylococcus aureus. Int. J. Med. Microbiol.304, 215–221. doi: 10.1016/j.ijmm.2013.10.009YuG.SmithD.ZhuH.GuanY.LamT. T. Y. (2017). Ggtree: an r package for visualization and annotation of phylogenetic trees with their covariates and other associated data. Methods Ecol. Evol.8 (1), 28–36. doi: 10.1111/2041-210X.12628ZareS.DerakhshandehA.HaghkhahM.NaziriZ.BroujeniA. M. (2019). Molecular typing of Staphylococcus aureus from different sources by RAPD-PCR analysis. Heliyon5, e02231. doi: 10.1016/j.heliyon.2019.e02231ZhaoH.HuF.JinS.XuX.ZouY.DingB.. (2016). Typing of panton-valentine leukocidin-encoding phages and lukSF-PV gene sequence variation in Staphylococcus aureus from China. Front. Microbiol.7. doi: 10.3389/fmicb.2016.01200