Transcriptomic Analysis of Staphylococcus aureus Under the Stress Condition Caused by Litsea cubeba L. Essential Oil via RNA Sequencing

Litsea cubeba L. essential oil (LCEO) is a natural essential oil with considerable antimicrobial activity, and it can gradually replace some chemical additives in the food industry. However, the genetic evidences of stress response of bacteria under sub-lethal treatment with LCEO is limited. To this end, transcriptomic analysis of Staphylococcus aureus 29213 under a low concentration of LCEO was performed. Bacterial RNA samples were extracted from 1/4 MIC (0.07 μL/mL) of LCEO-treated and non-treated S. aureus 29213. The transcriptional results were obtained by RNA sequencing (RNA-Seq). After treated with LCEO of S. aureus 29213, 300, and 242 genes were significantly up and down-regulated. Up-regulated genes were mainly related to cell membrane (wall) stress stimulon including genes related to two-component regulatory system (VraS), histidine metabolism (hisABCD etc.) and L-lysine biosynthesis (thrA, lysC, asd etc.). Significant differences were also founded between LCEO-treated and non-treated groups in peptidoglycan biosynthesis related pathways. Down-regulated genes were related to nitrogen metabolism (NarGHIJ etc.), carotenoid biosynthesis (all) and pyruvate metabolism (phdA, pflB, pdhC etc.) of S. aureus 29213 in an LCEO-existing environment compared to the control. At the same time, we confirmed that LCEO can significantly affect the staphyloxanthin level of S. aureus 29213 for the first time, which is closely related to the redox state of S. aureus 29213. These evidences expanded the knowledge of stress response of S. aureus 29213 strain under sub-lethal concentration of LCEO.


INTRODUCTION
Litsea cubeba L. (L. cubeba) is a plant of the Lauraceae family and distributed in many countries in East Asia, and it is mainly distributed in the Southern and western parts of China (Liu and Yang, 2012). All parts of the plant have a ginger-like smell (Ho et al., 2010). The fruits of L. cubeba are often used by people from Taiwan as a condiment for dishes (Liao et al., 2015). As an ethnobotanical material, L. cubeba has been used for headache, fatigue, chronic tracheitis and acute mastitis; the fresh leaves of it are grinded into a powder then used for skin problems, such as snake bite and furuncles (Chen et al., 2013;Nguyen et al., 2016). In the past decade, many studies focused on the pharmacological functions and the chemical compositions of L. cubeba essential oil (LCEO). The main component of L. cubeba in the fruit leaf, flower and twig oils is citral, 1,8-cineole, while in the stem oil there are limonene, citronellal, and citronellol (Ho et al., 2010). Research about the difference of main chemical composition in LCEO from various regions of China was also published (Si et al., 2012). The fruit oil of LC has been used for the flavor enhancer in cosmetics, cigarettes and food products (Wang et al., 2009;Liu and Yang, 2012). In some Chinese food factories, LCEO was used as food preservatives. Meanwhile, LCEO can be used as a material for citral (neral and geranial) refining in the market (Zhili et al., 2009). In pharmacological research, LCEO have been reported to have antioxidative (Hwang et al., 2005;She et al., 2019), anticancer [10,15] insecticidal (Seo et al., 2009) and antimicrobial activities (Wang and Liu, 2010;Liu and Yang, 2012).
Staphylococcus aureus (S. aureus) is a serious bacterial pathogen that causes infections (e.g., sepsis, meningitis, pneumonia, endocarditis, osteomyelitis, vomiting, nausea and abdominal cramping) and food poisoning in many animals and humans (Lowy, 1998;Argudin et al., 2010). Resent study has shown that the antibacterial ability of LCEO which was extracted in July on S. aureus was stronger than that of any other months (She et al., 2019). Research findings indicated the destructive effect of LCEO on the cell membrane of methicillinresistant S. aureus (MRSA) caused the leakage of intracellular macromolecules in cell . However, the deep impact of LCEO on S. aureus has not been well studied.
In recent years, various molecular tools and bioinformatics techniques have been widely used for identifying and analyzing how antibacterial agent interact with bacteria more accuracy and efficiency (Park et al., 2019). Thanks to the next generation sequencing, transcriptome analysis plays a significant role in understanding the mechanism of microbe which was treated by an antibacterial agent (Delpech et al., 2015). The aim of this study is to elucidate the stress response of S. aureus 29213 strain under sub-lethal concentration of LCEO at transcriptome level by RNA-Seq method.

