大竹蟶消化酶活性的響應面法分析

張雨 陳愛華 吳楊平

摘要：旨在利用響應面法研究溫度、鹽度和pH值3個因素及其聯合效應對大竹蟶3種消化酶（蛋白酶、淀粉酶和脂肪酶）活性的影響。結果表明，大竹蟶消化酶活性隨著鹽度、溫度及pH值的變化呈先增大后減小的變化規律。另外，通過Box-Behnken模型預測大竹蟶中3種消化酶活性最佳時的養殖條件，得出蛋白酶活性在鹽度為27.89‰、溫度為27.29 ℃、pH值為7.75時最佳，為0.479 U/mg;脂肪酶活性在鹽度為28.19‰、溫度為26.80 ℃、pH值為7.99時最佳，為0.291 U/mg;淀粉酶活性在鹽度為27.85‰、溫度為25.97 ℃、pH值為8時最佳，為0.958 U/mg?？紤]到實際養殖操作的便利性，將養殖條件調整為鹽度28‰、溫度27 ℃、pH值8，在該條件下大竹蟶的消化能力較好。

關鍵詞：大竹蟶;消化酶;響應面;養殖條件

中圖分類號： S968.3 ?文獻標志碼： A ?文章編號：1002-1302（2020）11-0166-06

收稿日期：2019-07-09

基金項目：江蘇省海洋水產研究所青年科學基金（編號：SQ201704）;南通市科技創新計劃（編號：JC2018136）。

作者簡介：張 雨（1986—），男，廣西玉林人，碩士，工程師，主要研究方向為貝類遺傳育種與增養殖技術。E-mail：zhangyutan610@163.com。

通信作者：陳愛華，研究員，主要研究方向為貝類遺傳育種與增養殖技術。E-mail：chenah540540@aliyun.com。 ?消化酶（digestive enzyme）是指由動物消化腺和消化系統分泌的起消化作用的酶類，主要包括糖類水解酶（淀粉酶、纖維素酶）、蛋白酶和脂肪酶，是衡量生物體生理健康狀態的重要指標[1]。消化酶不但隨動物種類、健康狀況和生長階段的不同而有所差異，還受溫度、鹽度、pH值等的諸多環境因子的影響[2]。目前，有關貝類消化酶的研究僅限于單一環境因素，如有關褶紋冠蚌（Cristaria plicata）[3]、三角帆蚌（Hyriopsis cumingii）[3]、紫貽貝（Mytilus edulis）[4]、濱螺（Littorina littorea）[4]、櫛孔扇貝（Chlamys Farreri）[5]及鮑（Abalone）[6-7]等消化酶的研究，較少考慮環境因素的聯合作用。大竹蟶俗稱蟶子王，是廣泛分布于我沿海的貝類之一。江蘇呂四漁場內如東、啟東一線歷來具有相當的大竹蟶資源量，但是近年來圍填海、港口建設等因素對海洋生境造成了不同程度的影響，使得大竹蟶的生理健康狀態受到了嚴重挑戰。因此，開展環境因素聯合作用對大竹蟶消化酶的影響研究具有重要意義。

響應面法是一種解決多變量問題的數學統計分析方法，主要通過提供合理的試驗設計和數據統計分析、評價多種變量之間的交互作用而建立數學模型和三維立體響應曲面圖形，最終得出最優參數組合的多元回歸方程[8]。目前，響應面法在食品工藝優化[9-11]和藥物研發[12-13]等領域有較多應用，在水產動物研究方面運用得較少。本研究利用響應面法研究溫度、鹽度和pH值3個因素及其聯合效應對大竹蟶消化酶活性的影響，并利用多元回歸方程模擬出大竹蟶消化酶活性最高時的養殖條件，以期為大竹蟶養殖條件優化、人工苗種繁育和增殖放流提供科學依據。

1 材料與方法

1.1 試驗材料

本試驗所用大竹蟶取自江蘇省海洋水產研究所文蛤良種場，選取相同規格、健康無損的大竹蟶作為試驗用貝，在聚乙烯水槽中暫養7 d后，投喂105 inds/mL的球等鞭金藻（Isochrysis galbana）并每天換1次水，試驗開始前停喂1 d。

1.2 試驗方法

