Daptomycin Resistant Staphylococcus aureus Clinical Strain With Novel Non-synonymous Mutations in the mprF and vraS Genes: A New Insight Into Daptomycin Resistance

Objectives: Daptomycin (DAP) resistance in Staphylococcus aureus is uncommon but there are increasing reports of the emergence of resistance during DAP therapy. Most clinical DAP-resistant S. aureus isolates investigated carried mutations in the mprF gene. The aim of this study was to identify mutations between a clinical pair of methicillin-susceptible S. aureus (MSSA) isolates (DAP-susceptible and DAP-resistant). Additionally, the activity of genes previously associated with DAP resistance was assessed. Materials and Methods: Two MSSA isolates from patient with left-sided endocarditis were analyzed by whole genome sequencing (WGS) and reverse transcription-quantitative real-time PCR (RT-qPCR). The first isolate, DAP-susceptible, was obtained before initiation of treatment and the second isolate, DAP-resistant, was recovered after 4 weeks of DAP therapy. Results: Comparison of complete genomes of DAP-susceptible and its DAP-resistant variant identified two non-synonymous and one synonymous mutations. The non-synonymous mutations consisted of a S829L substitution in mprF and a T331I substitution in vraS. The RT-qPCR experiments revealed an increased expression of vraS, dltA, mprF, and sceD genes in DAP-resistant variant. Strikingly, the expression of dltA and mprF genes was significantly downregulated by DAP. Conclusion: The mprF and vraS genes were previously associated with DAP resistance, however, none of the mutations described in this study had been previously identified and linked to DAP resistance. Moreover, we provide a new insight into the DAP action on S. aureus, in which the expression of key genes in DAP resistance is decreased by the antibiotic.


INTRODUCTION
Daptomycin (DAP) is an alternative to vancomycin for invasive methicillin-resistant Staphylococcus aureus (MRSA) infections or for serious methicillin-susceptible S. aureus (MSSA) infections in patients who are allergic to beta-lactams. The DAP nonsusceptibility in S. aureus (referred to as DAP resistance in this study for the ease of presentation) is an increasing problem and several reports have described the emergence of resistance during DAP therapy (Lee et al., 2010;Mammina et al., 2010). The current knowledge suggests that DAP resistance in S. aureus is complex and results from mutational changes in a number of different genes. Most clinical DAP-resistant S. aureus isolates (MICs of >1 µg/ml) investigated to date, harbored mutations in mprF, typically in the form of single-nucleotide polymorphisms (SNPs) (Jones et al., 2008;Ernst et al., 2009;Mishra and Bayer, 2013). The mprF gene encodes a bifunctional membrane protein that catalyzes the synthesis and translocation (flipping) of the positively charged phospholipid lysyl-phosphatidylglycerol within its cell membrane. The amino acid substitutions in the MprF protein identified in the strains showing resistance to DAP lead to altered cell membrane phospholipid profiles. It results in a cell membrane positive charge increase and changes in cell membrane fluidity . The dltABCD operon is involved in the addition of D-alanine to teichoic acids in many Gram-positive bacteria (Ernst et al., 2009). Mutations in the dlt operon and/or altered expression of its genes lead to a cell surface positive charge increase, as in the case of the mprF mutations. Data from numerous studies have suggested that charge repulsion arising from the mprF-and dlt-mediated enhancement of positive surface charge is mainly responsible for DAP-resistant phenotype among S. aureus strains (Jones et al., 2008;Patel et al., 2011;Peleg et al., 2012;Gasch et al., 2013;Yang et al., 2013;Cafiso et al., 2014;Kang et al., 2017;Ma et al., 2018). It has recently been discovered that deletion of the clpX gene caused a small reduction in DAP susceptibility (Baek et al., 2014). The highly conserved ClpX chaperone facilitates protein folding and with the ClpP protease, forming the ClpXP protease, controls cell size and is required for growth of S. aureus at low temperature (Stahlhut et al., 2017). Other determinants involved in DAP resistance include genes that encode enzymes associated with phospholipid metabolism, such as phosphatidylglycerol and cardiolipin synthetases (pgsA and cls, respectively) (Fischer et al., 2011). There are also two-component regulatory systems like WalKR also known as YycFG, VraSR or GraRS that directly or indirectly modulate transcription of a number of genes encoding proteins involved in wall metabolism and permeability (Tran et al., 2015). The emergence of mutations in these regulatory genes have been also associated with DAP resistance in S. aureus (Friedman et al., 2006;Mehta et al., 2012).
The main purpose of this study was to apply whole genome sequencing (WGS) to a clinical pair of MSSA isolates (DAP-susceptible and DAP-resistant) to detect genome-wide DNA sequence polymorphisms associated with DAP resistance. Additionally, the gene expression profiles of determinants previously linked with DAP resistance were investigated. Furthermore, the net cell-surface charge as the main mechanism responsible for the DAP-resistant phenotype in MSSA was assessed.

