黄曲条跳甲三个Cry3Aa杀虫相关基因的克隆及生物信息学分析

I引言黄曲条跳甲（Phyllotretastriolata）作为全球性分布在十字花科蔬菜范围内的害虫，在我国南方地区危害尤为严重，其成虫和幼虫分别危害叶片和根部，严重影响油菜、萝卜、白菜等作物的产量REF_Ref2338\r\h[1]。目前，该害虫的防治主要依赖化学农药，但长期使用导致其对辛硫磷、氯氰菊酯等药剂产生抗性，亟需开发新型高效、环境友好的防控策略REF_Ref29682\r\h[2]。Cry3Aa毒素作为苏云金芽孢杆菌（Bacillusthuringiensis）分泌的杀虫蛋白，对鞘翅目害虫具有特异性毒杀作用，但其作用机制及靶标受体在黄曲条跳甲中的分子基础尚不明确REF_Ref31282\r\h[3]REF_Ref31289\r\h[4]。Cry3Aa毒素的杀虫活性依赖于其与害虫中肠受体的结合，进而破坏细胞膜完整性，导致细胞裂解和害虫死亡REF_Ref31364\r\h[5]。已有研究表明，ABCB转运蛋白（如P-糖蛋白）作为ATP结合盒转运蛋白超家族成员，可通过外排毒素降低其胞内积累。叶甲科害虫的ABCB蛋白（如P-glycoprotein）在Cry3Aa抗性中起关键作用REF_Ref31429\r\h[6]；ADAM金属蛋白酶作为一类跨膜蛋白，其去整合素和金属蛋白酶结构域可能参与Cry3Aa的蛋白水解活化REF_Ref31547\r\h[7]；配体印迹实验证实，马铃薯甲虫（Leptinotarsadecemlineata）的ADAM同源蛋白能与Cry3Aa特异性结合REF_Ref31903\r\h[8]。碱性磷酸酶（ALP）作为GPI锚定蛋白，ALP可介导Cry3Aa与肠上皮细胞的初始结合，黄粉虫（Tenebriomolitor）中肠ALP被证实是Cry3Aa的功能性受体REF_Ref32014\r\h[9]。然而，上述研究多集中于其他鞘翅目害虫，黄曲条跳甲中与Cry3Aa互作的受体基因尚未被系统研究，其分子特征及功能机制仍不明确。本研究拟克隆黄曲条跳甲中与Cry3Aa毒力相关的ABCB、ADAM和ALP基因，并对其进行生物信息学分析，以揭示其结构特征及潜在功能。研究结果将有助于阐明Cry3Aa在黄曲条跳甲中的分子作用机制，为基于受体靶点的Bt生物农药优化及抗性治理提供理论依据，同时为开发新型靶向防控技术奠定分子基础。材料与方法2.1实验材料实验采用采自沙坪坝地区的黄曲条跳甲成虫，于实验室条件下饲养后，选取健康个体并按雌雄分组。为减少个体差异，每6头成虫随机混合为1个样本，每组处理设置3个独立的生物学重复（即3个混合样本），所有样本独立处理与分析。2.2主要试剂与仪器实验中所用的主要试剂包括：RNA提取使用Trizol试剂（Invitrogen）、氯仿、异丙醇和无水乙醇等；cDNA合成采用TaKaRa公司生产的PrimerScript™RT试剂盒；PCR扩增中使用普通Taq酶；克隆载体为pClone007（北京擎科）。主要实验仪器包括：组织研磨机（SCIENTZ-48）、冷冻离心机（Eppendorf5424R）、PCR仪（Bio-RadT100）及电泳系统。2.3实验方法2.3.1RNA提取与cDNA合成RNA提取：取6头黄曲条跳甲成虫与200μLTrizol试剂充分研磨，加入氯仿并进行分层离心，收集上清液。随后，利用异丙醇沉淀RNA，并使用75%乙醇进行洗涤，最后用DEPC水溶解沉淀的RNA。通过分光光度计测定RNA浓度及纯度，并利用1.5%琼脂糖凝胶电泳验证RNA的完整性。cDNA合成：使用TaKaRa试剂盒去除基因组DNA后，利用RNA作为模板进行反转录反应以合成cDNA，合成的cDNA产物保存于-80℃冷冻条件下以备后续实验使用REF_Ref32151\r\h[10]。2.3.2引物设计与PCR扩增根据鞘翅目赤拟谷盗、马铃薯甲虫等同源基因序列，从黄曲条跳甲的转录组中筛选出目标基因（包括ABCB、ADAM和ALP）。利用Primer5.0软件和NCBIPrimer-BLAST工具设计特异性引物（表1）。PCR反应体系构建为25μL，其中包含TaqMix12.5μL、上下游引物各1μL、cDNA模板2μL、ddH2O8.5μL。PCR反应程序设定为：94℃预变性3分钟；随后进行94℃变性30秒、50-60℃（根据引物Tm值调节）退火30秒和72℃延伸30秒至1分钟（根据片段长度调整），共进行35个循环；最后在72℃进行延伸10分钟。PCR产物通过1.2%琼脂糖凝胶电泳进行检测REF_Ref32220\r\h[11]。表1引物序列信息Table1Primersandtheirsequencesinformation引物名称Primer引物序列Primersequence(5'-3')ABCB-FTCGTTACTGAGTTTGCABCB-RACGGGACTGCTGACATADAM-FTGACTAACCACAATACAGTGADAM-RAACCCCTACAAGTACAAGATALP-FTCGCTCGCCATACCTAACTTALP-RTAATCCAAGCCTGGTCTAAATAA2.3.3克隆与测序PCR扩增的产物经过切胶纯化后，被连接到pClone007载体，随后又被转入Trelief-5a感受态细胞，含Amp抗性的LB平板上对其进行了培养操作，从平板挑出单菌落后进一步扩大培养，接着送去北京擎科公司安排测序，通用引物M13用来对插入片段加以验证确认。2.3.4生物信息学分析本研究对克隆获得的ABCB、ADAM和ALP基因进行生物信息学分析。首先，利用DNAman和DNAssist软件完成基因序列比对、开放阅读框（ORF）识别及基本序列特征分析。通过ExPASy平台在线工具预测蛋白理化性质与结构特征：ProtParam（/protparam/）计算分子量、等电点等参数；ProtScale（/protscale/）分析亲水性；TMHMM（http://www.cbs.dtu.dk/services/TMHMM/）预测跨膜结构域；SignalP（http://www.cbs.dtu.dk/services/SignalP）检测信号肽。基于氨基酸全长序列，采用MEGA7软件构建系统发育树：通过邻接法（Neighbor-Joining）生成进化关系拓扑结构，并以1000次自展重复评估分支置信度。此外，使用Pfam数据库（/）对蛋白功能结构域进行注释，明确各基因的潜在生物学功能REF_Ref32314\r\h[12]。结果与分析3.1ABCB基因的克隆与生物信息学分析3.1.1基因克隆与序列特征通过PCR扩增之后再经测序验证，黄曲条跳甲ABCB基因的完整序列便被成功取得，这一基因的开放阅读框（ORF）可达3843个碱基对，编码1280个氨基酸，从序列比对的结果来看，这个基因同鞘翅目昆虫玉米根萤叶甲（Diabroticavirgiferavirgifera）的ABCB基因达到了82%的同源程度，暗示出该基因在进化过程中具备比较高的保守属性。3.1.2蛋白质理化性质分析ABCB1蛋白的分子量为140.99kDa，等电点（pI）为8.80，表明其在碱性环境中更稳定。氨基酸组成分析显示，该蛋白富含亮氨酸（Leu）、丙氨酸（Ala）和甘氨酸（Gly），疏水性氨基酸占比较高（GRAVY指数：0.071），符合膜蛋白的特征。此外，该蛋白的脂肪族指数（Aliphaticindex）为102.12，表明其结构中脂肪族侧链（如Ala、Val、Ile、Leu）含量较高，这一特性通常与蛋白质的热稳定性及跨膜结构域的疏水环境相容性相关。结合不稳定系数（32.93，<40），说明ABCB1在体内外均具有较高的结构稳定性，适合作为跨膜转运蛋白长期行使功能。3.1.3跨膜结构域与信号肽预测跨膜结构预测显示，ABCB1蛋白含有12个明显的跨膜螺旋，其胞内和胞外区域交替分布（图3-1），后验概率接近1，表明预测结果高度可靠。这种典型的跨膜拓扑结构进一步证实了ABCB1作为膜转运蛋白的功能特性。信号肽预测结果表明，该蛋白不含任何已知类型的信号肽（Sec/SPI概率：0.0001）（图3-2），提示其可能通过非经典途径定位于细胞膜，或作为固有膜蛋白发挥作用。图3-1ABCB1蛋白跨膜结构预测Figure3-1PredictedTransmembraneStructureofABCB1Protein图3-2ABCB1蛋白信号肽预测Figure3-2PredictedSignalPeptideofABCB1Protein3.1.4系统进化分析基于ABCB基因构建的系统发育树显示（图3-3），黄曲条跳甲与鞘翅目昆虫玉米根萤叶甲（Diabroticavirgiferavirgifera）、杨树叶甲（Chrysomelatremula）聚为一支，Bootstrap支持率>80%，表明它们具有共同的进化起源。而与鳞翅目（如家蚕Bombyxmori）和双翅目（如果蝇Drosophilamelanogaster）的ABCB基因分化明显，反映了该基因在昆虫不同目中的功能特异性。图3-3基于ABCB基因的系统发育树Figure3-3PhylogeneticTreeBasedonABCBGene3.1.5功能结构域与生物学意义保守结构域分析发现，ABCB1蛋白含有典型的ABC转运蛋白ATP结合结构域（Pfam:PF00005），这一特征与其在药物外排和毒素抗性中的功能一致。结合已有研究，推测黄曲条跳甲ABCB1蛋白可能通过ATP水解提供能量，将Cry3Aa毒素泵出细胞外，从而降低其毒性作用REF_Ref31289\r\h[4]。这一机制可能与该害虫对Bt毒素的抗性演化密切相关。3.2ADAM基因的克隆与生物信息学分析3.2.1基因克隆与序列特征通过RT-PCR扩增技术还有测序工作的完成，顺利取得了黄曲条跳甲ADAM基因的全长cDNA序列，经过分析发现，这个基因开放阅读框（ORF）有长达3372bp的情况，能编码出1124个氨基酸，进一步做序列对比之后表明，这种基因与鞘翅目昆虫当中的科罗拉多马铃薯甲虫（Leptinotarsadecemlineata）ADAM基因同源性达到了78%左右，这种现象在鞘翅目昆虫群体当中，这基因存在较高的保守程度显现了出来。3.2.2蛋白质理化性质分析ADAM蛋白的分子量为130.27kDa，等电点（pI）为6.69，表明其在近中性环境中较为稳定。氨基酸组成分析显示，该蛋白富含亮氨酸（Leu，12.7%）、谷氨酸（Glu，6.9%）和丝氨酸（Ser，7.0%），具有较多的负电荷残基（138个）和正电荷残基（132个）。该蛋白的不稳定系数为42.87，属于不稳定型蛋白（>40），提示其在体内可能需要严格的调控机制来维持稳定性。脂肪族指数（Aliphaticindex）为98.43，表明其含有较高比例的脂肪族氨基酸（如Leu、Ile、Val），结合GRAVY值（-0.179）显示该蛋白整体呈微弱亲水性，但可能存在局部疏水结构域，这一特性与其作为膜锚定蛋白的功能特征相符。3.2.3跨膜结构域与信号肽预测DeepTMHMM跨膜区预测结果显示（图3-4），ADAM蛋白不含任何跨膜结构域，所有氨基酸残基均位于胞内（后验概率=1）。信号肽预测分析表明（图3-5），该蛋白不含任何已知类型的信号肽（Sec/SPI、Sec/SPII等概率均为0），综合预测结果为"其他"。这些结果共同表明ADAM蛋白是一个典型的胞内蛋白，不参与跨膜转运或分泌过程。图3-4ADAM蛋白跨膜区预测Figure3-4PredictedTransmembraneRegionofADAMProtein图3-5ADAM蛋白信号肽预测Figure3-5PredictedSignalPeptideofADAMProtein3.2.4系统进化分析基于ADAM基因构建的系统发育树显示（图3-6），黄曲条跳甲与鞘翅目昆虫科罗拉多马铃薯甲虫（Leptinotarsadecemlineata）、赤拟谷盗（Triboliumcastaneum）聚为一支，Bootstrap支持率>50%，证实了它们之间较近的亲缘关系。而与甘蔗黑蚜（Melanaphissacchari，半翅目）与次生隐白蚁（Cryptotermessecundus，蜚蠊目）的分支相距较远，反映了ADAM基因在跨目水平上的功能分化。图3-6基于ADAM基因的系统发育树Figure3-5PhylogeneticTreeBasedonADAMGene3.2.5功能结构域与生物学意义SMART数据库分析显示，ADAM蛋白含有典型的金属蛋白酶结构域（PF01421）和解整联蛋白结构域（PF00200）。结合已有研究，推测该蛋白可能通过其金属蛋白酶活性参与黄曲条跳甲中肠细胞的细胞外基质重塑，进而影响Cry3Aa毒素与中肠受体的结合过程REF_Ref32647\r\h[13]。虽然预测显示其为胞内蛋白，但其较高的不稳定性（42.87）提示可能在特定生理条件下发生转位或构象变化，从而参与毒素响应过程。3.3ALP基因的克隆与生物信息学分析3.3.1基因克隆与序列特征通过RT-PCR扩增及测序验证，成功获得黄曲条跳甲ALP基因的全长cDNA序列。该基因开放阅读框（ORF）长度为1617bp，编码539个氨基酸。序列比对分析显示，ALP基因与鞘翅目昆虫光肩星天牛（Anoplophoraglabripennis）的ALP基因同源性最高（75%），表明其在鞘翅目昆虫中具有较高的功能保守性。3.3.2蛋白质理化性质分析ALP蛋白的分子量为60.30kDa，等电点（pI）为5.82，表明其在弱酸性环境中更稳定。氨基酸组成分析显示，该蛋白富含丙氨酸（Ala，9.1%）、谷氨酸（Glu，7.6%）和亮氨酸（Leu，8.2%）。脂肪族指数（Aliphaticindex）为75.14，显著低于典型膜蛋白（通常>90），结合GRAVY值（-0.467）表明该蛋白具有显著亲水性特征，这一特性与其可能作为可溶性蛋白的功能定位相符。该蛋白的不稳定系数为41.65（>40），属于不稳定型蛋白，提示其在细胞内可能通过动态调控参与特定生理过程，这种结构不稳定性与其较低的脂肪族氨基酸含量可能相关。3.3.3跨膜结构域与信号肽预测DeepTMHMM跨膜区预测结果显示（图3-7），ALP蛋白不含典型跨膜结构域，其N端前20个氨基酸残基被预测为信号肽区域（Signal概率>0.8），而其余序列均位于胞内（Inside概率=1）。然而，SignalP6.0分析表明（图3-8），该区域未达到经典信号肽阈值（所有类型概率均为0），综合预测结果为"其他"。结合两种工具结果，推测ALP蛋白可能通过非经典途径短暂定位于细胞膜，但主要发挥胞内功能。图3-7ALP蛋白跨膜区预测Figure3-7PredictedTransmembraneRegionofALPProtein图3-8ALP蛋白信号肽预测Figure3-8PredictedSignalPeptideofALPProtein3.3.4系统进化分析基于ALP基因构建的系统发育树显示（图3-9），黄曲条跳甲与鞘翅目昆虫光肩星天牛（Anoplophoraglabripennis）、科罗拉多马铃薯甲虫（Leptinotarsadecemlineata）聚为一支，Bootstrap支持率>98%，证实其亲缘关系密切。而与鳞翅目昆虫棉铃虫（Helicoverpaarmigera）、家蚕（Bombyxmori）分化明显，反映了ALP基因在不同目昆虫中的功能特异性。图3-9基于ALP基因的系统发育树Figure3-9PhylogeneticTreeBasedonALPGene3.3.5功能结构域与生物学意义SMART数据库分析显示，ALP蛋白含碱性磷酸酶催化结构域（PF00245），这一特征与已报道的Cry3Aa受体蛋白功能一致REF_Ref32014\r\h[9]。结合跨膜预测结果，推测该蛋白可能通过GPI锚定短暂附着于中肠细胞膜表面，介导Cry3Aa毒素的结合。其较高的不稳定性（41.65）提示该结合过程可能具有动态调控特性，这可能与黄曲条跳甲对毒素的适应性响应相关。讨论在本研究中，我们成功克隆了黄曲条跳甲的三个与Cry3Aa毒素抗性相关的基因——ABCB、ADAM和ALP，并进行了生物信息学分析。通过这些分析，我们揭示了这些基因的序列特征、蛋白质性质以及它们在黄曲条跳甲抗性机制中的潜在作用。ABCB基因在本研究中表现出较高的保守性，其与玉米根萤叶甲的同源性高达82%，表明该基因在不同昆虫物种间具有较强的进化保守性。ABCB1蛋白属于ABC转运蛋白家族，已知在药物外排和毒素抗性中扮演重要角色REF_Ref31289\r\h[4]。在黄曲条跳甲中，ABCB1蛋白的跨膜结构特征和ATP结合域的存在使其可能通过能量驱动的方式将Cry3Aa毒素从细胞内泵出，从而降低毒素的细胞毒性REF_Ref121\r\h[14]。结合已有文献，推测ABCB1蛋白可能参与黄曲条跳甲对Bt毒素的抗性演化过程。此外，ABCB1蛋白的氨基酸组成分析显示其疏水性氨基酸比例较高，这符合其作为膜转运蛋白的功能特性。其较高的稳定性使其能够长期在细胞膜上发挥作用，这一特性为黄曲条跳甲在面对Cry3Aa毒素时的生存提供了重要支持。ADAM基因编码的蛋白在本研究中显示出与科罗拉多马铃薯甲虫等其他鞘翅目昆虫的较高同源性（78%），进一步证实了其在鞘翅目昆虫中具有保守性。与ABCB基因不同，ADAM蛋白被预测为典型的胞内蛋白，其不含跨膜结构域，且主要功能可能与细胞内基质的重塑和毒素的响应有关。基于其金属蛋白酶和解整联蛋白结构域，ADAM蛋白可能通过其金属蛋白酶活性参与黄曲条跳甲肠细胞中细胞外基质的调节。这种调节作用可能影响Cry3Aa毒素与中肠受体的结合，从而降低毒素的生物活性REF_Ref32647\r\h[13]。虽然ADAM蛋白