試驗采用Box-Behnken模型設計，根據預試驗結果，設計3因素3水平試驗，試驗設計及其結果見表1。本試驗所用鹽度、溫度及pH值均設3個梯度，其中鹽度分別設為22‰、28‰、34‰，溫度分別設為20、27、34 ℃，pH值分別設為6.5、8.0、9.5。海水鹽度通過海水、海水晶及淡水配制得到;溫度通過加熱板進行調節;pH值通過1 mol/L NaOH或1 mol/L HCl調至本試驗所需值。試驗在17個聚乙烯水槽中進行，每個水槽中放置20個大竹蟶，在試驗過程中應及時挑出死亡個體，脅迫3 d后任取3個大竹蟶進行解剖，取肝胰臟組織，用于消化酶[蛋白酶（PRS）、脂肪酶（LPS）和淀粉酶（AMS）]活性的測定。

蛋白酶、脂肪酶、淀粉酶活性及總蛋白含量均采用南京建成生物工程研究所的試劑盒測定并按說明書上的公式計算。

1.4 數據統計

使用Design-Expert 8.0軟件進行試驗設計和數據分析，以鹽度（A）、溫度（B）、pH值（C）為自變量，分別以蛋白酶活性、脂肪酶活性、淀粉酶活性為因變量建立模型方程，并對3種消化酶活性的最大值進行預測和條件優化?；貧w擬合大竹蟶消化酶模型曲線并用Origin 8.0軟件作相應的3D響應面圖。消化酶活性（y）的模型方程如下：

y=k+x1A+x2B+x3C+x4AB+x5AC+x6BC+x7A2+x8B2+x9C2+ε。

式中：k為截距;x1、x2、x3分別為鹽度與溫度、鹽度與pH值、溫度與pH值的一次效應;x4、x5、x6分別為鹽度與溫度、鹽度與pH值、溫度與pH值之間的交互作用;x7、x8、x9分別為鹽度、溫度、pH值之間的二次效應;ε為殘差，假定其服從均值為0的正態分布。

通過方差分析（analysis of variance，簡稱ANOVA）確定回歸方程模型及各試驗因素的顯著性，給出決定系數以考察模型的擬合優度，模型中各項效應采用最小二乘法進行估計。采用F統計量進行顯著性檢驗，當P<0.05時，表示差異顯著;當P<0.01時，表示差異極顯著。

2 結果與分析

2.1 蛋白酶活性模型分析及響應面分析

利用Design-Expert 8.0軟件對表2數據進行二次多元回歸擬合。如表3所示，F值為53.52，P<0.000 1，模型R2為0.980 7，修正系數為 0.967 3，表明該模型的擬合度較好，預測值與試驗值高度吻合，因此可用該模型來分析并預測大竹蟶蛋白酶活性。本試驗的精確度（Adeq precision）為19.917%，說明本試驗結果可信。另外，本試驗結果顯示，鹽度、溫度對大竹蟶蛋白酶活性的一次效應影響不顯著;pH值對大竹蟶蛋白酶活性的一次效應影響極顯著（P<0.01）;鹽度、溫度交互作用對大竹蟶蛋白酶活性的影響顯著（P<0.05），鹽度和pH值、溫度和pH值交互作用對大竹蟶蛋白酶活性的影響不顯著;鹽度、溫度和pH值的二次效應對大竹蟶蛋白酶活性的影響極顯著（P<0.01）。由上述回

根據上述結果，用Origin 8.0作出3D響應面圖。如圖1所示，當保持溫度、鹽度、pH值3個因子中的1個為中值時，隨著其他2個因子數值的升高，大竹蟶體內蛋白酶活性均呈現先上升后下降的趨勢;各響應曲面圖均顯示，在本試驗設置的溫度、鹽度、pH值范圍內存在蛋白酶活性的最大值。根據上述模型結果優化出蛋白酶活性最大時的養殖條件如下：鹽度為27.89‰，溫度為27.29 ℃，pH值為775，此時的PRS活性為0.479 U/mg。

2.2 脂肪酶活型模型的分析及響應面分析

