程薪宇，徐海軍，郭夢(mèng)橋，2，孫宇峰，曹焜，潘冬梅，魏國(guó)江，王曉飛，2，張旭，2，王曉楠*

（1.黑龍江省科學(xué)院大慶分院，黑龍江大慶163319；2.黑龍江省科學(xué)院高技術(shù)研究院，哈爾濱150020）

大麻（Cannabis sativa L.）又名寒麻、線麻、火麻等［1］，為大麻科（Cannabaceae）大麻屬一年生直立草本植物，原產(chǎn)錫金、不丹等國(guó)，現(xiàn)各國(guó)均有野生或栽培［2］。大麻是一種重要的經(jīng)濟(jì)作物，根據(jù)用途可分為纖維用、籽用和藥用3類(lèi)，現(xiàn)有品種150個(gè)左右［1］。但大麻同時(shí)也是一種具有制幻作用的毒品［3］，其主要毒性成分為四氫大麻酚（Tetrahydrocannabinol，THC）［4－7］。攝入過(guò)量的 THC可導(dǎo)致擬精神病癥狀，諸如人格解體、幻聽(tīng)、幻視、幻覺(jué)，對(duì)記憶力損害等［4］。國(guó)際上將THC含量低于0.3%的品種稱(chēng)為纖維大麻或工業(yè)大麻［8］，而將THC含量高于0.3%的品種稱(chēng)為藥用或毒品大麻（Marijuana、Hashish）［9］。目前，低毒大麻品種選育是工業(yè)大麻育種的重要方向之一［10－12］，而在低毒品種選育過(guò)程中，主要利用組織化學(xué)［9］或色譜法［13－14］來(lái)檢測(cè)毒性高低，但這些方法的高成本及低效性卻大大延長(zhǎng)了新品種選育周期。

植物的形態(tài)特征能在一定程度上反映作物的品質(zhì)，可作為經(jīng)濟(jì)作物品種選育的重要參考依據(jù)［15－17］。同時(shí)，形態(tài)學(xué)特征亦可作為大麻優(yōu)良品種選育的依據(jù)，如株高、莖粗常被用作高纖品種選育的參考依據(jù)［18］，穗緊湊程度及花性狀則主要作為低毒品種選育的依據(jù)［19］。以上形態(tài)學(xué)特征主要通過(guò)觀察快速生長(zhǎng)期及隨后發(fā)育時(shí)期的植株來(lái)獲得，大麻優(yōu)勢(shì)植株的首次篩選并非始于快速生長(zhǎng)期，而是始于對(duì)幼苗期植株的間苗過(guò)程。然而目前尚未見(jiàn)關(guān)于大麻幼苗期優(yōu)勢(shì)植株篩選依據(jù)方面的報(bào)道。Furr等［20］曾通過(guò)組織化學(xué)染色法對(duì)不同時(shí)期、不同部位的大麻植株進(jìn)行染色處理，發(fā)現(xiàn)大麻酚類(lèi)物質(zhì)主要集中于乳汁管及腺毛中，這一結(jié)論在后期的研究中也得到了證實(shí)［21－23］。由此可見(jiàn)，大麻腺毛及乳汁管數(shù)量很可能成為反映大麻酚類(lèi)物質(zhì)含量或毒性高低的間接判斷依據(jù)。但這些特征在田間實(shí)際生產(chǎn)中無(wú)法輕易獲得，需在找到與之相關(guān)性強(qiáng)且便于觀察的特征的前提下，方有可能將其應(yīng)用至田間工作中。因此，本研究選取4個(gè)大麻品種火麻1號(hào)、火麻3號(hào)、金刀－15和格列西亞以及兩個(gè)高世代種資源20170114和資源CA－4幼苗期（6葉期）完全展開(kāi)的第一真葉作為研究材料，對(duì)其外形、毛狀體特征和葉柄乳汁管數(shù)進(jìn)行觀察統(tǒng)計(jì)，并對(duì)其進(jìn)行相關(guān)分析及因子分析，旨在找到與大麻腺毛及乳汁管數(shù)量相關(guān)且能直接應(yīng)用于生產(chǎn)中的幼苗形態(tài)學(xué)性狀，為減輕大麻低毒品種育種工作量提供形態(tài)學(xué)依據(jù)。