Oil and Bacterial Strains
Litsea cubeba L. essential oil was came from a local agency (Nanning, China) and the chemical compositions of LCEO was shown in Table 1. The LCEO in present work was extracted from the fruit of LC. As essential oils do not dissolve in water, LCEO, tween-80 and distilled water (volume ratio of 1: 1: 5) was mixed to disperse the LCEO in water and the hydrophile-lipophile balance number of this system was 15.
S. aureus ATCC 29213 strain has the ability to form biofilm and to generate staphyloxanthin as well as it was often used in antimicrobial research (Musthafa and Voravuthikunchai, 2015;Sadrearhami et al., 2019;Singh et al., 2019). In this study, the S. aureus 29213 strain was provided by The First People's Hospital of Nanning, China.
One milliliter S. aureus 29213 was taken out from the refrigerator (−80 • C) then inoculated into tryptic soy broth (TSB) at 37 • C for 36 h with 250 RPM shake, then 10 µL cultured solution was incubated into TSB at 37 • C, 12 h with 250 RPM shake for viability recover. One microliter cultured solution was taken by inoculation hard loop to inoculate on a plate with tryptone soy agar (TSA) at 37 • C for 24 h. A single colony was incubated in a tube with 3.5 mL of TSB at 37 • C with shaking at 250 rpm for 8 h before subsequent trials.

Determination of Minimum Inhibitory Concentrations (MICs)
Minimum Inhibitory Concentration values of LCEO on S. aureus 29213 strain was determined by the serial-dilution culture method. TSB was used as the incubation medium. The LCEO solution was serially diluted by TSB. All tubes (13 × 100 mm) contained 1.75 mL LCEO solution and 1.75 mL of diluted S. aureus 29213 inoculum (approximately 10 4 CFU/mL), with a final concentration of LCEO solution of 71.4, 35.7, 17.8, 8.9, 4.5, 2.2, 1.1, 0.56, 0.28, 0.014, 0.070, and 0.035 µL/mL. After incubation at 37 • C for 24 h, MICs were measured by visual inspection of the turbidity of broth in tubes (Yang et al., 2020). That is to say, if the test tube was still clear and transparent (un-cloudy) after incubation, the cells cannot grow at that concentration and then the MIC value was obtained. Vancomycin was used as a positive anti-S. aureus control. This assay was carried out in triplicate.

Preparation of Bacterial Samples for RNA-Seq
A conical flask (150 mL) with TSB (45 mL) was inoculated with fresh S. aureus 29213 strain culture (500 µL) which mentioned above and then incubated at 37 • C with 250 rpm shaking. Cell samples were taken for LCEO stimulation as the optical density value (OD600) at 0.8 (1 × 10 6 CFU/mL) and subsequently, 3 samples of O_S29 group (O_S29_1, O_S29_2, O_S29_3) were treated with 1/4 MIC (0.07 µL/mL) of LCEO for 15 min. The control group (S29 group, 3 samples) was without incubated with LCEO and were also kept in the same conditions for 15 min.
RNA Extraction and Purification, cDNA Library Construction for RNA-Seq, RNA-Seq, and Data Analysis Total RNA was extracted from 6 bacterial samples (O_S29_1, O_S29_2, O_S29_3, S29_1, S29_2, S29_3) separately by a Total RNA Isolation System (Promega) according to the manufacturer's protocol. The quality of the RNA samples was examined using the Agilent 2100 Bioanalyzer. Library construction and Illumina sequencing was performed by Novogene China. An RNA-seq analysis was performed according to the protocol recommended by the manufacturer (Illumina Inc.) (Qiao et al., 2016;Liu et al., 2019).
For data analysis, RNA-Seq reads were mapped to the reference genome of S. aureus Newman (NC_009641). Differentially expressed genes (DEGs) were identified by the edgeR package based on Genes with two criterions: |log2(fold change)| > 0.58499 (|fold change| >1.50003) and p < 0.05

