薇甘菊頸盲蝽細胞色素P450 PmCYP4基因克隆及表達分析

張傳光， 澤桑梓， 朱家穎， 季 梅， 劉 凌， 張乃明

1云南農業大學植物保護學院，云南 昆明 650201； 2云南省林業科學院，云南 昆明650204； 3云南省林業有害生物防治檢疫局，云南 昆明 650051； 4西南林業大學，云南 昆明 650224

薇甘菊MikaniamicranthaH. B. K.是世界上最具侵略性的雜草之一(范志偉等,2016; Mingetal.,2017)，也是唯一一種全國林業檢疫有害植物(徐小偉等,2014)。物種在入侵過程中與入侵生境中其他物種相互競爭、協同適應，如入侵我國云南的薇甘菊揮發油內普遍存在β-cadinene、allo-aromadendrene、β-caryophyllene及5-(1,1-dimethylethyl)-2,3-1H-Inden-1-one 等多碳化合物(季梅等,2012; 孫盟等,2013)，而生長在原產地秘魯的薇甘菊檢測不到這些物質，說明入侵地薇甘菊的次生代謝物種類發生了改變，且含量也明顯增加(Nietal.,2007)。此外，薇甘菊在入侵的過程中其體內的化學防御素——縮合單寧也顯著高于本地物種(倪廣艷等,2014)。

Yiran & Liu (2017)首次報道了一個專食薇甘菊的新種——薇甘菊頸盲蝽PachypeltismicranthusMu et Liu，是控制云南省瑞麗市薇甘菊的一種本土天敵昆蟲。研究發現，被薇甘菊頸盲蝽取食后，薇甘菊葉片中的防御性酶POD和PAL的活性下降，PPO活性提高(季梅等,2014； 李勝等,2018)，說明在被頸盲蝽取食后薇甘菊很快做出了防御反應。植物次生物質能對昆蟲產生不利影響，甚至產生毒殺作用，而昆蟲則通過對植物次生物質忌避取食、解毒代謝等多種機制，對寄主植物產生適應性(陳澄宇等,2015)。其中，細胞色素P450 (cytochrome P450, CYP)作為最古老和最龐大的超基因家族，其介導的多功能氧化酶(microsomal mixed function oxidases，MFO)是昆蟲中涉及抗性的3大主要代謝解毒酶類之一，對多樣的外源化合物和內源化合物的氧化代謝起作用(Hrycay & Bandiera,2009; Omura,1999)，特別是CYP4家族基因在昆蟲對植物次生物質的解毒代謝及與寄主植物相互作用中發揮了重要作用(吳益東等,1997; Dingetal.,2013; Edietal.,2014; Feyereisen,1999)。

Feyereisenetal. (1989)克隆了第一個昆蟲CYP基因，此后，昆蟲CYPs逐漸成為研究熱點(艾均文等,2015)。對于微甘菊頸盲蝽而言，如何應對宿主植物的化學防御，其體內存在何種代謝解毒酶，相關基因在不同部位的表達量如何，都未有相關的報道。本文旨在克隆薇甘菊頸盲蝽CYP4基因，研究其在不同部位的表達情況，為深入研究其在薇甘菊頸盲蝽與薇甘菊互作中的作用奠定基礎。

1 材料與方法

1.1 材料

1.1.1 供試昆蟲 薇甘菊頸盲蝽采自云南省瑞麗市。

1.1.2 實驗試劑及主要儀器 Trizol、SMARTTMRACE cDNA amplification kit (Clontech)、DNA Marker DL2000均購自寶生物工程(大連)有限公司，PrimeScripttmRT reagent Kit with gDNA Eraser購自Takara公司。

高速冷凍離心機(5804R型，德國Eppendorf公司)，PCR儀(Veriti96型，美國Applied Biosystems公司)，熒光定量PCR儀(Rotor Gene-Q型，德國QIAGEN公司)，凝膠成像儀(GBOX iChemi型，英國Syngene公司)，超微量紫外分光光度計(NanoDrop 2000，美國Thermo公司)。

