<abstract xmlns="http://www.w3.org/1999/xhtml">Transcriptome changes of biofilm-forming <italic>Staphylococcus epidermidis</italic> response to total alkaloids of <italic>Sophorea alopecuroides</italic> was observed. Bioinformatic analyses were further used to compare the differential gene expression between control and the treated samples. It was found that 282 genes were differentially expressed, with 92 up-regulated and 190 down-regulated. These involved down-regulation of the sulfur metabolism pathway. It was suggested that inhibitory effects on <italic>Staphylococcus epidermidis</italic> and its biofilm formation of the total alkaloids of <italic>S. alopecuroides</italic> was mainly due to the regulation of the sulfur metabolism pathways of <italic>S. epidermidis</italic>.</abstract>

Transcriptome changes of biofilm-forming Staphylococcus epidermidis response to total alkaloids of Sophorea alopecuroides was observed. Bioinformatic analyses were further used to compare the differential gene expression between control and the treated samples. It was found that 282 genes were differentially expressed, with 92 up-regulated and 190 down-regulated. These involved down-regulation of the sulfur metabolism pathway. It was suggested that inhibitory effects on Staphylococcus epidermidis and its biofilm formation of the total alkaloids of S. alopecuroides was mainly due to the regulation of the sulfur metabolism pathways of S. epidermidis.

