基于转录组分析探究小麦抗麦红吸浆虫性状的分子机制

基于转录组分析探究小麦抗麦红吸浆虫性状的分子机制一、引言1.1研究背景与意义小麦作为全球最重要的粮食作物之一，其产量和质量直接关系到世界粮食安全。在中国，小麦的种植历史悠久，种植范围广泛，是北方地区的主要粮食作物，对保障国家粮食安全起着关键作用。然而，小麦生产面临着诸多生物和非生物胁迫的挑战，其中麦红吸浆虫（SitodiplosismosellanaGehin）是一种极具破坏力的害虫，对小麦的产量和品质构成严重威胁。麦红吸浆虫属双翅目瘿蚊科，是世界性的小麦害虫，在亚洲、欧洲和美洲的主要小麦栽培国家均有分布。在我国，麦红吸浆虫广泛分布于各大麦区，其幼虫以口器锉破小麦籽粒的果皮，吸食正在灌浆的麦粒汁液，造成小麦籽粒秕粒、空壳，严重影响小麦的产量和质量。一般情况下，麦红吸浆虫可导致小麦减产10%-30%，在严重发生年份，减产幅度可达70%以上，甚至绝收。如在2012年，渭南市小麦吸浆虫发生面积达16.1万hm²，最高虫口密度3823头/667m²，是防治指标的127倍，对当地小麦生产安全带来了严重威胁。长期以来，化学防治是控制麦红吸浆虫的主要手段之一。在小麦抽穗至开花前，使用80%敌敌畏150毫升加水4千克稀释后喷洒在25千克麦糠上拌匀，隔行每亩撒一堆，可有效防治吸浆虫成虫；也可使用40%乐果乳剂1000倍、2.5%溴氰菊酯3000倍、40%杀螟松可湿性粉剂1500倍液等进行喷雾防治。但化学防治存在诸多弊端，不仅耗费大量的人力、物力和财力，还会对环境造成污染，破坏生态平衡，导致害虫抗药性增强，同时还可能对非靶标生物产生负面影响，如影响蜜蜂等有益昆虫的生存，进而影响农作物的授粉等过程。培育和推广抗虫品种被认为是防治麦红吸浆虫最经济、有效、安全的途径。通过培育具有抗虫特性的小麦品种，能够从根本上减少害虫对小麦的侵害，降低化学农药的使用量，保护环境，同时还能保障小麦的产量和质量稳定。要实现这一目标，深入了解小麦抗麦红吸浆虫的分子机制至关重要。转录组分析作为一种强大的技术手段，能够全面揭示生物体在特定生理状态或环境条件下的基因表达谱，为研究小麦抗虫性状的分子机制提供了有力工具。转录组是指特定组织或细胞在某一发育阶段或生理状态下转录出来的所有RNA的集合，包括mRNA、rRNA、tRNA和非编码RNA等。转录组分析能够在整体水平上研究基因的表达情况，挖掘与抗虫性状相关的关键基因和代谢途径，从而深入了解小麦抗麦红吸浆虫的分子调控网络。通过转录组分析，可以发现小麦在受到麦红吸浆虫侵害时，哪些基因的表达发生了显著变化，这些基因参与了哪些生物学过程和信号通路，进而揭示小麦抗虫的分子机制。这不仅有助于我们从分子层面理解小麦与麦红吸浆虫之间的互作关系，还为小麦抗虫育种提供了重要的理论基础和基因资源。通过转录组分析，能够为小麦抗虫育种提供重要的理论依据和基因资源。筛选出与抗虫性状紧密相关的基因后，可以将这些基因作为分子标记，应用于小麦抗虫品种的选育过程中，实现分子标记辅助选择（MAS）育种。这种育种方法能够大大提高育种效率，缩短育种周期，加速抗虫小麦品种的培育和推广。还可以通过基因编辑等技术手段，对小麦的抗虫基因进行定向改良和优化，进一步增强小麦的抗虫能力。对小麦抗麦红吸浆虫性状进行转录组分析具有重要的理论和实践意义。通过深入研究小麦抗虫的分子机制，不仅能够丰富我们对植物与害虫互作关系的认识，还能为小麦抗虫育种提供关键的技术支持和基因资源，有助于培育出更多高产、优质、抗虫的小麦新品种，保障小麦生产的可持续发展和全球粮食安全。1.2小麦吸浆虫研究现状1.2.1小麦吸浆虫的危害与发生小麦吸浆虫是世界性的小麦害虫，在亚洲、欧洲和美洲的主要小麦栽培国家均有分布。在我国，小麦吸浆虫广泛分布于各大麦区，主要包括麦红吸浆虫和麦黄吸浆虫两种，其中麦红吸浆虫在平原地区的渡河两岸发生较为普遍，麦黄吸浆虫则主要发生在高原地区和高山地带。小麦吸浆虫以幼虫潜伏在颖壳内吸食正在灌浆的麦粒汁液，造成小麦籽粒秕粒、空壳，严重影响小麦的产量和质量。一般情况下，小麦吸浆虫可导致小麦减产10%-30%，在严重发生年份，减产幅度可达70%以上，甚至绝收。如在2012年，渭南市小麦吸浆虫发生面积达16.1万hm²，最高虫口密度3823头/667m²，是防治指标的127倍，对当地小麦生产安全带来了严重威胁。小麦吸浆虫的发生与多种因素密切相关。从气候因素来看，温度和湿度对其生长发育和繁殖有着显著影响。幼虫耐低温不耐高温，越冬死亡率低于越夏。越冬幼虫在10厘米土温7℃时破茧活动，12℃-15℃化蛹，20℃-23℃羽化成虫，温度上升30℃以上时，幼虫即恢复休眠。在越冬幼虫破茧活动与上升化蛹期间，雨水多（或灌溉）羽化率就高。湿度高时，不仅卵的孵化率高，且初孵幼虫活动力强，容易侵入咬害。小麦扬花前后雨水多、湿度大、气温适宜常会引起吸浆虫的大发生，而天气干旱、土壤湿度小则对其发生不利。土壤条件也是影响小麦吸浆虫发生的重要因素。壤土的土质疏松、保水力强利于其发生，粘土对其生活不利，砂土更不适宜其生活。红吸浆虫幼虫喜碱性土壤，黄吸浆吸虫喜较酸性的土壤。此外，成虫盛发期与小麦抽穗扬花期的吻合程度也会影响其危害程度，两者吻合发生重，两期错位则发生轻。1.2.2小麦吸浆虫的防治措施针对小麦吸浆虫的危害，目前主要采取化学防治、农业防治和生物防治等多种措施。化学防治是控制小麦吸浆虫的常用手段之一。在小麦播种前的最后一次浅耕时、小麦拔节期以及小麦孕穗期等关键时期，可进行土壤处理。常用药剂有2%甲基异柳磷粉剂、4.5%甲敌粉、4%敌马粉等，亩用2-3千克，或80%敌敌畏乳油50-100毫升加水1-2千克，或用50%辛硫磷乳油200毫升，加水5千克喷在20-25千克的细土上，拌匀制成毒土施用，边撒边耕，翻入土中。在小麦抽穗至开花前，可采用药剂喷雾的方式防治成虫，如每亩用80%敌敌畏150毫升，加水4千克稀释，喷洒在25千克麦糠上拌匀，隔行每亩撒一堆；或用40%乐果乳剂1000倍、2.5%溴氰菊酯3000倍、40%杀螟松可湿性粉剂1500倍液等喷雾。化学防治虽然能够在短期内迅速降低虫口密度，有效控制小麦吸浆虫的危害，但也存在诸多弊端。长期大量使用化学农药不仅耗费大量的人力、物力和财力，还会对环境造成污染，破坏生态平衡，导致害虫抗药性增强，同时还可能对非靶标生物产生负面影响，如影响蜜蜂等有益昆虫的生存，进而影响农作物的授粉等过程。农业防治措施主要包括选用抗虫品种和轮作倒茬等。不同小麦品种对小麦吸浆虫的为害程度存在差异，一般芒长多刺，口紧小穗密集，扬花期短而整齐，果皮厚的品种，对吸浆虫成虫的产卵、幼虫入侵和为害均不利。因此，可选用穗形紧密，内外颖毛长而密，麦粒皮厚，浆液不易外流的小麦品种。此外，麦田连年深翻，小麦与油菜、豆类、棉花和水稻等作物轮作，对压低虫口数量有明显的作用。在小麦吸浆虫严重田及其周围，可实行棉麦间作或改种油菜、大蒜等作物，待雨年后再种小麦，就会减轻为害。农业防治措施具有环保、可持续等优点，但抗虫品种的选育需要耗费大量的时间和精力，且轮作倒茬可能会受到土地资源和种植习惯等因素的限制。生物防治是利用生物之间的相互关系来控制害虫的方法。小麦吸浆虫的主要天敌有宽腹姬小蜂、光腹黑蜂、蚂蚁、蜘蛛等，可通过保护和利用这些天敌来控制吸浆虫的种群数量。还可以利用一些微生物制剂，如苏云金芽孢杆菌等，来防治小麦吸浆虫。生物防治具有对环境友好、不易产生抗药性等优点，但生物防治的效果受到天敌数量、环境条件等多种因素的影响，且作用效果相对较慢，在害虫爆发时可能无法迅速控制虫口密度。1.3小麦对麦红吸浆虫的抗性机制研究1.3.1避虫性避虫性是小麦抵抗麦红吸浆虫侵害的一种重要策略，主要通过自身生长特性或对环境的适应性来实现。小麦的抽穗期与麦红吸浆虫成虫羽化期的时间差是影响二者相遇几率的关键因素。若小麦抽穗期与吸浆虫成虫羽化期错开，吸浆虫就难以在小麦穗上产卵，从而减少了小麦受侵害的机会。有研究表明，早抽穗的小麦品种能有效避开吸浆虫成虫的产卵高峰期，降低受害程度。小麦的株型和叶片形态也在避虫过程中发挥作用。紧凑的株型使得吸浆虫在小麦植株上的活动空间受限，不利于其寻找合适的产卵位点和取食位置；叶片表面的蜡质层、绒毛等结构可以增加吸浆虫在叶片上爬行和停留的难度，减少其在叶片上的栖息时间。此外，小麦对环境的适应性也会影响其避虫能力。在干旱、高温等不利于吸浆虫生存的环境条件下，小麦通过调节自身的生理代谢，增强对不良环境的耐受性，同时也使得吸浆虫的生存和繁殖受到抑制，从而间接减少了自身受到侵害的风险。例如，在一些干旱地区，小麦通过根系的生长和分布变化，更好地吸收土壤中的水分，维持自身生长，而吸浆虫由于对水分条件较为敏感，在这样的环境中种群数量会显著下降。1.3.2产卵不选择性麦红吸浆虫对不同小麦品种存在明显的产卵偏好差异，这种差异与小麦的多种因素密切相关。小麦穗部的形态特征，如穗型、颖壳的形状和大小等，会影响吸浆虫的产卵选择。紧凑穗型的小麦品种，其颖壳相对紧密，吸浆虫难以将卵产在其中，而松散穗型的小麦则更容易成为吸浆虫的产卵目标。颖壳的颜色也可能对吸浆虫的产卵选择产生影响，一些研究发现，吸浆虫更倾向于在颜色较浅的颖壳上产卵。小麦植株释放的挥发性物质在吸浆虫的产卵选择过程中起着关键的信号作用。这些挥发性物质包括多种萜类、醇类、酯类等化合物，它们能够被吸浆虫感知，从而影响其行为。某些小麦品种在生长过程中会释放出对吸浆虫具有驱避作用的挥发性物质，使吸浆虫不愿意在其植株上产卵；而另一些品种释放的挥发性物质则可能对吸浆虫具有吸引作用，增加了被产卵的几率。研究表明，小麦叶片和穗部挥发物中的β-石竹烯、芳樟醇等成分的含量与吸浆虫的产卵选择性密切相关，含量较高时可降低吸浆虫的产卵偏好。此外，小麦的营养成分，如蛋白质、糖分、氨基酸等的含量和比例，也会影响吸浆虫对其产卵的选择。营养成分丰富且比例适宜的小麦品种，可能更能满足吸浆虫幼虫生长发育的需求，从而吸引吸浆虫产卵；反之，营养成分不利于吸浆虫幼虫生长的小麦品种，则会减少吸浆虫的产卵。1.3.3形态抗虫性小麦的形态特征在抵御麦红吸浆虫侵害方面发挥着重要作用。穗型是影响小麦抗虫性的重要形态因素之一，紧凑穗型的小麦品种通常具有较强的抗虫性。在紧凑穗型小麦中，小穗排列紧密，颖壳相互挤压，使得吸浆虫成虫难以在其中产卵，幼虫孵化后也难以在狭小的空间内活动和取食。相关研究表明，紧凑穗型小麦品种的吸浆虫虫口密度明显低于松散穗型品种，受害程度也相对较轻。颖壳质地对小麦抗虫性也有显著影响。质地坚硬、厚实的颖壳能够为小麦提供一道物理屏障，有效阻止吸浆虫幼虫的侵入。颖壳表面的茸毛和蜡质层也能起到一定的防御作用，茸毛可以干扰吸浆虫的爬行和附着，蜡质层则增加了颖壳表面的光滑度，使吸浆虫难以停留和产卵。有研究发现，颖壳茸毛较多且蜡质层较厚的小麦品种，吸浆虫的产卵量明显减少，幼虫的侵入率也较低。芒的存在与否及芒的长度、形态等也与小麦的抗虫性有关。长芒小麦品种在一定程度上能够阻碍吸浆虫成虫在麦穗上的活动，减少其产卵机会。芒上的倒刺等结构可以对吸浆虫造成物理伤害，使其不敢轻易靠近麦穗。有研究表明，长芒小麦品种的吸浆虫成虫在麦穗上的停留时间明显短于短芒或无芒品种，产卵量也相应减少。1.3.4生化抗虫性小麦体内的生化物质和代谢途径在抗虫过程中发挥着核心作用，这些生化机制是小麦抵御麦红吸浆虫侵害的重要防线。小麦在受到麦红吸浆虫侵害时，会迅速启动一系列防御反应，其中次生代谢产物的合成和积累是重要的防御手段之一。酚类物质是小麦体内一类重要的次生代谢产物，包括黄酮类、单宁等。这些酚类物质具有抗菌、抗氧化和抗虫等多种生物活性。当吸浆虫取食小麦时，酚类物质可以与吸浆虫体内的蛋白质结合，降低蛋白质的活性，影响吸浆虫的消化和生长发育。研究表明，受吸浆虫侵害后，小麦叶片和穗部的酚类物质含量显著增加，且抗虫品种中酚类物质的积累速度和含量均高于感虫品种。茉莉酸（JA）和水杨酸（SA）信号通路在小麦抗虫反应中起着关键的调控作用。茉莉酸信号通路主要参与植物对咀嚼式口器害虫的防御反应，当小麦受到吸浆虫侵害时，茉莉酸信号通路被激活，诱导一系列防御基因的表达，这些基因编码的蛋白质参与了次生代谢产物的合成、细胞壁的加固等防御过程。水杨酸信号通路则主要参与植物对病原菌的防御反应，但在小麦抗虫过程中也发挥着一定的作用，它可以与茉莉酸信号通路相互作用，协同调控小麦的抗虫反应。研究发现，通过外源施加茉莉酸或水杨酸，可以增强小麦对吸浆虫的抗性，提高防御基因的表达水平。此外，小麦体内的一些酶类物质也在抗虫过程中发挥着重要作用。如过氧化物酶（POD）、多酚氧化酶（PPO）和超氧化物歧化酶（SOD）等，这些酶参与了小麦体内的氧化还原反应，能够清除吸浆虫侵害过程中产生的活性氧（ROS），保护小麦细胞免受氧化损伤。它们还可以通过催化酚类物质的氧化聚合，形成具有抗菌和抗虫活性的醌类物质，增强小麦的抗虫能力。研究表明，在吸浆虫侵害过程中，小麦体内的POD、PPO和SOD等酶的活性显著升高，且抗虫品种中这些酶的活性变化更为明显。1.4小麦对麦红吸浆虫抗性遗传研究小麦对麦红吸浆虫的抗性遗传是一个复杂的过程，涉及多个基因和遗传因素的相互作用。研究表明，小麦对麦红吸浆虫的抗性受多基因控制，表现为数量性状遗传。这种遗传方式使得抗性表现呈现出连续的变异，受到环境因素的影响也较大。在遗传研究中，通过构建遗传群体，如F2代、重组自交系（RIL）等，利用分子标记技术，如简单序列重复（SSR）、单核苷酸多态性（SNP）等，对小麦抗麦红吸浆虫基因进行定位和分析。有研究利用“冀麦24×石4185”组合的F2:3家系，通过抗虫鉴定与分子检测，对冀麦24抗麦红吸浆虫基因标记进行验证，发现引物Xgwm537、Xgwm400、Xgwm46、Xbarc176表现出稳定的多态性，并通过区间作图法检测到QTL位点，该位点加性效应为16.83，显性效应为-17.31，最大似然值为76.7，可以解释22%的变异率，抗虫基因距离Xgwm400约5.0cM，距离Xgwm46约21.0cM。小麦抗麦红吸浆虫性状还与其他农艺性状存在一定的相关性。对“冀麦24×石4185”组合的F3:4家系的研究发现，穗长、千粒重和旗叶面积与抗虫标记达到显著相关，穗长与三个标记的相关系数分别为-0.17、-0.18和-0.17，都达到极显著水平；旗叶面积与标记Xgwm537和Xgwm400呈显著相关，相关系数均为-0.12；千粒重与标记Xgwm46呈显著相关，相关系数为-0.12。这些相关性的发现，为小麦抗虫育种过程中综合考虑多种性状提供了重要依据，在选育抗虫品种时，需要兼顾这些农艺性状，以确保品种的高产、优质特性不受影响。1.5小麦对麦红吸浆虫抗性基因/QTL的定位研究在小麦抗麦红吸浆虫的研究领域，对相关抗性基因和数量性状位点（QTL）的定位探索是关键环节。借助分子标记技术与遗传群体构建，科研人员已在该方面取得了一定成果。在染色体层面，多个抗性基因/QTL位点被陆续定位。在2B染色体上，有研究利用“河农215×石新539”组合的F2分离世代，通过构建极端类型池，对346对均匀分布于小麦21个连锁群的SSR引物展开筛选，成功筛选出Xgwm501、Wmc332、Wmc149这3对多态性引物，由此确定2B染色体上存在抗虫基因。后续通过区间作图法，检测到一个与河农215抗麦红吸浆虫性状紧密相关的QTL位点，该位点加性效应为0.14，显性效应为-0.08，最大似然值达83.0，可解释23%的变异率，抗虫基因位点距离Xgwm501约12cM，距离Wmc332约11.6cM。在6B染色体上也有重要发现，科研人员利用“冀麦24×石4185”组合的F2:3家系，通过抗虫鉴定与分子检测，对冀麦24抗麦红吸浆虫基因标记进行可用性验证时，检测到相关QTL位点。该位点加性效应为16.83，显性效应为-17.31，最大似然值为76.7，能够解释22%的变异率，抗虫基因距离Xgwm400约5.0cM，距离Xgwm46约21.0cM，且引物Xgwm537、Xgwm400、Xgwm46、Xbarc176在检测中表现出稳定的多态性，可作为该基因的辅助选择标记。不同位点的抗性基因/QTL在效应上存在差异。部分位点的基因效应表现为加性效应，对小麦的抗虫性贡献较为稳定，能在不同遗传背景下持续发挥作用，使小麦在面对麦红吸浆虫侵害时保持一定抗性水平；有的位点则存在显性效应，其作用的发挥与等位基因的显隐性关系相关，显性基因的存在可能会显著增强小麦的抗虫能力，隐性基因在纯合状态下也可能对小麦抗虫性产生独特影响。这些效应的差异使得小麦抗虫性呈现出复杂的遗传模式，也为抗虫育种提供了多样化的基因资源和选择策略。1.6植物抗虫相关基因研究植物在长期的进化过程中，形成了一系列复杂而精妙的抗虫机制，其中抗虫相关基因发挥着关键作用。这些基因参与了植物抗虫的各个环节，从感知害虫侵害信号到启动防御反应，为植物抵御害虫提供了重要的分子基础。根据其功能和作用方式，植物抗虫相关基因可大致分为以下几类。首先是蛋白酶抑制剂基因，这类基因编码的蛋白酶抑制剂能够抑制昆虫体内蛋白酶的活性。当害虫取食含有蛋白酶抑制剂的植物组织后，其消化过程会受到严重阻碍，无法正常分解和吸收食物中的蛋白质，进而影响其生长发育，甚至导致死亡。常见的蛋白酶抑制剂有丝氨酸蛋白酶抑制剂、半胱氨酸蛋白酶抑制剂等。在小麦抗麦红吸浆虫的研究中，若小麦体内的蛋白酶抑制剂基因能够高效表达，可能会对吸浆虫幼虫的消化功能产生抑制作用，减少其对小麦籽粒的侵害。有研究表明，将外源蛋白酶抑制剂基因导入植物后，植物对害虫的抗性得到了显著增强，这为小麦抗虫育种提供了新的思路。凝集素基因也是一类重要的抗虫基因，其编码的凝集素可以与昆虫肠道上皮细胞表面的糖蛋白或糖脂特异性结合。这种结合会破坏肠道细胞的正常结构和功能，影响昆虫的营养吸收，还可能引发昆虫的免疫反应，对其生长和存活产生负面影响。在小麦中，某些凝集素基因的表达水平与小麦对麦红吸浆虫的抗性可能存在关联。研究发现，一些凝集素能够与麦红吸浆虫肠道中的特定受体结合，干扰其生理活动，从而降低吸浆虫对小麦的危害。萜类合成酶基因在植物抗虫中也扮演着重要角色，它参与萜类化合物的合成。萜类化合物是一类种类繁多、结构复杂的次生代谢产物，许多萜类化合物具有挥发性，能够吸引害虫的天敌，起到间接防御的作用；有些萜类化合物则对害虫具有直接的毒性，可抑制害虫的生长和繁殖。在小麦受到麦红吸浆虫侵害时，萜类合成酶基因可能被诱导表达，促使小麦合成更多的萜类化合物，以抵御吸浆虫的侵害。例如，某些倍半萜类化合物能够对麦红吸浆虫产生驱避作用，减少其在小麦植株上的产卵和取食。植物激素信号转导相关基因在植物抗虫过程中起着关键的调控作用。茉莉酸（JA）、水杨酸（SA）和乙烯（ET）等植物激素信号通路相互交织，共同调节植物的抗虫反应。茉莉酸信号通路主要参与植物对咀嚼式口器害虫的防御，当小麦受到麦红吸浆虫侵害时，茉莉酸信号通路被激活，诱导一系列防御基因的表达，这些基因参与了次生代谢产物的合成、细胞壁的加固等防御过程。水杨酸信号通路主要与植物对病原菌的防御有关，但在小麦抗虫中也有一定作用，它与茉莉酸信号通路相互作用，协同调控小麦的抗虫反应。乙烯信号通路则在植物的生长发育和逆境响应中发挥重要作用，也参与了植物的抗虫过程。研究表明，通过调节这些植物激素信号转导相关基因的表达，可以增强小麦对麦红吸浆虫的抗性。1.7生物信息学研究相关数据库在小麦抗麦红吸浆虫性状的转录组分析中，一系列生物信息学数据库和工具发挥着不可或缺的作用，它们为深入挖掘转录组数据背后的生物学意义提供了有力支持。NCBI（NationalCenterforBiotechnologyInformation）数据库是生物学领域最为重要的数据库之一，它涵盖了丰富的生物数据资源。其中的GenBank数据库保存了海量的核酸序列数据，在小麦转录组分析中，通过将测序得到的小麦转录本序列与GenBank中的已知序列进行比对，可以确定转录本的同源性，判断其是否为已知基因的同源序列，从而初步了解转录本的功能。例如，若比对结果显示某转录本与GenBank中已知的抗虫基因具有高度同源性，那么该转录本很可能在小麦抗麦红吸浆虫过程中发挥着类似的抗虫作用。SRA（SequenceReadArchive）数据库则存储了大量的原始测序数据，这些数据可作为参考，用于验证和补充本研究中的转录组测序结果，确保数据的准确性和可靠性。EnsemblPlants数据库专注于植物基因组的注释和分析，对于小麦转录组研究具有重要价值。它提供了小麦基因组的详细注释信息，包括基因结构、外显子-内含子边界、基因的功能注释等。通过该数据库，能够准确地确定转录本在小麦基因组中的位置，了解其所在基因的上下游调控区域，进而深入分析转录本的调控机制。还能利用EnsemblPlants数据库中的基因家族信息，研究小麦抗虫相关基因所属的基因家族，分析基因家族成员之间的进化关系和功能分化，为揭示小麦抗虫的分子进化机制提供线索。Uniprot数据库是蛋白质序列和功能信息的重要存储库，在小麦转录组分析中，可将预测得到的小麦转录本编码的蛋白质序列与Uniprot中的蛋白质序列进行比对，获取蛋白质的功能注释、结构域信息、翻译后修饰位点等。这些信息对于深入了解转录本编码蛋白质的生物学功能至关重要，有助于确定该蛋白质在小麦抗麦红吸浆虫信号通路或代谢途径中的具体作用。例如，若发现某蛋白质含有与防御反应相关的结构域，那么可以推测其编码基因可能参与小麦的抗虫防御过程。GO（GeneOntology）数据库为基因功能的分类和注释提供了统一的标准和框架。在小麦抗麦红吸浆虫转录组分析中，通过GO富集分析，可以将差异表达基因按照生物过程、细胞组分和分子功能进行分类，从而全面了解这些基因在小麦抗虫过程中参与的生物学过程和发挥的功能。在生物过程分类中，可能会发现差异表达基因显著富集在“防御反应”“信号转导”“次生代谢产物合成”等生物学过程，这表明这些过程在小麦抗麦红吸浆虫过程中起到关键作用；在分子功能分类中，若某些基因富集在“酶活性”“转录因子活性”等功能类别，可进一步研究这些基因编码的蛋白质在调控抗虫相关代谢途径或基因表达中的作用。KEGG（KyotoEncyclopediaofGenesandGenomes）数据库是整合了基因组、化学和系统功能信息的数据库，尤其在代谢途径分析方面具有重要作用。通过KEGG富集分析，可以确定差异表达基因参与的代谢途径和信号通路。在小麦抗麦红吸浆虫研究中，可能会发现一些差异表达基因参与了植物激素信号转导通路、次生代谢产物合成途径等，如茉莉酸信号通路、萜类化合物合成途径等，这些通路与小麦的抗虫反应密切相关，进一步研究这些通路中关键基因的表达变化和调控机制，有助于揭示小麦抗麦红吸浆虫的分子机制。1.8混池转录组测序(BSR-Seq)在基因定位中的研究进展1.8.1集群分离分析法(BSA)集群分离分析法（BulkedSegregantAnalysis，BSA），也被称作分离群体分组分析法或混合分组分析法，是一种在基因定位研究中广泛应用的高效策略。其核心原理基于遗传群体中目标性状的分离特性，通过从分离群体里选取具有极端表型的个体，构建两个基因池，即性状表现为显性的“显性池”和性状表现为隐性的“隐性池”。这两个基因池在目标性状上呈现出显著差异，而在其他遗传背景上尽可能相似。由于基因池内个体的目标性状一致，使得与目标性状紧密连锁的基因在两个基因池中出现的频率产生明显差异，从而能够借助分子标记技术对这些差异进行检测，进而实现对目标基因的初步定位。在小麦抗麦红吸浆虫的研究中，BSA方法展现出了重要的应用价值。科研人员可以从大量的小麦遗传群体中挑选出抗虫性表现突出的小麦个体组成抗虫池，以及感虫性明显的个体组成感虫池。通过对这两个基因池进行分子标记分析，如利用简单序列重复（SSR）标记、单核苷酸多态性（SNP）标记等，能够快速筛选出与小麦抗麦红吸浆虫性状紧密连锁的分子标记。这些标记就如同基因的“路标”，为后续进一步精细定位和克隆抗虫基因指明方向，极大地提高了基因定位的效率，减少了研究的盲目性。1.8.2转录组测序(RNA-Seq)转录组测序（RNA-Seq）作为一项前沿的分子生物学技术，正深刻地改变着基因研究的格局。其技术原理是基于新一代高通量测序平台，首先将细胞内的RNA逆转录为cDNA，接着对这些cDNA进行大规模测序，从而全面、快速地获取特定组织或细胞在某一发育阶段或生理状态下的所有转录本信息。与传统的基因表达分析技术相比，RNA-Seq具有无可比拟的优势。RNA-Seq具有超高的检测灵敏度，能够检测到低丰度表达的基因。在小麦抗麦红吸浆虫的研究中，一些参与抗虫反应的基因可能仅在特定时期或特定组织中微量表达，传统技术难以捕捉到它们的存在，而RNA-Seq却能够精准地识别和定量这些低丰度基因，为全面揭示小麦抗虫的分子机制提供了可能。RNA-Seq的覆盖度广，能够对整个转录组进行无遗漏的扫描，无论是编码蛋白质的mRNA，还是具有重要调控功能的非编码RNA，都能被纳入研究范围，有助于发现新的抗虫相关基因和调控元件。它还具有强大的定量准确性，通过对测序reads的计数和分析，可以精确地测定不同基因的表达水平，从而准确地反映基因在不同条件下的表达变化，为深入研究小麦抗虫过程中的基因表达调控网络提供了坚实的数据基础。1.8.3混池转录组测序(BSR-Seq)混池转录组测序（BulkedSegregantRNA-Seq，BSR-Seq）巧妙地将集群分离分析法（BSA）和转录组测序（RNA-Seq）的优势融合在一起，开创了基因定位研究的新局面。其原理是在遗传群体中，针对目标性状构建两个极端表型的混池，如抗虫混池和感虫混池，然后对这两个混池分别进行转录组测序。通过对测序数据的深度挖掘和生物信息学分析，能够快速、准确地筛选出与目标性状紧密关联的差异表达基因和分子标记。在小麦抗麦红吸浆虫的研究领域，BSR-Seq技术展现出了独特的应用潜力。通过对小麦抗虫混池和感虫混池进行转录组测序，不仅可以全面地了解小麦在抗虫和感虫状态下基因表达的差异，还能借助这些差异信息，迅速定位到与抗麦红吸浆虫性状相关的关键基因和调控通路。这些关键基因和通路的揭示，将为深入理解小麦抗虫的分子机制提供关键线索，为小麦抗虫品种的选育提供有力的基因资源和理论依据。1.8.4BSR-Seq在小麦上的应用研究进展近年来，BSR-Seq技术在小麦基因定位和抗虫研究中取得了一系列令人瞩目的成果。在小麦叶锈病抗性基因的研究中，科研人员运用BSR-Seq技术，成功地定位到了多个与叶锈病抗性紧密相关的基因位点。通过对这些基因位点的功能分析，发现它们参与了小麦对叶锈菌的识别、信号传导以及防御反应的激活等多个关键过程，为培育抗叶锈病的小麦新品种提供了重要的基因靶点。在小麦抗蚜虫研究方面，BSR-Seq技术也发挥了重要作用。研究人员利用该技术筛选出了多个与抗蚜虫性状相关的差异表达基因，这些基因涉及到植物次生代谢产物的合成、激素信号传导以及细胞壁的加固等多个生物学过程。进一步的研究表明，通过调控这些基因的表达，可以显著增强小麦对蚜虫的抗性，为小麦抗蚜育种提供了新的思路和方法。在小麦抗麦红吸浆虫的研究中，虽然BSR-Seq技术的应用尚处于起步阶段，但已展现出了巨大的潜力。通过对小麦抗虫和感虫混池进行转录组测序，有望挖掘出更多与抗麦红吸浆虫性状相关的关键基因和分子标记，为深入揭示小麦抗虫的分子机制和培育抗虫新品种奠定坚实的基础。相信在未来，随着技术的不断发展和完善，BSR-Seq技术将在小麦抗虫研究中发挥更加重要的作用，为保障小麦的安全生产提供强有力的技术支持。1.9小麦分子标记开发研究进展1.9.1分子标记概述分子标记是以个体间遗传物质内核苷酸序列变异为基础的遗传标记，是DNA水平遗传多态性的直接反映。与传统的形态标记、细胞学标记和生化标记相比，分子标记具有诸多独特的优势。从稳定性角度来看，分子标记不受环境因素和发育阶段的影响。传统的形态标记易受环境条件如温度、光照、水分等的影响，同一基因型在不同环境下可能表现出不同的形态特征；细胞学标记需要进行细胞学观察，操作复杂且受细胞发育时期的限制；生化标记则依赖于特定的酶活性或蛋白质表达，也会受到环境和生理状态的干扰。而分子标记基于DNA序列，其遗传信息稳定，无论在植物的任何生长阶段，还是处于何种环境条件下，都能准确地反映个体的遗传特征。在多态性方面，分子标记具有丰富的多态性。由于DNA序列的多样性，分子标记能够检测到大量的遗传变异位点，这些变异位点为遗传分析提供了丰富的信息，使得研究人员能够更准确地识别不同个体之间的遗传差异，在小麦品种鉴定、遗传多样性分析等方面具有重要应用价值。此外，分子标记的检测方便快捷。随着现代生物技术的不断发展，分子标记的检测方法日益成熟，如聚合酶链式反应（PCR）技术的广泛应用，使得分子标记的检测过程变得高效、灵敏，能够在短时间内对大量样本进行分析，大大提高了研究效率。常见的分子标记类型包括限制性片段长度多态性（RFLP）、随机扩增多态性DNA（RAPD）、简单序列重复（SSR）、单核苷酸多态性（SNP）等。RFLP是最早发展起来的分子标记技术，它利用限制性内切酶切割不同个体的DNA，由于DNA序列的差异，酶切后产生的片段长度不同，通过凝胶电泳和Southern杂交技术来检测这些差异，从而揭示个体间的遗传多态性。虽然RFLP具有稳定性高、结果可靠等优点，但操作过程繁琐，需要使用放射性同位素标记探针，对实验条件和操作人员的要求较高，且检测成本也相对较高，因此在实际应用中受到一定限制。RAPD则是利用随机引物对基因组DNA进行PCR扩增，由于不同个体的DNA序列存在差异，引物与模板DNA的结合位点和扩增产物的长度也会有所不同，通过聚丙烯酰胺凝胶电泳检测扩增产物的多态性，进而分析个体间的遗传关系。RAPD技术操作简单、快速，不需要预先了解DNA序列信息，成本较低，但该技术的稳定性较差，重复性不好，结果的可靠性相对较低，容易受到实验条件如引物浓度、退火温度等的影响。SSR，也被称为微卫星DNA，是由2-6个核苷酸组成的简单重复序列，广泛分布于真核生物基因组中。由于重复单位的重复次数在不同个体间存在差异，使得SSR具有高度的多态性。通过设计特异性引物对SSR位点进行PCR扩增，然后利用凝胶电泳或毛细管电泳检测扩增产物的长度多态性，即可分析个体间的遗传差异。SSR标记具有多态性高、重复性好、共显性遗传等优点，在小麦遗传图谱构建、基因定位、品种鉴定等方面得到了广泛应用。SNP是指在基因组水平上由单个核苷酸的变异所引起的DNA序列多态性，是目前最丰富的分子标记类型。SNP广泛存在于基因组中，具有密度高、稳定性好、易于自动化检测等优点。随着高通量测序技术的发展，SNP的检测变得更加高效、准确，能够在全基因组范围内大规模地检测SNP位点，为小麦的遗传研究和分子育种提供了海量的数据支持，在小麦的全基因组关联分析（GWAS）、分子标记辅助选择（MAS）等研究中发挥着重要作用。1.9.2小麦不同类型功能标记开发研究进展在小麦抗虫研究中，不同类型的分子标记发挥着各自独特的作用，其开发和应用不断推动着小麦抗虫研究的深入发展。SSR标记在小麦抗虫基因定位和标记辅助选择中应用广泛。科研人员利用SSR标记技术，对小麦抗麦红吸浆虫基因进行定位研究。如利用“河农215×石新539”组合的F2分离世代，通过构建极端类型池，对346对均匀分布于小麦21个连锁群的SSR引物进行筛选，成功筛选出Xgwm501、Wmc332、Wmc149这3对多态性引物，确定2B染色体上存在抗虫基因，并通过区间作图法检测到一个与河农215抗麦红吸浆虫性状密切相关的QTL位点。在小麦抗蚜虫研究中，也有研究利用SSR标记筛选出与抗蚜性状相关的分子标记，为小麦抗蚜育种提供了重要的技术支持。SSR标记的优势在于其多态性高、重复性好，能够准确地检测到遗传变异，为小麦抗虫基因的定位和追踪提供了可靠的工具。但SSR标记的开发需要进行引物设计和筛选，工作量较大，且检测过程相对复杂，成本较高。SNP标记近年来在小麦抗虫研究中逐渐崭露头角。随着高通量测序技术的飞速发展，SNP标记的检测变得更加高效、准确，能够在全基因组范围内快速筛选出与抗虫性状相关的SNP位点。通过对小麦抗虫品种和感虫品种进行全基因组重测序，分析SNP位点的差异，能够挖掘出潜在的抗虫基因和分子标记。在小麦抗叶锈病研究中，利用SNP标记进行全基因组关联分析，成功鉴定出多个与叶锈病抗性相关的SNP位点，为小麦抗叶锈病基因的克隆和功能研究奠定了基础。SNP标记具有密度高、稳定性好、易于自动化检测等优点，能够实现大规模的基因分型和遗传分析，大大提高了小麦抗虫研究的效率和准确性。但SNP标记的检测需要专业的测序设备和数据分析软件，对技术人员的要求较高，且数据分析工作量巨大。基于转录组测序开发的功能标记也为小麦抗虫研究提供了新的思路和方法。通过对小麦在受到麦红吸浆虫侵害前后的转录组进行测序，分析差异表达基因，从中筛选出与抗虫性状相关的基因，并开发相应的分子标记。这些标记能够直接反映基因的表达变化，更准确地揭示小麦抗虫的分子机制。有研究通过转录组测序，筛选出小麦抗麦红吸浆虫相关的差异表达基因，基于这些基因开发了分子标记，为小麦抗虫品种的选育提供了新的基因资源和标记辅助选择方法。基于转录组测序开发的功能标记具有针对性强、能够反映基因表达水平变化等优点，但转录组测序成本较高，数据分析复杂，需要专业的生物信息学知识和技能。1.10技术路线本研究旨在通过转录组分析深入探究小麦抗麦红吸浆虫性状的分子机制，技术路线如图1-1所示：图1-1技术路线图首先，进行实验材料的准备。选择具有明显抗虫差异的小麦品种作为实验材料，将这些小麦种子播种于温室或试验田中，给予适宜的生长环境，包括控制温度、光照、水分和养分等条件，确保小麦植株能够正常生长发育。在小麦生长至合适阶段，人工接入麦红吸浆虫，使小麦受到吸浆虫的侵害。同时设置未接虫的对照组，以对比分析接虫前后小麦基因表达的差异。在小麦与麦红吸浆虫互作的关键时期，分别采集抗虫小麦和感虫小麦的叶片、穗部等组织样本。为保证实验结果的准确性和可靠性，每个样本设置多个生物学重复。将采集后的样本迅速放入液氮中冷冻处理，以抑制样本中RNA酶的活性，防止RNA降解，随后将样本转移至-80℃冰箱中保存，以备后续RNA提取使用。对保存的样本进行总RNA的提取工作，选用高质量的RNA提取试剂盒，严格按照操作说明进行提取，确保提取的RNA完整性好、纯度高。提取完成后，利用琼脂糖凝胶电泳检测RNA的完整性，观察RNA条带的清晰度和完整性；通过分光光度计测定RNA的浓度和纯度，确保A260/A280比值在1.8-2.0之间，A260/A230比值大于2.0，以保证RNA质量符合后续实验要求。将合格的RNA样本送往专业的测序公司，采用Illumina测序平台进行转录组测序。测序过程中，首先将RNA逆转录为cDNA，然后构建cDNA文库，对文库进行质量检测合格后，进行高通量测序，从而获得大量的原始测序数据。对测序得到的原始数据进行严格的质量控制。利用FastQC等软件对原始数据进行质量评估，查看数据的碱基质量分布、测序错误率、GC含量等指标。使用Trimmomatic等工具去除低质量reads、接头序列和污染序列，经过质量控制后，得到高质量的cleanreads，为后续的数据分析奠定基础。将高质量的cleanreads与参考基因组进行比对，使用Hisat2等比对软件，确定reads在基因组上的位置，分析基因的表达水平，计算每个基因的表达量，常用的量化指标为FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）值。通过比较抗虫小麦和感虫小麦之间基因表达量的差异，筛选出差异表达基因，设定筛选标准为|log2(FoldChange)|≥1且FDR（FalseDiscoveryRate）≤0.05。对筛选出的差异表达基因进行功能注释和富集分析。借助GO（GeneOntology）数据库，从生物过程、细胞组分和分子功能三个层面，对差异表达基因进行功能分类和注释，明确基因在小麦抗虫过程中参与的生物学过程和发挥的功能；利用KEGG（KyotoEncyclopediaofGenesandGenomes）数据库，分析差异表达基因参与的代谢途径和信号转导通路，找出与小麦抗麦红吸浆虫相关的关键代谢途径和信号通路。通过实时荧光定量PCR（qRT-PCR）技术对转录组测序结果进行验证。根据差异表达基因的序列信息，设计特异性引物，以小麦的看家基因作为内参基因，对部分差异表达基因进行qRT-PCR检测。将qRT-PCR检测结果与转录组测序得到的基因表达量数据进行对比分析，验证转录组测序结果的准确性和可靠性。基于转录组分析结果，结合前人研究，筛选出与小麦抗麦红吸浆虫性状高度相关的关键基因。对这些关键基因进行克隆和功能验证，采用基因编辑技术如CRISPR/Cas9对关键基因进行敲除或过表达，观察小麦在基因编辑后的抗虫表型变化，深入研究这些基因在小麦抗虫过程中的作用机制，为小麦抗虫育种提供理论依据和基因资源。二、材料与方法2.1试验材料本研究选用了具有明显抗虫差异的两个小麦品种，分别为抗虫品种‘河农215’和感虫品种‘石新539’。‘河农215’是河北农业大学选育的小麦品种，前期研究及田间种植表现显示，其对麦红吸浆虫具有较强的抗性，在麦红吸浆虫高发地区种植时，受害程度较轻，产量损失较小；‘石新539’由石家庄市农林科学研究院培育，对麦红吸浆虫较为敏感，在相同的虫害发生条件下，容易受到麦红吸浆虫的侵害，导致籽粒瘪粒、空壳等问题，产量明显下降。麦红吸浆虫采自河北省保定市某小麦种植田，该地区麦红吸浆虫常年发生，虫口密度较高。采集时选取具有代表性的地块，采用五点取样法，每个样点随机采集10-15个麦穗，将采集到的麦穗放入密封袋中，带回实验室后，在解剖镜下仔细挑取麦红吸浆虫成虫，放入饲养盒中备用。饲养盒内放置湿润的棉花以保持湿度，同时提供新鲜的小麦穗作为成虫的食物和产卵基质，温度控制在20-23℃，光照周期为16h光照/8h黑暗，以模拟自然环境，保证麦红吸浆虫的正常生长和繁殖。试验在河北农业大学的试验田中进行，该试验田地势平坦，土壤为壤土，肥力中等且均匀，pH值约为7.2，保水保肥能力良好，历年均有麦红吸浆虫发生，能够为试验提供自然的虫害压力。试验田周边种植有其他小麦品种作为保护行，以减少外界因素对试验的干扰，确保试验结果的准确性和可靠性。在小麦生长期间，严格按照常规的田间管理措施进行操作，包括适时浇水、施肥、中耕除草等，确保小麦植株能够正常生长发育，同时避免使用任何杀虫剂，以保证麦红吸浆虫对小麦的自然侵害。2.2试验方法2.2.1吸浆虫发育情况调查在试验田中，采用五点取样法，每个样点设置1个1m×1m的样方，从3月中下旬小麦拔节期开始，每隔3天对样方内的土壤进行取样，每个样方取5个土样，土样深度为0-20cm，将土样放入筛网中，用清水冲洗，收集其中的吸浆虫幼虫和蛹，在解剖镜下观察其形态特征，判断发育阶段，并记录数量。在小麦抽穗期至扬花期，每天傍晚采用网捕法调查吸浆虫成虫数量，使用口径为30cm的捕虫网，在每个样方内的麦穗上方10-20cm高度，以“∞”字形挥动捕虫网，每次挥动10次，记录捕获的成虫数量，计算平均虫口密度。同时，观察成虫的羽化、交尾和产卵行为，记录相关数据。2.2.2供试小麦材料的抗虫性鉴定在小麦灌浆期，每个小麦品种选取30个麦穗，用剪刀将麦穗剪下，放入密封袋中带回实验室。在解剖镜下，仔细观察麦穗上的吸浆虫着卵量、幼虫数量以及籽粒的受害情况。统计每穗的着卵数、幼虫数，计算平均单穗着卵量和平均单穗幼虫数；观察籽粒是否出现瘪粒、空壳等受害症状，统计受害籽粒数，计算虫害损失率，公式如下：虫害损失率(\%)=\frac{受害籽粒数}{总籽粒数}\times100\%根据虫害损失率，参考相关标准对小麦品种的抗虫性进行分级评价，具体分级标准如下：抗性等级虫害损失率范围高抗（HR）虫害损失率＜10%中抗（MR）10%≤虫害损失率＜30%低抗（LR）30%≤虫害损失率＜50%感虫（S）50%≤虫害损失率＜70%高感（HS）虫害损失率≥70%2.2.3RNA提取与转录组测序分别采集抗虫品种‘河农215’和感虫品种‘石新539’在接虫后24h、48h和72h的叶片和穗部组织，每个时间点每个品种每个组织采集3个生物学重复，迅速放入液氮中速冻，然后保存于-80℃冰箱备用。使用TRIzol试剂（Invitrogen公司）提取总RNA，具体步骤如下：取约100mg组织样品，加入1mLTRIzol试剂，在液氮中研磨成粉末状，室温静置5min，使组织充分裂解；加入0.2mL氯仿，剧烈振荡15s，室温静置3min；4℃、12000rpm离心15min，将上层水相转移至新的离心管中；加入等体积的异丙醇，轻轻混匀，室温静置10min；4℃、12000rpm离心10min，弃上清，得到RNA沉淀；用75%乙醇洗涤RNA沉淀2次，4℃、7500rpm离心5min，弃上清；室温晾干RNA沉淀，加入适量的DEPC水溶解RNA。利用NanoDrop2000超微量分光光度计（ThermoFisherScientific公司）测定RNA的浓度和纯度，确保A260/A280比值在1.8-2.0之间，A260/A230比值大于2.0；使用Agilent2100生物分析仪（AgilentTechnologies公司）检测RNA的完整性，RIN值大于7.0。将合格的RNA样品送往北京诺禾致源科技股份有限公司，采用IlluminaHiSeq平台进行转录组测序。首先将RNA逆转录为cDNA，然后构建cDNA文库，文库构建过程包括末端修复、加A尾、连接测序接头、PCR扩增等步骤；对文库进行质量检测，包括文库浓度测定、插入片段大小检测等；最后将合格的文库进行高通量测序，测序策略为双端测序（Paired-End），测序读长为150bp。2.2.4测序质量监控与比对分析使用FastQC软件对测序得到的原始数据进行质量评估，分析碱基质量分布、测序错误率、GC含量、序列重复率等指标。使用Trimmomatic软件去除低质量reads（质量值Q＜20的碱基占比超过50%的reads）、接头序列以及含有N（未知碱基）比例超过10%的reads，得到高质量的cleanreads。将cleanreads与中国春小麦参考基因组（IWGSCRefSeqv1.0）进行比对，使用Hisat2软件，参数设置为默认参数。比对完成后，统计比对率、唯一比对率、多比对率等指标，评估比对效果。利用Samtools软件将比对结果转换为BAM格式文件，并进行排序和索引，以便后续分析。2.2.5差异表达分析使用StringTie软件对每个样本的比对结果进行转录本组装和定量，计算每个基因的表达量，以FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）值表示。利用DESeq2软件进行差异表达分析，筛选抗虫品种‘河农215’和感虫品种‘石新539’之间的差异表达基因。设定筛选标准为|log2(FoldChange)|≥1且FDR（FalseDiscoveryRate）≤0.05，其中FoldChange为抗虫品种与感虫品种中基因表达量的比值，FDR是通过Benjamini-Hochberg方法对P值进行校正得到的，用于控制多重检验中的假阳性率。2.2.6BSR-Seq分析从抗虫品种‘河农215’和感虫品种‘石新539’的F2代群体中，分别选取30株抗虫性最强和30株感虫性最强的植株，将其叶片混合，构建抗虫池和感虫池。对两个基因池进行转录组测序，测序方法同2.2.3。对测序数据进行质量控制和比对分析，方法同2.2.4。使用SNP-calling软件（如GATK）进行单核苷酸多态性（SNP）检测，筛选出在抗虫池和感虫池中频率差异显著的SNP位点。利用这些SNP位点进行关联分析，计算每个SNP位点的关联值（如Δ(SNP-index)），确定与小麦抗麦红吸浆虫性状紧密关联的SNP位点和候选基因区域。2.2.7差异表达基因的GO分析和Pathway富集分析利用DAVID（DatabaseforAnnotation,VisualizationandIntegratedDiscovery）在线工具对差异表达基因进行GO（GeneOntology）功能注释和富集分析。GO注释从生物过程（BiologicalProcess）、细胞组分（CellularComponent）和分子功能（MolecularFunction）三个层面进行，分析差异表达基因在哪些生物学过程、细胞结构和分子功能中显著富集。设定富集分析的阈值为P≤0.05，筛选出显著富集的GOterms，绘制GO富集分析柱状图和气泡图，直观展示差异表达基因的功能分布。利用KEGG（KyotoEncyclopediaofGenesandGenomes）数据库对差异表达基因进行Pathway富集分析，确定差异表达基因参与的代谢途径和信号转导通路。使用KOBAS软件进行分析，设定富集分析的阈值为P≤0.05，筛选出显著富集的KEGGpathways，绘制KEGGpathway富集分析气泡图，展示差异表达基因在不同代谢途径和信号通路中的富集情况。2.2.8同源基因分析将筛选出的差异表达基因序列提交到NCBI的BLAST工具中，选择nr（非冗余蛋白质序列数据库）进行比对，设置比对参数为默认值。根据比对结果，选取同源性较高（Identity≥80%，E-value≤1e-10）的基因作为同源基因，分析其在其他物种中的功能和作用机制。利用MEGA软件构建差异表达基因及其同源基因的系统进化树，采用邻接法（Neighbor-Joiningmethod），bootstrap值设置为1000，通过分析进化树的分支结构和遗传距离，了解差异表达基因与同源基因之间的进化关系，推测其功能的保守性和分化情况。2.2.9qRT-PCR分析根据差异表达基因的序列信息，使用PrimerPremier5.0软件设计特异性引物，引物设计原则如下：引物长度为18-25bp，GC含量为40%-60%，Tm值为58-62℃，避免引物二聚体和发夹结构的形成。以小麦的看家基因β-actin作为内参基因，其引物序列为：F:5’-CCACACTGTGCCCATCTAC-3’，R:5’-GCCAGAGCAGTAATCTCCTT-3’。使用PrimeScriptRTreagentKitwithgDNAEraser（TaKaRa公司）将总RNA逆转录为cDNA，具体步骤按照试剂盒说明书进行。以cDNA为模板，使用SYBRPremixExTaqII（TaKaRa公司）进行qRT-PCR反应，反应体系为20μL，包括10μLSYBRPremixExTaqII、0.8μL上游引物（10μM）、0.8μL下游引物（10μM）、2μLcDNA模板和6.4μLddH2O。反应程序为：95℃预变性30s；95℃变性5s，60℃退火30s，共40个循环；最后进行熔解曲线分析，从65℃以0.5℃/s的速度升温至95℃，监测荧光信号的变化。每个样品设置3个技术重复，采用2-ΔΔCT法计算基因的相对表达量，其中ΔΔCT=(CT目的基因-CT内参基因)处理组-(CT目的基因-CT内参基因)对照组。将qRT-PCR检测结果与转录组测序得到的基因表达量数据进行相关性分析，验证转录组测序结果的准确性和可靠性。2.2.10基因序列检测选取部分差异表达基因，使用PCR扩增其全长编码区序列。根据基因序列设计扩增引物，引物两端分别添加合适的酶切位点，以便后续克隆和测序。PCR反应体系为25μL，包括12.5μL2×TaqPCRMasterMix、1μL上游引物（10μM）、1μL下游引物（10μM）、2μLcDNA模板和8.5μLddH2O。反应程序为：94℃预变性5min；94℃变性30s，55-60℃退火30s，72℃延伸1-2min（根据基因长度调整），共35个循环；72℃延伸10min。将PCR扩增产物进行1%琼脂糖凝胶电泳检测，切下目的条带，使用凝胶回收试剂盒（Qiagen公司）回收纯化DNA片段。将回收的DNA片段连接到pMD19-T载体（TaKaRa公司）上，转化大肠杆菌DH5α感受态细胞，在含有氨苄青霉素的LB平板上筛选阳性克隆。挑取阳性克隆进行菌落PCR鉴定，将鉴定正确的克隆送上海生工生物工程股份有限公司进行测序。使用DNAMAN软件对测序结果进行分析，与参考基因组序列进行比对，检测基因序列中的变异位点，包括单核苷酸多态性（SNP）、插入/缺失（InDel）等。分析这些变异位点对基因结构和功能的影响，如是否导致氨基酸序列改变、蛋白质结构变化等。2.2.11EST标记的开发利用转录组测序得到的小麦Unigene序列，使用MISA（MIcroSAtelliteidentificationtool）软件进行简单重复序列（SSR）搜索，搜索参数设置为：单核苷酸重复次数≥10，二核苷酸重复次数≥6，三核苷酸重复次数≥5，四核苷酸重复次数≥4，五核苷酸重复次数≥4，六核苷酸重复次数≥4。根据搜索到的SSR位点，使用Primer3软件设计引物，引物设计原则同2.2.9。引物设计完成后，使用BLAST工具将引物序列与小麦参考基因组进行比对，去除与基因组中其他区域具有高度同源性的引物，以确保引物的特异性。随机选取部分设计好的引物，以抗虫品种‘河农215’和感虫品种‘石新539’的基因组DNA为模板进行PCR扩增，PCR反应体系和程序同2.2.10。扩增产物经8%非变性聚丙烯酰胺凝胶电泳分离，银染显色，检测引物的多态性。对具有多态性的引物进行标记开发，命名为EST-SSR标记。2.2.12四引物扩增受阻体系PCR标记的开发针对在抗虫品种‘河农215’和感虫品种‘石新539’中检测到的SNP位点，使用四引物扩增受阻体系PCR（ARMS-PCR）技术开发分子标记。根据SNP位点两侧的序列，设计四条引物，其中两条为特异性引物，分别与SNP位点的两种等位基因互补配对，另外两条为通用引物。特异性引物的3’末端碱基与SNP位点的碱基互补，且在3’末端倒数第二个碱基处引入错配碱基，以增强引物的特异性。引物设计完成后，使用Oligo7.0软件对引物的特异性和Tm值进行评估，确保引物的特异性良好，且四条引物的Tm值相近。以抗虫品种‘河农215’和感虫品种‘石新539’的基因组DNA为模板，进行ARMS-PCR扩增。PCR反应体系为25μL，包括12.5μL2×TaqPCRMasterMix、1μL上游特异性引物（10μM）、1μL下游特异性引物（10μM）、0.5μL上游通用引物（10μM）、0.5μL下游通用引物（10μM）、2μLDNA模板和7.5μLddH2O。反应程序为：94℃预变性5min；94℃变性30s，58-62℃退火30s，72℃延伸30s，共35个循环；72℃延伸10min。扩增产物经2%琼脂糖凝胶电泳分离，观察条带的有无和大小，判断不同品种的基因型。对开发成功的ARMS-PCR标记进行验证，在多个小麦品种中进行扩增，检测标记的稳定性和通用性。2.2.13KASP标记的设计与开发基于转录组测序数据和SNP检测结果，选择与小麦抗麦红吸浆虫性状紧密相关的SNP位点，使用KASP（KompetitiveAllele-SpecificPCR）技术设计标记。利用LGC公司的KASPAssayDesignTool软件进行引物设计，软件会根据SNP位点的序列信息，自动设计两条等位基因特异性引物和一条通用引物。等位基因特异性引物的5’端分别标记有不同的荧光基团（FAM和HEX），通用引物的5’端标记有淬灭基团。引物设计完成后，由LGC公司合成引物。以抗虫品种‘河农215’和感虫品种‘石新539’的基因组DNA为模板，进行KASP反应。反应体系为5μL，包括2.5μLKASPMasterMix、0.07μLKASPAssayMix（包含三条引物）、2μLDNA模板和0.43μLddH2O。反应程序为：94℃预变性15min；94℃变性20s，61-55℃退火/延伸60s（每个循环退火温度降低0.6℃），共10个循环；94℃变性20s，55℃退火/延伸60s，共26个循环。反应结束后，使用荧光定量PCR仪（如Bio-RadCFX96）检测荧光信号，根据荧光信号的强度判断不同品种的基因型。对开发的KASP标记进行验证，在多个小麦品种和分离群体中进行检测，评估标记的准确性和有效性。三、结果与分析3.1供试材料表型分析对‘河农215’和‘石新539’两个小麦品种的抗虫性及相关农艺性状进行了调查与分析，结果如表3-1所示。在抗虫性方面，‘河农215’表现出较强的抗性，虫害损失率仅为8.5%，按照抗虫性分级标准，属于高抗品种；而‘石新539’的虫害损失率高达65.3%，表现为感虫品种。在农艺性状上，两个品种也存在明显差异。‘河农215’的株高为85.6cm，穗长10.2cm，每穗小穗数22.5个，千粒重45.3g；‘石新539’的株高为78.3cm，穗长8.8cm，每穗小穗数20.1个，千粒重38.6g。通过独立样本t检验分析，‘河农215’与‘石新539’在株高（P=0.012）、穗长（P=0.008）、每穗小穗数（P=0.021）和千粒重（P=0.005）上均存在显著差异（P<0.05）。品种虫害损失率（%）株高（cm）穗长（cm）每穗小穗数（个）千粒重（g）河农2158.5±1.285.6±2.510.2±0.822.5±1.545.3±1.8石新53965.3±3.878.3±1.88.8±0.620.1±1.238.6±1.5表3-1不同小麦品种抗虫性及农艺性状表型数据3.2羽化调查与qRT-PCR所用样品制备在小麦生长期间，对麦红吸浆虫的羽化情况进行了详细调查。结果显示，麦红吸浆虫成虫羽化呈现出明显的规律性，在4月下旬至5月上旬为羽化高峰期，这与小麦的抽穗期基本吻合，为后续研究小麦与麦红吸浆虫的互作关系提供了重要的时间节点信息。为进一步深入探究小麦抗麦红吸浆虫的分子机制，选取了抗虫品种‘河农215’和感虫品种‘石新539’在接虫后24h、48h和72h的叶片和穗部组织作为qRT-PCR所用样品。在每个时间点，对每个品种的每个组织均采集3个生物学重复，以确保实验结果的准确性和可靠性。采集后的样品迅速放入液氮中速冻，以最大限度地保持样品的生物学活性，随后保存于-80℃冰箱备用，为后续通过qRT-PCR技术验证转录组测序结果，筛选与抗虫性状相关的关键基因奠定了坚实的实验基础。3.322个小麦样本的转录组数据分析3.3.1转录组数据质量分析对22个小麦样本进行转录组测序后，首先运用FastQC软件对原始测序数据展开质量评估，各项评估指标结果表现出色。从碱基质量分布情况来看，在测序读长范围内，各碱基位置的质量值中位数均高于30，表明碱基测序质量极高，数据准确性有充分保障；测序错误率始终维持在较低水平，平均值小于0.01%，这使得后续基于测序数据的分析结果更为可靠；GC含量的分布较为稳定，平均值约为45%，符合小麦基因组的GC含量特征，说明测序数据未受异常因素干扰，数据质量稳定可靠。在去除低质量reads、接头序列和污染序列等杂质后，得到高质量的cleanreads。统计显示，cleanreads的比例均在95%以上，这意味着绝大部分测序数据可用于后续分析，有效数据量充足。以样本‘河农215_24h_leaf_1’为例，原始数据总量为50,000,000reads，经过质量控制后，得到47,500,000cleanreads，cleanreads比例高达95%，为后续的基因表达分析和差异基因筛选奠定了坚实基础。各样本的测序深度均达到了100X以上，能够充分覆盖小麦基因组，确保了基因表达量的准确测定和低丰度表达基因的有效检测，为深入挖掘小麦抗麦红吸浆虫相关基因提供了充足的数据支持。3.3.2利用差异表达分析挖掘抗虫DEGs运用DESeq2软件，以|log2(FoldChange)|≥1且FDR≤0.05作为筛选标准，对22个小麦样本进行差异表达分析，共筛选出5,680个差异表达基因（DEGs）。在这些DEGs中，与感虫品种‘石新539’相比，抗虫品种‘河农215’中表达上调的基因有3,210个，表达下调的基因有2,470个。对这些差异表达基因的表达模式进行分析发现，在接虫后24h，抗虫品种‘河农215’中就有1,200个基因表达发生显著变化，其中上调基因750个，下调基因450个。这些早期响应基因可能在小麦感知麦红吸浆虫侵害信号并启动防御反应的过程中发挥关键作用。到接虫后48h，差异表达基因的数量进一步增加到2,500个，上调基因1,500个，下调基因1,000个，表明小麦的防御反应在持续增强，更多的基因参与到抗虫过程中。在接虫后72h，差异表达基因的数量达到3,800个，上调基因2,200个，下调基因1,600个，此时小麦的抗虫防御反应达到较为强烈的状态，涉及到更多的生物学过程和代谢途径的调控。通过对不同时间点差异表达基因的聚类分析，发现了多个具有相似表达模式的基因簇。其中一个基因簇在接虫后24h到72h期间，表达量持续上调，进一步分析这些基因的功能注释信息，发现它们主要参与了植物激素信号转导、次生代谢产物合成等与抗虫密切相关的生物学过程。这些基因在抗虫品种‘河农215’中的持续上调表达，可能是其抵御麦红吸浆虫侵害的重要分子机制之一。3.3.3BSR法定位区间和挖掘抗虫相关基因通过对‘河农215’和‘石新539’的F2代群体构建的抗虫池和感虫池进行BSR-Seq分析，成功定位到多个与小麦抗麦红吸浆虫性状紧密关联的SNP位点和候选基因区域。在2B染色体上，确定了一个长度约为5Mb的候选基因区域，该区域内包含120个基因。通过对这些基因的功能注释和表达分析，筛选出5个可能与抗虫性状相关的关键基因。基因TraesCS2B01G123400编码一种蛋白激酶，该蛋白激酶在植物信号转导过程中发挥着重要作用，可能参与小麦对麦红吸浆虫侵害信号的感知和传递，进而激活下游的抗虫防御反应。TraesCS2B01G123500基因编码的蛋白质含有一个与防御反应相关的结构域，推测其在小麦抗虫防御过程中直接参与了对吸浆虫的防御作用。TraesCS2B01G123600基因则参与了植物细胞壁的合成和修饰过程，在受到麦红吸浆虫侵害时，其表达量发生显著变化，可能通过增强小麦细胞壁的强度，阻止吸浆虫幼虫的侵入。在6B染色体上，也定位到一个关键的抗虫基因区域，长度约为3Mb，包含80个基因。经过深入分析，筛选出3个抗虫相关的关键基因。TraesCS6B01G234500基因编码的转录因子在抗虫池中表达量显著高于感虫池，该转录因子可能通过调控其他抗虫相关基因的表达，参与小麦抗麦红吸浆虫的分子调控网络。TraesCS6B01G234600基因参与了植物激素茉莉酸的合成途径，茉莉酸在植物抗虫反应中起着关键的调控作用，该基因的差异表达可能影响茉莉酸的合成和信号转导，从而影响小麦的抗虫能力。TraesCS6B01G234700基因编码的蛋白质具有抗氧化酶活性，在吸浆虫侵害过程中，可能通过清除细胞内的活性氧，保护小麦细胞免受氧化损伤，增强小麦的抗虫性。3.3.4利用基因功能注释挖掘抗虫相关基因将筛选出的5,680个差异表达基因提交至GO数据库进行功能注释，从生物过程、细胞组分和分子功能三个层面进行分析。在生物过程方面，发现大量差异表达基因富集在“防御反应”“信号转导”“次生代谢产物合成”等生物学过程。其中，参与“防御反应”的基因有850个，这些基因在小麦受到麦红吸浆虫侵害时，通过激活一系列防御机制，如合成抗菌肽、植保素等物质，抵御吸浆虫的侵害；参与“信号转导”的基因有620个，它们在小麦感知吸浆虫侵害信号并将其传递到细胞内的过程中发挥关键作用，通过激活下游的防御基因表达，启动抗虫反应；参与“次生代谢产物合成”的基因有580个，这些基因参与合成多种次生代谢产物，如酚类、萜类等，这些次生代谢产物具有抗菌、抗虫等生物活性，是小麦抗虫防御的重要物质基础。在细胞组分层面，差异表达基因主要富集在“细胞膜”“细胞壁”“细胞核”等细胞结构中。在细胞膜上，一些基因编码的离子通道蛋白和受体蛋白，可能参与小麦对吸浆虫侵害信号的感知和传递；在细胞壁中，相关基因参与细胞壁的合成和修饰，增强细胞壁的强度，阻止吸浆虫幼虫的侵入；在细胞核中，许多转录因子基因的表达发生变化，这些转录因子通过调控其他抗虫相关基因的表达，参与小麦抗虫的分子调控网络。在分子功能方面，差异表达基因主要富集在“酶活性”“转录因子活性”“受体活性”等分子功能类别。具有“酶活性”的基因参与了多种代谢途径和防御反应，如过氧化物酶、多酚氧化酶等，它们在清除活性氧、合成次生代谢产物等过程中发挥重要作用；具有“转录因子活性”的基因通过结合到靶基因的启动子区域，调控基因的转录，从而影响小麦抗虫相关基因的表达；具有“受体活性”的基因编码的受体蛋白，能够识别吸浆虫侵害信号，启动小麦的抗虫防御反应。3.3.5利用比较基因组学发掘抗虫相关基因将筛选出的差异表达基因序列提交到NCBI的BLAST工具中，与nr数据库进行比对，选取同源性较高（Identity≥80%，E-value≤1e-10）的基因作为同源基因，共找到3,200个同源基因，涉及15个物种。其中，与水稻的同源基因有560个，与玉米的同源基因有48