飞蝗核受体基因LmE75的分子特性和功能分析

2020-07-31刘晓健郭俊张学尧马恩波张建珍

中国农业科学 2020年11期

刘晓健，郭俊，2，张学尧，马恩波，张建珍

飞蝗核受体基因的分子特性和功能分析

刘晓健1，郭俊1，2，张学尧1，马恩波1，张建珍1

（1山西大学应用生物学研究所，太原 030006；2山西大学生命科学学院，太原 030006）

【】核受体（nuclear receptor，NR）是一类配体依赖型的转录因子，在昆虫生长发育过程中发挥重要作用，本文以重要农业害虫飞蝗（）为研究对象，分析核受体基因的分子特性以及功能，为害虫防治提供新分子靶标。基于飞蝗转录组数据库获得3条的cDNA序列，运用Blast和ExPASy网站中相关软件分析各个基因的开放阅读框、预测氨基酸序列、理化特性及保守结构域；使用MEGA 6.0 软件中邻接法（neighbor-joining，NJ）与其他昆虫同源序列进行聚类分析。利用Primer 3.0软件设计3个亚型的特异性表达引物，运用 reverse transcription-quantitative PCR（RT-qPCR）技术分析3个亚型在5龄若虫不同组织及4—5龄不同发育日龄的表达特性。体内注射蜕皮激素（20-hydroxyecdysone，20E）诱导和干扰20E受体基因的表达后，RT-qPCR检测的表达，明确是否受20E调控。利用RNAi（RNA interference）技术，将体外合成的亚型共同区域的dsRNA注射至5龄飞蝗若虫体内，48 h后提取体壁RNA，反转录为cDNA，运用RT-qPCR技术检测亚型的沉默效率并观察表型；之后进一步利用H&E（hematoxylin and eosin）染色方法观察对飞蝗表皮结构的影响。基于飞蝗转录组数据库获得3条的cDNA序列，根据其结构特征命名为、和，GenBank登录号分别为MN584732、 MN584733和MN584734。3个LmE75亚型均具有典型的核受体结构域，其中DBD（DNA binding domain）和LBD（ligand binding domain）在各昆虫中高度保守；E75聚类分析结果显示不同目的昆虫类群E75亚型优先聚为一支；RT-qPCR结果显示3个亚型具有不同的组织和发育表达特性。在体壁和肌肉中高表达；在体壁中高表达，其次是前肠；特异性在体壁高表达；在4龄和5龄期均呈现先升高后降低的表达趋势；在N4D3和N5D5表达量最高，蜕皮前后表达量均较低；在4龄和5龄期均呈现逐渐升高的表达趋势，在蜕皮前表达量最高。20E诱导3 h和6 h后，处理组的表达量较对照组分别上调了4.2倍和4.5倍。利用RNAi技术干扰48 h后的表达量较对照组下调59%，表明受20E调控。在5龄第1天注射ds后的飞蝗具有明显的发育延迟现象，若虫无法蜕皮至成虫并最终100%死亡。体壁H＆E染色结果显示，注射ds组飞蝗没有发生皮层溶离现象。飞蝗中有3个亚型，三者具有不同的组织和发育表达特性；RNAi试验表明，在飞蝗蜕皮过程中发挥重要作用。

飞蝗；核受体；；蜕皮激素；RNA干扰

0 引言