1 材料與方法

1.1 試驗(yàn)材料

本研究以4個(gè)大麻品種火麻1號(hào)、火麻3號(hào)、金刀－15和格列西亞以及兩個(gè)高世代種資源20170114和資源CA－4為試驗(yàn)材料。其中金刀－15和格列西亞引自烏克蘭，其余材料由黑龍江省科學(xué)院大慶分院提供。試驗(yàn)于2017年5月在黑龍江省大慶市東風(fēng)農(nóng)場(chǎng)黑龍江省科學(xué)院大慶分院試驗(yàn)基地進(jìn)行，土壤為鹽堿化草甸土，其理化性質(zhì)為：pH 8.48，堿解氮220.13 mg／kg，速效磷24.57 mg／kg，速效鉀349.32 mg／kg，有機(jī)質(zhì) 17.6 mg／kg，測(cè)土配方施肥；播種后立即進(jìn)行封閉除草，視實(shí)際情況采取灌溉及病蟲(chóng)害綜合防治等措施。播種密度為450粒／m2。

1.2 試驗(yàn)方法

1.2.1 形態(tài)學(xué)試驗(yàn)材料采集

每種供試材料隨機(jī)采集5株幼苗期（6葉期）植株，重復(fù)3次，摘取第一真葉（第一對(duì)真葉中的一片），固定于FAA固定液（甲醛、冰醋酸、70%乙醇按1∶1∶18比例混合而成）中24 h以上。

1.2.2 觀察方法

待大麻葉完全固定后，按以下方法進(jìn)行觀察：

葉外形——將葉平展后，用相機(jī)拍照，并用Image J進(jìn)行葉柄長(zhǎng)度、葉片長(zhǎng)度、葉片寬度、葉面積的測(cè)量，柄葉比為葉柄長(zhǎng)度／葉片長(zhǎng)度。

葉柄乳汁管數(shù)——用雙面刀片制作大麻葉柄徒手切片，并將其制作成臨時(shí)裝片后，用萊卡DM4000B進(jìn)行觀察，Las X軟件計(jì)數(shù)。

毛狀體相關(guān)性狀——將葉片直接制成臨時(shí)裝片，用萊卡DM4000B在葉片最寬處附近隨機(jī)選取3個(gè)遠(yuǎn)離葉脈的區(qū)域進(jìn)行葉片腹背面（近軸面及遠(yuǎn)軸面）毛狀體特征的觀察，并用Las X軟件統(tǒng)計(jì)和測(cè)量非腺毛及腺毛的密度（每平方毫米中的毛狀體數(shù)量，重復(fù)觀察3個(gè)視野，取平均值）、總數(shù)（由密度值與葉面積值相乘所得）和腺毛直徑（視野中所有腺毛直徑的平均值，重復(fù)觀察3個(gè)視野，取平均值）。

1.2.3 統(tǒng)計(jì)分析

利用SPSS 20.0對(duì)大麻形態(tài)學(xué)觀察結(jié)果進(jìn)行相關(guān)性分析和因子分析。

2 結(jié)果與分析

2.1 葉形態(tài)學(xué)性狀分析

