基于黃酮合成通路基因表達特性的木薯對二斑葉螨防御機理初探

喬陽 梁曉 陳青 伍春玲 劉迎 劉小強 稅軍 姚曉文

關鍵詞：二斑葉螨；抗、感螨木薯品種；黃酮合成通路；基因表達

木薯（ManihotesculentaCrantz）是大戟科（Euphorbiaceae）木薯屬（Manihot）作物，原產于亞馬遜南部盆地，目前廣泛種植于全球亞熱帶和熱帶地區[1]，具有高產、高淀粉含量、耐貧耐旱等特點[2]，為世界上100多個國家近10億人口提供碳水化合物來源[3]。木薯還是重要的工業原料、能源作物、飼料作物和救荒作物，在淀粉糖合成、醫藥、化工、造紙和食品等領域發揮重要作用[4-5]。

二斑葉螨（TetranychusurticaeKoch）是世界危險性害螨，是危害木薯最嚴重的四大有害生物之一[6]，其嚴重為害會使木薯減產50%～70%，甚至絕收[7]。目前，化學殺螨劑依然是控制二斑葉螨最常用的方法，但該螨于木薯種植后6～8個月暴發成災，而此時木薯地已封行，不僅導致藥劑的有效利用率降低，還將造成害螨抗藥水平上升，產品和產地環境安全風險加劇等問題[8-9]。因此，尋求綠色、高效、環境友好的害蟲防治策略才能夠有效支撐木薯產業的可持續發展。

抗蟲品種的選育與利用是實現害蟲綠色綜合防控的重要手段[10-11]，研究作物抗蟲機理能為抗蟲品種的選育提供重要的理論依據。研究表明，植物能產生多種次生代謝物來防御生物和非生物脅迫[12-13]，其中，黃酮類化合物在作物形成抗蟲性中發揮重要作用[14]，它可以抑制昆蟲的攝食、發育和產卵[15]，黃酮類化合物及其合成調控基因的抗蟲作用已見于玉米抗棉鈴蟲（HelicoverpaarmigeraHübner）、茶樹抗茶尺蠖（EctropisobliquaProut）、水稻抗褐飛虱（NilaparvatalugensStal）等相關研究[16-18]。然而迄今為止，尚未見黃酮合成通路基因在木薯抗二斑葉螨中的功能報道。因此，本研究擬分析比較二斑葉螨取食不同抗性水平木薯品種前后黃酮合成通路基因的表達量差異，以初步明確與木薯抗二斑葉螨相關的黃酮合成通路關鍵基因，為深入開展基于黃酮代謝通路的抗螨木薯新種質的創制與分子設計育種提供重要的理論依據。

1材料與方法

1.1材料

1.1.1供試二斑葉螨參考伍牧鋒等[19]的方法進行二斑葉螨的室內飼養。選擇發育歷期相同、大小一致的雌成螨進行試驗。

1.1.2供試木薯遺傳穩定的抗螨木薯品種C1115、緬甸（Myanmar）、SC9，感螨木薯品種SC205、面包（Bread）、BRA900[20]由中國熱農業科學院熱帶作物品種資源研究所國家木薯種質資源圃提供。參考韓志玲等[21]的方法進行木薯的室內種植，待木薯生長約60d后，選擇長勢一致的健康木薯植株進行害螨接種。

1.1.3候選黃酮合成通路基因根據KEGG網站（https：//www.kegg.jp/kegg/）中已發布的黃酮類代謝物合成通路，繪制木薯黃酮類代謝物合成通路示意圖（圖1），選取植物黃酮合成調控研究中最為廣泛的7個基因作為候選基因，分別為CHS（查爾酮合酶）、PGT1（根皮苷合酶）、C3'H（5-O-香豆?；?D-喹啉3'-單加氧酶）、FLS（黃酮醇合酶）、F3H（柚皮素3-雙加氧酶）、LAR（白花青素還原酶）、CYP93B_16（黃酮合酶II）。

1.2方法

1.2.1二斑葉螨接種與取樣分別選取室內種植60d，長勢一致，健康的抗螨木薯品種C1115、緬甸（Myanmar）、SC9，感螨木薯品種SC205、面包（Bread）、BRA900進行害螨接種。每株木薯分別選取中部的3張葉片，用實驗室自制的吸蟲器小心吸取發育一致的二斑葉螨雌成螨，接種于葉背面（50頭/片），并以羊毛脂涂抹在葉柄基部以防止二斑葉螨逃逸。采集未受螨害（0d）、螨害1、4d后的葉片用于RNA提取。不同木薯品種不同螨害時間處理均設置3個生物學重復。

1.2.2RNA提取及cDNA第一條鏈的合成參照多糖多酚植物總RNA提取試劑盒（TIANGEN，美國）說明書提取木薯葉片RNA，取去除gDNA后的RNA樣品1.0μg進行cDNA第一鏈的合成，合成方法參照RTEasyMixforqPCR試劑盒（TOLOBIO，中國）。

