-
松材线虫病又称松树枯萎病(pine wilt disease),是由松材线虫(Bursaphelenchus xylophilus)引起的一种国际性重大检疫病害[1]。该病害具有发病部位隐蔽、传播蔓延迅速、致死率高、危害范围广和防治难度大等特点,因此又被称为松树中的“癌症”[2]。我国于1982年在南京中山陵首次发现松材线虫病[3],疫区范围已扩大到安徽、江西、云南、贵州、广东、湖南、陕西、四川、天津、辽宁、浙江、福建、湖北、广西以及重庆等19个省区市(国家林业和草原局2019年第4号公告),造成的经济损失达上千亿元。松材线虫病的危害和扩散趋势严重威胁着我国松林、众多重要的风景名胜区以及三峡库区的生态安全[4-5]。松材线虫分子致病机理还不清楚,导致松材线虫病防控难以取得突破性进展。
2011年Kikuchi等测定了松材线虫全基因组序列[6],这为松材线虫病致病机理及相关致病基因的研究奠定了良好基础。细胞色素P450(cytochrome p450,CYP450)广泛存在于生物体中,是一类结构和功能相似的血红素同工酶,催化生物体内一系列的氧化代谢反应,参与类固醇、脂肪酸、前列腺素、保幼激素和蜕皮激素等内源化合物的合成,同时在药物及杀虫剂等外源有害物质的解毒过程中具有重要作用[7-10]。松材线虫基因组中含有76个CYP450基因,研究表明这些基因在线虫代谢有害物质的解毒过程中发挥了重要作用[6]。然而,松材线虫细胞色素CYP450基因的具体功能目前还不清楚。
RNA干扰(RNA interference,RNAi)是由外源或内源双链RNA(double stranded RNA,dsRNA)启动,在核糖核酸内切酶的作用下将RNA分子切割为小分子干扰RNA(small interfering RNA,siRNA),然后与mRNA的同源序列特异性结合,在RNA聚合酶的催化下合成新的dsRNA,再由核酸内切酶剪切产生新的siRNA,与mRNA结合后产生新的siRNA,如此经过多次循环,沉默信号逐渐放大,从而引起功能基因表达沉默的现象[11-14]。RNA干扰技术是研究基因功能的重要工具,已被广泛应用于研究动物[15-16]、植物[17]、线虫[18]相关基因的功能。我们前期研究表明,松材线虫入侵黑松后其cyp-13A11基因表达量显著上调,说明该基因在松材线虫致病过程中发挥了重要作用[19]。然而,松材线虫cyp-13A11基因的具体功能目前还不明确。因此,本文以pMD19-T质粒构建了松材线虫cyp-13A11基因dsRNA表达载体,通过浸泡法检测cyp-13A11基因沉默对松材线虫的取食、繁殖和致病力的影响,为探明cyp-13A11基因的功能提供理论依据。
-
供试松材线虫(Bursaphelenchus xylophilus)分离自江西省九江市的黑松病木,将松材线虫接种至长满灰葡萄孢(Botrytis cinerea)的PDA(potato dextrose agar)培养基上,25 ℃恒温培养5 d,采用贝尔曼漏斗法分离线虫,用无菌水清洗3次后将线虫收集至1.5 mL的无菌离心管中,备用。
-
采用天根生化科技有限公司的动物组织总RNA提取试剂盒(RNAprep pure tissue kit)提取松材线虫总RNA,经紫外分光光度计测定浓度后电泳检测RNA质量。将提取的总RNA按照TaKaRa公司反转录试剂盒(Prime script 1st strand cDNA synthesis kit)说明书进行反转录,合成松材线虫cDNA。
-
以上述合成的cDNA为模板,cyp-13A11-F(5′-CAGGGTTACCCGTTCTTA-3′)和cyp-13A11-R(5′-TGATTCCCATCTTTCAGG-3′)为引物,进行PCR反应。反应条件为:94 ℃ 5 min预变性,94 ℃ 30 s,60 ℃ 34 s,72 ℃ 1 min,35个循环,72 ℃ 10 min延伸。扩增产物经1%琼脂糖凝胶检测分析后,按照AxyPrep DNA凝胶回收试剂盒的使用说明书进行切胶回收。回收片段与pEASY-T1载体(北京全式金生物技术有限公司)连接,取10 μL过夜连接产物与100 μL大肠杆菌Trans1-T1感受态细胞混匀,后均匀涂布于含有100 μL/mL Amp+的LB固体培养基上,37 ℃恒温培养过夜,挑取白色菌落进行菌液PCR。将PCR扩增产物送至上海生工公司测序,以确认cyp-13A11基因片段是否插入pEASY-T1载体。
-
以上述构建的表达载体为模板,使用带有T7启动子的引物(cyp-13A11-T7F:5′- TAATACGACTCACTATAGGGCAGGGTTACCCGTTCTTA-3′;cyp-13A11-R:5′- TGATTCCCATCTTTCAGG-3′;cyp-13A11-F:5′- CAGGGTTACCCGTTCTTA-3′;cyp-13A11- T7R:5′- TAATACGACTCACTATAGGGTGATTCCCATCTTTCAGG-3′)进行PCR反应,程序为:95 ℃预变性5 min;95 ℃变性30 s,53 ℃退火30 s,72 ℃延伸45 s,35个循环;72 ℃延伸10 min。PCR产物经1%琼脂糖凝胶电泳检测后,按照MEGAscript RNAi kit protocol试剂盒说明书合成dsRNA,使用紫外分光光度计测定合成的dsRNA浓度,凝胶电泳检测dsRNA的完整性。
-
松材线虫用无菌水清洗后浸泡于0.9%的硫酸链霉素溶液中,2 000 r/min离心3 min,将收集的松材线虫制成悬液。分别吸取80 μL(约3 200条)线虫悬液至等体积的cyp-13A11 dsRNA溶液和ddH2O中,置于转速为80 r/min的摇床中25 ℃浸泡48 h。
分别收集浸泡在ddH2O和dsRNA溶液中的松材线虫,用无菌水清洗3次后,从中挑取30对成虫接种至长满灰葡萄孢的培养皿上,每个处理接种3皿,25 ℃恒温培养6 d,拍照观察松材线虫取食灰葡萄孢的情况。培养6 d后用贝尔曼漏斗法分离不同处理的松材线虫,2 000 r/min离心3 min,弃上清收集线虫并在显微镜下统计每皿线虫的数量。
将cyp-13A11 dsRNA干扰后的线虫采用人工皮接法分别接种至3年生黑松上,每株接种50 μL(约2 000条)线虫悬液,分别将等体积ddH2O浸泡后的线虫悬液和cyp-13A11 dsRNA溶液(无线虫)接种至黑松作为对照。每个处理接种9株松苗,观察发病情况,持续30 d,定期统计发病率。
-
分别提取cyp-13A11 dsRNA溶液和ddH2O浸泡后松材线虫的总RNA,按照prime script 1st strand cDNA synthesis kit(TaKaRa)试剂盒的使用说明反转录合成cDNA的第一条链。依据GenBank上已发表的松材线虫cyp-13A11基因序列(登录号:NM_067303.1)设计qPCR引物:cyp-13A11-qPCR-F(5′-TGCTGCTCAATGTGCCTC-3′)和cyp-13A11-qPCR-R(5′-GACTATTGCGAATCCCACG-3′),以Actin基因为内参照,采用transstart green qPCR supermix(北京全式金生物技术有限公司)试剂盒进行qRT-PCR,测定松材线虫cyp-13A11基因的干扰效率,每个实验重复3次。
-
实验数据采用SPSS 13.0统计软件进行差异显著性分析。
-
以松材线虫总RNA反转录得到的cDNA为模板进行PCR扩增。扩增结束后,将PCR产物加入1%琼脂糖凝胶中检测分析,结果如图1所示,扩增片段大小介于1 200 ~ 1 300 bp之间,与预期大小1 269 bp相符。PCR产物按AxyPrep DNA凝胶回收试剂盒说明进行切胶回收,4 ℃保存备用。
-
将克隆的松材线虫cyp-13A11基因片段连接至pEASY-T1载体,连接产物转化至大肠杆菌Trans1-T1感受态细胞中,均匀涂布于含有Amp、IPTG、Xgal的固体LB培养基平板上,37 ℃培养12 h,挑取白色菌落进行菌液PCR。取10 μL PCR扩增产物进行琼脂糖凝胶电泳,结果如图2所示,电泳条带明亮清晰,大小与目的条带一致。同时将PCR产物送上海生工测序,将测序结果与NCBI网站的已知序列进行BLAST比对分析,发现二者之间的序列同源性为99.19%,表明松材线虫cyp-13A11基因已成功克隆至pEASY-T1载体。
-
以构建的RNA干扰表达载体为模板,按照MEGAscript RNAi kit protocol试剂盒说明书步骤合成松材线虫cyp-13A11基因dsRNA。经紫外分光光度计检测,合成的dsRNA浓度为2.6 μg/μL。取1 μL进行琼脂糖凝胶电泳,检测结果如图3所示,目的条带明亮、单一,并且与预期片段大小一致,说明合成的dsRNA可用于后续的基因沉默实验。
-
采用qRT-PCR检测RNA干扰对松材线虫cyp-13A11基因表达量的影响,结果如图4所示,将对照组的cyp-13A11基因的表达量设定为1,RNAi处理组松材线虫cyp-13A11基因的表达量为0.01,与对照组相比其表达量下调显著(P < 0.01),表明dsRNA干扰可有效抑制松材线虫cyp-13A11基因的表达。
-
松材线虫经ddH2O和dsRNA浸泡处理后接种至长满灰葡萄孢的培养基上,25 ℃恒温培养6天。松材线虫取食情况如图5所示,接种第3天ddH2O处理的线虫取食面积明显大于cyp-13A11 dsRNA处理的取食面积(图5A);接种第6天对照组松材线虫已将灰葡萄孢取食殆尽,而cyp-13A11 dsRNA处理组的培养基表面仍存在大量灰葡萄孢菌丝(图5B)。采用贝尔曼漏斗法收集培养6天后的线虫并统计其数量。由图6可知,cyp-13A11 dsRNA处理组与对照组的线虫数量分别为2 040和8 733条。上述结果表明对cyp-13A11基因的干扰影响了松材线虫取食和繁殖。
-
将干扰后和对照组的松材线虫接种至3年生的黑松上,观察黑松发病情况并统计其发病率。结合发病率和发病情况(表1和图7)可知,在接种10 d后,与ddH2O浸泡线虫处理相比,cyp-13A11 dsRNA浸泡线虫处理的黑松发病率减少了16.6%;在接种20 d后,ddH2O浸泡线虫处理黑松发病率为44.4%,而cyp-13A11 dsRNA浸泡线虫处理黑松发病率为33.3%;直至接种30 d,接种ddH2O和cyp-13A11 dsRNA浸泡线虫的黑松全部枯萎,发病率达到100%。与此同时,接种cyp-13A11 dsRNA溶液的黑松在整个实验期间生长状况良好,发病率为零。以上结果表明,将松材线虫cyp-13A11基因进行干扰后,线虫的致病能力减弱。
表 1 不同接种处理黑松的发病情况
Table 1. Disease status of Pinus thunbergii under different treatments
接种时间
Inoculating
time/d黑松发病率 Wilting rate of Pinus thunbergii/% 对照组1
Control group 1RNAi处理组
RNAi treatment group对照组2
Control group 210 0 5.6 22.2 20 0 33.3 44.4 30 0 100.0 100.0 注:对照组1:cyp-13A11dsRNA接种至黑松苗;RNAi:cyp-13A11 dsRNA干扰后的线虫悬液接种至黑松苗;对照组2:ddH2O浸泡后的线虫悬液接种至黑松苗。Notes: control group 1, cyp-13A11 dsRNA is inoculated into Pinus thunbergii; RNAi, the suspension of nematodes after interference with cyp-13A11 dsRNA is inoculated into Pinus thunbergii; control group 2, the ddH2O soaked nematode suspension is inoculated into Pinus thunbergii. -
松材线虫入侵松树后,寄主会产生一系列的防御反应,分泌出大量的蒎浠类和芳香族化合物等次生代谢产物,以抑制松材线虫的入侵[20]。在植物病原线虫的寄生过程中,对有害物质的解毒作用是线虫适应寄主体内环境的重要生理过程,也是线虫能否成功寄生的关键。因此,松材线虫为了能成功寄生于松树体内,必须通过解毒基因代谢转化松树产生的有毒代谢产物。有研究表明,CYP450基因在松树线虫的解毒过程中发挥了重要作用[6]。本研究成功构建松材线虫cyp-13A11基因干扰载体,通过RNAi方法对松材线虫cyp-13A11基因进行干扰,发现cyp-13A11基因沉默后松材线虫对灰葡萄孢的取食面积减小,线虫繁殖数量下降了76.6%。Benenati等通过显微注射方法对秀丽隐杆线虫(Caenorhabditis elegans)的cyp-31A3和cyp-31A2基因进行了RNA干扰,发现这两个基因沉默后对线虫卵壳的形成和胚胎的发育具有重要的影响[21]。我们推测松材线虫cyp-13A11基因干扰后不利于虫卵的发育,进而影响了线虫的繁殖和取食速度。
在接种前期(10 d),cyp-13A11基因干扰对松材线虫致病力的影响较为明显,接种干扰后松材线虫的黑松发病率仅为5.6%,而接种ddH2O浸泡松材线虫的黑松发病率为22.2%;在接种中期(20 d),接种RNA干扰后和ddH2O浸泡松材线虫的黑松发病率差距较小;在接种后期(30 d),接种松材线虫处理的黑松发病率均达到100%。由此可知,在接种前期,cyp-13A11基因干扰可以明显降低松材线虫的致病力,但随着接种时间的推移,干扰效果逐渐下降。我们的前期研究发现,RNAi对松材线虫F1 ~ F5代干扰效果明显,但对F6及以后代干扰效果较小[22]。因此,我们推测随着接种时间的推移,cyp-13A11基因RNA干扰效果逐渐消失,导致在接种后期cyp-13A11 dsRNA处理和ddH2O浸泡处理的黑松发病率趋于一致。大量研究显示,CYP450基因参与代谢过程和抗性机制主要通过表达量上调来实现。Maja等[23]研究表明,当秀丽隐杆线虫(Caenorhabditis elegans)受到外源有毒物质侵害时,其cyp-33D3基因的表达量上调显著。王璇等将松材线虫接种至5年生的马尾松上,发现松材线虫cyp-33C9和cyp-33C4基因的表达量上调,推测松材线虫CYP450基因可能参与了松材线虫蒎烯类物质代谢过程,可能是松材线虫致病相关基因之一[24]。目前关于线虫CYP450基因功能的报道相对较少,其研究主要集中在CYP450基因表达量方面。我们通过RNA干扰技术对松材线虫cyp-13A11基因的功能进行了研究,结果表明cyp-13A11基因表达被沉默后,减缓了松材线虫的取食速度,降低了松材线虫的繁殖能力和致病力,说明CYP450基因在松材线虫繁殖和致病过程中都具有重要作用,这为进一步阐明松材线虫分子致病机理提供了数据支撑,同时为松材线虫的生物防治提供了理论依据。然而,cyp-13A11基因对松材线虫繁殖和致病力的具体调节机制还不清楚,需要进一步深入研究。
-
本研究成功构建了松材线虫cyp-13A11基因RNA干扰载体,通过RNA干扰抑制了松材线虫cyp-13A11基因的表达。研究结果表明cyp-13A11基因表达沉默后减缓了松材线虫的取食速度,线虫的繁殖数量也显著下降(P < 0.001)。在松材线虫接种至黑松前期(10 d)和中期(20 d),cyp-13A11 dsRNA处理黑松的发病率比ddH2O处理分别降低了16.6%和11.1%,而在接种后期(30 d),RNA干扰对松材线虫致病力没有影响,两个处理的黑松发病率均达到了100%。
Construction and its function analysis of RNA interference vector of cytochrome cyp-13A11 gene in Bursaphelenchus xylophilus
-
摘要:
目的 构建松材线虫cyp-13A11基因RNA干扰载体,研究cyp-13A11基因的功能,为松材线虫病的生物防治提供理论依据。 方法 本文通过设计特异性引物,扩增松材线虫cyp-13A11基因片段,构建至pEASY-T1干扰载体上并转入Trans1-T1大肠杆菌菌株。采用浸泡法对松材线虫进行RNA干扰,通过qRT-PCR检测cyp-13A11基因的表达,将干扰后的线虫分别接种至长满灰葡萄孢菌丝的培养皿中和黑松苗上,测定RNA干扰后线虫的取食、繁殖、致病力和干扰效率等指标。 结果 成功构建了带有pEASY-T1干扰载体的Trans1-T1菌株并且能够合成cyp-13A11 dsRNA,通过RNA干扰抑制了松材线虫cyp-13A11基因的表达。松材线虫RNA干扰后接种至灰葡萄孢第3天,cyp-13A11 dsRNA处理松材线虫的取食面积较小,ddH2O处理松材线虫取食面积达到了整体的2/3;接种第6天,cyp-13A11 dsRNA处理松材线虫的取食面积为整体的1/3,而ddH2O处理松材线虫已将灰葡萄孢菌丝取食殆尽。ddH2O处理松材线虫的繁殖数量比cyp-13A11 dsRNA处理高4.28倍。松材线虫接种至黑松苗第10天,ddH2O处理黑松发病率为22.2%,而cyp-13A11 dsRNA处理黑松发病率为5.6%;接种至20天,ddH2O处理的黑松发病率为44.4%,cyp-13A11 dsRNA处理的黑松发病率为33.3%;直至接种30天,ddH2O和cyp-13A11 dsRNA处理的黑松全部枯萎,发病率达到100%。 结论 成功构建了松材线虫cyp-13A11基因RNA干扰载体,cyp-13A11基因表达沉默后影响了松材线虫的取食,降低了松材线虫的繁殖能力和致病力。 -
关键词:
- 松材线虫 /
- cyp-13A11基因 /
- 载体构建 /
- RNA干扰 /
- 功能分析
Abstract:Objective In order to reveal the function of cyp-13A11 gene in Bursaphelenchus xylophilus, the RNA interference vector was constructed and the function of cyp-13A11 gene of Bursaphelenchus xylophilus was analyzed to provide theoretical basis for biological control of pine wilt diseases. Method Specific primers were designed to amplify the cyp-13A11 gene fragment of Bursaphelenchus xylophilus. The products of amplification were constructed into pEASY-T1 interference vector, which was further transferred into Trans1-T1 escherichia coli strain. RNA interference efficiency was determined by qRT-PCR. The feeding reproduction and pathogenicity of pine wood nematodes were investigated after the nematodes soaked in cyp-13A11 dsRNA solution were inoculated on Botrytis cinerea and pine trees, respectively. Result Trans1-T1 strain with pEASY-T1 interference vector was successfully constructed and cyp-13A11 dsRNA was synthesized. RNA interference inhibited the expression of cyp-13A11 gene in Bursaphelenchus xylophilus. In addition, the feeding area was less in cyp-13A11 dsRNA treatment than that in ddH2O treatment, as two thirds of the Botrytis cinerea were fed by nematodes in ddH2O treatment on day 3. After 6 days, one third of Botrytis cinerea was fed by nematodes in cyp-13A11 dsRNA treatment, while almost all of the Botrytis cinerea were exhausted by nematodes in ddH2O treatment. The reproduction number of nematodes in ddH2O treatment was 4.28 times higher than that in cyp-13A11 dsRNA treatment. Furthermore, the wilting rates in both ddH2O and cyp-13A11 dsRNA treatments were 22.2% and 5.6% after Bursaphelenchus xylophilus being inoculated on pine trees for 10 days, respectively. The wilting rates of pine trees in both ddH2O and cyp-13A11 dsRNA treatments were 44.4% and 33.3% after 20 days, respectively. Obviously, the wilting rates were 100% in both treatments after 30 days. Conclusion The RNA interference vector of cyp-13A11 gene in pine wood nematode was successfully constructed. Silencing of cyp-13A11 gene exhibited an important effect on feeding of Bursaphelenchus xylophilus, as well as reduced the reproduction ability and pathogenicity of Bursaphelenchus xylophilus. -
Key words:
- Bursaphelenchus xylophilus /
- cyp-13A11 gene /
- vector construction /
- RNA interference /
- function analysis
-
图 5 松材线虫在灰葡萄孢上的取食情况
A. 接种3 d;B. 接种6 d。对照组:ddH2O浸泡松材线虫;RNAi处理组:cyp-13A11 dsRNA浸泡松材线虫。下同。A,inoculated for three days; B,inoculated for six days. Control group,Bursaphelenchus xylophilus soaked in ddH2O; RNAi treatment group,Bursaphelenchus xylophilus soaked in cyp-13A11 dsRNA. The same below.
Figure 5. Feeding condition of Bursaphelenchus xylophilus on Botrytis cinerea
表 1 不同接种处理黑松的发病情况
Table 1. Disease status of Pinus thunbergii under different treatments
接种时间
Inoculating
time/d黑松发病率 Wilting rate of Pinus thunbergii/% 对照组1
Control group 1RNAi处理组
RNAi treatment group对照组2
Control group 210 0 5.6 22.2 20 0 33.3 44.4 30 0 100.0 100.0 注:对照组1:cyp-13A11dsRNA接种至黑松苗;RNAi:cyp-13A11 dsRNA干扰后的线虫悬液接种至黑松苗;对照组2:ddH2O浸泡后的线虫悬液接种至黑松苗。Notes: control group 1, cyp-13A11 dsRNA is inoculated into Pinus thunbergii; RNAi, the suspension of nematodes after interference with cyp-13A11 dsRNA is inoculated into Pinus thunbergii; control group 2, the ddH2O soaked nematode suspension is inoculated into Pinus thunbergii. -
[1] Xia Y, Qi Y, Yu X, et al. Nematicidal effect against Bursaphelenchus xylophilus of harmine quaternary ammonium derivatives, inhibitory activity and molecular docking studies on acetylcholinesterase[J]. European Journal of Plant Pathology, 2019, 153(1): 239−250. doi: 10.1007/s10658-018-1559-8. [2] 朱克恭, 姚仕义. 松材线虫病在我国蔓延的严重性[J]. 林业科技开发, 1992(3):6−7. Zhu K G, Yao S Y. The spread of pine wood nematode disease in China is serious[J]. China Forestry Science and Technology, 1992(3): 6−7. [3] 程瑚瑞, 林茂松, 黎伟强. 南京松树线虫研究[J]. 南京农业大学学报, 1983, 6(4):84. Cheng H R, Lin M S, Li W Q. Study on Bursaphelenchus xylophilus in Nanjing[J]. Journal of Nanjing Agricultural College, 1983, 6(4): 84. [4] Gao R, Shi J, Huang R, et al. Effects of pine wilt disease invasion on soil properties and Masson pine forest communities in the Three Gorges Reservoir Region, China[J]. Ecology and Evolution, 2015, 5(8): 1702−1716. doi: 10.1002/ece3.1326 [5] 理永霞, 张星耀. 松材线虫病致病机理研究进展[J]. 环境昆虫学报, 2018, 40(2):231−241. Li Y X, Zhang X Y. Research advance of pathogenic mechanism of pine wilt disease[J]. Journal of Environmental Entomology, 2018, 40(2): 231−241. [6] Kikuchi T, Cotton J A, Dalzell J J, et al. Genomic insights into the origin of parasitism in the emerging plant pathogen Bursaphelenchus xylophilus[J]. PLoS Pathogens, 2011, 7(9): e1002219. doi: 10.1371/journal.ppat.1002219. [7] Coon M J. Cytochrome P450: nature’s most versatile biological catalyst[J]. Annual Reviews Pharmacol Toxicol, 2005, 45: 1−25. doi: 10.1146/annurev.pharmtox.45.120403.100030. [8] Thomas J H. Rapid birth-death evolution specific to xenobiotic cytochrome P450 genes in vertebrates[J]. PLoS Genetics, 2007, 3(5): e67. doi: 10.1371/journal.pgen.0030067. [9] Scott J G, Liu N, Wen Z. Insect cytochromes P450: diversity, insecticide resistance and tolerance to plant toxins[J]. Comparative Biochemistry and Physiology Part C: Pharmacology, Toxicology and Endocrinology, 1998, 121(1−3): 147−155. doi: 10.1016/S0742-8413(98)10035-X. [10] 郭亭亭, 姜辉, 高希武. 昆虫细胞色素P450基因的多样性, 进化及表达调控[J]. 昆虫学报, 2009, 52(3):301−311. Guo T T, Jiang H, Gao X W. Diversity, evolution and expression regulation of cytochrome P450 monooxygenase genes in insects[J]. Acta Entomologica Sinica, 2009, 52(3): 301−311. [11] 邹丰才, 宋慧群, 陈宁, 等. RNA干扰技术及其在寄生线虫基因功能研究中的应用[J]. 中国兽医科学, 2006, 36(1):80−84. Zou F C, Song H Q, Chen N, et al. RNA interference technique and its application to studies of gene functions of parasitic nematodes[J]. Veterinary Science in China, 2006, 36(1): 80−84. [12] 陈贵元, 谭德勇. RNA干扰的研究现状与应用前景[J]. 大理学院学报, 2012, 11(10):52−55. Chen G Y, Tan D Y. The research status and application prospect of RNA interference[J]. Journal of Dali University, 2012, 11(10): 52−55. [13] Fire A, Xu S Q, Montgomery M K, et al. Potent and specific genetic interference by double-stranded RNA in Caenorhabditis elegans[J]. Nature, 1998, 391: 806. doi: 10.1038/35888. [14] 郭葆玉. RNA干扰技术[J]. 药物生物技术, 2008, 15(2):137−141. Guo B Y. RNA interference technology[J]. Pharmaceutical Biotechnology, 2008, 15(2): 137−141. [15] 夏欢, 邹声雷, 李超程, 等. 家兔Zfx基因shRNA干扰载体构建及验证[J]. 现代畜牧兽医, 2019(5): 14−19. Xia H, Zou S L, Li C H, et al. The shRNA interference vector construction and validation from Lepus Zfx gene[J]. Modern Journal of Animal Husbandry and Veterinary Medicine, 2019 (5): 14−19. [16] 王越, 张苏芳, 徐瑶, 等. 美国白蛾几丁质酶细菌表达的RNA干扰载体构建及其介导的RNA干扰[J]. 林业科学研究, 2019, 32(2):1−8. Wang Y, Zhang S F, Xu Y, et al. Construction of the expression vector and RNAi mediated by bacteria expressed dsRNA of chitinase gene from Hyphantria cunea[J]. Forest Research, 2019, 32(2): 1−8. [17] 陈鸿鹏, 谭晓风, 谢耀坚, 等. 油茶FAD2 基因的植物表达载体和RNA干扰载体构建及其功能分析[J]. 植物研究, 2015, 35(3):347−354. Chen H P, Tan X F, Xie Y J, et al. Constructing and function of plant expression vector and RNA interference vector of a FAD2 gene from Camellia oleifera[J]. Bulletin of Botanical Research, 2015, 35(3): 347−354. [18] 王守先, 牛宝龙, 沈卫锋, 等. 松材线虫RNA聚合酶基因的RNA干扰研究[J]. 浙江农业科学, 2007, 49(6):690−693. Wang S X, Niu B L, Shen W F, et al. RNA interference of the RNA polymerase gene in Bursaphelenchus xylophilus[J]. Journal of Zhejiang Agricultural Sciences, 2007, 49(6): 690−693. [19] Qiu X, Wu X, Huang L, et al. Specifically expressed genes of the nematode Bursaphelenchus xylophilus involved with early interactions with pine trees[J]. PLoS ONE, 2013, 8(10): e78063. doi: 10.1371/journal.pone.0078063. [20] Futai K. Pine wood nematode, Bursaphelenchus xylophilus[J]. Annual Review of Phytopathology, 2013, 51: 61−83. doi: 10.1146/annurev-phyto-081211-172910. [21] Benenati G, Penkov S, Müller-Reichert T, et al. Two cytochrome P450s in Caenorhabditis elegans are essential for the organization of eggshell, correct execution of meiosis and the polarization of embryo[J]. Mechanisms of Development, 2009, 126(5−6): 382−393. doi: 10.1016/j.mod.2009.02.001. [22] Kikuchi T, Shibuya H, Aikawa T, et al. Cloning and characterization of pectate lyases expressed in the esophageal gland of the pine wood nematode Bursaphelenchus xylophilus[J]. Molecular Plant-Microbe Interactions, 2006, 19(3): 280−287. doi: 10.1094/MPMI-19-0280. [23] Studencka M, Wesołowski R, Opitz L, et al. Transcriptional repression of Hox genes by C. elegans HP1/HPL and H1/HIS-24[J]. PLoS Genetics, 2012, 8(9): e1002940. doi: 10.1371/journal.pgen.1002940. [24] 王璇, 理永霞, 刘振宇, 等. 松材线虫CYP450基因与松树蒎烯类物质代谢的相关性[J]. 林业科学, 2017, 53(6):105−110. Wang X, Li Y X, Liu Z Y, et al. Relationship between the cytochrome P450 gene of pine wood nematode and the accumulation of pine pinene[J]. Scientia Silvae Sinicae, 2017, 53(6): 105−110. -