<div xmlns="http://www.w3.org/1999/xhtml"><italic>S. epidermidis</italic> is an opportunistic pathogen; however, it is capable of storing and transmitting drug-resistant genes (Bloemendaal <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_001">2010</a>; Otto, <a ref-type="bibr" href="#j_pjm-2018-024_ref_012">2013</a>). The ability of <italic>S. epidermidis</italic> to form biofilms <italic>in vivo</italic> allows it to be highly resistant to chemotherapeutics and can ultimately lead to chronic disease (Cerca <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_002">2005</a>). Some Chinese medicinal herbal products have been found to prevent biofilm formation in bacteria (Guan <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_004">2013</a>; Wang <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_015">2015</a>; Yang <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_016">2016</a>). Total alkaloids of <italic>S. alopecuroides</italic> (TASA) are an alkaloid mixture that has demonstrated a better inhibitory effect on the late stage of <italic>S. epidermidis</italic> biofilm thickening than ciprofloxacin (CIP), and erythromycin (ERY) (Li <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_008">2016</a>).In order to investigate the inhibitory mechanism of TASA, the global transcriptome changes of biofilm-forming <italic>S. epidermidis</italic> (ATCC 35984) in response to TASA were analysed in this study. Some pre-experiments were performed to determine the transcriptome sequencing conditions. In the pre-experiments, the relative expression quantity of four main biofilm formation related genes (<italic>ica, sigB, agr</italic> and <italic>fbe</italic>) were compared after the <italic>S. epidermidis</italic> treatment by TASA in different sub-MIC (1/2 MIC and 1/4 MIC) and time periods (6 h, 12 h, and 24 h). The results showed that these four genes were down-regulated significantly when the <italic>S. epidermidis</italic> was treated by 1/2 MIC (12.5 mg/ml) of TASA for 12 h. Based on the pilot experiment results, the transcriptome sequencing was used to investigate the transcriptional differences of <italic>S. epidermidis</italic> ATCC 35984 in a blank-medium control group and the treatment group, cultured with 1/2 MIC of TASA for 12 h. The inhibitory mechanism of Chinese medicine TASA was examined through the bioinformatic analyses.<italic>S. epidermidis</italic> ATCC 35984 was cultured overnight at 37°C, 120 rpm. The resulting suspension was diluted to 0.5 McFarland with TSB medium and 3 ml of the diluted suspension was seeded into a six-well plate. The plate was sealed with an aseptic glass cover and placed in the incubator upside down. After 12 h of incubation at 37°C, the liquid in each well was removed and washed twice with PBS. The biofilm bacteria adhered to the glass cover in the blank control samples were collected with 1 ml TSB, and centrifuged at 12000 r/min for 10 min at 4°C. For TASA-treated samples, the biofilm bacteria on the glass covers were incubated for an additional 12 h in the presence of TASA at 12.5 mg/ml (1/2 MIC) prior to the bacteria harvesting, which followed the same methodology as the blank control. Total mRNA extraction was performed with RNAprep pure Cell/Bacteria Kit (TianGen), according to the manufacturer’s instructions, then stored at –80°C, and prepared for transcriptome sequencing.Standard sequencing procedure and data analysis were performed by using Illumina HiSeq 2500 high-throughput sequencing platform from Biomarker Technologies Co., Ltd (China). Separate sequence read dataset was used as inputs into DESeq2 package (Love <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_009">2014</a>) to analysis the unigenes expression based on RPKM (reads per kilobase transcriptome per million mapped reads). A fold change of ≥ 2 and a minimum false discovery rate (FDR) of &lt; 0.01 were accepted as indicators of the differentially expressed genes (DEGs) after Benjamini-Hochberg <italic>post hoc</italic> correction. The DEGs were BLASTX against the Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Group (COG) and Gene Ontology (GO) to identify their predicted biological function. In order to verify the RNA-Seq results, the expression levels of selected DEGs were quantified by RT-qPCR. First, Nanodrop 2000 (Thermo Scientific, Wilmington, USA) was used to measure the concentration of total mRNA and 500 ng of RNA was reverse-transcribed into cDNA with PrimeScriptTM RT Master Mix (TaKaRa). RT-qPCR was performed with an iQ5 light cycler (Bio-Rad) by using SYBR® Premix EX TaqTM II (TaKaRa) in a 20 μl reaction volume, which consisted of 0.25 mM of each primer, 10 μl of SYBR Premix Ex Taq II and 1 μL of template cDNA. PCR conditions were as follows: 30 s at 95°C followed by 40 cycles of 5 s at 95°C and 30 s at 58°C. The reference gene was 16S rRNA, data were acquired through Bio-Rad and analyzed by using the 2–ΔΔ<italic>CT</italic> method (Livak and Schmittgen, <a ref-type="bibr" href="#j_pjm-2018-024_ref_006">2001</a>). All experiments were performed in triplicate.282 unigenes were differentially expressed under the screening criteria (fold change ≥ 2 and FDR &lt; 0.01). Within these DEGs, 92 unigenes were significantly up-regulated and 190 unigenes were significantly down-regulated. The DEGs were searched against the KEGG database to identify their biological pathways. The 20 greatest enriched pathways were listed in Fig. <a ref-type="fig" href="#j_pjm-2018-024_fig_001">1</a>. Pathways (represented by symbols) in the upper left quadrant of Fig. <a ref-type="fig" href="#j_pjm-2018-024_fig_001">1</a> contain DEGs with more significant and reliable enrichment levels, therefore the sulfur metabolism pathway ranked first. In sulfur metabolism pathway, there are nine DEGs which were all down-regulated significantly, corresponding to the encoded enzymes included serine acetyltransferase (CysE, EC 2.3.1.30), phosphoadenosine phosphosulfate reductase (CysH, EC 1.8.1.2), sulfite reductase (CysI and CysJ, EC 1.8.4.8), cystathionine beta-lyase (Cbl, EC 2.5.1.48), cysteine synthase (CysK and CysM, EC 2.5.1.47), sulfate adenylyl transferase (Sat, EC 2.7.7.4) and adenylylsulfate kinase (CysC, EC 2.7.1.25) (Fig. <a ref-type="fig" href="#j_pjm-2018-024_fig_002">2</a>). These enzymes directly affect the content of sulfur metabolites such as cysteine (Cys), methionine (Met), glutathione (GSH), etc. RT-qPCR results further demonstrated that all nine genes showed similar expression patterns to those of RNA-seq, confirming the RNA-seq based transcriptome datasets were accurate and robust (Fig. <a ref-type="fig" href="#j_pjm-2018-024_fig_003">3</a>).<figure id="j_pjm-2018-024_fig_001" fig-type="figure"><h2>Fig. 1.</h2><figCaption>KEGG pathway scatter diagram of differentially expressed genes. Each pattern represents a KEGG pathway with its name listed in the right. X-coordinate represents the enrichment factor showing the ratio of all genes annotated in the pathway, in contrast to the differentially expressed genes annotated in the same pathway. The smaller the enrichment factor, the enrichment level of the differentially expressed genes was more significant in this pathway. Y-coordinate represents log10(Q_value). The higher the absolute Q_value, the enrichment level of the differentially expressed genes was more significant and reliable in this pathway.</figCaption><img xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="graphic/j_pjm-2018-024_fig_001.jpg" src="https://sciendo-parsed.s3.eu-central-1.amazonaws.com/64727a9b215d2f6c89dc9d0a/j_pjm-2018-024_fig_001.jpg?X-Amz-Algorithm=AWS4-HMAC-SHA256&amp;X-Amz-Date=20240505T094914Z&amp;X-Amz-SignedHeaders=host&amp;X-Amz-Expires=18000&amp;X-Amz-Credential=AKIA6AP2G7AKP25APDM2%2F20240505%2Feu-central-1%2Fs3%2Faws4_request&amp;X-Amz-Signature=dfd501fc322dfea70b8a673eafc4c444b753aafdb5feec14e0a1d1a9b3614ad7" class="mw-100"/></figure><figure id="j_pjm-2018-024_fig_002" fig-type="figure"><h2>Fig. 2.</h2><figCaption>Distribution of the differentially expressed genes in a sulfur metabolism pathway. Compared to the control group, the enzyme marked by a green frame is related to a drown-regulated gene product. The number in the frame represents the enzyme code.</figCaption><img xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="graphic/j_pjm-2018-024_fig_002.jpg" src="https://sciendo-parsed.s3.eu-central-1.amazonaws.com/64727a9b215d2f6c89dc9d0a/j_pjm-2018-024_fig_002.jpg?X-Amz-Algorithm=AWS4-HMAC-SHA256&amp;X-Amz-Date=20240505T094914Z&amp;X-Amz-SignedHeaders=host&amp;X-Amz-Expires=18000&amp;X-Amz-Credential=AKIA6AP2G7AKP25APDM2%2F20240505%2Feu-central-1%2Fs3%2Faws4_request&amp;X-Amz-Signature=338b897f85809b2be2397b9750964e2eb04e582e54adb6717944aeffc1cc3ce0" class="mw-100"/></figure><figure id="j_pjm-2018-024_fig_003" fig-type="figure"><h2>Fig. 3.</h2><figCaption>Comparison of folds changes detected by RNA-seq (RPKM) and q-PCR.Bars represent the log2 value of the fold changes of the gene expressions between the control and TASA treatment groups. On the top of the bars are the abbreviations of enzymes encoded by the nine genes.</figCaption><img xmlns:xlink="http://www.w3.org/1999/xlink" xlink:href="graphic/j_pjm-2018-024_fig_003.jpg" src="https://sciendo-parsed.s3.eu-central-1.amazonaws.com/64727a9b215d2f6c89dc9d0a/j_pjm-2018-024_fig_003.jpg?X-Amz-Algorithm=AWS4-HMAC-SHA256&amp;X-Amz-Date=20240505T094914Z&amp;X-Amz-SignedHeaders=host&amp;X-Amz-Expires=18000&amp;X-Amz-Credential=AKIA6AP2G7AKP25APDM2%2F20240505%2Feu-central-1%2Fs3%2Faws4_request&amp;X-Amz-Signature=20a658761862b87016f917474a9442ca463810098cc90c4f3e9261c971bf5db6" class="mw-100"/></figure>In this research, according to the analysis of transcriptome data, the sulfur metabolism pathway in <italic>S. epidermidis</italic> was greatly influenced by TASA. Sulfur metabolism is an important metabolic pathway in bacteria and its metabolites are involved in many physiological and biochemical processes in cells (Zeng <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_017">2013</a>). Cysteine (Cys), methionine (Met) and glutathione (GSH) are associated with bacterial activity and biofilm formation (Murillo <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_010">2005</a>; Gales <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_003">2008</a>; Soutourina <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_014">2009</a>; Li <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_006">2013</a>; Solis <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_013">2016</a>). A link between sulfate assimilation, Met biosynthesis and biofilm formation have been found in <italic>Candida albicans</italic> (Murillo <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_010">2005</a>). Deletion of the <italic>ecm17</italic> gene encoding the sulfite reductase beta subunit resulted in reduced adhesion and poor biofilm formation in <italic>C. albicans</italic> (Li <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_006">2013</a>). Additionally, deletion of the master regulator of Cys metabolism <italic>cymR</italic> in <italic>S. aureus</italic> also results in diminished biofilm formation (Soutourina <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_014">2009</a>). A recent research has provided further evidence of a role for sulfate assimilation and Cys/Met biosynthesis in <italic>S. epidermidis</italic> ATCC 35984 (RP62A) biofilm formation (Solis <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_013">2016</a>). According to the comparison of non-biofilm forming <italic>S. epidermidis</italic> ATCC 12228 and biofilm-forming ATCC 35984, Solis <italic>et al</italic>. (<a ref-type="bibr" href="#j_pjm-2018-024_ref_013">2016</a>) showed that the sulfate assimilation and cysteine/methionine biosynthesis pathways in ATCC 35984 contained elevated levels (~25% increase) of methionine that were likely linked to biofilm formation. GSH, another important sulfur metabolite, which is associated with intracellular reactive oxygen species (ROS) and hydrogen sulfide (H2S) also played an important role in biofilm formation (Gales <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_003">2008</a>; Klare <italic>et al.</italic>, <a ref-type="bibr" href="#j_pjm-2018-024_ref_005">2016</a>; Ooi and Tan, <a ref-type="bibr" href="#j_pjm-2018-024_ref_011">2016</a>). In this study, TASA significantly influenced the sulfur metabolism by down-regulation of nine important genes in this pathway. The concentration of Cys, Met and GSH were directly affected, and most probably finally disrupted the biofilm formation of <italic>S. epidermidis</italic>. This result also validated the above-mentioned relationship between sulfur metabolism and biofilms.In the pre-experiments, the relative expression of biofilm-related genes including <italic>ica, sigB, agr</italic> and <italic>fbe</italic> in <italic>S. epidermidis</italic> was measured. Results showed that these four genes were differentially expressed after treatment with TASA at 1/2 MIC for 12 h. However, these four genes did not appear in the 282 DEGs obtained by transcriptome sequencing due to the DEGs screening criteria (fold change ≥ 2 and FDR &lt; 0.01). It was suggested that genes related to biofilm formation did have differential expression, but they may not be directly regulated by TASA. Instead, their expression may be indirectly affected by the changes of other metabolic pathways (<italic>i.e</italic>. sulfur metabolism).For a long time, the formation of biofilm has been widely recognized as a dynamic process and regulated by some genes including <italic>ica, sigB, agr</italic> and <italic>fbe etc</italic>. Now, it was shown that there is a relationship between sulfur metabolism and biofilm formation in <italic>S. epidermidis</italic>. Based on the results of transcriptome analysis, it was suggested that TASA’s inhibitory effects on <italic>S. epidermidis</italic> and its biofilm formation is mainly due to the regulation of the sulfur metabolism pathway.The regulatory mechanism of TASA has not been analyzed due to its multi-channeled and multi-targeted actions. In the study, the effect of TASA on <italic>S. epidermidis</italic> was comprehensively analyzed by using RNA-seq. The enrichment of DEGs showed that metabolism, genetic information processing and environmental information processing were greatly influenced by TASA. The sulfur metabolism pathway was the most significant with all of the key enzymes in this pathway being down-regulated. Given the important physiological role of sulfur metabolism in bacteria and its effect on biofilm formation, it is concluded that the inhibitory effect of TASA on <italic>S. epidermidis</italic> and its biofilm formation are mainly due to its actions on this pathway.In conclusion, according to the analysis of differential expression and metabolic pathway enrichment, the molecular mechanism of TASA regulation of pathogenicity, virulence and metabolism of <italic>S. epidermidis</italic> was further characterized. The sulfur metabolism pathway was identified and results could provide valuable information for follow-up studies to examine the regulatory mechanism of TASA on <italic>S. epidermidis</italic>, and could serve as the basis for exploring potential drug targets.</div>

Global Transcriptome Changes of Biofilm-Forming <italic>Staphylococcus epidermidis</italic> Responding to Total Alkaloids of <italic>Sophorea alopecuroides</italic>

Polish Journal of Microbiology

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Global Transcriptome Changes of Biofilm-Forming  Responding to Total Alkaloids of 

right-running

Transcriptome changes of biofilm-forming Staphylococcus epidermidis response to total alkaloids of Sophorea alopecuroides was observed. Bioinformatic analyses...

Global Transcriptome Changes of Biofilm-Forming Staphylococcus epidermidis Responding to Total Alkaloids of Sophorea alopecuroides

Article Category: short-communication

Published Online: Jun 30, 2018

Page range: 223 - 226

Received: Jul 06, 2017

Accepted: Dec 03, 2017

DOI: https://doi.org/10.21307/pjm-2018-024

Keywords
sulfur metabolism, total alkaloids of, transcriptome

© 2018 Cui-Ping Guan et al., published by Sciendo

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Global Transcriptome Changes of Biofilm-Forming Staphylococcus epidermidis Responding to Total Alkaloids of Sophorea alopecuroides

Article Category: short-communication

Published Online: Jun 30, 2018

Page range: 223 - 226

Received: Jul 06, 2017

Accepted: Dec 03, 2017

DOI: https://doi.org/10.21307/pjm-2018-024

Keywordssulfur metabolism, total alkaloids of, transcriptome

© 2018 Cui-Ping Guan et al., published by Sciendo

This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License.

Keywords
sulfur metabolism, total alkaloids of, transcriptome