Next Article in Journal
Classical Microbiological Diagnostics of Bacteremia: Are the Negative Results Really Negative? What is the Laboratory Result Telling Us About the “Gold Standard”?
Previous Article in Journal
Strategies to Prevent Biofilm Infections on Biomaterials: Effect of Novel Naturally-Derived Biofilm Inhibitors on a Competitive Colonization Model of Titanium by Staphylococcus aureus and SaOS-2 Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transporters and Efflux Pumps Are the Main Mechanisms Involved in Staphylococcus epidermidis Adaptation and Tolerance to Didecyldimethylammonium Chloride

1
Department of Food Science and Technology, Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, SI-1000 Ljubljana, Slovenia
2
Department of Agronomy, Biotechnical Faculty, University of Ljubljana, Jamnikarjeva 101, SI-1000 Ljubljana, Slovenia
3
Omega d.o.o., Dolinškova 8, SI-1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Microorganisms 2020, 8(3), 344; https://doi.org/10.3390/microorganisms8030344
Submission received: 5 February 2020 / Revised: 27 February 2020 / Accepted: 28 February 2020 / Published: 29 February 2020
(This article belongs to the Section Antimicrobial Agents and Resistance)

Abstract

:
Staphylococcus epidermidis cleanroom strains are often exposed to sub-inhibitory concentrations of disinfectants, including didecyldimethylammonium chloride (DDAC). Consequently, they can adapt or even become tolerant to them. RNA-sequencing was used to investigate adaptation and tolerance mechanisms of S. epidermidis cleanroom strains (SE11, SE18), with S. epidermidis SE11Ad adapted and S. epidermidis SE18To tolerant to DDAC. Adaptation to DDAC was identified with up-regulation of genes mainly involved in transport (thioredoxin reductase [pstS], the arsenic efflux pump [gene ID, SE0334], sugar phosphate antiporter [uhpT]), while down-regulation was seen for the Agr system (agrA, arC, agrD, psm, SE1543), for enhanced biofilm formation. Tolerance to DDAC revealed the up-regulation of genes associated with transporters (L-cysteine transport [tcyB]; uracil permease [SE0875]; multidrug transporter [lmrP]; arsenic efflux pump [arsB]); the down-regulation of genes involved in amino-acid biosynthesis (lysine [dapE]; histidine [hisA]; methionine [metC]), and an enzyme involved in peptidoglycan, and therefore cell wall modifications (alanine racemase [SE1079]). We show for the first time the differentially expressed genes in DDAC-adapted and DDAC-tolerant S. epidermidis strains, which highlight the complexity of the responses through the involvement of different mechanisms.

Graphical Abstract

1. Introduction

Staphylococcus epidermidis is a member of the coagulase-negative Staphylococcus group, and it is the most common bacterial species of human skin and mucosa [1,2,3]. S. epidermidis is a leading cause of infections related to medical devices, and as highlighted by the European Centre for Disease Control, the resistance to multiple antimicrobial agents among S. epidermidis strains is a concerning trend that might limit treatment options for infections that are already difficult to treat [4]. A reduction in S. epidermidis prevalence would eliminate the reservoir of resistance genes and decrease the risk of further transfer of resistance genes and resistance-associated mobile genetic elements [5]. S. epidermidis strains are also known to easily adapt to different environmental changes. Thus, the presence of low concentrations of disinfectants will provide the opportunity for these bacteria to adapt and/or to develop tolerance [6,7].
S. epidermidis strains are one of the most common cleanroom isolates. Cleanrooms are places that are considered to be microbially reduced environments [8] with defined and controlled air quality (e.g., generally using high-efficiency particulate air-filtration systems), particles in the air, temperature and humidity. Cleanrooms ensure the quality and safety of many production areas, where contamination with microorganisms and particles is critical, such as in surgical units, manufacture of pharmaceuticals, food production and biosafety laboratories [8,9]. Therefore, disinfection of cleanrooms is a common procedure that is very important for maintenance of low levels of microbial burden.
Quaternary ammonium compounds (QACs) and didecyldimethylammonium chloride (DDAC) are commonly used as disinfectants, as they have broad-spectrum activities [10,11,12,13]. However, regular surface disinfection can result in low levels of QACs maintained in the environment. Therefore, bacteria might adapt towards an individual antimicrobial it is exposed to, or even develop tolerance to the disinfectants used [14]. It is also known that adaptation or tolerance to QACs can lead to cross-tolerance to other antimicrobial substances, such as other biocides and antibiotics [15].
Adaptation is usually a reversible biological property, such that it is lost with the removal of the disinfectant. However, bacterial tolerance to a disinfectant is a permanent change that allows tolerant cells to grow at higher concentrations of the disinfectant than non-tolerant cells. Here, tolerance remains even after the removal of the disinfectant, and it can be transferred across generations [16]. Several bacterial adaptation and tolerance mechanisms are known to be involved, such as modulation of cell growth rate, biofilm formation, efflux pump activity, fatty-acid profile, and the cell size, cell membranes and cell wall [16,17,18,19]. However, in general, tolerance development to QACs is still not fully understood [13].
RNA-sequencing (RNA-seq) is a common tool for analysis of differentially expressed genes (DEGs) [20] and we used this method to identify candidate genes and the underlying S. epidermidis biological pathways involved in adaptation and development of tolerance to DDAC. The transcriptomes of staphylococci (S. epidermidis, S. aureus) have been analyzed after exposure to environmental stresses [7,21,22], but none of these studies have been directly associated to cleanroom strains and adaptation or tolerance to DDAC. Adaptation of the OJ82 seafood strain of Staphylococcus spp. to high salt concentrations resulted in higher expression of genes involved in basic cellular processes, such as translation (ribosomal proteins, tRNA synthetase, translation initiation factor genes) and ATP production, together with membrane synthesis, transport systems, and metabolic pathways [22]. Similar responses in nitrogen metabolic pathways in S. aureus and S. epidermidis human strains were observed as a response to sub-inhibitory concentrations of sapienic acid [21]. Investigation of the food-strain transcriptome of S. aureus grown in the presence of disinfectants (ethanol, chloramine T, their combination) revealed significant up-regulation of genes associated with biofilm formation, transport systems, metabolism pathways and membrane composition [7]. To our knowledge, the precise mechanisms involved in the adaptation and tolerance of cleanroom S. epidermidis strains to DDAC has not been studied.
In the present study, we used the adapted and tolerant S. epidermidis cleanroom strains that we had previously obtained by their exposure to sub-inhibitory concentrations of DDAC [19]. We aimed to gain insight into the transcriptomic changes for the biological processes involved in the adaptation and tolerance mechanisms of these S. epidermidis cleanroom strains. Using RNA-seq, we analyzed the differential expression of genes through a comparison of the transcriptomes of these adapted and tolerant strains (SE11Ad, SE18To) with the transcriptomes of their original strains before adaptation to DDAC as control strains (SE11, SE18, respectively). To the best of our knowledge, this is the first time that the adaptation and tolerance response mechanisms of S. epidermidis strains have been investigated using RNA-seq to study adaptation and tolerance to DDAC, as one of the commonly used QAC disinfectants.

2. Materials and Methods

2.1. Staphylococcus epidermidis Strains and Growth Conditions

The S. epidermidis SE11 and S. epidermidis SE18 cleanroom strains were used as the control strains (henceforth referred to as SE11, SE18). These strains had been exposed to DDAC (Merck, Darmstadt, Germany) in our previous study [19]. The DDAC-adapted SE11Ad and DDAC-tolerant SE18To S. epidermidis strains obtained were used in the present study (henceforth referred to as SE11Ad, SE18To). The DDAC minimal inhibitory concentration (MIC) against SE11Ad was 18.0 µg/mL, while for the control SE11 this was lower, at 4.5 µg/mL. The DDAC MIC against SE18To was 9.0 µg/mL, and against SE18 it was 0.6 µg/mL [19].
For RNA isolation, SE11 and SE18 (1 mL) were grown in Tryptic Soya Broth (TSB; Oxoid, Hampshire, UK) (9 mL), while SE11Ad and SE18To were grown in TSB with 18.0 µg/mL and 9.0 µg/mL DDAC, respectively. All of these S. epidermidis strains were incubated at 37 °C for 4.5 h with shaking (120 rpm), to a concentration of approximately 109 cells/mL.

2.2. Next-Generation Sequencing

2.2.1. RNA Isolation and Quantification

The total RNA of SE11, SE11Ad, SE18 and SE18To were isolated using PureLink RNA mini kit (Thermo Fisher Scientific, Carlsbad, CA, USA), following the manufacturer protocol, with some minor modifications. Cell lysis was performed in lysis solution with 1% 2-mercaptoethanol (98%; Sigma-Aldrich Chemie, Taufkirchen, Germany), and the cells were disrupted using a laboratory mill (Retsch MM400; Retsch GmBH, Haan, Germany) at the maximum speed (1800 rpm) for 35 s, followed by cooling on ice (2 min). During the extraction step, DNAse treatment was performed using PureLink DNAse sets (Thermo Fisher Scientific). The total RNA yield and quality were investigated with Qubit RNA HS assay kits (Thermo Fisher Scientific), according to manufacturer protocol, using a fluorometer (Qubit v4; Thermo Fisher Scientific). The messenger RNA (mRNA) fraction from the total RNA was enriched using beads (NEXTflex PolyA; Perkin Elmer, Waltham, MA, USA), according to the manufacturer protocol.