1.1.3 培養基 LB培養基：酵母提取物、胰蛋白胨、NaCl分別為5、10、10 g·L-1，用NaOH調節pH為7；加入15 g·L-1瓊脂為LB固體培養基。

1.2 實驗方法

1.2.1 樣品采集方法 將捕捉到的雌、雄活體切分為觸角、殘體、翅膀、足4個部分，迅速保存至液氮中帶回實驗室，每部分3個重復，每個重復100頭。

1.2.2 總RNA的提取及CYP基因的克隆 參照Trizol試劑說明書，以樣品研磨提取總RNA，提取到的總RNA，采用紫外分光光度計檢測其質量，儲存在-80 ℃冰箱備用。

以提取的總RNA為模板，用SMARTTMRACE cDNA amplification kit (Clontech)反轉錄試劑盒合成cDNA第一鏈。根據實驗室前期對薇甘菊頸盲蝽cDNA文庫測序獲得的CYP片段序列，設計5′RACE和3′RACE特異性引物(5′-GGCAAAGGTCTGCTGACAAGTA-3′和5′-TGAATGGTACGTAGGCAA

ACGG-3′)，以反轉錄合成的cDNA為模板，參照RACE試劑盒說明書，PCR擴增該基因的3′端和5′端序列。PCR反應條件為：94 ℃預變性2 min；94 ℃變性30 s，55 ℃退火30 s，72 ℃延伸40 s，35個循環；72 ℃延伸10 min。用1%的瓊脂糖凝膠電泳對PCR產物進行檢測。利用凝膠回收試劑盒對PCR產物進行回收純化，TA克隆接入pGEM?-T-easy載體(Promega)，藍白斑篩選，挑取陽性克隆送至昆明碩擎生物科技有限公司進行測序。

1.2.3 序列分析 用ORFfinder (https:∥www.ncbi.nlm.nih.gov/orffinder/)將CYP的核苷酸序列翻譯成氨基酸序列；用ProtParam (https:∥web.expasy.org/protparam/)對蛋白進行基本進化性質分析；用SignalP 4.1 Server (http:∥www.cbs.dtu.dk/services/SignalP-4.1/)對信號肽進行預測；用ClustalX 1.83進行多序列比對分析進化樹(Thompsonetal.,1997)，MEGA 7.0軟件構建NJ (neighbor-joining)分子系統進化樹(Kumaretal.,2016)。

1.2.4 熒光定量PCR 用Trizol試劑提取各樣品總RNA，將各樣品總RNA定量至1 μg，用PrimeScripttmRT reagent Kit with gDNA Eraser試劑盒對提取的各樣品總RNA先進行基因組DNA的去除，定量后再反轉錄成cDNA模板。根據實驗前期克隆獲得的CYP cDNA序列，設計特異性引物(5′-CGTTTTTCTCCGAGCGTTC-3′和5′-TTGAACTTCGA

CGGGTTGG-3′)，對該基因在薇甘菊頸盲蝽不同性別、不同部位中的表達量進行熒光定量PCR分析。以18S RNA為內參基因(5′-TTTCAAATGTCTGCCTTATC-3′和5′-TGTGGTAGCCGTTTCTCA-3′)，每個樣品重復3次。PCR反應條件為：95 ℃預變性2 min；95 ℃變性10 s，55 ℃退火20 s，39個循環。熒光定量PCR結果采用E-??CT法進行計算，根據該方法把表達量最低的樣品表達量值定義為1(Livak & Schmittgen,2001)。最后，用Excel對數據進行整理、IBM SPSS Statistics進行方差分析和S-N-K多重比較、Origin Pro 2016進行圖形繪制。

2 結果與分析

2.1 PmCYP4生物信息學分析

2.1.1PmCYP4編碼氨基酸理化性質PmCYP4 cDNA全長1713 bp，ORF為1500 bp，編碼500個氨基酸，5′端非編碼區27 bp，3′端非編碼區183 bp，無信號肽(圖1)。編碼氨基酸的基本理化性質見表1。

