群居型与散居型飞蝗在真菌感染前后可变剪切的差异及机制探究

群居型与散居型飞蝗在真菌感染前后可变剪切的差异及机制探究一、引言1.1研究背景飞蝗（Locustamigratoria）作为世界上分布最为广泛的蝗虫种类，长期以来一直是农业领域的重大威胁。其分布范围覆盖了东半球的温带和热带地区，包括非洲、亚洲、中东以及澳大利亚等区域，给当地的农业生产带来了严重的破坏。在我国，飞蝗主要包括东亚飞蝗（Locustamigratoriamanilensis）、亚洲飞蝗（Locustamigratoriamigratoria）和西藏飞蝗（LocustamigratoriatibetnsisChen）三个亚种，其中东亚飞蝗广泛分布于东部季风区，是我国分布最广、危害最重的历史性害虫。飞蝗具有群居型和散居型两种生活型，这两种生活型在形态、生理和行为等方面均存在显著差异。群居型飞蝗通常体色较深，如灰黄褐色，它们具有较强的聚集性和迁飞能力，能够远距离飞行，形成庞大的群体，对农作物造成毁灭性的危害。当群居型飞蝗大规模爆发时，所到之处农作物被迅速啃食，常常导致颗粒无收，严重影响农业经济和粮食安全。相比之下，散居型飞蝗体色较浅，头、胸、后足常带绿色，它们较为分散，活动范围相对较小，对农作物的危害程度相对较轻。然而，在适宜的环境条件下，散居型飞蝗可以转变为群居型，从而引发蝗灾。这种生活型的转变与多种因素有关，包括种群密度、食物资源、环境温度等。例如，当种群密度增加时，飞蝗个体之间的相互接触和信息交流增多，会促使它们逐渐聚集并表现出群居行为。在农业生产中，蝗虫灾害与水、旱灾害常交替发生，给农业生产带来了极大的挑战。据统计，全世界常年发生蝗虫的面积达4680万km²，全球约1/8的人口经常受到蝗灾的袭扰。在我国，历史上就有多次飞蝗大爆发的记载，如《史记》中曾记载公元前243年发生的一次蝗灾：“蝗虫从东方来，蔽天。”这些蝗灾不仅对当时的农业生产造成了巨大破坏，还引发了饥荒等一系列社会问题。传统的蝗虫防治方法主要依赖化学杀虫剂，然而，化学杀虫剂的大量使用带来了诸多问题。一方面，化学杀虫剂会对环境造成严重污染，破坏生态平衡，影响非靶标生物的生存和繁衍；另一方面，长期使用化学杀虫剂还会导致蝗虫产生抗药性，使得防治效果逐渐下降。此外，化学杀虫剂的残留还可能对食品安全构成威胁，影响人类健康。利用昆虫病原真菌来防治蝗虫是一种绿色、可持续的生物防治方法，近年来受到了广泛关注。昆虫病原真菌能够通过感染蝗虫，在其体内生长繁殖，最终导致蝗虫死亡。与化学防治相比，真菌防治具有环境友好、不易产生抗药性等优点。例如，绿僵菌（Metarhiziumanisopliae）已进化为一种专门的蝗虫病原体，能够有效地控制蝗虫种群数量。然而，蝗虫在长期的进化过程中也形成了一系列防御机制来应对真菌的侵染。研究蝗虫在真菌感染前后的免疫应答机制，对于深入了解蝗虫与真菌之间的相互作用关系，提高真菌防治蝗虫的效果具有重要意义。可变剪切是一种重要的基因表达调控机制，在昆虫的生长、发育和免疫等过程中发挥着关键作用。通过可变剪切，一个基因可以产生多种不同的mRNA转录本，进而翻译出不同的蛋白质异构体，增加了蛋白质组的复杂性和生物功能的多样性。在昆虫免疫方面，可变剪切参与了免疫信号通路的调控、免疫相关基因的表达以及抗菌肽的合成等过程。例如，在果蝇中，一些免疫相关基因的可变剪切能够影响其对病原体的抗性。然而，目前关于群居型和散居型飞蝗在真菌感染前后可变剪切的差异研究还相对较少。深入探究这两种生活型飞蝗在真菌感染前后的可变剪切变化，有助于揭示飞蝗免疫应答的分子机制，为开发更加有效的蝗虫生物防治策略提供理论依据。1.2研究目的与意义本研究旨在深入探究群居型和散居型飞蝗在真菌感染前后可变剪切的差异，揭示可变剪切在飞蝗免疫应答过程中的作用机制。通过对两种生活型飞蝗在真菌感染前后的转录组进行测序和分析，全面鉴定可变剪切事件，筛选出与免疫相关的差异可变剪切基因，并进一步验证这些基因的功能。这不仅有助于我们从分子层面深入理解飞蝗的免疫防御机制，还能为开发高效、绿色的蝗虫生物防治技术提供理论支持。蝗虫灾害长期以来严重威胁全球农业生产，给粮食安全带来了巨大挑战。据统计，全球每年因蝗虫灾害造成的农业损失高达数十亿美元。传统化学防治方法虽然在一定程度上能够控制蝗虫数量，但带来的环境污染、抗药性产生等问题日益突出，迫切需要寻找更加环保、可持续的生物防治策略。昆虫病原真菌作为一种天然的生物防治剂，具有环境友好、特异性强等优点，然而蝗虫对真菌侵染的免疫防御机制尚不完全清楚，限制了真菌防治技术的进一步发展。可变剪切作为基因表达调控的重要方式，在昆虫免疫过程中发挥着关键作用。研究群居型和散居型飞蝗在真菌感染前后的可变剪切差异，有助于揭示飞蝗免疫应答的分子调控网络，为开发针对蝗虫免疫弱点的生物防治策略提供新的思路和靶点。例如，通过干扰关键免疫基因的可变剪切，可能削弱蝗虫的免疫防御能力，从而增强真菌对蝗虫的感染效果。此外，深入了解可变剪切在飞蝗免疫中的作用，还能为其他害虫的生物防治研究提供借鉴，推动整个农业害虫防治领域的发展，对于保障全球农业生产和生态平衡具有重要的现实意义。1.3研究现状1.3.1飞蝗的生物学特性研究飞蝗作为一种典型的具有群居和散居两种生活型的昆虫，其生物学特性一直是研究的热点。在形态特征方面，群居型飞蝗通常体色较深，体型相对较大，具有发达的飞行肌肉和较长的翅膀，便于远距离迁飞；而散居型飞蝗体色较浅，体型稍小，飞行能力相对较弱。例如，东亚飞蝗的群居型成虫体色多为灰黄褐色，前胸背板中隆线两侧的黑色带纹明显，前翅较长，超过后足股节2倍以上；散居型成虫头、胸、后足常带绿色，前胸背板中隆线两侧的黑色带纹不明显，前翅较短，不到后足股节的2倍。飞蝗的生活习性也与其生活型密切相关。群居型飞蝗具有高度的聚集性，它们常常形成庞大的群体进行迁飞和取食，对农作物造成严重的危害。当群居型飞蝗大规模迁飞时，遮天蔽日，所到之处农作物被迅速啃食殆尽。散居型飞蝗则较为分散，单独活动，对农作物的危害相对较小。研究表明，飞蝗生活型的转变受到多种因素的影响，其中种群密度是关键因素之一。当种群密度较低时，飞蝗表现为散居型；随着种群密度的增加，飞蝗会逐渐转变为群居型。此外，食物资源、环境温度、光照等因素也会对飞蝗生活型的转变产生影响。在繁殖特性上，飞蝗的繁殖能力较强，产卵量较大。以东亚飞蝗为例，雌虫一般产4-5个卵块，每个卵块有卵约50-75粒。飞蝗的卵通常产在土壤中，形成卵囊，卵囊具有一定的保护结构，能够抵御外界环境的不利影响，保证卵的正常孵化。1.3.2昆虫病原真菌防治蝗虫的研究利用昆虫病原真菌防治蝗虫是一种重要的生物防治方法，具有广阔的应用前景。目前，研究较多的用于防治蝗虫的昆虫病原真菌主要有绿僵菌、白僵菌等。绿僵菌能够通过体表入侵蝗虫体内，在蝗虫体内生长繁殖，消耗蝗虫的营养物质，最终导致蝗虫死亡。研究发现，绿僵菌对不同生活型的飞蝗都具有一定的致病性，但致病性可能存在差异。昆虫病原真菌防治蝗虫的效果受到多种因素的影响。首先，真菌的种类和菌株特性对防治效果有重要影响。不同种类的真菌以及同一真菌的不同菌株，其致病力、生长速度、抗逆性等方面存在差异。其次，环境因素如温度、湿度、光照等也会影响真菌的生长和侵染能力。例如，绿僵菌在适宜的温度（25-30℃）和湿度（80%-90%）条件下，能够更好地生长和侵染蝗虫。此外，蝗虫的生理状态、免疫防御能力等也会对真菌的防治效果产生影响。在实际应用中，昆虫病原真菌防治蝗虫还面临一些挑战。例如，真菌的大规模生产技术还不够成熟，生产成本较高；真菌制剂的稳定性和储存条件要求较高，在实际应用中容易受到环境因素的影响而降低防治效果。因此，需要进一步加强对昆虫病原真菌的研究，优化生产工艺，提高真菌制剂的质量和稳定性，以推动其在蝗虫防治中的广泛应用。1.3.3可变剪切在昆虫领域的研究进展可变剪切是一种重要的基因表达调控机制，在昆虫的生长、发育、免疫等过程中发挥着关键作用。在昆虫生长发育方面，可变剪切参与了昆虫蜕皮、变态等重要过程的调控。例如，在果蝇中，一些基因的可变剪切能够影响其蜕皮激素信号通路，从而调控果蝇的蜕皮和变态发育。在蚊子中，特定基因的可变剪切与蚊子的生殖发育密切相关。在昆虫免疫领域，可变剪切同样发挥着重要作用。研究发现，昆虫在受到病原体侵染时，免疫相关基因会发生可变剪切，产生不同的转录本和蛋白质异构体，这些异构体可能具有不同的免疫功能。例如，在果蝇中，Toll信号通路中的一些基因在受到真菌感染时会发生可变剪切，产生的不同异构体能够调节Toll信号通路的激活程度，从而影响果蝇对真菌的免疫应答。在蜜蜂中，免疫相关基因的可变剪切也参与了蜜蜂对病毒和细菌感染的免疫防御过程。目前，关于可变剪切在昆虫中的研究主要集中在模式昆虫如果蝇、蚊子等，对于飞蝗等农业害虫的可变剪切研究相对较少。特别是群居型和散居型飞蝗在真菌感染前后可变剪切的差异研究还处于起步阶段，深入探究这一领域，有助于揭示飞蝗免疫应答的分子机制，为蝗虫的生物防治提供新的理论依据。二、材料与方法2.1实验材料2.1.1飞蝗样本采集本研究选取东亚飞蝗作为实验对象，分别采集群居型和散居型飞蝗样本。群居型飞蝗样本采集于山东省东营市黄河三角洲地区，该地区是东亚飞蝗的重要发生地之一，具有典型的蝗虫栖息地特征，如广阔的滩涂、丰富的芦苇等植物资源，为飞蝗的生存和繁殖提供了适宜的环境。散居型飞蝗样本采集于河南省新乡市的农田周边，这些区域农作物种植较为分散，飞蝗密度相对较低，符合散居型飞蝗的生活环境特点。在采集飞蝗样本时，使用捕虫网进行捕捉。为确保样本的代表性，每个类型的飞蝗随机采集50只，涵盖不同龄期和性别。采集后的飞蝗样本立即装入透气的昆虫饲养盒中，并标记采集地点、时间和类型等信息。将采集到的飞蝗带回实验室后，分别在不同的养殖箱中进行养殖。养殖箱尺寸为长50cm×宽40cm×高30cm，底部铺设一层5cm厚的湿润土壤，以模拟飞蝗的自然栖息环境。在养殖箱内放置足量的新鲜小麦苗作为食物，小麦苗每日更换，以保证飞蝗能获取充足的营养。养殖箱放置在温度为28±2℃、相对湿度为60±5%、光照周期为14h光照/10h黑暗的人工气候箱中，为飞蝗提供适宜的生长环境。在养殖过程中，定期观察飞蝗的生长状态，及时清理养殖箱中的粪便和剩余食物，确保养殖环境的清洁。2.1.2真菌菌株选择与培养选用金龟子绿僵菌（Metarhiziumanisopliae）作为感染飞蝗的真菌菌株，该菌株购自中国农业微生物菌种保藏管理中心，编号为ACCC30533。金龟子绿僵菌是一种广泛应用于害虫生物防治的昆虫病原真菌，对蝗虫具有较强的致病性。将金龟子绿僵菌接种到马铃薯葡萄糖琼脂（PDA）培养基上，培养基配方为：马铃薯200g、葡萄糖20g、琼脂15g、水1000mL。先将马铃薯去皮切块，煮沸20min后过滤取汁，然后加入葡萄糖和琼脂，加热溶解后分装到三角瓶中，121℃高压灭菌20min。待培养基冷却至50℃左右时，在无菌操作台上将金龟子绿僵菌接种到培养基上，置于25℃恒温培养箱中培养7天。培养过程中，观察真菌的生长情况，待菌落长满培养基表面且孢子成熟后，用于后续实验。在进行飞蝗感染实验前，将培养好的金龟子绿僵菌用无菌水冲洗，制成浓度为1×10^8个/mL的孢子悬浮液，使用血球计数板进行孢子浓度的精确计数。为确保孢子悬浮液的均匀性，在使用前需充分振荡。2.2实验方法2.2.1真菌感染实验设计采用体表涂抹法对群居型和散居型飞蝗进行真菌感染。将饲养至5龄的飞蝗饥饿处理12h，使其处于空腹状态，以便更好地吸收真菌孢子。在无菌操作台上，用移液器吸取5μL浓度为1×10^8个/mL的金龟子绿僵菌孢子悬浮液，均匀涂抹于飞蝗的胸部背板和腹部体表。对照组则涂抹等量的无菌水。每组设置3个重复，每个重复包含20只飞蝗。感染后的飞蝗继续饲养在相同的环境条件下，分别在感染后0h、12h、24h、36h、48h、60h和72h进行观察，记录飞蝗的行为变化、感染症状等信息。例如，观察飞蝗是否出现行动迟缓、食欲不振、体色改变等症状，以及这些症状出现的时间和发展趋势。2.2.2样本收集与处理在感染后的不同时间点（0h、12h、24h、36h、48h、60h和72h），每组随机选取5只飞蝗进行样本收集。收集的样本部位包括脂肪体、血淋巴和中肠，这些部位在昆虫免疫过程中发挥着重要作用。收集脂肪体时，将飞蝗用75%酒精消毒后，在冰上解剖，迅速取出脂肪体，放入预冷的RNase-free离心管中，加入1mLTRIzol试剂，立即用匀浆器匀浆，使脂肪体充分裂解。血淋巴的收集则是用无菌注射器从飞蝗的腹部节间膜处抽取，将抽取的血淋巴注入含有抗凝剂（0.1%肝素钠）的离心管中，轻轻混匀。中肠的收集方法与脂肪体类似，解剖飞蝗后取出中肠，用预冷的PBS缓冲液冲洗干净，去除表面的杂质，然后放入含有TRIzol试剂的离心管中匀浆。匀浆后的样本在4℃下以12000rpm的转速离心10min，取上清液转移至新的离心管中，用于后续的RNA提取。2.2.3RNA提取与测序使用TRIzol试剂提取飞蝗样本中的总RNA。具体步骤如下：将上述匀浆后的上清液在室温下静置5min，使核酸蛋白复合物充分解离。加入0.2mL氯仿，剧烈振荡15s，室温下静置3min。然后在4℃下以12000rpm的转速离心15min，此时溶液分为三层，上层为无色透明的水相，含有RNA；中间层为白色的蛋白层；下层为红色的有机相。小心吸取上层水相转移至新的离心管中，加入0.5mL异丙醇，轻轻混匀，室温下静置10min，使RNA沉淀。再次在4℃下以12000rpm的转速离心10min，弃上清液，此时可见离心管底部有白色的RNA沉淀。用1mL75%乙醇洗涤RNA沉淀2次，每次在4℃下以7500rpm的转速离心5min，弃上清液。最后将RNA沉淀在室温下晾干5-10min，加入适量的RNase-free水溶解RNA。使用NanoDrop2000超微量分光光度计检测RNA的浓度和纯度，要求RNA的OD260/OD280比值在1.8-2.0之间，OD260/OD230比值大于2.0。同时，用1%的琼脂糖凝胶电泳检测RNA的完整性，观察28S和18SrRNA条带的亮度和清晰度，28SrRNA条带的亮度应约为18SrRNA条带的2倍。将提取的高质量RNA送往专业的测序公司进行高通量测序。采用IlluminaHiSeq平台进行测序，构建链特异性文库，测序读长为150bp，双端测序。测序过程严格按照标准操作流程进行，以确保数据的准确性和可靠性。2.2.4可变剪切分析方法利用生物信息学工具对测序数据进行可变剪切分析。首先，使用TopHat软件将测序得到的reads比对到飞蝗的参考基因组上，参数设置为默认值。然后，利用Cufflinks软件对转录本进行组装和定量，识别出可变剪切事件。在组装转录本时，设置最小转录本长度为200bp，最小外显子长度为50bp。根据Cufflinks软件的分析结果，将可变剪切事件分为5种类型：外显子跳跃（exonskipping）、内含子保留（intronretention）、可变5'端剪接位点（alternative5'splicesite）、可变3'端剪接位点（alternative3'splicesite）和互斥外显子（mutuallyexclusiveexons）。通过计算不同类型可变剪切事件在不同样本中的发生频率，分析群居型和散居型飞蝗在真菌感染前后可变剪切模式的差异。为了筛选出差异可变剪切基因，采用Cuffdiff软件进行差异分析，设置FDR（falsediscoveryrate）<0.05和|log2(foldchange)|>1作为筛选标准。对筛选出的差异可变剪切基因进行功能注释和富集分析，使用DAVID在线工具对基因进行GO（GeneOntology）功能富集分析和KEGG（KyotoEncyclopediaofGenesandGenomes）通路富集分析，了解这些基因参与的生物学过程和信号通路。2.2.5验证实验设计为了验证高通量测序结果的准确性，采用qRT-PCR和RT-PCR方法对部分差异可变剪切基因进行验证。对于qRT-PCR实验，根据差异可变剪切基因的序列设计特异性引物，引物设计原则为：引物长度为18-25bp，GC含量在40%-60%之间，退火温度在58-62℃之间。以飞蝗的β-actin基因作为内参基因，对每个样本进行3次重复检测。使用SYBRGreen荧光染料法进行qRT-PCR反应，反应体系为20μL，包括10μLSYBRGreenMasterMix、0.5μL上游引物（10μM）、0.5μL下游引物（10μM）、2μLcDNA模板和7μLRNase-free水。反应程序为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；最后进行熔解曲线分析，以确保扩增的特异性。根据qRT-PCR结果，计算差异可变剪切基因在不同样本中的相对表达量，采用2^-ΔΔCt法进行计算。对于RT-PCR实验，同样根据差异可变剪切基因的序列设计引物，扩增包含可变剪切位点的片段。反应体系为25μL，包括12.5μL2×TaqPCRMasterMix、0.5μL上游引物（10μM）、0.5μL下游引物（10μM）、2μLcDNA模板和9.5μLddH2O。反应程序为：95℃预变性5min；95℃变性30s，58℃退火30s，72℃延伸1min，共35个循环；最后72℃延伸10min。PCR产物用1.5%的琼脂糖凝胶电泳进行检测，观察不同样本中扩增条带的大小和亮度，判断可变剪切事件的发生情况。三、群居型和散居型飞蝗的基础特征3.1形态特征差异群居型与散居型飞蝗在形态特征上存在显著不同。在体型方面，一般而言，群居型飞蝗的体型相对较大，这使其在长途迁飞过程中更具优势，能够更好地适应远距离飞行所消耗的能量需求。例如，东亚飞蝗的群居型成虫体长普遍比散居型更长，雄成虫体长可达33-48毫米，雌成虫体长39-52毫米；而散居型成虫的体型相对较小，在长度上略逊一筹。体色也是两者的重要区别之一。群居型飞蝗的体色通常为赤褐、黑褐色，其黑褐色斑纹较多且固定，颜色较深。这种较深的体色在其大规模聚集飞行时，能够形成醒目的视觉效果，可能有助于群体内个体之间的识别和交流。相比之下，散居型飞蝗的体色随周围环境而异，多呈现草绿色或淡褐色，黑褐色斑纹少而淡，头、前胸背板及后足腿节常为绿色。这种与环境相似的保护色有助于散居型飞蝗在自然环境中隐藏自己，避免被捕食者发现。翅长的差异在两者之间也较为明显。群居型飞蝗的前翅较长，通常超过后足股节2倍以上，发达的翅膀赋予了它们强大的飞行能力，使其能够进行远距离的迁飞，寻找更适宜的生存和繁殖环境。而散居型飞蝗的前翅相对较短，不到后足股节的2倍，这也导致其飞行能力相对较弱，活动范围较为有限，多在相对较小的区域内活动。此外，头胸高度比在群居型和散居型飞蝗之间也有所不同。群居型飞蝗头顶稍高于前胸背板，而散居型飞蝗头顶低于前胸背板。后足腿节的特征同样存在差异，群居型飞蝗后足腿节稍短于或约等于前翅长度的一半，且通常等于或稍大于2（雄2.0-2.17，雌1.78-2.22）；散居型飞蝗后足腿节通常稍长于前翅长度的一半，且通常不大于2（雌雄均为1.8-1.96)。这些细微的形态差异反映了两种生活型飞蝗在长期进化过程中对不同生存策略的适应。3.2生态习性差异在栖息环境方面，群居型飞蝗偏爱地势低洼、易涝易旱或水位不稳定的海滩、湖滩、河滩、荒地或耕作粗放的农田。这些区域通常生长着大量芦苇、盐蒿、稗草、荻草、莎草等飞蝗嗜食植物，为其提供了丰富的食物来源。例如，在中国的黄淮海平原，滨湖蝗区、沿海蝗区、内涝蝗区和河泛蝗区等多种类型的蝗区，都是群居型飞蝗的适宜栖息地。在干旱年份，这些区域的宜蝗荒滩、荒地面积增大，有利于群居型飞蝗的繁衍，容易酿成蝗灾。散居型飞蝗则相对较为分散，常栖息于农作物种植较为分散、飞蝗密度相对较低的区域，如农田周边、草地边缘等。这些地方虽然食物资源相对不如群居型飞蝗栖息地丰富，但也能满足散居型飞蝗的生存需求。由于其分散的特性，散居型飞蝗对单一区域的危害相对较小。在取食偏好上，两者都偏好禾本科和莎草科植物，像麦类、玉米、高粱、谷子、水稻、糜子等禾本科作物，以及芦苇、狗尾草、稗草等莎草科杂草，都是它们喜爱的食物。在食物充足的情况下，群居型飞蝗由于群体数量庞大，取食更为集中，往往会对农作物造成大面积的破坏，所到之处农作物被迅速啃食，常常导致颗粒无收。而散居型飞蝗由于个体分散，取食相对分散，对农作物的危害程度相对较轻。但当食物资源短缺时，散居型飞蝗也可能会聚集起来，寻找更多的食物，此时其取食行为可能会向群居型飞蝗靠拢。繁殖行为方面，东亚飞蝗无滞育现象，发生代数由北向南递增，各地都以卵囊在土壤中越冬。大致北京以北地区每年发生1代，黄淮流域一年2代，长江中下游2-3代，华南3代，海南4代。在2代区，越冬卵于4月底至5月上中旬孵化，蝗蝻在6月中旬至7月上旬羽化。这一代蝗虫称为夏蝗。夏蝗寿命55-60天，7月上中旬进入产卵盛期，卵在7月中旬至8月中旬孵化，蝗蝻25-30天后羽化为成虫，称为秋蝗。秋蝗在9月上中旬进入产卵盛期，大部分卵直接越冬，少数在高温干旱年份，在9月中旬前后又孵出第三代蝗蝻，羽化后不能产卵，多在冬季被冻死。群居型飞蝗和散居型飞蝗在繁殖过程中也存在差异。成虫羽化后，需经1-2周时间补充营养和飞翔，才能交配产卵。产卵多选择植被稀疏，土面坚实，土壤湿度和含盐量适宜的向阳地带。夏蝗旱季多在水畔荒滩和低洼农田产卵，秋蝗在汛期多选择未被水淹的河畔海滩或附近高燥农田产卵。雌虫用产卵瓣钻土成孔，深入土面下4-6厘米产卵。每只雌蝗一般产4-5个卵块，每个卵块有卵约50-75粒。群居型飞蝗体内含脂肪量多，水分少，活动力强，但卵巢管数少，产卵量低；而散居型则相反，卵巢管数较多，产卵量相对较高。这可能是由于群居型飞蝗需要将更多的能量用于迁飞和群体活动，而散居型飞蝗则将更多能量投入到繁殖后代上。3.3免疫相关基础差异飞蝗作为昆虫，主要依赖固有免疫来抵御病原体的入侵。固有免疫是飞蝗在长期进化过程中形成的一种天然防御机制，能够对病原体做出快速响应。它包括物理防御和免疫细胞介导的防御以及体液免疫等多个方面。在物理防御层面，飞蝗的表皮和肠道上皮等屏障结构起着重要作用。表皮作为飞蝗与外界环境接触的第一道防线，具有多层结构，包括表皮层、真皮细胞层等。表皮层中的几丁质和蛋白质等成分构成了坚韧的物理屏障，能够阻止真菌等病原体的入侵。同时，表皮上还分布有蜡质层，进一步增强了其防水和防病原体的能力。肠道上皮则是飞蝗抵御经口感染病原体的重要防线。肠道上皮细胞紧密排列，形成了一道物理屏障，能够防止病原体穿透肠道进入体内。此外，肠道上皮细胞还能分泌黏液，黏液中含有多种抗菌物质，如抗菌肽、溶菌酶等，能够抑制病原体的生长和繁殖。在免疫细胞介导的防御方面，飞蝗的血细胞在识别和清除病原体过程中发挥着关键作用。飞蝗的血细胞主要包括原血细胞、浆血细胞、颗粒血细胞等。原血细胞是血细胞的前体细胞，具有分化为其他类型血细胞的能力。浆血细胞能够通过吞噬作用摄取病原体，将其包裹在吞噬体中，然后与溶酶体融合，利用溶酶体中的各种酶类将病原体降解。颗粒血细胞则含有丰富的颗粒物质，这些颗粒中包含多种抗菌物质和免疫调节因子。当颗粒血细胞受到病原体刺激时，会释放这些颗粒物质，对病原体进行杀伤和免疫调节。在体液免疫方面，飞蝗在受到真菌感染后，会激活一系列免疫信号通路，如Toll信号通路、Imd信号通路等。Toll信号通路在飞蝗抵御真菌感染中起着至关重要的作用。当飞蝗的免疫细胞识别到真菌细胞壁上的特定分子模式，如β-1,3-葡聚糖等时，会激活Toll受体，进而引发一系列的信号转导过程。Toll受体激活后，会招募下游的信号分子，如MyD88、Tube等，形成信号复合物。这个信号复合物会进一步激活蛋白激酶Pelle，Pelle会磷酸化并激活转录因子Dorsal和Relish。激活后的Dorsal和Relish会进入细胞核，与相关基因的启动子区域结合，调控免疫相关基因的表达，如抗菌肽基因等。抗菌肽是飞蝗体液免疫中的重要效应分子，它们具有广谱的抗菌活性，能够对真菌、细菌等病原体产生杀伤作用。不同种类的抗菌肽具有不同的作用机制，有些抗菌肽能够破坏病原体的细胞膜，导致细胞内容物泄漏；有些抗菌肽则能够抑制病原体的蛋白质合成或核酸复制。群居型和散居型飞蝗在固有免疫的各个环节可能存在差异。在物理防御方面，两者表皮和肠道上皮的结构和组成可能有所不同。例如，群居型飞蝗由于其聚集生活的特性，可能面临更高的病原体传播风险，因此其表皮的结构可能更加致密，蜡质层更厚，以增强对病原体的防御能力。而散居型飞蝗由于生活相对分散，病原体传播风险相对较低，其表皮结构可能相对较疏松。在免疫细胞介导的防御方面，群居型和散居型飞蝗的血细胞数量、种类和功能可能存在差异。研究表明，群居型飞蝗的血细胞数量可能相对较多，且浆血细胞和颗粒血细胞的活性可能更高，这可能与其在高密度群体中需要更强的免疫防御能力有关。在体液免疫方面，群居型和散居型飞蝗的免疫信号通路的激活程度和抗菌肽的表达水平可能存在差异。当受到真菌感染时，群居型飞蝗可能更快地激活Toll信号通路和Imd信号通路，从而更迅速地启动免疫应答。同时，群居型飞蝗可能表达更高水平的抗菌肽，以应对群体中可能爆发的病原体感染。而散居型飞蝗由于个体之间接触较少，病原体传播相对较慢，其免疫信号通路的激活和抗菌肽的表达可能相对较弱。四、真菌感染对飞蝗的影响4.1感染后的生理变化在感染初期，飞蝗的行为会发生明显改变。飞蝗开始表现出行动迟缓的症状，其正常的爬行和跳跃能力受到显著影响，移动速度明显减慢。在饲养箱中观察发现，感染后的飞蝗常常长时间静止不动，不再像健康个体那样活跃地四处觅食和活动。部分飞蝗还会出现食欲减退的现象，对原本喜爱的小麦苗等食物表现出明显的拒食行为。这种食欲的下降可能是由于真菌感染引发的生理不适，导致飞蝗的摄食欲望降低。同时，飞蝗的反应能力也有所下降，对外部刺激的敏感度降低，例如当用小木棍触碰感染后的飞蝗时，其逃避反应明显迟缓。随着感染时间的延长，飞蝗的生长发育受到严重阻碍。在体型发育方面，与未感染的飞蝗相比，感染真菌的飞蝗生长速度明显减缓，体型较小。通过定期测量飞蝗的体长和体重发现，感染后36h的飞蝗体长增长速度比对照组降低了约30%，体重增长也明显滞后。在蜕皮过程中，感染真菌的飞蝗常常出现蜕皮困难的情况。蜕皮是昆虫生长发育的重要阶段，正常情况下，飞蝗能够顺利完成蜕皮，实现体型的增长和形态的转变。然而，感染真菌后，部分飞蝗在蜕皮时会出现旧表皮难以脱落的现象，导致蜕皮过程延长，甚至因蜕皮失败而死亡。据统计，感染真菌的飞蝗蜕皮成功率比对照组降低了约40%。此外，真菌感染还可能影响飞蝗的生殖器官发育，导致生殖能力下降。在对感染后的雌性飞蝗解剖观察中发现，其卵巢发育不完整，卵巢管数量减少，卵细胞的发育也受到抑制。在感染后期，飞蝗的体色会发生显著变化。原本绿色或褐色的体色逐渐变为黑色或深褐色，这种颜色的改变在飞蝗的体表各个部位都较为明显，尤其是胸部和腹部。飞蝗的身体表面还会出现一些白色的菌丝体，这些菌丝体从飞蝗的气孔、节间膜等部位生长出来，逐渐覆盖飞蝗的体表。随着菌丝体的生长，飞蝗的身体逐渐变得僵硬，失去了正常的柔韧性和活动能力。最终，飞蝗因真菌的持续侵染和体内营养物质的大量消耗而死亡。在感染后72h，大部分感染真菌的飞蝗已经死亡，其尸体干瘪，附着着大量的菌丝体。4.2感染后的免疫反应飞蝗免疫系统对真菌感染的应答是一个复杂而有序的过程，涉及多种免疫机制的协同作用。当飞蝗受到金龟子绿僵菌感染时，其免疫系统迅速启动，首先通过模式识别受体（PRRs）识别真菌细胞壁上的特定分子模式，如β-1,3-葡聚糖等。这些PRRs广泛存在于飞蝗的免疫细胞表面以及一些组织细胞表面，能够特异性地结合真菌的特征分子，从而触发免疫信号通路的激活。在免疫细胞变化方面，飞蝗的血细胞在真菌感染后发挥着关键作用。研究发现，感染真菌后，飞蝗体内的血细胞数量会发生明显变化。在感染初期，血细胞数量迅速增加，这是机体为了增强免疫防御能力而做出的反应。通过显微镜观察发现，在感染后12h，飞蝗血淋巴中的血细胞数量相比未感染组增加了约50%。这些增加的血细胞主要包括浆血细胞和颗粒血细胞。浆血细胞具有较强的吞噬能力，能够通过伸出伪足将真菌孢子包裹起来，形成吞噬体，然后与溶酶体融合，利用溶酶体中的各种酶类将真菌孢子降解。颗粒血细胞则能够释放多种抗菌物质和免疫调节因子，如抗菌肽、活性氧等。这些物质能够直接杀伤真菌，或者调节其他免疫细胞的活性，增强飞蝗的免疫防御能力。抗菌肽的产生是飞蝗免疫应答的重要环节。当飞蝗受到真菌感染时，体内的免疫信号通路被激活，其中Toll信号通路和Imd信号通路在抗菌肽的诱导表达中发挥着关键作用。在Toll信号通路中，模式识别受体识别真菌细胞壁成分后，激活Toll受体，进而招募下游的信号分子MyD88、Tube等，形成信号复合物。该复合物激活蛋白激酶Pelle，Pelle磷酸化并激活转录因子Dorsal和Relish。激活后的Dorsal和Relish进入细胞核，与抗菌肽基因的启动子区域结合，促进抗菌肽基因的转录和表达。在Imd信号通路中，免疫受体识别病原体相关分子模式后，激活Imd蛋白，Imd蛋白招募下游的信号分子FADD、Dredd等，形成信号复合物。该复合物激活转录因子Relish，Relish进入细胞核，调控抗菌肽基因的表达。通过实时定量PCR技术检测发现，在感染金龟子绿僵菌后，飞蝗体内多种抗菌肽基因的表达水平显著上调。例如，防御素（Defensin）基因在感染后24h的表达量相比未感染组增加了约10倍，天蚕素（Cecropin）基因的表达量在感染后36h增加了约8倍。这些抗菌肽具有广谱的抗菌活性，能够对真菌、细菌等病原体产生杀伤作用。它们通过破坏病原体的细胞膜、抑制病原体的蛋白质合成或核酸复制等方式，发挥免疫防御功能。除了血细胞和抗菌肽，飞蝗的脂肪体在免疫应答中也扮演着重要角色。脂肪体是飞蝗体内重要的代谢和免疫器官，在真菌感染后，脂肪体会合成和分泌多种免疫相关蛋白。研究表明，脂肪体中一些免疫相关基因的表达在感染后发生显著变化。例如，免疫球蛋白超家族蛋白（Immunoglobulinsuperfamilyproteins）基因在感染后表达上调，这些蛋白可能参与了免疫细胞的识别和黏附过程，增强了免疫细胞对真菌的清除能力。此外，脂肪体还能够合成和分泌一些蛋白酶抑制剂，这些抑制剂可以抑制真菌分泌的蛋白酶活性，从而阻止真菌对飞蝗组织的破坏。4.3群居型与散居型飞蝗感染差异对比在感染率方面，研究结果显示，散居型飞蝗的感染率相对较高。在相同的感染条件下，即使用浓度为1×10^8个/mL的金龟子绿僵菌孢子悬浮液进行体表涂抹感染，感染后48h，散居型飞蝗的感染率达到了约70%，而群居型飞蝗的感染率仅为约40%。这表明群居型飞蝗对金龟子绿僵菌的感染具有更强的抵抗力，可能是由于群居型飞蝗在长期的进化过程中，形成了更为有效的防御机制来应对真菌的侵染。死亡率的差异也较为显著。随着感染时间的延长，散居型飞蝗的死亡率迅速上升。在感染后72h，散居型飞蝗的死亡率高达85%，而群居型飞蝗的死亡率为55%。这进一步说明了群居型飞蝗对真菌的耐受性更强，能够在更长的时间内抵御真菌的侵害。从免疫反应强度来看，通过对血细胞数量变化和抗菌肽基因表达水平的检测，发现群居型飞蝗在感染后的免疫反应更为强烈。在血细胞数量变化方面，感染后12h，群居型飞蝗血淋巴中的血细胞数量相比未感染组增加了约70%，而散居型飞蝗血细胞数量增加了约50%。这表明群居型飞蝗能够更快地动员血细胞参与免疫防御，增强对真菌的清除能力。在抗菌肽基因表达水平上，以防御素（Defensin）基因和天蚕素（Cecropin）基因的表达为例，感染后24h，群居型飞蝗体内防御素基因的表达量相比未感染组增加了约15倍，而散居型飞蝗增加了约10倍；感染后36h，群居型飞蝗天蚕素基因的表达量增加了约12倍，散居型飞蝗增加了约8倍。这说明群居型飞蝗在受到真菌感染时，能够更有效地激活免疫信号通路，诱导抗菌肽基因的表达，从而增强免疫防御能力。此外，在免疫相关基因的表达模式上，群居型和散居型飞蝗也存在差异。通过对Toll信号通路和Imd信号通路中相关基因的表达分析发现，群居型飞蝗在感染后，Toll信号通路中的关键基因如Toll、MyD88等的表达上调更为迅速和显著。在感染后8h，群居型飞蝗脂肪体中Toll基因的表达量就开始显著上调，而散居型飞蝗在感染后12h才出现明显的上调。这表明群居型飞蝗能够更快地启动Toll信号通路，从而更早地开始免疫应答。在Imd信号通路中，群居型飞蝗的Relish基因在感染后的表达水平也明显高于散居型飞蝗，这进一步说明群居型飞蝗在免疫反应中对Imd信号通路的激活更为强烈。五、飞蝗可变剪切的数据分析5.1可变剪切事件鉴定通过对测序数据进行深入分析，本研究共鉴定出了多种类型的可变剪切事件。其中，外显子跳跃（exonskipping）事件是最为常见的类型之一，共鉴定出5628个，占总可变剪切事件的35.4%。外显子跳跃是指在mRNA前体剪接过程中，某个或某些外显子被跳过，不参与成熟mRNA的形成。这种可变剪切方式可以产生不同长度和结构的mRNA转录本，进而翻译出具有不同功能的蛋白质异构体。例如，在飞蝗的免疫相关基因中，外显子跳跃可能导致蛋白质结构域的缺失或增加，从而影响蛋白质与病原体的相互作用能力。内含子保留（intronretention）事件鉴定出2764个，占比17.4%。内含子保留是指在剪接过程中，某些内含子没有被完全切除，而是保留在成熟mRNA中。保留的内含子可能会影响mRNA的稳定性、翻译效率以及蛋白质的功能。在飞蝗中，内含子保留可能与免疫信号通路的调控有关，通过调节相关基因的表达水平来应对真菌感染。可变5'端剪接位点（alternative5'splicesite）事件有2456个，占15.4%。可变5'端剪接位点是指mRNA前体在剪接时，5'端的剪接位点发生变化，导致不同的外显子被连接在一起。这种可变剪切方式可以改变蛋白质的N端序列，进而影响蛋白质的功能。在飞蝗免疫过程中，可变5'端剪接位点的变化可能会影响免疫相关蛋白与其他分子的相互作用，从而调节免疫应答的强度。可变3'端剪接位点（alternative3'splicesite）事件数量为2212个，占13.9%。可变3'端剪接位点与可变5'端剪接位点类似，只是剪接位点的变化发生在3'端。这种可变剪切方式会改变蛋白质的C端序列，对蛋白质的功能产生影响。在飞蝗真菌感染前后，可变3'端剪接位点的变化可能与免疫相关基因的表达调控以及免疫效应分子的活性有关。互斥外显子（mutuallyexclusiveexons）事件鉴定出1920个，占12.0%。互斥外显子是指在mRNA前体剪接过程中，多个外显子之间存在互斥关系，只能选择其中一个外显子参与成熟mRNA的形成。这种可变剪切方式可以产生多种不同的mRNA转录本，增加了蛋白质组的复杂性。在飞蝗中，互斥外显子可能参与了免疫细胞的分化和免疫功能的调节，通过产生不同的蛋白质异构体来适应不同的免疫需求。其他类型的可变剪切事件相对较少，共鉴定出620个，占3.9%。这些可变剪切事件虽然数量较少，但也可能在飞蝗的生理过程中发挥着重要作用，需要进一步深入研究。通过对不同类型可变剪切事件的鉴定和分析，我们对飞蝗在真菌感染前后的基因表达调控机制有了更深入的了解。这些可变剪切事件的发生可能与飞蝗的免疫应答、生长发育等过程密切相关，为后续的研究提供了重要的线索。5.2群居型飞蝗感染前后可变剪切变化在基因表达水平上，对群居型飞蝗感染金龟子绿僵菌前后的转录组数据分析发现，共有1256个基因发生了可变剪切变化。其中，892个基因的可变剪切在感染后显著上调，这意味着这些基因产生了更多不同的转录本形式。例如，基因LOC105372345在感染前主要产生一种转录本，而在感染后通过可变剪切产生了三种不同的转录本。这些新产生的转录本可能编码具有不同功能的蛋白质异构体，从而增强群居型飞蝗对真菌感染的免疫应答能力。通过对这些上调可变剪切基因的功能分析，发现它们主要参与了免疫防御、细胞代谢和信号转导等生物学过程。在免疫防御方面，一些基因编码的蛋白质可能具有更强的抗菌活性，能够更有效地抑制真菌的生长和繁殖；在细胞代谢方面，相关基因的可变剪切可能调节细胞内的代谢途径，为免疫应答提供更多的能量和物质基础。另一方面，364个基因的可变剪切在感染后显著下调。以基因LOC105371892为例，感染前该基因存在多种可变剪切形式，产生多个转录本，但在感染后，其可变剪切事件减少，转录本种类也相应减少。这些下调可变剪切的基因可能在正常生理状态下发挥重要作用，但在真菌感染时，其功能可能受到抑制，以避免对免疫应答产生干扰。对这些下调基因的功能研究发现，部分基因与飞蝗的生长发育相关，在感染期间，飞蝗可能通过减少这些基因的可变剪切，将更多的能量和资源投入到免疫防御中。从功能角度分析，对发生可变剪切变化的基因进行GO功能富集分析，结果显示，在生物过程（biologicalprocess）类别中，“免疫应答”（immuneresponse）、“抗菌体液应答”（antimicrobialhumoralresponse）和“活性氧代谢过程”（reactiveoxygenspeciesmetabolicprocess）等功能条目显著富集。在“免疫应答”过程中，多个免疫相关基因发生可变剪切，如Toll样受体基因（Toll-likereceptorgenes），其不同的转录本可能具有不同的配体结合能力，从而影响免疫信号的传递和免疫应答的强度。在“抗菌体液应答”中，抗菌肽基因的可变剪切可能产生具有不同抗菌活性的抗菌肽异构体，增强飞蝗对真菌的抵抗能力。“活性氧代谢过程”的富集表明，可变剪切可能调节飞蝗体内活性氧的产生和清除，活性氧在免疫防御中具有重要作用，适量的活性氧可以杀伤病原体，但过多的活性氧会对细胞造成损伤，可变剪切可能通过调节相关基因的表达，维持活性氧的平衡。在分子功能（molecularfunction）类别中，“抗菌肽活性”（antimicrobialpeptideactivity）、“氧化还原酶活性”（oxidoreductaseactivity）和“蛋白结合”（proteinbinding）等功能条目富集明显。具有“抗菌肽活性”的基因发生可变剪切，可能改变抗菌肽的结构和功能，使其对真菌的杀伤作用更加有效。“氧化还原酶活性”相关基因的可变剪切，可能影响飞蝗体内的氧化还原平衡，进而影响免疫细胞的活性和免疫应答的进程。“蛋白结合”功能的富集说明，可变剪切可能产生的蛋白质异构体具有不同的蛋白质相互作用能力，通过与其他免疫相关蛋白结合，调节免疫信号通路的激活和免疫应答的调控。在细胞组成（cellularcomponent）类别中，“细胞外区域”（extracellularregion）、“细胞膜”（cellmembrane）和“细胞质”（cytoplasm）等功能条目富集。在“细胞外区域”，可变剪切产生的蛋白质可能分泌到细胞外，参与免疫防御，如一些抗菌肽和免疫调节因子。“细胞膜”相关基因的可变剪切可能改变细胞膜的结构和功能，影响免疫细胞对病原体的识别和吞噬作用。“细胞质”中相关基因的可变剪切可能调节细胞内的免疫信号转导和免疫相关物质的合成。通过KEGG通路富集分析发现，“Toll和Imd信号通路”（TollandImdsignalingpathway）、“吞噬体”（phagosome）和“氧化磷酸化”（oxidativephosphorylation）等通路显著富集。在“Toll和Imd信号通路”中，多个关键基因发生可变剪切，如Toll受体基因、MyD88基因和Relish基因等。这些基因的可变剪切可能导致信号通路的激活程度发生改变，从而调节免疫相关基因的表达。例如，Toll受体基因的可变剪切可能产生具有不同活性的受体异构体，影响其对真菌细胞壁成分的识别和信号传递能力。在“吞噬体”通路中，可变剪切可能影响免疫细胞对真菌的吞噬和降解过程，相关基因的可变剪切可能改变吞噬体的形成和成熟过程，增强免疫细胞对真菌的清除能力。“氧化磷酸化”通路的富集表明，可变剪切可能调节飞蝗细胞的能量代谢，为免疫应答提供充足的能量。5.3散居型飞蝗感染前后可变剪切变化在基因表达层面，经深入分析转录组数据可知，散居型飞蝗在感染金龟子绿僵菌后，有1028个基因出现可变剪切变化。其中，654个基因的可变剪切水平在感染后显著上调，这表明这些基因在转录后加工过程中产生了更多种类的mRNA转录本。例如，基因LOC105370987在感染前仅产生一种主要转录本，感染后通过可变剪切产生了另外两种具有不同结构的转录本。这些新增转录本可能编码不同功能的蛋白质异构体，从而对散居型飞蝗的免疫应答产生影响。通过对这些上调可变剪切基因的功能分析，发现它们广泛参与细胞代谢、信号传导以及免疫防御等生物学过程。在细胞代谢方面，部分基因可能通过改变代谢途径，为免疫细胞提供更多的能量和物质支持；在信号传导过程中，相关基因的可变剪切可能影响免疫信号的传递效率，进而调控免疫应答的强度。另一方面，374个基因的可变剪切在感染后显著下调。以基因LOC105371235为例，感染前该基因存在多种可变剪切形式，产生多个转录本，但在感染后，其可变剪切事件减少，转录本种类也相应减少。这些下调可变剪切的基因可能在正常生理状态下发挥重要作用，但在真菌感染时，其功能可能受到抑制，以避免对免疫应答产生干扰。对这些下调基因的功能研究发现，部分基因与飞蝗的生长发育相关，在感染期间，散居型飞蝗可能通过减少这些基因的可变剪切，将更多的能量和资源投入到免疫防御中。从功能角度分析，对发生可变剪切变化的基因进行GO功能富集分析，结果显示，在生物过程（biologicalprocess）类别中，“免疫应答”（immuneresponse）、“抗菌体液应答”（antimicrobialhumoralresponse）和“活性氧代谢过程”（reactiveoxygenspeciesmetabolicprocess）等功能条目显著富集。在“免疫应答”过程中，多个免疫相关基因发生可变剪切，如免疫球蛋白超家族基因（Immunoglobulinsuperfamilygenes），其不同的转录本可能具有不同的抗原结合能力，从而影响免疫细胞对病原体的识别和清除。在“抗菌体液应答”中，抗菌肽基因的可变剪切可能产生具有不同抗菌活性的抗菌肽异构体，增强散居型飞蝗对真菌的抵抗能力。“活性氧代谢过程”的富集表明，可变剪切可能调节散居型飞蝗体内活性氧的产生和清除，活性氧在免疫防御中具有重要作用，适量的活性氧可以杀伤病原体，但过多的活性氧会对细胞造成损伤，可变剪切可能通过调节相关基因的表达，维持活性氧的平衡。在分子功能（molecularfunction）类别中，“抗菌肽活性”（antimicrobialpeptideactivity）、“氧化还原酶活性”（oxidoreductaseactivity）和“蛋白结合”（proteinbinding）等功能条目富集明显。具有“抗菌肽活性”的基因发生可变剪切，可能改变抗菌肽的结构和功能，使其对真菌的杀伤作用更加有效。“氧化还原酶活性”相关基因的可变剪切，可能影响散居型飞蝗体内的氧化还原平衡，进而影响免疫细胞的活性和免疫应答的进程。“蛋白结合”功能的富集说明，可变剪切可能产生的蛋白质异构体具有不同的蛋白质相互作用能力，通过与其他免疫相关蛋白结合，调节免疫信号通路的激活和免疫应答的调控。在细胞组成（cellularcomponent）类别中，“细胞外区域”（extracellularregion）、“细胞膜”（cellmembrane）和“细胞质”（cytoplasm）等功能条目富集。在“细胞外区域”，可变剪切产生的蛋白质可能分泌到细胞外，参与免疫防御，如一些抗菌肽和免疫调节因子。“细胞膜”相关基因的可变剪切可能改变细胞膜的结构和功能，影响免疫细胞对病原体的识别和吞噬作用。“细胞质”中相关基因的可变剪切可能调节细胞内的免疫信号转导和免疫相关物质的合成。通过KEGG通路富集分析发现，“Toll和Imd信号通路”（TollandImdsignalingpathway）、“吞噬体”（phagosome）和“氧化磷酸化”（oxidativephosphorylation）等通路显著富集。在“Toll和Imd信号通路”中，多个关键基因发生可变剪切，如Toll受体基因、MyD88基因和Relish基因等。这些基因的可变剪切可能导致信号通路的激活程度发生改变，从而调节免疫相关基因的表达。例如，Toll受体基因的可变剪切可能产生具有不同活性的受体异构体，影响其对真菌细胞壁成分的识别和信号传递能力。在“吞噬体”通路中，可变剪切可能影响免疫细胞对真菌的吞噬和降解过程，相关基因的可变剪切可能改变吞噬体的形成和成熟过程，增强免疫细胞对真菌的清除能力。“氧化磷酸化”通路的富集表明，可变剪切可能调节散居型飞蝗细胞的能量代谢，为免疫应答提供充足的能量。5.4两型飞蝗可变剪切差异对比在可变剪切事件类型和数量上，群居型飞蝗在感染后共鉴定出5种主要类型的可变剪切事件，包括外显子跳跃、内含子保留、可变5'端剪接位点、可变3'端剪接位点和互斥外显子，总计15568个事件。其中，外显子跳跃事件数量最多，为5628个，占比35.4%。散居型飞蝗在感染后同样鉴定出这5种类型的可变剪切事件，总数为13288个。外显子跳跃事件数量为4892个，占比36.8%。从整体数量上看，群居型飞蝗的可变剪切事件总数略多于散居型飞蝗。在各类可变剪切事件的占比上，两者较为相似，但外显子跳跃事件在散居型飞蝗中的占比相对更高，而内含子保留事件在群居型飞蝗中的占比相对较高。例如，内含子保留事件在群居型飞蝗中占比17.4%，在散居型飞蝗中占比16.2%。在相关基因功能方面，对群居型和散居型飞蝗感染后发生可变剪切的基因进行功能富集分析，结果显示，两者在免疫相关功能上存在一些差异。在GO功能富集分析中，在生物过程类别下，群居型飞蝗中与“免疫细胞活化”（immunecellactivation）相关的基因显著富集，这些基因的可变剪切可能通过调节免疫细胞的活性，增强群居型飞蝗的免疫防御能力。例如，一些编码免疫细胞表面受体的基因发生可变剪切，可能改变受体的结构和功能，影响免疫细胞对病原体的识别和应答。而散居型飞蝗中与“免疫信号转导负调控”（negativeregulationofimmunesignaltransduction）相关的基因富集明显，这可能表明散居型飞蝗在感染后通过调节免疫信号转导的强度，避免过度的免疫反应对自身造成损伤。例如，某些信号通路中的抑制因子基因发生可变剪切，可能产生具有更强抑制活性的蛋白质异构体，从而调控免疫信号的传递。在分子功能类别下，群居型飞蝗中具有“趋化因子活性”（chemokineactivity）的基因发生可变剪切，趋化因子能够吸引免疫细胞向感染部位聚集，增强免疫细胞对病原体的清除能力。而散居型飞蝗中与“蛋白酶抑制剂活性”（proteaseinhibitoractivity）相关的基因可变剪切显著，蛋白酶抑制剂可以抑制真菌分泌的蛋白酶活性，阻止真菌对飞蝗组织的破坏。在KEGG通路富集分析中，群居型飞蝗中“细胞因子-细胞因子受体相互作用”（Cytokine-cytokinereceptorinteraction）通路显著富集，该通路中的基因可变剪切可能影响免疫细胞之间的通讯和协作，增强群居型飞蝗的免疫应答。散居型飞蝗中“吞噬体”（phagosome）通路的基因可变剪切更为突出，这可能与散居型飞蝗免疫细胞对真菌的吞噬和降解过程密切相关。例如，参与吞噬体形成和成熟的相关基因发生可变剪切，可能改变吞噬体的结构和功能，影响免疫细胞对真菌的清除效率。六、关键可变剪切基因的功能分析6.1筛选关键可变剪切基因通过对群居型和散居型飞蝗感染前后可变剪切数据分析，依据差异显著性及功能重要性，筛选出多个关键可变剪切基因。在群居型飞蝗中，基因LOC105372345在感染后可变剪切显著上调，其编码的蛋白可能参与免疫细胞的活化过程，对增强免疫防御具有重要作用。基因LOC105371892在感染后可变剪切显著下调，该基因在正常生理状态下可能参与飞蝗的生长发育过程，但在真菌感染时，其功能可能受到抑制，以避免对免疫应答产生干扰。在散居型飞蝗中，基因LOC105370987感染后可变剪切上调明显，可能通过调节细胞代谢途径，为免疫细胞提供更多的能量和物质支持，从而增强免疫应答。基因LOC105371235可变剪切下调显著，可能与散居型飞蝗在感染后对免疫信号转导的调控有关，通过减少该基因的可变剪切，避免过度的免疫反应对自身造成损伤。这些筛选出的关键可变剪切基因，在两型飞蝗的免疫应答过程中具有重要功能，为进一步研究飞蝗免疫机制提供了关键靶点。6.2基因功能预测利用生物信息学工具对筛选出的关键可变剪切基因进行功能预测，发现这些基因参与了多个重要的信号通路和生物学过程。在信号通路方面，部分基因参与了Toll和Imd信号通路。例如，基因LOC105372345编码的蛋白含有与Toll受体相似的结构域，推测其可能在Toll信号通路中发挥作用。当飞蝗受到真菌感染时，Toll受体识别真菌细胞壁上的β-1,3-葡聚糖等分子模式，激活下游的信号转导过程。LOC105372345基因的可变剪切可能影响其编码蛋白与Toll受体的相互作用，从而调节Toll信号通路的激活程度，进而影响免疫相关基因的表达和免疫应答的强度。另一关键基因LOC105370987与Imd信号通路相关。Imd信号通路在飞蝗抵御革兰氏阴性菌和某些真菌的感染中起着重要作用。该基因编码的蛋白可能参与了Imd信号通路中信号复合物的形成，通过可变剪切产生的不同转录本可能影响蛋白的结构和功能，从而改变信号复合物的活性，调控免疫相关基因的表达。例如，可变剪切可能导致蛋白的某些结构域发生变化，影响其与其他信号分子的结合能力，进而影响Imd信号通路的传递效率。在生物学过程方面，基因LOC105371892参与了细胞代谢过程。通过对其基因序列和编码蛋白的结构分析，发现该蛋白含有多个与代谢酶相关的结构域。在飞蝗正常生长发育过程中，该基因可能通过可变剪切产生不同的转录本，编码具有不同催化活性的代谢酶，参与细胞内的物质合成与分解代谢，为细胞的正常生理活动提供能量和物质基础。而在真菌感染后，其可变剪切模式发生改变，可能导致参与代谢的酶的种类和活性发生变化，从而影响细胞代谢过程，为免疫应答提供更多的能量和物质支持。基因LOC105371235则与免疫调节过程密切相关。该基因编码的蛋白可能作为免疫调节因子，通过可变剪切产生的不同转录本，调控免疫细胞的活性和免疫应答的进程。例如，某些转录本编码的蛋白可能具有抑制免疫细胞过度活化的功能，避免免疫反应对飞蝗自身组织造成损伤。在真菌感染时，可变剪切可能使该基因产生具有更强免疫调节活性的转录本，从而维持免疫平衡，增强飞蝗对真菌的抵抗能力。6.3功能验证实验为了进一步验证关键可变剪切基因在飞蝗免疫应答和真菌感染过程中的功能，本研究设计并实施了基因敲除和过表达实验。在基因敲除实验中，采用CRISPR/Cas9技术对关键可变剪切基因进行敲除。以群居型飞蝗中的基因LOC105372345为例，首先根据该基因的序列设计特异性的sgRNA。通过生物信息学分析，选择了位于基因编码区的两个位点作为sgRNA的靶向位点，以确保能够有效地敲除基因。将设计好的sgRNA与Cas9蛋白混合，形成RNP复合物。然后，利用显微注射技术将RNP复合物注射到飞蝗的受精卵中。在注射过程中，严格控制注射量和注射压力，以保证受精卵的存活率。注射后的受精卵在适宜的条件下孵化，待孵化出的飞蝗幼虫生长至3龄时，提取其基因组DNA，通过PCR扩增和测序验证基因敲除的效果。结果显示，在成功敲除LOC105372345基因的飞蝗中，该基因的编码序列发生了缺失或突变，导致基因功能丧失。对基因敲除后的飞蝗进行真菌感染实验，观察其免疫应答和感染情况。将敲除基因的飞蝗和野生型飞蝗同时用金龟子绿僵菌进行体表涂抹感染，感染后定期观察飞蝗的行为变化和感染症状。结果发现，敲除LOC105372345基因的飞蝗在感染后行动迟缓、食欲不振等症状出现的时间更早，且感染率和死亡率明显高于野生型飞蝗。在感染后48h，敲除基因的飞蝗感染率达到了80%，而野生型飞蝗感染率为40%；在感染后72h，敲除基因的飞蝗死亡率高达90%，而野生型飞蝗死亡率为55%。通过检测免疫相关指标发现，敲除基因的飞蝗在感染后血细胞数量的增加幅度明显低于野生型飞蝗，抗菌肽基因的表达水平也显著降低。例如，防御素（Defensin）基因的表达量在敲除基因的飞蝗中仅为野生型飞蝗的30%。这表明LOC105372345基因在飞蝗的免疫应答中发挥着重要作用，敲除该基因会削弱飞蝗对真菌感染的抵抗能力。对于散居型飞蝗中的基因LOC105370987，同样采用CRISPR/Cas9技术进行基因敲除。设计并合成针对该基因的sgRNA，将其与Cas9蛋白形成RNP复合物后注射到散居型飞蝗的受精卵中。经过验证，成功获得了基因敲除的飞蝗。对基因敲除的散居型飞蝗进行真菌感染实验，结果表明，敲除LOC105370987基因的飞蝗在感染后的生长发育受到更严重的阻碍，体型明显小于野生型飞蝗。在感染后36h，敲除基因的飞蝗体长增长速度比野生型飞蝗降低了约50%。同时，敲除基因的飞蝗免疫细胞对真菌的吞噬能力下降，免疫相关基因的表达受到抑制。这说明LOC105370987基因在散居型飞蝗的免疫应答和生长发育过程中具有重要功能。在过表达实验中，构建关键可变剪切基因的过表达载体。以群居型飞蝗中的基因LOC105371892为例，通过PCR扩增获得该基因的全长编码序列，将其克隆到pUAST过表达载体中。利用脂质体转染法将构建好的过表达载体导入飞蝗的细胞系中，筛选出稳定表达该基因的细胞克隆。将过表达细胞克隆注射到飞蝗的血淋巴中，使基因在飞蝗体内实现过表达。对过表达基因的飞蝗进行真菌感染实验，观察其免疫应答和感染情况。与野生型飞蝗相比，过表达LOC105371892基因的飞蝗在感染后免疫应答更为迅速和强烈。在感染后12h，过表达基因的飞蝗血淋巴中的血细胞数量相比野生型飞蝗增加了约30%，抗菌肽基因的表达水平也显著上调。例如，天蚕素（Cecropin）基因的表达量在过表达基因的飞蝗中为野生型飞蝗的2倍。同时，过表达基因的飞蝗感染率和死亡率明显低于野生型飞蝗。在感染后72h，过表达基因的飞蝗死亡率为35%，而野生型