馮欣，登巴達(dá)吉，孔四新，2，李海奎，Shuya K.，仲格嘉*

1.中國藏學(xué)研究中心 藏醫(yī)藥研究所，北京 100094；2.世紀(jì)神農(nóng)(北京)國際生物技術(shù)有限公司，北京 100073；3.中國林業(yè)科學(xué)研究院 資源信息研究所，北京 100091；4.加利福尼亞大學(xué)，美國 圣克魯茲 CA 95064

波棱瓜Herpetospermumpedunculosum(Ser.) Baill.，又名塞季美朵、塞美塞古、色吉美朵，屬葫蘆科波棱瓜屬一年生攀援藤本植物，是藏藥材常用品種，在西藏主要分布于藏東南地區(qū)海拔2300～3500 m的山坡及林緣地帶，以種子入藥[1-3]；味苦，臨床用于清熱解毒、護(hù)肝，治黃疸性傳染性肝炎、消化不良等[1-4]。新近研究證實(shí)，其治療肝炎、肝纖維化的作用機(jī)制包括調(diào)節(jié)TGF-β信號(hào)通路[4]和改善能量代謝、雙向調(diào)節(jié)支鏈氨基酸等[5]。近年來，波棱瓜子市場需求劇增，受掠奪式采收等人為因素影響，野生資源日顯不足，亟需進(jìn)一步發(fā)展人工栽培。波棱瓜為雌雄異株，雌雄比嚴(yán)重失調(diào)；花果期較長，種子成熟期不一致，造成種子質(zhì)量不均一、發(fā)芽和出苗不整齊，從而影響其產(chǎn)量和質(zhì)量，嚴(yán)重制約波棱瓜產(chǎn)業(yè)發(fā)展。

本研究通過在西藏自治區(qū)林芝市設(shè)立定點(diǎn)試驗(yàn)小區(qū)，開展波棱瓜的人工栽培技術(shù)試驗(yàn)，在取得波棱瓜種子產(chǎn)量與相對應(yīng)技術(shù)措施對應(yīng)數(shù)據(jù)的基礎(chǔ)上，運(yùn)用數(shù)量化理論Ⅰ方法對定性因子量化分析，建立波棱瓜子產(chǎn)量(Y)與相關(guān)技術(shù)措施(催芽方法、播期、雌雄比、密度、施肥、疏花打頂)關(guān)系的回歸模型，探討多因子條件下篩選影響波棱瓜子產(chǎn)量的關(guān)鍵因子，減少波棱瓜人工栽培試驗(yàn)的工作量，為波棱瓜產(chǎn)量預(yù)測及生產(chǎn)調(diào)控提供科學(xué)依據(jù)。

1 材料

以中國藏學(xué)研究中心從西藏采集的波棱瓜種子為試驗(yàn)材料。經(jīng)中國醫(yī)學(xué)科學(xué)院藥用植物研究所鑒定，為葫蘆科植物波棱瓜Herpetospermumpedunculosum(Ser.) Baill.的種子。

2 方法

2.1 試驗(yàn)設(shè)計(jì)

2.1.1 試驗(yàn)地點(diǎn) 試驗(yàn)地設(shè)在西藏林芝市八一鎮(zhèn)核桃村，尼洋河畔，中國藏學(xué)研究中心藏藥研究基地。位于29°50′N，93°25′E，海拔2900 m。年降雨量650 mm，年均氣溫8.7 ℃，年均日照2 022.2 h，無霜期180 d，屬于溫暖半濕潤氣候。試驗(yàn)小區(qū)選擇在地勢高崗、沙質(zhì)土壤的向陽地塊，土壤pH 6.5。

2.1.2 試驗(yàn)區(qū)整地與小區(qū)設(shè)置 試驗(yàn)地每667 m2施腐熟有機(jī)肥2500 kg做底肥，對種植地翻耕30 cm、平整，整理為0.8 m寬的高畦，溝寬0.3 m。播種前灌透水，畦面表層10 cm土壤相對濕度達(dá)70% 時(shí)播種。播種方法：穴播，每穴播7粒種子，覆土2～3 cm，輕壓。點(diǎn)穴距離距高壟邊緣為0.25 m，溝寬0.3 m。畦面覆蓋地膜。每畦中，每6 m間隔1 m作為小區(qū)間的隔離帶，每2畦并列為1個(gè)處理小區(qū)。分別設(shè)不同催芽方法、不同播期、不同密度、不同施肥方案種植小區(qū)。每個(gè)小區(qū)12 m2(2 m×6 m)，小區(qū)間的隔離帶為1 m。隨機(jī)排列，共73個(gè)小區(qū)。