【研究意义】核受体（nuclear receptor，NR）是一类在生物体中广泛分布的配体依赖型的转录因子[1]，可通过多个信号途径调控机体的胚胎及胚后发育、衰老、昼夜节律等生命过程[2]。飞蝗（）是重要的农业害虫，围绕其核受体基因开展研究工作，对阐明核受体在飞蝗蜕皮过程中的作用机制以及基于RNA干扰（RNAi）的飞蝗防治具有重要意义。【前人研究进展】昆虫经过变态发育获得飞翔和生殖能力[3]。前胸腺分泌的蜕皮激素（20-hydroxyecdysone，20E）与咽侧体分泌的保幼激素（juvenile hormone，JH）是昆虫变态的主要因素[4]。蜕皮激素20E调控昆虫蜕皮的分子机制已在完全变态昆虫中有较多研究，20E与蜕皮激素受体（EcR）和超气门蛋白（USP）的复合物结合成为三聚体，启动一系列下游转录因子的表达，进而调控昆虫蜕皮等生理活动[5]。参与20E信号传导的基因包括转录因子E74、Broad和一些核受体。核受体含有5个典型结构域：A/B域（转录激活域，trans-activation domain）、C域（DNA结合域，DNA binding domain，DBD）、D域（铰链域，hinge domain）、E域（配体结合域，ligand binding domain，LBD）和F域。可分为7个亚家族：NR0、NR1、NR2、NR3、NR4、NR5和NR6[6]。E75核受体属于NR1亚家族，命名为NR1D3。在黑腹果蝇（）中研究发现，有4个亚型，是20E信号传导通路的早期响应基因，其纯合的缺失品系致使1龄幼虫期发育延迟，不能蜕皮为2龄幼虫而死亡[7]。迄今为止，已在埃及伊蚊（）[8]、烟草天蛾（）[9]、大蜡螟（）[10]、云杉芽卷蛾（）[11]、家蚕（）[12]、印度谷螟（）[13]、西方蜜蜂（）[14]、赤拟谷盗（）[15]、马铃薯甲虫（）[16]、德国小蠊（）[17]等多种昆虫中克隆获得。已有研究表明在昆虫生长发育过程中发挥重要作用，如德国小蠊中有5个亚型，注射共同区域的dsRNA后大部分个体无法正常蜕皮而死亡[17]。马铃薯甲虫中有3个亚型，各亚型均受20E调控，利用RNAi技术将沉默后可影响蜕皮激素20E和保幼激素JH的滴度从而调控蜕皮激素和保幼激素信号途径，最终导致马铃薯甲虫化蛹困难而死亡[16]。【本研究切入点】目前在完全变态昆虫的研究中取得了较大进展[7-16]，但该基因在不完全变态昆虫飞蝗中的作用机制有待深入解析。【拟解决的关键问题】通过RNAi结合H&E染色方法分析飞蝗核受体基因5的生物学功能，丰富蜕皮激素调控昆虫蜕皮的理论体系，为基于RNAi的害虫控制提供新的靶标基因。

1 材料与方法

试验于2018—2019年在山西大学应用生物学研究所完成。

1.1 供试昆虫

飞蝗虫卵由河北省沧州市建民虫业公司提供。将虫卵置于孵化盒中进行孵化，期间喷洒适量水以保持沙土潮湿，孵化出的1龄若虫转移至25 cm×25 cm×25 cm的纱笼中，置于温度为（30±2）℃，相对湿度为（40±10）%，光照14 h、黑暗10 h的养虫室中，饲喂新鲜的小麦苗，待其生长发育至3龄若虫期开始辅以麦麸饲喂。

1.2 LmE75 cDNA序列分析

1.2.1cDNA序列搜索及分析通过搜索飞蝗转录组数据库获得3条的cDNA序列。利用NCBI中的Open Reading Frame Finder（ORF Finder）（http://cms.ncbi.nlm.nih.gov/orf finder/）软件找到的全长开放阅读框。利用ExPASy网站中（https://www.expasy.org/）相关软件对各个基因的开放阅读框序列进行翻译、预测各蛋白分子量（Mw）及等电点（pI）。SMART（http://smart. embl-heidelberg. de/）分析预测其蛋白结构域。运用GENEDOC软件将、和特异性的A/B域核苷酸进行比对，分析三者的核苷酸一致度；运用GENEDOC软件将、和氨基酸序列与马铃薯甲虫（、和）和德国小蠊（、和）氨基酸序列进行比对，分析其DNA结合域（DBD）和配体结合域（LBD），而两个结构域通过不保守的绞链结构域连接。

1.2.2 昆虫E75系统发育分析在NCBI网站上下载其他物种E75蛋白的同源序列，运用MEGA 6.0软件中NJ（neighbor-joining）法构建不同昆虫E75的系统进化树。不同昆虫E75的GenBank登录号见表1。

1.3 LmE75亚型在不同组织及发育时期的表达特性分析

1.3.1 样品收集为了研究亚型在不同组织和发育时间的表达特性，选用5龄第2天（N5D2）飞蝗若虫，在体视显微镜下解剖触角、脑、前肠、胃盲囊、中肠、后肠、精巢、卵巢、体壁、脂肪体、肌肉和翅12个组织；选取4—5龄每天（N4D1—N4D5和N5D1—N5D7）若虫第2—3腹节处体壁立即冻存于液氮中，-80℃保存备用。以上所有样品均取4个生物学重复，每个生物学重复至少4头若虫。