Quantitative Real-Time PCR
To ensure the accuracy of RNA-seq data, 8 DEGs were selected for quantitative real-time PCR (RT-PCR) verification and they were based on two standards: (a) Functional classification. The 8 DEGs were selected from different GO terms and KEGG pathways which were significantly enriched. (b) The value of fold change. 4 DEGs were up-regulated and the other 4 DEGs down-regulated and with |fold change| > 1.8. All experiments were performed in triplicate with at least three independent experiments.
Quantitative real-time PCR assay was carried out as described previously with some modifications (Park et al., 2019). The cDNA for qRT-PCR was synthesized using a quanti tect reverse transcription kit [Sangon Biotech (Shanghai) Co., Ltd.]. QRT-PCR was performed using specific primer pairs and iQTM SYBR Green supermix [Sangon Biotech (Shanghai) Co., Ltd.]. The primers were designed with the Primer-BLAST (NCBI) and purchased from Sangon Biotech (Shanghai) Co., Ltd., as shown in Table 2. The amplification and detection of PCR products were performed using CFX Connect TM Real-Time System (Bio-Rad). The thermal cycling conditions were as follows. After activation of the polymerase and a DNA denaturation step at 95 • C, 40  amplification cycles were performed with a denaturation step at 95 • C followed by an annealing and extension step at 55 • C. The cDNA values were normalized with the value of 16s rRNA, which was constant in different conditions (data not shown).

Staphyloxanthin Levels Assay
Staphyloxanthin production was measured as described previously with some modifications (Singh et al., 2018). Briefly, 8 h S. aureus 29213 strain cultures were diluted 1:100 in tubes with 3.5 mL TSB at 250 RPM shaken and 37 • C for 24 h with different lever of LCEO: 1/2 (0.14 µL/mL), 1/4 (0.07 µL/mL), and 1/8 (0.035 µL/mL) MIC. The control samples (without LCEO) were incubated in the same culture conditions. All the culture in tubes were centrifugated and then washed with distilled water for two times. Each tube was added into 2 mL methanol for the ultrasonic extraction of staphyloxanthin under 600 watt at 55 • C for 20 min, and the supernatant with staphyloxanthin was obtained by centrifugation. The value of OD465 was measured as the biosynthesis level of staphyloxanthin. Three independent experiments were run in triplicate.

Statistical Analysis of Quantitative Real-Time PCR and Staphyloxanthin Levels Assay
Statistical analysis was performed using SPSS (version 16.0, SPSS Inc., Chicago, IL, United States). Differences between groups were analyzed using one-way analysis of variance (ANOVA) followed by post hoc analysis using Tukey test. Differences were considered significant at P < 0.05.

General Features of the Transcriptome Profile
RNA-Seq generated total 9,494,858 to 11,617,170 clean reads from LCEO-treated (O_S29) and control (S29) cDNA libraries, respectively. These clean reads were mapped to the reference genome of S. aureus Newman (NC_009641) (Supplementary Data Sheet 1). RNA-seq data were available for further transcriptome analysis as the gene expression levels of three independent sample in each group were highly consistent. The expression level of genes were determined by average Per Kilobase of transcript per Million mapped reads (RPKM, full information of reads quantification) and these values showed that 2989 total genes were expressed in LCEO-treated and control groups.
The raw sequence data reported in this paper have been deposited in the Genome Sequence Archive (Genomics, Proteomics and Bioinformatics 2017) in BIG Data Center (Nucleic Acids Res 2019), Beijing Institute of Genomics Frontiers in Microbiology | www.frontiersin.org (BIG), Chinese Academy of Sciences, under accession numbers CRA002405 that is publicly accessible at https://bigd.big.ac. cn/gsa.

GO Functional Enrichment Analysis
The 542 DEGs were processed by Gene Ontology (GO) enrichment analysis for understanding the stress response at genetic level of S. aureus 29213 which was treated by 1/4 MIC LCEO. As shown in Supplementary Table 2, 910, 180, and 451 specific GO terms in biological process, cellular FIGURE 4 | Significantly enriched KEGG pathway "Histidine metabolism" (A) and "Lysine biosynthesis" (B) (from KEGG database, green rectangles represent proved Organism-specific gene product and green frames represent down-regulation and red frames represent up-regulation). component and molecular function were reported, respectively. In Figure 2, 30 significantly enriched GO terms from three specific categories were showed.
Ten GO terms were enriched in the biological process, including amino acid transport, anion transport, organic anion transport, organic acid transport, carboxylic acid transport, amino acid transmembrane transport, anion transmembrane transport, organic acid transmembrane transport, membrane fusion, and cell adhesion.
Ten GO terms were enriched in the cellular component, including extracellular space, viral envelope, viral membrane, extracellular region, cell wall, extracellular region part, integral component of membrane, virion, intrinsic component of membrane, and immunological synapse.
In the category of molecular function, Ten GO terms (amino acid transmembrane transporter activity, solute:sodium symporter activity, hydro-lyase activity, organic anion transmembrane transporter activity, iron ion binding, symporter activity, organic acid transmembrane transporter activity, solute:cation symporter activity, carboxylor carbamoyltransferase activity, and carboxylic acid transmembrane transporter activity) were significantly enriched.

