類黃酮合成通路介導的木薯對木薯綿粉蚧的防御機理

耿越 陳青 梁曉 劉迎 伍春玲 陳銀華

摘??要：木薯是世界重要的糧食作物、能源作物和工業原料。木薯綿粉蚧（Phenacoccus?manihoti?Matile-Ferrero）是一種世界危險性檢疫性害蟲，培育和利用抗蟲木薯品種可以有效阻斷其定殖為害。挖掘具有抗蟲作用的次生代謝物質及其調控基因是開展抗蟲育種的重要策略之一。而類黃酮是植物抵御生物與非生物脅迫的特有次生代謝物質，但類黃酮及其合成通路基因在木薯抗木薯綿粉蚧中的功能尚不清楚。據此，本研究分析了抗蟲木薯品種（C1115、SC9、緬甸）和感蟲木薯品種（KU50、SC205、面包）被木薯綿粉蚧為害不同時間（0、1、4、8?d）后，葉片中類黃酮合成通路相關基因（CHS、CHI、FLS、LAR、DFR、F3H、CCoAOMT、C4H、C3'H、ANR）表達量以及類黃酮含量的變化情況。結果表明：木薯綿粉蚧取食后，葉組織中CCoAOMT、C3'H、ANR、C4H雖然上調表達，但與為害前相比并無顯著差異，且抗、感木薯品種間的表達量也無顯著差異；與之相對，CHS、CHI、FLS、F3H、DFR、LAR基因表達量均比為害前顯著提高，并且在相同的為害時間內，這6個基因在抗蟲木薯品種中的表達量也顯著高于感蟲木薯品種。進一步通過相關性分析發現，CHS、CHI、FLS、F3H、LAR基因的表達量與木薯品種的抗蟲性呈顯著正相關。此外，總黃酮含量測定結果表明：木薯綿粉蚧為害1?d時總黃酮含量與為害前相比均顯著上升，為害4?d后抗蟲木薯品種總黃酮含量顯著高于感蟲木薯品種。相關性分析顯示，總黃酮含量也與木薯品種的抗蟲性呈顯著正相關。因此推測，類黃酮含量的上升及其調控基因CHS、CHI、FLS、F3H、LAR在抗蟲木薯品種中的顯著上調表達與其對木薯綿粉蚧的抗性有關。本研究為深入解析類黃酮合成基因調控木薯對木薯綿粉蚧的抗蟲防御反應分子機制，以及木薯抗蟲分子設計育種提供重要的前期基礎。

關鍵詞：抗、感木薯品種；木薯綿粉蚧；類黃酮合成途徑；基因表達；抗蟲機理中圖分類號：S435.33??????文獻標識碼：A

Flavonoid?Synthesis?Pathway-mediate?Defense?Mechanism?of?Cassava?Against?Phenacoccus?manihoti

GENG?Yue1，2，3，?CHEN?Qing2，3*，?LIANG?Xiao2，3，?LIU?Ying?2，3，?WU?Chunling2，3，?CHEN?Yinhua1

1.?School?of?Tropical?Agriculture?and?Forestry，?Hainan?University，?Haikou，?Haian?570228，?China;?2.?Environment?and?Plant?Protection?Institute，?China?Academy?of?Tropical?Agricultural?Sciences?/?Key?Laboratory?of?Integrated?Pest?Management?on?Tropical?Crops，?Ministry?of?Agriculture?and?Rural?Affairs?/?Hainan?Engineering?Research?Center?for?Biological?Control?of?Tropical?Crops?Diseases?and?Insect?Pests，?Haikou，?Hainan?571101，?China;?3.?Sanya?Research?Academy，?Chinese?Academy?of?Tropical?Agriculture?Science?/?Hainan?Key?Laboratory?for?Biosafety?Monitoring?and?Molecular?Breeding?in?Off-Season?Reproduction?Regions，?Sanya，?Hainan?572000，?China

Abstract：?Cassava?is?an?important?food?crop，?energy?crop?and?industrial?raw?material.?Cassava?mealybug?（Phenacoccus?manihoti?Matile-Ferrero）?is?a?dangerous?quarantine?pest?in?the?world，?the?cultivation?and?utilization?of?insect-resistant?cassava?cultivars?can?effectively?block?its?colonization.?Mining?insect-resistant?secondary?metabolites?and?the?regulatory?genes?is?one?of?the?important?strategies?for?insect-resistant?breeding.?Flavonoids?are?unique?secondary?metabolites?for?plants?to?resist?biotic?and?abiotic?stresses，?but?the?function?of?flavonoids?and?the?synthesis?pathway?genes?in?cassava?resistance?to?cassava?mealybug?is?still?unclear.?Based?on?this，?this?study?analyzed?the?expression?levels?of?flavonoid?synthesis?pathway-related?genes?（CHS，?CHI，?FLS，?LAR，?DFR，?F3H，?CCoAOMT，?C4H，?C3'H?and?ANR）?and?flavonoid?content?in?the?leaves?of?insect-resistant?（C1115，?SC9，?Myanmar）?and?insect-susceptible?（KU50，?SC205?and?Bread）?cassava?cultivars?after?being?damaged?by?cassava?mealybug?for?different?time?（0，?1，?4，?8?d）.?Although?the?expression?of?CCoAOMT，?C3'H，?ANR?and?C4H?in?the?leaf?was?up-regulated?after?feeding，?there?was?no?significant?difference?compared?with?that?before?damage，?and?there?was?no?significant?difference?in?the?expression?level?between?the?resistant?and?susceptible?cassava?cultivars.?In?contrast，?the?expression?levels?of?CHS，?CHI，?FLS，?F3H，?DFR?and?LAR?genes?were?significantly?higher?than?those?before?damage，?and?in?the?same?damage?time，?the?expression?levels?of?the?six?genes?in?the?insect-resistant?cassava?cultivars?were?also?significantly?higher?than?those?in?the?insect-susceptible?cassava?cultivars.?Further?correlation?analysis?showed?that?the?expression?levels?of?CHS，?CHI，?FLS，?F3H?and?LAR?genes?were?significantly?positively?correlated?with?the?insect?resistance?of?cassava?cultivars.?In?addition，?the?results?of?total?flavonoid?content?determination?showed?that?the?total?flavonoid?content?increased?significantly?after?one?day?of?damage?compared?with?that?before?damage，?and?the?total?flavonoid?content?of?insect-resistant?cassava?cultivars?was?significantly?higher?than?that?of?the?insect-susceptible?cassava?cultivars?after?4?days?of?damage.?Correlation?analysis?showed?that?the?total?flavonoid?content?was?also?significantly?positively?correlated?with?the?insect?resistance?of?cassava?cultivars.?Therefore，?it?is?speculated?that?the?increase?of?flavonoid?content?and?the?significantly?up-regulation?of?CHS，?CHI，?FLS，?F3H?and?LAR?in?the?insect-resistant?cassava?cultivars?may?be?related?to?the?resistance?to?cassava?mealybug.?This?study?would?provide?an?important?preliminary?basis?for?the?in-depth?analysis?of?the?molecular?mechanism?of?flavonoid?synthesis?genes?regulating?the?insect-resistant?defense?response?of?cassava?to?mealybug，?as?well?as?the?molecular?design?and?breeding?of?cassava?insect-resistant.