1.3.2 总RNA提取、cDNA合成及RT-qPCR 提取上述组织与发育时间的样品总RNA，具体步骤均按照 RNAisoTMPlus（TaKaRa）试剂说明书操作，用1.5%的琼脂糖凝胶电泳检测所提取RNA的质量，利用NanoDrop 2000定量各样品的RNA浓度，并检验A260/A280和A260/A230的值。以1 μg总RNA为模板，根据Reverse Transcriptase M-MLV（RNase H-）（TaKaRa）反转录试剂盒说明书，合成cDNA。运用Primer 3.0软件设计3个亚型的特异性表达引物，以作为内参[18-19]，引物信息见表2。利用RT-qPCR方法检测和在不同组织和发育时间的表达特性。在Applied Biosystems 7300 real-time PCR system上进行RT-qPCR，PCR反应体系按照SYBR ® Green real-time PCR Master Mix（TaKaRa）说明书进行配置，扩增程序第一步为预变性10 min，第二阶段为95℃ 15 s，退火温度60℃ 1 min，40个循环；DNA熔解曲线为65—95℃，每步上升1℃。每个样品均设置3个生物学重复和2 个技术重复，基因相对表达量采用2-ΔCt法[20]进行分析，每个组织与发育时间点的数据以平均数±标准差表示。采用SPSS中的 Tukey’s HSD 进行差异显著性检验，柱形图上不同字母代表基因表达差异显著。

表1 本研究中用于构建系统进化树的物种和E75的 GenBank 登录号

表2 所用RT-qPCR和RNAi的引物信息

T7序列T7 sequence: taatacgactcactataggg

1.4 20E诱导表达分析

1.4.1 20E诱导为了检测对20E诱导的响应，配置10%的乙醇并将20E（Sigma）稀释至浓度为5 μmol·L-1保存备用，具体参考文献[21]的方法；挑取健康的5龄第2天（N5D2）若虫注射10 μg 20E，同时注射等体积的10%的乙醇作为对照组，分别在注射3、6、12 h后取飞蝗第2—3体节的体壁提取RNA，并反转录为cDNA。每个时间点的对照组和处理组均设置5个生物学重复，4头试虫为1个生物学重复，共计20头试虫。采用RT-qPCR方法检测在对照组和20E处理组各时间点中的表达量。

1.4.2 干扰20E核受体后的表达变化为了进一步检测响应20E诱导表达情况，首先合成的dsRNA，具体参考文献[22]的方法。在N5D2若虫注射10 μg双链 RNA（ds），以注射的双链RNA（ds）作为对照组。注射48 h后取飞蝗第2—3体节的体壁提取RNA，并反转录为cDNA。对照和注射ds组均设置5个生物学重复，4头试虫为1个生物学重复，共计20头试虫。采用RT-qPCR方法检测和的表达量。1.4.1和1.4.2的数据采用student-test方法进行差异表达分析，*、**和***分别表示在＜0.05、＜0.01、＜0.001水平差异显著。

1.5 基于RNAi的LmE75功能分析

1.5.1 dsRNA合成与注射根据3个亚型共同区域的核苷酸序列，在E-RNAi（http://www.dkfz.de/ signaling/e-rnai3/idseq.php）网站设计dsRNA引物，引物信息见表2。以稀释后的和质粒DNA分别为模板进行PCR扩增以获得用于体外转录dsRNA的模板，按照2×Taq PCR MasterMix（TIANGEN）说明书进行PCR反应体系的配置及反应程序的设置。PCR产物的回收纯化按照FastPure Gel DNA Extraction Mini Kit（Vazyme）说明书操作，之后用1.5%的琼脂糖凝胶电泳检测扩增条带的单一性。利用纯化后的DNA产物为模板，参照T7 RiboMAXTMExpress RNAi System（Promega）试剂盒说明书步骤分别合成和的dsRNA，利用NanoDrop 2000定量各样品的dsRNA浓度和纯度。用ddH2O将合成的dsRNA稀释至终浓度为2 μg·μL-1，置于-80℃冰箱保存备用。挑取N5D1的若虫进行ds和ds的注射。每头试虫的dsRNA注射量均为10 μg，设置3个生物学重复，20头若虫为1个生物学重复，共计60头试虫。注射完毕后，将所有试虫放置于养虫室内细心饲养用于后续试验。

1.5.2 RNAi后表达量检测及表型观察注射ds和ds48 h后分别收集对照组和处理组的第2—3体节体壁提取RNA，反转录为cDNA，采用 RT-qPCR方法检测的沉默效果，设置3个生物学重复，每组3头试虫，2个技术重复，RT-qPCR方法与1.3的检测方法相同，采用2-ΔCt法对的相对表达量进行计算，数据处理及统计方法参考1.4，剩余试虫用来观察表型并拍照。