KEGG Pathway Enrichment Analysis
The 542 DEGs were mapped to the KEGG database, and then analyzed by KEGG pathway enrichment analysis to understand the involved pathways (Supplementary Table 3). As shown in Figure 3, KEGG pathways: histidine metabolism, fatty FIGURE 5 | Significantly enriched KEGG pathway "Two-component system" (A) and "Peptidoglycan biosynthesis" (B) (from KEGG database, green rectangles represent proved Organism-specific gene product and green frames represent down-regulation and red frames represent up-regulation). acid degradation, nitrogen metabolism, pyruvate metabolism, carotenoid biosynthesis were significantly enriched.
DEGs Involved in "Histidine Metabolism," "Lysine Biosynthesis," "Two-Component System" and "Peptidoglycan Biosynthesis" Pathways Study has shown that the leakage of intracellular macromolecules was caused by the destructive effect of 0.5 mg/mL (1 MIC) LCEO on the cell membrane of MRSA . Actually, based on the transmission electron microscopy (TEM) micrographs provided by this reference above, the cell wall of the bacterium was also damaged by LCEO. Therefore, LCEO has the effect of disintegrating cell walls and membranes of S. aureus. In this experiment, 1/4 mic of LCEO also triggered cell wall (membranes) stress to S. aureus 29213.
Histidine protein kinases (HPKs) are a large family of signaltransduction enzymes. The typical HPK is a transmembrane receptor and it can autophosphorylate the conserved histidine residue of itself (Wolanin et al., 2002). Bacterial two-component signaling system allows organisms to perceive and respond to change in many different environmental conditions. Some signal-transduction were achieved by the sensor histidine kinase (HK) and phosphorylation response regulator (RR) of twocomponent system (Harms et al., 2001). Well, it is easy to know that histidine is urgently needed in the repair process once the cell wall (membranes) is damaged. There are ten steps of biosynthesis from PRPP to L-histidine, and the 3 key enzymes (hisB, hisC, and hisI) are the utmost importance in this histidine biosynthesis pathway of S. aureus (Henriksen et al., 2011) in LCEO treated S. aureus 29213 and were significantly up-regulated as comparison with the control. L-histidine can be metabolized into L-glutamate through a series of reactions. Simultaneously, 4 (all the) annotated organism-specific genes in these steps FIGURE 6 | Significantly enriched KEGG pathway "Fatty acid degradation" (from KEGG database, green rectangles represent proved Organism-specific gene product and green frames represent down-regulation and red frames represent up-regulation).
FIGURE 7 | Significantly enriched KEGG pathway "Nitrogen metabolism" (from KEGG database, green rectangles represent proved Organism-specific gene product and green frames represent down-regulation and red frames represent up-regulation).
FIGURE 8 | Significantly enriched KEGG pathway "Carotenoid biosynthesis" (from KEGG database, green rectangles represent proved Organism-specific gene product and green frames represent down-regulation and red frames represent up-regulation).
were also down-regulated. These results suggested that the histidine biosynthesis of S. aureus 29213 was increased under the stimulation of 1/4 MIC of LCEO ( Figure 4A).
FIGURE 9 | Staphyloxanthin level assay. OD465 value represent the level of staphyloxanthin production in different groups. The data represent the means ± standard deviation of the results. Lowercase marked without the same superscripts (a,b) differ significantly (P < 0.05).
Like alanine and glutamate, L-lysine is not only used to synthesize proteins but also a component of peptidoglycan of cell wall in gram-positive bacteria. Meso-2,6-diaminopimelate can be utilized in the biosynthesis process of peptidoglycan and L-lysine (Grundy et al., 2003). Results in Figure 4B showed that more than half of the genes encoding the enzymes which are regulating lysine biosynthesis were up-regulated, and they were mainly in the early stages involved in this pathway. LysA [EC:4.1.1.20] is a decarboxylase enzyme to catalyze meso-2,6diaminopimelate into L-lysine and the fold change of its encoded gene is + 1.437 (slightly below the baseline of this analysis, and |fold change| <1.50003 will not be highlighted). These results suggested that the lysine biosynthesis of S. aureus 29213 has increased under the stress of LCEO.
The gene (vraS) encoding the sensor histidine kinase VraS in two-component system was up-regulated in LCEO treated S. aureus 29213 ( Figure 5A). Once VraS is stimulated by cell wall stress stimulon, it phosphorylates the histidine residues of itself into VraR, then VraR sends signals to downstream for peptidoglycan biosynthesis (Kuroda et al., 2003;Muthaiyan et al., 2008). So, LCEO is a VraSR-dependent cell wall stress stimulon.
Gene encoding pbp4 [EC:3.4.16.4], one of the enzymes in the final steps of peptidoglycan biosynthesis, was up-regulated by LCEO stress (Figure 5B), which suggests that the biosynthesis of peptidoglycan may be increased. Recent study also showed that there is a close relationship between pbp4 and peptidoglycan fragments (Maya-Martine et al., 2019).
From these results, we can conclude that LCEO is a cell wallactive agent , and the stress response of LCEO treated S. aureus 29213 could be explained by the transcriptional profiling evidences of "Histidine metabolism, " "Lysine biosynthesis, " "Two-component system, " and "Peptidoglycan biosynthesis" pathways.
DEGs Involved in "Fatty Acid Degradation" and "Nitrogen Metabolism" Pathways Fatty acids are required to maintain cell viability and S. aureus ingest exogenous fatty acids into its phospholipid bilayer (Delekta et al., 2018). The fatty acids in the cytoplasm are synthesized by fadE [EC:6.2.1.3] into hexadecanoyl-CoA then transported to the mitochondria by carnitine palmitoyl transferase I (CPT1) and CPT2. Hexadecanoyl-CoA will be degraded into acetyl-CoA after seven cycles. Each cycles is consists of four steps (dehydrogenation, hydration, dehydrogenation and transformation) and then total 4 mol ATP will be generated by the 2 dehydrogenation processes in a cycle (Parekh, 1977). In present results (Figure 6) were up-regulated, which indicated the fatty acid degradation activity of LCEO treated S. aureus 29213 was increased. Research showed that the hexose monophophate pathway (HMP) and its key enzyme (glucose-6-phosphate dehydrogenase) activity of methicillin-resistant S. aureus were inhibited by LCEO treatment . However, the |fold change| of genes involved in HMP is under 1.5 in present results. That may be because our S. aureus 29213 was treated with LCEO for a short time (only 15 min). Although there was no significantly enriched, present data showed that almost all the screened genes FIGURE 10 | Significantly enriched KEGG pathway "Pyruvate metabolism" (from KEGG database, green rectangles represent proved Organism-specific gene product and green frames represent down-regulation and red frames represent up-regulation, yellow frames represent both down and up-regulation).
The existence of dissimilatory nitrate reduction makes the microbes with multiple environmental adaptabilities (such as the metabolism in both oxic and anoxic status) and resource independence (like anaerobic metabolism under indirect nitrate supplyment) (Kamp et al., 2015). Denitrification is a type of respiration. When bacteria are exposed to low or poor oxygen environment, nitrate or nitrite is reduced as a terminal electron acceptor in this respiration (Cabello et al., 2004). Both of these methods can also provide energy to microorganisms (higher than that of fermentation but lower than that of aerobic respiration). The reductases NarGHIJ can reduce nitrate into nitrite, and the reductases NirBD can reduce nitrite into ammonia. As the genes encoding NarGHIJ and NirBD were donw-regulated (Figure 7), we assumed that LCEO inhibits the nitrate respiration function of S. aureus 29213. Related studies have shown that the nitrogen metabolism is closely related to the antibiotic resistance (Borriello et al., 2004;Yan et al., 2013). Therefore, we assumed that the combination of LCEO and antibiotics can reduce the use of antibiotics on S. aureus.
From this part, we assumed that LCEO inhibited the metabolism of carbohydrate and nitrogen, so that in order to survive, the bacteria has to increase activity of fatty acid degradation.