Isolates
The isolates were obtained from a 50-year-old male, with a history of alcohol abuse and several comorbidities, who had a mitral valve replacement and in the following 6 months experienced two episodes of left-sided endocarditis, diagnosed according to current guidelines (Habib et al., 2009). In the last episode, Enterococcus faecalis and MSSA (isolate IT1-S) were obtained from blood cultures. As the patient had previously shown an allergic reaction to penicillins, treatment was instituted with gentamycin (discontinued after 2 weeks) and DAP at the dose of 500 mg daily. After 4 weeks the patient's conditions remained serious, with no improvement. Blood culture yield an MSSA (isolate IT4-R) that was resistant to DAP (Table 1).

Next-Generation Sequencing (NGS)
Next-generation sequencing (NGS) was carried out on the Illumina MiSeq system using 300 × 2 paired-end reads. Sequence libraries were created using the standard Illumina Nextera XT library creation kit. The SeqMan NGen software version 11.2.1 (DNASTAR) was used for de novo assembly of the reads and the resulting contigs were ordered by Mauve Contig Mover. The remaining gaps between contigs were closed by PCR amplification and Sanger sequencing allowing for the analysis of fully closed chromosomes and plasmids. Manual sequence editing was conducted using the SeqBuilder software (DNASTAR). The DNA sequences were aligned using the MegAlign (DNASTAR) and BLASTn software.
spa Typing and MLST The procedure was conducted as previously described (Aires-de-Sousa et al., 2006). The spa types were assigned using Ridom StaphType software version 1.4.6 (Ridom GmbH, Würzburg, Germany). The MLST STs were assigned through the publicly available MLST server 2 on the basis of WGS data. Each isolate was grown in triplicate in Mueller Hinton Broth (MHB) in the absence or presence of DAP. The overnight cultures were diluted to produce an inoculum concentration of approximately 10 8 cfu/ml. The inoculum was then diluted 1:100 in MHB supplemented with 50 µg/ml CaCl 2 or in MHB supplemented with 0.08 µg/ml DAP and 50 µg/ml CaCl 2 . It resulted in starting number of the cells approximately 10 6 cfu/ml in each culture. Cells were then harvested from the IT1-S and IT4-R cultures grown to mid exponential (∼10 8 cfu/ml, after 5.5 h and 4 h, respectively), late exponential (∼10 9 cfu/ml, after 8 h and 6 h, respectively) and stationary (∼10 10 cfu/ml, after 14 h and 12 h, respectively) phase. The cell concentration was confirmed by colony counting after plating and incubation onto blood agar plates. Total RNA was isolated using the RNeasy Mini kit (Qiagen) according to the manufacturer's instructions. RNA was quantified using Qubit (Thermo Fisher Scientific) and the quality of the RNA extracted was assessed by TapeStation 2200 (Agilent Technologies). The Agilent TapeStation 2200 system, which is an automated instrument for nucleic acid gel electrophoresis, assigns RNA Integrity Number (RIN) values ranging from 1 to 10, with 10 being the highest quality.
Only samples with preserved 16S and 23S peaks and RIN values >8 were selected for gene expression analyses. The RIN values >8 indicate intact, high quality RNA samples for downstream applications (Fleige and Pfaffl, 2006). Total RNA was further treated with Baseline-ZERO DNase (Epicentre) followed by the RNeasy MinElute Cleanup kit (Qiagen) according to the manufacturer's instructions. The absence of contaminating DNA was verified by PCR. For RT-qPCR analyses, an iTaq Universal SYBR Green One-Step Kit (Bio-Rad Laboratories) was used. Each reaction mix with a volume of 20 µl was prepared with 300 nM each primer (final concentration) and 20 ng of template RNA. A CFX96 Touch Real-Time PCR detection system (Bio-Rad Laboratories) was used for the measurements using a protocol with the following thermal cycling conditions: reverse transcription reaction at 50 • C for 10 min and polymerase activation and DNA denaturation at 95 • C for 1 min, followed by 35 cycles of denaturation at 95 • C for 10 s and annealing/extension at 60 • C for 30 s. After the last amplification cycle, a melting curve analysis was carried out by heating from 65 to 95 • C in increments of 0.5 • C/s. Negative controls (without template or reverse transcriptase enzyme) were included in each run. Gene-specific primers ( Table 2) for the genes used in RT-qPCR experiments were designed using SeqBuilder software (version 15.0.1) from DNASTAR. Calculations of primer efficiencies were performed using the software CFX Manager version 3.1 (Bio-Rad Laboratories). Fold changes in the expression levels of the investigated genes were normalized in relation to the levels of gyrB mRNA. The relative changes in gene expression were quantified using the Pfaffl method (Pfaffl, 2001): gene expression ratio = (E target ) Ct target(control − sample) /(E reference ) Ct reference(control − sample) , where E target is the amplification efficiency of target (gene of interest), E reference is the amplification efficiency of reference (gyrB), Ct is the point at which the fluorescence rises above the background fluorescence, Ct target is the Ct deviation of the control minus the sample of the target gene transcript, and Ct reference is the Ct deviation of the control minus the sample of the reference gene transcript. The statistical significance of differences in the results was analyzed with the Student t-test and P-values of ≤0.05 were considered significant.