1.5.3 H＆E 染色为了观察干扰后对飞蝗表皮形成的影响，分别解剖对照组脊线开裂时若虫和处理组相应时间点以及最终死亡时的飞蝗体壁组织。取其第4—5腹节体壁，用3%戊二醛固定，4℃过夜以制备石蜡切片。固定后的体壁样品经过脱水和石蜡包埋后进行5 μm切片，苏木精和伊红（H＆E）进行染色，具体方法参照文献[23]，封片后利用Olympus BX51显微镜（Japan）观察注射ds和ds的体壁结构。

2 结果

2.1 LmE75序列特性

根据其结构特征，将获得的3个飞蝗命名为、和，GenBank登录号分别为MN584732、MN584733和MN584734。、和的cDNA序列长度分别为2 238、2 280和2 271 bp，编码745、759和756个氨基酸；使用ExPaSy网站上的相关软件预测蛋白分子量及等电点，LmE75A、LmE75B和LmE75C的相对分子量分别为81.60、82.55和82.37 kD，其pI均为碱性，分别为8.43、8.08和8.32。3个LmE75亚型具有典型的核受体结构域：A/B域、C域（DBD）、D域、E域（LBD）和F域（图1-A）。飞蝗3个E75亚型N端A/B域有较小的差异，和以及的A/B域核苷酸一致度为45.0%和47.7%，和的A/B域核苷酸一致度为64.1%（图1-B）。LmE75A、LmE75B和LmE75C的DBD和LBD结构域与其他昆虫的DBD和LBD结构域具有高度一致度，其中DBD结构域包含两个锌指区域，每一个具有4个半胱氨酸残基，LBD结构域中包含4个参与血红素结合的重要氨基酸残基（图1-C）。

2.2 系统发育分析

E75的系统发育分析表明，发育树分别形成了鞘翅目类群、半翅目类群、蜚蠊目类群、膜翅目类群、直翅目类群、鳞翅目类群和双翅目类群，其中半翅目类群分为两支，绿盲蝽（）和温带臭虫（）的E75聚为一支，褐飞虱（）两个E75亚型单独聚为一支，鳞翅目类群和双翅目类群进而形成长翅亚群总支，飞蝗中的3个LmE75亚型属于直翅目类群（图2）。

A：3个LmE75亚型结构域分析Domain structure analysis of the three LmE75 nuclear receptor isoforms；B：3个LmE75亚型A/B域核苷酸序列比对Nucleotide alignment of the A/B domain of three LmE75 isoforms；C：LmE75亚型氨基酸序列与马铃薯甲虫和德国小蠊E75多重比对Multiple alignment of LmE75 isoforms with E75 of L. decemlineata and B. germanica。在序列的上方标记DNA结合域、铰链域和配体结合域。序列下方的圆圈指示锌指中的4个Cys残基，三角形标记参与血红素结合的重要氨基酸残基The DNA binding domain, hinge domain and ligand binding domain are tagged above the sequence, the circles indicate that there are four Cys residues in the zinc finger, and the important amino acid residues involved in heme binding are marked with triangles

图2 昆虫E75系统发育分析

2.3 LmE75亚型基因组织和时期表达特性

RT-qPCR分析结果表明3个亚型具有不同的组织和发育表达特性。在体壁和肌肉中高表达；在体壁中高表达，其次是前肠；特异性在体壁高表达；在4龄和5龄均呈现先升高后降低的表达趋势，在N4D3—N4D4以及N5D4—N5D6表达量最高；在4龄和5龄的表达趋势一致，均在龄期中一天表达量最高（N4D3和N5D5），蜕皮前后表达量均较低；在4龄和5龄均呈现逐渐升高的表达趋势，在蜕皮前表达最高（图3）。

2.4 20E诱导表达及LmEcR RNAi后LmE75的表达变化

为了检测是否受20E诱导，选取N5D2的飞蝗若虫进行了20E和ds注射，RT-qPCR结果如图所示，与对照组相比，在20E诱导3 h和6 h后表达量显著上调，分别上调了4.2倍和4.5倍（图4-A）；注射ds48 h后的沉默效率为83%，的表达量亦显著被抑制，与对照组相比其表达量下调了59%（图4-B）。

2.5 基于RNAi的LmE75的生物学功能

RT-qPCR的结果显示，注射ds的飞蝗与注射ds对照组相比被显著沉默（图5-A）。对照组7 d羽化为成虫，而注射ds组7 d后飞蝗进食较对照组明显减少，活动力下降，且100%出现严重的发育延迟现象，5龄若虫期长达13—18 d，并未出现蜕皮迹象，死亡时外形无显著变化（图5-B）。体壁H＆E染色的结果显示，对照组脊线开裂时有皮层溶离现象，旧表皮降解，新表皮明显形成。ds处理组则未观察到皮层溶离，旧表皮不降解，亦无新表皮的形成（图5-C）。

