- Research article
- Open Access
Metabolic characterization of the chitinolytic bacterium Serratia marcescens using a genome-scale metabolic model
BMC Bioinformatics volume 20, Article number: 227 (2019)
Serratia marcescens is a chitinolytic bacterium that can potentially be used for consolidated bioprocessing to convert chitin to value-added chemicals. Currently, S. marcescens is poorly characterized and studies on intracellular metabolic and regulatory mechanisms would expedite development of bioprocessing applications.
In this study, our goal was to characterize the metabolic profile of S. marcescens to provide insight for metabolic engineering applications and fundamental biological studies. Hereby, we constructed a constraint-based genome-scale metabolic model (iSR929) including 929 genes, 1185 reactions and 1164 metabolites based on genomic annotation of S. marcescens Db11. The model was tested by comparing model predictions with experimental data and analyzed to identify essential aspects of the metabolic network (e.g. 138 essential genes predicted). The model iSR929 was refined by integrating RNAseq data of S. marcescens growth on three different carbon sources (glucose, N-acetylglucosamine, and glycerol). Significant differences in TCA cycle utilization were found for growth on the different carbon substrates, For example, for growth on N-acetylglucosamine, S. marcescens exhibits high pentose phosphate pathway activity and nucleotide synthesis but low activity of the TCA cycle.
Our results show that S. marcescens model iSR929 can provide reasonable predictions and can be constrained to fit with experimental values. Thus, our model may be used to guide strain designs for metabolic engineering to produce chemicals such as 2,3-butanediol, N-acetylneuraminic acid, and n-butanol using S. marcescens.
With recent advances in high-throughput data and bioprospecting, the breadth of novel and interesting biochemistry continues to expand. In terms of the large amount of data and resources available, one current challenge is to interpret raw data into knowledge that can help better understand complex biological systems. Genome-scale metabolic modeling (GSMM) provides gene-protein-reaction level specificity and can combine metabolic network with fundamental genomic and biochemical information. However, there are still a number of limitations . For instance, because metabolic pathway fluxes are underdetermined, which means always existence of alternative flux states according to different pathway usage that generate an indistinguishable physiological phenotype. This is an underlying issue with GSMMs that limits the level of detailed predictions these models can make when running simulations (presence of alternate optimal solutions) if not additional data/information is used. In this study, we investigate the poorly-characterized chitinolytic bacterium, Serratia marcescens, and utilize RNAseq data to achieve a better understanding of its metabolic network.
S. marcescens is unique among enteric bacteria in many aspects. It secretes extracellular DNase, gelatinase, lipase, several proteases, a red pigment (prodigiosin), chitinases and a chitin binding protein. It is believed to be one of the most efficient chitin-degrading bacteria in the environment [2, 3]. S. marcescens Db11 contains ten chitinase-related proteins  leading to studies focused on various facets of the S. marcescens chitinolytic mechanisms [5,6,7,8]. Due to the high efficiency of processing chitin, several individual S. marcescens chitinase genes have been cloned into model bacterial species (e.g. Escherichia coli) [9,10,11, 46, 47]. The cloned enzymes were isolated in good concentrations but failed to show similar level of chitinolytic activity as is found in S. marcescens. This may be due to the complexity of chitin degradation systems and synergy between multiple enzymes [4, 8, 12, 13]. Hence, it is suggested that the best way to fully utilize the chitinolytic capabilities of S. marcescens may be developing S. marcescens rather than moving its chitinases into other systems by heterologous expression.
The ability to produce chemicals of industrial importance using inexpensive chitinolytic biomass has been a recent focus [14, 49, 50]. Microbial conversion of chitin waste into value-added chemicals, a similar concept as consolidated bioprocessing (CBP), can be socially and economically beneficial [14,15,16, 48]. For S. marcescens Db11, the sequencing of its genome in 2014 sets up a milestone towards understanding this industrially applicable microbe . The initial developments in the characterization of S. marcescens are listed in Table 1. This is a promising development toward making use of the chitinolytic capabilities of this microbe to reduce the complex multi-step bioprocess to CBP.
In studying S. marcescens Db11, we believe that it can be developed into a facile chitinolytic system for consolidated bioprocessing. Recently, we established molecular tools for genetic engineering this organism using: 1) an in-frame gene deletion approach based on endogenous exonucleases by introducing linear DNA fragments and 2) a shuttle vector that contains functional replication origin and expression elements (e.g. promoter and RBS) . These approaches allowed us to study the function of a chitinase regulator protein (ChiR) by characterizing physiological differences in a ChiR overeexpression strain and a ChiR deletion strain. ChiR is characterized as a positive regulator protein and is essential for chitinase production in S. marcescens. In the same study, the ChiR overexpression S. marcescens produced 1.13 g/L 2,3-butanediol from 2% crystal chitin. The ability to convert chitin into biochemicals makes this organism an interesting candidate for a chitin-based CBP [15, 18, 49]. Thus, our group aims to present a system level understanding of the metabolic network of S. marcescens to facilitate future engineering applications.
A reconstructed GSMM can provide a framework to understand cellular processes and can be used to find target for focused metabolic engineering to yield products of biotechnological value . The most widely-used algorithms for design and simulation of genome-scale constraint-based metabolic models such as OptKnock , OptForce , EMiLio  are based on flux distribution and flow through chemically balanced reactions (FBA). Flux balance analysis (FBA) uses linear algorithm to optimize an objective function:
where, Z is the flux towards objective function, e.g. biomass production and production optimization, S is stoichiometry of the reactions as matrix, v is the reaction flux vector, ai and bi are the constraints based on the flux vi of the reaction i .
After compiling biochemical information into a draft model, a second step of model curation can be done by integrating high-throughput omics data to realize “content in context” . A variety of applications using incorporating omics data with the framework of GSMM have been reported for many prokaryotic microbes to more closely match cellular processes [24,25,26]. In this study, we aim at understanding metabolite profile of S. marcescens Db11 and its chitinolytic regulation systems. Therefore, a genome-scale metabolic model, iSR929, consisting of 1185 reactions, 1164 metabolites and 929 genes was constructed. The constructed model was analyzed for 1) accuracy compared to experimental results; 2) prediction of essential genes; 3) metabolic differences under different carbon growth conditions. We also proposed three potential chemical compounds for metabolic engineering implementation.
Construction of S. marcescens GSMM
Based upon the genome sequence of S. marcescens Db11  and available physiological evidence, a genome-scale constraint-based metabolic model for S. marcescens, hereafter denoted iSR929, was constructed. The reaction database used for drafting the model includes but is not limited to sources such as KEGG , BiGG , rBioNet , UniProt  and MBRole . iSR929 contains 1185 reactions representing the function of 929 genes and 1164 metabolites (see Table 2). Metabolic gaps were identified and addressed guided by a computational gap-filling algorithm [19, 24, 25, 32, 33]. Among the 1100 gene associated reactions of the total 1185 reactions, there are 795 reactions that are associated with only one reaction and there are 305 reactions associated with more than one enzyme, meaning either isozymes or enzyme complexes. Of the 1164 metabolites, 1099 metabolites are intracellular, 43 metabolites are extracellular, and 22 metabolites are boundary. We have added gene-protein-reactions (ChiA, ChiB, ChiC, SMDB11_1994, Chb, NagZ, SMDB11_1542, and SMDB11_4602) that lead to converting chitin into N-acetylglucosamine. In addition, we have added ChiB as a deacetylase that can convert N-acetylglucosamine to glucosamine, and added NagZ as a glucosamine kinase that phosphorylates glucosamine to glucosamine-6-P. Thus, the model can indicate the biotransformation of chitin to glucosamine, and glucosamine can be a carbon source that enters into carbon central metabolism. Furthermore, we added chitobiose degradation reactions to the model based upon the annotation and experimental evidence (EC 22.214.171.124, SMDB11_0477, SMDB11_1190 SMDB11_1542 and SMDB11_4602).
The breakdown of iSR929 by functional categories is represented in Fig. 1. Amino acid metabolism was the largest subsystem with 253 reactions. Experimentally, S. marcescens appears to have the capacity to synthesize all 20 amino acids , and iSR929 reflects this. Other large groups of reactions include subsystems related to carbohydrate metabolism (glycolysis, pentose phosphate pathway, pyruvate metabolism, and TCA cycle), cofactor and vitamin metabolism (nicotinate and nicotinamide metabolism, folate biosynthesis, and porphyrin and chlorophyll metabolism) and nucleotide metabolism (including reactions related to purine and pyrimidine synthesis). The model also includes reactions relating to the synthesis of lipids, including fatty acid, glycerophospholipid, and glycerolipid metabolism.
To test the predictions of the iSR929 model, we simulated growth of S. marcescens by applying FBA, assuming minimal media conditions with one of three possible carbon sources (glucose, N-acetylglucosamine, or glycerol). We then compared model predictions to experimentally observed growth rates and fermentation product secretion profiles of S. marcescens grown in batch culture. In addition to the carbon source, the in silico minimal medium used for simulations contained water (H2O), ammonia (NH4), sulfate (SO4), phosphate (pi), calcium (Ca2), ferrous iron (Fe3), hydrogen sulfide (H2S), potassium (K), magnesium (Mg2), pantothenate (pnto-r), and nicotinate D-ribonucleotide (nmn). In each of the three simulation conditions (glucose, N-acetylglucosamine, and glycerol), we applied progressively more experimentally determined constraints associated with by-product secretion rates to determine how closely the computational results could match the experimental results given the possibility of alternate optimal solutions. Figure 2 shows the simulation results for each growth condition when exchange rates for the carbon sources, acetic acid, and succinic acid were constrained to match experimental observations.
Computational prediction of gene knockout targets
A comprehensive analysis of in silico single gene deletions were conducted with iSR929 using glucose as a carbon source and the other minimal media components (as described above), and constraining the glucose uptake rate to its experimentally observed value for batch growth (10.53 mmol/gDW/h). Gene essentiality results are shown in Fig. 1. In the case of growth on glucose, we found that 138 (14.9%) of S. marcescens genes included in iSR929 were predicted to be essential. We also examined which subsystems of iSR929 contained the highest percentage of essential reactions (‘vulnerable subsystems’). Among the most vulnerable subsystems are amino acid metabolism, nucleotide metabolism, lipid metabolism, and metabolism of cofactors and vitamins.
mRNA sequencing of S. marcescens under growth on different carbon sources
mRNA samples of S. marcescens Db11 were prepared from cell cultures under four conditions: M9 medium with glucose, M9 medium with N-acetylglucosamine, M9 medium with glycerol, and LB medium. Overall, the average coverage per gene (4831 genes) (see Additional file 1: Table S1) implies a high-depth-coverage (generally 250).
Integrating RNASeq data into the GSMM
The integration of RNAseq data of S. marcescens grown in M9 medium with glucose, N-acetylglucosamine, or glycerol was conducted using the iSR929 model and a mixed integer linear programming algorithm approach . In order to establish a cutoff value of gamma (lower bound cutoff), we initially categorized gene expression level into 58 cutoffs based on an exponential of 1.2, and a distribution of gene numbers against each gene expression level range was then plotted, a representative figure of M9 glucose condition was shown in Additional file 1: Figure S1. The lower bound is 3.00, which excludes around 10% of the total number of genomic genes across all the conditions, and the upper bound is 850.56, shown in Additional file 1: Table S2.
Simulations were then run on the S. marcescens model after transcriptomic data integration for the three growth conditions (M9 glucose, M9 N-acetlyglucosamine, M9 glycerol) separately, where a reaction flux is computationally constrained to 1000 if the corresponding gene’s RPKM is over 850.56 and a reaction flux is 0 if the gene’s RPKM is less than 3.00. Raw data of reaction flux can be found in Additional file 1: Table S3 and a summary of the active metabolic reactions for the three growth conditions were categorized in different metabolic pathway modules (Table 3). When growing on glucose, there were 150 active reactions. Among these 150 reactions, the most highly represented functional categories were: carbohydrate metabolism (42), amino acid metabolism (38), lipid metabolism (12), cofactor metabolism (12), and nucleic acid metabolism (6). When growing on N-acetylglucosamine, there were 132 active reactions. Among these 132 active reactions, the majority of the reactions were: carbohydrate metabolism (39), amino acid metabolism (28), lipid metabolism (6), cofactor metabolism (9), and nucleic acid metabolism (10). When growing on glycerol, there were 146 active reactions. Among these 146 active reactions, the functional categories were: carbohydrate metabolism (44), amino acid metabolism (41), lipid metabolism (6), cofactor metabolism (8), and nucleic acid metabolism (7).
The main active pathways for the metabolism of glucose as a sole carbon source were shown in Fig. 3a. After transport, glucose input is routed through glycolysis and pentose phosphate pathway to yield pyruvate. Then, the main flux goes to the TCA cycle, butanoate metabolism (yielding 2,3-butanediol), amino acid synthesis, lipoprotein synthesis, and pyruvate metabolism (yielding ethanol, formate, lactate, and acetate).
The main active pathway for the metabolism of N-acetylglucosamine as a sole carbon source was shown in Fig. 3b. For the metabolism of N-acetylglucosamine, S. marcescens is believed to phosphorylate transported N-acetylglucosamine by a membrane type II phosphotransferase (SMDB11_0473). The SMDB11_0473 gene expression level in the N-acetylglucosamine condition (RPKM 967.20) was 12.09- and 5.13-fold than those under glucose (RPKM 80.23) and glycerol (RPKM 188.36) conditions, respectively, supporting that N-acetylgluosamine is likely phosphorylated intracellular. The carbon flows through glycolysis metabolism from fructose-6P to regenerate glucose-6P. The glucose-6P flows to the pentose phosphate pathway to yield pyruvate. Then, the majority of the carbon flux was distributed into purine and pyrimidine metabolism, butanoate metabolism (secreted 2,3-butanediol as an end-product), pyruvate metabolism (yielding ethanol, formate, lactate, and acetate as end-products), amino acid synthesis, and lipoprotein synthesis.
The main active pathways for the metabolism of glycerol as a sole carbon source were shown in Fig. 3c. Glycerol metabolism is thought to be converted to pyruvate by three steps: a membrane facilitator protein (SMDB11_4021), an ATP-dependent glycerol kinase (SMDB11_4022), and a glycerol-3-phosphate dehydrogenase (SMDB11_3886) . The SMDB11_4021 expression level in the glycerol growth condition (RPKM 1772.81) was 89.5- and 70.8-fold higher than those of glucose and N-acetylglucosamine condition. The SMDB11_4022 expression level in the glycerol condition (RPKM 4535.25) was 76.6- and 63.6-fold higher than those of glucose (RPKM 59.20) and N-acetylglucosamine (RPKM 71.28) condition. The gene expression level of SMDB11_3886 in the glycerol condition (RPKM 1683.89) was 59.9- and 54.7-fold higher than those of glucose (RPKM 28.10) and N-acetylglucosamine (RPKM 30.78) condition. The main active flux goes to glycolysis, pentose phosphate pathway, amino acid synthesis, and pyruvate metabolism.
S. marcescens exhibits low activity in TCA cycle with N-acetylglucosamine as a sole carbon source
According to model simulation results, we have examined metabolic fluxes differences in the carbon central metabolism, namely, glycolysis, pentose phosphate pathway, and TCA cycles, shown in Additional file 1: Figure S2A, S2B, and Fig. 4, respectively. The main metabolic differences among the three carbon sources are found in the TCA cycle. When glucose is the sole carbon source, the TCA cycle is fully active; under glycerol condition, the TCA cycle is partially active (R00361_smac, R01082_smac, R02164_smac, R00405_smac, R02570_smac, and R00621_smac); while under N-acetylglucosamine condition, the TCA cycle is not active. The gene expression level of the TCA cycle pathway under N-acetylglucosamine, glucose or glycerol condition were mapped in Fig. 4a. In general, the mRNA level of the genes in the TCA cycle support the model prediction results: compared to mRNA level of the TCA cycle genes under N-acetylglucosamine condition, those corresponding genes under glucose or glycerol condition are generally higher. In particular, the SMDB11_0505 gene, which encodes a citric acid synthase that converts acetyl-CoA to citric acid, may play the key role in determining the model prediction since under glucose or glycerol condition, the SMDB11_0505 mRNA level are both above the high threshold, while under N-acetylglucosamine condition, the SMDB11_0505 mRNA level is below the high threshold.
To test for differences in TCA cycle flux, measurement of intracellular citric acid and isocitric acid levels can be used as indicators of TCA cycle activity . Experimental measurement of intracellular citric acid and isocitric acid concentrations for growth on N-acetylglucosamine were lower than those for growth on glucose or glycerol, providing some support for the predicted lower TCA cycle flux results (Fig. 4b and c). Additionally, when N-acetylglucosamine is the sole carbon source, S. marcescens appears to utilize intermediates (e.g. oxaloacetate, fumarate, and 2-oxoglutarate) as precursors for amino acids synthesis, and synthesis of nucleotides (e.g. purine and pyrimidine) comes from the ribose-5P generated from the pentose phosphate pathways since N-acetylglucosamine is an ideal carbon source for amino-sugar and nucleotides synthesis (e.g. C:N ratio 8:1).
Cofactor metabolism and energy ATP metabolism
Based on metabolic flux analysis predicted by the model, the common consumption/re-generation of NADH under the three carbon source conditions are supported in the amino acid metabolism, butanoate metabolism and pentose phosphate metabolism; the common consumption/re-generation of NADPH under the three carbon source conditions are supported from the glycolysis and gluconeogenesis metabolism, one carbon pool by folate, and transport; common pathways associated with ATP consumption and generation under the three-carbon source conditions are supported from amino acid metabolism, one carbon pool by folate, pyruvate metabolism and transportation of molecules, shown in Additional file 1: Table S4.
When growing S. marcescens under N-acetylglucosamine as a sole carbon source, NAD+/NADH tends to be a main cofactor compared to NADP/NADPH. It is worth noting that under N-acetylglucosamine as a sole carbon source condition, additional NADH cylces can be achieved in the pyrimidine metabolism, and pentose and glucuronate interventions by a nucleotide phosphorylase and a 3-dehydro-L-gulonate 2-dehydrogenase. This is in agreement with the observance shown in the carbon metabolism map that nucleotide metabolism and pentose phosphate pathways are highly active when using N-acetylglucosamine. Particular ATP consumption and generation pathways in S. marcescens under N-acetylglucosamine condition include glycine, serine and threonine metabolism and pyrimidine metabolism.
Metabolic engineering applications of S. marcescens using the model iSR929
S. marcescens possesses interesting native metabolic capabilities (active chitinolytic system and native production of 2,3-butanediol) and can potentially be engineered for biochemical production. Based on the aforementioned experimental observations and computational simulation results, we selected three chemical compounds (2,3-butanediol, n-butanol, and N-acetylneuraminic acid) as potential production targets and illustrate pathway designs for the production of each.
S. marcescens Db11 is a native 2,3-butanediol producer. Direct production of 2,3-butanediol from crystal chitin by S. marcescens Db11 have been experimentally validated and the model iSR929 shows 2,3-butanediol production across all growth conditions . The proposed pathways for production of 2,3-butanediol and the corresponding annotated genes are shown in Fig. 5a. Unlike other reported S. marcescens 2,3-butanediol producers [37, 38], S. marcescens Db11 harbors only two 2,3-butanediol dehydrogenases (meso and (2S,3S)) and relatively small amount of (2S,3S)-2,3-butanediol dehydrogenase was expressed compared to the meso-2,3-butanediol dehydrogenase, indicating meso-2,3-butanediol is the major product.
The second target of interest is a biofuel chemical, n-butanol. The iSR929 incorporates most of the reactions present in the butanoate metabolism as inferred from genomic and biochemical data. Besides the 2,3-butanediol synthetic pathway, two active reactions (R01171 and R01975) were found under glucose growth conditions and one active reaction (R01171) was observed under glycerol growth conditions. One consideration of the computational simulation is that butanoate is a product of a secondary pathway and S. marcescens may not overproduce genes associated with this pathway during normal exponential growth. Thus, it is not necessarily surprising that growth simulations of a wild-type strain with no genetic manipulation may not show flux through the butanoate metabolism pathways. The proposed n-butanol synthetic pathway and associated gene expression levels in S. marcescens Db11 are shown in Fig. 5b. The annotated pathway for n-butanol production is partially complete but is missing the enzyme required for the last reaction (EC 126.96.36.199 alcohol dehydrogenase). A bifunctional alcohol dehydrogenase (AdhE2) from Clostridium acetobutylicum is reported with highest activity of alcohol production  and can be an appropriate enzyme candidate for heterologous expression in S. marcescens. Addition of an alcohol dehydrogenase would theoretically enable production of n-butanol, but it is also likely that additional pathway optimization would be necessary to improve flux through the necessary reactions as expression of the remaining endogenous butanoate genes is relatively low (e.g. SMDB11_4664 acetoacetate-CoA transferase).
The third chemical target of interest is N-acetylneuraminic acid which is not natively produced by S. marcescens Db11. Compared to glucose-based N-acetylneuraminic acid synthesis (demonstrated in E. coli ), use of N-acetylglucsaomine/chitin as a carbon source has the following advantages: 1) it requires only two steps for N-acetylneuraminic acid synthesis from N-acetylglucosamine compared to glucose and 2) no additional nitrogen source is needed [15, 41]. Growth of S. marcescens on N-acetylglucosamine showed that the enzyme (SMDB11_0473) associated with uptake of N-acetylglucosamine is highly active (RPKM 967.20). After transport, SMDB11_0473 (N-acetylglucosamine transferase) is a reversible N-acetylglucosamine transferase that can generate intracellular N-acetylglucosamine , shown in Fig. 5c. Intracellular N-acetylglucosamine can be converted to N-acetylneuraminic acid by additional two biochemical reactions (N-acetylmannosamine 2-epimerase and N-acetylneuraminic acid synthase) . Further overproduction of N-acetylneuraminic acid may require carefully balance the heterologous pathways since heterologous expression may generate metabolic burden .
In this study, we constructed a genome-scale metabolic model, iSR929, for the chitinolytic bacterium Serratia marcescens. The model contains 929 genes, 1185 reactions and 1164 metabolites. Initial experimental testing showed that the iSR929 model is able to reasonably predict quantitative phenotypes for S. marcescens. Integrating transcriptomic data for growth on different carbon sources (e.g. glucose, N-acetylglucosamine and glycerol) to the model provides insight into metabolic difference that are condition-dependent including different utilization of the TCA cycle for growth on N-acetylglucosamine. Three biochemical targets (2,3-butanediol, n-butanol, and N-acetylneuraminic acid) were proposed as potential metabolic engineering implementation targets and pathway activities were given. Overall, our results provide insight into the metabolic function of S. marcescens and reveal potential targets for engineering S. marcescens for biochemical production.
Serratia marcescens Db11 was purchased from the Caenorhabditis Genetics Center (Twin City, USA http://www.cbs.umn.edu/CGC) . The S. marcescens strain was grown in M9 medium or LB medium supplemented with various carbon sources. Carbon sources in the M9 medium were used at the following concentrations: glucose, 5 g/L; N-acetylglucosamine, 5 g/L; glycerol 5%. S. marcescens pre-cultures were grown in LB medium at 30 °C and 250 rpm overnight. 2.5% seed culture was inoculated at 50 mL M9 minimum medium with a 250 mL Erlenmeyer flask at 30 °C, initial pH 7.5, and 220 rpm for growth experiments. Stock cultures of S. marcescens were maintained at − 80 °C in 26% (v/v) glycerol.
Construction of the model iSR929
The core S. marcescens GSMM is a compile of metabolic reactions occurring in S. marcescens, compiled based on biochemical information from genome annotations  and experimental evidence. A draft model of initial biochemical reactions was constructed based on gene reading frames that encode enzymes with predicted functions in the genomic annotations available from IMG, UniProt, and KEGG [27, 30, 44]. Specifically, Enzyme Commission (EC) numbers of annotated S. marcescens genes were adapted to select reactions from a set of database reactions.
Flux balance analysis (FBA)
To this end, no Serratia biomass equation is available; therefore, we modified E. coli biomass equation by including components that are essential for Serratia growth [51,52,53]. Detailed biomass equation can be found as follows:
Constraint-based models can be used to analyze in silico gene essentiality through running FBA . Briefly, using the GPR relationships, the effect of an in silico gene knockout on reaction activities is assessed by assuming a gene is deleted or non-functional. If a gene is in silico deleted, that reaction is constrained to “zero” flux to simulate the effect of the gene deletion. Then, the model with the new constraint is analyzed by FBA, maximizing flux using the growth objective for batch growth on glucose. If the maximum flux of the biomass reaction is “zero”, the deleted gene is determined as an essential gene. On the contrary, if the maximum flux of the biomass reaction is over “zero”, the deleted gene is determined as a non-essential gene. The process is repeated for the next cycle to analyze another gene by resetting all the constraints back to default values. The essentiality analysis was conducted using M9 minimal medium with glucose as carbon source.
The total RNA of S. marcescens Db11 was extracted from cells harvested during mid-log phase. A QIAGEN (Venlo, Netherlands) RNeasy Mini kit along with the proper RNA protect reagent was used. Then, the mRNA was isolated, enriched, and reverse transcribed into cDNA. The sequencing of cDNA was conducted using paired-end reads using Illumina Hiseq 2500 (San Diego, CA). The detailed information of the RNASeq data can be found in a previous publication . In addition, the completed gene transcriptional levels of S. marcescens are listed in Additional file 1: Table S1.
Measurement of cell density
The measurement of S. marcescens culture density was generally quantified at OD600 using a Biomate3 UV/VIS spectrophotometer (Thermo Fisher Scientific, Waltham, MA).
Quantification of secreted metabolites
The concentration of N-acetylglucosamine, glucose, glycerol and all metabolic end-products (acetic acid, succinic acid, 2,3-butanediol, and ethanol) were analyzed by HPLC described in details in our previous publication .
Integration of RNASeq data with genome-scale reconstruction
The integration of RNASeq data with the iSR929 model was analyzed using a mixed integer linear programming algorithm approach , In brief, gene expression states were determined according to a gene expression value (g) vs. an arbitrary threshold value (γ), shown as follows:
Where g is the relative expression level as determined by RNAseq (e.g. RPKM) and γ is the threshold value that was manually set. The resulting gene states were then integrated into iSR929 using the gene-protein-reaction (GPR) relationships to generate lists of reactions predicted to be high or low fluxes .
Determination of intracellular citric acid and isocitric acid concentration
Intracellular levels of citric acid and isocitric acid were evaluated using a citrate assay kit (Sigma-Aldrich, St. Louis, MO) and an isocitric acid kit (BioVision, Milpitas, CA), respectively. In brief, S. marcescens cell growth was stopped at 2.5 h (OD600 about 0.4) and cells were collected by centrifuging at 4 °C, 15,000 rpm, and 10 min. Cell pellets were lysed using a commercial detergent BugBuster (Millipore Sigma, Temecula, CA). Citric acid and isocitric acid assays were conducted according to the manual’s instruction. 10 μL supernatant samples were added in a total volume of 100 μL reaction. Citric acid and isocitric acid concentrations were determined at 570 nm and 450 nm by a VERSAmax microplate reader (Molecular Devices, Sunnyvale, CA), respectively.
Genome-scale metabolic model
Peptidoglycan subunit of Serratia
Long M, Ong W, Reed J. Computational methods in metabolic engineering for strain design. Curr Opin Biotechnol. 2015;34:135–41.
Vaaje-Kolstad G, Houston D, Sorlie M, Eijsink V. The chitinolytic machinery of Serratia marcescens-- a model system for enzymatic degradation of recalcitrant polysaccharides. FEBS J. 2013;280(13):3028–49.
Monreal J, Reese E. The chitinase of Serratia marcescens. Can J Microbiol. 1969;15:689–96.
Yan Q, Hong E, Fong SS. Study of ChiR function in Serratia marcescens and its application for improving 2,3-butanediol from crystal chitin. Appl Microbiol Biotechnol. 2017;101:7567–78.
Kim K, Creagh A, Haynes C. Effective production of N-acetyl-β-D-glucosamine by Serratia marcescens using chitinaceous waste. Biotechnol Bioprocess Eng. 1998;3:71–7.
Mekasha S, Byman IR, Lynch C, Toupalova H, Andera L, Nas T, Vaaje-Kolstad G, Eijsink V. Development of enzyme cocktail for complete saccharification of chitin using mono-component enzymes from Serratia marcescens. Process Biochem. 2017;56:132–8.
Vaikuntapu PR, Rambabu S, Madhuprakash J, Podile AR. A new chitinase-D from a plant growth promoting Serratia marcescens GPS5 for enzymatic conversion of chitin. Bioresour Technol. 2016;220:200–7.
Hamilton JJ, Marlow VL, Owen RA, Costa Mde A, Guo M, Buchanan G, Chandra G, Trost M, Coulthurst SJ, Palmer T, et al. A holin and an endopeptidase are essential for chitinolytic protein secretion in Serratia marcescens. J Cell Biol. 2014;207(5):615–26.
Suzuki K, Sugawara N, Suzuki M, Uchiyama T, Katouno F, Nikaidou N, Watanabe T. Chitinases a, B, and C1 of Serratia marcescens 2170 produced by recombinant Escherichia coli: enzymatic properties and synergism on chitin degradation. Biosci Biotechnol Biochem. 2002;66(5):1075–83.
Brurberg M, Eijsink V, Haandrikman A, Venema G, Nes I. Chitinase B from Serratia marcescens BJL200 is exported to the periplasm without processing. Microbiology. 1995;141:123–31.
Tuveng TR, Hagen LH, Mekasha S, Frank J, Arntzen MO, Vaaje-Kolstad G, Eijsink VGH. Genomic, proteomic and biochemical analysis of the chitinolytic machinery of Serratia marcescens BJL200. Biochim Biophys Acta. 1865;2017:414–21.
Suzuki K, Uchiyama T, Suzuki M, Nikaidou N, Regue M, Watanabe T. LysR-type transcriptional regulator ChiR is essential for production of all chitinases and a chitin-binding protein, CBP21, in Serratia marcescens 2170. Biosci Biotechnol Biochem. 2001;65(2):338–47.
Suzuki K, Shimizu M, Sasaki N, Ogawa C, Minami H, Sugimoto H, Watanabe T. Regulation of the chitin degradation and utilization system by the ChiX small RNA in Serratia marcescens 2170. Biosci Biotechnol Biochem. 2016;80:376–85.
Yan N, Chen X. Sustainability: don’t waste seafood waste. Nature. 2015;524:155–7.
Yan Q, Fong SS. Bacterial chitinase: nature and perspectives for sustainable bioproduction. Bioresour Bioprocess. 2015;2:31–9.
Dahiya N, Tewari R, Hoondai G. Biotechnological aspects of chitinolytic enzymes: a review. Appl Microbiol Biotechnol. 2006;71(6):773–82.
Iguchi A, Nagaya Y, Pradel E, Ooka T, Ogura Y, Katsura K, Kurokawa K, Oshima K, Hatton M, Parkhill J, et al. Genome evolution and plasticity of Serratia marcescens, an important multidrug-resistant nosocomial pathogen. Genome Biol Evo. 2014;6(8):2096–110.
Yan Q, Fong SS. Challenges and advances for genetic engineering of non-model bacteria and uses in consolidated bioprocessing. Front Microbiol. 2017;8:2060.
Roberts SB, Gowen CM, Brooks JP, Fong SS. Genome-scale metabolic analysis of Clostridium thermocellum for bioethanol production. BMC Syst Biol. 2010;4:31–47.
Burgard A, Pharkya P, Maranas C. Optknock: a bilevel programming framework for identifying gene knockout strategies for microbial strain optimization. Biotechnol Bioeng. 2003;84:647–57.
Ranganathan S, Suthers P, Maranas C. OptForce: an optimization procedure for identifying all genetic manipulations leading to targeted overproductions. PLoS Comput Biol. 2010;6:e1000744.
Yang L, Cluett W, Mahadevan R. EMILiO: a fast algorithm for genome-scale strain design. Metab Eng. 2011;13:272–81.
Palsson BO. Systems biology: constraint-based reconstruction and analysis. New York: Cambridge University Press; 2015.
Gowen CM, Fong SS. Genome-scale metabolic model integrated with RNAseq data to identify metabolic states of Clostridium thermocellum. Biotechnol J. 2010;5:759–67.
Vanee N, Brooks JP, Spicer V, Shamshuin D, Krokhin O, Wilkins JA, Deng Y, Fong SS. Proteomics-based metabolic modeling and characterization of the cellulolytic bacterium Thermobifida fusca. BMC Syst Biol. 2014;8:86–99.
Vanee N, Brooks JP, Fong SS. Metabolic profile of the cellulolytic industrial actinomycete Thermobifida fusca. Metabolites. 2017;7(4):57.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of gees and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Schellenberger J, Park J, Conrad T, Palsson B. BiGG: a biochemical genetic and genomic knowledgebase of large scale metabolic reconstructions. BMC Bioinf. 2010;11:213.
Thorleifsson S, Thiele I. rBioNet: a COBRA toolbox extension for reconstructing high-quality biochemical networks. Bioinformatics. 2011;27:2009–10.
Consortium TU. The universal protein resource (UniProt). Nucleic Acids Res. 2008;36:D190–5.
Chagoyen M, Pazos F. MBRole: enrichment analysis of metabolomic data. Bioinformatics. 2011;27:730–1.
Brooks JP, Burns W, Fong SS, Gowen CM, Roberts SB. Gap detection for genome-scale constriant-based models. Adv Bioinforma. 2012;2012:323472.
Vanee N, Roberts SB, Fong SS, Manque P, Buck G. A genome-scale metabolic model of Cryptosporidium hominis. Chem Biodivers. 2010;7:1026–39.
Shlomi T, Cabili MN, Herrgard MJ, Passon BO, Ruppin E. Network-based prediction of human tissue-specific metabolism. Nat Biotech. 2008;26:1003–10.
Murarka A, Dharmadi Y, Yazdani SS, Gonzalez R. Fermentation utilization of glycerol by Escherichia coli and its implications for the production of fuels and chemicals. Appl Environ Microbiol. 2008;74(4):1124–35.
Sadykov MR, Olson ME, Halouska S, Zhu Y, Fey PD, Powers R, Somerville GA. Tricarboxylic acid cycle-dependent regulation of Staphylococcus epidermidis polysaccharide intercellular adhesin synthesis. J Bacteriol. 2008;190:7621–32.
Bai F, Dai L, Fan J, Truong N, Rao B, Zhang L, Shen Y. Engineered Serratia marcescens for efficient (3R)-acetoin and (2R,3R)-2,3-butanediol production. J Ind Microbiol Biotechnol. 2015;42:779–86.
Zhang L, Guo Z, Chen J, Xu Q, Lin H, Hu K, Guan X, Shen Y. Mechanism of 2,3-butanediol stereoisomers formation in a newly isolated Serratia sp. T241. Sci Rep. 2016;6:19257.
Deng Y, Fong SS. Metabolic engineering of Thermobifida fusca for direct aerobic bioconversion of untreated lignocellulosic biomass to 1-propanol. Metab Eng. 2011;13:570–7.
Kang J, Gu P, Wang Y, Li Y, Yang F, Wang Q, Qi Q. Engineering of an N-acetylneuraminic acid synthetic pathway in Escherichia coli. Metab Eng. 2012;13:623–9.
Steiger MG, Mach-Aigner AR, Gorsche R, Rosenberg EE, Mihovilovic MD, Mach RL. Synthesis of an antiviral drug precursor from chitin using a saprophyte as a whole-cell catalyst. Microbial Cell Fact. 2011;10:102.
Wu G, Yan Q, Jones JA, Tang YJ, Fong SS, Koffas MA. Metabolic burden: cornerstones in synthetic biology and metabolic engineering applications. Trends Biotechnol. 2016;34:652–64.
Flyg C, Kenne K, Boman H. Insect pathogenic properties of Serratia marcescens: phage-resistant mutants with a decreased resistance to Cecropia immunity and a decreased virulence to Drosophila. J Gen Microbiol. 1980;120:173–81.
Markowitz V, Szeto E, Palaniappan K, Grechkin Y, Chu K, Chen I, Dubchak I, Anderson I, Lykidis A, Mavromatis K, et al. The integrated microbial genomes (IMG) system in 2007: data content and analysis tool extensions. Nucleic Acids Res. 2008;36:D528–33.
Oh Y, Palsson BO, Park S, Schilling C, Mahadevan R. Genome-scale reconstruction of metabolic network in Bacillus subtilis based on high-throughput phenotyping and gene essentiality data. J Biol Chem. 2007;282(39):28791–9.
Vaaje-Kolstad G, Westereng B, Horn SJ, Liu Z, Zhai H, Sorlie M, Eijsink VG. An oxidative enzyme boosting the enzymatic conversion of recalcitrant polysaccharides. Science. 2010;330(6001):219–22.
Watanabe T, Kimura K, Sumiya T, Nikaidou N, Suzuki K, Suzuki M, Taiyoji M, Ferrer S, Regue M. Genetic analysis of the chitinase system of Serratia marcescens 2170. J Bacteriol. 1997;179(22):7111–7.
Yan Q, Fong SS. Increasing carbon source uptake rates to improve chemical productivity in metabolic engineering. Curr Opin Biotechnol. 2018;53:254–63.
Yan Q, Fong SS. Design and modularized optimization of one-step production of N-acetylneuraminic acid from chitin in Serratia marcescens. Biotechnol Bioeng. 2018;115:2255–67.
Yan Q, Fong SS. Cloning and characterization of a chitinase from Thermobifida fusca reveals Tfu_0580 as a thermostable and acidic endochitinase. Biotechnol Rep. 2018;19:e00274.
Hertle R, Brutsche S, Groeger W, Hobbie S, Koch W, Konninger U, Braun V. Specific phosphatidylethanolamine depedence of Serratia marcescens cytotoxin activity. Mol Microbiol. 1997;26(5):853–65.
Poole K, Braun V. Influence of growth temperature and lipopolysaccharide on hemolytic activity of Serratia marcescens. J Bacteriol. 1988;170(11):5146–52.
Yichoy M, Nakayasu ES, Shpak M, Aguilar C, Aley SB, Almeida IC, Das S. Lipidomic analysis reveals that phosphatidylglycerol and phosphatidylethanolamine are newly generated phospholipids in an early-divergent protozoan, Giardia lamblia. Mol Biochem Parasitol. 2009;165(1):67–78.
No specific funding was received for this study.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its supplementary information files. The in-house scripts such as model file (Additional file 2), biomass objective function file (Additional file 3), exchange file (Additional file 4), and gpr file (Additional file 5) are also included in the supplementary data.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Average RPKM values of S. marcescens gene expression level at three different growth conditions. Table S2. Number of genes excluded/included after the lower/upper bound cutoff. Table S3. S.marcescens iSR929 reaction flux values after running Ruppin algorithm. Table S4. NAD+/NADH associated reactions under three carbon source conditions. Figure S1. Representative figure of gene expression level distribution under M9 glucose medium growth condition. Figure S2. Relative gene expression levels of reactions in the metabolic map of (A) Glycolysis (B) Pentose Phosphate Pathway of S. marcescens. (ZIP 898 kb)
Model file. (TXT 116 kb)
Biomass objective function file. (TXT 1 kb)
Exchange file. (TXT 2 kb)
GPR file. (TXT 199 kb)
About this article
Cite this article
Yan, Q., Robert, S., Brooks, J.P. et al. Metabolic characterization of the chitinolytic bacterium Serratia marcescens using a genome-scale metabolic model. BMC Bioinformatics 20, 227 (2019). https://doi.org/10.1186/s12859-019-2826-1
- Genome-scale metabolic model
- Serratia marcescens
- N-acetylneuraminic acid
- Flux balance analysis