球形節(jié)桿菌A152產透明質酸酶的發(fā)酵優(yōu)化及動力學模型的建立

姜言暉，張京良，江曉路,3,4

?

球形節(jié)桿菌A152產透明質酸酶的發(fā)酵優(yōu)化及動力學模型的建立

姜言暉1，張京良2,3，江曉路1,3,4

1 中國海洋大學食品科學與工程學院，山東青島 266003 2 中國海洋大學醫(yī)藥學院，山東青島 266003 3 青島海洋生物醫(yī)藥研究院，山東青島 266071 4 青島明月海藻集團有限公司，海藻活性物質國家重點實驗室，山東青島 266400

姜言暉, 張京良, 江曉路.球形節(jié)桿菌A152產透明質酸酶的發(fā)酵優(yōu)化及動力學模型的建立. 生物工程學報, 2017, 33(11): 1883–1888.Jiang YH, Zhang JL, Jiang XL. Establishment of kinetics digital model for hyaluronate lyase productionbased on fermentation optimization of Arthrobacter globiformis A152. Chin J Biotech, 2017, 33(11): 1883–1888.

以制備高產量的透明質酸酶為出發(fā)點，利用5 L發(fā)酵罐對球形節(jié)桿菌A152發(fā)酵生產透明質酸酶的條件進行優(yōu)化，并對其發(fā)酵動力學模型進行研究。研究結果表明，轉速400 r/min、通氣量3.5 L/min時，生物量和透明質酸酶酶活力最優(yōu)，分別為5 g/L、16.4 U/mL；同時對發(fā)酵過程中菌體生長、產物生成及基質消耗的規(guī)律進行研究，應用Logistic方程、Luedeking-Piret方程和底物消耗的物料平衡方程建立了球形節(jié)桿菌A152發(fā)酵過程的動力學模型，并通過MATLAB軟件進行最優(yōu)參數(shù)估計和非線性擬合。模型的計算值與實驗值能較好地擬合，表明所建的模型能較好地反映發(fā)酵過程，可為發(fā)酵過程的在線控制和預測提供理論基礎。

透明質酸酶，發(fā)酵優(yōu)化，動力學模型，非線性擬合

透明質酸酶 (Hyaluronidase，HAase) 是能夠使透明質酸產生低分子化作用酶的總稱。根據(jù)作用機制不同，分為哺乳動物及動物毒液來源的內切β-N-乙酰氨基葡糖苷酶，水蛭、十二指腸蟲來源的內切-β-葡糖苷酸酶及細菌來源的透明質酸酶3類[1]。

透明質酸酶可改善藥物的組織滲透性，促進藥物擴散和吸收[2]，它降解生成的低分子透明質酸寡糖具有促進骨和血管的生成[3]、促進創(chuàng)傷愈合[4]、抗腫瘤[5]及抗炎[6]等功效，在醫(yī)藥方面具有良好的應用前景?，F(xiàn)階段，透明質酸酶多從動物睪丸組織提取[7]，由于原料有限、含量低等原因限制了它的開發(fā)應用，而微生物來源的透明質酸酶可避免這些缺陷。目前微生物來源的透明質酸酶的研究主要集中在菌株的選育方面[8]，對產酶菌株發(fā)酵動力學方面的研究未見報道。

發(fā)酵罐發(fā)酵是實驗室搖瓶發(fā)酵的放大，為工業(yè)化放大提供依據(jù)，有利于實現(xiàn)成果從實驗室水平向工業(yè)化水平的轉變。通過前期研究發(fā)現(xiàn)球形節(jié)桿菌A152發(fā)酵生產的透明質酸酶具有良好的溫度和pH穩(wěn)定性，在醫(yī)藥和生物化工方面具有良好應用價值。本文在搖瓶優(yōu)化培養(yǎng)基配方的基礎上，利用5 L發(fā)酵罐對其發(fā)酵條件進行優(yōu)化，通過對通氣量、轉速及補料的研究，確定最適的發(fā)酵條件及控制參數(shù)。在發(fā)酵條件優(yōu)化的基礎上，建立A152生產透明質酸酶的動力學模型，探索發(fā)酵過程變化規(guī)律，為發(fā)酵過程的預測和控制提供理論方程，為其在線控制和自動化提供理論依據(jù)。

1 材料與方法

1.1 菌種及培養(yǎng)基