3 讨论

本研究基于飞蝗转录组数据库搜索获得的3条cDNA序列，根据其结构特征分别命名为、和。在黑腹果蝇[7]、家蚕[24]、赤拟谷盗[15]、马铃薯甲虫[16]和德国小蠊[17]等昆虫中也发现具有不同数量的亚型，如马铃薯甲虫中有3个亚型，黑腹果蝇有4个亚型，德国小蠊有5个亚型。的3个亚型具有典型的核受体结构域：一个非依赖配体的激活结构域（A/B结构域）、DNA结合域（DBD，C结构域）、铰链连接域（D结构域）、配体结合域（LBD、E结构域）和F结构域。飞蝗、和中，A/B结构域分别由40、54和51个氨基酸组成，A/B结构域在核受体成员间高度可变，表明其来源于同一个基因的不同剪切方式。飞蝗的3个亚型与其他昆虫的DBD和LBD结构域具有高度一致度，DBD结构域中包含两个锌指区域，有研究报道其中的4个半胱氨酸残基可参与单个锌原子与激素响应原件DNA序列上的结合[25]，也有一些昆虫如马铃薯甲虫中只有一个锌指区域。果蝇E75是血红素传感器，位于LBD结构域的血红素中心能与双原子分子（CO或NO）结合[26]，并且含有4个参与血红素结合的重要氨基酸残基，其中第1和第3个半胱氨酸残基以及第4个组氨酸在飞蝗中是保守的。

本文采用RT-qPCR技术研究了3个亚型在5龄飞蝗不同组织中的表达特性，结果表明3个亚型具有不同的组织表达特性。在体壁和肌肉中高表达；在体壁中高表达，其次是前肠；特异性在体壁高表达；在体壁中特异性高表达，推测该基因可能参与了飞蝗表皮的发育过程。的多组织分布特性在家蚕[24]、黑腹果蝇[27]、烟草天蛾[28]和马铃薯甲虫[16]中也有报道，但不同昆虫亚型的表达特性有所不同，如马铃薯甲虫中3个亚型均在体壁低表达，而在马氏管、中肠、腹神经节、脑-心侧体-咽侧体复合体等多个组织高表达[16]。时期表达分析发现，3个亚型在飞蝗若虫期的表达量并不一致，和在4龄和5龄均为蜕皮前后表达量均较低，但表达量高于；在4龄和5龄蜕皮前表达量最高。家蚕中3个亚型也表现出不同的发育表达模式[24]，而在马铃薯甲虫中，3个亚型具有相似的发育表达特性[16]，表明不同昆虫转录本在各龄期可能发挥着不同的功能，从而共同完成的功能。在5龄若虫中期表达量较高，这一表达模式与20E滴度[22]正相关，推测的表达可能受20E的调控。事实上，体内注射20E诱导3 h和6 h后均显著上调，沉默20E受体基因后表达显著下调。类似结果在其他昆虫中也有报道，如家蚕中3个亚型均在mRNA和蛋白水平表现出对20E的响应[24]。马铃薯甲虫中，体内注射20E和体外中肠培养均表明20E可直接诱导和表达[16]。德国小蠊中，解剖脂肪体进行体外组织培养，加入20E后可诱导的4个亚型（、、）表达[17]。

图4 20E诱导及干扰LmEcR后LmE75的表达

OC：旧表皮old cuticle；NC：新表皮new cuticle；EC：皮细胞层epidermal cell；NA：无皮层溶离 No apolysis

目前仅在黑腹果蝇[7]德国小蠊[17]和褐飞虱[29]中报道了亚型的功能，结果发现亚型功能存在冗余现象，如德国小蠊不同转录本RNAi均能下调各自靶基因的表达量，但其他转录本上调，最终没有可见表型出现。因此本文选择3个亚型共同区域设计dsRNA，以全部沉默3个转录本进一步研究在飞蝗中的生物学功能，ds注射至5龄第1天后飞蝗出现发育延迟现象，死亡率达100%，且体壁无皮层溶离现象。类似表型也见于其他昆虫，如赤拟谷盗在幼虫阶段沉默后，一部分幼虫在静息阶段死亡，另一些幼虫在蜕皮时死亡[15]。在德国小蠊中，注射ds的6龄若虫可存活长达80 d，最终不能蜕皮为成虫而死亡[17]。马铃薯甲虫中通过调节20E和JH的滴度影响其信号传导通路，沉默后导致个体不能正常化蛹，在入土18 d后，大多数个体死亡[16]。以上研究表明功能在完全变态和不完全变态昆虫中是保守的。

