Abstract Shewanella putrefaciens is a Gram-negative bacterium that can cause seafood spoilage under low-temperature conditions. The bacterium easily forms biofilms to enhance its survival in challenging environments. Our previous research revealed that the biofilm formed by S. putrefaciens WS13 under the low temperature (4 °C) has larger biomass and tighter structure than at an optimum growth temperature (30 °C). In this study, comparative transcriptome analysis was further performed to get insights into the global-level of gene expression in the biofilm formed by S. putrefaciens WS13 under the refrigerating and optimal temperatures using Illumina RNA-Sequencing technique. The results revealed that a total of 761 genes were differentially expressed, of which 497 were significantly up-regulated and 264 were significantly down-regulated (p<0.05). The qRT-PCR results of randomly selected differentially expressed genes (DEGs) confirmed the RNA sequencing results. Comparison of transcriptome data revealed 28 significantly changed metabolic pathways under the cold stress, including the down-regulated chemotaxis, and motility, and up-regulated tryptophan metabolism, histidine biosynthesis, and quorum sensing, which benefited the biofilm formation of S. putrefaciens WS13 under the adverse circumstance. This study provided useful data for better understanding of the biofilm formation of S. putrefaciens, and also laid a theoretical foundation for novel vaccine and drug targets against the severe spoilage bacterium under the cold stress. Keywords: Shewanella putrefaciens WS 13, biofilm, cold stress, transcriptome, specific spoilage organism Introduction Shewanella putrefaciens is a severe spoilage bacterium in seafood, particularly under low temperature conditions. The bacterium inhabits ubiquitously in the environment because of its excellent environmental adaptability ([27]Xie et al., 2018; [28]Zhen-Quan et al., 2018). The amount of S. putrefaciens cells has often been used as an index to evaluate the quality of seafood ([29]Gram et al., 2002). S. putrefaciens can reduce trimethylamine oxide, the umami taste substance in seafood, to trimethylamine, and generate histamine and other harmful volatile substances ([30]Vogel et al., 2005), which poses a serious threat to the seafood processing industry ([31]Xie et al., 2011; [32]Hou et al., 2013). S. putrefaciens is easy to adhere to the surface of food processing equipment such as stainless steel to form a biofilm that is composed of polysaccharides, proteins, nucleic acids, and lipids ([33]Bagge et al., 2001; [34]Flemming and Wingender, 2010). Biofilm can enhance the stress tolerance of bacterial cells ([35]Bagge et al., 2001; [36]Flemming et al., 2016; [37]Yan and Xie, 2021). Recently, our prior study revealed that the biofilm formed by S. putrefaciens under 4 °C had 1.61-fold larger biomass and tighter structure than that at 30°C ([38]Yan and Xie, 2020). Therefore, it was speculated that the formation of biofilm could enhance the survival ability of S. putrefaciens under the cold stress. However, few information on transcriptome profiles and regulatory factors of S. putrefaciens biofilm cells under the cold stress is available so far. With increased breakthrough of sequeneing technology, RNA-Sequencing (RNA-Seq) technique has been use to study bacterial differential express genes (DEGs), transcript structures, new transcripts and isomers, and alternative splicing and allele-specific expression under adverse circumstances, such as Clostridium acetobutylicum in a high-salt environment ([39]Ao et al., 2020), and Escherichia coli under different heating methods. In recent years, gene expression between biofilm cells and planktonic cells under stress conditions ([40]Charlebois et al., 2016; [41]Ao et al., 2020), has been investigated at the transcriptomic level in Clostridium acetobutylicum ([42]Dong et al., 2015), Porphyromonas gingivalis ([43]Lo et al., 2009), and Gardnerella vaginalis ([44]Castro et al., 2017). However, current literature on the molecular mechanism of biofilm formation of spoilage bacteria under the cold stress is still rare. Therefore, in this study, we aimed to determine DEGs during the formation of S. putrefaciens WS13 biofilm under the cold stress by transcriptomics analysis. The results in this study will provide crucial clues for the targeted inhibition of the biofilm of S. putrefaciens under the low temperature. Materials and Methods Bacterial Strain and Growth Conditions S. putrefaciens WS13 strain was isolated from spoilage shrimp Litopenaeus vanname in refrigerator ([45]Chen et al., 2019). The isolate was maintained in Luria Broth (LB, Land Bridge Technology, Beijing, China) with 50% (v/v) glycerol at -80 °C freezer in our laboratory at Shanghai Ocean University, Shanghai, China. S. putrefaciens WS 13 was inoculated in 9 mL LB medium (pH 7.4) and incubated at 30°C with shaking at 200 rpm for 12 h, and repeated twice for further analysis. Biofilm Assay Biofilm assay was performed according to the method described by [46]Yan and Xie (2020). S. putrefaciens WS13 were incubated overnight to approximately 8 log colony forming units (CFU) mL^-1 (OD[600nm]≈0.8), and diluted with fresh LB medium (1:100, v/v). A 1 mL of diluted culture was added to each well of sterile 24-well polystyrene microtiter plates. Each sample was tested in six replicates. The samples were incubated at 4°C, and 30°C statically to form biofilms for 24h, and 84 h, respectively. Plastic wraps were used to minimize evaporative loss ([47]Yan and Xie, 2020). RNA Extraction, Library Preparation, and RNA Sequencing Mature biofilm cells of S. putrefaciens WS13 at 4°C and 30°C were harvested at 24h, and 84h, respectively. RNA extraction, cDNA library preparation, and RNA sequencing were carried out as described previously ([48]Ao et al., 2020). Quality Control and Mapping Raw paired end reads were trimmed using Fastp v 0.20.0 software ([49]https://github.com/OpenGene/fastp), and low-quality reads and removing reads with size inferior to 50 bp. The clean reads were separately aligned to the reference genome of S. putrefaciens WS 13 (GenBank: [50]CP028435.1) with orientation mode using HISAT2 v 2.1.0 software ([51]http://daehwankimlab.github.io/hisat2/). Next, the mapped reads of each sample were assembled using StringTie v 1.3.6 software in a reference-based approach. All sequences were quantified as Fragments Per Kilobase Million Mapped Reads (FPKM) by StringTie. The formula was defined as FPKM = 10^6×F/(NL×10^-3), where F is the number of fragments assigned to a certain gene in a certain sample, N is the total number of mapped reads in the certain sample, and L is the length of the certain gene. Differential Expression Gene Analysis DEGs were determined by DEseq2, and genes with FDR <0.05 and |log2 fold change| >1 were identified as DEGs. DEGs were aligned against Gene Ontology (GO, [52]http://geneontology.org/) and Kyoto Encyclopedia of Genes and Genomes (KEGG) ([53]https://www.genome.jp/kegg/) databases. The R package cluster Profiler ([54]http://bioconductor.org/packages/release/bioc/html/clusterProfiler .html) was used to identify enriched GO terms and KEGG pathways with a cut-off of P-value < 0.05. DGEs in biofilm formation of S. putrefaciens WS 13 were further analyzed using hierarchical clustering. FDR<0.05 is the standard for screening genes with significant differential expression. The gene expression level of S. putrefaciens mature biofilm cells at 4°C was used as a reference, whether a gene was up-regulated or down-regulated was determined by comparing its expression level with that at 4°C. Quantitative Real-Time Polymerase Chain Reaction Assay To validate the transcriptome data, ten DEGs were selected randomly for qRT-PCR assay, and 16S rRNA gene was used as the internal reference ([55] Table 1 ). qRT-PCR was using an ABI Stepone Plus Platform (Thermo, USA). Each gene was analyzed in three biological samples, and three reaction repeats were performed for each biological sample as described previously. Table 1. Genes and primers used in the qRT-PCR assay. Gene ID Gene name Forward Reverse Fragment size (bp) TM [56]AVV84311.1 rpsD CTAACACACCGTAAATACGACGAA TAAACTGGAAACTGCACCTGGA 110 60 [57]AVV84994.1 GAGTGGTAATAAGGTTGGCGTC GGTGTATCTGGGCAAGTAGGGT 250 60 [58]AVV84995.1 GTAACTCACCCATACCGGAAATAA CCCAAGTCTAAAGCAGACCAAG 130 60 [59]AVV82000.1 sucA TGAAGCGGTTGCTTTTGTGT ATAGATCTTACGTGGTGTAGGGTGT 181 60 [60]AVV84336.1 fusA GGAAAAACGCCGTAAAGAAAA GTTGAAAATCAGCCAAAGCAA 107 60 [61]AVV85624.1 gabD TAGATGATGTACAGACCTGTCCCG ACTTTCCCTACTACCAAATTGCG 405 60 [62]AVV86072.1 TTCTGTTATACCCGCTTTGCTTT CTGTTTAGTCTGTCACGGTTCTGT 298 60 [63]AVV84898.1 speC AGAAGCCTGCTTGTTGTTTGTGT GGTTGATCGTATTGGTCATCTATGT 173 60 [64]AVV82314.1 uspE AGCATTATTAACCACGCCATC AATTCAGCATCTAACTGAGCAGC 497 60 [65]AVV84205.1 katG TCGAGCGTTTTAAATGCTTCG CATGGTGGTAATACCTCCGTCAC 143 60 16S CGGTGAATACGTTCYCGG GGWTACCTTGTTACGACTT 128 60 [66]Open in a new tab Statistical Analyses All the experiments were conducted in six independent biological replicates. Related data to biofilm formation were tested using Duncan’s multiple range test in SPSS 22.0 software (IBM, New York, USA). All data were reported as mean ± standard deviation. Differences with a p value < 0.05 were regarded as statistically significant. Results and Discussion Determination of Transcriptomes of Biofilm Cells Formed by S. putrefaciens WS13 Based on our prior research ([67]Yan and Xie, 2020), the biofilm of S. putrefaciens WS13 grown in LB medium (pH=7.2) reached maturity at 24 h and 84 h at 4°C and 30°C, respectively (Figures not shown). The cells of mature biofilm at both temperatures were collected, and transcriptomes at a global gene expression level were obtained using Illumina RNA sequencing technique for the further analysis. Identification of DEGs in S. putrefaciens WS13 Induced by the Cold Stress DEGs of the biofilm cells formed by S. putrefaciens WS13 at 4°C and 30°C were identified, and the results are shown in [68]Figure 1 . A total of 761 DEGs were discovered, among which the expression of 497 DEGs was significantly up-regulated, and 264 DEGs were significantly down-regulated (p <0.05). Figure 1. [69]Figure 1 [70]Open in a new tab The volcano plots of DEGs between biofilms formed by S. putrefaciens WS13 under 4°C and 30°C. All DEGs were classified into three major functional categories in the Gene Ontology (GO) database, including biological process (BP), cellular component (CC), and molecular function (MF). The GO enrichment analysis of the DEGs revealed that the most abundant GO function was the peptide metabolic process (11.51%, 64/556), followed by protein-containing complex subunit organization (10.61%, 59/556), translation (10.25%, 57/556), peptide biosynthetic process (10.25%, 57/556), and ion transport (9.89%, 55/556) in BP, whereas cellular respiration (5.04%, 28/556), ATP hydrolysis coupled transmembrane transport (2.52%, 14/556), and ATP hydrolysis coupled ion transmembrane transport (2.52%, 14/556) showed an opposite pattern ([71] Figure 2 ). The protein-containing complex (23.20%, 129/556) was the most enriched DEGs in CC, while the percentages of the DEGs in the structural constituent of ribosome (6.83%, 38/556), and structural molecule activity (6.83%, 38/556) was the highest, followed by the rRNA binding (3.60%, 20/556) in MF ([72] Figure 2 ). Figure 2. [73]Figure 2 [74]Open in a new tab GO enrichment analyses of DEGs. The KEGG pathway enrichment analysis was further performed on the identified DEGs in the obtained transcriptomes of S. putrefaciens biofilm cells, and the results revealed 28 significantly changed metabolic pathways, including ribosome, oxidative phosphorylation, citrate cycle, transporters, photosynthesis proteins, histidine metabolism, photosynthesis, tryptophan metabolism, pyruvate metabolism, glyoxylate and dicarboxylate metabolism, carbon fixation pathways in prokaryotes, lysine degradation, MAPK signaling pathway-plants, propanoate metabolism, butanoate metabolism, mitochondrial biogenesis, inositol phosphate metabolism, β-alanine metabolism, longevity regulating pathway-worm, amyotrophic lateral sclerosis, pathways of neurodegeneration-multiple diseases, translation factor, and others, antimicrobial resistance genes, membrane trafficking, phenylalanine metabolism, and longevity regulating pathway-multiple species ([75] Figure 3 ). Figure 3. [76]Figure 3 [77]Open in a new tab KEGG enrichment analysis of DEGs. Identified DEGs Involved in Biofilm Formation of S. putrefaciens WS13 at the Cold Stress Four significantly altered metabolic pathways associated with biofilm formation were found in S. putrefaciens WS13 at 4°C, including the biofilm formation, amino acid metabolism, two-component system (TCS), and quorum sensing (QS [78]Table 2 ). Table 2. The major DEGs in biofilm cells of S. putrefaciens WS13 induced by the cold stress. Matabolic pathway Gene ID Name Log[2]FC Definition Biofilm formation - Pseudomonas aeruginosa [79]AVV85673.1 -2.77 Hypothetical protein SPWS13_3995 [80]AVV82533.1 -3.75 Crp/Fnr family transcriptional regulator [81]AVV82532.1 -5.68 Crp/Fnr family transcriptional regulator [82]AVV85155.1 nrfG -2.66 Nitrite reductase [83]AVV83548.1 -1.96 GntR family transcriptional regulator [84]AVV85962.1 flgM -1.01 Anti-sigma-28 factor FlgM [85]AVV85989.1 flrA -1.01 Fis family transcriptional regulator [86]AVV85720.1 -1.10 Membrane protein [87]AVV85731.1 -1.87 Serine/threonine protein phosphatase Tryptophan metabolism [88]AVV82001.1 sucB 2.79 Dihydrolipoamide succinyltransferase [89]AVV85789.1 amiE 1.60 Amidase [90]AVV82915.1 katE 2.84 Catalase [91]AVV84271.1 5.68 Aldehyde dehydrogenase [92]AVV82993.1 1.69 Enoyl-CoA hydratase [93]AVV85788.1 amiE 2.36 Amidase [94]AVV85097.1 pdhD 1.53 Dihydrolipoamide dehydrogenase [95]AVV86128.1 fadJ 2.12 Multifunctional fatty acid oxidation complex subunit alpha [96]AVV82394.1 bkdB 1.15 Dihydrolipoamide acetyltransferase [97]AVV82994.1 1.17 Enoyl-CoA hydratase Histidine metabolism [98]AVV82365.1 -3.60 TonB-dependent receptor [99]AVV82366.1 -2.85 TonB-dependent receptor [100]AVV82367.1 -1.83 TonB-dependent receptor [101]AVV84622.1 hutI -2.21 Imidazolonepropionase [102]AVV84105.1 -1.24 Arginase [103]AVV84398.1 urdA -1.38 Cytochrome C [104]AVV82368.1 -1.34 TonB-dependent receptor [105]AVV84082.1 -1.21 Hypothetical protein SPWS13_2316 [106]AVV86163.1 hutH -1.08 Histidine ammonia-lyase [107]AVV82091.1 hisC 2.88 Histidinol-phosphate aminotransferase [108]AVV84271.1 5.68 Aldehyde dehydrogenase [109]AVV82088.1 2.89 Hypothetical protein SPWS13_0228 [110]AVV82089.1 hisH 2.58 Imidazole glycerol phosphate synthase [111]AVV82090.1 hisB 2.52 Imidazoleglycerol-phosphate dehydratase [112]AVV82092.1 hisD 1.84 Histidinol dehydrogenase [113]AVV82087.1 hisF 2.35 Imidazole glycerol phosphate synthase [114]AVV82086.1 hisI 1.94 Phosphoribosyl-AMP cyclohydrolase [115]AVV84178.1 hutH 1.00 Histidine ammonia-lyase [116]AVV82085.1 hisI 1.83 Phosphoribosyl-ATP pyrophosphatase Lysine degradation [117]AVV82001.1 sucB 2.79 Dihydrolipoamide succinyltransferase [118]AVV84271.1 5.68 Aldehyde dehydrogenase [119]AVV82993.1 1.69 Enoyl-CoA hydratase [120]AVV85097.1 pdhD 1.53 Dihydrolipoamide dehydrogenase [121]AVV86128.1 fadJ 2.12 Multifunctional fatty acid oxidation complex subunit alpha [122]AVV82394.1 bkdB 1.15 Dihydrolipoamide acetyltransferase [123]AVV82994.1 1.17 Enoyl-CoA hydratase Two-component system [124]AVV84741.1 nrfC -6.34 4Fe-4S ferredoxin 4Fe-4S ferredoxin [125]AVV85250.1 petA -9.12 Ubiquinol-cytochrome C reductase [126]AVV83964.1 -3.21 Chemotaxis protein [127]AVV82196.1 hydA -3.63 Quinone-reactive Ni/Fe-hydrogenase small chain [128]AVV82533.1 -3.75 Crp/Fnr family transcriptional regulator [129]AVV84073.1 -2.62 Histidine kinase [130]AVV83344.1 glnA -2.11 Glutamine synthetase [131]AVV82195.1 hydB -3.36 Hydrogenase 2 large subunit [132]AVV85596.1 -2.04 Chemotaxis protein [133]AVV83291.1 -1.60 Sensor histidine kinase [134]AVV83573.1 frdA -1.82 Fumarate reductase [135]AVV86289.1 cheV -1.50 Chemotaxis protein CheW [136]AVV85217.1 CheC -1.57 Chemotaxis protein CheC [137]AVV83574.1 frdA -1.66 Fumarate reductase [138]AVV82532.1 -5.68 Crp/Fnr family transcriptional regulator [139]AVV85476.1 norR -1.52 Transcriptional regulator [140]AVV84252.1 glnL -1.94 Nitrogen regulation protein NR(II) [141]AVV85586.1 -1.79 Amino acid ABC transporter substrate-binding protein [142]AVV82193.1 hyaC -2.37 Hydrogenase Ni/Fe-hydrogenase, b-type cytochrome subunit [143]AVV82194.1 hydB -2.43 Hydrogenase 2 large subunit [144]AVV84209.1 glnA -1.56 Glutamine synthetase [145]AVV83483.1 -1.28 Chemotaxis protein [146]AVV84098.1 -1.32 Chemotaxis protein [147]AVV83827.1 -1.14 Peptidase [148]AVV84924.1 -1.35 Cell division protein ZapB [149]AVV83129.1 motA -1.45 Flagellar motor protein PomA [150]AVV85069.1 frdC -1.56 Fumarate reductase [151]AVV82265.1 ampC -1.20 Beta-lactamase [152]AVV85068.1 frdC -1.15 Fumarate reductase [153]AVV82166.1 -1.00 Chemotaxis protein [154]AVV84289.1 -1.85 Histidine kinase [155]AVV84253.1 glnL -1.47 Nitrogen regulation protein NR(II) [156]AVV84398.1 urdA -1.38 Cytochrome C [157]AVV83971.1 psrB -1.96 Polysulfide reductase subunit B [158]AVV85962.1 flgM -1.01 Anti-sigma-28 factor FlgM [159]AVV84053.1 maeB -1.05 Malate dehydrogenase [160]AVV84712.1 cusB -2.01 RND transporter MFP subunit [161]AVV84254.1 glnL -1.32 Nitrogen regulation protein NR(II) [162]AVV84075.1 -1.54 Response regulator receiver protein [163]AVV83905.1 -1.15 Chemotaxis protein [164]AVV85989.1 flrA -1.01 Fis family transcriptional regulator [165]AVV85931.1 -1.09 Membrane protein cyd operon protein YbgT [166]AVV83697.1 glnK -1.66 Nitrogen regulatory protein P-II 1 [167]AVV85720.1 -1.10 Membrane protein [168]AVV85990.1 flrB -1.38 Sensor histidine kinase [169]AVV83549.1 cydA -2.29 Cytochrome D ubiquinol oxidase subunit I [170]AVV85070.1 frdA -1.23 Fumarate reductase [171]AVV84287.1 frdC -4.57 Fumarate reductase [172]AVV83377.1 glnB -1.11 Nitrogen regulatory protein P-II [173]AVV84286.1 frdB -4.81 Fumarate reductase iron-sulfur subunit [174]AVV83550.1 cydB -1.80 Ubiquinol oxidase subunit II, cyanide insensitive [175]AVV83918.1 ttrB -1.97 Tetrathionate reductase subunit B [176]AVV84284.1 frdA -4.55 Fumarate reductase flavoprotein subunit [177]AVV86161.1 -2.76 Cytochrome C flavocytochrome c [178]AVV85991.1 flrC -1.11 Fis family transcriptional regulator [179]AVV84288.1 frdD -3.58 Fumarate reductase [180]AVV84265.1 cpxA -1.06 Sensor histidine kinase [181]AVV83919.1 ttrC -1.95 Polysulfide reductase [182]AVV84285.1 frdA -1.65 Fumarate reductase flavoprotein subunit Quorum sensing [183]AVV84313.1 secY 2.27 Preprotein translocase subunit SecY [184]AVV86104.1 secF 2.06 Preprotein translocase subunit SecF [185]AVV84314.1 secY 2.13 Preprotein translocase, SecY subunit [186]AVV84540.1 yidC 1.61 Membrane protein insertase [187]AVV85368.1 2.07 RND transporter [188]AVV84996.1 1.41 Outer membrane adhesin-like protein [189]AVV83999.1 1.41 Hypothetical protein SPWS13_2218 [190]AVV86108.1 4.64 Peptidase S8 [191]AVV82086.1 hisI 1.94 Phosphoribosyl-AMP cyclohydrolase [192]AVV82884.1 dppB 1.22 Peptide ABC transporter permease [193]AVV82009.1 1.72 Hypothetical protein SPWS13_0149 [194]AVV86102.1 yajC 1.14 Preprotein translocase subunit YajC [195]AVV86191.1 trpG 1.13 Anthranilate synthase subunit II [196]AVV82085.1 hisI 1.83 Phosphoribosyl-ATP pyrophosphatase [197]AVV84384.1 1.11 Cytochrome B561 [198]AVV82919.1 gadB 1.76 Glutamate decarboxylase Carbohydrate metabolism [199]AVV83139.1 ykgG -3.15 L-lactate dehydrogenase complex protein LldG [200]AVV83137.1 Nan -3.32 L-lactate dehydrogenase complex protein LldE [201]AVV83138.1 Nan -3.01 L-lactate dehydrogenase complex protein LldF [202]AVV85349.1 Nan -1.53 Acetoin utilization protein AcuB [203]AVV85417.1 sfsA -1.26 Sugar fermentation stimulation protein A Energy metabolism [204]AVV84740.1 nrfB -6.79 Cytochrome c nitrite reductase small subunit [205]AVV84741.1 nrfC -6.35 Polysulfide reductase chain B [206]AVV84742.1 nrfD -5.78 Protein NrfD [207]AVV85151.1 ccmF -3.50 Cytochrome c-type biogenesis protein NrfE [208]AVV85149.1 ccmH -3.58 Formate-dependent nitrite reductase complex subunit NrfG [209]AVV85155.1 nrfG -2.67 Formate-dependent nitrite reductase complex subunit NrfG [210]AVV83971.1 psrB -1.96 Polysulfide reductase chain B [211]AVV82007.1 Nan -1.06 Ferredoxin/flavodoxin—NADP+ reductase [212]AVV83970.1 psrC -1.72 Polysulfide reductase chain C [213]AVV85150.1 ccmF -3.15 Cytochrome c-type biogenesis protein NrfE [214]Open in a new tab Some genes encoding transcripitional regulators were slightly down-regulated (p<0.05), such as flgM, and flrA genes. The former that encodes an anti-sigma-28 factor FlgM can regulate the expression of flagellar genes in a complex regulatory network controlling chemotaxis, swimming and biofilm formation in Rhodobacter sphaeroides ([215]Wilkinson et al., 2011). It has been reported that the flrA gene that encodes Fis family transcriptional regulator was strongly sensitive to environmental stress ([216]Gang et al., 2016). The silencing of the flrA gene led to deficiencies in adhesion, motility, flagellar assembly, biofilm formation and exopolysaccharide (EPS) production in Vibrio alginolyticus ([217]Gang et al., 2016). These results suggested that the reduction of flagellar synthesis and motility may have enhanced the biofilm formation of S. putrefaciens WS13 at 4°C. Identified DEGs Involved in Amino Acid Metabolism in S. putrefaciens Biofilm at the Cold Stress The biofilm is composed of extracellular substances secreted by S. putrefaciens WS13, such as polysaccharides, proteins, lipids, and other substances. Proteins play a crucial role in maintaining the structural stability of biofilms ([218]Yan and Xie, 2020). Amino acids are key intermediates in both carbon and nitrogen metabolisms in microorganisms. Bacterial amino acid metabolism is usually sensitive to environmental stress. In this study, comparative transtriptomics analyses revealed remarkedly changed DEGs in the amino acid metabolism in biofilm cells of S. putrefaciens WS13 under the cold stress. These DEGs were significantly enriched in the amide biosynthetic process, peptide metabolic process, translation, peptide biosynthetic process, and other protein-related GO functions. For example, interestingly, all the DEGs involved in the tryptophan metabolism were significantly up-regulated in the S. putrefaciens WS13 biofilm at 30°C, including the sucB, amiE, katE, pdhD, fadJ, and bkdB genes (2.23 fold to 7.15 fold) (p<0.05). It has been reported that the tryptophan biosynthesis genes were up-regulated in the biofilms of Escherichia coli and Salmonella enterica ([219]Domka et al., 2007; [220]Hamilton et al., 2009). It was shown that exogenous tryptophan significantly accelerated the biofilm formation of S. enterica and Fusobacterium nucleatum and completely restored the deleted mutant of S. enterica in biofilm formation ([221]Hamilton et al., 2009; [222]Sasaki-Imamura et al., 2010). In the amino acid metabolism, some DEGs involved in histidine metabolism were significantly up-regulated in S. putrefaciens WS13 at 4°C (2.23 fold to 61.60 fold) (p<0.05). It has been reported that under the acetate, butyrate, or butanol stress, the genes involved in histidine biosynthesis were up-regulated in C. acetobutylicum ([223]Alsaker et al., 2010; [224]Wang et al., 2013). The induction of histidine biosynthesis genes was also observed under the acid tolerance in Lactobacillus casei ([225]Broadbent et al., 2010). Amino acid metabolism played a vital role in bacterial adaptation to certain circumstances including metabolite stress ([226]Alsaker et al., 2010; [227]Wang et al., 2013), oxygen tolerance ([228]Hillmann et al., 2009), and sporulation ([229]Jones et al., 2008). The highly increased amino acid metabolism may enable bacterial cells to adjust the structure and function of biofilm in response to the cold stress. Identified DEGs Involved in TCSs in S. putrefaciens Biofilm at the Cold Stress Bacteria have a variety of signal transduction systems, which can sense external signal stimuli and respond adaptively to changes (e.g., osmolarity, light, temperature, and oxygen) in the surrounding environment ([230]Liu et al., 2013). TCS that compose of histidine protein kinases (HKs) and response regulators (RRs) are widely present in gram-negative bacteria. TCSs are very important signaling pathways that coordinate responses to environmental stimulus, and regulates bacterial sporulation, biofilm formation, competence, and chemotaxis ([231]Sieuwerts et al., 2010). In this study, comparative transcriptomic analyses revealed 59 DEGs involved in the TCSs, which were all significantly up-regulated in S. putrefaciens WS13 at 4°C, e.g., cheV (2.84 fold), motA (2.75 fold). For example, bacterial chemotaxis system is a typical coupling protein-dependent signal transduction system and play a crucial role in bacterial colonigetion and adhesion. It has been reported that the up-regulated genes involved in chemotaxis significantly reduced the ability of adhesion of the organism, motility, chemotaxis, and biofilm formation in Vibrio harveyi. Environmenal factors such as temperatures, salinities, and pH values affected the chemotactic gene expression involved in the regulation of adhesion ability ([232]Xu et al., 2021). In this study, the cheW and cheC genes were significantly down-regulated in the biofilm formation by S. putrefaciens at the cold stress. The cheV gene encodes a linker protein, while the cheC encoding protein has phosphatase activity ([233]Moon et al., 2016). Moreover, expression of six additional genes coding for chemotaxis proteins were all reduced at the transcriptional level in S. putrefaciens at 4°C, including the [234]AVV85596.1 (0.24 fold), [235]AVV83483.1 (0.41 fold), [236]AVV84098.1 (0.40 fold), and [237]AVV83905.1 (0.45 fold). These results suggested inhibital chemotaxis and/or motility of the bacterium. The movement of flagella via the flagellar motor complex affected mature biofilm architecture ([238]Wood et al., 2006). In this study, the motA gene that encodes a key component of the motor compelx was significantly down-regulated (-1.45 fold) (p<0.05) suggesting poor motility forming flatter microcolony structures of S. putrefaciens biofilm at 4°C. Additionally, the DEGs enriched in TCSs may sense external signal stimulation and regulate the movement ability of microbial cells, thereby affecting the biofilm formation. Identified DEGs Involved in QS in S. putrefaciens Biofilm at the Cold Stress Biofilm formation is essentially coordinated through a cell density-dependent gene regulation system known as QS ([239]Fuqua et al., 1994; [240]Kjelleberg and Molin, 2002; [241]Parsek and Greenberg, 2005). In this study, remarkly, the DEGs involved in the QS were significantly down-regulated in S. putrefaciens at 4°C, including secY (4.83 fold), secF (4.20 fold), hisI (3.83 fold), gadB (3.38 fold), yidC (3.06 fold), dppB (2.34 fold), yajC (2.21 fold), trpG (2.19 fold). For example, YidC has been recognized as a drug to inhibit biofilm formation in Staphylococcus aureus ([242]Hofbauer et al., 2018). Expression of the yidC gene was significantly influenced bu pH and starvation stress in Vibrio alginolyticus, and the bacterial adhesion was significantly decreased after silencing of the yidC gene ([243]Lemos et al., 2004). YidC function as integral membrane chaperone/insertase associated with the classical SecYEG translocon, which could contribute to inhibit biofilm formation ([244]Hofbauer et al., 2018). The upregulated expression of the secY and yidC genes involved in extracellular polymeric substances was also observed in the stage of biofilm maturation of Bifidobacterium Longun FGSZY16M3 ([245]Liu et al., 2021). In this study, the secF gene expression was also highly increased by 2.06 fold, which encodes a component of the Sec translocon. Additionally, expression of the yajC (1.14 fold), dppB (1.22fold), and [246]AVV84996.1 (1.41 fold) were significantly increased in S. putrefaciens WS13 biofilm at the cold stress, which encode a potein translocase subunit YajC, a peptide ABC transporter permease, and an outer membrane adhesin-like protein, respectively. These results suggested QS plays an important role in the biofilm formation of S. putrefaciens at the cold stress. Identified DEGs Involved in the Other Key Metabolic Pathways in S. putrefaciens Biofilm at the Cold Stress Energy metabolism such as carbohydrate metabolism, pyruvate metabolism, oxidative phosphorylation, and citrate cycle helps bacteria adapt to the changing environment. In this study, approximately 5 DEGs involved in the carbohydrate metabolism were significantly down-regulated (0.10 fold to 0.41 fold) in S. putrefaciens WS13 biofilm formed at 30°C (p<0.05). For energy metabolism, approximately 10 DEGs were significantly down-regulated (0.01 fold to 0.48 fold) at 30°C. For example, expression of the ppsA gene was increased by 2.62 fold in S. putrefaciens at 4°C which can promote biofilm formation by enhancing bacterial adhesion ([247]Ao et al., 2020). Moreover, approximately 11 DEGs involved in aminoacyl-tRNA biosynthesis were highly down-regulated (2.02 fold to 5.15 fold) in S. putrefaciens WS13 biofilm at 4°C, which may have released amino acids to feed the energy-providing pathways ([248]Bénédicte et al., 2009), and benefited the survival of S. putrefaciens WS13 under the cold stress. Confirmation of the DEGs by the qRT-PCR Assay The relative expression levels of randomly selected DEGs were determined and calculated using 16S rRNA as the internal reference gene, including: [249]AVV84311.1, [250]AVV84994.1, [251]AVV84995.1, [252]AVV82000.1, [253]AVV84336.1, [254]AVV85624.1, [255]AVV86072.1, [256]AVV84898.1, [257]AVV82314.1, and [258]AVV84205.1. The obtained qRT-PCR results ([259] Table 3 ) confirmed the transcription changes of these DEGs in the comparative transcriptome analyses. Table 3. The results of the qRT-PCR assay. Index Sample Name Assay Name mean ΔCT ΔΔCT Fold Change Up/Down Control 1 S30 [260]AVV84311.1 -11.36 -2.87 7.29 Up 16S WS4 [261]AVV84311.1 -8.49 2 S30 [262]AVV84994.1 -8.63 -5.23 37.46 Up 16S WS4 [263]AVV84994.1 -3.40 3 S30 [264]AVV84995.1 -10.56 -3.96 15.52 Up 16S WS4 [265]AVV84995.1 -6.60 4 S30 [266]AVV82000.1 -8.12 -4.33 20.06 Up 16S WS4 [267]AVV82000.1 -3.79 5 S30 [268]AVV84336.1 -8.30 -3.54 11.67 Up 16S WS4 [269]AVV84336.1 -4.75 6 S30 [270]AVV85624.1 -4.72 2.30 0.20 Down 16S WS4 [271]AVV85624.1 -7.02 7 S30 [272]AVV86072.1 -5.59 2.57 0.17 Down 16S WS4 [273]AVV86072.1 -8.16 8 S30 [274]AVV84898.1 -5.43 5.29 0.03 Down 16S WS4 [275]AVV84898.1 -10.72 9 S30 [276]AVV82314.1 -3.06 2.26 0.21 Down 16S WS4 [277]AVV82314.1 -5.32 10 S30 [278]AVV84205.1 -5.54 1.85 0.28 Down 16S WS4 [279]AVV84205.1 -7.39 [280]Open in a new tab Conclusions This study was the first to characterize the global-level gene expression of biofilm cells of S. putrefaciens WS13 under the cold stress. Distinct transcriptomic profiles were obtained using Illumina RNA sequencing technique. Comparative transcriptomic analyses revealed a total of 761 DEGs in the biofilm formed at 4°C, among which the expression of 497 DEGs was significantly up-regulated, and 264 DEGs were significantly inhibited (p<0.05). Although carbohydrate and energy metabolisms were repressed in the biofilm cells at the harsh condition, S. putrefaciens WS13 reduced chemotaxis, and mobility, but enhanced histidine biosynthesis, tryptophan metabolism, and QS to construct the biofilm and survive at the cold stress. This work provides valuable insights into the transcriptiomic regulation in biofilm cells under cold stress and laid a theoretical foundation for the targeted inhibition of the biofilm formation of the severe spoilage S. putrefaciens WS13. Data Availability Statement Raw data of the transcriptomes was deposited in NCBI database under the accession number PRJNA759975 ([281]http://www.ncbi.nlm.nih.gov/bioproject/759975). Author Contributions JX, JY conceived the idea. JY carried out the laboratory work and wrote the paper. ZY assisted in the data analysis. JX revised the manuscript. All authors have read and approved the final manuscript. Funding This work was supported by the National Natural Science Foundation of China (31972142), Shanghai Engineering Research Center Construction Special Fund from Shanghai Municipal Science and Technology Commission (19DZ2284000), the Startup Foundation for Young Teachers of Shanghai Ocean University. Conflict of Interest The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest. Publisher’s Note All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher. References