Relative Positive Surface Charge Determination
The cells were grown overnight (16 h) in Tryptic Soy Broth (TSB) medium in the absence of DAP, washed twice with 20 mM morpholinepropanesulfonic acid (MOPS) buffer (pH 7.0), and resuspended in the same buffer at an optical density at 600 nm (OD 600 ) of 1.0. The cells were incubated with cytochrome c (Sigma) dissolved in 20 mM MOPS buffer (pH 7.0) using three different concentrations: 50 µg/ml, 100 µg/ml, and 500 µg/ml. Moreover, the cells incubated in 20 mM MOPS buffer (pH 7.0) without cytochrome c were used as a negative control. After 15 min of incubation, the cells were pelleted by centrifugation and absorbance of the supernatants was measured at 530 nm. The amount of remaining (unbound) cytochrome c in the supernatant was determined by comparison to a standard curve (serially diluted cytochrome c). Cytochrome c is a cationic peptide and there is a direct correlation between the amount of unbound cytochrome c detected in the supernatant and the positive charge of the bacterial surface (Peschel et al., 1999). The data were converted and shown as the percent amount of bound cytochrome c. Independent runs were performed in triplicate on separate days. Changes in cytochrome c binding were compared by the Student t-test. A P-value of ≤0.05 was considered significant. Statistical analysis was performed using the GraphPad software.

Membrane Potential Assay
DAP-induced bacterial membrane depolarization was performed essentially as described previously (Silverman et al., 2003) by using the membrane potential-sensitive fluorescent dye 3,3dipropylthiacarbocyanine [DiSC3(5); Thermo Fisher Scientific] with the following modifications. Briefly, isolates were grown at 37 • C to early exponential phase (OD 600 , 0.2 to 0.3) in 50 ml of MHB. Cells were harvested by centrifugation, washed twice with 5 mM HEPES buffer (pH 7.2) supplemented with 50 µg/ml CaCl 2 and resuspended in the same buffer (containing 50 µg/ml CaCl 2 ) to an OD 600 of 0.2. The dye DiSC3(5) was added to cell suspension to make a final concentration of 0.18 µM. The cells were incubated with the dye DiSC3(5) for 15 min at room temperature. Then KCl was added (100 mM final concentration) to equilibrate the cytoplasmic and external potassium ions (K + ) concentrations. The 200-µl cell suspension aliquots were transferred to the wells of a white 96-well microplate. The desired DAP concentrations were subsequently added to the microplate wells and an increase in the fluorescence intensity at a wavelength of 670 nm was measured over 60 min using a Synergy HTX multimode microplate reader (BioTek). The fluorescence leakage (F L ) was determined using equation as previously described by Cheng et al. (2014): where F F was the fluorescence intensity of the cell suspension over 60 min of treatment with DAP, F I was the initial fluorescence intensity of the cell suspension with DAP, and F B was the fluorescence intensity of the blank (only cells and the dye). Data were normalized relative to the maximum fluorescence leakage (expressed as 100%). Results were shown as the mean of triplicate measurements. Changes in membrane depolarization were compared by the Student t-test. A P-value of ≤0.05 was considered significant. Statistical analysis was performed using the GraphPad software.