蜕皮是昆虫等特有的生命现象，表皮是昆虫体壁的典型结构，脂类、蛋白质以及几丁质是其主要成分，并且昆虫蜕皮受激素调控。解析昆虫蜕皮发育关键基因的功能是RNAi技术应用于害虫防治的关键环节，目前关于飞蝗蜕皮发育相关基因的研究主要集中在以下几个方面：（1）蜕皮激素和保幼激素信号通路相关基因，如蜕皮激素受体EcR[22]、核受体HR3[30]和HR39[31]等；（2）表皮脂代谢关键基因，如表皮碳氢化合物合成的基因[32]、参与表皮脂转运的基因[33]等；（3）几丁质代谢以及排布关键基因，如几丁质合成酶1[34]、几丁质酶5[21]、几丁质脱乙酰基酶基因[35]等；（4）表皮蛋白基因，如ACP7[36]等。本文研究的E75是蜕皮激素信号通路中的关键核受体基因，基因沉默可导致飞蝗100%不蜕皮死亡，比较这些基因作用机制以及RNAi效果的差异可为筛选更合适的害虫控制靶标提供基础数据。

4 结论

飞蝗核受体有3个亚型，均具有典型的核受体结构域。3个亚型基因具有不同的组织和发育表达特性。体内20E诱导及RNAi试验表明，受介导的20E信号通路调控。基于RNAi功能分析发现，显著沉默后飞蝗发育延迟导致100%不蜕皮死亡。

[1] ARANDA A, PASCUAL A. Nuclear hormone receptors and gene expression., 2001, 81(3): 1269-1304.

[2] HELSEN C,CLAESSENS F. Looking at nuclear receptors from a new angle., 2014, 382(1): 97-106.

[3] 周树堂, 郭伟, 宋佳晟. 昆虫变态的激素与基因调控. 生物学通报, 2012, 47(9): 1-6.

ZHOU S T, GUO W, SONG J S. Insect metamorphosis hormones and gene regulation., 2012, 47(9): 1-6. (in Chinese)

[4] 何倩毓, 张原熙, 裴泽华, 李美鑫, 张旭. 保幼激素对昆虫变态发育调控的分子机制. 昆虫学报, 2017, 60(5): 594-603.

HE Q Y, ZHANG Y X, PEI Z H, LI M X, ZHANG X. Molecular mechanisms of juvenile hormone regulation on insect metamorphosis., 2017, 60(5): 594-603. (in Chinese)

[5] GUO X P, HARMON M A, LAUDET V, MANGELSDORF D J, PALMER M J. Isolation of a functional ecdysteroid receptor homologue from the ixodid tick(L.)., 1997, 27(11): 945-962.

[6] FAHRBACH S E, SMAGGHE G, VELARDE R A. Insect nuclear receptors., 2012, 57:83-106.

[7] BIALECKI M, SHILTON A, FICHTENBERG C, SEGRAVES W A, THUMMEL C S. Loss of the ecdysteroid-inducible E75A orphan nuclear receptor uncouples molting from metamorphosis in., 2002, 3(2): 209-220.

[8] PIERCEALL W E, LI C, BIRAN A, MIURA K, RAIKHEL A S, SEGRAVES W A. E75 expression in mosquito ovary and fat body suggests reiterative use of ecdysone-regulated hierarchies in development and reproduction., 1999, 150(1/2): 73-89.

[9] HIRUMA K, RIDDIFORD L M. Differential control of MHR3 promoter activity by isoforms of the ecdysone receptor and inhibitory effects of E75A and MHR3., 2004, 272(2): 510-521.

[10] JINDRA M, SEHNAL F, RIDDIFORD L M. Isolation, characterization and developmental expression of the ecdysteroid-induced E75 gene of the wax moth., 1994, 221(2): 665-675.

[11] RETNAKARAN A, MACDONALD A, TOMKINS W L, DAVIS C N, BROWNWRIGHT A J, PALLI S R. Ultrastructural effects of a non-steroidal ecdysone agonist, RH-5992, on the sixth instar larva of the spruce budworm,., 1997, 43(1): 55-68.

[12] SWEVERS L, ITO K, IATROU K. The BmE75 nuclear receptors function as dominant repressors of the nuclear receptor BmHR3A., 2002, 277(44): 41637-41644.

[13] SIAUSSAT D, BOZZOLAN F, QUEGUINER I, PORCHERON P, Debernard S. Effects of juvenile hormone on 20-hydroxyecdysone- inducible,,gene expression in imaginal wing cells oflepidoptera., 2004, 271(14): 3017-3027.