球形節(jié)桿菌A152，由食品科學與工程學院應用微生物學實驗室提供。

種子培養(yǎng)基：K2HPO41 g/L、MgSO40.5 g/L、NaCl 0.5 g/L、CaCl20.01 g/L、牛肉膏5 g/L、蛋白胨1 g/L、HA 0.5 g/L，pH調至7.0，115 ℃滅菌30 min。

發(fā)酵培養(yǎng)基：K2HPO41.5 g/L、MgSO40.5 g/L、NaCl 0.5 g/L、CaCl20.01 g/L、牛肉膏5 g/L、蛋白胨1 g/L、HA 4 g/L，pH調至7.0，115 ℃滅菌30 min。

補料培養(yǎng)基：K2HPO41.5 g/L、MgSO40.5 g/L、NaCl 0.5 g/L、CaCl20.01 g/L、牛肉膏25 g/L、蛋白胨5 g/L、HA 20 g/L，pH調至7.0，115 ℃滅菌30 min。

1.2 儀器與設備

BLBIO-5GJ-4-H 5L四聯(lián)自動發(fā)酵罐，上海百倫生物科技有限公司；LDZX-75KB立式壓力蒸汽滅菌器，上海申安醫(yī)療器械廠；GL-20G-II高速冷凍離心機，上海安亭科學儀器廠；UV-6000PC紫外型可見分光光度計，上海元析儀器有限公司。

1.3 培養(yǎng)方法

取一環(huán)菌接種于裝有25 mL種子培養(yǎng)基的250 mL三角瓶中，160 r /min、28 ℃培養(yǎng)24 h，得一級種子培養(yǎng)液。將一級種子培養(yǎng)液以5%的接種量接入裝有25 mL種子培養(yǎng)基的250 mL三角瓶中，同種方式培養(yǎng)獲得二級種子液。將二級種子液以5%的接種量接種至裝有3.5 L發(fā)酵培養(yǎng)基的5 L發(fā)酵罐中，控制適宜的轉速和通氣量，于24 ℃條件下發(fā)酵培養(yǎng)。

1.3.1 發(fā)酵培養(yǎng)條件優(yōu)化

控制其他條件相同，分別考察轉速 (200、400、600、800 r /min)、通氣量 (1.75、3.50、4.75 L/min)、補料對發(fā)酵體系的溶氧量 (DO)、透明質酸酶的產量和生物量的影響，篩選適宜的轉速、通氣量及補料條件。

1.4 分析方法

1.4.1 生物量的測定

在發(fā)酵條件優(yōu)化實驗中，以660 nm的吸光值 (660) 表征生物量[9]；在動力學模型構建實驗中，以菌體干重表示生物量[10]。

菌體干重 ()：取10 mL發(fā)酵液于預先稱重的離心管 (0) 中，8 000 r/min離心10 min，棄上清，加入10 mL蒸餾水離心洗滌2次，棄上清，沉淀于80 ℃烘箱烘干至恒重 (n)。

=n–0

1.4.2 酶活力的測定

酶活力單位定義[11]：37 ℃條件下，每分鐘催化形成1 μmol不飽和雙糖消耗的酶量。酶活力的測定：取20 μL經6 000 r/min離心10 min的發(fā)酵上清液、3.98 mL 0.02 mol/L Tris-HCL (pH 7.5) 和4 mL 0.2%透明質酸 (0.02 mol/L Tris-HCL配制，pH 7.5) 加入15 mL的比色管中，37 ℃水浴反應20 min后，置于沸水浴中煮沸5 min，對照管以相同的條件加入滅活的發(fā)酵上清液，測定232 nm處光吸收值。

式中，指232 nm處的吸光值；?指毫摩爾消光系數(shù)，取5.5；d指比色皿的厚度，1 cm；指反應的時間 (min)；t指反應的總體積 (mL)；s指樣品的體積 (mL)。

1.4.3 殘?zhí)橇康臏y定

采用DNS法[12]測定發(fā)酵液中的殘?zhí)恰?/p>

1.4.4 溶氧量的測定

采用溶氧電極測定發(fā)酵液中的溶氧量。

1.5 動力學模型

1.5.1 菌體生長動力學模型的建立

利用Logistic方程[9]建立A152的菌體生長動力學模型，描述菌體生長速率的Logistic微分方程為

