1. Introduction
Staphylococcus aureus has shown a high capacity for developing resistance to each newer anti- staphylococcal antibiotic developed. Antibiotic resistance in
S.
aureus results in higher morbidity, mortality, length of hospitalization and health expenditure [
1]. The cyclic anionic lipopeptide antibiotic daptomycin (DAP) [
2] is becoming a main resource in therapy against multidrug-resistant
S.
aureus, especially methicillin-resistant
S. aureus (MRSA), due to its activity against glycopeptide and linezolid-resistant MRSA [
3,
4]. The mechanism of action of DAP differs from conventional glycopeptides (vancomycin, teicoplanin), and is more like cationic antimicrobial peptides produced by the immune system. Vancomycin and teicoplanin impair the bacterial cell wall synthesis. They bind to the terminal dipeptide D-alanyl-D-alanine (D-Ala-D-Ala), within the disaccharide pentapeptide attached to lipid II and the forming PGN chain. This way, transglycosylation and transpeptidation are prevented, thus hindering proper wall synthesis [
5,
6]. In addition, vancomycin has other secondary actions, probably less important in the lysis of the microorganism, such as inhibition of RNA synthesis and alteration of wall permeability.
The action of DAP is based on binding to the bacterial cell membrane, through its hydrophobic end, in the presence of calcium ions [
7,
8,
9]. This binding occurs both in exponential growth phase and stationary phase, causing membrane depolarization due to the loss of potassium ions from the cytoplasm [
10]. This process leads to the disruption of multiple functions of the bacterial cell membrane, such as synthesis of proteins, DNA and RNA without the need for the antimicrobial to penetrate the cytoplasm. Unlike other antimicrobials, the mechanism of action does not involve lysis of bacterial cells [
11] but rather bacterial cell death due to alterations in cellular homeostasis [
7,
12].
Intrinsic resistance to DAP by hydrolytic cleavage of ester bonds inside the DAP molecule has been described, in absence of contact with this or other antibiotics, among low percentage of GC gram-positive bacteria, such as staphylococci [
13].
Streptomyces strains highly resistant to daptomycin (MIC≥256 µg/ml) have been found to inactivate the antibiotic through two primary modes of daptomycin inactivation: ring hydrolysis, leading to linearization of the cyclic compound and deacylation of the lipid tail [
14]. Additionally,
Paenibacillus lautus has been shown to inactivate daptomycin. The production of the inactivating activity in
P. lautus was inducible by exposure to daptomycin and some data suggest might be caused by metalloesterases (14).
This “natural” resistance mechanism has been suggested to be associated with natural production of DAP-like molecules by different microorganisms. Though this type of resistance has not been shown in clinically important bacteria yet, resistance genes coding for these DAP inactivating enzymes might be a potential source of resistance if they could be captured by pathogenic bacteria, as has happened before for other resistance determinants, such as CTX-M enzymes, originated in
Kluyvera spp. and now present in most clinically important enterobacteria [
15].
Nevertheless, DAP-non susceptibility observed in isolates obtained both
in vitro and from therapeutic failures in
S. aureus, seems to be mainly associated with electrostatic repulsion of the DAP–calcium complex from the cell surface, due to an increase in the positive charge of the bacterial surface [
16].
DAP-non susceptible isolates obtained both
in vitro and from therapeutic failures, usually harbor mutations in genes associated with cytoplasmic membrane (
mprF) and cell wall homeostasis (
yycFG, also known as
walK), and mutations in RNA polymerase subunits
rpoB/
rpoC [
7,
8]. DAP-non susceptible isolates have been described, both
in vitro and
in vivo, lacking mutations in some of these genes, but not in all of them [
7,
9,
10,
11]. Other authors have shown that DAP-non susceptibility may also involve the upregulation of genes associated with cell wall synthesis and turnover, such as the two-component regulator
vraSR [
12]. Otherwise, whole genome sequencing studies have shown that DAP-non susceptible isolates may harbor mutations in other genes, such as
agrA,
pnpA,
pgsA,
prs,
cls2 and
clpP, some of them associated with the cytoplasmic membrane and cell wall metabolism [
10,
11].
Therefore, the present study had the following aims: identification of S. aureus strains (MRSA and MSSA) non-susceptible to DAP isolated from 3 different hospitals in Spain; test sensitivity to different drugs such as oxacillin, vancomycin, linezolid, dalbavancin and telavancin; amplification and sequencing of the main genes associated with DAP resistance.
2. Materials and Methods
2.1. Selection of bacterial strains
Nine S. aureus clinical isolates identified as DAP non-susceptible S. aureus (MIC >1 mg/L) by microdilution conventional susceptibility systems (WalkAway, Perkin Elmer, USA) were included in this study. Seven of them had been isolated in hospitals in the province of Badajoz (Southwest of Spain), and two in the Hospital Universitario de Salamanca (Middle West of Spain). All isolates except MSa9 were from patients who had received daptomycin treatment, although data on duration of treatment and dosage were not accessible. To our knowledge, the patient from whom the MSa9 isolate was obtained did not receive treatment with daptomycin at any time.
Microorganisms were identified by MALDI-TOF mass spectrometry (MicroflexTM, Bruker Daltonics, Germany). All the isolates used were preserved and stored at -80°C in glycerol broth with 10% skimmed milk.
2.2. Antibiotic susceptibility
Susceptibility to DAP, vancomycin (VAN), linezolid (LZD) and oxacillin (OXA) were tested by using a broth microdilution method (MicroScan WalkAway, Perkin Elmer, USA), according manufacturer instructions. Telavancin (TLV) and dalbavancin (DLV) were studied by E-test (Liofilchem, Italy). The microdilution sensitivity study was carried out using an automated commercial method (MicroScan WalkA-way, Perkin Elmer, USA), in which the antimicrobials are lyophilised so that a predetermined concentration is reached by adding the set amount of culture medium. The microorganism is inoculated in cation-adjusted Mueller Hinton broth, supplemented with 2% NaCl for oxacillin sensitivity and 50 microg/ml calcium for DAP. An inoculum of the microorganism equivalent to 0.5 of the McFarland (McF) standard is inoculated and incubated at 37°C for 20-24 hours. The reading is done in an automated way in a MicroScan equipment.
For the E-test study, a Mueller Hinton agar plate was inoculated with an inoculum equivalent to 0.5 McF, strips of the different antimicrobials are deposited on the inoculated agar, incubated at 37°C for 24 hours and read visually. For the DAP E-tests it is not necessary to supplement the agar with Ca2+, as the DAP strips have it incorporated.
Methicillin resistance, in strains with e-test sensitive to oxacillin, was verified by PCR specific for
mecA and
mecC genes (the primers used appear in
Table 1). The isolates were subjected to the E-test according to the manufacturer's instructions. Susceptibility criteria for VAN, LZD, DLV, TLV and OXA were those proposed by the CLSI (VAN: < 2 mg/L, susceptible; 4-8 mg/L, intermediate; > 16 mg/L, resistant. LZD: < 4 mg/L, susceptible; > 8 mg/L, resistant. DLV: < 0.25 mg/L, susceptible; intermediate and resistant categories not defined. TLV: < 0.12 mg/L, susceptible; intermediate and resistant categories not defined. OXA: < 0.5 mg/L, susceptible; > 1 mg/L, resistant) [
13]. The only categories considered for DAP were "susceptible" and "non-susceptible", as CLSI reports only the sensitivity criterion (MIC < 1 mg/L).
2.3. DNA extraction and amplification of genes associated with DAP non-susceptibility in S. aureus
A suspension with turbidity equivalent to 1 McF was obtained from colonies grown for 24 hours, at 37 ºC on blood agar. Then, 500 μl of the suspension were transferred to the automated nucleic acid extraction system NucliSENS® easyMAG (bioMérieux, France). Nucleic acids extraction was performed according to manufacturer instructions. DNA was preserved at -80°C.A series of genes previously associated with DAP non-susceptibility [
7,
14] were amplified and sequenced. Amplification was carried out in an Eppendorf Master Cycler® thermal cycler (Eppendorf, Germany). All primers used in the study were designed specifically (
Table 1), based on the
S. aureus ATCC 25923 complete genome sequence deposited in GenBank (gi: 685631213) using the Primer Blast® application (NCBI, USA). The primers were manufactured by Sigma-Genosys (USA), and the MasterMix® PCR amplification mix (Promega, USA9, a premixed 2X solution of Taq DNA Polymerase, dNTPs and reaction buffer, was used for PCR reaction. PCR conditions were as follows: 95°C x 5 minutes (denaturation); 35 cycles of: 95°C x 30 seconds, 1 minute at the temperature reported in supplementary section (
Table 1S) specific for each primer (hybridization), 72°C x 1 minute (elongation); and finally, 72°C x 5 minutes (elongation).
2.4. Visualization of amplicons and purification of PCR product
The products were subjected to electrophoresis in 2% agarose gel together with a size control (DNA Ladder 100 bp, Promega, USA). Blue/Orange 6x Loading Dye (Promega, USA) was used as loading buffer, and RedSafe® Nucleic Acid Staining Solution (Invitrogen, USA) was used as a fluorescent dye for DNA detection.
As a general rule, amplified DNA was cleaned for sequencing by using Exo-SAP It® (Chang Biosciences, China). When DNA was purified from the gel, the MEGAquick-Spin Total Fragment DNA Purification kit (Intron Biotechnology, South Korea) was used.
Amplicons were sequenced by the Sanger method in a 3500 Genetic Analyzer sequencer (Applied Biosystems, USA) and visualized by using the software Chromas 2.5.1.0. Sequences were compared with the S. aureus ATCC 25923 complete genome sequence deposited in GenBank (gi: 685631213). S. aureus ATCC 25923 is commonly used as a control in standard laboratory tests. Comparison and analysis were carried out by using the BLAST® software, freely available in the NIH website.
4. Discussion
DAP resistance seems to be based on complex mechanisms and has been associated with the presence of mutations in different genes, most of them involved in the metabolism and homeostasis of the cell membrane, and, in some cases, of the bacterial wall [
7,
14,
15]. A cause-effect relationship between specific mutations in specific genes and a specific increase in the MICs of DAP has hardly been described, and the repulsion hypothesis does not explain resistance to DAP in all
S. aureus isolates. Only some studies have shown that the replacement of the mutated
mprF gene with a wild one reverses, at least partially, the MIC of DAP [
16,
17]. In the present study, the MIC interval of DAP ranged between 2 and 4 mg/L in all cases, except MSa2, MIC of 1 mg/L. Although non-sensitivity to daptomycin in
S. aureus is MIC > 1 mg/L, due to its proximity to the range, MSa2 was not ruled out from the study.
Mutations in
rpoB and
mprF genes have been reported in most staphylococcal clinical isolates showing DAP MICs increase [
7,
10,
11].
mprF codes for the bifunctional enzyme MprF, that contributes to the positive charge of the cell membrane surface trough the lysinylation of phosphatidylglycerol (PG) and the translocation of lysinylated PG (L-PG) from the inner to the outer leaflet of the cell membrane. The positively charged cell membrane surface resulting from this process, helps repel the DAP molecule from the surface [
18,
19,
20]. The association between
mprF mutations and DAP non-susceptibility is supported by the finding that the inactivation of the mutated
mprF can reverse, at least partially, DAP non-susceptibility [
21].
MprF protein is composed of 14 transmembrane domains and a cytosolic C-terminal domain [
11]. The first eight N-terminal transmembrane domains are involved in translocation of L-PG to the outer CM. The next four transmembrane “central” are involved both in L-PG synthesis and flipping, while the cytosolic C-terminal domain is involved only in L-PG synthesis. Single nucleotide polymorphism (SNP) in
mprF genes have been reported in most staphylococcal clinical isolates showing DAP MICs increase [
7,
10,
11], and are associated with a gain of function. Nevertheless, some authors have shown that SNPs can also be found in 30% of DAP- susceptible isolates, and that some DAP-susceptible isolates with MICs of 1 mg/L can show up to 30 amino acid changes in
mprF [
22], suggesting a high level of heterogeneity within
mprF in some staphylococcal isolates. Nevertheless, the DAP-non susceptible isolates studied by these authors [
22], showed only one amino acid substitution in MprF (L341S or L826F). Changes in these or in very near positions have been previously reported associated with DAP-non susceptibility (S829L, S295L, P314L, S337L, T345I, T345A, I420N, [
7,
8,
11,
23,
24].
Unlike results reported by other authors, we found one SNPs in
mprF only in one isolate (MSa5). Mutations at the 314 position had been reported previously associated with DAP-non susceptibility, but as P314L [
7] instead of P314T, as appears in MSa5. P314L was found in other strain (MSa6,
Table 2). The other isolates had a high number of mutations affecting the three areas. They have between 6 and 8 mutations in the N-terminal “flippase” domains (mode: 6 mutations), 17-20 mutations (mode: 18 mutations) in the central, bifunctional domains, and 8- 9 mutations (mode: 8) in the phosphatidylglycerol lysinylation domain. Six isolates (MSa1, MSa4, MSa5, MSa6, MSa7 and MSa8) included in this wide group of mutations some amino acid changes in positions previously associated with DAP-non susceptibility (S295, P314, T345) [
7,
11,
16], but other isolates (MSa2, MSa3 and MSa9,) show mutations previously not reported in association with DAP- non susceptibility. Moreover, they show no mutations in the area comprised between S295 and I420, where most DAP-non susceptibility associated mutations reported are located, though they show more than 30 amino acid changes throughout the gene. The study reflects a large genetic variability between our clinical isolates, and the strain used as reference strain (
S. aureus ATCC 25923), especially in some of the genes studied. Some of the mutations observed in the different genes have been previously described associated with daptomycin resistance in previous studies [
7,
8,
11,
16,
23], but others have not. As already mentioned in the results section, the absence of MLST typing forces us to consider the possibility that, at least some of these previously undescribed changes are not associated with daptomycin resistance, but with the genetic variability existing between the different MLST types.
All the isolates showed
rpoB mutations. The
rpoB and
rpoC genes encode for bacterial RNA polymerase β and β’ subunits [
7,
8]. Mutations described in
rpoB and
rpoC associated with DAP- non susceptibility are different from those affecting to other antibiotics, such as rifampicin [
25]. Unlike
mprF mutations that, at least
in vitro, usually emerge at relatively early times during the selection process,
rpoB and
rpoC mutations emerge later [
7]. Some studies have reported
rpoB changes in mutants selected
in vitro, but not in clinical isolates [
7], while other studies [
11] did not detect
rpoB changes in strains selected
in vitro and detected them only in clinical isolates. Some authors [
25,
26] have shown that single point mutations in
rpoB, such as A477D or A621E, can reduce the susceptibility both to daptomycin and vancomycin. Such mutations have been shown to cause cell wall thickening and reduction of the negative charge of the outer layer in
S.
aureus. All the isolates tested in this study showed a
rpoB mutation (F737Y) that had not been previously described in non-susceptible isolates, whether clinical or obtained
in vitro. Changes in the
rpoC gene were also detected in all isolates, excepting MSa5. Single
rpoC mutations (N341D) have also associated with treatment failures in staphylococcal infections [
27]. Mutations in the walkR sequence, a regulator of cell wall metabolism and virulence in
S. aureus, B. subtilis and
S. pneumoniae [
28] have been associated with vancomycin resistance [
29], and also to reduced DAP susceptibility, in association or not with
mprF and
agrA mutations [
10,
29].
fakA had only been associated with DAP-non susceptibility in one study, and only in mutants selected
in vitro [
10]. In that study, authors find a L133H mutations associated with DAP-non susceptibility. In our study, this protein showed different changes in 8 out of 9 clinical isolates. Since the function of this protein is unknown, it is not possible to know its role in DAP susceptibility, or even if there is really a cause-effect relationship.
cls genes encode for a cardiolipin synthase, and mutations in these genes have been associated with DAP low susceptibility in staphylococci [
11]. Cls proteins are membrane-bound enzymes that synthesize cardiolipin, an important anionic membrane phospholipid, from the phosphatidyl moiety of two PG molecules [
30]. While in logarithmic growth PG is the main membrane phospholipid, in stationary growth conditions, and under conditions of stress, such as unfavorable growth conditions or cell-wall acting antibiotics, cardiolipin can accumulate up to 25%–30% of membrane phospholipid.
cls2 has been the most frequently reported
cls gene in this aspect. In this study, all isolates excepting MSa5 showed mutations in this gene. Most isolates shared two mutations (V135I, H205R), and two of them added a third mutation (A459L, A471E). None of these mutations had been described previously. Mutations previously reported appear between the amino acids 20 and 60 (A23V, T33N, L52F, F60S) [
11], in the two transmembrane domains described in this protein, while our mutations appear outside these domains and the cardiolipin synthase active domains.
Mutations in
cls1 have been described with a much lower frequency in association to DAP reduced susceptibility. Nevertheless, we found mutations in this gene in the same isolates in which we found mutations in
cls2, and the number of mutations in this gene/isolate ranged from 6 to 16. Most mutations found in the isolates tested in this study appear outside the transmembrane and cardiolipin synthase active domains, but V18A (in
cls1), that appears in all isolates excepting MSa5, and G20A, that appears in MSa6 and MSa9 (
Table 3), are in the first transmembrane domain, and I238V, that also appears in all isolates excepting MSa5, and I421M (isolates MSa1, MSa4 and MSa9) appear in cardiolipin synthase domains. Results concerning
cls1 and
cls2 suggest that, though DAP-non susceptibility associated mutations are more frequently described in
cls2, in the isolates tested in this study it seems more likely that the most transcendent mutations are those observed in
cls1 since, unlike
cls2, in almost all isolates we found mutations that specifically affect the transmembrane and cardiolipin synthase functional domains. Mutations in these areas may have a summation effect with those found in
mprF, since the negative charge provided by cardiolipin would be reduced or even suppressed, reducing the attraction between DAP and the bacterial membrane [
31].