[14] PAUL R K, TAKEUCHI H, KUBO T. Expression of two ecdysteroid-regulated genes,and, in the brain and ovary of the honeybee (L.)., 2006, 23(12): 1085-1092.

[15] PARTHASARATHY R, TAN A, BAI H, PALLI S R. Transcription factor broad suppresses precocious development of adult structures during larval-pupal metamorphosis in the red flour beetle,., 2008, 125(3/4): 299-313.

[16] GUO W C, LIU X P, FU K Y, SHI J F, LÜ F G, Li G Q. Nuclear receptor ecdysone-induced protein 75 is required for larval-pupal metamorphosis in the colorado potato beetle(say)., 2016, 25(1): 44-57.

[17] MANE-PADROS D, CRUZ J, VILAPLANA L, PASCUAL N, BELLES X, MARTIN D. The nuclear hormone receptor BgE75 links molting and developmental progression in the direct-developing insect., 2008, 315(1): 147-160.

[18] YANG, Q P, LI Z, CAO J J, ZHANG S D, ZHANG H J, WU X Y,ZHANG Q W, LIU X X. Selection and assessment of reference genes for quantitative PCR normalization in migratory locust(Orthoptera: Acrididae)., 2014, 9(6):e98164.

[19] LIU X J, LI F,LI D Q, MA E B, ZHANG W Q,ZHU K Y, ZHANG J Z. Molecular and functional analysis ofUDP--acetylglucosamine pyrophosphorylases from themigratory locust,,2013, 8(8): e71970.

[20] LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCtmethod., 2001, 25(4): 402-408.

[21] LI D Q, ZHAGN J Q, WANG Y, LIU X J, MA E B, SUN Y, LI S, ZHU K Y, ZHANG J Z. Two chitinase 5 genes from: Molecular characteristics and functional differentiation., 2015, 58:46-54.

[22] LIU X J, SUN Y W, LI D Q, LI S, MA E B, ZHANG J Z. Identification ofas a 20-hydroxyecdysone response gene in the chitin biosynthesis pathway from the migratory locust,., 2018, 25(2): 211-221.

[23] SONG T Q, YANG M L, WANG Y L, LIU Q, WANG H M, ZHANG J, LI T. Cuticular protein LmTwdl1 is involved in molt development of the migratory locust., 2016, 23(4): 520-530.

[24] LI K, GUO E E, HOSSAIN M S, LI Q R, CAO Y, TIAN L, DENG X J, LI S.E75 isoforms display stage- and tissue-specific responses to 20-hydroxyecdysone.2015, 5: 12114.

[25] ESCRIVA H, BERYRAND S, LAUDET V. The evolution of the nuclear receptor superfamily., 2004, 40(2): 11-26.

[26] CACERES L, NECAKOV A S, SCHWARTZ C, KIMBER S, ROBERTS I J, KRAUSE H M. Nitric oxide coordinates metabolism, growth, and development via the nuclear receptor E75., 2011, 25(14): 1476-1485.

[27] JOHNSTON D M, SEDKOV Y, PETRUK S, RILEY K M, FUJILKA M, JAYNES J B, MAZO A. Ecdysone- and no-mediated gene regulation by competing EcR/Usp and E75A nuclear receptors duringdevelopment., 2011, 44(1): 51-61.

[28] KESHAN B, HIRUMA K, RIDDIFORD L M. Developmental expression and hormonal regulation of different isoforms of the transcription factor E75 in the tobacco hornworm., 2006, 295(2): 623-632.

[29] 李凯龙. 褐飞虱蜕皮及变态信号途径相关基因的功能分析[D]. 北京: 中国农业科学院, 2016.

LI K L. Functional analysis of the genes involved in molting and metamorphosis signal pathways in[D]. Beijing: Chinese Academy of Agricultural Sciences,2016. (in Chinese)

[30] ZHAO X M, QIN Z Y, LIU W M, LIU X J, Moussian B, MA E B, LI S, ZHANG J Z. Nuclear receptor HR3 controls locust molt by regulating chitin synthesis and degradation genes of., 2018, 92: 1-11.

[31] ZHAO X M, QIN Z Y, ZHANG J, YANG Y, JIA P, YANG Q, MA E B, ZHANG J Z. Nuclear receptor HR39 is required for locust molting by regulating the chitinase and carboxypeptidase genes., 2019, 28(4): 537-549.

[32] YU Z T, ZHANG X Y, WANG Y W, Moussian B, ZHU K Y, LI S, MA E B, ZHANG J Z.: An oenocyte-specific cytochrome P450 gene required for cuticular waterproofing in the migratory locust,., 2016, 6: 29980.