Nucleotide Sequence Accession Numbers
Chromosome and plasmid sequence data of the two S. aureus isolates were annotated using the NCBI Prokaryotic Genome Annotation Pipeline and deposited in GenBank 3 under accession numbers CP028468-CP028471.

Antibiotic Resistance
The first aim of our investigations was to determine the antibiotic susceptibility of the isolates recovered before and after treatment with DAP. The MIC values for selected antibiotics are presented in Table 1. Antibiotic susceptibility testing showed that isolate IT1-S (recovered before the DAP treatment) was susceptible to all the antibiotics tested with the exception of clindamycin that showed intermediate resistance (MIC = 0.5 µg/ml). MIC value for DAP was 0.5 µg/ml. Isolate IT4-R (recovered after the DAP treatment) was resistant to DAP (MIC = 2 µg/ml) and susceptible to all other antibiotics tested. Therefore, the mechanism conferring DAP resistance in IT4-R did not confer resistance to other antibiotics, including vancomycin, for which the two isolates showed the same MIC ( Table 1).

Molecular Characterization
This part of the study aimed at molecular typing of the isolates and comparing their virulence potential. The characteristics of the isolates are shown in Table 1. spa typing assigned spa type t223 to both isolates. From the WGS data, in silico MLST identified the isolates belonged to single sequence type (ST) 22. Moreover, the isolates were positive for toxic shock syndrome toxin-1, enterotoxin N, enterotoxin O and enterotoxin P. The full virulence profile of the IT1-S and IT4-R isolates is presented in Table 3. Both isolates showed exactly the same virulence potential as they carried the same set of the genes.

The S. aureus IT4-R Genome Carries Mutations That Can Be Linked to DAP Resistance
Subsequently, the genome sequences of the two isolates were compared to identify mutational changes responsible for the DAP resistance of the IT4-R isolate. The analyzed S. aureus isolates, IT1-S and IT4-R, had identical genomes in length, which included a 2,774,523-bp circular chromosome (32.8% G+C content) and a 34,104-bp circular plasmid (30.3% G+C content). The NCBI annotation pipeline revealed 2,895 genes in total in each genome. Whole-genome comparison of both isolates, DAPsusceptible and its DAP-resistant variant, identified three single nucleotide polymorphisms (SNPs), one in each of three genes. Two of the mutations were non-synonymous (leading to a change in the amino acid sequence in the translated gene product) while the third mutation was synonymous (without change in amino acid sequence). The non-synonymous mutations consisted of (i) a S829L substitution in mprF and (ii) a T331I substitution in vraS. The synonymous mutation was found in sceD (A→G at the nucleotide 105), the gene involved in cell wall turnover, growth and cell separation (Howden et al., 2010). Although mprF and vraS genes were previously associated with DAP resistance, none of the SNPs that we identified had been previously described and linked to DAP resistance.