式中，d/d為菌體的生長速率 (g/(L·h))；為菌體濃度 (g/L)；為發(fā)酵時間 (h)；為比生長速率 (h–1)；max為最大比生長速率 (h–1)；max為最大菌體濃度 (g/L)。

初始條件=0時，得=0，積分式為

1.5.2 產物合成動力學模型的建立

采用Luedeking-Piret方程[9]建立A152的產物合成動力學模型，描述透明質酸酶合成速率的Luedeking-Piret微分方程為

式中，d/d為產物合成速率 (g/(L·h))；為產物的濃度 (g/L)；為與菌體生長速率有關聯(lián)的產物合成常數(shù) (g/(L·h))，為與菌體濃度關聯(lián)的產物合成濃度 (g/(L·h))?！?、=0時，產物合成與菌體生長相偶聯(lián)；≠0、≠0時，產物合成與菌體生長部分偶聯(lián)；=0、≠0時，產物合成與菌體生長非偶聯(lián)。

初始條件=0時，得=0，積分式為

1.5.3 底物消耗動力學模型的建立

在發(fā)酵過程中，底物消耗主要用于菌體的生長、產物合成及細胞代謝的維持。根據(jù)基質平衡原理，建立底物消耗的物料平衡微分方程[9]為

式中，–/為底物消耗速率(g/(L·h))，為底物濃度(g/L)，x/s為菌體對底物的得率(g/g)，p/s為產物對底物的得率(g/g)，m為菌體的維持因子(h–1)。

初始條件=0時，得=0，積分式為

1.5.4 數(shù)據(jù)處理與分析

采用EXCEL 2010及MATLAB R2012a進行實驗數(shù)據(jù)處理分析與模型擬合分析。

2 結果分析

2.1 發(fā)酵優(yōu)化結果

2.1.1 轉速對透明質酸酶發(fā)酵的影響

微生物發(fā)酵過程培養(yǎng)基組分的均勻和溶氧都與攪拌相關。由圖1可知，200 r/min時，溶氧始終保持極低的水平，溶氧不足導致菌體生長緩慢，透明質酸酶的合成延滯，發(fā)酵周期延長。而在400 r/min時， 18 h達到產酶最高點，為8.68?U/mL，轉速600 r/min與400 r/min的規(guī)律一致，在18?h達到產酶最高點，為8.66 U/mL，二者酶產量相差不大，而在800 r/min時，在14 h即可達到產酶最高點，酶產量較400 r/min和600 r/min低，為8.02 U/mL。綜上，選擇400 r/min作為發(fā)酵生產透明質酸酶的最適轉速。

圖1 轉速對Arthrobacter globiformis A152產透明質酸酶發(fā)酵過程的影響(A：DO、OD660；B：酶活力)

2.1.2 通氣量對透明質酸酶發(fā)酵的影響

通氣量可以改變培養(yǎng)基中菌體細胞之間的氧傳遞系數(shù)而間接影響菌體對培養(yǎng)基組分的利用效率，適當?shù)耐饬坑欣诰w細胞內外養(yǎng)分平衡、菌體濃度的積累以及目的產物的釋放[13]。由圖2可以看出，通氣量越大，溶氧低水平期越短，通氣量1.75 L/min時，生物量總體偏低，產酶最高點也略微延遲；通氣量為3.50 L/min時，酶產量在18?h達到最大，為8.16 U/mL，通氣量為4.75 L/min時，酶產量相對較低，為7.98 U/mL。通氣量低影響菌體生長，生物量和酶產量偏低；通氣量過高，培養(yǎng)基中的氧氣處于過飽和狀態(tài)，菌體代謝異常，酶產量略有下降。與通氣量為1.75 L/min相比，通氣量3.50 L/min和4.75 L/min生物量和酶活力均有提高，說明氧氣供應充足，可以充分滿足菌體生長，但通氣過量，會抑制菌體產酶。從節(jié)省能源和酶產量兩方面考慮，選擇3.50 L/min作為發(fā)酵生產透明質酸酶的最適通氣量。

2.1.3 補料對透明質酸酶發(fā)酵的影響

根據(jù)菌體的生長和初始培養(yǎng)基的特點，在發(fā)酵的某些階段進行補加培養(yǎng)基，使菌體或其代謝產物的生產時間延長，使產物生成增加[14]。由圖3可知，發(fā)酵前期，菌體生長迅速，底物消耗增加，殘?zhí)橇繙p少，在12?h最低，限制菌體的生長。因此，選擇在12 h補料600 mL。與未補料相比，補料后殘?zhí)橇吭黾?，為菌體生長提供充足的營養(yǎng)，生物量增加，但酶活力卻有所降低，不利于產酶，后續(xù)發(fā)酵選擇不補料。