圖1 PmCYP4基因的開放閱讀框和推測的氨基酸序列Fig.1 ORF and the deduced amino acid sequences of PmCYP4

指標Index預測結果 Predicted outcome分子式 Molecular formulaC2611H4061N673O739S23分子質量 Molecular weight/(ku)57.44氨基酸數 Number of amino acids500等電點 Isoelectric point6.93負電荷殘基 Asp+Glu63正電荷殘基 Arg+Lys62平均親水系數Average hydrophilicity-0.279不穩定系數 Instability index28.84脂肪系數 Aliphatic index84.02

2.1.2PmCYP4基因序列分析PmCYP4編碼的氨基酸序列中的2個保守區域見圖1和圖2(下劃線部分)：(1)CYP4家族成員保守區EVDTFMFEGHDTT(周夏,2015);(2)FXXGXXXCXG/A區域，為紅素結合環，C是保守的半胱氨酸，X代表任意密碼子，其為血紅素鐵第5個配體，也是在450 nm處與CO結合產生吸收峰特征性光譜的原因(邱星輝和冷欣夫,1998; 鄭明奇等,2001)。

不同物種、不同來源的CYP基因均有高度的序列同源性，序列比對發現，PmCYP4與地中海實蠅Ceratitiscapitata、家蠅Muscadomestica、辣椒實蠅Bactroceralatifrons、銅綠蠅LuciliacuprinaCYP4的氨基酸同源性分別為41%、42%、43%、45%(圖2)。系統進化樹顯示，PmCYP4與溫帶臭蟲CimexlectulariusCPY4基因編碼的蛋白同源性處于同一分支上，進化關系最近(圖3) 。

圖2 PmCYP4與其他昆蟲CYP多序列比對Fig.2 Multiple alignment of PmCYP4 and other insects CRY：砂白蟻,PNF27066.1;ZOO：內華達古白蟻,KDR12230.1;NIL：褐飛虱,AIW79997.1;CAM：佛羅里達弓背蟻，XP_025266601.1;PAC：薇甘菊頸盲蝽。下劃線部分為PmCYP4編碼的氨基酸序列中的2個保守區域。黑色和灰色陰影分別表示保守性極高和較高的殘基位點。 CRY： Cryptotermes secundus, PNF27066.1; ZOO： Zootermopsis nevadensis, KDR12230.1; NIL： Nilaparvata lugens, AIW79997.1;CAM： Camponotus floridanu, XP_025266601.1; PAC： Pachypeltis micranthus. The underlined parts were two conserved domains of the amino acid sequences of PmCYP4. Residues identical or similar are highlighted in black and grey respectively.

圖3 PmCYP4基因系統發育樹Fig.3 Phylogenetic tree of PmCYP4CER:地中海實蠅,XP_004521346.1;BAC:辣椒實蠅,XP 018799661.1;MUS:家蠅,XP_005186030.1;LUC:銅綠蠅, XP_023297259.1;FOP:阿里山潛蠅繭蜂,XP_011300548.1;NEO:紅頭松葉蜂,XP_015512314.1;ZOO:內華達古白蟻,XP_021934792.1;CRY:砂白蟻,XP_023705844.1;ONT:嗡蜣螂,XP_022906724.1;HAL:茶翅蝽,XP_014284338.1;CIM:溫帶臭蟲,XP_014249034.1。CER: Ceratitis capitata, XP_004521346.1; BAC: Bactrocera latifrons, XP 018799661.1; MUS: Musca domestica, XP_005186030.1; LUC: Lucilia cuprina, XP_023297259.1; FOP: Fopius arisanus, XP_011300548.1; NEO: Neodiprion lecontei, XP_015512314.1; ZOO: Zootermopsis nevadensis, XP_021934792.1; CRY: Cryptotermes secundus, XP_023705844.1; ONT: Onthophagus taurus, XP_022906724.1;HAL: Halyomorpha halys, XP_014284338.1; CIM: Cimex lectularius, XP_014249034.1.