Comparison of Gene Expression Between the IT1-S and IT-4 Isolates
We wanted to know if the mutations, which arose in genomic DNA during DAP therapy had an influence on gene expression. The gene expression levels of IT4-R were compared to that of its parent strain IT1-S during mid-exponential, late exponential and stationary growth phases using RT-qPCR (Figure 1). All three genes, in which the SNPs were identified in this study (mprF, sceD and vraS) and the genes previously linked with DAP resistance (dltA, walK, graR, cls, pgsA and clpX) were investigated. The vraS gene showed significantly increased expression (P value ≤ 0.05) in IT4-R relative to IT1-S in all growth phases in both the absence and presence of DAP (Figure 1 and Supplementary Table S1).
Likewise, dltA, mprF and sceD were upregulated in IT4-R but the differences between the two isolates were not always statistically significant (Supplementary Table S1). Other genes showed no difference in regulation between the two isolates. In conclusion, the increased expression of the vraS, dltA, mprF, and sceD genes was observed in the DAP-resistant variant compared to its parent DAP-susceptible strain.

Influence of DAP on Gene Expression
It was shown previously that concentrations of some antibiotics below the MIC were able to modulate the expression of virulenceassociated genes in S. aureus (Bernardo et al., 2004;Stevens et al., 2007;Turner and Sriskandan, 2015). However, DAP seemed to have no significant effects on virulence factor expression by S. aureus (Otto et al., 2013). In the current study, we wanted to explore if subinhibitory concentration of DAP could influence the expression of the genes linked with DAP resistance. Therefore, the expression levels of mprF, sceD, vraS, dltA, walK, graR, cls, pgsA and clpX were compared in IT1-S or in IT4-R in the absence and presence of DAP during different stages of growth (Figure 1). DAP significantly downregulated the expression of the dltA gene.
With the exception of the stationary phase in IT4-R, the mprF gene was also downregulated in the presence of DAP. However, the difference was not significant for the mid-exponential phase in IT1-S (Supplementary Table S2). Based on these observations we could conclude that DAP downregulates the expression of dltA and mprF genes.

The Effect of Growth Phase on Gene Expression
Throughout bacterial growth, the cell density, nutrient conditions, pH, and other factors are changing. Therefore, we analyzed the growth phase-dependent gene expression profiles (Figure 1). In both isolates in the absence or presence of DAP, the sceD gene showed the highest transcript level in the mid-exponential phase and during further growth its transcriptional activity was decreasing, while the highest and lowest transcriptional activity of cls was found in late exponential and mid-exponential growth phases, respectively (Figure 1). However, the differences between the growth phases were not always statistically significant (Supplementary Table S3). Moreover, dltA and graR were always the most active in the late exponential growth phase and clpX in the mid-exponential phase.
In case of dltA and graR the differences were always statistically significant. In conclusion, independently on mutational changes during DAP therapy as well as presence or absence of DAP in medium, strong impact of the growth phase on gene expression was observed in case of sceD and cls.

Relative Positive Surface Charge
The emergence of DAP-resistant S. aureus strains can occur by way of several different mechanisms involving the cell membrane and/or cell wall. Recent analysis of a panel of clinical DAPresistant MSSA isolates showed that enhancement of positive surface charge, reducing DAP binding, may be the main FIGURE 1 | Expression analysis of genes previously linked with DAP resistance. Transcript levels of the analyzed genes were determined by RT-qPCR in relation to gyrB expression. Values represent the mean of results obtained for three replicates independently grown cultures for each isolate. Error bars indicate the standard deviation of comparisons between three replicates. IT1-S -isolate IT1-S grown in the absence of DAP, IT1-S-DAP -isolate IT1-S grown in the presence of 0.08 µg/ml daptomycin, IT4-R -isolate grown in the absence of daptomycin, IT4-R-DAP -isolate IT1-S grown in the presence of 0.08 µg/ml daptomycin. Mid-log, mid-exponential growth phase; Late Exp, late exponential growth phase; Plateau, stationary growth phase. mechanism of DAP resistance among the MSSA strains (Kang et al., 2017). To test this hypothesis the cytochrome c binding analysis was performed. Cytochrome c binding assay revealed that DAP-resistant isolate IT4-R had significantly decreased positive surface charge compared to its DAP-susceptible parental isolate IT1-S (Figure 2).