Keywords：?insect-resistant?and?insect-susceptible?cassava?cultivars;?Phenacoccus?manihoti;?flavonoid?biosynthetic?pathway;?gene?expression;?resistance?mechanism

DOI：?10.3969/j.issn.1000-2561.2023.12.010

木薯（Manihot?esculenta?Crantz）為大戟科木薯屬灌木或亞灌木作物，在全球一百多個國家包括我國的海南、廣西、廣東、福建、湖南、云南等省（區）大面積種植[1-2]，是重要的糧食作物、能源作物、工業原料和飼料作物[3-4]，在保障世界糧食安全和能源安全中具有不可替代的地位，并作為“先鋒作物”在“一帶一路”倡議中發揮重要作用。

木薯綿粉蚧（Phenacoccus?manihoti?Matile-?Ferrero）是半翅目、粉蚧科、綿粉蚧屬的一種世界危險性檢疫性害蟲[5]，最初在南美洲被發現，2011年該蟲被我國農業農村部列為高風險性檢疫有害生物[6]，目前已在世界木薯主產區普遍嚴重發生[7-8]。該蟲主要通過若蟲和成蟲直接刺吸植株組織汁液為害，造成頂芽、葉片卷曲和束頂，節間減少，發育遲緩，莖桿和根系變弱，常導致減產50%以上，嚴重危害時導致整株死亡和絕收[9-10]。此外，該蟲還會誘發煤煙病加重為害。目前對于該蟲的防治主要以化學防治和生物防治[11]為主，但目前登記用于防治該蟲的藥劑非常少，而生物防治因木薯產地生態環境多樣性低導致天敵種群數量不穩定，生物防治效果差。因此，尋求新的有效綠色防控木薯綿粉蚧的策略和途徑成為目前國內外木薯產業發展中亟待解決的重大課題。

作物抗蟲性是使植物通過內稟的抗蟲機制免于受害、利用內因治蟲所特有的手段，培育和利用抗蟲品種是目前國內外公認的最經濟、簡便、有效的一項帶有方向性的蟲害綠色防控途徑[12-13]，挖掘具有抗蟲作用的次生代謝物質是目前開展抗蟲育種的重要策略。黃酮類化合物是植物體內一類特有的次生代謝物質，在酚類化合物中數量最多，目前已知的黃酮類化合物已超過10?000種[14]。已報道與植物抗蟲相關的多種黃酮類化合物及其調控基因。研究發現無色花青素還原酶（LAR）、花青素還原酶（ANR）基因的過表達可提高單寧含量，使木薯對二斑葉螨具有抗性[15]；在水稻中添加櫻花素作為激發子可以有效提高水稻對稻飛虱的抗性[16]；將查爾酮異構酶（CHI）基因導入抗蟲棉中可增強其對病害和棉鈴蟲的抗性[17]；松茸毒蛾和松毛蟲在為害馬尾松后，馬尾松的類黃酮合成通路基因類黃酮3-羥化酶（F3H）、肉桂酸4-羥化酶（C4H）和類黃酮-3'5'-羥化酶（F3'5'H）顯著上調表達，提高了馬尾松的抗蟲性[18]。然而，迄今為止尚無類黃酮調控基因作用于木薯抗木薯綿粉蚧中的相關研究報道。因此，本研究以抗（C1115、SC9、緬甸）、感（KU50、SC205、面包）木薯品種為參試材料，系統開展木薯綿粉蚧取食前后抗、感木薯品種類黃酮含量及其合成途徑基因表達的差異研究，以期明確類黃酮合成通路在木薯抗木薯綿粉蚧中的重要作用。

1??材料與方法

1.1??材料

1.1.1??供試木薯品種??選用經抗蟲性鑒定且表性穩定的抗蟲木薯品種C1115、SC9、緬甸和感蟲木薯品種KU50、SC205、面包為參試材料[19-21]，上述品種均由中國熱帶農業科學院環境與植物保護研究所提供。種植土比例（土壤∶泥炭土∶珍珠巖=1∶1∶1）。

1.1.2??供試木薯綿粉蚧??木薯綿粉蚧為用感蟲木薯品種KU50在中國熱帶農業科學院環境與植物保護研究所室內長期繼代飼養的實驗種群，飼養溫度設置為（28±2）℃，相對濕度為75%±5%。

