Isolation and identification of a psychrotolerant dimethyl phthalate-degrading bacterium from selected frozen soil of high-latitude areas in China and optimization of its fermentation conditions using response surface methodology

Abstract Large quantities of phthalic acid esters (PAEs) are released and migrate into the environment continuously and become ubiquitous in soil, surface water, groundwater and other areas, eventually leading to serious environmental pollution. Therefore, the methods of PAEs mineralization have become a hot research topic. In this study, a highly efficient dimethyl phthalate (DMP)-degrading bacterium capable of utilizing DMP as the sole source of carbon and energy was isolated from frozen soil collected from a refuse disposal station. The DMP-degrading bacterium was identified as Pseudomonas sp. QDF12 according to the phylogeny of 16S rDNA sequence. The degradation process well fitted a first-order kinetic model, with a half-life of 9.8 h and the highest degradation rate of 20.19 mg/L•h at 15.0 °C. The fermentation conditions that enhanced the fermentation of biomass were optimized via Plackett-Burman design and response surface methodology. Finally, the carbon-nitrogen ratio, CaCO3, temperature and shaking speed were identified as key factors. Verification experiments indicated that the maximal biomass was 1.85 × 109 cfu/mL under the optimized conditions, which was 134.18% and 66.67% higher than those obtained under the Luria broth medium and the basic fermentation conditions, respectively. The time and temperature of fermentation were reduced to 12 h and 12.0 °C respectively. In conclusion, the strain QDF12 is more suitable for application in low-temperature areas owing to its fast, efficient fermentation and low-temperature resistance. Because of the shorter fermentation period and lower fermentation temperature, the large-scale commercial application of strain QDF12 will reduce the production costs and accelerate the fermentation speed.


Introduction
Phthalic acid esters (PAEs) are a class of refractory synthetic organic plasticizer compounds which are widely utilized as softeners and plasticizers in various polymeric products, including plastics, coatings, paints, daily chemical products and pesticides [1]. Some PAEs are easily released in the process of using plastics because their interactions with plastic molecules are not a chemical bond [2]. Therefore, these compounds become global environmental organic pollutants found in soil, water, the atmosphere and even the human body [3]. PAEs can be divided into high-molecular-weight (HMW) and low-molecularweight (LMW) PAEs, and the latter include dimethyl phthalate (DMP), diethyl phthalate (DEP), dibutyl phthalate (DBP), diisobutyl phthalate (DiBP) and diethylhexyl phthalate (DEHP). LMW PAEs are used predominantly as solvents in adhesives, waxes, inks, cosmetics, insecticides and pharmaceuticals [4]. LMW PAEs entering the bodies through food chains act as anti-androgens in male mammals to significantly injure reproductive functions, such as damaging the reproductive system, causing testicular atrophy and sperm loss, and changing the ultrastructure of germ cells. Furthermore, they are moderately toxic to the embryonic development of female mammals [5][6][7]. LMW PAEs are classified as environmental pollutants of priority by the U.S. Environmental Protection Agency (EPA), European Union and the China National Environmental Monitoring Center [8].
In China, PAEs account for 90% of the plasticizers usage in polyvinyl chloride productions. The current consumption of PAEs in China is more than 2.2 million tons per year and is predicted to increase [9]. Highload PAEs have been found in 2.3 billion acres of cultivated soil in China [10]. The top three provinces with the highest concentrations of P 15PAEs are Fujian, Guangdong and Xinjiang, and the concentration of DEHP particularly exceeds those of all other organic pollutants. The highest soil concentration of PAEs in some parts of Hubei Province of China is 157,620 ng/g, which is 24 times that of the safety standard of the U.S. EPA (6500 ng/g) [11]. The soil concentrations of DMP, DBP and DEHP are all significantly different among provinces. The high LMW PAEs concentration seriously reduces the quality and yield of crops, indicating that LMW PAEs can seriously hinder the agricultural productions [12]. Moreover, LMW PAEs can inhibit the growth and reduce the diversity and activity of soil microorganisms, and weaken soil resistance against soil-borne diseases. This could lead to a reduction in soil fertility, crop quality and yield and eventually to serious soil degradation [13].
In order to reduce the high residues of LMW PAEs in the environment and improve soil fertility, researchers are looking for ways to mineralize PAEs. Because the hydrolysis and photolysis rates of PAEs are very slow, biodegradation is the main pathway for the decomposition of environmental PAEs [14]. As reported, PAEs with shorter ester hydrocarbon chains are more easily biodegraded and mineralized by microorganisms [12]. Therefore, searching for efficient PAE-degrading bacteria is considered the primary step toward the removal of LMW PAE residues from soil. The DMP-degrading bacteria isolated so far are most prokaryotic microbes, including Acinetobacter, Arthrobacter, Bacillus, Paracoccus, Brevibacterium, Comamonas, Chryseobacterium, Enterobacter, Gordonia, Lactobacillus, Novosphingobium, Pasteurella, Xanhtomonas, Pseudomonas, Bukholderia, Rhodococcus, Pandoraea and Sphingomonas [15][16][17][18]. The optimum degradation temperatures of these bacteria ranging from 20 to 40 C will not be applicable to high-latitude environments with low annual average temperature, since low-temperature significantly inhibits the degradation activity of DMP-degrading bacteria. However, few studies are available on psychrotolerant DMP-degrading bacteria. In this study, we isolated a psychrotolerant DMPdegrading bacterium from selected frozen soil collected from high-latitude areas in China. The discovery of this bacterium enriches the family of PAEs-degrading bacteria and provides a strain resource for DMP biodegradation in low-temperature environment.