DAP-Induced Membrane Depolarization
The target of DAP is the bacterial cell membrane (Silverman et al., 2003). Once inside the membrane, DAP aggregates (Pogliano et al., 2012) leading to rapid depolarization and loss of membrane potential, and thus to bacterial cell death. We compared depolarization of the cell membrane of the IT1-S and IT4-R isolates to identify differences, which could explain higher resistance of the mutant to DAP. Changes in membrane depolarization were determined by membrane potential assay using a fluorescence dye, DiSC3(5). The fluorescence of DiSC3(5) decreases as it incorporates into polarized membranes because at high concentration the dye aggregates and self-quenches. With the addition of a membrane-disrupting agent, such as DAP, the dye is released from cells into the media, which in turn leads to fluorescence dequenching. DAP concentrations that were used in the membrane potential assay were 1 µg/ml and 4 µg/ml that corresponded to twice the DAP MIC for IT1-S and IT4-R, respectively. The addition of DAP into the assay medium gradually depolarized the membrane of the isolates (Figure 3). The rate of membrane depolarization was reduced in DAPresistant isolate compared to that of DAP-susceptible isolate. However, statistically significant differences (p < 0.05) between IT1-S and IT4-R at the same time points were observed only between 10 and 40 min at the DAP concentration of 1 µg/ml, and between 10 and 30 min at the DAP concentration of 4 µg/ml. At the DAP concentration of 1 µg/ml, a maximum depolarization (100%) of the IT1-S membrane was observed after 50 min, while the antibiotic depolarized the membrane of IT4-R by 85.6% (compared to the maximum depolarization of IT1-S) after 60 min. At the DAP concentration of 4 µg/ml, the time to reach the maximum of depolarization was shorter for IT1-S (30 min), than for IT4-R (40 min). However, the peak of depolarization for IT4-R was lower than that of IT1-S (91.2% versus 100%). This assay demonstrated that both the total amount of membrane depolarization and the rate of depolarization were reduced in DAP-resistant isolate, which can explain higher resistance of the mutant to the antibiotic.