供試材料16個(gè)形態(tài)學(xué)性狀的平均值及變異情況見(jiàn)表1。由表1可知，葉片長(zhǎng)度、葉片寬度、葉面積、腹面腺毛直徑和背面腺毛直徑的變異系數(shù)相對(duì)較小，分別為11.99%、12.50%、22.28%、9.87%、9.48%，表明與大麻葉片外形及腺毛直徑相關(guān)的性狀相對(duì)穩(wěn)定。而葉柄外形及毛狀體數(shù)量相關(guān)的性狀在大麻個(gè)體間的差異卻較大，其中葉柄長(zhǎng)度、柄葉比、葉柄乳汁管數(shù)、腹面非腺毛密度、腹面非腺毛總數(shù)、腹面腺毛密度、腹面腺毛總數(shù)、背面非腺毛密度、背面腺毛密度、背面腺毛總數(shù)的變異系數(shù)分別 42.27%、40.91%、54.07%、52.63%、60.15%、77.98%、82.42%、50.86%、102.17%、103.33%，表明這些性狀在研究大麻個(gè)體差異及變化規(guī)律方面具有更高的實(shí)用價(jià)值。

2.2 大麻葉形態(tài)學(xué)性狀相關(guān)性分析

對(duì)大麻幼苗葉外形、毛狀體及乳汁管性狀進(jìn)行相關(guān)性分析（詳見(jiàn)表2），結(jié)果顯示葉片外形性狀間存在極顯著正相關(guān)關(guān)系，其中葉面積與葉片長(zhǎng)度和葉片寬度的相關(guān)系數(shù)分別為0.892和0.851；與葉柄外形相關(guān)的兩個(gè)性狀，葉柄長(zhǎng)度和柄葉比之間的相關(guān)系數(shù)為0.950；與毛狀體數(shù)量有關(guān)的性狀間的相關(guān)性均達(dá)到極顯著水平，其中腹面非腺毛密度與腹面腺毛密度、腹面腺毛總數(shù)、背面腺毛密度和背面腺毛總數(shù)的相關(guān)系數(shù)依次為0.893、0.864、0.872和0.880。腺毛直徑相關(guān)的兩個(gè)性狀腹面腺毛直徑和背面腺毛直徑與本研究觀察的其它大麻葉性狀相關(guān)性較弱，而二者彼此間卻呈極顯著正相關(guān)，相關(guān)系數(shù)0.734。由此可見(jiàn)，選擇一些指征性良好的指標(biāo)來(lái)代替絕大多數(shù)觀測(cè)指標(biāo)是可行的。

2.3 葉形及毛狀體性狀因子分析

2.3.1 葉形性狀因子分析

由上述結(jié)果可知，大麻葉各類(lèi)特征內(nèi)部性狀彼此間常存在較強(qiáng)的相關(guān)性，經(jīng)主成分進(jìn)一步分析，可知：葉形相關(guān)性狀可提取為兩個(gè)主因子，貢獻(xiàn)率分別為64.739%和27.019%，累計(jì)貢獻(xiàn)率為91.758%，說(shuō)明這兩個(gè)主因子可以對(duì)葉形性狀做出很好的解釋。但葉片長(zhǎng)度、柄葉比及葉柄長(zhǎng)在兩個(gè)主因子中的載荷量均較大，為了更好的定義因子并解釋向量，將成分矩陣進(jìn)行正交旋轉(zhuǎn)，將上述3個(gè)變量分別歸入2個(gè)公因子中，結(jié)果表明，第一公因子主要包括葉面積、葉片長(zhǎng)度和葉片寬度，其荷載系數(shù)分別為0.966、0.927、0.811，主要反映的是葉片外形特征，將其命名為“葉片外形因子”；柄葉比及葉柄長(zhǎng)度在公因子2上有較大的載荷，載荷量分別為0.996和0.930，主要反應(yīng)的是葉柄外形相關(guān)特征，因此命名為“葉柄外形因子”。

?

表3 葉形特征因子分析Tab.3 The factor analysis of leaf shape characteristics

2.3.2 毛狀體性狀因子分析