2.2.2. Ion Torrent Library Preparation and Sequencing

RNA-sequencing next-generation sequencing libraries were prepared using Ion Total RNA-seq kits v2 (Thermo Fisher Scientific), following the manufacturer instructions and including four different barcodes. The amounts and size distributions of the prepared library fragments were determined (Labchip GX; Perkin Elmer). Strand-specific sequencing was performed using the Ion Proton system (Thermo Fisher Scientific). Emulsion PCR and the enrichment steps for the four pooled libraries were carried out using Ion P Hi-Q OT2 200 kits (Thermo Fisher Scientific), as described by the manufacturer. Library templates were sequenced with Ion P Hi-Q Sequencing 200 kits (Thermo Fisher Scientific) using a single Ion PI Chip kit v3 (Thermo Fisher Scientific). Signal processing, base calling, and adapter trimming were performed with the Torrent Suite software version 5.10.1 (Thermo Fisher Scientific). The sequences of the raw data have been deposited with the NCBI Sequence Read Archive (SRA) under BioProject SRA accession number PRJNA596062 (SE11, SRR10715330; SE11Ad, SRR10715329; SE18, SRR10715328; SE18To, SRR10715327).

2.3. Data Analysis

The bioinformatics analysis was performed using the CLC Genomics Workbench (version 12.0.3) and the CLC Genomics Server (version 11.0.2) tools. The basic quality control of the reads was performed using the ‘QC for Sequencing Reads’ tool (Qiagen, Redwood City, CA, USA). The reference genome used for the mapping sequences was the S. epidermidis ATCC12228 genome from NCBI (NC_004461), including its six plasmids (NC_005008, NC_005007, NC_005006, NC_005005, NC_005004, NC_005003) with the gene coordinates. For screening the qac genes, the S. epidermidis pHOB1 plasmid reference sequence was used (NZ_CP018843.1). The mapping of the sequencing reads to the reference sequences was done with the ‘RNA-seq analysis’ tool using the default settings with the strand-specific option set as ‘Forward’. The corresponding mapping profiles were used in a statistical differential expression test implemented in the ‘Differential expression for RNA-seq tool’, which uses multi-factorial statistics based on a negative binomial generalized linear model. The DEGs were chosen based on fold change (FC) >1 and p <0.05. SE11Ad was analyzed in comparison to its control, SE11, and SE18To to its control, SE18. The gene ontologies for biological processes, cellular components and molecular functions for these genes were searched in the UniProt database (https://www.uniprot.org/), based on gene model accession numbers and according to their functional annotations using the Blast2GO command line (version 1.3.0). These DEGs were further analyzed according to the Kyoto Encyclopaedia of Genes and Genome (KEGG) pathways, in the KEGG mapper tool (https://www.genome.jp/kegg/mapper.html).

3. Results and Discussion

Transcriptomic analyses were performed for SE11Ad (DDAC-adapted S. epidermidis SE11 strain) and SE18To (DDAC-tolerant S. epidermidis SE18 strain). Strand-specific RNA-seq was analyzed using the Ion Torrent technology. Altogether, after default Torrent Suite quality filtering and trimming, 42.26 M reads were obtained for the four transcriptomes studied (Table 1). The majority of the reads mapped to the reference S. epidermidis ATCC12228 genome with six plasmids and the S. epidermidis pHOB1 plasmid (Table 1). The data are available in the SRA under BioProject accession number PRJNA596062.
The mechanisms of S. epidermidis that are involved in adaptation and tolerance to DDAC were studied for the first time according to the DEGs and their involvement in biological pathways and processes. SE11Ad showed 720 DEGs, with 385 up-regulated and 335 down-regulated. SE18To showed 1013 DEGs, with 519 up-regulated and 494 down-regulated (Supplementary Table S1).

3.1. Response Mechanisms to DDAC of the DDAC-adapted S. epidermidis Strain, SE11Ad

Detailed analysis of SE11Ad revealed 111 significant DEGs when analyzed in comparison to its control SE11. Among these, 37 genes were up-regulated and 74 genes were down-regulated (Supplementary Table S2). According to gene ontology analysis, different metabolic and biosynthetic processes, transport systems, and DNA transformations were the most significant biological processes (Figure 1A). ‘Integral component of membrane’ was the main cellular component (Figure 1B), and ‘ATP-binding’, ‘Transporter activity’, ‘DNA-binding’ and ‘Metal-ion binding’ were the main molecular function activities (Figure 1C) that characterized these DEGs in SE11Ad.
The KEGG pathway analysis defined the most representative cellular pathways of these DEGs involved in SE11Ad adaptation to DDAC (Supplementary Table S2). The highest number of DEGs, along with their higher fold changes, were identified among the metabolic and biosynthetic pathways. Glycerol ester hydrolase is involved in lipid metabolism, and showed the highest fold change (gehD [gene ID, SE0185], -495.5-fold). Due to its lipase activity, it was defined as an important virulence factor that enables S. epidermidis to grow on a lipid-rich skin environment [23,24]. This gene was also down-regulated in an S. epidermidis RNAIII mutant with a reduced quorum sensing (Agr) system, which suggested that gehD is Agr-dependent [25]. Seventeen other genes were also involved in different metabolic pathways and in biosynthesis of secondary metabolites and antibiotics (Table 2). Three genes involved in the shikimic acid pathway were down-regulated (chorismate mutase [SE1599], -26.1-fold; shikimate kinase [SE1224], -6.2-fold; 3-phosphoshikimate 1-carboxyvinyltransferase [SE1153], -5.8-fold). This metabolic pathway leads to the biosynthesis of chorismic acid, a precursor of aromatic amino acids, folate, ubiquinone, and other isoprenoid quinones [26]. Furthermore, dihydrolipoamide dehydrogenase was down-regulated (SE1199, -11.1-fold), similar to that seen as the manuka honey effect in S. aureus [27]. Reduced expression of these metabolism-associated enzymes might be related to the SE11Ad adaptation response as decreased energy consumption of the cells for metabolism.
The down-regulated fused signal recognition particle receptor (ftsY [SE0910]) is part of the S. epidermidis Sec secretion system, and is involved in cell division, as a highly energy demanding process. The expression of ftsY was also reduced when S. epidermidis cells were exposed to α-mangostin, a natural xanthone with antimicrobial activity [28]. This indicates another energy saving process as an adaptation response.
Another adaptation response identified here for SE11Ad was reduced production of secondary metabolites, as in this way cells can save energy to ‘streamline’ their primary metabolism [29,30]. The down-regulated Ipk gene (SE2288, -21.7-fold) is involved in terpenoid biosynthesis, which is essential for cell survival. Terpenoid biosynthesis enzymes have already been investigated as potential drug targets in staphylococci, and decreased expression of this biosynthesis pathway is known to lead to down-regulation of primary metabolism and a shift in the regulation of virulence factors and cell-wall biosynthesis genes [31]. The down-regulated uroporphyrinogen III synthase (SE1344, -5.9-fold) is involved in hem biosynthesis [32], together with cytochrome c oxidase assembly protein subunit 15, hem synthase (SE0814, -8.5-fold) and hemC (SE1345, -12.8-fold). Hem biosynthesis reduction has already been shown for biofilm and persister cells (small colony variants) of S. aureus [33]. MiaA was also down-regulated (SE0981, -8.9-fold), and it encodes tRNA dimethyl allyl transferase, which is responsible for tRNA biosynthesis. Strong repression of genes responsible for DNA synthesis and its metabolism was shown in S. aureus exposed to the antimicrobial mupirocin, along with the same for those for translation processes [34]. Similar effects were reported in S. aureus exposed to ethanol, where genes for transcription, translation and nucleotide biosynthesis were generally down-regulated [35]. It has been suggested that exposure to this antimicrobial agent triggers a global adaptation response that is intended to adjust the cell biosynthetic pathways to the availability of energy and the required precursors.
Genes involved in quorum sensing were highly down-regulated (Table 2). Among these, four were from the accessory gene regulator (Agr) system (agrA [SE1638], -81.3-fold; agrC [SE1637], -77.0-fold; agrD [SE1636], -391.5-fold) and phenol soluble modulin β1/β2 (SE0484, -285.1-fold). The Agr system is considered as the main quorum-sensing mechanism in Gram-positive bacteria, including S. epidermidis strains [36]. The reduced expression of these genes might be directly associated with the adaptation of SE11Ad to DDAC, as down-regulation of the Agr system is associated with increased biofilm formation and decreased formation of phenol-soluble modulins (low-molecular-weight toxins) and other exoenzymes [25,36]. Through the Agr system, cells maintain energy for the vital activities connected to biofilm formation, which can protect them against DDAC. This regulation is usually described as initial bacterial colonization and development of chronic human infections [25,36]. However, for SE11Ad, the reduced expression of genes involved in the Agr system was identified as the mechanism that might have a large contribution in adaptation to DDAC, as reduced expression of the genes of the Agr system enhance biofilm formation as a response to DDAC. Indeed, stronger biofilm formation was detected in SE11Re, in comparison to the control SE11 [19]. The Agr system also regulates some extracellular proteases, such as V8 protease (SE1543, -534.6-fold). This protease is involved in extracellular protein cleavage and has been associated with protein-based biofilm removal [37,38]. As reported previously, in S. epidermidis, a reduced Agr system is connected with reduced production of extracellular enzymes, like proteases [25]. Down-regulation of V8 protease appears to contribute to greater biofilm formation, as an adaptation response to DDAC.
Other highly down-regulated genes included putative transposases for insertion-sequence-like elements (SE0079, -1795.91-fold; SE0090, -636.85-fold; BUM85_RS13075, -1447.36-fold) (Supplementary Table S2). Transposases are normally responsible for genetic recombination and the repair of genetic material, along with another down-regulated gene that encodes the recombination protein recR (recR, -13.86-fold). S. epidermidis has many mobile genetic elements that enable mobilisation and horizontal gene transfer. Insertion sequences are usually small DNA sequences that include one or two transposase genes [39]. Often, these are associated with resistance genes, which can randomly move to new genomic locations. This is one of the most concerning phenomena of the spread of resistance genes among bacterial cells of the same, or even different, species, and this genome plasticity provides bacteria with rapid adaptation to changing environmental conditions [39,40,41]. This mechanism is known to occur in stressful environments [42], and it might increase the adaptation of SE11Ad to DDAC. As high repression of putative transposase for insertion-sequence-like elements was seen for SE11Ad, they appear to be involved in adaptation.
There were eight highly up-regulated genes for SE11Ad among the transporters, which are responsible for the transport of different substances, such as mineral and organic ions, oligosaccharides, lipids, sterols, peptides, proteins, metals, amino acids and drugs and other antimicrobials (Table 2). Transporters have been suggested to be defined as virulence factors, as they are responsible for the transport of numerous substrates that are essential for virulence and survival in stress environments [43]. Thioredoxin reductase (pstS [SE1070]) was the most up-regulated gene for SE11Ad, with the highest fold-change (124.9-fold). It is part of a transporter involved in the thiol-specific redox system. This system is responsible for maintenance of cell homeostasis and protection of the cell from oxidising compounds, and therefore it is involved in cell defense against toxic compounds [44]. A phosphate transporter gene (pstC [SE1069]) was also highly up-regulated (38.5-fold). Together with PstA, this PstC protein forms a channel for the transport of phosphate across the cytoplasmic membrane [45]. Additionally, pstB (SE1067) was up-regulated (8.00-fold). Usually, this phosphate transport is more expressed in environments with limited phosphate. When S. aureus was exposed to the antimicrobial peptide ranalexin, similar effects were observed [46], and this was thus suggested as a Staphylococcus adaptation mechanism. Among the ABC transporters, genes for amino-acid transport were also up-regulated, such as the D-methionine transport system permease (metI [SE2321]), the methionine transport ATP-binding protein (metN [SE2322]), and the L-cystine transport system permease (tcyB [SE1992]). A supply of these particular amino acids is needed for growth and survival of SE11Ad [47]. The oligopeptide transport system permease (oppC [SE0681]) was also up-regulated, which is associated with the transport of small peptides for cell-wall biosynthesis [48]. This is part of the opp operon, which is necessary for survival for S. aureus [49].
As well as these transporters, the efflux systems were involved in the response of SE11Ad to this DDAC adaptation (Table 2). The up-regulation of genes involved in different efflux systems have also been reported previously after exposure and adaptation to antimicrobials, such as for S. aureus after exposure to silver nanoparticles [50], for Pseudomonas aeruginosa adapted to benzalkonium chloride [51], and for S. epidermidis adapted to ethidium bromide [52]. An important efflux system for the arsenic efflux pump (SE0334) showed an increase in 11.8-fold (Table 2). This efflux pump actively exports toxic arsenic compounds from cells [53]. The arsenic efflux system has multidrug-resistance members [54], and according to these data, we suggest that this system contributes to the active export of DDAC in SE11Ad, although further supporting evidence is needed. Among the two component systems, the uhpT gene that encodes a sugar-phosphate antiporter was highly up-regulated (11.5-fold). This is involved in phosphomycin resistance in methicillin-resistant S. aureus strains [55].
According to the metabolic pathways, the most up-regulated gene in metabolism encodes carbamate kinase (arc [SE0228], 20.7-fold). This is involved in arginine metabolism through the arginine deiminase pathway and is known as an important pathogenesis mechanism as its final product is ammonia, which increases cytoplasmic pH and therefore protects cells from toxic acidic conditions and maintains cytosolic pH homeostasis [56].

3.2. Response Mechanisms to DDAC of the DDAC-tolerant S. epidermidis Strain, SE18To

Detailed analysis of SE18To identified 42 DEGs (FC >5, p <0.05), among which 21 were up-regulated and 21 were down-regulated (Supplementary Table S2). Gene ontology analysis revealed different metabolism and biosynthetic processes and transport systems as the most significant biological processes associated with these DEGs (Figure 2A). According to the cellular component, the majority of the DEG products were ‘Integral components of membrane’ (Figure 2B), and ‘Metal-ion binding’, ‘Transporter activity’ and ‘DNA-binding’ were the main molecular functions that characterized these DEGs in SE18To (Figure 2C).
According to KEGG analysis, the most down-regulated genes were identified as among those involved in amino-acid biosynthesis. The dapE gene (SE2220, 3-4.0-fold) encodes succinyl-diaminopimelate desuccinylase (Table 3). This enzyme is involved in the biosynthesis of lysine, an important amino acid in bacterial cell-wall synthesis [57]. SE1079 was highly down-regulated (-16.9-fold) and encodes alanine racemase, which is involved in cell-wall synthesis and growth in S. aureus [58]. It converts L-alanine to D-alanine, which is a precursor of the pentapeptide moiety that links the glycan chains in peptidoglycan. A reduction in these peptidoglycan linkages results in greater elasticity of the S. aureus cell wall [59]. Alanine racemase has also been related to enterococci resistance against vancomycin, an antibiotic that targets the cell wall [18]. Modification of the peptidoglycan and bacterial cell wall are highly regulated mechanisms, and these can occur as responses to different environmental stimuli, with the basic aim being to protect the inner content of cells [16,48]. Such cell-membrane and cell-wall modifications usually result in thickening of the staphylococcal cell wall as a resistance mechanism, to help protect the cell against antimicrobials [18]. Other genes involved in amino-acid biosynthetic pathways were down-regulated, including (SE0275, -7-fold); this gene encodes 1-(5-phosphoribosyl)-5-[(5-phosphoribosylamino) methylideneamino] imidazole-4-carboxamide isomerase, which is involved in L-histidine biosynthesis. Down-regulation of hisA was also seen for Staphylococcus xylosus after exposure to aspirin [60]. Ketol-acid reductoisomerase (SE1657, -8.1-fold) is involved in branched-chain amino-acid biosynthesis [61]. Cystathionine β-lyase (metC [SE2380], -7.8-fold) is an enzyme that primarily converts L-cystathionine to homocysteine, pyruvate and ammonia, and it is involved in methionine biosynthesis [62]. A reduction in the basic metabolic processes and growth of cells is a survival strategy [63]. This down-regulation of the biosynthesis of amino acids and secondary metabolites indicates major changes in the amino-acid profile. Modifications to amino acids have been shown to be a response mechanism for a changing environment. This was shown by an analysis of the amino-acid profile of S. aureus after exposure to hydrogen peroxide [64]. Modifications to the amino-acid profile influence many biological pathways in cells, as amino acids are required for protein synthesis, and used for synthesis of the cell-wall and nucleic acids, and they can be used as energy substrates [64]. Reductions in cytoplasmic amino acids has been associated with slower rates of metabolism, to conserve energy for basal metabolism [65]. Modifications to cell membranes and the cell wall are known to be important to maintain cellular homeostasis, by the protection of cells from changes in the environment. This is achieved by maintenance of cell shape and turgor pressure, and is involved in cell division, energy production and permeability regulation [18].
The most up−regulated gene for SE18To was SE0875, which encodes uracil permease (27.6−fold), although this does not have a defined KEGG status (Supplementary Table S2). This is an electrochemically driven transporter that is responsible for transport of the nucleotide uracil. In a vancomycin−resistant S. aureus strain, this gene increased purine biosynthesis [66]. A shift in purine metabolism is one of the most common resistance responses in methicillin−resistant S. aureus [67]. Overall, up−regulation of transport, and particularly of the ABC transporters, was also reported in an analysis of the transcriptomes of ciprofloxacin−resistant S. aureus strains [68]. The transporter gene for L−cysteine, tcyB, was among the most up−regulated transporters in SE18To (SE1992, 16.0−fold). Together, tcyA, tcyB and tcyC encode the ABC transporter, which is responsible for uptake of cysteine, as required for cell growth. It has been suggested that in an environment with limited cysteine, the genes for this cysteine transporter are up−regulated [69]. As the biosynthesis of amino acids in SE18To was reduced, import from the environment provides the source for the cell needs. Other transport systems were also up−regulated, including the zinc transport system permease (znuB [SE1242], 6.8−fold), a part of the zinc transporter that is responsible for transporting zinc and other metals into and out of the cell, and therefore for the regulation of cell metal homeostasis. The up−regulated lmrP gene (SE0242, 7.0−fold; Supplementary Table S2) encodes a transporter protein that is an orthologue of the Lactococcus lactis LmrP multidrug efflux protein. This actively exports cationic amphiphilic substrates out of the cell, and is predicted to also export drugs in S. epidermidis [70,71]. Another multidrug efflux pump that was up−regulated was the arsenic efflux pump (arsB [SE0135], 5.97−fold). This energy−dependent efflux system belongs to the multi−drug−resistance protein family, which is involved in the active export of toxic compounds from cells [53,54].
The best described mechanism of bacterial tolerance to QACs involves the overexpression of efflux systems [16,72,73,74]. It is known that S. epidermidis has nearly 30 different drug efflux pumps, some of which are predicted to actively export QACs out of the cell [75]. The nature of DDAC and its non−specific membrane activity mean that it can enter cells by interacting with the phospholipid bilayer, which disrupts the membrane diffusive forces and ion gradients. This damage often leads to membrane lysis and cell death [13,76]. Active export of DDAC out of cells is a common mechanism of tolerance that eliminates the toxic DDAC from cells, to prevent cell damage [13]. This might be the major tolerance mechanism in SE18To.

3.3. L−cystine Transport and Arsenic Efflux Systems in DDAC−adapted and −tolerant S. epidermidis Strains

Results of the detailed analysis of SE11Ad and SE18To revealed that L−cystine transport and arsenic efflux are common responses to increased DDAC concentrations. The gene that encodes L−cystine permease (tcyB [SE1992]), which is part of the L−cystine amino−acid transport, was overexpressed in both SE11Ad and SE18To. In SE11Ad, tcyB was up−regulated 10.5−fold, and in SE18To, 16.0−fold. Up to 12 different amino acids are required for the growth of staphylococci, which include cysteine [77]. L−cysteine is also part of the low−molecular−weight thiols that have protective role in cells when they are exposed to reactive oxygen species, reactive electrophilic species, antibiotics and heavy metals [78]. Exposure of S. aureus to diamide leads to S−thiolation, and formation of disulphide bonds among proteins and low−molecular−weight thiols, and these products protect cytoplasmic proteins from irreversible oxidation [79]. Along with diamide, cysteine metabolism is connected to hydrogen peroxide, superoxide, nitric oxide, thiol−reactive electrophiles and metal ions [80]. We suggest that up−regulation of the tcyB gene that is involved in the uptake of L−cysteine is an important adaptation and tolerance mechanism for maintenance of cellular homeostasis in both SE11Ad and SE18To. These cells have greater need for L−cysteine to form low−molecular−weight thiols, and uptake of L−cysteine from the environment is preferred over de−novo synthesis [18,81].
The gene that encodes the arsenic efflux system (SE0334 for SE11Ad; arsB [SE0135] for SE18To) was also up−regulated in both of these strains. In SE11Ad, SE0334 was up−regulated 11.8−fold, and in SE18To, arsB was up−regulated 6.0−fold (Supplementary Table S2). There are many different arsenic efflux systems in cell membranes, and the ars operon is the most studied and known of these, although other novel arsenic efflux pumps have also been described [53]. Usually, these are considered as multidrug−resistance members of efflux systems [54]. These efflux systems represent another common mechanism when these S. epidermidis strains are adapted or tolerant to DDAC.

3.4. RNA−sequencing Results and Phenotypic Properties of These S. epidermidis Strains

The RNA−seq results (Table 2 and Table 3) show that genes involved in transport systems and efflux pumps were highly up−regulated in SE11Ad and SE18To. These results confirm our previous study, where we showed increased activity of the efflux pumps in SE11Ad and SE18To, based on an acridine−orange assay [19].
SE11Ad and SE18To showed down−regulation of genes for different metabolic pathways, including the biosynthesis of amino acids and secondary metabolites, whereby the cells will be saving energy for their basic metabolic processes through the reduction in processes that are not crucial for survival. The SE11Ad and SE18To cells showed decreased sizes after the DDAC treatments, in comparison to the control SE11 and SE18 cells, respectively, when examined using scanning electron microscopy [19]. Modifications to the metabolic pathways are related to decreasing cell size, although many different pathways were detected using the RNA−seq analysis here. The reduction in growth and regulation of the basal metabolic activity of cells represent adaptation responses and survival strategies, along with other metabolic changes, to reduce the harmful effects of antimicrobials [63,65,67].
With this RNA−seq analysis, we also identified the significant down−regulation of genes involved in the Agr system, and thus these cells increased biofilm formation, which is a known adaptation mechanism in S. epidermidis strains [82]. In our previous study, we showed that at the phenotypic level, SE11Ad increased biofilm formation compared to the control SE11 [19]. Biofilm formation protects cells from the harmful effects of antimicrobials as this limits the diffusion into the cells, and so a shift to biofilm production is one of the adaptation response mechanisms.

4. Conclusions

This analysis of the transcriptomes of these DDAC−adapted (SE11Ad) and DDAC−tolerant (SE18To) S. epidermidis strains revealed for the first time the most significant genes, and therefore mechanisms, that are involved in the adaptation and tolerance responses of S. epidermidis to this commonly used disinfectant.
The transcriptomic analysis of SE11Ad showed that adaptation to DDAC comprises different mechanisms, among which the most relevant were up−regulation of the different transporters (pstS, pstC, pstB, MetI, MetN) and efflux pump systems (arsenic efflux pump, tcyB), and down−regulation of the Agr system (AgrA, ArC, AgrD, psm, V8 protease), the consequence of which was greater biofilm formation. Tolerance to DDAC for SE18To was identified with the down−regulation of genes involved in the metabolic and biosynthetic pathways that are connected to peptidoglycan and cell−wall synthesis, while the most up−regulated genes were for transporters, which suggests that this represents another resistance mechanism. There was increased expression of genes associated with efflux systems (multidrug transporters, LmrP homologue, arsenic efflux pump) as one of the main mechanisms in this S. epidermidis tolerance to DDAC.
The continuing emergence of resistant bacterial strains is a well−known phenomenon, although tolerance to disinfectants is a relatively new aspect. Adaptation and tolerance development to DDAC (as a commonly used disinfectant) in S. epidermidis (the most common bacterium associated with cleanrooms) had not been studied to date. Using RNA−seq here, we have revealed that many genes are involved in these complex adaptation and tolerance mechanisms. We have shown that multiple mechanisms have co−developed during adaptation and tolerance to DDAC in these cleanroom S. epidermidis strains. These mechanisms of adaptation and tolerance to DDAC that we have identified for cleanroom S. epidermidis strains will help in the targeting of S. epidermidis strains in terms of the development of disinfectants and disinfection regimes, to control the incidence and spread of bacterial adaptation and tolerance to DDAC.

Supplementary Materials

The following Supplementary Materials are available online at https://www.mdpi.com/2076-2607/8/3/344/s1. Table S1. Differentially expressed genes (FC >1, p <0.05) in SE11Ad and SE18To, the didecyldimethylammonium chloride (DDAC)−adapted and DDAC−tolerant, respectively, S. epidermidis strains. Table S2. Analysis of the differentially expressed genes (FC >5, p <0.05) in SE11Ad and SE18To, the didecyldimethylammonium chloride (DDAC)−adapted and DDAC−tolerant, respectively, S. epidermidis strains, according to the gene ontology and KEGG pathways.

Author Contributions

Conceptualization, U.R., B.J., N.T. and S.K.; methodology, N.T. and S.K.; software, J.J. and S.K.; formal analysis, U.R., B.J., N.T. and S.K.; investigation, U.R., B.J., N.T. and S.K.; writing, reviewing and editing, U.R., B.J., N.T., S.K., M.K., J.J. and A.K.; visualization, J.J., U.R. and B.J. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Slovenian Research Agency Programme (P4−0116, J4−9299).

Acknowledgments

The authors acknowledge the bioinformatics hardware and software support by the Slovenian Research Agency Program P4−0077.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hellmark, B.; Söderquist, B.; Unemo, M.; Nilsdotter−Augustinsson, Å. Comparison of Staphylococcus epidermidis isolated from prosthetic joint infections and commensal isolates in regard to antibiotic susceptibility, agr type, biofilm production, and epidemiology. Int. J. Med. Microbiol. 2013, 303, 32–39. [Google Scholar] [CrossRef] [PubMed]
  2. Pinheiro, L.; Brito, C.I.; Cataneli Pereira, V.; de Oliveira, A.; Camargo, C.H.; de Lourdes Ribeiro de Souza da Cunha, M. Reduced susceptibility to vancomycin and biofilm formation in methicillin−resistant Staphylococcus epidermidis isolated from blood cultures. Mem. Inst. Oswaldo Cruz. 2014, 109, 871–878. [Google Scholar] [CrossRef] [PubMed]
  3. Balraadjsing, P.P.; de Jong, E.C.; van Wamel, W.J.B.; Zaat, S.A.J. Dendritic cells internalize Staphylococcus aureus more efficiently than Staphylococcus epidermidis, but do not differ in induction of antigen−specific T cell proliferation. Microorganisms 2020, 9, 19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. European Centre for Disease Prevention and Control. Multidrug−resistant Staphylococcus epidermidis. In Proceedings of the ECDC, Sweden, Stockholm, 8 November 2018. [Google Scholar]
  5. Xu, Z.; Misra, R.; Jamrozy, D.; Paterson, G.K.; Cutler, R.R.; Holmes, M.A.; Gharbia, S.; Mkrtchyan, H.V. Whole genome sequence and comparative genomics analysis of multi−drug resistant environmental Staphylococcus epidermidis ST59. Genes Genomes Genet. 2018, 8, 2225–2230. [Google Scholar] [CrossRef] [Green Version]
  6. Onyango, L.A.; Dunstan, R.H.; Roberts, T.K.; Macdonald, M.M.; Gottfries, J. Phenotypic variants of Staphylococci and their underlying population distributions following exposure to stress. PLoS ONE 2013, 8, e77614. [Google Scholar] [CrossRef] [Green Version]
  7. Slany, M.; Oppelt, J.; Cincarova, L. Formation of Staphylococcus aureus biofilm in the presence of sublethal concentrations of disinfectants studied via a transcriptomic analysis using transcriptome sequencing (RNA−seq). Appl. Environ. Microbiol. 2017, 83, e01643-17. [Google Scholar] [CrossRef] [Green Version]
  8. Mahnert, A.; Vaishampayan, P.; Probst, A.J.; Auerbach, A.; Moissl−Eichinger, C.; Venkateswaran, K.; Berg, G. Cleanroom maintenance significantly reduces abundance but not diversity of indoor microbiomes. PLoS ONE 2015, 10, e0134848. [Google Scholar] [CrossRef]
  9. Sandle, T.; Leavy, C.; Jindal, H.; Rhodes, R. Application of rapid microbiological methods for the risk assessment of controlled biopharmaceutical environments. J. Appl. Microbiol. 2014, 116, 1495–1505. [Google Scholar] [CrossRef] [Green Version]
  10. Jansen, A.C.; Boucher, C.E.; Coetsee, E.; Kock, J.L.F.; van Wyk, P.W.J.; Swart, H.C.; Bragg, R.R. The influence of didecyldimethylammonium chloride on the morphology and elemental composition of Staphylococcus aureus as determined by NanoSAM. Sci. Res. Essays 2013, 8, 152–160. [Google Scholar]
  11. Wessels, S.; Ingmer, H. Modes of action of three active substances: A review. Regul. Toxicol. Pharmacol. 2013, 67, 456–467. [Google Scholar] [CrossRef]
  12. Buffet−Bataillon, S.; Tattevin, P.; Maillard, J.−Y.; Bonnaure−Mallet, M.; Jolivet−Gougeon, A. Efflux pump induction by quaternary ammonium compounds (QAC) is a first step in the evolution of fluoroquinolone−resistance in bacteria. Future Microbiol. 2016, 11, 81–92. [Google Scholar]
  13. Morrison, K.R.; Allen, R.A.; Minbiole, K.P.C.; Wuest, W.M. More QACs, more questions: Recent advances in structure activity relationships and hurdles in understanding resistance mechanisms. Tetrahedron Lett. 2019, 60, 150935. [Google Scholar] [CrossRef]
  14. Smith, K.; Gemmell, C.G.; Hunter, I.S. The association between biocide tolerance and the presence or absence of qac genes among hospital−acquired and community−acquired MRSA isolates. J. Antimicrob. Chemother. 2008, 61, 78–84. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Capita, R.; Alonso−Calleja, C. Antibiotic−resistant bacteria: A challenge for the food industry. Crit. Rev. Food Sci. Nutr. 2013, 53, 11–48. [Google Scholar] [CrossRef] [PubMed]
  16. Jennings, M.C.; Minbiole, K.P.C.; Wuest, W.M. Quaternary ammonium compounds: An antimicrobial mainstay and platform for innovation to address bacterial resistance. ACS Infect. Dis. 2015, 1, 288–303. [Google Scholar] [CrossRef] [PubMed]
  17. Gadea, R.; Fernández Fuentes, M.Á.; Pérez Pulido, R.; Gálvez, A.; Ortega, E. Effects of exposure to quaternary−ammonium−based biocides on antimicrobial susceptibility and tolerance to physical stress in bacteria from organic foods. Food Microbiol. 2017, 63, 58–71. [Google Scholar] [CrossRef]
  18. Onyango, L.A.; Alreshidi, M.M. Adaptive metabolism in Staphylococci: Survival and persistence in environmental and clinical settings. J. Pathog. 2018, 2018, 1–11. [Google Scholar] [CrossRef] [Green Version]
  19. Ribič, U.; Polak, T.; Lušnic Polak, M.; Klančnik, A.; Jeršek, B. Adaptation response mechanisms of Staphylococcus epidermidis strains exposed to increasing concentrations of didecyldimethylammonium chloride. Microb. Drug Resist. 2020. ahead of print. [Google Scholar]
  20. Tian, X.; Zhang, L.; Feng, S.; Zhao, Z.; Wang, X.; Gao, H. Transcriptome analysis of apple leaves in response to powdery mildew (Podosphaera leucotricha) infection. Int. J. Mol. Sci. 2019, 20, 2326. [Google Scholar] [CrossRef] [Green Version]
  21. Moran, J.C.; Alorabi, J.A.; Horsburgh, M.J. Comparative transcriptomics reveals discrete survival responses of S. aureus and S. epidermidis to sapienic acid. Front. Microbiol. 2017, 8, 33. [Google Scholar] [CrossRef] [Green Version]
  22. Choi, S.; Jung, J.; Jeon, C.O.; Park, W. Comparative genomic and transcriptomic analyses of NaCl−tolerant Staphylococcus sp. OJ82 isolated from fermented seafood. Appl. Microbiol. Biotechnol. 2014, 98, 807–822. [Google Scholar] [CrossRef]
  23. Pourmand, M.R.; Clarke, S.R.; Schuman, R.F.; Mond, J.J.; Foster, S.J. Identification of antigenic components of Staphylococcus epidermidis expressed during human infection. Infect. Immun. 2006, 74, 4644–4654. [Google Scholar] [CrossRef] [Green Version]
  24. Longshaw, C.M.; Farrell, A.M.; Wright, J.D.; Holland, K.T. Identification of a second lipase gene, gehD, in Staphylococcus epidermidis; comparison of sequence with those of other staphylococcal lipases. Microbiology 2000, 146, 1419–1427. [Google Scholar] [CrossRef] [Green Version]
  25. Olson, M.E.; Todd, D.A.; Schaeffer, C.R.; Paharik, A.E.; Van Dyke, M.J.; Büttner, H.; Dunman, P.M.; Rohde, H.; Cech, N.B.; Fey, P.D.; et al. Staphylococcus epidermidis agr quorum−sensing system: Signal identification, cross talk, and importance in colonization. J. Bacteriol. 2014, 196, 3482–3493. [Google Scholar] [CrossRef] [Green Version]
  26. Zucko, J.; Dunlap, W.C.; Shick, J.M.; Cullum, J.; Cercelet, F.; Amin, B.; Hammen, L.; Lau, T.; Williams, J.; Hranueli, D.; et al. Global genome analysis of the shikimic acid pathway reveals greater gene loss in host−associated than in free−living bacteria. BMC Genom. 2010, 11, 628. [Google Scholar] [CrossRef] [Green Version]
  27. Packer, J.M.; Irish, J.; Herbert, B.R.; Hill, C.; Padula, M.; Blair, S.E.; Carter, D.A.; Harry, E.J. Specific non−peroxide antibacterial effect of manuka honey on the Staphylococcus aureus proteome. Int. J. Antimicrob. Agents 2012, 40, 43–50. [Google Scholar] [CrossRef]
  28. Sivaranjani, M.; Leskinen, K.; Aravindraja, C.; Saavalainen, P.; Pandian, S.K.; Skurnik, M.; Ravi, A.V. Deciphering the antibacterial mode of action of alpha−mangostin on Staphylococcus epidermidis RP62A through an integrated transcriptomic and proteomic approach. Front. Microbiol. 2019, 10, 150. [Google Scholar] [CrossRef]
  29. De Carvalho, C.C.C.R.; Fernandes, P. Production of metabolites as bacterial responses to the marine environment. Mar. Drugs 2010, 8, 705–727. [Google Scholar] [CrossRef]
  30. Raymond−Bouchard, I.; Tremblay, J.; Altshuler, I.; Greer, C.W.; Whyte, L.G. Comparative transcriptomics of cold growth and adaptive features of a aury− and steno−psychrophile. Front. Microbiol. 2018, 9, 1565. [Google Scholar] [CrossRef]
  31. Heuston, S.; Begley, M.; Gahan, C.G.M.; Hill, C. Isoprenoid biosynthesis in bacterial pathogens. Microbiology 2012, 158, 1389–1401. [Google Scholar] [CrossRef]
  32. Anzaldi, L.L.; Skaar, E.P. Overcoming the heme paradox: Heme toxicity and tolerance in bacterial pathogens. Infect. Immun. 2010, 78, 4977–4989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Cao, S.; Huseby, D.L.; Brandis, G.; Hughes, D. Alternative evolutionary pathways for drug−resistant small colony variant mutants in Staphylococcus aureus. MBIO 2017, 8, e00358-17. [Google Scholar] [CrossRef] [Green Version]
  34. Reiß, S.; Pané−Farré, J.; Fuchs, S.; François, P.; Liebeke, M.; Schrenzel, J.; Lindequist, U.; Lalk, M.; Wolz, C.; Hecker, M.; et al. Global analysis of the Staphylococcus aureus response to mupirocin. Antimicrob. Agents Chemother. 2012, 56, 787–804. [Google Scholar] [CrossRef] [Green Version]
  35. Pando, J.M.; Pfeltz, R.F.; Cuaron, J.A.; Nagarajan, V.; Mishra, M.N.; Torres, N.J.; Elasri, M.O.; Wilkinson, B.J.; Gustafson, J.E. Ethanol−induces stress response of Staphylococcus aureus. Can. J. Microbiol. 2017, 63, 745–757. [Google Scholar] [CrossRef]
  36. Le, K.Y.; Otto, M. Quorum−sensing regulation in staphylococci−an overview. Front. Microbiol. 2015, 6, 1174. [Google Scholar] [CrossRef] [Green Version]
  37. Martí, M.; Pilar Trotonda, M.; Ángeles Tormo−Más, M.; Vergara−Irigaray, M.; Cheung, A.L.; Lasa, I.; Penadés, J.R. Extracellular proteases inhibit protein−dependent biofilm formation in Staphylococcus aureus. Microbes Infect. 2010, 12, 55–64. [Google Scholar] [CrossRef] [Green Version]
  38. Loughran, A.J.; Atwood, D.N.; Anthony, A.C.; Harik, N.S.; Spencer, H.J.; Beenken, K.E.; Smeltzer, M.S. Impact of individual extracellular proteases on Staphylococcus aureus biofilm formation in diverse clinical isolates and their isogenic sarA mutants. Microbiologyopen 2014, 3, 897–909. [Google Scholar] [CrossRef] [Green Version]
  39. Partridge, S.R.; Kwong, S.M.; Firth, N.; Jensen, S.O. Mobile genetic elements associated with antimicrobial resistance. Clin. Microbiol. Rev. 2018, 31, e00088-17. [Google Scholar] [CrossRef] [Green Version]
  40. Winstel, V.; Liang, C.; Sanchez−Carballo, P.; Steglich, M.; Munar, M.; Bröker, B.M.; Penadés, J.R.; Nübel, U.; holst, O.; Dandekar, T.; et al. Wall teichoic acid structure governs horizontal gene transfer between major bacterial pathogens. Nat. Commun. 2013, 4, 2345. [Google Scholar] [CrossRef]
  41. Vandecraen, J.; Chandler, M.; Aersten, A.; Van Houdt, R. The impact of insertion sequences on bacterial genome plasticity and adaptability. Crit. Rev. Microbiol. 2017, 43, 709–730. [Google Scholar] [CrossRef]
  42. Matic, I.; Taddei, F.; Radman, M. Survival versus maintenance of genetic stability: A conflict of priorities during stress. Res. Microbiol. 2004, 155, 337–341. [Google Scholar] [CrossRef]
  43. Tanaka, K.J.; Song, S.; Mason, K.; Pinkett, H.W. Selective substrate uptake: The role of ATP−binding cassette (ABC) importers in pathogenesis. Biochim. Biophys. Acta 2018, 1860, 868–877. [Google Scholar] [CrossRef]
  44. Gómez−Gómez, B.; Pérez−Corona, T.; Mozzi, F.; Pescuma, M.; Madrid, Y. Silac−based quantitative proteomic analysis of Lactobacillus reuteri CRL 1101 response to the presence of selenite and selenium nanoparticles. J. Proteom. 2019, 195, 53–65. [Google Scholar] [CrossRef]
  45. Wang, X.; Han, H.; Lv, Z.; Lin, Z.; Shang, Y.; Xu, T.; Wu, Y.; Zhang, Y.; Qu, D. PhoU2 but not PhoU1 as an important regulator of biofilm formation and tolerance to multiple stresses by participating in various fundamental metabolic processes in Staphylococcus epidermidis. J. Bacteriol. 2017, 199, e00219-17. [Google Scholar] [CrossRef] [Green Version]
  46. Overton, I.M.; Graham, S.; Gould, K.A.; Hinds, J.; Botting, C.H.; Shirran, S.; Barton, G.J.; Coote, P.J. Global network analysis of drug tolerance, mode of action and virulence in methicillin−resistant S. aureus. BMC Syst. Biol. 2011, 5, 68. [Google Scholar] [CrossRef] [Green Version]
  47. Krismer, B.; Liebeke, M.; Janek, D.; Nega, M.; Rautenberg, M.; Hornig, G.; Unger, C.; Weidenmaier, C.; Lalk, M.; Peschel, A. Nutrient limitation governs Staphylococcus aureus metabolism and niche adaptation in the human nose. PLoS ONE 2014, 10, e1003862. [Google Scholar] [CrossRef]
  48. Vermassen, A.; Leroy, S.; Talon, R.; Provot, C.; Popowska, M.; Desvaux, M. Cell wall hydrolases in bacteria: Insight on the diversity of cell wall amidases, glycosidases and peptidases toward peptidoglycan. Front. Microbiol. 2019, 10, 331. [Google Scholar] [CrossRef]
  49. Chaffin, D.O.; Taylor, D.; Skerrett, S.J.; Rubens, C.E. Changes in the Staphylococcus aureus transcriptome during early adaptation to the lung. PLoS ONE 2012, 7, e41329. [Google Scholar] [CrossRef] [Green Version]
  50. Singh, N.; Rahwade, J.; Paknikar, K.M. Transcritpome analysis of silver nanoparticles treated Staphylococcus aureus reveals potential targets for biofilm formation. Colloids Surf. B Biointerface 2019, 175, 487–497. [Google Scholar] [CrossRef]
  51. Kim, M.; Hatt, J.K.; Weigand, M.R.; Krishnan, R.; Pavlostathis, S.G.; Konstantinidis, K.T. Genomic and transcriptomic insights into how bacteria withstands high concentrations of benzalkonium chloride biocides. Appl. Environ. Microbiol. 2018, 84, e00197. [Google Scholar] [CrossRef] [Green Version]
  52. Santos Costa, S.; Viveiros, M.; Pomba, C.; Couto, I. Active antimicrobial efflux in Staphylococcus epidermidis: Building up of resistance to fluoroquinolones and biocides in a major opportunistic pathogen. J. Antimicrob. Chemother. 2017, 73, 320–324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Ben Fekih, I.; Zhang, C.; Ping Li, Y.; Zhao, Y.; Alwathnani, H.A.; Saquib, Q.; Rensing, C.; Cervantes, C. Distribution of arsenic resistance genes in prokaryotes. Front. Microbiol. 2018, 9, 2473. [Google Scholar] [CrossRef] [PubMed]
  54. Yang, H.−C.; Fu, H.−L.; Lin, Y.−F.; Rosen, B.P. Pathways of arsenic uptake and efflux. Curr. Top. Membr. 2012, 69, 325–358. [Google Scholar]
  55. Fu, Z.; Ma, Y.; Chen, C.; Guo, Y.; Hu, F.; Liu, Y.; Xu, X.; Wang, M. Prevalence of Fosfomycin resistance and mutations in murA, glpT, and uhpT in methicillin−resistant Staphylococcus aureus strains isolated from blood and cerebrospinal fluid samples. Front. Microbiol. 2016, 6, 1544. [Google Scholar] [CrossRef]
  56. Xiong, L.; Teng, J.L.L.; Botelho, M.G.; Lo, R.C.; Lau, S.K.P.; Woo, P.C.Y. Arginine metabolism in bacterial pathogenesis and cancer therapy. Int. J. Mol. Sci. 2016, 17, 363. [Google Scholar] [CrossRef] [Green Version]
  57. Heath, T.K.; Lutz, M.R.; Reidl, C.T.; Guzman, E.R.; Herbert, C.A.; Nocek, B.P.; Holz, R.C.; Olsen, K.W.; Ballicora, M.A.; Becker, D.P. Practical spectrophotometric assay for the dapE−encoded N−succinyl−L,L−diaminopimelic acid desuccinylase, a potential antibiotic target. PLoS ONE 2018, 13, e0196010. [Google Scholar] [CrossRef] [Green Version]
  58. Ciustea, M.; Mootien, S.; Rosato, A.E.; Perez, O.; Cirillo, P.; Yeung, K.R.; Ledizet, M.; Cynamon, M.H.; Aristoff, P.A.; Koski, R.A.; et al. Thiadiazolidinones: A new class of alanine racemase inhibitors with antimicrobial activity against methicillin−resistant S. aureus. Biochem. Pharmacol. 2012, 83, 368–377. [Google Scholar] [CrossRef] [Green Version]
  59. Loskill, P.; Pereira, P.M.; Jung, P.; Bischoff, M.; Herrmann, M.; Pinho, M.G.; Jacobs, K. Reduction of the peptidoglycan crosslinking causes a decrease in stiffness of the Staphylococcus aureus cell envelope. Biophys. J. 2014, 107, 1082–1089. [Google Scholar] [CrossRef] [Green Version]
  60. Xu, C.−G.; Yang, Y.−B.; Zhou, Y.−H.; Hao, m.−Q.; Ren, Y.−Z.; Wang, X.−T.; Chen, J.−Q.; Muhammad, I.; Wang, S.; Liu, D.; et al. Comparative proteomic analysis provides insight into the key proteins as possible targets involved in aspirin inhibiting biofilm formation in Staphylococcus xylosus. Front. Microbiol. 2017, 8, 543. [Google Scholar] [CrossRef] [Green Version]
  61. Patel, K.M.; Teran, D.; Zheng, S.; Kandale, A.; Garcia, M.; Lv, Y.; Schembri, M.A.; McGeary, R.P.; Schenk, G.; Guddat, L.W. Crystal structures of Staphylococcus aureus ketol−acid reductoisomerase in complex with two transition state analogues that have biocidal activity. Chem. Eur. J. 2017, 23, 18289–18295. [Google Scholar] [CrossRef]
  62. Troccaz, M.; Benattia, F.; Borchard, G.; Clark, A.J. Properties of recombinant Staphylococcus haemolyticus cystathionine ß−lyase (metC) and its potential role in the generation of volatile thiols in axillary malodor. Chem. Biodivers. 2008, 5, 2372–2385. [Google Scholar] [CrossRef]
  63. Alreshidi, M.M.; Dunstan, R.H.; Macdonald, M.M.; Smith, N.D.; Gottfries, J.; Roberts, T.K. Metabolomic and proteomic responses of Staphylococcus aureus to prolonged cold stress. J. Proteom. 2015, 121, 44–55. [Google Scholar] [CrossRef]
  64. Murphy, G.R.; Dunstan, R.H.; Macdonald, M.M.; Gottfries, J.; Roberts, T.K. Alterations in amino acid metabolism during growth by Staphylococcus aureus following exposure to H2O2− a multifactorial approach. Helyon 2018, 4, e00620. [Google Scholar] [CrossRef] [Green Version]
  65. Alreshidi, M.M.; Dunstan, R.H.; Macdonald, M.M.; Smith, N.D.; Gottfries, J.; Roberts, T.K. Amino acids and proteomic acclimation of Staphylococcus aureus when incubated in a defined minimal medium supplemented with 5% sodium chloride. Microbiologyopen 2019, 8, e772. [Google Scholar] [CrossRef] [Green Version]
  66. Mongodin, E.; Finan, J.; Climo, M.W.; Rosato, A.; Gill, S.; Archer, G.L. Microarray transcription analysis of clinical Staphylococcus aureus isolates resistant to vancomycin. J. Bacteriol. 2003, 185, 4638–4643. [Google Scholar] [CrossRef] [Green Version]
  67. Schelli, K.; Zhong, F.; Zhu, J. Comparative metabolomics revealing Staphylococcus aureus metabolic response to different antibiotics. Microb. Biotechnol. 2017, 10, 1764–1774. [Google Scholar] [CrossRef]
  68. Matsuo, M.; Hiramatsu, M.; Singh, M.; Sasaki, T.; Hishinuma, T.; Yamamoto, N.; Morimoto, Y.; Kirikae, T.; Hiramatsu, K. Genetic and transcriptomic analyses of ciprofloxacin−tolerant Staphylococcus aureus isolated by the replica plating tolerance isolation system (REPTIS). Antimicrob. Agents Chemother. 2019, 63, e02019-18. [Google Scholar]
  69. Kim, J.; Senadheera, D.B.; Lévesque, C.M.; Cvitkovich, D.G. TcyR regulated L−cystine uptake via the TcyABC transporter in Streptococcus mutans. FEMS Microbiol. Lett. 2012, 328, 114–121. [Google Scholar] [CrossRef] [Green Version]
  70. Gill, S.R.; Fouts, D.E.; Archer, G.L.; Mongodin, E.F.; DeBoy, R.T.; Ravel, J.; Paulsen, I.T.; Kolonay, J.F.; Brinkac, L.; Beanan, M.; et al. Insights of evolution of virulence and resistance from the complete genome analysis of an early methicillin−resistant Staphylococcus aureus strain and a biofilm−producing methicillin−resistant Staphylococcus epidermidis strain. J. Bacteriol. 2005, 187, 2426–2438. [Google Scholar] [CrossRef] [Green Version]
  71. Schaedler, T.A.; Tong, Z.; van Veen, H.W. The multidrug transporter LmrP protein mediates selective calcium efflux. J. Biol. Chem. 2012, 287, 27682–27690. [Google Scholar] [CrossRef] [Green Version]
  72. Huet, A.A.; Raygada, J.L.; Mendiratta, K.; Seo, S.M.; Kaatz, G.W. Multidrug efflux pump overexpression in Staphylococcus aureus after single and multiple in vitro exposures to biocides and dyes. Microbiology 2008, 154, 3144–3153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Santos Costa, S.; Viveiros, M.; Amaral, L.; Couto, I. Multidrug efflux pumps in Staphylococcus aureus: An update. Open Microbiol. J. 2013, 7, 59–71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Wales, A.D.; Davies, R.H. Co−selection of resistance to antibiotics, biocides and heavy metals, and its relevance to foodborne pathogens. Antibiotics 2015, 4, 567–604. [Google Scholar] [CrossRef] [Green Version]
  75. Hassan, K.A.; Skurray, R.A.; Brown, M.H. Active export proteins mediating drug resistance in Staphylococci. J. Mol. Microbiol. Biotechnol. 2007, 12, 180–196. [Google Scholar] [CrossRef]
  76. Jennings, M.C.; Forman, M.E.; Duggan, S.M.; Minbiole, K.P.C.; Wuest, W.M. Efflux pumps may not be the major drivers of QAC resistance in methicillin−resistant Staphylococcus aureus. Chembiochem 2017, 18, 1573–1577. [Google Scholar] [CrossRef]
  77. Somerville, G.A.; Prostor, R.A. At the crossroads of bacterial metabolism and virulence factor synthesis in Staphylococci. Microbiol. Mol. Biol. Rev. 2009, 73, 233–248. [Google Scholar] [CrossRef] [Green Version]
  78. Chandrangsu, P.; Van Loi, V.; Antelmann, H.; Helmann, J.D. The role of bacillithiol in gram−positive Firmicutes. Antioxid. Redox Signal. 2018, 28, 445–462. [Google Scholar] [CrossRef] [Green Version]
  79. Pöther, D.−C.; Liebeke, M.; Hochgräfe, F.; Antelmann, H.; Becher, D.; Lalk, M.; Lindequist, U.; Borovok, I.; Cohen, G.; Aharonowitz, Y.; et al. Diamide triggers mainly S thiolation in the cytoplasmic proteomes of Bacillus subtilis and Staphylococcus aureus. J. Bacteriol. 2009, 191, 7520–7530. [Google Scholar]
  80. Soutourina, O.; Dubrac, S.; Poupel, O.; Msadek, T.; Martin−Verstraete, I. The pleiotropic CymR regulator of Staphylococcus aureus plays an important role in virulence and stress response. PLoS Pathog. 2010, 6, e1000894. [Google Scholar] [CrossRef] [Green Version]
  81. Liebeke, M.; Dörries, K.; Zühlke, D.; Bernhardt, J.; Fuchs, S.; Pané−Farré, J.; Engelmann, S.; Völker, U.; Bode, R.; Dandekar, T.; et al. A metabolomics and proteomics study of the adaptation of Staphylococcus aureus to glucose starvation. Mol. Biosyst. 2011, 7, 1241–1253. [Google Scholar] [CrossRef]
  82. Forbes, S.; Dobson, C.B.; Humphreys, G.J.; McBain, A.J. Transient and sustained bacterial adaptation following repeated sublethal exposure to microbicides and a novel human antimicrobial peptide. Antimicrob. Agents Chemother. 2014, 58, 5809–5817. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Gene ontology analysis of the significant differentially expressed genes (FC >5, p <0.05) of SE11Ad (DDAC-adapted S. epidermidis strain) for ‘Biological processes’ (A), ‘Cellular component’ (B) and ‘Molecular function’ (C).
Figure 1. Gene ontology analysis of the significant differentially expressed genes (FC >5, p <0.05) of SE11Ad (DDAC-adapted S. epidermidis strain) for ‘Biological processes’ (A), ‘Cellular component’ (B) and ‘Molecular function’ (C).
Microorganisms 08 00344 g001
Figure 2. Gene ontology analysis of the significant differentially expressed genes (FC >5, p <0.05) of SE18To (DDAC-tolerant S. epidermidis strain) for ‘Biological processes’ (A), ‘Cellular component’ (B) and ‘Molecular function’ (C).
Figure 2. Gene ontology analysis of the significant differentially expressed genes (FC >5, p <0.05) of SE18To (DDAC-tolerant S. epidermidis strain) for ‘Biological processes’ (A), ‘Cellular component’ (B) and ‘Molecular function’ (C).
Microorganisms 08 00344 g002
Table 1. RNA sequencing statistics.
Table 1. RNA sequencing statistics.
S. epidermidisConditionLibrary
Strain N°. of ReadsData Yield (bp)Average Length (bp)Mapped to Reference Genome (%)
SE11Control6,748,314796,550,256118.0485.46
SE11Ad Adapted6,939,6981,284,519,656185.0196.02
SE18Control13,150,2011,826,332,974138.8891.12
SE18To Tolerant15,425,0271,945,548,307126.1395.63
Table 2. The most significant differentially expressed genes (FC >5, p <0.05) in SE11Ad, the didecyldimethylammonium chloride (DDAC)-adapted S. epidermidis strain, when exposed to the adapted concentration of 18.0 µg/mL DDAC, according to the KEGG pathways.
Table 2. The most significant differentially expressed genes (FC >5, p <0.05) in SE11Ad, the didecyldimethylammonium chloride (DDAC)-adapted S. epidermidis strain, when exposed to the adapted concentration of 18.0 µg/mL DDAC, according to the KEGG pathways.
Gene ExpressionKEGG PathwayGene IDGeneFold ChangeProduct
Down-TransportersSE0405 −767.0Lipoprotein; ABC transporter
regulation SE0407 −250.9ABC transporter ATP-binding protein
SE0910 −5.6Signal recognition particle
Quorum sensingSE1543 −534.6V8 protease
SE1636agrD−391.5AgrD protein
SE0484psm β1/β2−285.1Phenol soluble modulin β1/ β2
SE1638agrA−81.3Accessory gene regulator A
SE1637agrC−77.0Accessory gene regulator C
SE0910 −5.6Signal recognition particle
Metabolism and SE0185gehD−495.5Glycerol ester hydrolase
biosynthesis ofSE1599 −26.1Chorismate mutase
secondary SE2288ipk−21.74-Diphosphocytidyl-2C-methyl-D-erythritol kinase
metabolitesSE0857mraY−15.1Phospho-N-acetylmuramoyl-pentapeptide transferase
SE2081 −14.3Beta-subunit of L-serine dehydratase
SE2311gltD−14.1Glutamate synthase subunit β
SE1345hemC−12.8Porphobilinogen deaminase
SE0349 −11.9Phosphomethylpyrimidine kinase
SE1199 −11.1Dihydrolipoamide dehydrogenase
SE0981miaA−8.9tRNA d(2)-isopentenylpyrophosphate transferase
SE0514nrdF−8.8Ribonucleotide-diphosphate reductase subunit β
SE1305 −8.7Iron–sulphur cofactor synthesis protein-like protein
SE0814 −8.5Heme synthase
SE1889murQ−7.6N-acetylmuramic acid-6-phosphate etherase
SE1596nadE−7.1NAD synthetase
SE1224 −6.2Shikimate kinase
SE1344 −5.9Uroporphyrinogen III synthase
SE1153 −5.83-Phosphoshikimate 1-carboxyvinyltransferase
Up-Transporters and SE1070pstS124.9Thioredoxin reductase
regulationefflux pumpSE1069pstC38.5Phosphate ABC transporter
SE2321metI19.6ABC transporter permease
SE0334 11.8Arsenic efflux pump protein
SE1992tcyB10.5ABC transporter permease
SE0681oppC8.1Oligopeptide transport system permease OppC
SE1067pstB8.0Phosphate transporter ATP-binding protein
SE2322metN5.4ABC transporter ATP-binding protein
Two-component systemSE0164uhpT11.5Sugar phosphate antiporter
Metabolism and SE0228arcC20.7Carbamate kinase
biosynthesis ofSE2175sat17.3Sulphate adenylyltransferase
secondarySE2197 13.4Alkaline phosphatase
metabolitesSE0319ispD9.02-C-methyl-D-erythritol 4-phosphate cytidylyltransferase
SE1654 8.0Dihydroxy-acid dehydratase
SE2174 7.1Adenylyl-sulphate kinase
SE2333 7.0NADH dehydrogenase subunit 5
SE2323 6.5Cystathionine gamma-synthase
Drug resistanceSE0681oppC8.1Oligopeptide transport system permease OppC
Table 3. The most significant differentially expressed genes (FC >5, p <0.05) in SE18To, the didecyldimethylammonium chloride (DDAC)-tolerant S. epidermidis strain, when exposed to the adapted concentration of 9.0 µg/mL DDAC, according to the KEGG pathways.
Table 3. The most significant differentially expressed genes (FC >5, p <0.05) in SE18To, the didecyldimethylammonium chloride (DDAC)-tolerant S. epidermidis strain, when exposed to the adapted concentration of 9.0 µg/mL DDAC, according to the KEGG pathways.
Gene ExpressionKEGG PathwayGene IDGeneFold ChangeProduct
Down−Metabolism and SE2220dapE−34.0Succinyl−diaminopimelate desuccinylase
regulationbiosynthesis ofSE1079 −16.9Alanine racemase
secondarySE0181 −16.38−Amino−7−oxononanoate synthase
metabolitesSE2600 −8.8Thiazole synthase
SE0486 −8.56−Pyruvoyl tetrahydrobiopterin synthase
SE1657 −8.0Ketol−acid reductoisomerase
SE2380metC−7.8Cystathionine beta−lyase
SE0275hisA−7.0Imidazole−4−carboxamide isomerase1−(5−phosphoribosyl)−5−[(5− phosphoribosylamino)methylideneamino] imidazole−4−carboxamide isomerase
TransportersSE2400 −9.9vraD protein
SE2255 −7.7Transport system protein
Up−TransportersSE1992tcyB16.0ABC transporter permease
regulation SE2248secY−211.8Accessory Sec system protein translocase subunit SecY2
SE1242znuB6.8ABC transporter
SE0871lspA5.2Lipoprotein signal peptidase
Metabolism and SE0274 8.0Imidazole glycerol phosphate synthase subunit HisH
biosynthesis ofSE0510 8.07−Cyano−7−deazaguanine reductase
secondary metabolitesSE0873 7.0Bifunctional pyrimidine regulatory protein PyrR uracil phosphoribosyltransferase
RibosomeSE1017rpmG25.350S ribosomal protein L33 2
SE0913rpsP5.630S ribosomal protein S16
Drug resistanceSE1606 5.5Penicillinase repressor

Share and Cite

MDPI and ACS Style

Ribič, U.; Jakše, J.; Toplak, N.; Koren, S.; Kovač, M.; Klančnik, A.; Jeršek, B. Transporters and Efflux Pumps Are the Main Mechanisms Involved in Staphylococcus epidermidis Adaptation and Tolerance to Didecyldimethylammonium Chloride. Microorganisms 2020, 8, 344. https://doi.org/10.3390/microorganisms8030344

AMA Style

Ribič U, Jakše J, Toplak N, Koren S, Kovač M, Klančnik A, Jeršek B. Transporters and Efflux Pumps Are the Main Mechanisms Involved in Staphylococcus epidermidis Adaptation and Tolerance to Didecyldimethylammonium Chloride. Microorganisms. 2020; 8(3):344. https://doi.org/10.3390/microorganisms8030344

Chicago/Turabian Style

Ribič, Urška, Jernej Jakše, Nataša Toplak, Simon Koren, Minka Kovač, Anja Klančnik, and Barbara Jeršek. 2020. "Transporters and Efflux Pumps Are the Main Mechanisms Involved in Staphylococcus epidermidis Adaptation and Tolerance to Didecyldimethylammonium Chloride" Microorganisms 8, no. 3: 344. https://doi.org/10.3390/microorganisms8030344

APA Style

Ribič, U., Jakše, J., Toplak, N., Koren, S., Kovač, M., Klančnik, A., & Jeršek, B. (2020). Transporters and Efflux Pumps Are the Main Mechanisms Involved in Staphylococcus epidermidis Adaptation and Tolerance to Didecyldimethylammonium Chloride. Microorganisms, 8(3), 344. https://doi.org/10.3390/microorganisms8030344

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

Article Metrics

Back to TopTop