催芽方法包括：直播；45 ℃溫水浸種48 h；青霉素200 mg·L-1浸種；赤霉素200 mg·L-1浸種；ABT生根粉6號(hào)200 mg·L-1浸種。浸種12 h后，自來水沖洗2次。播期：2016年4月13日、5月10日、6月10日。密度：分別按穴距35 cm×70 cm和50 cm×70 cm進(jìn)行穴播，折合密度為40 816穴/hm2和28 571穴/hm2。施肥方案共3種，分別為：對照組常規(guī)施用復(fù)合肥(N∶P∶K=18∶12∶5，總質(zhì)量分?jǐn)?shù)35%)45 g·m-2、EM益生菌肥葉面噴施、追施磷肥(過磷酸鈣，P2O5質(zhì)量分?jǐn)?shù)12%)45 g·m-2。

2.2 采收和產(chǎn)量測定

按小區(qū)分別于7—10月份果實(shí)成熟期(以果實(shí)開始變黃褐色、裂口為標(biāo)準(zhǔn))陸續(xù)采收果實(shí)，自然陰干脫粒、凈化后分別測定產(chǎn)量。

2.3 數(shù)據(jù)

本研究共73組數(shù)據(jù)，觀測總數(shù)n=73。獲得種子產(chǎn)量與相應(yīng)技術(shù)措施(播期、催芽方法、雌雄比、密度、施肥方案、農(nóng)藝措施)間變化數(shù)據(jù)。

其中Y為波棱瓜種子產(chǎn)量，kg·hm-2；X1為播期，4月、5月、6月分別賦值為1、 2、3；X2為催芽方法，直播賦值為0，溫水浸種賦值為1，青霉素浸種賦值為2，赤霉素浸種賦值為3，ABT生根粉6號(hào)浸種賦值為4；X3為密度，設(shè)定為：40 816穴/hm2，28 571穴/hm2；X4為施肥方案，常規(guī)施肥賦值為0，EM益生菌肥賦值為1，追施磷肥賦值為2；X5為農(nóng)藝措施，常規(guī)賦值為0；疏花打頂賦值為1；X6為雌株比例(%)，閾值15～50。

2.4 計(jì)算分析方法

采用數(shù)量化方法Ⅰ進(jìn)行計(jì)算分析。數(shù)量化方法Ⅰ屬于一般線性模型，既可以處理定性因子，也可以處理定量因子，通過F檢驗(yàn)來篩選對因變量有顯著性影響的因子，建立因變量對自變量回歸模型[6-8]。

3 結(jié)果與分析

3.1 波棱瓜子產(chǎn)量預(yù)測數(shù)學(xué)模型建立

3.1.1 波棱瓜子產(chǎn)量與相應(yīng)因子變化匯總及定性因子數(shù)量化

根據(jù)實(shí)測數(shù)據(jù)和試驗(yàn)設(shè)計(jì)，將各因子按數(shù)量化方法Ⅰ進(jìn)行數(shù)據(jù)轉(zhuǎn)換。

表1 波棱瓜子產(chǎn)量各對應(yīng)因子水平量化賦值

3.1.2 顯著性因子篩選

1)方差分析

方差分析結(jié)果見表2。

表2 波棱瓜子產(chǎn)量多因子方差分析

注：R2=0.924 9(調(diào)整R2=0.911 4)。

2)剔出截距、X2(催芽方法)、X3(密度)和X5(農(nóng)藝措施)等對產(chǎn)量影響不顯著的因子，進(jìn)一步方差分析見表3。

表3 波棱瓜子產(chǎn)量篩選因子方差分析

注：R2=0.917 3(調(diào)整R2=0.911 1)。