經(jīng)主成分進(jìn)一步分析，大麻毛狀體相關(guān)性狀可提取兩個(gè)主因子，貢獻(xiàn)率分別為73.691%和14.823%，累計(jì)貢獻(xiàn)率為88.514%，說(shuō)明這兩個(gè)主因子可以對(duì)毛狀體相關(guān)性狀做出很好的解釋。但腹面腺毛直徑在兩個(gè)主因子中的載荷量均較大，為了更好地定義因子并解釋向量，將成分矩陣進(jìn)行正交旋轉(zhuǎn)，結(jié)果表明，腹面腺毛總數(shù)、腹面腺毛密度、腹面非腺毛總數(shù)、腹面非腺毛密度、背面腺毛總數(shù)、背面非腺毛總數(shù)、背面腺毛密度和背面非腺毛密度8個(gè)變量在公因子1上的載荷量較高，依次為0.949、0.944、0.938、0.934、0.930、0.898、0.892和0.873。這些變量均與毛狀體數(shù)量相關(guān)，因此將公因子1命名為“毛狀體數(shù)量因子”。公因子2中僅有背面腺毛直徑和腹面腺毛直徑兩個(gè)變量，在公因子2上的載荷量分別為0.930和0.894，因此將其命名為“腺毛直徑因子”。

表4 毛狀體特征因子分析Tab.4 The factor analysis of trichome characteristics

2.4 葉柄乳汁管數(shù)與新因子間的相關(guān)性分析

將葉柄乳汁管數(shù)性狀與之前因子分析提取的4個(gè)公因子變量進(jìn)行相關(guān)性分析（詳見(jiàn)表5）。結(jié)果顯示葉柄乳汁管數(shù)與葉柄外形因子及毛狀體數(shù)量因子均存在極顯著的相關(guān)性，相關(guān)系數(shù)分別為0.697和0.721。而葉柄外形因子與毛狀體數(shù)量因子彼此間的相關(guān)性也達(dá)到了極顯著水平，相關(guān)系數(shù)0.771。葉片外形因子和腺毛直徑因子與葉柄乳汁管數(shù)和其它兩個(gè)公因子的相關(guān)性較弱。這表明可通過(guò)觀察葉柄外形來(lái)估測(cè)葉柄乳汁管數(shù)和毛狀體數(shù)量或通過(guò)統(tǒng)計(jì)毛狀體數(shù)量來(lái)估測(cè)葉柄乳汁管數(shù)。但卻無(wú)法通過(guò)觀察葉片外形和腺毛直徑來(lái)估測(cè)其它性狀情況。

表5 葉柄乳汁管數(shù)與新因子間的相關(guān)性分析Tab.5 The correlation analysis between petiole laticifer number and new factors

3 討論

植物形態(tài)特征除作為系統(tǒng)學(xué)及分類(lèi)學(xué)研究［24］的依據(jù)外，在作物優(yōu)良品種選育過(guò)程中也發(fā)揮著重要作用［15－17］。目前大麻低毒品種幼苗期主要保留目測(cè)長(zhǎng)勢(shì)較好且整齊的植株［20，25］。但此種篩選方式并不能有效地降低大麻低毒品種選育的周期。而現(xiàn)有研究結(jié)果表明，大麻酚類(lèi)物質(zhì)主要分布于大麻的腺毛及乳汁管中［20－23］。因此，了解腺毛及乳汁管數(shù)量與其它大麻形態(tài)學(xué)性狀間的關(guān)系，可為大麻低毒品種選育過(guò)程中育種中間材料的篩選提供新的參考依據(jù)。

本研究結(jié)果顯示腺毛數(shù)量相關(guān)性狀和非腺毛數(shù)量有關(guān)的性狀共同組成了毛狀體數(shù)量因子。毛狀體數(shù)量因子中腹面腺毛總數(shù)載荷量最大，為0.949，其次為腹面腺毛密度（0.944）、腹面非腺毛總數(shù)（0.938）和腹面非腺毛密度（0.934），這4個(gè)變量在毛狀體數(shù)量因子中的載荷量較為接近。而這4個(gè)變量彼此間的相關(guān)系數(shù)均大于0.850。綜合以上信息后，認(rèn)為在實(shí)際生產(chǎn)過(guò)程中，觀察這4個(gè)變量中任一變量，均可以在一定程度上反映其它3個(gè)變量的情況。再考慮到腺毛在葉上的總數(shù)和密度遠(yuǎn)不及非腺毛多（前者在葉片腹面的平均總數(shù)和密度分別為3154.76個(gè)和5.63個(gè)／mm2，后者則分別為31253.99個(gè)和56.81個(gè)／mm2），以及腹面非腺毛密度與所有腺毛數(shù)量相關(guān)性狀的相關(guān)系數(shù)均大于0.860，從實(shí)用性角度出發(fā)，認(rèn)為可通過(guò)觀察腹面非腺毛密度這一性狀來(lái)推測(cè)腺毛總數(shù)及密度的高低，并可由其來(lái)代表毛狀體數(shù)量因子。