圖2 通氣量對Arthrobacter globiformis A152產透明質酸酶發(fā)酵過程的影響(A：DO、OD660；B：酶活力)

2.2 Arthrobacter globiformis A152生產透明質酸酶的發(fā)酵過程曲線及分析

A152產透明質酸酶的發(fā)酵過程曲線如圖4所示，殘?zhí)橇肯染徛陆岛罂焖傧陆底詈筅呌谄椒€(wěn)，在發(fā)酵4 h內緩慢下降，4–12 h迅速下降，12 h之后維持在0.5 mg/mL以下，保持穩(wěn)定。生物量與其變化趨勢相反，12 h后穩(wěn)定在5 g/L，生物量與殘?zhí)橇康淖兓竭M行。產酶量的變化略有不同，0–6 h緩慢增加，6–18 h迅速升高，在18 h后維持在16.4 U/mL左右，達到最高。以上數(shù)據(jù)表明，發(fā)酵過程中菌體生長消耗底物，使基質中的殘?zhí)橇砍掷m(xù)降低，生物量逐漸增加，透明質酸酶酶活力逐漸升高，當殘?zhí)橇吭?.5 mg/mL以下時，不足以維持菌體的生長，此時生物量達到最大，pH降至最低，產酶量繼續(xù)增加，在4?h后達到產酶最高點。pH在整個發(fā)酵的過程中呈先降低后升高的趨勢，表明菌體在合成透明質酸酶的過程中產酸，在12 h生物量達到最大，此時pH和殘?zhí)橇烤陆抵磷畹忘c，12 h后pH逐漸升高然后趨于穩(wěn)定，pH的變化可能是由于殘?zhí)橇坎蛔阋鸬?。酶產量在pH降至最低點后6?h達到最高點，可通過pH的變化對整個發(fā)酵終點進行控制，以獲得高產量的透明質酸酶。

圖3 補料對Arthrobacter globiformis A152生產透明質酸酶發(fā)酵過程的影響

圖4 Arthrobacter globiformis A152生產透明質酸酶的發(fā)酵曲線

2.3 模型參數(shù)求解及數(shù)據(jù)擬合

采用MATLAB軟件對球形節(jié)桿菌生產透明質酸酶的發(fā)酵動力學方程 (2)、(4)、(6) 進行非線性擬合，得到的動力學參數(shù)如表1所示。

將模型參數(shù)值代入 (2)、(4)、(6) 式，得到菌體生長、產物合成和底物消耗動力學方程分別為：

2.4 擬合曲線分析

A152發(fā)酵生產透明質酸酶的動力學曲線如圖5所示，從圖中可以看出，擬合值與實驗值較為接近，曲線擬合度2分別為0.983 8、0.973 6和0.921 2，表明模型的擬合效果良好，能夠較好地反映實際的發(fā)酵過程。

表1 發(fā)酵動力學模型參數(shù)值

3 討論

發(fā)酵工程的理論和方法已廣泛應用于工業(yè)微生物產品的生產及工藝優(yōu)化，使現(xiàn)代工業(yè)微生物發(fā)酵取得了高速發(fā)展[15]。本文利用5?L發(fā)酵罐生產透明質酸酶，通過優(yōu)化其培養(yǎng)條件，得到穩(wěn)定的發(fā)酵工藝，在轉速為400 r/min，通氣量為3.5 L/min的條件下發(fā)酵培養(yǎng)18 h，其產酶高達16.4 U/mL，與搖瓶 (發(fā)酵周期為24 h，酶產量為12.8 U/mL) 相比，發(fā)酵周期縮短了25%，酶產量提高了28%。且通過濁度法測得透明質酸酶酶活力達1.36×104U/mL，遠高于報道的輕型鏈球菌MTCC 2695[16]和噬尼古丁節(jié)桿菌[17]產的透明質酸酶的酶活力，且本文的發(fā)酵周期縮短，可降低生產成本，提高生產效率。此外，本研究的菌株易于培養(yǎng)，營養(yǎng)要求簡單，可規(guī)?；瘧糜诠I(yè)生產。