1.1.3??類黃酮合成途徑基因??根據已發表的相關文獻以及在KEGG（https：//www.kegg.jp/kegg/）網站上發布的類黃酮合成通路，選取10個研究廣泛、挖掘潛力最大的類黃酮合成基因，用于木薯綿粉蚧取食前后類黃酮合成基因在抗、感蟲木薯品種間表達差異及其與木薯對木薯綿粉蚧的抗蟲性相關性分析。10個類黃酮合成通路基因為：查爾酮合酶（CHS）、CHI、F3H、黃酮醇合成酶（FLS）、二氫黃酮醇還原酶（DFR）、LAR、ANR、咖啡酰輔酶A-O-甲基轉移酶（CCOAOMT）、C4H、香豆酸-3羥基化酶（C3'H）。

1.2??方法

1.2.1??木薯品種對木薯綿粉蚧的抗蟲性鑒定??選擇前期研究中對二斑葉螨和木瓜秀粉蚧分別表現出抗性的木薯品種C1115、SC9、緬甸和感性的木薯品種KU50、SC205、面包[19-21]，同時通過室內鑒定上述品種對木薯綿粉蚧的抗蟲性。根據《木薯種質資源抗蟲性鑒定技術規程》（NY/T?2445—2013），木薯綿粉蚧為害高峰期時，每個品種調查12株木薯，每株木薯選取受蟲害最嚴重的3片葉，即樣本量N=36，根據葉片枯黃和卷曲程度將葉片受木薯綿粉蚧為害后分為0、1、2、3、4級，根據以下公式計算蟲害指數（I）：

式中，I為蟲害指數，S為葉片受害級別，Ns為該受害級別葉片數，N為調查葉片總數。據此，可將田間木薯的抗蟲性水平分為免疫（I=0）、高抗（0

1.2.2??木薯綿粉蚧接種與取樣??選取生長周期約50?d、長勢相近的健康木薯，用細毛筆輕輕蘸取大小相近長期繼代飼養生長周期15?d、二齡左右的雌粉蚧接種于木薯中部葉片的背面，以30頭/葉接種于木薯中部3片葉上，并均勻涂抹羊毛脂于葉柄基部以防止粉蚧遷移。隨后在木薯綿粉蚧為害0、1、4、8?d時對葉片取樣，輕掃下葉片上的粉蚧后迅速將葉片用液氮處理，隨后于–80?℃下冷凍保存，備用。每株木薯取樣1次，每個蟲害時間點設3個生物學重復。

1.2.3??RNA提取及合成cDNA第一條鏈??使用多糖多酚植物總RNA提取試劑盒（TIANGEN，美國）對木薯葉片進行RNA提取，并通過1%瓊脂糖凝膠電泳檢測其完整性和純度，并用超微量紫外分光光度計（Thermo）檢測其濃度。確定RNA可使用后，進一步通過Thermoscript?1st?Strand?cDNA?Synthesis?Kit（Vazyme，中國）進行基因組DNA的清除以及逆轉錄反應，進行cDNA第一條鏈的合成，其中RNA樣品加入1.0?μg，選擇Oligo?（dT）作為逆轉錄引物。

1.2.4??實時熒光定量PCR??以木薯的Meactin基因作為內參基因，相關基因的引物序列信息見表1。通過2×Q3?SYBR?qPCR?Master?Mix試劑盒（TOLOBIO，中國）說明書配置qPCR反應體系，qPCR反應條件為：首先在95?℃下進行60?s預變性，循環反應為：95?℃變性15?s，58?℃下退火30?s，72?℃延伸20?s，使用LightCycler?96儀器（Roche，瑞士）默認程序采集擴增曲線和溶解曲線。將未受木薯綿粉蚧取食的類黃酮合成通路基因的表達量統一設置為1.0，粉蚧為害后的類黃酮合成通路基因的表達量變化情況以為害前的相對倍數表示，以2–ΔΔCt法[22]計算分析，每個處理設3個生物性重復，并進行3個技術性重復。

1.2.5??木薯葉片總黃酮含量的測定??在堿性亞硝酸鹽溶液中的類黃酮和鋁離子在510?nm處可形成有特征吸收峰的紅色絡合物，據此采用微量法測定木薯葉片類黃酮含量，參照植物類黃酮試劑盒（Mlbio，中國上海）說明書（https：//www.?mlbio.cn/）進行測定。

1.3??數據處理

采用Gradphad?Prism?9軟件進行數據匯總和整理，采用完全隨機設計（completely?randomized?design）進行試驗設計，并采用SPSS?26.0軟件對木薯綿粉蚧取食不同木薯品種后基因表達量變化倍數進行方差分析和多重比較。首先對原始數據的正態性和方差齊性進行檢驗，通過峰度和偏度的計算確認數據是否服從正態分布；其次若方差齊性則采用Duncans新復極差法方法進行多重比較，若方差不齊時則對數據進行平方根轉換，并采用Dunnetts?T3方法進行數據間的多

重比較（P=0.05）。通過Spearman相關性分析計算基因表達量或類黃酮含量和蟲害指數之間相關系數（ρ）。

2??結果與分析

2.1??木薯品種對木薯綿粉蚧的抗蟲性鑒定

由表2可知抗蟲木薯品種C1115、SC9、緬甸和感蟲木薯品種KU50、SC205、面包的蟲害指數（I）分別為7.90%、9.60%、11.50%、93.25%、91.80%、95.20%，因此，C1115、SC9、緬甸對木薯綿粉蚧的抗蟲水平均為高抗，而KU50、SC205、面包對木薯面粉蚧的抗蟲水平均為高感。

2.2??木薯綿粉蚧為害抗、感木薯品種前后葉片CHS基因的表達差異分析

