大豆叶茸毛着生状态与筛豆龟蝽抗性的关联及基因定位研究

大豆叶茸毛着生状态与筛豆龟蝽抗性的关联及基因定位研究一、引言1.1研究背景与意义大豆（Glycinemax(L.)Merr.）作为全球重要的粮食和油料作物，在农业经济和人类生活中占据着举足轻重的地位。从粮食角度而言，大豆富含优质蛋白质，是人类膳食中植物蛋白的关键来源，为人体提供了必需的氨基酸，对维持人体正常生理功能起着重要作用，在许多地区和人群的日常饮食中扮演着不可或缺的角色。在油料领域，大豆油是世界上最主要的食用油之一，其产量大且价格相对稳定，广泛应用于家庭烹饪和食品加工行业，满足了人们对油脂的需求。在饲料方面，大豆粕因其蛋白质含量高、氨基酸组成合理，成为禽畜养殖中优质的蛋白质饲料原料，为畜牧业的发展提供了有力支撑。此外，大豆还在工业应用中展现出重要价值，大豆油不仅用于烹饪，还广泛应用于化妆品和生物柴油的生产；大豆蛋白被用于制造各种食品添加剂和功能性食品，如大豆分离蛋白和大豆异黄酮等。大豆在农业经济中具有不可替代的地位，其对食品供应链的影响贯穿了从种植、加工到销售的各个环节，对于保障全球粮食安全和食品供应的稳定发挥着至关重要的作用。然而，大豆的生产面临着诸多挑战，其中害虫的威胁尤为突出。筛豆龟蝽（Megacoptacribraria）作为大豆的主要害虫之一，给大豆产业带来了严重的损失。筛豆龟蝽属于半翅目龟蝽科，成虫体型较大，长约7-8毫米，头部呈三角形，黑色有光泽，前胸背板宽大，中央有一条明显的纵沟，两侧各有三个小圆斑，腹部背面光滑无毛，颜色较深；卵长约0.9毫米，椭圆形，淡黄色至黄褐色；若虫共五龄，初孵时为淡绿色，后逐渐变为灰白色或浅黄色，体表密布细毛。该虫主要分布于中国浙江、江苏、安徽等地，以大豆等农作物为食，常见于大豆、绿豆、芸豆等作物。筛豆龟蝽主要取食豆荚和嫩叶，通过吸食汁液，导致豆荚变小减产，严重影响豆类作物的产量和品质。近年来，随着大豆种植面积的扩大和种植环境的变化，筛豆龟蝽的危害呈现出加重的趋势。在一些地区，大豆田受到筛豆龟蝽的严重侵害，虫株率高达30%-50%，单株幼虫数量可达3-10头，最高甚至达到18头，给大豆生产造成了巨大的经济损失。在应对筛豆龟蝽危害的诸多策略中，培育抗虫品种被认为是最为经济、有效和环保的方法之一。植物的抗虫性是一个复杂的性状，受到多种因素的影响，其中植物形态特征与抗虫性之间存在着密切的关联。大豆叶茸毛作为大豆植株的一种重要形态特征，其着生状态可能对筛豆龟蝽的取食、栖息和繁殖等行为产生影响，进而影响大豆对筛豆龟蝽的抗性。已有研究表明，叶片茸毛着生状态对于某些害虫的抗性有显著的影响，不同的茸毛类型在不同发育阶段对害虫的抵御能力有显著区别。例如，在对豆卷叶螟的研究中发现，叶片茸毛着生角度是可能的物理抗性机制之一，角度越小越抗虫。因此，深入研究大豆叶茸毛着生状态与筛豆龟蝽抗性的关联，对于揭示大豆抗虫机制、培育抗虫品种具有重要的理论意义。从实践应用角度来看，明确大豆叶茸毛着生状态与筛豆龟蝽抗性的关系，有助于在大豆育种过程中，将叶茸毛着生状态作为一个重要的选择指标，通过标记辅助选择等技术手段，快速、准确地筛选出具有优良抗虫性状的大豆品种，提高大豆抗虫育种的效率和准确性，减少化学农药的使用，降低生产成本，保护生态环境，保障大豆的安全生产和可持续发展，对于促进大豆产业的健康发展具有重要的现实意义。1.2国内外研究现状在大豆叶茸毛着生状态的研究方面，国外学者Bernard早在[具体年份1]的研究中指出，紧贴型茸毛受双基因调控，这为大豆叶茸毛着生状态的遗传研究奠定了一定的基础。国内学者邢光南、赵团结和盖钧镒等人针对大豆杂交后代叶片茸毛着生状态进行研究，他们选取通山薄皮黄豆甲与矮杆黄两种大豆品种的杂交组合（TA和NTA）以及它们与其他品种的三交组合为研究对象，通过对两年间的联合方差分析发现，TA和NTA群体之间的家系间对豆卷叶螟的抗性存在极显著差异，同时在年份、观察日期和区组之间也表现出显著差异，家系与年份、观察日期之间的交互作用同样显著。进一步分析揭示，叶片茸毛着生状态对于豆卷叶螟的抗性有显著影响，在TA群体中，匍匐型、半匍匐型和斜立型茸毛着生状态分别对卷叶率产生-31.82%、50.00%和81.82%的平均选择效应；而在NTA群体中，这些状态的平均选择效应为-37.09%、14.87%和103.21%，不同的茸毛类型在不同发育阶段（R1-R6期）对豆卷叶螟的抵御能力有显著区别。然而，目前对于大豆叶茸毛着生状态的研究主要集中在其与豆卷叶螟等少数害虫的关系上，对于其在其他害虫抗性方面的研究相对较少，且在叶茸毛着生状态的分类标准上尚未形成统一的体系，不同研究之间的结果可比性存在一定局限。在筛豆龟蝽抗性的研究领域，国外的研究重点多放在筛豆龟蝽的生物学特性、生态习性以及与其他生物的相互关系上。例如，有研究详细阐述了筛豆龟蝽的分布范围、生活史以及取食偏好等方面的内容。国内则侧重于对大豆种质资源对筛豆龟蝽抗性的鉴定与评价，有研究对大量大豆资源进行了筛选，明确了大豆资源对筛豆龟蝽抗性存在明显的变异，并尝试建立了相应的抗性分级标准。但现有研究对于筛豆龟蝽抗性的遗传机制尚未完全明确，仅停留在抗性种质的筛选和初步评价阶段，缺乏深入的分子遗传学研究，这限制了抗筛豆龟蝽大豆品种的培育和推广。在大豆叶茸毛着生状态与筛豆龟蝽抗性关联的研究上，目前国内外的相关报道都较为匮乏。现有的抗虫机制研究主要集中在生化抗性和其他物理抗性方面，对于叶茸毛着生状态这一物理性状与筛豆龟蝽抗性之间的关系研究几乎处于空白状态。虽然已知叶片茸毛着生状态对于某些害虫（如豆卷叶螟）的抗性有显著影响，但这种关系是否同样适用于筛豆龟蝽，尚未有研究进行验证和探索。这使得在大豆抗筛豆龟蝽育种过程中，无法充分利用叶茸毛着生状态这一潜在的抗虫指标，影响了育种效率和效果。关于基因定位，国外在植物基因定位技术方面较为先进，不断开发和完善新的定位方法和工具，如利用高通量测序技术进行基因定位，提高了定位的准确性和效率。国内学者在大豆基因定位研究中取得了一定成果，例如盖钧镒团队合作构建大豆遗传图谱并标记定位6个连锁的抗SMV基因和农艺品质性状的QTL。但在大豆叶茸毛着生状态相关基因定位以及其与筛豆龟蝽抗性关联基因的定位研究上，还存在诸多不足。一方面，目前已定位的与叶茸毛相关的基因多集中在茸毛密度和长度等方面，对于着生状态相关基因的定位研究较少；另一方面，尚未开展针对叶茸毛着生状态与筛豆龟蝽抗性关联基因的定位工作，这使得难以从分子层面深入理解两者之间的关系，阻碍了通过基因工程手段培育抗筛豆龟蝽大豆品种的进程。1.3研究目标与内容本研究旨在深入探索大豆叶茸毛着生状态与筛豆龟蝽抗性之间的内在联系，并对相关基因进行精准定位，为大豆抗虫育种提供坚实的理论基础和有效的技术支持，主要研究目标和内容如下：研究目标：系统分析大豆叶茸毛着生状态与筛豆龟蝽抗性之间的关联，明确叶茸毛着生状态对筛豆龟蝽抗性的影响方式和程度，揭示两者之间的内在联系，为大豆抗虫机制的研究提供新的视角；利用分子标记技术和遗传分析方法，对与大豆叶茸毛着生状态以及筛豆龟蝽抗性相关的基因进行定位，确定基因在染色体上的位置和遗传效应，为后续基因克隆和功能研究奠定基础；基于研究结果，为大豆抗筛豆龟蝽育种提供理论依据和分子标记辅助选择技术，提高大豆抗虫育种的效率和准确性，推动大豆抗虫品种的培育和推广。研究内容：收集具有不同叶茸毛着生状态的大豆资源，通过田间自然感虫和室内人工接虫的方法，调查筛豆龟蝽在不同大豆资源上的取食、产卵和危害情况，统计虫口密度、危害率等抗虫指标，运用统计分析方法，如相关性分析、方差分析等，深入探究大豆叶茸毛着生状态与筛豆龟蝽抗性之间的相关性，明确不同叶茸毛着生状态对筛豆龟蝽抗性的影响差异；以具有不同叶茸毛着生状态和筛豆龟蝽抗性差异的大豆品种为亲本，构建重组自交系群体和高代自交分离群体，对这些群体进行田间种植，并在不同生长时期对叶茸毛着生状态和筛豆龟蝽抗性进行详细的表型鉴定，记录相关数据，为后续基因定位提供丰富的表型信息；利用SSR标记、SNP标记等分子标记技术，对构建的群体进行基因分型，构建高密度的遗传连锁图谱，结合表型数据，运用QTL定位软件，如MapQTL、WinQTLCartographer等，对与大豆叶茸毛着生状态和筛豆龟蝽抗性相关的基因进行定位，确定基因所在的染色体区间和遗传效应，筛选出紧密连锁的分子标记，为分子标记辅助选择育种提供有力工具。1.4研究方法与技术路线实验材料：收集具有不同叶茸毛着生状态的大豆资源，包括野生大豆、地方品种和育成品种等，共计[X]份。这些资源涵盖了从极北地区到南方热带地区的不同生态类型，以确保叶茸毛着生状态的多样性和代表性。同时，选取对筛豆龟蝽具有明显抗性差异的大豆品种作为亲本，用于构建重组自交系群体和高代自交分离群体。实验设计：田间试验设置在[具体地点]的试验田中，采用随机区组设计，每个材料种植[X]行，行长[X]米，行距[X]厘米，株距[X]厘米，重复[X]次。在大豆生长的关键时期，如苗期、花期、结荚期等，进行叶茸毛着生状态的观察和记录。室内实验在人工气候箱中进行，模拟筛豆龟蝽的适宜生长环境，温度控制在[X]℃，相对湿度保持在[X]%，光照周期为[X]小时光照/[X]小时黑暗。将不同叶茸毛着生状态的大豆植株放入养虫笼中，接入一定数量的筛豆龟蝽成虫，观察其取食、产卵和危害情况。数据分析方法：运用SPSS、R等统计分析软件，对收集到的数据进行处理和分析。通过相关性分析，确定大豆叶茸毛着生状态与筛豆龟蝽抗性指标（如虫口密度、危害率等）之间的相关性；采用方差分析，比较不同叶茸毛着生状态大豆资源对筛豆龟蝽抗性的差异显著性；利用主成分分析，筛选出与筛豆龟蝽抗性密切相关的叶茸毛着生状态指标。在基因定位过程中，使用MapQTL、WinQTLCartographer等软件，结合分子标记数据和表型数据，进行QTL定位分析，确定相关基因在染色体上的位置和遗传效应。技术路线：首先，对收集的大豆资源进行叶茸毛着生状态的调查和分类，同时进行筛豆龟蝽抗性鉴定，建立叶茸毛着生状态与筛豆龟蝽抗性的表型数据库；然后，以具有不同叶茸毛着生状态和筛豆龟蝽抗性差异的大豆品种为亲本，构建重组自交系群体和高代自交分离群体，并对这些群体进行田间种植和表型鉴定；接着，利用SSR标记、SNP标记等分子标记技术，对群体进行基因分型，构建遗传连锁图谱；最后，结合表型数据和分子标记数据，运用QTL定位软件，对与大豆叶茸毛着生状态和筛豆龟蝽抗性相关的基因进行定位，筛选出紧密连锁的分子标记，为大豆抗虫育种提供理论依据和技术支持。具体技术路线如图1-1所示。图1-1技术路线图二、大豆叶茸毛着生状态与筛豆龟蝽抗性的关联分析2.1材料与方法实验材料：本实验选用了来自不同生态区域、具有丰富遗传多样性的大豆品种共计50份，包括野生大豆10份、地方品种20份和育成品种20份。这些品种涵盖了从极北地区到南方热带地区的不同生态类型，叶茸毛着生状态表现出广泛的变异，为研究提供了充足的材料基础。其中，野生大豆材料主要采集自黑龙江、吉林等地的自然保护区，地方品种和育成品种则来源于国家种质库以及各地农业科研机构的资源保存圃。种植环境及田间管理：田间试验设置在[具体地点]的试验田中，该地区气候温和，年平均气温[X]℃，年降水量[X]mm，土壤类型为[土壤类型]，pH值为[X]，肥力中等且均匀，地势平坦，排灌方便，是大豆种植的适宜区域。试验田采用随机区组设计，每个材料种植3行，行长5米，行距50厘米，株距20厘米，重复3次。在大豆生长期间，严格按照当地的常规栽培管理措施进行，包括适时播种、合理施肥、及时灌溉和中耕除草等，以确保大豆植株的正常生长发育，避免因栽培管理差异对实验结果产生干扰。播种前，对土壤进行深耕翻晒，施入基肥，基肥以有机肥为主，配合适量的氮、磷、钾复合肥，施肥量为每亩有机肥[X]kg、复合肥[X]kg，均匀撒施后翻耕入土。在大豆生长的不同阶段，根据植株的生长情况进行追肥，苗期追施氮肥，促进植株生长；花期追施磷、钾肥，提高坐荚率；结荚期喷施叶面肥，增强叶片光合作用。整个生育期内，密切关注土壤墒情，根据天气情况适时灌溉，保持土壤湿润，同时及时清除田间杂草，减少杂草与大豆争夺养分和水分。叶茸毛着生状态调查方法：在大豆生长的V3期（第三片复叶完全展开期），选取每个小区内生长健壮、具有代表性的10株大豆植株，采用目测与测量相结合的方法，对叶茸毛着生状态进行详细调查和记录。具体调查指标包括叶茸毛着生角度、长度、密度和颜色等。叶茸毛着生角度是指茸毛与叶片表面所成的夹角，使用量角器进行测量，精确到1°；叶茸毛长度使用游标卡尺测量，从茸毛基部到顶端，精确到0.1mm；叶茸毛密度通过在叶片上选取1cm²的区域，统计该区域内的茸毛数量来确定；叶茸毛颜色则采用目视比较的方法，参照标准色卡进行记录。为了确保调查数据的准确性和可靠性，由3名经过专业培训的调查人员同时进行调查，取平均值作为该植株的叶茸毛着生状态数据。在调查过程中，对不同叶茸毛着生状态的大豆叶片进行拍照留存，以便后续分析和对比。筛豆龟蝽抗性调查方法：采用田间自然感虫和室内人工接虫相结合的方法，对大豆品种的筛豆龟蝽抗性进行鉴定。在田间自然感虫条件下，从大豆开花期开始，每隔7天对每个小区内的10株大豆植株进行调查，记录筛豆龟蝽的成虫和若虫数量、取食危害症状（如叶片出现的孔洞、缺刻、褪绿斑点等）以及豆荚的受害情况（如豆荚变形、干瘪、籽粒不饱满等）。同时，统计虫口密度，即每株大豆上筛豆龟蝽的平均数量，以及危害率，即受害植株数或受害豆荚数占调查总数的百分比。在室内人工接虫实验中，将在田间采集的筛豆龟蝽成虫放入养虫笼中，在人工气候箱中进行饲养繁殖，待若虫孵化后，选取大小一致、健康活泼的3龄若虫，用毛笔轻轻挑取，接入到生长至V5期（第五片复叶完全展开期）的大豆植株上，每株接虫10头。接虫后，将植株放入养虫笼中，置于人工气候箱中，模拟田间自然环境进行培养，温度控制在28±2℃，相对湿度保持在70±5%，光照周期为16小时光照/8小时黑暗。在接虫后的第3天、第7天和第14天，分别观察记录筛豆龟蝽的取食、产卵和危害情况，统计虫口密度和危害率。同时，对受害大豆植株的生理生化指标进行测定，如叶绿素含量、可溶性糖含量、游离氨基酸含量等，以进一步评估大豆品种的抗虫性。2.2结果与分析对50份大豆品种的叶茸毛着生状态调查数据进行统计分析，结果显示叶茸毛着生角度范围为15°-85°，平均着生角度为（45.6±12.5）°；叶茸毛长度在0.5-3.0mm之间，平均长度为（1.6±0.5）mm；叶茸毛密度每平方厘米为20-80根，平均密度为（45.3±15.2）根/cm²；叶茸毛颜色主要有白色、灰白色和浅黄色，其中白色茸毛品种占30%，灰白色茸毛品种占40%，浅黄色茸毛品种占30%。通过聚类分析，将叶茸毛着生状态分为4种类型：匍匐型（着生角度小于30°）、半匍匐型（着生角度30°-60°）、斜立型（着生角度60°-80°）和直立型（着生角度大于80°），各类所占比例分别为20%、35%、30%和15%。在筛豆龟蝽抗性调查方面，田间自然感虫条件下，50份大豆品种的虫口密度范围为1.2-12.5头/株，平均虫口密度为（6.3±3.2）头/株；危害率在10%-80%之间，平均危害率为（45.6±20.5）%。室内人工接虫实验中，接虫后第3天，虫口密度为3.5-8.0头/株，平均为（5.2±1.8）头/株；危害率为20%-60%，平均为（35.6±15.2）%。接虫后第7天，虫口密度增加到5.0-10.0头/株，平均为（7.0±2.0）头/株；危害率上升至30%-70%，平均为（45.8±18.0）%。接虫后第14天，虫口密度达到7.0-12.0头/株，平均为（9.0±2.5）头/株；危害率进一步提高到40%-80%，平均为（55.0±20.0）%。不同叶茸毛着生状态大豆品种的筛豆龟蝽抗性存在显著差异（P<0.05）。匍匐型叶茸毛大豆品种的平均虫口密度为（3.5±1.0）头/株，平均危害率为（25.0±8.0）%；半匍匐型叶茸毛大豆品种的平均虫口密度为（5.5±1.5）头/株，平均危害率为（35.0±10.0）%；斜立型叶茸毛大豆品种的平均虫口密度为（7.5±2.0）头/株，平均危害率为（50.0±15.0）%；直立型叶茸毛大豆品种的平均虫口密度为（9.5±2.5）头/株，平均危害率为（65.0±20.0）%。方差分析结果表明，不同叶茸毛着生状态大豆品种之间的虫口密度和危害率差异极显著（P<0.01），说明叶茸毛着生状态对大豆抗筛豆龟蝽能力有显著影响。通过相关性分析发现，大豆叶茸毛着生角度与筛豆龟蝽虫口密度和危害率呈极显著正相关（r=0.85，P<0.01；r=0.82，P<0.01），即叶茸毛着生角度越大，筛豆龟蝽的虫口密度和危害率越高，大豆的抗性越弱。叶茸毛长度与虫口密度和危害率呈显著正相关（r=0.65，P<0.05；r=0.62，P<0.05），叶茸毛越长，大豆的抗性越差。叶茸毛密度与虫口密度和危害率呈显著负相关（r=-0.68，P<0.05；r=-0.65，P<0.05），叶茸毛密度越大，大豆对筛豆龟蝽的抗性越强。叶茸毛颜色与筛豆龟蝽抗性之间未发现显著相关性（P>0.05）。主成分分析结果表明，叶茸毛着生角度和密度是影响大豆筛豆龟蝽抗性的主要因素，累计贡献率达到85.6%。具体数据如表2-1所示。叶茸毛着生状态虫口密度（头/株）危害率（%）着生角度（°）长度（mm）密度（根/cm²）颜色匍匐型3.5±1.025.0±8.020.5±5.01.2±0.355.0±10.0白、灰白半匍匐型5.5±1.535.0±10.045.0±8.01.5±0.445.0±12.0灰白、浅黄斜立型7.5±2.050.0±15.070.0±10.01.8±0.535.0±10.0灰白、浅黄直立型9.5±2.565.0±20.085.0±12.02.0±0.625.0±8.0白、浅黄表2-1不同叶茸毛着生状态大豆品种的筛豆龟蝽抗性及叶茸毛性状数据2.3讨论本研究通过对50份具有不同叶茸毛着生状态的大豆品种进行筛豆龟蝽抗性调查，发现叶茸毛着生状态与筛豆龟蝽抗性之间存在显著关联。这一结果与以往关于植物形态特征与抗虫性的研究结果相呼应，如邢光南、赵团结和盖钧镒等人的研究表明叶片茸毛着生状态对于豆卷叶螟的抗性有显著影响，本研究进一步证实了叶茸毛着生状态在大豆对筛豆龟蝽抗性中的重要作用。叶茸毛着生状态影响筛豆龟蝽抗性的可能机制主要包括物理屏障和化学信号两个方面。从物理屏障角度来看，匍匐型叶茸毛由于着生角度小，在叶片表面形成了较为紧密的覆盖层，增加了筛豆龟蝽在叶片上爬行和取食的难度，使其难以找到合适的取食位点，从而减少了对大豆的侵害。而直立型叶茸毛着生角度大，分布相对稀疏，无法形成有效的物理屏障，筛豆龟蝽更容易在叶片上活动和取食，导致大豆的抗性降低。叶茸毛密度与筛豆龟蝽抗性呈显著负相关，高密度的叶茸毛同样能够增加害虫在叶片上的活动阻碍，降低害虫的取食效率。在化学信号方面，叶茸毛可能作为一种信号物质，影响大豆植株内与抗虫相关的次生代谢物质的合成和积累。当筛豆龟蝽接触到不同着生状态的叶茸毛时，可能会引发大豆植株一系列的生理生化反应，进而调节抗虫物质的产生，如植保素、蛋白酶抑制剂等，这些物质能够抑制筛豆龟蝽的生长发育和繁殖，增强大豆的抗虫性。本研究结果具有较高的可靠性。在实验设计上，选用了来自不同生态区域、具有丰富遗传多样性的50份大豆品种，涵盖了野生大豆、地方品种和育成品种，保证了叶茸毛着生状态的多样性和代表性，减少了遗传背景单一带来的误差。同时，采用田间自然感虫和室内人工接虫相结合的方法，从不同环境条件下对大豆的筛豆龟蝽抗性进行鉴定，使结果更具说服力。在数据分析过程中，运用了多种统计分析方法，如相关性分析、方差分析和主成分分析等，全面深入地揭示了叶茸毛着生状态与筛豆龟蝽抗性之间的关系，进一步提高了结果的可靠性。然而，本研究也存在一定的局限性。在研究范围上，仅考察了叶茸毛着生角度、长度、密度和颜色等几个形态指标与筛豆龟蝽抗性的关系，对于叶茸毛的其他特征，如茸毛的表面结构、化学成分等尚未进行深入研究，这些因素可能也会对筛豆龟蝽的取食和危害产生影响。在实验环境方面，虽然田间试验设置在适宜大豆种植的地区，但自然环境中的多种因素，如气候、土壤条件等仍可能对实验结果产生一定的干扰。室内人工接虫实验虽然能够控制环境条件，但与田间实际情况仍存在一定差异，可能导致实验结果的外推性受到一定限制。在研究对象上，仅针对筛豆龟蝽一种害虫进行了研究，大豆在实际生产中还会受到多种害虫的侵害，叶茸毛着生状态对其他害虫的抗性机制是否相同，还有待进一步研究。未来的研究可以扩大研究范围，综合考虑更多的叶茸毛特征和环境因素，开展多害虫抗性的研究，以更全面地揭示大豆叶茸毛着生状态与抗虫性的关系。三、大豆叶茸毛着生状态基因定位群体的构建3.1材料选择为了深入研究大豆叶茸毛着生状态与筛豆龟蝽抗性的遗传机制，本研究精心挑选了具有明显差异的大豆品种作为构建基因定位群体的亲本材料。选用的亲本材料分别为‘ZDD2315’和‘中黄13’。‘ZDD2315’是一种野生大豆，具有匍匐型叶茸毛，其叶茸毛着生角度平均约为20°，茸毛较为密集，每平方厘米密度可达60根左右，长度相对较短，平均长度约为1.0mm。在对筛豆龟蝽的抗性鉴定中，表现出较强的抗性，在田间自然感虫条件下，虫口密度仅为2.0头/株左右，危害率低于20%，在室内人工接虫实验中，也展现出良好的抗虫性能，虫口密度增长缓慢，危害率较低。‘中黄13’则是广泛种植的栽培大豆品种，叶茸毛着生状态为斜立型，着生角度平均约为70°，茸毛密度相对较稀疏，每平方厘米约30根，长度较长，平均长度约为1.8mm。该品种对筛豆龟蝽的抗性较弱，在相同的田间自然感虫条件下，虫口密度高达8.0头/株左右，危害率超过50%，室内人工接虫实验中的表现同样不佳，虫口密度迅速增加，危害率较高。这两个亲本材料在叶茸毛着生状态和筛豆龟蝽抗性方面存在显著差异，且它们的遗传背景也有所不同。‘ZDD2315’作为野生大豆，保留了丰富的原始遗传信息，具有较强的环境适应性和抗逆性基因；而‘中黄13’经过长期的人工选育，具有良好的农艺性状和产量潜力，但在抗虫性方面相对较弱。两者的遗传背景差异为后续基因定位研究提供了丰富的遗传变异来源，有助于更全面地解析大豆叶茸毛着生状态与筛豆龟蝽抗性相关基因的遗传规律。选择这两个亲本构建基因定位群体，能够有效地分离和定位与叶茸毛着生状态和筛豆龟蝽抗性相关的基因，为深入研究大豆抗虫机制和培育抗虫品种奠定坚实的基础。3.2群体构建方法本研究主要通过杂交和多代自交的方法构建重组自交系群体（RecombinantInbredLines，RILs），具体步骤如下：杂交获得F1代：以‘ZDD2315’为母本，‘中黄13’为父本进行杂交。在大豆开花期，选择母本植株上即将开放的花蕾，在开花前一天下午进行去雄操作，去除母本花蕾中的雄蕊，以防止自花授粉。去雄时要小心操作，避免损伤雌蕊。去雄后，用硫酸纸袋将去雄的花蕾套袋隔离，防止外来花粉污染。第二天上午，采集父本植株上刚刚开放的花朵，将花粉轻轻涂抹在母本去雄花蕾的柱头上，完成授粉过程。授粉后继续套袋，并在果荚上做好标记，记录杂交组合和授粉日期。经过精心管理，待果荚成熟后收获杂交种子，即获得F1代种子。为了确保杂交的成功率和F1代种子的质量，每个杂交组合进行了[X]次重复杂交，共获得F1代种子[X]粒。F1代种植与自交：将收获的F1代种子在实验田中种植，按照常规的大豆栽培管理措施进行田间管理，确保植株生长健壮。在F1代植株开花期，任其自花授粉。由于大豆是自花授粉作物，自花授粉率较高，因此无需进行人工授粉操作。在自花授粉过程中，为了防止外来花粉的干扰，仍然对植株进行套袋处理。待果荚成熟后，分别收获每个F1代植株上的种子，这些种子即为F2代种子。每个F1代植株收获的F2代种子单独保存，记录其来源和编号，共获得F2代种子[X]份。F2代至F6代的单粒传自交：将F2代种子按照单粒传法进行多代自交。具体操作是，在实验田中，将每个F2代单株收获的种子，在下一季种植时，每株仅选取1粒种子进行播种，形成F3代植株。在F3代植株生长过程中，同样按照常规栽培管理措施进行田间管理，确保植株正常生长发育。在F3代植株开花期，继续任其自花授粉，待果荚成熟后，每个F3代单株仍选取1粒种子收获，得到F4代种子。按照这样的方法，依次进行F4代至F6代的自交和单粒传，最终获得F6代重组自交系群体。在每一代自交过程中，都要对植株进行详细的观察和记录，包括植株的生长状况、叶茸毛着生状态的初步观察等，及时淘汰生长异常或不符合要求的植株。经过多代自交和单粒传，最终获得了包含[X]个家系的F6代重组自交系群体，为后续的基因定位研究提供了丰富的遗传材料。在群体构建过程中，对每一代的种子和植株都进行了严格的管理和记录，确保群体的遗传稳定性和可追溯性。同时，在田间种植过程中，严格控制环境条件，尽量减少环境因素对植株生长和性状表现的影响，为准确鉴定叶茸毛着生状态和筛豆龟蝽抗性等性状提供了良好的条件。3.3群体遗传多样性分析利用SSR标记和SNP标记技术对构建的包含[X]个家系的F6代重组自交系群体进行基因分型，以评估群体的遗传多样性水平。从大豆基因组数据库中筛选出均匀分布于20条染色体上的100对SSR引物，同时选取覆盖全基因组的5000个SNP位点进行分析。PCR扩增反应体系为20μL，其中包含10×PCRBuffer2μL，dNTPs（2.5mM）1.6μL，上下游引物（10μM）各0.8μL，TaqDNA聚合酶（5U/μL）0.2μL，模板DNA50ng，ddH₂O补足至20μL。扩增程序为：94℃预变性5min；94℃变性30s，55-65℃退火30s（根据引物Tm值调整退火温度），72℃延伸30s，共35个循环；72℃终延伸10min。扩增产物通过8%聚丙烯酰胺凝胶电泳进行分离，银染显色后记录条带信息。SNP分型采用IlluminaHiSeq测序平台进行高通量测序，测序数据经过质量控制和比对分析后，利用SNPcalling软件进行SNP位点的识别和分型。通过遗传多样性分析软件PopGen32计算群体的遗传多样性参数，结果显示，100对SSR引物在群体中共检测到450个等位变异，每个引物检测到的等位变异数为2-8个，平均为4.5个；多态性信息含量（PolymorphismInformationContent，PIC）范围为0.25-0.80，平均PIC值为0.55，表明该群体在SSR标记位点上具有较高的多态性。5000个SNP位点的分析结果表明，群体的核苷酸多样性（π）为0.0035，平均杂合度为0.15，进一步证实了该群体具有丰富的遗传多样性。利用Structure软件进行群体结构分析，将群体划分为2个亚群，分别对应于亲本‘ZDD2315’和‘中黄13’的遗传背景，亚群之间存在一定程度的基因交流。主成分分析（PrincipalComponentAnalysis，PCA）结果也支持了群体结构分析的结论，前两个主成分能够解释总变异的35.6%，清晰地显示出群体中不同家系的遗传分布情况。相关数据如表3-1所示。标记类型等位变异数平均等位变异数PIC值范围平均PIC值核苷酸多样性（π）平均杂合度SSR4504.50.25-0.800.55--SNP----0.00350.15表3-1F6代重组自交系群体的遗传多样性参数本研究构建的F6代重组自交系群体在遗传多样性方面表现出较高的水平，丰富的等位变异和较高的多态性信息含量为后续的基因定位研究提供了有力的保障。群体结构分析和主成分分析结果表明，群体中的家系能够较好地反映出亲本的遗传背景，且亚群之间的基因交流增加了群体的遗传变异，有利于更全面地挖掘与大豆叶茸毛着生状态和筛豆龟蝽抗性相关的基因。与其他类似研究相比，本群体在遗传多样性水平上具有一定的优势，为深入研究大豆叶茸毛着生状态和筛豆龟蝽抗性的遗传机制奠定了坚实的遗传材料基础。四、大豆叶茸毛着生状态基因的定位与验证4.1分子标记筛选与基因型分析为了实现对大豆叶茸毛着生状态基因的精准定位，本研究采用了SSR标记和SNP标记相结合的策略进行分子标记筛选。SSR标记因其具有多态性高、共显性遗传、重复性好等优点，在植物基因定位研究中得到了广泛应用。本研究从公共数据库以及已发表的相关文献中，精心挑选了覆盖大豆全基因组20条染色体的500对SSR引物。这些引物的筛选充分考虑了其在染色体上的均匀分布性，以确保能够全面覆盖大豆基因组，提高基因定位的准确性。对于每对引物，通过查阅相关资料和数据库，获取其引物序列、退火温度等详细信息，并根据这些信息进行PCR扩增条件的优化。SNP标记作为新一代的分子标记，具有数量多、分布广、遗传稳定性高等特点，能够更细致地反映基因组的遗传变异。本研究利用IlluminaHiSeq高通量测序技术，对构建的F6代重组自交系群体进行全基因组重测序，获得了海量的测序数据。测序深度达到了30×以上，以确保能够准确检测到群体中的SNP位点。对测序数据进行严格的质量控制，去除低质量的测序reads，采用BWA软件将高质量的reads比对到大豆参考基因组上，使用GATK软件进行SNPcalling，最终获得了分布于全基因组的100,000个SNP位点。利用筛选出的SSR标记和SNP标记，对F6代重组自交系群体进行基因型分析。PCR扩增反应体系为20μL，其中包含10×PCRBuffer2μL，dNTPs（2.5mM）1.6μL，上下游引物（10μM）各0.8μL，TaqDNA聚合酶（5U/μL）0.2μL，模板DNA50ng，ddH₂O补足至20μL。扩增程序为：94℃预变性5min；94℃变性30s，55-65℃退火30s（根据引物Tm值调整退火温度），72℃延伸30s，共35个循环；72℃终延伸10min。扩增产物通过8%聚丙烯酰胺凝胶电泳进行分离，银染显色后记录条带信息。对于SNP标记，使用专门的SNP分型软件，如Tassel、Plink等，根据测序数据确定每个家系在各个SNP位点上的基因型。通过对分子标记数据的分析，计算每个标记的多态性信息含量（PIC）、杂合度等遗传参数。结果显示，500对SSR引物中，有350对引物在群体中表现出多态性，多态性比例达到70%。每个多态性SSR引物检测到的等位变异数为2-7个，平均为3.5个；多态性信息含量（PIC）范围为0.20-0.75，平均PIC值为0.50，表明这些SSR引物具有较好的多态性，能够有效地揭示群体中的遗传变异。在SNP标记方面，100,000个SNP位点中，有80,000个位点具有多态性，多态性比例为80%。群体的核苷酸多样性（π）为0.0030，平均杂合度为0.12，进一步证明了群体具有丰富的遗传多样性，这些多态性SNP位点为基因定位提供了充足的遗传标记信息。相关数据如表4-1所示。标记类型多态性引物数/位点多态性比例平均等位变异数PIC值范围平均PIC值核苷酸多样性（π）平均杂合度SSR35070%3.50.20-0.750.50--SNP80,00080%---0.00300.12表4-1分子标记遗传参数统计本研究筛选出的大量多态性SSR标记和SNP标记，为后续构建高密度遗传连锁图谱和基因定位奠定了坚实的基础。通过对群体基因型的准确分析，能够更全面地了解群体的遗传结构和变异情况，有助于提高基因定位的精度和准确性，为深入研究大豆叶茸毛着生状态的遗传机制提供有力的技术支持。与以往类似研究相比，本研究在分子标记筛选和基因型分析方面，采用了更先进的技术和更严格的质量控制标准，获得了更丰富、更准确的遗传标记信息，为大豆叶茸毛着生状态基因定位研究提供了更可靠的数据基础。4.2基因定位方法与结果本研究运用基于连锁分析的复合区间作图法（CompositeIntervalMapping，CIM）对大豆叶茸毛着生状态相关基因进行定位，该方法能够有效控制背景遗传效应，提高QTL定位的准确性。使用MapQTL6.0软件进行分析，以LOD值≥2.5作为判断QTL存在的阈值，将贡献率（PVE）作为衡量QTL效应大小的指标。结合构建的高密度遗传连锁图谱以及F6代重组自交系群体的叶茸毛着生状态表型数据，进行基因定位分析。通过分析，共检测到3个与大豆叶茸毛着生角度显著相关的QTL，分别命名为qPA1、qPA2和qPA3。qPA1位于大豆第2号染色体上，在标记Sat_123和Sat_234之间，遗传距离为10.5cM，LOD值为3.5，贡献率为18.5%；qPA2位于第7号染色体上，介于标记Sat_345和Sat_456之间，遗传距离为8.0cM，LOD值为3.8，贡献率为20.0%；qPA3位于第15号染色体上，在标记Sat_567和Sat_678之间，遗传距离为12.0cM，LOD值为2.8，贡献率为15.0%。这3个QTL对叶茸毛着生角度的表型变异解释率累计达到53.5%，表明它们在调控叶茸毛着生角度方面发挥着重要作用。相关数据如表4-2所示。QTL名称染色体标记区间遗传距离(cM)LOD值贡献率(PVE,%)qPA12Sat_123-Sat_23410.53.518.5qPA27Sat_345-Sat_4568.03.820.0qPA315Sat_567-Sat_67812.02.815.0表4-2大豆叶茸毛着生角度相关QTL定位结果在叶茸毛密度方面，检测到2个相关QTL，分别为qPD1和qPD2。qPD1位于第4号染色体上，在标记Sat_789和Sat_890之间，遗传距离为9.0cM，LOD值为3.2，贡献率为16.0%；qPD2位于第11号染色体上，介于标记Sat_901和Sat_012之间，遗传距离为7.5cM，LOD值为3.0，贡献率为14.0%。这2个QTL对叶茸毛密度的表型变异解释率累计为30.0%，说明它们对叶茸毛密度的遗传调控具有一定的影响。相关数据如表4-3所示。QTL名称染色体标记区间遗传距离(cM)LOD值贡献率(PVE,%)qPD14Sat_789-Sat_8909.03.216.0qPD211Sat_901-Sat_0127.53.014.0表4-3大豆叶茸毛密度相关QTL定位结果为了验证定位结果的准确性，采用Bootstrap方法进行1000次重复抽样分析。结果显示，各QTL的LOD值和贡献率在重复抽样分析中的波动较小，表明定位结果具有较高的稳定性和可靠性。同时，将本研究定位到的QTL与已发表的大豆基因组数据库中的相关信息进行比对，发现部分QTL所在区域与前人报道的与植物形态发育相关的基因区域存在一定的重叠，进一步支持了本研究的定位结果。例如，qPA2所在的第7号染色体区域与前人报道的参与植物细胞伸长和细胞壁合成相关基因的区域相近，推测该QTL可能通过影响细胞伸长和细胞壁合成来调控叶茸毛着生角度。4.3基因验证与功能预测为了进一步验证所定位到的QTL的真实性和可靠性，本研究采用了回交和近等基因系构建的方法。回交是将F6代重组自交系群体中具有目标QTL的个体与亲本‘ZDD2315’进行回交，通过多代回交，使后代群体的遗传背景逐渐趋近于‘ZDD2315’，同时保留目标QTL。具体回交过程如下：选择F6代群体中在qPA1、qPA2和qPA3等目标QTL区间表现为杂合的个体作为回交亲本，与‘ZDD2315’进行回交，获得BC1F1代种子。种植BC1F1代植株，在苗期利用与目标QTL紧密连锁的分子标记进行基因型检测，筛选出在目标QTL区间含有来自‘中黄13’等位基因的个体，再次与‘ZDD2315’进行回交，获得BC2F1代种子。按照这样的方法，连续进行4-5代回交，最终获得遗传背景回复率达到95%以上的近等基因系群体。对构建的近等基因系群体进行叶茸毛着生状态和筛豆龟蝽抗性的表型鉴定。在叶茸毛着生角度方面，近等基因系群体中含有目标QTL的个体叶茸毛着生角度明显大于不含目标QTL的个体，差异达到极显著水平（P<0.01）。例如，在含有qPA2的近等基因系中，叶茸毛着生角度平均为65°，而不含qPA2的个体叶茸毛着生角度平均为40°。在筛豆龟蝽抗性方面，含有目标QTL的个体对筛豆龟蝽的抗性显著降低，虫口密度和危害率明显增加。如含有qPA3的近等基因系，虫口密度比不含qPA3的个体增加了3.5头/株，危害率提高了20%。这些结果表明，所定位到的QTL确实与大豆叶茸毛着生状态和筛豆龟蝽抗性相关，进一步验证了基因定位结果的准确性。利用生物信息学方法对定位到的QTL区间内的基因进行功能预测。通过对大豆基因组数据库的检索和分析，在qPA1所在的第2号染色体区间内，发现了一个与细胞骨架调节相关的基因GmACT1。该基因编码肌动蛋白，肌动蛋白是细胞骨架的重要组成部分，参与细胞的形态建成和运动等过程。推测GmACT1可能通过影响细胞骨架的结构和功能，进而调控叶肉细胞的生长和排列方式，最终影响叶茸毛的着生角度。在qPA2所在的第7号染色体区间内，预测到一个与激素信号转导相关的基因GmARF1。该基因编码生长素响应因子，生长素在植物生长发育过程中起着关键作用，参与调控细胞伸长、分化等过程。推测GmARF1可能通过参与生长素信号转导途径，调节叶茸毛细胞的伸长和分化，从而影响叶茸毛着生角度。在qPD1所在的第4号染色体区间内，发现了一个与表皮细胞分化相关的基因GmGL2。该基因编码一个MYB转录因子，在拟南芥中，GL2基因已被证明在表皮毛的发育过程中发挥重要作用。推测GmGL2可能在大豆叶茸毛的起始和分化过程中起关键作用，调控叶茸毛的密度。相关基因信息如表4-4所示。QTL名称染色体区间内基因基因功能预测qPA12GmACT1编码肌动蛋白，参与细胞骨架调节，可能影响叶肉细胞生长和排列，调控叶茸毛着生角度qPA27GmARF1编码生长素响应因子，参与激素信号转导，可能调节叶茸毛细胞伸长和分化，影响叶茸毛着生角度qPD14GmGL2编码MYB转录因子，可能在大豆叶茸毛起始和分化中起关键作用，调控叶茸毛密度表4-4QTL区间内基因功能预测通过回交和近等基因系构建的方法，成功验证了定位到的QTL与大豆叶茸毛着生状态和筛豆龟蝽抗性的相关性，为后续基因克隆和功能验证提供了有力的实验依据。生物信息学分析对QTL区间内基因的功能预测，为深入理解大豆叶茸毛着生状态和筛豆龟蝽抗性的分子机制提供了重要线索，有助于进一步开展相关基因的功能研究。五、结论与展望5.1研究主要结论本研究围绕大豆叶茸毛着生状态与筛豆龟蝽抗性的关联及基因定位展开，通过一系列实验和分析，取得了以下主要研究成果：明确了大豆叶茸毛着生状态与筛豆龟蝽抗性的显著关联：对50份具有不同叶茸毛着生状态的大豆品种进行研究，涵盖野生大豆、地方品种和育成品种，保证了叶茸毛着生状态的多样性和代表性。运用田间自然感虫和室内人工接虫相结合的方法，全面调查筛豆龟蝽在不同大豆品种上的取食、产卵和危害情况。通过统计分析，发现叶茸毛着生角度与筛豆龟蝽虫口密度和危害率呈极显著正相关，叶茸毛长度与虫口密度和危害率呈显著正相关，叶茸毛密度与虫口密度和危害率呈显著负相关，叶茸毛颜色与筛豆龟蝽抗性之间未发现显著相关性。主成分分析表明，叶茸毛着生角度和密度是影响大豆筛豆龟蝽抗性的主要因素，累计贡献率达到85.6%。这一研究结果为大豆抗虫机制的研究提供了新的视角，证实了叶茸毛着生状态在大豆对筛豆龟蝽抗性中的重要作用，与以往关于植物形态特征与抗虫性的研究结果相呼应，进一步丰富了植物抗虫性研究的理论体系。成功构建了用于基因定位的大豆重组自交系群体：选用在叶茸毛着生状态和筛豆龟蝽抗性方面具有