發(fā)酵動力對發(fā)酵工藝的分析、優(yōu)化、放大和控制具有重要的指導意義。王瑩等[18]對羊肚菌液體深層發(fā)酵動力學進行研究，對菌體生長、多糖產生和基質消耗模型進行了參數(shù)估計。本文對A152發(fā)酵產透明質酸酶的動力學進行了研究，通過建立發(fā)酵動力學模型，可以實時監(jiān)控發(fā)酵過程中菌體生長、產物合成、底物消耗的變化規(guī)律。模型擬合結果表明，三條曲線的擬合度分別為0.983 8、0.973 6和0.921 2，擬合值與實驗值相近，模型能較好地反映發(fā)酵過程規(guī)律，對A152發(fā)酵過程的在線控制及工業(yè)化生產具有重要指導意義。考慮到實驗過程中不能進行全程自動化控制，人工抽樣檢測可能會對實驗結果帶來一定影響，需要進一步的實驗對其進行深入研究。

[1] Meyer K, Rapport MM. Hyaluronidases. Adv Enzymol, 1952, 13: 199–236.

[2] Su K, Ji AG. Research of Hyaluronidases and Future Development. Biotechnol Bull, 2014, (3): 15–21 (in Chinese).蘇康, 吉愛國. 透明質酸酶的研究進展. 生物技術通報, 2014, (3): 15–21.

[3] Wang YF, Wang FS, Guo XP. Preparation of low molecular weight hyaluronic acid and its angiogenesis effects. Chin J Biochem Pharm, 2007, 28(2): 107–109 (in Chinese).王彥厚, 王鳳山, 郭學平. 低相對分子質量透明質酸的制備及其促血管生成作用. 中國生化藥物雜志, 2007, 28(2): 107–109.

[4] Hu ZL. Effects of hyaluronic acid in healing of periodontal wound. China Pharm, 2005, 14(9): 23–24 (in Chinese).胡澤利. 透明質酸在牙周組織創(chuàng)傷愈合過程中的作用. 中國藥業(yè), 2005, 14(9): 23–24.

[5] Alaniz L, Rizzo M, Malvicini M, et al. Low molecular weight hyaluronan inhibits colorectal carcinoma growth by decreasing tumor cell proliferation and stimulating immune response. Cancer Lett, 2009, 278(1): 9–16.

[6] Zhang HW, Wei WR, Cheng XF, et al. To investigate the effect of biological activity of hyaluronic acid brushing treatment of gingivitis. J Clin Med, 2015, 2(23): 4808–4809 (in Chinese).張華偉, 魏旺榮, 程錫芳, 等. 探討生物活性透明質酸刷牙治療牙齦炎的療效. 臨床醫(yī)藥文獻電子雜志, 2015, 2(23): 4808–4809.

[7] Yu J, Jin HN, Bai YJ, et al. A stabilizer and liquid preparation for hyaluronidase: CN, 104768535A. 2015-07-08 (in Chinese).禹具, 金荷娜, 白永俊, 等. 用于透明質酸酶的穩(wěn)定劑和包括透明質酸酶的液體制劑: 中國, 104768535A. 2015-07-08.

[8] Xue W. Strain Screening and Optimization of Fermentation Medium for Higher Production of Hyaluronidase. Jinan: Shandong University, 2015 (in Chinese).薛蔚. 產透明質酸酶菌株的篩選及發(fā)酵條件的優(yōu)化.濟南: 山東大學, 2015.

[9] Yu XM, Zhang JL, Sun YC, et al. Establishment of digital model for levan production based on fermentation optimization and kinetics. Food Fermen Ind, 2016, 42(12): 20–25 (in Chinese).于曉萌, 張京良, 孫永超, 等. 基于發(fā)酵優(yōu)化和動力學建立Levan果聚糖生產的數(shù)字化模型. 食品與發(fā)酵工業(yè), 2016, 42(12): 20–25.

[10] Zhang Y, Feng XH, Li S, et al. Fed-batch fermentation kinetics of ε-polylysine by. MY 5-36. Chin J Bioprocess Eng, 2012, 10(3): 23–27 (in Chinese).張揚, 馮小海, 李莎, 等.. MY 5-36產ε-聚賴氨酸分批補料發(fā)酵動力學. 生物加工過程, 2012, 10(3): 23–27.