在木薯綿粉蚧為害后，木薯類黃酮合成通路CHS基因的表達量在抗、感木薯品種中均有顯著變化，呈先上升后下降的趨勢，其中抗蟲木薯品種中基因表達量的變化更為顯著。在木薯綿粉蚧為害1?d時，C1115、SC9、緬甸葉片中的CHS基因表達量分別上升到為害前的1.81倍、2.12倍、1.66倍；粉蚧為害4?d時，C1115、SC9葉片中的CHS基因表達量進一步上升到為害前的4.29倍、2.51倍，而感蟲木薯品種KU50、面包、SC205中無顯著變化；粉蚧為害8?d時，6個品種的CHS基因表達量均顯著降低（圖1A）。進一步橫向對比同一時間不同木薯品種的CHS基因表達量，發現緬甸在粉蚧為害0、1?d時的表達量遠遠大于其他木薯品種，抗蟲木薯品種在粉蚧為害4?d后顯著高于感蟲品種，為害8?d時抗感差異最顯著（圖1B）。相關性分析表明，粉蚧為害4?d后CHS基因的表達水平和木薯品種的抗蟲性呈顯著正相關（ρ=0.8285，?P=0.0415）。

2.3??木薯綿粉蚧為害抗、感木薯品種前后葉片CHI基因的表達差異分析

在木薯綿粉蚧為害后，CHI基因的表達量在抗蟲木薯品種C1115、SC9、緬甸以及感蟲木薯品種面包葉片中變化顯著，均呈先增后降的趨勢，而在SC205、KU50中基因表達量的變化不顯著或稍有所下降。其中在木薯綿粉蚧為害1?d后，CHI基因表達量在面包、SC9、緬甸中分別上升到為害前的1.74倍、1.98倍、1.61倍；粉蚧為害4?d后，抗蟲品種C1115、緬甸分別上升到為害前的2.22倍、4.41倍；在粉蚧為害8?d時，CHI基因表達量在各個品種中均下降至較低水平（圖2A）。對比同一時間不同木薯品種CHI基因的表達量，發現抗蟲木薯品種在粉蚧為害4?d后開始顯著高于感蟲木薯品種（圖2B）。相關性分析表明，粉蚧為害4?d時CHI基因表達量與木薯品種的抗蟲性呈顯著正相關（ρ=0.9428，?P=0.0048）。

2.4??木薯綿粉蚧為害抗、感木薯品種前后葉片FLS基因的表達差異分析

木薯綿粉蚧為害后，FLS基因的表達量在抗蟲木薯品種C1115、SC9、緬甸以及感蟲品種KU50、SC205中均呈先升后降的趨勢，其中抗蟲品種的表達量變化較大。在粉蚧為害1?d后，KU50、SC205、SC9葉片中的FLS基因表達量分別上升到為害前的2.09倍、2.19倍、4.36倍，隨后在為害4?d時基因表達量下降。粉蚧為害4?d

后，C1115、緬甸分別上升到為害前的2.07倍、3.54倍；面包在第8天上升至為害前的1.66倍（圖3A）。對比同一時間不同木薯品種的FLS基因表達量，發現在木薯綿粉蚧為害4?d后抗蟲木薯品種的FLS基因表達量顯著高于感蟲木薯品種，其中抗蟲品種C1115的表達量最高（圖3B）。相關性分析表明，粉蚧為害4?d后FLS基因的表達量與木薯品種的抗蟲性呈顯著正相關（ρ=0.9428，?P=0.0048）。

2.5??木薯綿粉蚧為害抗、感木薯品種前后葉片LAR基因的表達差異分析

在木薯綿粉蚧為害后，抗、感木薯品種葉片中LAR基因的表達量變化顯著，表現出先升后降的趨勢，且抗蟲木薯品種升高更為顯著。其中在粉蚧為害1?d時，SC9葉片中的LAR基因表達量上升到為害前的1.48倍，在面包中上升到1.34倍，而KU50、SC205無顯著變化；粉蚧為害4?d后，C1115、緬甸上升到為害前的1.79倍、1.55倍（圖4A）。對比同一時間不同木薯品種的LAR基因表達量后發現，粉蚧為害4?d時抗蟲木薯品種顯著高于感蟲木薯品種，且抗蟲木薯品種SC9在粉蚧為害0?d和1?d時也表現出高表達量（圖4B）。相關性分析表明，粉蚧為害4?d后LAR基因的表達量與木薯品種的抗蟲性呈顯著正相關（ρ=0.8285，?P=0.0415）。

2.6??木薯綿粉蚧為害抗、感木薯品種前后葉片DFR基因的表達差異分析

在木薯綿粉蚧為害后，除KU50、SC205葉片中DFR基因的表達量未升高外，C1115、SC9、面包均呈先升后降的趨勢，且抗蟲木薯品種比感蟲品種升高更為顯著。其中SC9、面包在木薯綿粉蚧為害1?d后DFR基因的表達量上升到為害前

的1.83倍、1.29倍，SC9、C1115在粉蚧為害4?d后上升到為害前的1.92倍、3.73倍（圖5A）。進一步對比同一時間不同木薯品種DFR基因的表達量后發現，在粉蚧為害8?d后，抗蟲木薯品種顯著大于感蟲木薯品種（圖5B）。相關性分析表明，粉蚧為害8?d時DFR基因的表達量與木薯品種的抗蟲性呈不顯著正相關（ρ=0.7142，?P=0.1107）。

2.7??木薯綿粉蚧為害抗、感木薯品種前后葉片F3H基因的表達差異分析

木薯綿粉蚧為害后，F3H基因在感蟲木薯品種KU50葉片中的表達量下降，在SC205、面包、C1115、SC9、緬甸中顯著上升后下降。粉蚧為害1?d后，抗蟲木薯品種C1115、SC9、緬甸葉片中的F3H基因表達量分別上升到為害前的1.55倍、2.44倍、3.64倍，感蟲品種SC205、面包上升到為害前的1.30倍、1.70倍。粉蚧為害4?d后，C1115