由表5可知，葉柄乳汁管數(shù)與葉柄外形因子及毛狀體數(shù)量因子間均為極顯著正相關(guān)，相關(guān)系數(shù)分別為0.697和0.721。認(rèn)為可通過(guò)這兩個(gè)因子來(lái)推測(cè)葉柄乳汁管數(shù)。而葉柄外形因子中柄葉比對(duì)公因子的貢獻(xiàn)最大，其因子載荷量達(dá)到0.996，且與公因子中的另一個(gè)變量葉柄長(zhǎng)的相關(guān)系數(shù)也達(dá)到了0.950，因此在生產(chǎn)中可直接用柄葉比代表葉柄外形因子使用。再結(jié)合之前提到的用腹面非腺毛密度來(lái)代表毛狀體數(shù)量因子，最終認(rèn)為，可通過(guò)觀察大麻葉的柄葉比或腹面非腺毛密度來(lái)推測(cè)葉柄乳汁管數(shù)的高低。

根據(jù)表5的分析結(jié)果可知，葉片外形因子和腺毛直徑因子與葉柄乳汁管數(shù)和毛狀體數(shù)量間無(wú)較強(qiáng)的相關(guān)性，這意味著無(wú)法通過(guò)葉片外形和腺毛直徑特征來(lái)反映乳汁管及毛狀體的相關(guān)特征。

通過(guò)本次研究，最終篩選出兩個(gè)較為容易觀察的性狀——柄葉比和葉腹面非腺毛密度。這兩個(gè)性狀可在一定程度上間接反映與大麻腺毛及葉柄乳汁管數(shù)量的多寡，為進(jìn)一步研究大麻形態(tài)學(xué)特征與其所含大麻酚類(lèi)物質(zhì)含量間的關(guān)系提供了新的形態(tài)學(xué)依據(jù)。

參考文獻(xiàn)：

［1］王德珠，陳建，李宏俊.大麻纖維及其應(yīng)用［J］.中國(guó)纖檢，2012（5）：81－83.

［2］中國(guó)科學(xué)院中國(guó)植物志.中國(guó)植物志，23卷，第一分冊(cè)［M］.北京：科學(xué)出版社，1998，223.

［3］旺堆，旦增旺扎，汪荀.簡(jiǎn)述大麻毒品的社會(huì)危害性及其利用GC／MS檢驗(yàn)大麻成分的分析技術(shù)［J］.西藏科技，2015（4）：19－20.

［4］何錦風(fēng)，杜軍強(qiáng)，陳天鵬.漢麻葉的生物活性成分研究現(xiàn)狀［J］.中國(guó)食品學(xué)報(bào)，2011，11（8）：133－140.

［5］Aizpurua－Olaizola O，Omar J，Navarro P，etal.Identification and quantification of cannabinoids in Cannabissativa L.plants by high performance liquid chromatography－mass spectrometry［J］.Analytical and Bioanalytical Chemistry，2014，406（29）：7549－7560.

［6］Benjamin D B，Kevin M，Verstraete A G，etal.Evolution of the content of THC and othermajor cannabinoids in drugmilype cannabis cuttings and seedlings during growth of plants［J］.Journal of Forensic Sciences，2012，57（4）：918－922.

［7］Elsohly M A，Slade D.Chemical constituents ofmarijuana：the complexmixture of natural cannabinoids［J］.Life Sciences，2005，78（5）：539－548.