[11] Allen AG, Lindsay H, Seilly D, et al. Identification and characterisation of hyaluronate lyase from. Microb Pathog, 2004, 36(6): 327–335.

[12] Yan FH, Zhou JM, Wu RC. Colorimetric determination of reducing sugar by methanoic acid pretreatment bagasse with 3,5-dinitrosalicylic acid. Food Res and Dev, 2015, 36(2): 126–128 (in Chinese).彥繁鶴, 周金梅, 吳如春. DNS法測定甘蔗渣中還原糖含量. 食品研究與開發(fā), 2015, 36(2): 126–128.

[13] Roukas T, Liakopoulou-Kyriakides M. Production of pullulan from beet molasses byin a stirred tank fermentor. J Food Eng, 1999, 40(1/2): 89–94.

[14] Hu B, Liu H, Wu ZL, et al. Study on D ribose production by fed-batch fermentation. Chem React Eng and Technol, 2006, 22(2): 151–155 (in Chinese).胡濱, 劉虎, 吳兆亮, 等. 補料發(fā)酵生產D-核糖工藝的研究. 化學反應工程與工藝, 2006, 22(2): 151–155.

[15] Xia Y, Wang B, Wang GX, et al. Establishment of batch fermentation kinetics model of polysaccharide ofgroup C. Chin J Biol, 2014, 27(11): 1473–1476 (in Chinese).夏永, 王斌, 王公孝, 等. C群流行性腦膜炎球菌多糖分批發(fā)酵動力學模型的建立. 中國生物制品學雜志, 2014, 27(11): 1473–1476.

[16] Mahesh N, Balakumar S, Parkavi R, et al. Optimization and production of hyaluronidase by streptococcus mitis MTCC 2695. Biomolecul, 2012, 1: 101.

[17] Su K. Study on optimal cultivation condition and characterization of hyaluronidase from[D]. Jinan: Shandong University, 2014 (in Chinese).蘇康. 噬尼古丁節(jié)桿菌產透明質酸酶條件優(yōu)化及酶學性質研究[D]. 濟南: 山東大學, 2014.

[18] Wang Y, Piao MZ, Sun YH. Genetic algorithm for fermentation kinetics of submerged fermentation by. Chin J Biotech, 2008, 24(8): 1454–1457 (in Chinese).王瑩, 樸美子, 孫永海. 基于遺傳算法的羊肚菌液體發(fā)酵動力學模型的建立. 生物工程學報, 2008, 24(8): 1454–1457.

(本文責編 郝麗芳)

Establishment of kinetics digital model for hyaluronate lyase productionbased on fermentation optimization ofA152

Yanhui Jiang1, Jingliang Zhang2,3, and Xiaolu Jiang1,3,4

1 College of Food Science and Engineering, Ocean University of China, Qingdao 266003, Shandong, China2 College of Medicine and Pharmacy, Ocean University of China, Qingdao 266003, Shandong, China3 Marine Biomedical Research Institute of Qingdao, Qingdao 266071, Shandong, China4 State Key Laboratory of Bioactive Seaweed Substances, Qingdao Brightmoon Seaweed Group Co. Ltd, Qingdao 266400, Shandong, China

In order to produce hyaluronate lyase of high yield, we optimized the fermentationA152 in quadruple fermentation of 5 L, and studied the kinetics of fermentation. Both the highest biomass and enzyme activity could be achieved when the rotation speed was 400 r/min and the ventilation volume was 3.5 L/min. In addition, digital models of cell growth, product synthesis and substrate consumption were built by equation of logistic, luedeking-piret, product synthesis and substrate consumption. Nonlinear fitting and estimation of optimal parameters were obtained by MATLAB. The model correlated well between prediction and experimental data, and reflected the change rules of cell growth, hyaluronidase synthesis and substrate consumption during the process of producing hyaluronate lyase. The establishment of fermentation kinetics digital models could provide basis for controlling and prediction of the production process.

hyaluronate lyase, fermentation optimization, kinetics model, nonlinear fitting

December 13, 2016;

June 8, 2017

Xiaolu Jiang. Tel: +86-532-82032290; E-mail: jiangxl@ouc.edu.cn

Supported by:Scientific and Technological Project in Shandong Province (No. 2015GSF115002).

山東省科技廳 (No. 2015GSF115002) 資助。

2017-09-15

網絡出版地址：http://kns.cnki.net/kcms/detail/11.1998.Q.20170915.1039.001.html