葉片的F3H基因表達量進一步上升到為害前的2.24倍，而其他木薯品種中呈下降趨勢。粉蚧為害8?d時，各木薯品種的F3H基因表達量進一步降低（圖6A）。進一步對比相同時間不同木薯品種葉片中F3H基因的表達量后發現，抗蟲木薯品種在粉蚧為害4?d后顯著高于感蟲品種（圖6B）。相關性分析表明，粉蚧為害4?d后F3H基因的表達量與木薯品種的抗蟲性呈顯著正相關（ρ=0.8285，?P=0.0415）。

2.8??木薯綿粉蚧為害抗、感木薯品種前后葉片CCoAOMT基因的表達差異分析

木薯綿粉蚧為害后，KU50、SC205、緬甸葉片中的CCoAOMT基因表達量逐漸升高，而SC9無顯著變化，其他品種均呈先升后降的趨勢。粉蚧為害1?d后，木薯品種SC205、面包、C1115葉片中的CCoAOMT基因表達量上升至為害前的2.51倍、1.91倍、1.92倍；粉蚧為害4?d后，感

蟲品種KU50上升至為害前的4.40倍，抗蟲品種C1115、緬甸分別上升至為害前的9.06倍、2.23倍；粉蚧為害8?d后緬甸中的表達量上升至為害前的3.59倍（圖7A）。進一步對比相同時間不同木薯品種葉片中CCoAOMT基因的表達量后發現，在粉蚧為害前與為害1?d時，SC9與面包有較高的表達量，總體上抗蟲品種中的CCoAOMT基因表達量并不一定高于感蟲木薯品種（圖7B）。

2.9??木薯綿粉蚧為害抗、感木薯品種前后葉片C4H基因的表達差異分析

木薯綿粉蚧為害后，C4H基因的表達量在抗蟲木薯品種C1115、SC9葉片中呈先升后降的趨勢，在感蟲木薯品種KU50中表現為先降后升的趨勢，其他木薯品種中無顯著變化或稍有下降。其中粉蚧為害1?d后，抗蟲品種C1115、SC9葉片中C4H基因的表達量上升到為害前的2.16倍、2.72倍；為害4?d后，所有木薯品種中的表達量均呈下降趨勢或無顯著變化；為害8?d后，KU50上升至為害前的1.93倍（圖8A）。對比相同時間不同木薯品種葉片中C4H基因的表達量后發現，SC9中的表達量在粉蚧為害前以及為害4?d時遠大于其他品種，但總體抗、感木薯品種無顯著差異，抗蟲木薯品種并不一定大于感蟲品種（圖8B）。

2.10??木薯綿粉蚧為害抗、感木薯品種前后葉片C3'H基因的表達差異分析

木薯綿粉蚧為害后，KU50、SC205葉片中的C3'H基因表達量呈下降趨勢，C1115、SC9、緬甸、面包中均呈先升后降的趨勢。在粉蚧為害1?d后，C3'H基因在抗蟲品種C1115、SC9、緬甸中分別上升到為害前的1.69倍、2.43倍、1.88倍，在感蟲木薯品種面包中上升至1.68倍，而KU50、SC205中無顯著變化。粉蚧為害4?d時，C3'H基因在抗蟲品種C1115、SC9中分別下降到為害前的27.5%、44.9%，在感蟲品種KU50、SC205中下降到為害前的43%、26.7%（圖9A）。對比相同時間不同木薯品種葉片中C3'H基因的表達量后發現，在粉蚧為害前，抗蟲木薯品種葉片中的C3'H基因表達量大于感蟲木薯品種，粉蚧為害1?d時C1115、SC9的表達量遠遠大于其他品種，而在粉蚧為害4?d和8?d時，抗、感木薯品種之間的表達差異不大或無顯著差異（圖9B）。

2.11??木薯綿粉蚧為害抗、感木薯品種前后葉片ANR基因的表達差異分析

木薯綿粉蚧為害后，抗蟲木薯品種SC9葉片

中的ANR基因表達量呈先升后降再升的趨勢，緬甸中的表達量隨粉蚧為害時間的延長逐漸上升，C1115中無顯著變化；感蟲木薯品種KU50、SC205、面包中均無顯著變化或稍有下降。其中SC9葉片中的ANR基因表達量在粉蚧為害1?d時上升到為害前的1.87倍，4?d時下降到為害前的水平，8?d后升高到為害前的2.12倍；緬甸在為害1?d時上升至為害前的1.53倍，為害8?d后上升至為害前的1.86倍（圖10A）。橫向對比相同時間不同木薯品種葉片中ANR基因的表達量后發現，抗蟲木薯品種SC9在粉蚧為害1、8?d時的表達量高于其他品種，而其他品種葉片中的表達量在各個時間均無顯著差異（圖10B）。

2.12??木薯綿粉蚧為害抗、感木薯品種前后葉片總黃酮含量的差異分析

與木薯綿粉蚧為害前相比，為害后抗、感木薯品種葉片的總黃酮含量變化顯著，為害1?d后，抗、感木薯品種總黃酮含量均顯著上升；而在為害4?d后，抗蟲木薯品種C1115、SC9、緬甸的總黃酮含量仍然維持在較高水平，而感蟲品種KU50、面包、SC205中則呈下降趨勢；粉蚧為害8?d后，6個品種的總黃酮含量無顯著變化或稍有下降，但抗蟲木薯品種中的含量顯著高于感蟲品種（表3）。相關性分析顯示，粉蚧為害4?d時，總黃酮含量與木薯品種的抗蟲性呈顯著正相關（ρ=0.9428，?P=0.0048）。

3??討論