1.2.3實時熒光定量PCR分析根據GenBank中已經發布的木薯黃酮類代謝物合成通路7個關鍵基因序列設計qPCR引物（表1）。cDNA樣品經RNase-freeddH2O稀釋160倍后作為qPCR的模板，以木薯Meactin為內參基因[7]。qPCR反應體系的配制參照2×Q3SYBRqPCRMasterMix試劑盒（TOLOBIO，中國）。qPCR反應條件為95℃30s；95℃10s，60℃30s，40個循環。使用LightCycler?96儀器（Roche，瑞士）默認溶解曲線程序采集溶解曲線。根據LIVAK等[22]的2-ΔΔCt（ΔΔCt=ΔCt2-ΔCt1，ΔCt1為2個樣本內參基因Ct值的差值，ΔCt2為2個樣本目的基因Ct值的差值）方法計算而得，每個處理均設置3個重復，目的基因在2個樣本間的相對表達量為3個重復的平均值。

1.3數據處理

采用Excel和GradphadPrism軟件進行數據匯總整理，使用SPSSStatistics24軟件進行數據分析，采用Duncans新復極差法進行數據間的多重比較，采用皮爾遜相關性分析（Personscorrelationanalysis）方法分析基因表達量和木薯品種螨害指數間的相關性，其中螨害指數用于表征木薯品種的抗螨水平，前期研究結果表明，抗螨木薯品種C1115、SC9和Myanmar的螨害指數分別為12.50、26.60和14.40，感螨木薯品種Bread、SC205和BRA900的螨害指數分別為95.40、93.60和95.70[20]。

2結果與分析

2.1二斑葉螨取食前后抗、感螨木薯品種葉組織中CHS基因表達差異分析

比較螨害前后同一木薯品種的CHS基因相對表達量發現，與螨害前相比，感螨木薯品種Bread、SC205和BAR900的CHS基因表達量呈先顯著降低再顯著升高的趨勢，而抗螨木薯品種C1115、SC9和Myanmar的CHS相對表達量呈先顯著升高再顯著降低的趨勢（圖2A）。

比較螨害前后抗、感螨品種間CHS的表達量發現，螨害4d后，抗螨木薯品種中CHS的表達量總體上高于感螨木薯品種，而螨害1d后，抗螨木薯品種中CHS的表達量均顯著高于感螨木薯品種（圖2B）。進一步相關性分析表明（以螨害1d后數據進行分析），螨害后CHS基因的相對表達量與木薯品種抗螨性呈正相關（r=0.8684，P=0.012）。

2.2二斑葉螨取食前后抗、感木薯品種葉組織中的PGT1基因表達差異分析

比較螨害前后同一木薯品種的PGT1基因相對表達量發現，與螨害前相比，感螨木薯品種Bread、SC205和BAR900的PGT1基因表達量呈先顯著降低再顯著升高的趨勢，而抗螨木薯品種C1115、SC9和Myanmar的PGT1相對表達量呈先升高再顯著降低的趨勢（圖3A）。

比較螨害前后抗、感螨品種間PGT1的表達量發現，螨害前（0d），抗螨木薯品種中PGT1的表達量總體上高于感螨木薯品種，而螨害1d后，抗螨木薯品種中PGT1的表達量均顯著高于感螨木薯品種（圖3B），進一步相關性分析結果表明（以螨害1d后數據進行分析），螨害后PGT1基因相對表達量與木薯品種抗螨性呈正相關（r=0.9891，P=0.0001）。

比較螨害前后抗、感螨品種間F3H的表達量發現，螨害前（0d），抗、感木薯品種間F3H基因表達量差異不顯著，但抗螨木薯品種C1115的F3H表達量顯著高于其他品種，螨害4d后，抗螨木薯品種中F3H的表達量總體上高于感螨木薯品種。而螨害1d后，抗螨木薯品種C1115、SC9和Myanmar的F3H表達量均顯著高于感螨木薯品種（圖4B）。進一步相關性分析結果表明（以螨害1d后數據進行分析），螨害后F3H基因相對表達量與木薯品種抗螨性呈正相關（r=0.6438，P=0.084）。

2.3二斑葉螨取食前后抗、感木薯品種葉組織中的F3H基因表達差異分析

比較螨害前后同一木薯品種的F3H基因相對表達量發現，與螨害前相比，感螨木薯品種Bread、SC205和BAR900的F3H基因表達量呈先顯著降低再升高的趨勢，而抗螨木薯品種SC9和Myanmar的F3H基因相對表達量呈先顯著升高再顯著降低的趨勢，C1115的F3H基因相對表達量則呈顯著降低的趨勢（圖4A）。