2.2 mCYP4基因熒光定量PCR分析

PmCYP4在薇甘菊頸盲蝽雌、雄蟲中的表達量都是足中最高，顯著高于翅膀、觸角、殘體。雌蟲殘體中表達量最低，觸角、翅膀、足中的表達量分別是其1.15、1.88和11.90倍，殘體、觸角和翅膀中的表達量無顯著差異；雄蟲觸角表達量最低，殘體、翅膀、足中的表達分別是其1.02、6.07和16.32倍，且翅膀中的表達量顯著高于殘體和觸角，而殘體和觸角中的表達量無顯著差異；雌蟲翅膀、殘體及觸角與雄蟲相應組織中的表達量無顯著差異，但雄蟲翅膀中的表達量顯著高于雌蟲翅膀，是其2.37倍(圖4)。

圖4 PmCYP4基因在薇甘菊頸盲蝽不同部位中的表達量Fig.4 Relative expression levels of PmCYP4 gene in different parts of P. micranthusFA:雌蟲觸角；FB:雌蟲殘體；FL:雌蟲足；FW:雌蟲翅膀；MA:雄蟲觸角；MB:雄蟲殘體； ML:雄蟲足；MW:雄蟲翅膀。不同字母表示在 5%水平上差異顯著。FA: Antenna of female; FB: Residual body of female; FL: Leg of female; FW: Wing of female; MA: Antenna of male; MB: Residual body of male; ML: Leg of male; MW: Wing of male. Different letters mean significant differences at 5% level．

3 討論

CYP基因家族包括36個基因族(李顯春等,1999; 張宇宏等,2015)，昆蟲中已鑒定的CYP分屬于CYP4、CYP6、CYP9、CYP12、CYP18、CYP28等27個家族，除CYP4家族外全是昆蟲特異的家族(Feyereisen,1999)。通過序列比對，此次克隆到的薇甘菊頸盲蝽PmCYP4基因與CYP4家族基因編碼的蛋白同源性最高,且存在2個保守區域。Blast結果表明，薇甘菊頸盲蝽CYP基因與其他昆蟲的CYP4家族的同源性都大于40%，所以PmCYP4應屬于CYP4家族。

CYP4基因起源于1000多萬年前，早于脊椎動物和無脊椎動物分化(Bradfieldetal.,1991)，其功能包括促進多種激素及甾醇的生物合成和降解，調節昆蟲發育和形態改變(楊帆和王進軍,2008)，代謝植物有毒物質、多種殺蟲劑,誘導誘變劑等(Gu & Knipple,2013)。薇甘菊頸盲蝽的CYP4家族基因的主要作用可能是調節內源性物質和代謝食物中的有毒物質，PmCYP4應主要分布于中腸、脂肪體和馬氏管中，表達部位及量應無性別差別,在殘體中有幾乎同水平的表達量，但表達量很低。因此，PmCYP4基因可能參與了食物源有毒化合物的代謝，但不是其主要功能。CYP4基因也參與了信息素的失活(Lazardetal.,1990; Ma?bèchecoisneetal.,2005)。在繁殖期，薇甘菊頸盲蝽雌蟲通過氣味分子來引誘雄蟲交尾(王大偉等,2014; 澤桑梓等,2017)，而PmCYP4基因在雌、雄蟲觸角中的表達量很低且無差異，所以PmCYP4基因沒有參與性信息素的失活，但可能參與了環境信息素的失活。

本研究發現，PmCYP4在雌、雄蟲的足中表達量顯著大于其他部位，雄蟲翅膀中的表達量顯著大于雌蟲，GO數據庫對此基因的功能注釋之一是作為配對供體，結合或減少氧分子來增加氧化還原酶活性(GO:0016705)。前期對薇甘菊頸盲蝽的生物特性研究中發現，該蟲的爬行能力較強，飛翔能力較弱，且雄蟲飛翔能力遠強于雌蟲，據此推測，PmCYP4的主要功能可能是編碼與運動相關的酶。