RNA‑seq analyses of antibiotic resistance mechanisms in Serratia marcescens

The present study aimed to further clarify the genetic mechanisms responsible for the antimicrobial resistance of Serratia marcescens (S. marcescens) using RNA sequencing. Three drug‑susceptible S. marcescens strains (named MYQT1, MYQT2, and MYQT3) and three multidrug‑resistant S. marcescens strains (named MYQT4, MYQT5, and MYQT6) were isolated from six different patients and subjected to RNA sequencing. Differentially expressed genes (DEGs) between the multidrug‑resistant S. marcescens strains and drug‑susceptible strains were screened and compared, followed by functional enrichment analysis. In addition, a protein‑protein interaction (PPI) network was constructed, and significant modules were extracted from it. Genes enriched in the significant modules were subjected to further enrichment analysis. MYQT3 had very a different expression pattern from MYQT1 and MYQT2, and thus, MYQT3 was excluded from the following analysis. A total of 225 DEGs were identified, of which SMDB11_RS09300 (GTP cyclohydrolase FolE2) was the most significantly upregulated with a log2 FC of 6.4; these DEGs were enriched in different GO terms, including hydrogen sulfide biosynthetic process, sulfur compound transmembrane transporter activity, and ABC transporter complex. Additionally, several genes were identified to be important genes in the PPI network, including SMDB11_RS17755 (upregulated; glutamate synthase large subunit), SMDB11_RS00590 (upregulated; sulfite reductase subunit α), and SMDB11_RS04505 (upregulated; cystathionine β‑synthase). Thus, SMDB11_RS09300, SMDB11_RS17755, SMDB11_RS00590, and SMDB11_RS04505 may play significant roles in the antimicrobial resistance of S. marcescens by participating in folate metabolism or the integrity of cell membranes. However, further experiments are required to clarify these findings.