DEGs Involved in "Carotenoid
Biosynthesis" and "Two-Component System" Pathways Staphyloxanthin is a golden pigment of S. aureus and it distinguishes S. aureus from other Staphylococci even grampositive cocci. Staphyloxanthin is biosynthesized by a series of steps and is the only C 30 golden carotenoid of membraneembedded. Staphyloxanthin is an important antioxidant as well as an anti-host neutrophil-based killing substance. Also, the biosynthesis of staphyloxanthin is under the oxygen-sensing and redox-signaling regulation of the two-component system (Gao et al., 2017;Hall et al., 2017). CrtM and CrtN are the key enzymes of the golden staphyloxanthin biosynthetic process and they could be important targets of novel drug design to handle the infections of S. aureus (Yang et al., 2019). According to the results of transcriptome analysis, all the DEGs involved in carotenoid biosynthesis pathway were downregulated (Figure 8). Furthermore, the down-regulation of CydA which is involved in the redox related part of two-component system pathway as well as relevant to aerobic respiration was also recorded (Figure 5A). To further confirm the transcriptomic reports of present analysis, we treated S. aureus 29213 strains with LCEO and found that LCEO inhibited the staphyloxanthin lever of S. aureus 29213 in a dose-dependent manner (Figure 9). As staphyloxanthin is related to the antioxidant activity of LCEO, this result suggests that S. aureus may be more vulnerable to antibacterial agent or environment after treated with LCEO (Wang et al., 2012;Bi et al., 2017;She et al., 2019). If it is true, LCOE will have a broader application prospect.
DEGs Involved in "Staphylococcus aureus Infection" Pathway S. aureus has ability to colonize to the host and to form biofilm, which is likely to be partially determined by its ability of adhesion (Foulston et al., 2014;Scherr et al., 2014). Vitro research showed that the proteins ClfB, IsdA, SdrC/D, and SasG are able to promote adhesion of staphylococci (Foster, 2009). Related studies have shown that the formation of biofilm of S. aureus is closely related to the genes which encoding these proteins (O'Brien et al., 2002;Feuillie et al., 2017;Liu et al., 2018). Recent report has shown that 1 mg/mL was the minimal biofilm inhibitory concentration of LCEO on S. aureus and the MIC is 0.5 mg/mL (Mei, 2019). Thus, the up-regulated genes encoding surface proteins ClfB and SdrC/D (data not shown) which promote the adhesion of S. aureus suggested the enhanced biofilm viability to adapt to the unsuitable conditions.

DEGs Involved in "Pyruvate Metabolism" Pathways
Pyruvate is the final product of glycolysis pathway, which is reduced to lactic acid for energy supply in the cytoplasm, or oxidized to acetyl CoA in the mitochondria, and then oxidized to carbon dioxide and water to complete the aerobic energy supply process of glucose. Pyruvate can also convert into sugars, fats and amino acids through the acetyl CoA and tricarboxylic acid cycles. As a consequence, pyruvate plays an important pivotal role in the metabolism of the three major nutrients (Hills, 1938;Sevag et al., 1950). Pyruvate dehydrogenase complex (PDHC) is a group of speed-limiting enzymes that catalyzes the irreversible oxidative decarboxylation of pyruvate into acetyl-CoA. In LCEO treated S. aureus 29213, the genes encoding PDHC which consists of pyruvate dehydrogenase E1 component alpha subunit  were down-regulated compared with the control (Figure 10). The results may suggest the reduced lever of acetyl-CoA which comes from pyruvate of LCEO treated S. aureus 29213. In addition, this may be related to the inhibition effect of carbohydrate metabolism mentioned above.

RT-PCR Validation
Results of qRT-PCR assay showed that the mRNA expression levels of 8 selected DEGs were highly consistent with the results of RNA-Seq analysis (Figure 11). This result indicated that RNA-Seq assay was properly performed.

CONCLUSION
Evidences of genetic lever for the stress response of S. aureus 29213 strain under sub-lethal concentration of LCEO by transcriptomic analysis were founded in this study. 300 and 242 genes were significantly up and down-regulated in LCEO-treated S. aureus 29213. Up-regulated genes were mainly involved in cell membrane and down-regulated genes were mainly involved in oxidation-reduction process. These findings indicated that the cell membrane stress of S. aureus 29213 was occurred by LCEO. At the same time, we reported that LCEO can significantly affect the staphyloxanthin synthesis of S. aureus 29213 for the first time, which was closely related to the redox state of S. aureus 29213. LCEO may also inhibit the nitrogen metabolism and promote fatty acid degradation of S. aureus 29213. These evidences expanded the knowledge of stress response of S. aureus 29213 strain under sub-lethal concentration of LCEO.

DATA AVAILABILITY STATEMENT
The raw data of this manuscript have been deposited in the Genome Sequence Archive (Genomics, Proteomics & Bioinformatics 2017) in BIG Data Center (Nucleic Acids Res 2019), Beijing Institute of Genomics (BIG), Chinese Academy of Sciences, under accession number CRA002405 that is publicly accessible at https://bigd.big.ac.cn/gsa.

AUTHOR CONTRIBUTIONS
WZ and HS designed this work. All authors reviewed and approved the final manuscript.