2.4二斑葉螨取食前后抗、感木薯品種葉組織中的FLS基因表達差異分析

比較螨害前后同一木薯品種的FLS基因相對表達量發現，與螨害前相比，感螨木薯品種Bread、SC205和BAR900的FLS基因表達量呈逐漸降低的趨勢，而抗螨木薯品種C1115、SC9和Myanmar的FLS相對表達量呈先顯著升高再降低的趨勢（圖5A）。

比較螨害前后抗、感螨品種間FLS的表達量發現，螨害4d后，抗螨木薯品種中FLS的表達量總體上高于感螨木薯品種，而螨害1d后，抗螨木薯品種的FLS表達量均顯著高于感螨木薯品種（圖5B）。進一步相關性分析結果表明（以螨害1d后數據進行分析），螨害后FLS基因相對表達量與木薯品種抗螨性呈正相關（r=0.8073，P=0.026）。

2.5二斑葉螨取食前后抗、感木薯品種葉組織中的LAR基因表達差異分析

比較螨害前后同一木薯品種的LAR基因相對表達量發現，與螨害前相比，感螨木薯品種Bread和SC205的LAR基因表達量呈先降低再升高的趨勢，BAR900的LAR基因相對表達量呈升高的趨勢?？跪臼砥贩NSC9和Myanmar的LAR基因相對表達量呈先顯著升高再降低的趨勢，而C1115的LAR基因相對表達量呈降低趨勢（圖6A）。

比較螨害前后抗、感螨品種間LAR的表達量發現，螨害1d后，抗螨木薯品種SC9和Myanmar的LAR基因相對表達量均顯著高于各感螨木薯品種，而C1115的LAR基因相對表達量與感螨木薯之間無顯著差異（圖6B）。進一步相關性分析結果表明（以螨害1d后數據進行分析），螨害后LAR基因相對表達量與木薯品種抗螨性呈正相關（r=0.5722，P=0.118）。

2.6二斑葉螨取食前后抗、感木薯品種葉組織中的C3'H基因表達差異分析

比較螨害前后同一木薯品種的C3'H基因相對表達量發現，與螨害前相比，感螨木薯品種Bread和BAR900的C3'H基因相對表達量呈先顯著降低再顯著升高的趨勢，SC205的C3H基因相對表達量呈升高趨勢，抗螨木薯品種SC9和Myanmar的C3'H基因相對表達量呈先顯著升高再顯著降低的趨勢，而C1115的C3'H基因相對表達量無顯著變化（圖7A）。

比較螨害前后抗、感螨品種間C3'H的表達量發現，螨害前（0d），抗螨木薯品種中C3'H的表達量總體高于感螨木薯品種，而螨害1d后，抗螨木薯品種中C3'H的表達量均高于感螨木薯品種（圖7B）。進一步相關性分析結果表明（以螨害1d后數據進行分析），螨害后C3'H基因相對表達量與木薯品種抗螨性呈正相關（r=0.7704，P=0.037）。

2.7二斑葉螨取食前后抗、感木薯品種葉組織中的CYP93B_16基因表達差異分析

比較螨害前后同一木薯品種的CYP93B_16基因相對表達量發現，與螨害前相比，感螨木薯品種Bread、SC205和BAR900的CYP93B_16基因表達量呈先顯著降低再顯著升高的趨勢，而抗螨木薯品種SC9和Myanmar的CYP93B_16相對表達量呈先升高再顯著降低的趨勢，C1115的CYP93B_16基因相對表達量呈顯著下降的趨勢（圖8A）。

比較螨害前后抗、感螨品種間CYP93B_16的表達量發現，螨害1d后，抗螨木薯品種SC9和Myanmar的CYP93B_16基因相對表達量均顯著高于各感螨木薯品種，而C1115的CYP93B_16基因相對表達量與感螨木薯之間無顯著差異（圖8B）。進一步相關性分析結果表明（以螨害1d后數據進行分析），螨害后CYP93B_16基因相對表達量與木薯品種抗螨性呈正相關（r=0.6250，P=0.092）。

3討論

黃酮類化合物廣泛存在于植物體內，并且在抵御蟲害脅迫中發揮著重要作用[23]。體內分析研究表明，植物中的黃酮類化合物，如兒茶素、山奈酚和蘆丁可以影響害蟲的取食偏好和產卵[24-26]。落葵中的芹菜素和牡荊素可以抑制斜紋夜蛾（SpodopteralituraFabricius）幼蟲的生長發育[27]。綠原酸、槲皮素和阿魏酸等提高了花生對棉鈴蟲（Helicoverpaarmigera）和蚜蟲（AphiscraccivoraKoch）的抗性[28]。紫杉葉素和兒茶素在云杉防御樹皮甲蟲（Ipstypographus）中發揮重要作用[29]。槲皮素和蘆丁可以影響蚜蟲的口針刺吸植物的汁液[30]。通過對39個水稻品種對褐飛虱的抗性分析結果表明，黃酮類化合物含量的提高與水稻抗褐飛虱能力的增強有關[31]。體外分析結果表明，在人工飼料中摻入黃酮類化合物也會極大地影響昆蟲的生長和發育。如在較高濃度的槲皮素和綠原酸的飼喂處理下，美國白蛾（Hyphantriacunea）的死亡率分別達到了70.00%和67.50%[32]。將柚皮素和槲皮素添加到豌豆蚜（Acyrthosiphonpisum）的人工飼料中，顯著降低了蚜蟲的繁殖力，并顯著增加了其發育歷期和成蟲的死亡率[33]。