［8］Salentijn EM J，Zhang Q，Amaducci S，et al.New developments in fiber hemp（Cannabis sativa L.）breeding［J］.Industrial Crops＆Products，2015，68（3）：32－41.

［9］盧延旭，董鵬，崔曉光，等.工業(yè)大麻與毒品大麻的區(qū)別及其可利用價(jià)值［J］.中國(guó)藥理學(xué)通報(bào)，2007，23（8）：1112－1114.

［10］宋憲友，李江，吳廣文，等.低毒雌雄同株大麻新品種USO－31及高產(chǎn)栽培技術(shù)［J］.中國(guó)麻業(yè)科學(xué)，2007（4）：201－203.

［11］宋憲友，張利國(guó)，房郁妍，等.大麻新品種“龍大麻1號(hào)”的選育［J］.中國(guó)麻業(yè)科學(xué)，2011，33（3）：109－111.

［12］粟建光，陳基權(quán)，謝小美.大麻育種現(xiàn)狀與前景［J］.中國(guó)麻業(yè)科學(xué)，2006，28（4）：212－217.

［13］翟晚?xiàng)?，張春水，高利?高效液相色譜檢測(cè)大麻中△～9－四氫大麻酚的分析方法研究［J］.分析測(cè)試學(xué)報(bào)，2012，31（11）：1379－1384.

［14］崔廣東.低毒工業(yè)大麻葉莖化學(xué)成分的分離以及抑菌活性物質(zhì)的跟蹤和篩選［D］.蘭州：蘭州理工大學(xué)，2008.

［15］佟屏亞.從植株形態(tài)指標(biāo)評(píng)價(jià)玉米育種方向［J］.玉米科學(xué)，2006，14（6）：1－3.

［16］陳亮.水稻根系育種的意義與前景［J］.南方農(nóng)業(yè)，2016，10（24）：151－152.

［17］眭書(shū)祥，李?lèi)?ài)國(guó)，朱青竹，等.黃河流域棉花形態(tài)育種探討［J］.河北農(nóng)業(yè)科學(xué)，2012，16（7）：75－77.

［18］郭媛，王玉富，邱財(cái)生，等.干旱脅迫對(duì)不同大麻品種生理特性和生長(zhǎng)的影響研究初報(bào)［J］.中國(guó)麻業(yè)科學(xué)，2011，33（5）：235－239.

［19］郭鴻彥，胡學(xué)禮，陳裕，等.早熟籽用型工業(yè)大麻新品種云麻2號(hào)的選育［J］.中國(guó)麻業(yè)科學(xué)，2009，31（5）：285－287.

［20］Furr M，Mahlberg PG.Histochemical analyses of laticifers and glandular trichomes in Cannabis sativa［J］.Journal of Natural Products，1981，44（2）：153－159.

［21］Sirikantaramas S，Taura F，Tanaka Y，et al.Tetrahydrocannabinolic acid synthase，the enzyme controllingmarijuana psychoactivity，is secreted into the storage cavity of the glandular trichomes［J］.Plant＆Cell Physiology，2005，46（9）：1578－82.

［22］陳璇，楊明，郭鴻彥.大麻植物中大麻素成分研究進(jìn)展［J］.植物學(xué)報(bào)，2011，46（2）：197－205.

［23］Paul G.Mahlberg，Eun Soo Kim.Accumulation of Cannabinoids in Glandular Trichomes of Cannabis（Cannabaceae）［J］.Journal of Industrial Hemp，2004，9（1）：15－36.

［24］楊永紅.大麻屬植物分類(lèi)研究史［J］.中國(guó)麻業(yè)科學(xué)，2003，25（1）：9－11.

［25］胡學(xué)禮，郭鴻彥，劉旭云，等.云南工業(yè)大麻品種在黑龍江大興安嶺地區(qū)的適應(yīng)性研究［J］.西南農(nóng)業(yè)學(xué)報(bào)，2012，25（3）：838－841.