類黃酮通路合成基因在植物應對生物和非生物脅迫過程中發揮重要作用。研究發現，類黃酮通路糖苷轉移酶基因（UGT79B2/B3）的過量表達可以顯著增強擬南芥對低溫、干旱和鹽脅迫的耐受性[23]；大豆中突變黃烷酮-3-羥化酶（F3H）和黃酮合酶（FNS）基因，防止其與異黃酮合酶（IFS）基因競爭底物，從而提高異黃酮的含量以及對大豆花葉病毒（SMV）的抗性[24]；異黃酮合成酶基因的表達量在大豆受高濃度鎘處理時顯著提高，表明其與耐受鎘脅迫密切相關[25]，MURATA等[26]研究發現，柚皮素7-O-甲基轉移酶（NOMT）基因在抗真菌品種日本晴和敏感品種秈稻中的表達量存在顯著差異，從而導致抗真菌毒素物質累積及抗菌水平的顯著差異。有研究發現，在大豆胞囊線蟲的脅迫下，抗病品種灰皮支黑豆和感病品種Williams?82的CHS基因表達量存在顯著差異，抗病品種顯著上調表達，初步認為CHS基因在抗大豆胞囊線蟲的抗性機制中發揮重要作用[27]。本研究表明，木薯綿粉蚧為害4?d后，抗蟲木薯品種類黃酮通路合成基因（CHS、CHI、FLS、DFR、LAR、F3H）的表達量顯著高于感蟲木薯品種，其中CHS、CHI、FLS、F3H、LAR與品種抗蟲性呈顯著正相關，可以推測上述類黃酮通路合成基因可能在木薯抗木薯綿粉蚧中發揮作用，這與上述研究結果具有相似性。研究發現類黃酮合成通路基因的表達水平會影響類黃酮含量累積進而影響植物的抗蟲性。如LAR和ANR是調控黃酮醇與

花青素合成的關鍵基因[28]，過表達LAR、ANR基因提高了原花青素類物質縮合單寧的含量，從而提高了木薯對二斑葉螨的抗性[15]。類黃酮通路基因的表達還存在級聯效應，如JAN等[29]通過干預F3H基因可以影響下游基因FLS、DFR等的表達，進而影響花青素、山奈酚、槲皮素等多種類黃酮物質含量，并改變了水稻的抗蟲性。此外，對轉錄因子的表達進行調控，同樣可以影響合成通路上基因的表達及類黃酮水平。如在擬南芥中過表達MYB75轉錄因子，將同時上調CHS、DFR、FLS等基因的表達，提高花青素和黃酮醇的水平，最終提高了植物對歐洲粉蝶幼蟲的抗性[30]。研究表明，植物類黃酮含量與其抗蟲性密切相關。在番茄防御煙粉虱的研究中發現，類黃酮高含量品系比低含量品系更能降低煙粉虱對其韌皮部的攝食效率和番茄黃葉卷曲病毒（TYLCV）的傳播[31]，此外，黃酮類化合物還對害蟲具有體外活性，如異黃酮類物質高麗槐素、牛蒿素等在濃度為100?mmol/L時可以有效阻止棉鈴蟲幼蟲對鷹嘴豆的取食，而加入綠原酸又可進一步增強其拒食活性[32]。本研究發現，與木薯綿粉蚧為害前相比，抗蟲木薯品種中類黃酮含量以及類黃酮合成通路基因CHS、CHI、FLS、LAR、DFR和F3H的表達量均顯著升高，并顯著高于感蟲木薯品種的含量和表達水平，推測抗蟲木薯品種被粉蚧為害后，誘導了上述黃酮合成基因的上調表達，從而使具有抗蟲功能的黃酮類物質含量相較于感蟲品種顯著提高，進而增強木薯抗蟲水平。

抗蟲木薯品種受木薯綿粉蚧為害后，葉片中總黃酮含量顯著提高，類黃酮合成通路基因CHS、CHI、FLS、LAR、DFR、F3H顯著上調表達，其中CHS、CHI、FLS、F3H、LAR基因的表達量與木薯品種抗蟲性呈顯著正相關，表明上述5個基因可能參與類黃酮合成途徑介導的木薯抗木薯綿粉蚧的分子機制。本研究分析了木薯綿粉蚧為害后抗、感木薯品種葉片中類黃酮合成通路基因表達量以及總黃酮含量的差異，后續研究需開展相應合成酶活性的測定，明確基因-蛋白-代謝物的變化趨勢是否具有一致性，并進一步通過體外抗蟲活性試驗明確具體發揮抗蟲作用的類黃酮物質，從而進一步揭示特定的類黃酮合成基因介導的木薯抗蟲功能。

參考文獻

[1]?ROY?S，?MARIE?A.?Food?and?Agriculture?Organization?of?the?United?Nations[J].?Australian?Veterinary?Journal，?2006，?35（3）：?105-110.

[2]?歐文軍，?羅秀芹，?安飛飛，?李開綿.?氣候變化與我國木薯北移的可能性分析[J].?中國熱帶農業，?2014，?59（4）：?4-8.OU?W?J，?LUO?X?Q，?AN?F?F，?LI?K?M.?Analysis?of?climate?change?and?the?possibility?of?cassava?moving?northward?in?China[J].?China?Tropical?Agriculture，?2014，?59（4）：?4-8.?（in?Chinese）

• EL-SHARKAWY?M?A.?Cassava?biology?and?physiology?cassava：?a?crop?for?sustainable?agriculture?and?food?security?in?developing?countries[J].?Plant?Molecular?Biology，?2004，?56：?481-501.
• PARMAR?A，?STURM?B，?HENSEL?O.?Crops?that?feed?the?world：?production?and?improvement?of?cassava?for?food，?feed，?and?industrial?uses[J].?Food?Security，?2017，?9（5）：?907-927.
• 盧輝，?盧芙萍，?梁曉，?陳青.?基于生態位模型的木薯綿粉蚧潛在適生性分析研究[C]//中國熱帶作物學會第九次全國會員代表大會暨2015年學術年會論文摘要集，?2015：?69.LU?H，?LU?F?P，?LIANG?X，?CHEN?Q.?Potential?suitability?analysis?of?mealybug?in?cassava?based?on?niche?model[C].?Summary?of?Papers?from?the?9th?National?Congress?of?the?Chinese?Tropical?Crop?Society?and?the?2015?Academic?Annual?Conference，?2015：?69.?（in?Chinese）
• 冬曉.?我國將木薯綿粉蚧和異株莧亞屬雜草列為進境檢疫性有害生物[J].?農業技術與裝備，?2011（14）：?77.DONG?X.?Phenacoccus?manihoti?and?Amaranthus?weeds?are?listed?as?quarantine?pests?in?China[J].?Agricultural?Technology?&?Equipment，?2011（14）：?77.?（in?Chinese）