黃酮類化合物在植物體內的積累及轉化受其合成基因的調控，因此，黃酮合成通路基因的表達與植物抗蟲性有關。黃酮合成通路基因F3H、C4H和F3'5'H在馬尾松被松茸毒蛾（DasychiraaxuthaCollenette）和松毛蟲（Dendrolimus）為害時快速表達，提高了馬尾松的抗蟲性[34]。擬南芥的花青素和黃酮醇的合成受MYB轉錄因子的調控，在擬南芥中過表達MYB75導致花青素和黃酮醇水平增加，可以增強擬南芥對鱗翅目害蟲的抗性[15]。轉F3H基因水稻對褐飛虱（Nilaparvatalugens）和白背飛虱（SogatellafurciferaHorváth）的抗性顯著高于野生型水稻，遺傳分析發現F3H基因可以通過影響FLS和DFR的表達量，調控山奈酚、槲皮素和花青素的含量進而影響水稻的抗蟲性[35]。前期研究發現，二斑葉螨為害后，抗螨木薯中黃酮合成基因MeLAR和MeANR的表達量，LAR和ANR酶活以及縮合單寧含量均顯著高于感螨木薯品種，進一步在感螨野生型木薯中過表達MeLAR和MeANR基因，可以顯著提高木薯對二斑葉螨的抗性，從分子水平證明了這2個黃酮合成基因的抗螨功能[36]。然而，木薯中涉及黃酮合成的基因有多個，篩選并驗證具有調控黃酮類化合物含量的關鍵基因及其功能，能夠有效支撐抗螨木薯種質分子精準設計育種工作，而明確抗、感螨木薯品種中黃酮合成通路基因的表達差異，能夠為候選基因的篩選提供理論依據。

本研究初步揭示了較短的螨害時間內，木薯黃酮類代謝物合成通路調控基因CHS、PGT1、F3H、FLS、LAR、C3'H和CYP93B_16在感螨木薯品種中，受螨害后黃酮合成基因的表達量總體呈先降低（1d）再升高（4d）至螨害前（0d）水平的趨勢，其螨害期間的基因表達水平較螨害前并未提高，而抗螨木薯品種中黃酮合成基因的表達量總體呈先顯著升高（1d）再降低（4d）的趨勢，但均顯著高于螨害前（0d）的水平，研究表明，抗螨木薯中基因誘導上調表達的能力強于感螨木薯。進一步比較螨害前后抗、感螨木薯品種黃酮合成通路基因表達量發現，抗螨木薯品種顯著高于感螨木薯品種。害蟲為害后，抗蟲品種中黃酮合成通路基因的防御響應速度和表達量均顯著高于感蟲作物品種的現象在大豆與二斑葉螨互作[37]，高粱與玉米蚜（RhopalosiphummaidisFitch）互作[38]和水稻與飛虱互作[35]的研究中均有報道，這與本研究結果具有一致性。黃酮合成通路C3'H基因可調控綠原酸的合成，而外施綠原酸可以提高害蟲的死亡率[32]，CHS、PGT1和CYP93B_16基因分別調控根皮素、根皮苷和木犀草素的合成，因此，推測抗螨木薯品種受螨害后，可能通過誘導上述黃酮合成基因的上調，從而增加具有抗蟲功能的黃酮類代謝物含量以增強其抗螨水平。后續研究可通過干預上述基因的表達，構建轉基因木薯株系，以期直接驗證上述基因及其調控的黃酮類物質在木薯抗蟲性中的功能。

本研究一方面初步闡明了木薯中黃酮合成通路基因可能具有抗螨性功能，另一方面表明上述基因可能具有作為鑒定評價木薯抗螨性水平的分子指標的潛力。但必須指出的是，本研究僅從黃酮類代謝物合成通路基因轉錄水平初步闡述了木薯對二斑葉螨的抗性機制，而二斑葉螨為害前后，不同抗螨木薯品種的黃酮類代謝物合成通路關鍵酶和類黃酮物質含量差異尚未進行測定，其轉錄水平、酶學水平和代謝水平是否具有一致性，以及其抗螨分子功能仍需后續研究進一步證實。