Microorganism and medium conditions
Frozen soil polluted by domestic waste was collected at a depth of 0.5 m from a refuse disposal station (47 22'43" N, 124 0'42" E) in Tiefeng District, Qiqihar City. In each time, 5 g of frozen soil was loaded in a 250 mL Erlenmeyer glass flask containing 100 mL of minimum salt medium (MSM)-DMP solution. The flask was then incubated on an SHA-C rotary shaker (Guohua Electrical, China) at 15 C and 100 rpm in the dark for 8 d.
A complete MSM was prepared by adding 2 mL of fresh Solution B to 1000 mL of Solution A. Solution A consisted of 2.0 g/L NH 4 NO 3 , 3.0 g/L NaCl, 1.64 g/L K 2 HPO 4 3H 2 O and 0.39 g/L KH 2 PO 4 in 1000 mL of ddH 2 O; Solution B (trace elements) consisted of 0.246 g/L MgSO 4 7H 2 O, 0.166 g/L CoCl 2 6H 2 O, 0.074 g/ L MnSO 4 H 2 O, 0.288 g/L ZnSO 4 7H 2 O, 0.08 g/L CuSO 4 , 0.055 g/L CaCl 2 and 0.25 g/L FeSO 4 7H 2 O in 1000 mL of ddH 2 O. The MSM was initially adjusted to pH 7.6 with 1.0 mol/L NaOH or HCl. MSM-DMP was made by diffusing a DMP solution into MSM. A standard Luria broth (LB) medium (pH 7.6) was prepared by dissolving 3.0 g/L Beef extract, 5.0 g/L peptone and 5.0 g/L NaCl in 1000 mL of ddH 2 O. A domestication medium was prepared by mixing MSM-DMP and LB at appropriate proportions. The basic fermentation medium initially at pH 7.6 was composed of 5.0 g/L Beef extract, 2.0 g/L NH 4 NO 3 , 0.1 g/L CaCO 3 , 2.0 mL of Solution B, 100 lL of Tween-80, 7.0 g/L macro-elements (NaCl: K 2 HPO 4 3H 2 O: KH 2 PO 4 ¼30:16.4:3.6) and 1000 mL of ddH 2 O. The volume of the Erlenmeyer flask used in the basic fermentation experiment was 150 mL. The final DMP concentrations in the MSM-DMP, the domestication medium and the basic fermentation medium were all 1000 mg/L. DMP solution and Tween-80 were filter-sterilized through a 0.22-lm membrane before use. All other media and tools were sterilized in an autoclave at 121 C for 25 min.

Culture enrichment and isolation
In each experiment, 1 mL of a soil suspension was transferred into a 250-mL Erlenmeyer flask containing 100 mL of the initial domestication medium (MSM-DMP: LB ¼ 1:1, v/v). The flask was then incubated at 15 C and 100 rpm for 3 d in the dark. When the flask became turbid, 1 mL of the inoculum from the initial domestication medium was inoculated for up to 4 d into the second domestication medium, in which the proportion of MSM to LB increased to 2:1. In this way, the LB concentration in the second domestication medium decreased, whereas the MSM-DMP increased continuously. When the seventh domestication experiment (which lasted 9 d) was started, LB was no longer added to the domestication medium.
The seventh domestication culture (0.1 mL) was spread onto an MSM-DMP solid plate. After incubation at 15 C in the dark for 5-8 d, the colonies with different morphologies were transferred into the agar-free MSM-DMP and then incubated at the same conditions. When the liquid medium was turbid, 0.1 mL of the liquid culture was spread onto a new MSM-DMP solid plate. Three circulations of 'liquid medium-solid medium' alternative experiments were carried out according to the above processes. Finally, nine bacterial colonies (numbered QDA-QDI) were streaked onto slant test-tubes of MSM-DMP and then incubated at 15 C for 3 d in the dark.

Identification of DMP-degrading bacteria
The purified strain QDF12 was identified by its colour, morphology, size and Gram staining. The morphology of the colonies on agar plates was observed by the naked eye. Cell morphology of strain QDF12 was investigated by using optical microscopy (Olympus BX51, Japan), and its cell size was measured by a microscopic micrometer. The Gram reaction was examined by using the Gram stain.
Further identification was performed by 16S rDNA gene sequencing. The bacterial strain QDF was enriched overnight in the MSM-DMP, and then collected by centrifugation at 10,000 rpm (Sigma 1-14, Germany) for 5 min and were washed twice with sterile water. Cell-free supernatant was transferred into fresh sterile Eppendorf tubes and was used as template DNA for the polymerase chain reaction (PCR). Genomic DNA was extracted by a Takara Genomic DNA extraction kit. Then about 1500 bp length of 16S rDNA was amplified through PCR using the bacterial universal primers 27 F (5 0 -AGAGTTTGATCCTGGCTCAG-3 0 ) and 1492 R (5 0 -TACGGYTACCTTGTTACGACTT-3 0 ). The PCR amplification conditions for each PCR mixture (20 lL) were composed of 2.0 lL of 10 Â Taq buffer solution, 1.6 lL of deoxyribonucleoside triphosphate (dNTP)(2.5 mmol/L), 0.8 lL of each Primer (20 lmol/L), 0.5 lL of template DNA, 0.2 lL of Taq polymerase (5 U/ lL), and 14.1 lL of ddH 2 O. PCR was performed as follows [19]: preheating at 95 C for 5 min; then denaturation at 95 C for 30 s, annealing at 55 C for 30 s, extension at 72 C for 1.5 min for 24 cycles, and a final extension at 72 C for 10 min. Purified PCR products were directly sequenced using an ABI Prism 3730XL Sequencer by Shanghai Majorbio Bio-pharm Technology Co., Ltd. and compared with the referenced sequences from the National Center for Biotechnology Information (NCBI) using BLAST search. The closest references were retrieved from NCBI GenBank and aligned with CLUSTAL W with all parameters set at their default values. Phylogenetic trees were then generated by a neighbour-joining algorithm [20] and Maximum Composite Likelihood model of base substitution as implemented in MEGA7.0.14 [21]. The trees were validated using bootstrap analysis performed with 1000 replicates.

Degradation experiments
The MSM-DMP (100 mL) was inoculated by transferring a loopfull of bacteria from the slant culture into a 250-mL Erlenmeyer flask on a rotary shaker at 15 C and 100 rpm for 3 d. The medium was harvested on a centrifuge (Sigma 1-14, Germany) in a 50 mL Eppendorf centrifuge tube at 8000 rpm for 10 min. The cell pellets were washed twice with a 0.05 mol/L phosphate buffer solution (pH8.0) under sterile conditions, and resuspended in the same buffer to optical density at 600 nm (OD 600nm ) of 0.2. The bacterial suspension was prepared for the subsequent degradation and optimization experiments.
All degradation experiments were carried out in 500-mL glass flasks, each containing 1.0 mL of the bacterial suspension and 200 mL of MSM-DMP. The flasks were then incubated on the rotary shaker at 15 C and 100 rpm in the dark for 84 h. At the interval of 12 h, 2 mL of the solution in each flask was withdrawn and preserved at À20 C for residual DMP determination by high-performance liquid chromatography (HPLC) (Elite P230II, Dalian China) [22]. A blank control was set up in the same way except that no cells were added, which eliminated the possibility of DMP decomposition by chemical reactions, such as hydrolysis and oxidation. All the experiments were performed in triplicate. A high-performance DMP-degrading bacterium named QDF12 was obtained from the degradation experiments.

Effects of fermentation conditions on the growth of QDF12
Eight fermentation conditions were studied, including carbon-nitrogen ratio (C/N, beef extract/NH 4 NO 3 ), CaCO 3 concentration (g/L), culture temperature ( C), trace elements (mL), macro-elements (g/L), Erlenmeyer flask volume (mL), Tween-80 (lL/L) and shaking speed (rpm) ( Table 1). When one fermentation condition was tested, the other conditions were kept constant and at the same levels as those of the basic fermentation medium. The fermentation biomass was determined by enumerating viable cells with the 'drop count method' [23]. The fermentation solutions in decimal serial dilutions were prepared with sterile water. The 3-day-old MSM-DMP solid plates were each dropped with the fermentation solution (0.1 mL) and then incubated at 15 C for 48 h in the dark. The fermentation biomass was calculated as colony forming units (CFU) per milliliter.

Optimization of fermentation conditions for maximum biomass
The medium components and culture conditions for the fermentation of strain QDF12 were optimized in two steps: Plackett-Burman design (PBD) and response surface methodology (RSM). All experiments were repeated three times and the results were analyzed on Design Expert 8.0.5 and Minitab 16.0 [24].

Plackett-Burman design
PBD is a two-factorial design that identifies the critical physiochemical parameters required to elevate the fermentation biomass of strain QDF12 by screening n variables in n þ 1 experiments. All the variables were investigated at two widely-spaced intervals of Negative values (low level, À1) and Positive values (high level, þ1) [25]. The details of the design matrix and the results of variable screening were denoted in Table 2. The effects of individual parameters on the fermentation biomass, E(X i ), were determined as follows: where M þ i and M À i represent the response (the fermentation biomass) at which the parameter is at its higher or lower levels respectively; N is the total number of trials, which is 12.
Response surface methodology. Four significant factors were screened out, including x 1 (C/N), x 3 (CaCO 3 ), x 6 (temperature) and x 8 (shaking speed). Then the levels of these factors and the interaction effects were analyzed and optimized by the central composite design (CCD) [26]. Each variable was studied at five levels: Àa, À1, 0, þ1 and þ a. The full experimental plan consisted of 30 trials (Table 3). The fermentation biomass of strain QDF12 was taken as the dependent variable (Response Y). The second-order polynomial coefficients were calculated and analyzed on Design Expert software 8.0.5 (State-Ease Inc., USA) as follows [27]: where Y i is the predicted response (the fermentation biomass), i and j both range from 1 to the number of  Table 2. Experimental design matrix with response value of PBD for screening significant variables.

Runs
Variables a Biomass (Â10 7 cfu) is the ijth interaction coefficient, x i and x j are the levels of the independent variables (C/N, CaCO 3 , temperature and shaking speed) that influence the response variable Y. The significance of the model was determined via analysis of variance (ANOVA) performed by Fisher's test (F-value). The model adequacy was evaluated by the coefficient of determination R 2 , adjusted R 2 (which adjusts the number of explanatory terms in a model), and predicted R 2 (which in regression analysis indicates how well the model predicts the responses for new observations). For each variable, the quadratic models were represented as contour plots and response surface curves were generated on Design Expert.

Determination of DMP and HPLC conditions
For analysis of extracellular residual DMP in the solution, firstly the suspension of QDF12 (2 mL) was centrifuged at 10,000 rpm for 10 min. Then the supernatant was collected and filtered through a 0.22-lm filter for DMP determination by HPLC (Elite P230II, Dalian China) equipped with a diode array detector and a Hypersil ODS-C18 column (200 Â 4.6 Â 5lm). The mobile phase consisted of methanol and 0.01 mol/L phosphate solution (pH 7.0) (50:50, v/v) at 25 C and was added at a rate of 0.5 mL/min. The sample volume (20 lL) was injected into the column and detected at 254 nm. The retention time of DMP was identified according to reference standards [28]. The chemical concentration was quantified using external standards and calibration.

Data analysis
Statistical analysis of the data was done using SPSS 22.0 software and statistically significant differences were determined by applying Student's t-test. Values p < 0.05 and p < 0.001 were considered to indicate statistical significance. All data were expressed as mean values ± standard deviation (±SD).

DMP biodegradation kinetics
The degradation kinetics of nine bacterial strains were investigated at initial DMP concentration of 1000 mg/L. The first-order kinetics equation, which is frequently used to describe the biodegradation process, was used here to describe the DMP biodegradation by the DMPdegrading bacteria QDA-QDI [29]: where C is the residual concentration of DMP (mg/ L), k is the degradation rate constant, t is time (d) and A is a constant. The half-life of DMP degradation (t 1=2 , h) and the decreasing rate of DMP concentration (v,  mg/Lh) were calculated as follows: where C 0 is the initial concentration of DMP (mg/L). The experimental data of DMP biodegradation were fitted using Equation (3) (Figure 1 and Table 4). The degradation ability of the nine DMP-degrading bacteria can be divided into three levels of high, medium and low efficiency. Only strain QDF was classified as 'high efficiency', whereas strains QDC and QDE showed the lowest efficiency, and the other strains were of 'medium efficiency'. The degradation rates of the nine strains ranged from 9.44 to 20.19 mg/Lh, and the correlation coefficients R were all higher than 98%. The shortest half-life of degradation at 15 C was 9.8 h (strain QDF). In the MSM-DMP initially containing 1000 mg/L DMP, the degradation rate within 48 h by strain QDF reached 96.89%, and the DMP degradation rate was up to 20.19 mg/Lh, which both surpassed the results of all other strains. The DMP-degrading bacterium Cyanothece sp. PCC7822 screened out by Zhang et al. [30] consumed 40% of 50 mg/L DMP after 72 h and completely degraded it within 120 h, with the maximum degradation rate of only 0.67 mg/Lh. Another DMP-degrading bacterium, Rhodococcus sp. L4, screened and identified by Lu et al. [31] took 96 h to completely degrade 100 mg/L DMP. In other studies, DMP was degraded by mixed bacteria. For example, the multiple microorganisms collected from anaerobic sludge degraded 10 mg/L DMP completely after 6 d, and the half-life of degradation was up to 23.9 h [32]. In conclusion, we screened out a high-efficiency DMPdegrading bacterium that had a higher degradation rate compared with some DMP-degrading bacteria reported previously.

HPLC method validation
The analytical figures of merit were obtained using the implemented HPLC conditions used for DMP determination. The validation parameters studied were: sensitivity (linearity and linear range), limits of detection (LOD) and quantification (LOQ), and repeatability. A representative analytical curve equation was described as follows: where S peak is the peak area and C DMP is the concentration of DMP (mg/L). The correlation coefficient R of 0.9985 indicated the linear model was adequate to describe the proportionality between peak areas and DMP concentrations. LOD is the lowest detectable concentration of the analyte and provides a signal-to-noise ratio of 3, whereas LOQ is the lowest quantifiable concentration of the analyte and gives the signal-to-noise ratio of 10. In the present study, LOD and LOQ were 0.11 mg/L and 0.37 mg/L, respectively, from which it was inferred that the HPLC method was sensitive to determine the concentrations of DMP. The repeatability of the HPLC method was estimated by comparing the relative standard deviation (RSD%) obtained after three consecutive measurements of nine DMP solutions treated and analyzed by HPLC. The RSD% for the lowest levels was 0.9% (confidence interval of 95%), although it was below 2.8% for the other levels. These results suggested a very high precision of the HPLC method even in the lowest portion of the analytical curve. Table 5 summarizes the analytical parameters of merit of the analytical method described above. Based on these data, we concluded that the HPLC method can be used for the determination of DMP solutions under the experimental conditions.

Characterization and phylogenetic analysis of QDF12
After analysis of the degradation experiment results, we screened out a highly efficient DMP-degrading bacterium QDF12. The strain had cell size of 0.2-0.6 lm Â 1.1-2.1 lm and was gram-negative and short-rod-shaped. The strain formed colonies with a diameter of 0.5-1.5 mm after culturing for 24-36 h on the MSM-DMP agar medium; they were ivory, translucent and morphologically round. The colonies turned white after 84 h of incubation. Obvious clear water-like substances were found around the colonies, which may be an extracellular secretion of the bacterium containing a component capable of degrading extracellular DMP. The 16S rDNA gene was amplified via PCR using a universal bacterial primer set (27 F and 1492 R), which was obtained from strain QDF. The 16S rDNA partial gene sequences of strain QDF were aligned and compared with the 16S rDNA bacterial gene sequences in the GenBank on MEGA 7.0. Based on the phylogenetic analysis, the 16S rDNA sequence similarity search -BLAST showed the highest similarity with Pseudomonas (95%) and query cover of 100%. The phylogenetic tree was constructed ( Figure 2) using the neighbour-joining method with bootstrapping and rooted by Pseudomonas putida. It showed that the strains fell within a cluster of gammaproteobacteria Pseudomonas. The tree shows the relative relationship between strain QDF12 and related species. We compared the characteristics of strain QDF12 with various Pseudomonas species listed in the phylogenetic tree. However, we cannot clearly identify it at the species level. So far, the acceptable positional standard is that if the similarity of sequences between the tested strain and a reference strain is higher than 97-98%, they are regarded as the same genus [33]. Therefore, based on the results of the phylogenetic analysis and phenotypic analysis, the strain QDF12 was designated as Pseudomonas sp.QDF12. The partial nucleotide sequence of this strain was deposited in GenBank under the accession number MK491530. This strain was preserved in the China Center for Type Culture Collection (CCTCC) with the reference number of CCTCC NO. M2018272.

Effects of fermentation conditions on the growth of Pseudomonas sp. QDF12
The effects of fermentation conditions (C/N, CaCO 3 , culture temperature, trace elements, macro-elements, Erlenmeyer flask volume, Tween-80 and shaking speed) on fermentation biomass at different levels are shown in Figure 3. The main nutrients in the fermentation medium were of carbon (beef extract) and nitrogen (NH 4 NO 3 ) sources. The fermentation biomass depended on the carbon/nitrogen (C/N) ratio and was generally promoted by a higher C/N ratio [34,35]. The experiments showed that Pseudomonas sp. QDF12 had the highest fermentation biomass when the C/N ratio was 4/1 (Figure 3(a)), which was similar with another study [36]. The effects of trace elements on fermentation biomass were quite different from those of macro-elements (Figure 3(b,d)). Trace elements even at low concentrations (<3.5 mL/L) obviously inhibited the fermentation, whereas macro-elements only showed obvious inhibitory effect at high concentrations (>11.0 g/L). This indicates that low concentrations of trace elements will lead to a decline in the fermentation yield. In general, trace elements are involved in many metabolic pathways in cells, e.g. Mg 2þ , Zn 2þ and Cu 2þ , which are activators of many enzymes; Fe 2þ , Co 2þ , Mn 2þ and Mo 2þ , which are structural components of many special molecules;   Note: Error bars indicate averaged values ± standard deviation from three independent experiments. Statistically significant differences (p < 0.05) were determined by applying Student's t-test and denoted with different letters above each column.
Ca 2þ and K þ , which are cofactors of some enzymes and can keep cell stability and osmotic pressure. Sufficient trace elements should be maintained in order to keep the normal growth and metabolism of cells. Thus, the concentrations of trace elements should be strictly controlled.
CaCO 3 was added as a pH stabilizer to the fermentation medium. When abundant H þ existed in the fermentation medium, CO 3 2þ quickly combined with it to form HCO 3 þ or CO 2 and H 2 O. This kept the fermentation medium at neutral or weak alkaline conditions for a long time, thereby benefiting the growth of prokaryotic cells [37]. Moreover, the maximum biomass with the presence of 0.25 g/L CaCO 3 was 22.17% higher than that in the basic fermentation medium (p < 0.05) (Figure 3(c)), further indicating that the stable pH of the fermentation medium significantly promoted the fermentation. Tween-80 is a commonly-used nonionic surfactant that enhances the permeability of cell membranes and increases the cellular ability to absorb nutrients or secrete metabolites [38]. In this experiment, the addition of 50-100 lL/L Tween-80 enhanced the fermentation (Figure 3(e)). As reported similarly, the non-ionic surfactant below the concentration of 5 mL/L in the culture medium could significantly improve the permeability of cell membranes and increase the fermentation yield [39]. At present, most of the identified DMP-degrading bacteria are mesophilic bacteria, which grow at an optimum temperature between 20 and 40 C [16,[27][28][29][30]. This temperature range is not suitable for the growth of these bacteria in northern areas, especially in the environment where the annual average temperature is below 20 C. Thus, the lower temperatures in the studied areas were more favourable for the growth of Pseudomonas sp. QDF12. The highest fermentation biomass at 15.0 C was 1.14 Â 10 9 cfu/mL (Figure 3(f)), and the corresponding DMP degradation rate was up to 20.19 mg/ Lh ( Table 4), indicating that Pseudomonas sp. QDF12 was suitable for production and application in lowtemperature areas.
We innovatively chose the shaking bottle volume as an optimization factor for the first time. The gas/ liquid contact area per unit volume of a fermenter is decided by its fermentation capacity, which affects the volumetric oxygen transfer coefficient, thereby affects the oxygen transfer rate (OTR) and the fermentation efficiency [40]. The results showed that the fermentation biomass maximized when the flask volume, the gas/liquid contact area of the fermentation medium and the specific surface area of the fermentation liquid were 300 mL, 46.57 cm 2 and 0.1552 cm 2 /mL, respectively (Figure 3(g)). The fermentation biomass began to decrease when the flask volume was enlarged. Reportedly, the fermentation medium shakes more intensely in a larger-volume shaking bottle under a constant shaking speed (100 rpm), and cells cannot grow in a relatively unstable environment, resulting in slow growth. Especially in the process of producing a metabolite, a faster shaking speed will greatly reduce the fermentation yield [41].

Optimization of fermentation conditions
Identifying the significant factors using Plackett-Burman design The important fermentation factors with respect to their main effects were screened out by PBD. The effects of the screened factors on the fermentation biomass were studied through 12 runs of experiments ( Table 2). Statistical analysis of the responses was also performed ( Table 6). The obtained F-value of 33.92 indicates that the model is significant (p < 0.05). The coefficient of determination R 2 (0.9890) and adjustment coefficient R (0.9797) indicate that the data variability can be well explained by the models [42].
The p-value of ANOVA helped to identify the significance of each term at the designed confidence level. Any variable with confidence level above 95% was considered as significant. Accordingly, x 1 , x 3 , x 6 and x 8 were significant terms, and the other factors in the tested ranges did not cause significant variation in the fermentation process. The effects of the eight factors on the fermentation biomass in the tested ranges ranked as follows x 1 >x 3 >x 8 >x 6 >x 2 and x 5 >x 7 >x 4 . Based on the ANOVA results, a first-order polynomial equation representing the production of fermentation biomass of QDF12 against eight independent variables was derived: The positive or negative sign represents the direction of the effect that a certain factor has on the fermentation biomass. The coefficients of C/N (À17.58), CaCO 3 (þ11.75), culture temperature (À0.79) and shaking speed (À0.22) indicate that the high levels of CaCO 3 concentrations and the low levels of C/N, culture temperature and shaking speed are all beneficial for the fermentation biomass.
In conclusion, the fermentation conditions of C/N, CaCO 3 , culture temperature and shaking speed all significantly affected the fermentation biomass. The optimal levels of these four factors were further optimized by statistical optimization methods.

Optimization by RSM
The interactions among the significant factors and their optimal levels were studied via CCD [43]. When a significant factor was studied, the other factors were maintained constant, which yielded the maximal fermentation biomass (Figure 3). The results of the fermentation biomass are presented in Table 3. According to the results of multiple regression analysis, a second-order polynomial equation was derived as follows: Y ¼ 175:50 þ 17:87x 1 þ 6:79x 3 À5:13x 6 þ2:04x 8 À2:81x 1 x 3 À4:06x 1 x 6 À1:94x 1 x 8 À2:06x 3 x 6 À3:69x 3 x 8 þ 4:06x 6 x 8 À21:14x 1 2 À4:51x 3 2 À5:89x 6 2 À7:26x 8 2 (8) where Y is the predicted fermentation biomass (Â10 7 cfu); x 1 , x 3 , x 6 and x 8 are actual values of C/N, CaCO 3 , culture temperature and shaking speed, respectively. The model adequacy was tested using ANOVA through Fisher's test ( Table 7). The F-value of 58.97 indicates that the model is significant and the Lack of Fit (LOF) F-value of 2.56 is insignificant. A coefficient of determination (R 2 ) closer to 1 indicates a better correlation between the predicted and observed values, and a smaller coefficient of variation (C.V.%) implies a higher reliability of the experiment [35,44]. In this study, the R 2 of 0.9823, predicted R of 0.9538 and adjusted R of 0.9826 were consistent, and the predicted R and adjusted R were within 0.20 of each other. This indicates the consistency between the experimental and predicted data. The C.V.% of 3.69 indicates the experiments were highly reliable. Diagnostic plots were drawn to assess how well the model approximated the true response surface [45], and some selected plots are shown in Figure 4. The normal probability of residuals (Figure 4(a)) was normally distributed and the residuals were approximated along a straight line, which signified an acceptable model [46,47]. The studentized residuals plotted against the response of fermentation biomass were scattered in a desired random way (Figure 4(b)). The actual responses plotted against the predicted responses showed good agreement (Figure 4(c)). Diagnostic plots, including the Box-Cox plot ( Figure  4(d)), indicated that data transformation was unnecessary and therefore was not performed here. All the statistical results prove the second-order polynomial equation is highly accurate and generally applicable and adequate to describe the responses observed in the experiments. The p-values are of vital importance for understanding the patterns of mutual interactions among the experimental variables. The p-value of less than 0.001 indicates high significance [48]. The linear terms of the independent variables including x 1 , x 3 , and x 6 significantly affected the response Y ( Table 7). The quadratic terms of the four variables and the interactions between x 1 and x 6 , between x 3 and x 8 , and between x 6 and x 8 also significantly affected the fermentation biomass. Notably, the effects of x 1 , x 3 , x 6 , x 1 2 , x 3 2 , x 6 2 and x 8 2 (p < 0.001) were higher than other effects, and thus were the most significant variables influencing the fermentation biomass. The response surface contour plots are generally the graphical representation of the regression equation. Based on the contour plots and response surface curves ( Figure 5(a-f)), the fermentation biomass at different levels of the four factors can be predicted [49]. Each contour ellipse plot represents an infinitive number of combinations of two factors, with the other two factors maintained at their central levels [50]. The maximum predicted value is indicated by the response surface confined in the smallest ellipse in the contour diagram. The maximal fermentation biomass under the conditions of 4.70 of C/N and 0.28 g/L of CaCO 3 was 1.81 Â 10 9 CFU/mL, whereas the temperature and shaking speed were both held at the zero level ( Figure 5(a)). The highest fermentation biomass at 4.94 of C/N and 12.1 C was 1.81 Â 10 9 CFU/mL, whereas CaCO 3 and shaking speed were held at 0.25 g/L and 110 rpm respectively ( Figure  5(b)). With C/N at 4.84 and shaking speed at 112 rpm, the fermentation biomass maximized to 1.79 Â 10 9 CFU/mL with CaCO 3 and temperature held at 0.25 g/L and 15.0 C respectively ( Figure 5(c)). At the conditions of 0.29 g/L CaCO 3 and 12.2 C (Figure 5(d)), 0.29 g/L CaCO 3 and shaking speed of 109 rpm ( Figure 5(e)), and 12.9 C and shaking speed of 110 rpm ( Figure 5(f)), the maximum fermentation biomass was 1.80 Â 10 9 , 1.78 Â 10 9 , and 1.77 Â 10 9 CFU/mL respectively. According to Equation (8), it can be predicted that the fermentation biomass will reach a maximum of 1.84 Â 10 9 cfu/mL when the four factors are C/N at 4.71:1 (beef extract 8.249 g/L, NH 4 NO 3 1.751 g/L), 0.29 g/L CaCO 3 , 12.0 C and shaking speed at 99 rpm. The other non-significant fermentation conditions were 3.5 mL/L of solution B, 7.0 g/L of

Validation and application of the experimental model
An optimized fermentation medium was prepared by combining all the optimized fermentation factors at the required levels. The maximum fermentation biomass peaked to 1.85 Â 10 9 cfu/mL under the optimized fermentation conditions (Figure 6), which agrees well with the prediction of 1.84 Â 10 9 cfu/mL. This comparison validates that the statistical model is successful, accurate and reliable for predicting the fermentation biomass. The optimized fermentation medium can enhance the fermentation biomass and accelerate the growth compared with the basic fermentation medium. The maximum fermentation biomass was 66.67% higher compared with that obtained in the basic fermentation medium (1.11 Â 10 9 cfu/mL) and 134.18% higher compared with the LB medium (0.79 Â 10 9 cfu/mL). The fermentation time was shortened from 72 to 12 h, greatly improving the fermentation efficiency, and the optimum fermentation temperature dropped to 12.0 C. Overall, the DMPdegrading bacterium Pseudomonas sp. QDF12 is characterized by high efficiency, low energy consumption and high yield under the optimized conditions.

Conclusions
A psychrotolerant DMP-degrading bacterium QDF12 was isolated from selected frozen soil by gradient domestication and solid-liquid mixed separation and identified as Pseudomonas sp. The strain QDF12, which used DMP as the sole carbon and energy source, could efficiently degrade DMP at low temperature. The fermentation conditions of strain QDF12 were optimized through Plackett-Burman design and response surface methodology in order to enhance the fermentation yield. Compared with the basic fermentation medium and LB medium, the optimized fermentation conditions significantly increased the fermentation biomass, shortened the fermentation time and reduced the fermentation temperature. Therefore the application of the optimized fermentation conditions will reduce costs, increase speeds and improve benefits for large-scale productions of the DMP-degrading bacterium Pseudomonas sp. QDF12. Note: Data are averaged values from three independent experiments. Error bars indicate standard deviation (±SD). Statistically significant differences (p < 0.05) were determined by applying Student's t-test and were denoted with different letters above each column.