[7]?BELLOTTI?A?C，?BRAUN?A?R，?ARIAS?B，?CASTILLO?J?A，?GUERRERO?J?M.?Origin?and?management?of?Neotropical?cassava?arthropod?pests[J].?Louvain?Medical，?1994，?2（4）：?S41-42.

[8]?PARSA?S，?KONDO?T，?WINOTAI?A.?The?cassava?mealybug?（Phenacoccus?manihoti）?in?Asia：?first?records，?potential?distribution，?and?an?identification?key[J].?Public?Library?of?Science，?2012，?7（10）：?343.

[9]?NWANZE?F?K.?Relationships?between?cassava?root?yields?and?crop?infestations?by?the?mealybug，?Phenacoccus?manihoti[J].?Tropical?Pest?Management，?1982，?28（1）：?27-32.

[10]?MUNIAPPAN?R，?SHEPARD?B?M，?WATSON?G?W，?CARNER?G?R，?RAUF?A，?SARTIAMI?D，?HIDAYAT?P，?AFUN?J?V?K，?GOERGEN?G，?RAHMAN?A?K?M，?ZIAUR.?New?records?of?invasive?insects?（Hemiptera：?Sternorrhyncha）?in?Southeast?Asia?and?West?Africa[J].?Journal?of?Agricultural?&?Urban?Entomology，?2009，?26（4）：?167-174.

[11]?LOHR?B，?VARELA?AM，?SANTOS?B.?Exploration?for?natural?enemies?of?the?cassava?mealybug，?Phenacoccus?manihoti?（Homoptera：?Pseudococcidae），?in?South?America?for?the?biological?control?of?this?introduced?pest?in?Africa[J].?Bulletin?of?Entomological?Research，?1990，?80（4）：?417-425.

[12]?WAR?A?R，?PAULRAJ?M?G，?AHMAD?T，?BUHROO?A?A，?HUSSAIN?B，?IGNACIMUTHU?S，?SHARMA?H?C.?Mechanisms?of?plant?defense?against?insect?herbivores[J].?Plant?Signaling?Behavior，?2012，?7（10）：?1306-1320.

• SMITH?C?M.?Conventional?breeding?of?insect-resistant?crop?plants：?still?the?best?way?to?feed?the?world?population[J].?Current?Opinion?in?Insect?Science，?2021，?45：?7-13.
• 王軍妮，?黃艷紅，?牟志美，?王彥文.?植物次生代謝物黃酮類化合物的研究進展[J].?蠶業科學，?2007，?33（3）：?499-505.WANG?J?N，?HUANG?Y?H，?MU?Z?M，?WANG?Y?W.?Research?progress?on?flavonoid?of?the?plant?secondary?metabolites[J].?Science?of?Sericulture，?2007，?33（3）：?499-505.?（in?Chinese）
• CHEN?Q，?LIANG?X，?WU?C，?LIU?Y，?LIU?X?Q，?ZHAO?H?P，?LI?K?M，?CHEN?S?B，?WANG?H?Y，?HAN?Z?L，?WU?M?F，?YAO?X?W，?SHUI?J，?QIAO?Y，?ZHAN?X，?ZHANG?Y.?Overexpression?of?leucoanthocyanidin?reductase?or?anthocyanidin?reductase?elevates?tannins?content?and?confers?cassava?resistance?to?two-spotted?spider?mite[J/OL].?Frontiers?in?Plant?Science，?2022，?13.?[2023-08-14].?https：//doi.org/10.3389/fpls.?2022.994866.
• 顏健，?楊忠艷，?黎平，?白根祥.?櫻花素作為激發子在誘導農作物抗蟲中的應用：?CN202110313675.1[P].?2023-08-14.YAN?J，?YANG?Z?Y，?LI?P，?BAI?G?X.?The?application?of?sakuranin?as?an?elicitor?in?inducing?crop?insect?resistance：?CN202110313675.1[P].?2023-08-14.?（in?Chinese）

[17]?簡桂良，?盧美光，?王志興，?賈士榮，?鄒亞飛，?程紅梅，?孫文姬.?Chi基因導入對轉Bt棉花抗病蟲性的影響[J].?農業生物技術學報，?2004，?12（4）：?478-479.JIAN?G?L，?LU?M?G，?WANG?Z?X，?JIA?S?R，?ZOU?Y?F，?CHENG?H?M，?SUN?W?J.?Influence?of?Chi?gene?on?cotton?insect?and?diseace?resistance[J].?Journal?of?Agricultural?Biotechnology，?2004，?12（4）：?478-479.?（in?Chinese）

[18]?YANG?Z?Q，?CHEN?H，?TAN?J?H，?Xu?H?L，?Jia?J，?Feng?Y?H.?Cloning?of?three?genes?involved?in?the?flavonoid?metabolic?pathway?and?their?expression?during?insect?resistance?in?Pinus?massoniana?Lamb[J].?Genetics?and?Molecular?Research，?2016，?15（4）：?15049332.

[19]?YAO?X?W，?LIANG?X，?CHEN?Q，?LIU?Y，?WU?C?L，?WU?M?F，?SHUI?J，?QIAO?Y，?ZHANG?Y，?GENG?Y.?MePAL6?regulates?lignin?accumulation?to?shape?cassava?resistance?against?two-spotted?spider?mite[J].?Frontiers?in?Plant?Science，?2023，13：?1067695.