[33] YU Z T, WANG Y W, ZHAO X M, LIU X J, MA E B, Moussian B, ZHANG J Z. The ABC transporteris needed for cuticle barrier construction in.2017, 87: 90-99.

[34] ZHENG J Z, LIU X J, ZHANG J Q, LI D Q, SUN Y, GUO Y P, MA E B, ZHU K Y. Silencing of two alternative splicing-derived mrna variants of chitin synthase 1 gene by RNAi is lethal to the oriental migratory locust,(meyen)., 2010, 40(11): 824-833.

[35] YU R R, LIU W M, LI D Q, ZHAO X M, DING G W, ZHANG M, MA E B, ZHU K Y, LI S, MOUSSIAN B, ZHANG J Z. Helicoidal organization of chitin in the cuticle of the migratory locust requires the function of the chitin deacetylase 2 enzyme (LmCDA2)., 2016, 291(47): 24352-24363.

[36] ZHAO X M, GOU X, LIU W M, MA E B, MOUSSIAN B, LI S, ZHU K Y, ZHANG J Z. The wing-specific cuticular protein LmACP7 is essential for normal wing morphogenesis in the migratory locust., 2019, 112: 103206.

Molecular Characteristics and Function Analysis of Nuclear Receptor Gene

LIU XiaoJian1, Guo Jun1,2, Zhang Xueyao1, MA Enbo1, ZHANG Jianzhen1

(1Research Institute of Applied Biology, Shanxi University, Taiyuan 030006;2College of Life Science, Shanxi University, Taiyuan 030006)

【】Nuclear receptors (NRs) act as ligand-inducible transcription factors, which play an important role during growth and development in insects.is an important agricultural pest. The objective of this study is to analyze the molecular characteristics and biological function of the nuclear receptor gene, and to provide a new molecular target for pest control.【】The cDNA sequences of threegenes were obtained based on the transcriptome database of. The open reading frames, predicted amino acid sequences, physicochemical properties and conserved domains of each gene were analyzed using Blast and ExPASy websites. Using MEGA 6.0 software, the neighbor-joining (NJ) method was used to construct a phylogenetic tree with the homologous sequences of E75 from other insects. Specific primers of threeisoforms were designed using Primer 3.0 software. Reverse transcription-quantitative PCR (RT-qPCR) was used to analyze the expression characteristics of threeisoforms in different tissues and developmental days from 4th to 5th instar nymphs. To further test whetherwas regulated by 20-hydroxyecdysone (20E), the expression ofwas detected by RT-qPCR after injected with 20Eandinterfered with the 20E receptor geneby RNAi. The double-stranded RNA (dsRNA) of the common region of threeisoforms was synthesized,and then was injected into the 5th instar nymphs. The total RNA of integument was extracted after 48 h, and the first strand cDNA was synthesized as the template of RT-qPCR. The silencing efficiency ofwas detected by RT-qPCR technology and the phenotype was observed. The effect on the structure of cuticle after injected with dswas analyzed by H&E staining method.【】Three cDNA sequences ofwere obtained based on the transcriptome database of. According to structural characteristics, they were named,and, respectively. The GenBank accession numbers were MN584732, MN584733 and MN584734. Each of the three LmE75 isoforms had typical nuclear receptor conserved domains, and DBD (DNA binding domain) and LBD (ligand binding domain) were highly conserved among insects. Phylogenetic tree analysis showed that E75 isoforms of different orders of insects were preferentially clustered into one group. RT-qPCR results showed that threeisoforms had different tissue and developmental expression characteristics.was highly expressed in the integument and muscle.was highly expressed in the integument, followed by the foregut.was specifically expressed in the integument. The expression level ofincreased first and then decreased at the 4th and 5th instar nymphs. The expression levels ofwere the highest in N4D3 and N5D5, and lower before and after molting. The expression level ofincreased gradually at both the 4th and 5th instar nymphs, with the highest expression before molting. Compared with the control group, the expression level ofsignificantly increased by 4.2 and 4.5 times after 20E induction for 3 h and 6 h, respectively. The expression level ofwas significantly down-regulated by 59% after interference withfor 48 h by RNAi. After injection of dson the 1st day of 5th instar nymphs,exhibited obvious developmental delay phenotype and allcould not molt to adults. The H&E staining of the integument found that the cuticle of ds-injectedshowed no apolysis.【】There are threeisoforms in, and they have different tissue and developmental expression characteristics. RNAi results show thatplay an important role in the molting of

; nuclear receptor (NR);; 20-hydroxyecdysone (20E); RNA interference (RNAi)