用Design-Expert 8.0軟件對表2中的脂肪酶活性進行二次多元回歸擬合。如表4所示，F值為5.43，P<0.05，模型的R2為0.874 7，修正系數為0713 5，表明該模型的擬合度較好，預測值與試驗值高度吻合，因此可以用該模型來分析并預測大竹蟶的脂肪酶活性。本試驗的精確度為6.319%，說明本試驗結果較為可信。另外，本試驗結果顯示，鹽度、溫度和pH值的一次效應及交互效應對大竹蟶脂肪酶活性的影響不顯著;溫度、鹽度的二次效

應對大竹蟶脂肪酶活性的影響顯著（P<0.05），pH值的二次效應對大竹蟶脂肪酶活性的影響極顯著（P<0.01）。由上述回歸結果擬合出的鹽度（A）、溫度（B）和pH值（C）與脂肪酶活性的編碼值回歸方程如下：

根據上述結果，用Origin 8.0作出3D響應面圖。如圖2所示，當保持溫度、鹽度、pH值3個因子中的1個為中值時，隨著其他2個因子數值的升高，大竹蟶體內的脂肪酶活性均呈現先上升后下降的

趨勢;各響應曲面圖均顯示，在本試驗設置的溫度、鹽度、pH值范圍內存在脂肪酶活性的最大值。根據上述模型結果優化出脂肪酶活性最大時的養殖條件如下：鹽度為28.19‰，溫度為26.80 ℃，pH值為7.99，此時的LPS活性為0.291 U/mg。

2.3 淀粉酶活性模型分析及響應面分析

用Design-Expert 8.0軟件對表2中的淀粉酶活性進行二次多元回歸擬合。如表5所示，F值為20.41，P<0.01，模型的R2為0.963 3，修正系數為0.916 1，表明該模型的擬合度較好，預測值與試驗值高度吻合，因此可以用該模型來分析和預測大竹蟶淀粉酶活性。本試驗的精確度為12002%，說明本試驗結果可信。另外，本試驗結果顯示，鹽度和pH值對大竹蟶淀粉酶活性的一次效應影響不顯著，溫度對大竹蟶淀粉酶活性的一次效應影響顯著（P<0.05）;鹽度、溫度和pH值交互作用對大竹蟶淀粉酶活性的影響不顯著;鹽度、溫度和pH值的二次效應對大竹蟶淀粉酶活性的影響極顯著（P<001）。由上述回歸結果擬合的鹽度（A）、溫度（B）和pH值（C）與淀粉酶活性的編碼值回歸方程如下：

根據上述結果，用Origin 8.0作出3D響應面圖。由圖3可見，當保持溫度、鹽度、pH值3個因子中的1個為中值時，隨著其他2個因子數值的升高，大竹蟶體內的淀粉酶活性均呈現先上升后下降的趨勢，各響應曲面圖均顯示，在本試驗設置的溫度、鹽度、pH值范圍內存在淀粉酶活性的最大值。根據上述模型結果優化出淀粉酶活性最大時的養殖條件如下：鹽度為2785‰，溫度為25.97 ℃，pH值為8，此時的AMS活性為0.958 U/mg。

3 討論

鹽度、溫度及pH值均是影響濾食性貝類消化酶活性的關鍵因素[1]。貝類是變滲透壓動物，當海水鹽度過高或過低時，貝類機體會主動排除水分或鹽分以適應環境[14]。在調節滲透壓的過程中，貝類會消耗大量能量，消化能力相應減小。本試驗結果顯示，大竹蟶消化酶活性在鹽度為28‰左右時較大。鹽度過高或過低，消化酶活性均有所降低，符合上述規律。酶促反應動力學認為，在一定溫度范圍內，隨著溫度的升高，酶促反應的速度也加快，超出一定溫度后便開始下降。本試驗結果表明，當溫度約為27 ℃時，消化酶活性達到較大值，符合上述規律。pH值可以通過影響底物的解離狀況、酶分子活性部位上有關基團的解離、中間絡合物（ES）的解

離狀態來影響酶活性[15]。本試驗結果表明，大竹蟶在pH值約為8時，其消化酶活性達到較大值，可作為養殖的最佳參考條件。