[20]?陳青，?梁曉，?伍春玲，?陳謙，?王亞茹.?不同木薯品種對木瓜秀粉蚧生殖和發育的影響[J].?基因組學與應用生物學，?2020，?39（5）：?2177-2182.?CHEN?Q，?LIANG?X，?WU?C?L，?CHEN?Q，?WANG?Y?R.?Effects?of?different?cassava?cultivars?on?reproduction?and?development?of?mealybug[J].?Genomics?and?Applied?Biology，?2020，?39（5）：?2177-2182.?（in?Chinese）

[21]?陳謙.?木瓜秀粉蚧對不同木薯種質的寄主選擇性研究[D].?海口：?海南大學，?2022.CHEN?Q.?Study?on?the?host?selectivity?of?Paracoccus?marginatus?to?different?cassava?germplasms[D].?Haikou：?Hainan?University，?2022.?（in?Chinese）

[22]?LIVAK?K?J，?SCHMITTGEN?T?D.?Analysis?of?relative?gene?expression?data?using?real-time?quantitative?PCR?and?the?2?ΔΔCT?method[J].?Methods，?2001，?25（4）：?402-408.

[23]?LI?P，?LI?Y?J，?ZHANG?F?J，?ZHANG?G?Z，?JIANG?X?Y，?YU?H?M，?HOU?B?K.?The?Arabidopsis?UDP-glycosyltransferases?UGT79B2?and?UGT79B3，?contribute?to?cold，?salt?and?drought?stress?tolerance?via?modulating?anthocyanin?accumulation[J].?Plant?Journal，?2017，?89（1）：?85-103.

[24]?ZHANG?P?P，?DU?H?Y，?WANG?J，?YU?D?Y.?Multiplex?CRISPR/Cas9-mediated?metabolic?engineering?increases?soya?bean?isoflavone?content?and?resistance?to?Soyabean?mosaic?virus[J].?Plant?Biotechnology?Journal，?2020，?18（6）：?1384-1395.

[25]?劉曉慶，?張大勇，?徐照龍，?許玲，?黃益洪，?何曉蘭，?馬鴻翔，?易金鑫.?Cd脅迫對大豆幼苗異黃酮合成關鍵酶基因表達的影響[J].?江蘇農業學報，?2013，?29（5）：?974-978.LIU?X?Q，?ZHANG?D?Y，?XU?Z?L，?XU?L，?HUANG?Y?H，?HE?X?L，?MA?H?X，?YI?J?X.?Gene?expression?of?isoflavonoid?synthesis?related?enzymes?in?soybean?[J].?Jiangsu?Journal?of?Agricultural?Sciences，?2013，?29（5）：?974-978.?（in?Chinese）

• MURATA?K，?KITANO?T，?YOSHIMOTO?R，?TAKATA?R，?UBE?N，?UENO?K，?UENO?M，?YABUTA?Y，?TERAISHI?M，?HOLLAND?C，?JANDER?G，?OKUMOTO?Y，?MORI?N，?ISHIHARA?A.?Natural?variation?in?the?expression?and?catalytic?activity?of?a?naringenin?7-O-methyltransferase?influences?antifungal?defenses?in?diverse?rice?varieties[J].?Plant?Journal，?2020，?101（5）：?1103-1117.
• 王爽，?張順斌，?王旭東，?于佰雙，?陳井生，?王惠，?段玉璽.?大豆胞囊線蟲脅迫下GmCHS基因表達及異黃酮含量變化分析[J].?大豆科學，?2021，?40（1）：?21-27.WANG?S，?ZHANG?S?B，?WANG?X?D，?YU?B?S，?CHEN?J?S，?WANG?H，?DUAN?Y?X.?Expression?analysis?of?GmCHS?gene?and?isoflavone?content?changes?under?soybean?cyst?nematode?stress[J].?Soybean?Science，?2021，?40（1）：?21-27.?（in?Chinese）
• CAO?Y，?XING?M，?XU?C，?LI?X.?Biosynthesis?of?flavonol?and?its?regulation?in?plants[J].?Acta?Horticulturae?Sinica，?2018，?45（1）：?177-192.
• JAN?R，?KHAN?M?A，?ASAF?S，?LEE?I?J，?KIM?K?M.?Overexpression?of?OsF3H?modulates?WBPH?stress?by?alteration?of?phenylpropanoid?pathway?at?a?transcriptomic?and?metabolomic?level?in?Oryza?sativa[J].?Scientific?Reports，?2020，?10（1）：?14685.
• ONKOKESUNG?N，?REICHELT?M，?VAN?D?A，?SCHUURINK?R?C，?VAN?L?J?J?A，?MARCEL?D.?Modulation?of?flavonoid?metabolites?in?Arabidopsis?thaliana?through?overexpression?of?the?MYB75?transcription?factor：?role?of?kaempferol-3，7-dirhamnoside?in?resistance?to?the?specialist?insect?herbivore?Pieris?brassicae[J].?Journal?of?Experimental?Botany，?2014，?65（8）：?2203-2217.
• YAO?Q，?PENG?Z，?TONG?H，?YANG?F?B，?XING?G?S，?WANG?L?J，?ZHENG?J?J，?ZHANG?Y?J，?SU?Q.?Tomato?plant?flavonoids?increase?whitefly?resistance?and?reduce?spread?of?Tomato?yellow?leaf?curl?virus[J].?Journal?of?Economic?Entomology，?2019，?112（6）：?2790-2796.
• SIMMONDS?M，?STEVENSON?P?C.?Effects?of?isoflavonoids?from?cicer?on?larvae?of?Heliocoverpa?armigera[J].?Journal?of?Chemical?Ecology，?2001，?27（5）：?965-77.