由此表明：對波棱瓜子產(chǎn)量有顯著性影響的因子主要有播期(X1)、施肥(X4)和雌株比例(X6)，此3個(gè)因子對產(chǎn)量差異的解釋(擬合度)達(dá)到0.917 3。

3)確定波棱瓜子產(chǎn)量預(yù)測模型(解集)

以上分析表明，經(jīng)多因子篩選，確定影響(Y)的貢獻(xiàn)值最大的關(guān)鍵因子分別為：播期(X1)、施肥(X4)和雌株比例(X6)。

由此建立的預(yù)測擬合模型(解集)為：

Y=X1([X1=1]=226.0；[X1=2]=311.83；[X1=3]=-207.13)+X4([X4=0]=-189.33；[X4=1]=41.84；[X4=2]=0)+20.53X6

3.2 模型分析

按上述預(yù)測模型推算，西藏林芝地區(qū)海拔2900 m左右的沙質(zhì)壤土大田種植波棱瓜，在X1=2，X4=1的條件下，X6值越高，產(chǎn)量(Y)值越大。即選擇5月播期，花果期施用EM益生菌肥，雌株比例(15%～50%)越大，波棱瓜子產(chǎn)量越高，當(dāng)雌株比例達(dá)到50%時(shí)，波棱瓜子產(chǎn)量最高，可達(dá)1 380.17 kg·hm-2。

4 討論與結(jié)論

本研究通過運(yùn)用線性回歸分析模型(數(shù)量化理論Ⅰ方法)將定性因子賦值，實(shí)現(xiàn)量化分析，建立產(chǎn)量(Y)與各栽培技術(shù)措施(催芽方法、播期、雌株比、密度、施肥、農(nóng)藝措施)間的關(guān)系模型，并篩選出影響波棱瓜子產(chǎn)量的關(guān)鍵因子，提出了基于數(shù)量化Ⅰ理論的波棱瓜子產(chǎn)量預(yù)測模型的建模方法。

數(shù)量化理論Ⅰ最早由日本的林知己夫提出，其核心是將定性因子引入統(tǒng)計(jì)分析研究，通過實(shí)現(xiàn)1個(gè)或多個(gè)定性指標(biāo)進(jìn)行數(shù)據(jù)定量轉(zhuǎn)換，建立要素間的回歸數(shù)學(xué)模型[6-8]，解決相應(yīng)技術(shù)問題。近年來，數(shù)量化理論Ⅰ在生物學(xué)相關(guān)領(lǐng)域的定性因子研究方面得以應(yīng)用，通過生物數(shù)學(xué)模型的建立，使得基于多因子的篩選更為效率、科學(xué)，減少試驗(yàn)工作量。波棱瓜子產(chǎn)量受多項(xiàng)栽培因子影響，利用數(shù)量化理論Ⅰ的方法，經(jīng)對多因子回歸分析，可以篩選出均達(dá)顯著水平的關(guān)鍵性栽培因子。應(yīng)用該方法，本研究獲得對波棱瓜子產(chǎn)量貢獻(xiàn)值最大關(guān)鍵栽培因子依次為：雌株比例(X6)、播期(X1)和施肥(X4)，為波棱瓜中藥材生產(chǎn)質(zhì)量管理規(guī)范(GAP)栽培技術(shù)要點(diǎn)的確定提供了科學(xué)依據(jù)。對于在生長管理過程中的限制性措施應(yīng)在種植時(shí)充分考慮，特別是雌雄比的調(diào)控。

本研究建立的波棱瓜產(chǎn)量預(yù)測模型，可對波棱瓜因技術(shù)措施變化帶來的產(chǎn)量變化實(shí)現(xiàn)科學(xué)預(yù)測和評估。應(yīng)用該模型，還可對生產(chǎn)調(diào)控的關(guān)鍵技術(shù)措施進(jìn)一步篩選，如：優(yōu)化提高雌株比例的技術(shù)方案、EM菌肥使用量和噴施技術(shù)，明確于波棱瓜標(biāo)準(zhǔn)生產(chǎn)規(guī)程中，為波棱瓜高效生產(chǎn)體系的建立提供依據(jù)。