不同貝類體內的消化酶活性與其食性有關，喜食大型藻的貝類一般其纖維素酶活性相對較高，如濱螺等腹足綱貝類;而食單細胞藻類的貝類淀粉酶、蛋白酶活性相對較高，如紫貽貝等雙殼貝類[4]。在本試驗中，投喂大竹蟶的餌料為球等鞭金藻，沒有纖維素構成的細胞壁，鑒于食性對貝類消化酶種類的影響，本試驗僅測定大竹蟶3種功能消化酶（蛋白酶、脂肪酶、淀粉酶）活性而未考察纖維素酶活性。本研究結果顯示， 大竹蟶體內消化酶活性排序淀粉酶>蛋白酶>脂肪酶， 這與馬氏珠母貝的研究結果[16]相似。

目前針對貝類消化酶的研究較少考察環境因素的聯合作用及曲面效應。本試驗構建的Box-Behnken模型擬合度很高，說明該模型能夠很好地解釋大竹蟶消化酶隨環境因素的變化。本試驗得出大竹蟶消化酶活性最佳時的養殖條件如下：鹽度為28‰左右，溫度為27 ℃左右，pH值約為8，該條件可用于實際生產。

參考文獻：

[1]楊蕙萍，童圣英，王子臣. 國內外關于水產動物消化酶研究的概況[J]. 大連水產學院學報，1998，13（3）：64-71.

[2]劉 敏. 不同溫度和鹽度對施氏獺蛤消化酶和免疫酶活力的影響[D]. 湛江：廣東海洋大學，2015.

[3]費志良，喬 慧，郝 忱，等. 3種淡水蚌消化酶活力的初步研究[J]. 淡水漁業，2006，36（4）：3-6.

[4]劉萬順，李 瀕. 海洋無脊椎動物消化酶的研究——Ⅰ、紫貽貝、日本蟳、濱螺消化酶的初步分析及其應用[J]. 山東海洋學院學報（自然科學版），1988（1）：54-62.

[5]張 碩，趙 艷. 櫛孔扇貝蛋白酶和淀粉酶活力的初步研究[J]. 大連水產學院學報，1997，12（1）：15-20.

[6]楊蕙萍，童圣英，王子臣. 皺紋盤鮑淀粉酶和褐藻酸酶的研究[J]. 水產學報，1998，22（4）：345-351.

[7]黃 勃，王桂林，李二超. 雜色鮑與九孔鮑消化酶活力的比較[J]. 水產學報，2003，27（2）：119-123.

[8]Xu D P，Zhou Y，Zheng J，et al. Optimization of ultrasound-assisted extraction of natural antioxidants from the flower of Jatropha integerrima by response surface methodology[J]. Molecules，2016，21（1）：1-12.

[9]涂行浩，張帥中，唐景華，等. 澳洲堅果油微乳體系的構建[J]. 熱帶作物學報，2019，40（2）：359-367.

[10]高丹丹，郭鵬輝，祁高展，等. 響應面法優化薄荷全草總黃酮的提取工藝[J]. 食品工業科技，2015，36（2）：299-303，322.

[11]高丹丹，常 通，曹郁生，等. 響應面法優化中性蛋白酶水解棉籽蛋白工藝[J]. 食品科學，2009，30（16）：153-157.

[12]Vuong Q V，Nguyen V T，Thanh D T，et al. Optimization of ultrasound-assisted extraction conditions for euphol from the medicinal plant，Euphorbia tirucalli，using response surface methodology[J]. Industrial Crops and Products，2015，63：197-202.

[13]Goudarzi M，Madadlou A，Mousavi M E，et al. Optimized preparation of ACE-inhibitory and antioxidarive whey protein hydrolysate using response surface method （RSM）[J]. Dairy Science and Technology，2012，92（6）：641-653.

[14]呂昊澤，劉 健，陳錦輝，等. 鹽度對長江口3種濾食性貝類濾水率、攝食率、同化率的影響[J]. 海洋科學，2016，40（8）：10-17.

[15]范德朋，潘魯青，肖國強，等. 溫度、pH對縊蟶（Sinonovacula constricta）消化酶活力的影響[J]. 海洋湖沼通報，2003（4）：69-73.

[16]朱曉聞. 溫度、鹽度、pH對馬氏珠母貝稚貝生理學特征的聯合效應研究[D]. 湛江：廣東海洋大學，2012.魯 珊，閻旭東，肖荷霞，等. 春播青貯玉米不同種植模式的增產效應[J]. 江蘇農業科學，2020，48（11）：172-175.