DISCUSSION
The MprF protein consists of an N-terminal transmembrane flippase domain, C-terminal catalytic synthase domain, and a central bifunctional domain bridging the flippase and synthase domains. The changes in amino acid sequence linked to DAP resistance of S. aureus were identified in all domains of MprF: G61V in the flippase domain (Peleg et al., 2012); S295L, P314L, S337L, T345A, and T345I in a central bifunctional domain (Friedman et al., 2006;Mehta et al., 2012;Peleg et al., 2012); and I420N and L826F in synthase domain (Peleg et al., 2012). Taking into account a frequency of above mutations they occur most often in a central bifunctional domain both in clinical and in vitro obtained isolates. Our study revealed a novel amino acid substitution in C-terminal catalytic synthase domain of the MprF protein at amino acid residue 829 (S829L), which was found in a clinical isolate with the DAP-resistant phenotype after prolonged treatment of a patient with endocarditis initially infected with a fully susceptible strain.
Although several molecular mechanisms of resistance to DAP have been proposed, mainly associated with phenotypic changes in cell wall and cell membrane, the basis of DAP resistance in S. aureus is still incompletely understood (Tran et al., 2015). This results from the fact that DAP resistance in S. aureus appears to be multi-factorial and strain specific. It has been proposed that charge repulsion of the DAP antibiotic molecule from the cell surface, which is associated with the increased positive charge of cell membrane and cell wall may be the main mechanism to explain DAP resistance in S. aureus strains (Jones et al., 2008;Patel et al., 2011;Peleg et al., 2012;Gasch et al., 2013;Yang et al., 2013;Cafiso et al., 2014;Kang et al., 2017;Ma et al., 2018). During the emergence of DAP resistance in S. aureus, the first SNPs in the genomic DNA appear most often in the mprF gene and are associated with a gain of enzymatic function, resulting in an increase in the positive charge of the cell membrane . The mutant in this study displayed an amino acid substitution in the MprF protein and increased expression of the mprF gene. These two factors possibly reduced dissipation of membrane potential of the mutant in the presence of DAP, promoting resistance to the antibiotic molecules by electrostatic repulsion.
Another strategy utilized by S. aureus to enhance the positive surface charge is increased expression of the dlt operon (Bertsche et al., 2013;Cafiso et al., 2014;Mishra et al., 2014). The dltABCD gene products are involved in the incorporation of D-alanine into cell wall teichoic acids. Although in the RT-qPCR experiments we observed increased expression of the dltA gene in IT4-R, this isolate had higher capacity to bind cytochrome c than IT1-S. However, the test with cytochrome c measures the overall cell surface charge and not the membrane charge itself. Moreover, it has been previously shown that some DAP-resistant isolates with mprF mutations did not show significantly altered net surface charge (Bayer et al., 2015). It will be interesting to compare more clinical DAP-resistant S. aureus strains with non-increased surface charge relative to their respective isogenic susceptible strains by DAP-induced cell membrane potential assay to explore further the mechanisms of DAP resistance in S. aureus.
Second non-synonymous mutation identified in the genome of IT4-R was found in vraS, a gene of the two-component system VraSR, which positively modulates the regulation of cell-wall biosynthesis pathway in S. aureus. Moreover, the expression levels of vraS were significantly increased in IT4-R compared to that of DAP-susceptible counterpart. In a study recently published by Chen et al. (2018), the authors showed results indicating a causal relationship between the L431F substitution in the MprF protein and increased expression of the vraSR genes in the DAPresistant strain leading concurrently to vancomycin resistance. But how the mutant MprF protein affected expression of vraSR and the cell wall-related genes was not elucidated. We can assume that in the IT4-R isolate the amino acid substitution in the MprF synthase domain (S829L) can lead not only to changes in the cell membrane of but also to pleiotropic effects, including upregulation of the VraSR system and the transglycosylase gene sceD, a feature involved in cell-wall turnover, resulting in increased resistance to DAP.
A previous study showed that the mprF sequence variations within the same clonal complex were not only associated with the DAP-resistant S. aureus strains but also with DAP-susceptible strains . There is an increasing interest to use WGS for genotypic prediction of antimicrobial resistance, which has the potential to reduce bed-to-diagnosis time by eliminating culture-based antimicrobial susceptibility testing. However, accuracy of predictions depends on a priori knowledge of mutations in mprF as well as in other genes that are conferring the DAP-resistant phenotype in S. aureus. Our study identified two novel SNPs leading to the amino acid substitutions (S829L in MprF and T331I in VraS), which can be linked to S. aureus DAP resistance, improving our recognition of resistance signatures within core genomic determinants.
This study provided a new insight into the DAP action on S. aureus by way of decreasing the expression of mprF and dltA. The dltABCD operon, like the mprF gene, contributes to the staphylococcal net positive surface charge. It will be important to identify the mechanism by which DAP reduces transcriptional activity of the genes, which contribute to positive surface charge since our study showed that the antibiotic did not influence the expression of two-component systems, VraSR, GraSR, and WalKR, associated with resistance of S. aureus to DAP.

CONCLUSION
We identified the new point mutation in the mprF gene leading to the amino acid substitution in the MprF protein that counteracts DAP antibiotic activity. Our results support the suggestion that vraSR contributes directly or indirectly to DAP resistance in S. aureus. Description of new mutations in the core genome, which can be linked with S. aureus resistance to DAP will allow development of better and more rapid diagnostic methods for identifying of drug resistance. The membrane potential assay demonstrated that both the total amount of membrane depolarization and the rate of depolarization were reduced in DAP-resistant isolate, possibly resulting in resistance to the antibiotic molecules by electrostatic repulsion. However, the reduced effect of DAP on membrane depolarization in DAPresistant isolate, compared to its DAP-susceptible parent, was not associated with increased positive net cell surface charge, which warrants further investigation of a higher number of such isolates to better understand the mechanisms of DAP resistance in S. aureus. Finally, we shed a new light on the DAP action in S. aureus, in which the expression of key genes in DAP resistance is decreased by the antibiotic.

AUTHOR CONTRIBUTIONS
AS, MT, HG, AP, and AF designed the project. MT and AP provided the isolates, clinical and epidemiological data. VA, MM, MDG, and GE performed the experiments. AS wrote the manuscript. All authors interpreted the data and reviewed the manuscript.

FUNDING
This project was financed in part by funds granted by the European Regional Development Fund within the EurHealth-1Health project (EU/INTERREG VA-681377 to AS, VA, and AF).