基于黄瓜酸黄瓜渐渗系的抗霜霉病QTLs定位及早期互作表达谱解析

基于黄瓜酸黄瓜渐渗系的抗霜霉病QTLs定位及早期互作表达谱解析一、引言1.1研究背景与意义黄瓜（CucumissativusL.）作为全球广泛种植的重要蔬菜作物，在蔬菜生产和人们日常生活中占据着举足轻重的地位。中国作为黄瓜种植和消费大国，种植面积和产量均位居世界首位。据FAO数据显示，2022年我国黄瓜收获面积达131.15万hm²，总产量高达7730.73万t，为蔬菜的稳产保供、农民增收以及乡村振兴做出了卓越贡献。黄瓜不仅营养丰富，富含多种维生素和矿物质，而且口感清脆爽口，既可生食、凉拌，也能用于烹饪和加工，深受消费者喜爱。然而，在黄瓜的种植过程中，霜霉病（Pseudoperonosporacubensis）的威胁不容小觑。黄瓜霜霉病是由鞭毛菌亚门真菌古巴假霜霉菌侵染所致，是黄瓜最常见且危害严重的病害之一。该病主要危害叶片，也能对茎和花序造成影响，从苗期至成株期均可发病。幼苗发病时，子叶正面会出现不规则水渍状褪绿黄斑，在潮湿条件下，子叶背面会产生灰黑色霉层，严重时导致子叶干垂、幼苗死亡。成株多在开花至结瓜期发病较重，常从植株下部叶片向上部叶片发展蔓延。染病初期，叶片背面形成水渍状小点，后病斑逐渐扩大，受叶脉限制呈多角形、水渍状，潮湿时病斑上生紫黑色霉层，严重时多个病斑连接成片，全叶变为黄褐色，最终干枯、死亡。黄瓜霜霉病具有传播速度快、流行范围广的特点，一旦爆发，若防治不及时，极易导致黄瓜大幅度减产，甚至绝收。据相关研究表明，在适宜的发病条件下，霜霉病可使黄瓜减产30%-50%，严重影响了黄瓜的产量和品质，给黄瓜产业带来了巨大的经济损失。目前，针对黄瓜霜霉病的防治，主要采取化学药剂防治、农业防治和生物防治等措施。化学药剂防治虽然在一定程度上能够有效控制病害的发生，但长期大量使用化学农药，不仅会导致病原菌产生抗药性，降低防治效果，还会造成环境污染，危害食品安全，对生态平衡和人体健康构成潜在威胁。农业防治措施如轮作、合理密植、加强田间管理等，虽然有助于减少病害的发生，但难以从根本上解决问题。生物防治作为一种绿色环保的防治手段，具有广阔的应用前景，但目前生物防治产品的效果还不够稳定，且成本较高，限制了其大规模应用。随着分子生物学技术的飞速发展，利用抗病分子育种技术培育抗病品种，已成为防治黄瓜霜霉病的最经济有效且可持续的方法。通过抗病分子育种，可以将抗病基因导入优良黄瓜品种中，提高黄瓜对霜霉病的抗性，减少化学农药的使用，实现黄瓜的绿色安全生产。而黄瓜-酸黄瓜渐渗系作为一种重要的抗病资源，在黄瓜抗病分子育种中具有巨大的潜力。酸黄瓜（Cucumissativusvar.hardwickiiAlef.）是黄瓜的近缘野生种，具有丰富的遗传多样性和多种优良性状，如对多种病害具有较强的抗性。通过将酸黄瓜的优良基因渐渗到栽培黄瓜中，构建黄瓜-酸黄瓜渐渗系，能够为黄瓜抗病育种提供新的种质资源和基因来源。本研究旨在对黄瓜-酸黄瓜渐渗系进行抗霜霉病QTLs定位，并深入分析其与霜霉菌早期互作的表达谱，揭示黄瓜抗霜霉病的分子机制。这不仅有助于丰富黄瓜抗霜霉病的遗传理论，为黄瓜抗病分子育种提供重要的理论依据，还能够筛选出与抗霜霉病相关的关键基因和分子标记，为黄瓜抗霜霉病新品种的培育提供有力的技术支持，推动黄瓜产业的可持续发展。1.2国内外研究现状1.2.1黄瓜霜霉病研究概况黄瓜霜霉病的病原菌为古巴假霜霉菌（Pseudoperonosporacubensis），属于鞭毛菌亚门真菌。其孢子囊梗无色，单生或2-4根束生，从气孔伸出，上部呈3-5次锐角分枝，分枝末端着生一个孢子囊；孢子囊卵形或柠檬形，顶端具乳状突起，单胞，淡褐色，可萌发长出芽管，或释放出游动孢子，变为圆形休止孢子后萌发芽管，侵入寄主。卵孢子球形，黄色，表面有瘤状突起。黄瓜霜霉病的症状较为典型，主要危害叶片。苗期发病时，子叶上起初出现褪绿斑，逐渐变为黄色不规则形斑，潮湿时子叶背面产生灰黑色霉层，随着病情发展，子叶很快变黄、枯干。成株期发病，叶片上初现浅绿色水浸斑，扩大后受叶脉限制，呈多角形，黄绿色转淡褐色，后期病斑汇合成片，全叶干枯，由叶缘向上卷缩，潮湿时叶背面病斑上生出灰黑色霉层，严重时全株叶片枯死。抗病品种病斑少而小，叶背霉层也稀疏。该病的流行规律与环境条件密切相关。黄瓜霜霉病最适宜发病温度为16-24℃，低于10℃或高于28℃较难发病，低于5℃或高于30℃基本不发病；适宜的发病湿度为85%以上，特别在叶片有水膜时，最易受侵染发病，湿度低于70%，病菌孢子难以发芽侵染，低于60%，病菌孢子不能产生。在北方，黄瓜霜霉病的传播途径主要是从温室传到大棚，再传到春季露地黄瓜上，然后传到秋季露地黄瓜上，最后又传回到温室黄瓜上，主要靠气流和雨水传播。一旦田间出现中心病株，在适宜的温湿度条件下，只需3-4次的扩大再侵染，即可酿成大灾。针对黄瓜霜霉病的防治，目前主要有以下几种手段：一是选用抗病品种，这是防治霜霉病最经济有效的方法之一，但目前高抗霜霉病的黄瓜品种相对较少；二是加强栽培管理，如采用高畦栽培、合理密植、控制浇水、及时通风排湿等，创造不利于病害发生的环境条件；三是物理防治，如选用无滴膜，预防棚膜滴水增加棚室湿度，发病严重药剂控制较难的温棚，可采用高温闷棚的方法；四是化学防治，在发病初期及时喷药保护，可选用烯酰吗啉、霜脲・锰锌、甲霜灵・锰锌、霜霉威盐酸盐、代森锰锌等药剂，但长期使用化学药剂易导致病原菌产生抗药性；五是生物防治，利用有益微生物或其代谢产物来抑制病原菌的生长和繁殖，但生物防治的效果受环境因素影响较大，且作用速度相对较慢。1.2.2黄瓜酸黄瓜渐渗系研究进展渐渗系（Introgressionline）是指通过杂交、回交等手段，将野生种或近缘种的染色体片段导入到栽培种中，从而形成的一系列含有不同外源染色体片段的后代群体。在黄瓜遗传改良中，黄瓜-酸黄瓜渐渗系具有重要作用。酸黄瓜作为黄瓜的近缘野生种，携带了许多栽培黄瓜所不具备的优良基因，如对多种病虫害的抗性基因、抗逆基因等。通过构建黄瓜-酸黄瓜渐渗系，可以将这些优良基因引入栽培黄瓜中，拓宽黄瓜的遗传基础，为黄瓜新品种的选育提供丰富的种质资源。目前，国内外在黄瓜-酸黄瓜渐渗系的研究方面已取得了一定成果。在种质资源创新方面，科研人员通过杂交和回交等方法，成功构建了多个黄瓜-酸黄瓜渐渗系群体，并对这些渐渗系的农艺性状、抗病性等进行了鉴定和评价。例如，有研究利用酸黄瓜与栽培黄瓜杂交，获得了一系列具有不同抗性表现的渐渗系材料，为黄瓜抗病育种提供了新的种质资源。在基因挖掘方面，借助分子标记技术和基因组测序技术，对渐渗系中导入的酸黄瓜染色体片段进行了定位和分析，成功挖掘出了一些与抗病、抗逆等性状相关的基因。例如，通过对黄瓜-酸黄瓜渐渗系进行抗蔓枯病基因定位，发现了多个与抗蔓枯病相关的QTLs。然而，目前对于黄瓜-酸黄瓜渐渗系抗霜霉病的研究还相对较少，尤其是在QTLs定位和与霜霉菌互作的分子机制方面，仍有待深入探究。1.2.3QTLs定位技术在植物抗病研究中的应用QTLs（QuantitativeTraitLoci）即数量性状位点，是指控制数量性状的基因在基因组中的位置。QTLs定位的原理是基于遗传连锁分析，通过构建遗传图谱，利用分子标记与目标性状之间的连锁关系，确定控制目标性状的QTLs在染色体上的位置和效应。常用的QTLs定位方法包括区间作图法、复合区间作图法、多区间作图法等。区间作图法是利用相邻标记构建区间，通过比较区间内标记基因型与表型的关系来确定QTLs的位置；复合区间作图法在区间作图的基础上，进一步考虑了其他QTLs的影响，提高了定位的准确性；多区间作图法则可以同时分析多个区间，更全面地检测QTLs。在黄瓜抗病研究中，QTLs定位技术已被广泛应用于挖掘抗病基因。例如，有研究通过对黄瓜白粉病抗性材料进行QTLs定位，在黄瓜7条染色体上均检测到了与叶片白粉病抗性相关的位点。在其他植物抗病研究中，QTLs定位也取得了众多成果。在番茄抗叶霉病研究中，利用QTLs定位技术定位到了多个与抗叶霉病相关的基因位点，并通过分子标记辅助选择，将这些抗病基因聚合到优良品种中，培育出了高抗叶霉病的番茄新品种。在水稻抗稻瘟病研究中，通过QTLs定位鉴定出了多个抗稻瘟病的QTLs，并对其中一些关键基因进行了克隆和功能验证，为水稻抗稻瘟病分子育种提供了重要的基因资源和理论依据。1.2.4植物与病原菌互作的表达谱分析研究进展植物与病原菌互作是一个复杂的生物学过程，涉及到植物和病原菌双方基因的表达调控。当植物受到病原菌侵染时，会启动一系列的防御反应，包括细胞壁加厚、产生植保素、激活防御相关基因的表达等。而病原菌为了成功侵染植物，也会分泌各种效应分子，抑制植物的防御反应。表达谱分析技术是研究植物与病原菌互作分子机制的重要手段之一。常用的表达谱分析技术包括基因芯片技术和RNA-Seq技术。基因芯片技术是将大量的基因探针固定在芯片上，通过与样本中的RNA杂交，检测基因的表达水平；RNA-Seq技术则是利用高通量测序技术，对样本中的RNA进行测序，从而全面、准确地分析基因的表达情况。通过表达谱分析，可以筛选出在植物与病原菌互作过程中差异表达的基因，进而深入研究这些基因的功能和调控机制。例如，在小麦与条锈菌互作的研究中，利用RNA-Seq技术分析了互作过程中基因的表达变化，发现了许多与小麦抗病相关的基因，包括一些参与信号转导、防御反应等过程的基因。在拟南芥与病原菌互作的研究中，通过基因芯片技术筛选出了一系列受病原菌诱导表达的基因，并对这些基因的启动子区域进行分析，揭示了其调控机制。1.3研究目标与内容本研究的目标是通过对黄瓜-酸黄瓜渐渗系进行抗霜霉病QTLs定位，明确黄瓜抗霜霉病的遗传基础，并分析其与霜霉菌早期互作的表达谱，筛选出与抗霜霉病相关的关键基因，为黄瓜抗霜霉病分子育种提供理论依据和基因资源。具体研究内容如下：构建黄瓜-酸黄瓜渐渗系群体，并对其进行抗霜霉病鉴定，筛选出高抗霜霉病的渐渗系材料。利用酸黄瓜与栽培黄瓜进行杂交和回交，构建包含多个世代的渐渗系群体。采用人工接种霜霉菌的方法，对渐渗系群体进行抗霜霉病鉴定，根据发病症状和病情指数，筛选出抗霜霉病表现优异的材料，为后续研究提供材料基础。利用筛选出的抗霜霉病渐渗系材料，构建遗传图谱，进行抗霜霉病QTLs定位。提取渐渗系群体及双亲的DNA，利用SSR、SNP等分子标记技术，构建高密度的遗传图谱。结合抗霜霉病表型数据，运用QTL定位软件，如MapQTL、QTLIciMapping等，进行抗霜霉病QTLs定位，确定QTLs在染色体上的位置、效应和遗传贡献率。分析抗霜霉病渐渗系与霜霉菌早期互作的表达谱，筛选差异表达基因。在霜霉菌接种后的不同时间点，采集抗霜霉病渐渗系和感病对照的叶片样本，提取RNA，利用RNA-Seq技术进行表达谱分析。通过生物信息学分析，筛选出在抗霜霉病渐渗系与霜霉菌互作过程中差异表达的基因，并对这些基因进行功能注释和富集分析，初步揭示黄瓜抗霜霉病的分子机制。对筛选出的差异表达基因进行验证和功能分析。采用实时荧光定量PCR技术，对RNA-Seq结果进行验证，确保差异表达基因的可靠性。选取部分关键的差异表达基因，利用基因编辑技术（如CRISPR/Cas9）进行功能验证，明确其在黄瓜抗霜霉病过程中的作用机制。1.4研究方法与技术路线1.4.1研究方法黄瓜-酸黄瓜渐渗系群体构建：选择具有代表性的酸黄瓜和栽培黄瓜品种作为亲本，进行杂交获得F1代，然后以栽培黄瓜为轮回亲本，与F1代进行连续回交，并结合自交，构建包含多个世代的渐渗系群体。抗霜霉病鉴定：采用人工接种霜霉菌的方法，在温室或人工气候箱中对渐渗系群体进行抗霜霉病鉴定。接种后，定期观察发病症状，记录病情指数，病情指数计算公式为：病情指数=∑（各级病株数×相对级数值）/（调查总株数×最高级数值）×100。DNA提取与分子标记分析：采用CTAB法提取渐渗系群体及双亲的基因组DNA，利用SSR、SNP等分子标记技术，对DNA进行多态性分析。SSR标记可通过查阅相关文献或数据库获得引物序列，进行PCR扩增后，利用聚丙烯酰胺凝胶电泳或毛细管电泳检测扩增产物的多态性；SNP标记可利用高通量测序技术进行检测，如SNP芯片技术或全基因组重测序技术。QTLs定位：利用JoinMap软件构建遗传图谱，结合抗霜霉病表型数据，运用MapQTL、QTLIciMapping等软件进行QTLs定位分析。采用复合区间作图法或多区间作图法，确定QTLs在染色体上的位置、效应和遗传贡献率，以LOD值大于2.5作为QTLs存在的阈值。RNA提取与表达谱分析：在霜霉菌接种后的0h、6h、12h、24h、48h等时间点，采集抗霜霉病渐渗系和感病对照的叶片样本，采用Trizol法提取总RNA，利用RNA-Seq技术进行表达谱分析。通过与黄瓜参考基因组进行比对，分析基因的表达水平，筛选出差异表达基因。差异表达基因的筛选标准为：|log2（foldchange）|≥1且FDR＜0.05。基因功能验证：采用实时荧光定量PCR技术对RNA-Seq结果进行验证，根据验证结果选取部分关键的差异表达基因，利用CRISPR/Cas9基因编辑技术构建基因敲除载体，通过农杆菌介导转化法将载体导入黄瓜植株中，获得基因敲除突变体。对突变体进行抗霜霉病鉴定，分析其抗病性变化，从而验证基因的功能。1.4.2技术路线本研究的技术路线如图1所示：材料准备：选择合适的酸黄瓜和栽培黄瓜品种作为亲本，进行种植和管理，为后续杂交和回交做准备。同时，准备霜霉菌菌种，进行活化和扩繁，用于抗霜霉病鉴定。渐渗系群体构建与抗霜霉病鉴定：按照上述方法进行杂交、回交和自交，构建渐渗系群体。对渐渗系群体进行人工接种霜霉菌，鉴定抗霜霉病能力，筛选出高抗霜霉病的渐渗系材料。DNA提取与分子标记分析：提取渐渗系群体及双亲的DNA，利用SSR、SNP等分子标记进行多态性分析，构建遗传图谱。QTLs定位：结合抗霜霉病表型数据和遗传图谱，运用QTL定位软件进行QTLs定位分析，确定抗霜霉病QTLs的位置、效应和遗传贡献率。RNA提取与表达谱分析：在霜霉菌接种后的不同时间点采集叶片样本，提取RNA，进行RNA-Seq分析，筛选差异表达基因。基因功能验证：采用实时荧光定量PCR技术验证差异表达基因，选取关键基因利用CRISPR/Cas9技术进行基因敲除，获得突变体，对突变体进行抗霜霉病鉴定，验证基因功能。[此处插入技术路线流程图]图1技术路线图[此处插入技术路线流程图]图1技术路线图图1技术路线图二、材料与方法2.1实验材料本研究所用的黄瓜-酸黄瓜渐渗系材料，是以栽培黄瓜（CucumissativusL.）品种‘津优2号’作为母本，酸黄瓜（Cucumissativusvar.hardwickiiAlef.）作为父本，通过杂交、回交和自交等一系列育种手段构建而成。‘津优2号’是生产上广泛应用的优良黄瓜品种，具有生长势强、产量高、品质好等优点，但对霜霉病的抗性较弱。酸黄瓜作为黄瓜的近缘野生种，在长期的自然选择过程中，积累了丰富的遗传多样性，对多种病害，尤其是霜霉病表现出较强的抗性，同时还具备一些其他优良性状，如抗逆性强等。本研究构建的黄瓜-酸黄瓜渐渗系群体包含多个世代，涵盖了丰富的遗传变异，为后续的抗霜霉病研究提供了充足的材料基础。实验中使用的霜霉菌菌株（Pseudoperonosporacubensis）采自天津市西青区黄瓜种植基地发病严重的黄瓜植株。在田间，仔细挑选具有典型霜霉病症状的叶片，这些叶片上病斑呈现多角形，水渍状，背面有紫黑色霉层。将采集到的病叶迅速装入无菌自封袋中，低温保存，尽快带回实验室。在实验室中，采用活体植株保存法对霜霉菌菌株进行保存。具体操作如下：选择感病的黄瓜品种‘津优2号’幼苗，在温室中培养至3-4片真叶期。将采集的霜霉病菌病叶上的孢子囊，用无菌水冲洗下来，制成孢子囊悬浮液，浓度调整为1×10⁴-1×10⁵个孢子囊/mL。使用小型手持喷雾器，将孢子囊悬浮液均匀地喷洒在‘津优2号’幼苗叶片的正反两面，至叶片表面欲滴为止。将接种后的幼苗置于温度为18-22℃、相对湿度90%-100%的培养箱中，黑暗保湿24h，之后每日光照14h，继续培养4-5d，待叶片上出现明显的霜霉病症状后，即可作为保存的菌株材料。每隔3-4周，将保存的菌株转接至新的感病黄瓜幼苗上，以保持菌株的活性和致病性。2.2实验方法2.2.1黄瓜酸黄瓜渐渗系群体构建以栽培黄瓜‘津优2号’为母本，酸黄瓜为父本进行杂交。在杂交过程中，于黄瓜开花前一天下午，选取‘津优2号’母本植株上即将开放的雌花，去雄后套袋隔离。次日上午，采集酸黄瓜父本植株上刚开放的雄花，将花粉轻轻涂抹在母本雌花的柱头上，授粉后再次套袋，并做好标记。待杂交果实成熟后，收获杂交种子，即F1代种子。将F1代种子播种于温室中，待植株生长至开花期，以‘津优2号’为轮回亲本，与F1代进行回交。回交过程同样进行严格的去雄、授粉和套袋操作，收获回交一代（BC1F1）种子。之后，对BC1F1代植株进行自交，获得BC1F2代种子。按照此方法，连续进行回交和自交，构建包含多个世代的黄瓜-酸黄瓜渐渗系群体。本研究构建的渐渗系群体规模为300个单株。群体规模的确定主要依据统计学原理和前人的研究经验。从统计学角度来看，较大的群体规模能够更全面地涵盖双亲的遗传变异，提高QTL定位的准确性和可靠性。参考相关研究，在植物QTL定位研究中，一般认为300个左右单株的群体规模能够满足定位分析的要求，可有效降低假阳性和假阴性结果的出现概率。在田间种植管理方面，将渐渗系群体及双亲种植于天津市农业科学院蔬菜研究所试验基地。土壤为壤土，肥力中等，排灌方便。播种前，对土壤进行深耕、耙平，并施足基肥，基肥以有机肥为主，配合适量的复合肥。采用高畦栽培，畦宽1.2m，沟宽0.3m，每畦种植两行，株距0.3m。播种后，及时浇水，保持土壤湿润，促进种子发芽和幼苗生长。在生长期间，根据黄瓜的生长需求，合理进行追肥、浇水、中耕除草等管理措施。定期喷施杀虫剂和杀菌剂，防治病虫害的发生，确保植株正常生长。2.2.2霜霉病抗性鉴定采用人工接种霜霉菌的方法对黄瓜-酸黄瓜渐渗系群体进行霜霉病抗性鉴定。在温室中，将保存的霜霉菌菌株转接至感病黄瓜品种‘津优2号’幼苗上进行扩繁。待叶片上产生大量新鲜的孢子囊后，用无菌水冲洗叶片，收集孢子囊，制成孢子囊悬浮液。使用血球计数板对孢子囊悬浮液进行计数，并将浓度调整为5×10⁴个孢子囊/mL，同时加入0.05%的吐温-20作为表面活性剂，以增强孢子囊悬浮液在叶片表面的附着性。当渐渗系群体及双亲的黄瓜幼苗长至3-4片真叶期时，进行接种。采用喷雾接种法，使用小型手持喷雾器将配制好的孢子囊悬浮液均匀地喷洒在黄瓜幼苗叶片的正反两面，至叶片表面欲滴为止。接种后，将黄瓜幼苗置于温度为18-22℃、相对湿度90%-100%的培养箱中，黑暗保湿24h，以促进孢子囊的萌发和侵染。之后，将幼苗转移至温室中，保持温度为20-25℃，每日光照14h，正常管理。抗性评价分级标准参考前人的研究并结合本实验的实际情况进行制定，具体如下：0级，叶片无病斑，表现为免疫；1级，叶片病斑面积占叶面积的10%以下，病情轻微，为高抗；3级，叶片病斑面积占叶面积的10%-25%，有少量病斑，属于抗病；5级，叶片病斑面积占叶面积的25%-50%，病斑较多，为中抗；7级，叶片病斑面积占叶面积的50%-75%，病斑密布，属于感病；9级，叶片病斑面积占叶面积的75%以上，或叶片枯死，为高感。在接种后的第5d、7d、10d分别进行调查，记录每株黄瓜幼苗叶片的发病情况，按照上述分级标准进行病情分级，并计算病情指数。病情指数计算公式为：病情指数=∑（各级病株数×相对级数值）/（调查总株数×最高级数值）×100。以3次调查的病情指数平均值作为最终的抗性评价依据。2.2.3DNA提取与分子标记筛选采用CTAB法从黄瓜叶片中提取基因组DNA。具体步骤如下：取0.2g黄瓜幼嫩叶片，置于预冷的研钵中，加入适量液氮，迅速研磨成粉末状。将研磨好的叶片粉末转移至1.5mL离心管中，加入700μL预热至65℃的CTAB提取缓冲液（100mmol/LTris-HCl，pH8.0；20mmol/LEDTA，pH8.0；1.4mol/LNaCl；2%CTAB；0.2%β-巯基乙醇，使用前加入），轻轻颠倒混匀，使叶片粉末与提取缓冲液充分接触。将离心管置于65℃水浴锅中保温30-60min，期间每隔10min轻轻颠倒混匀一次，以促进DNA的释放和溶解。保温结束后，取出离心管，冷却至室温，加入等体积的氯仿-异戊醇（24:1）混合液，轻轻颠倒混匀10-15min，使溶液充分乳化。12000rpm离心15min，将上清液转移至新的1.5mL离心管中。加入2/3体积预冷的异丙醇，轻轻颠倒混匀，可见白色絮状DNA沉淀析出。12000rpm离心10min，弃上清液，用70%乙醇洗涤DNA沉淀2-3次，每次洗涤后12000rpm离心5min，弃去乙醇。将DNA沉淀在室温下晾干，加入50μLTE缓冲液（10mmol/LTris-HCl，pH8.0；1mmol/LEDTA，pH8.0）溶解DNA，置于-20℃冰箱中保存备用。利用SSR（SimpleSequenceRepeat）和SNP（SingleNucleotidePolymorphism）分子标记对黄瓜-酸黄瓜渐渗系群体及双亲的DNA进行多态性筛选。SSR标记引物序列参考相关文献和数据库，由上海生工生物工程有限公司合成。PCR扩增体系为20μL，包括10×PCRbuffer2μL，2.5mmol/LdNTPs1.6μL，10μmol/L上下游引物各0.8μL，5U/μLTaqDNA聚合酶0.2μL，模板DNA2μL，ddH₂O12.6μL。PCR扩增程序为：94℃预变性5min；94℃变性30s，55-65℃退火30s（根据引物Tm值调整退火温度），72℃延伸30s，共35个循环；72℃延伸10min。扩增产物在6%的聚丙烯酰胺凝胶上进行电泳分离，银染显色，观察并记录多态性条带。SNP标记筛选采用高通量测序技术，将渐渗系群体及双亲的DNA进行全基因组重测序。测序数据经过质量控制和比对分析后，利用生物信息学软件（如GATK、SAMtools等）进行SNP位点的检测和筛选。筛选标准为：质量值大于20，深度大于10，最小等位基因频率大于0.05。2.2.4QTLs定位分析利用筛选出的多态性SSR和SNP分子标记，采用JoinMap4.0软件构建黄瓜-酸黄瓜渐渗系群体的遗传图谱。在构建遗传图谱时，将分子标记按照连锁群进行排列，计算标记之间的遗传距离，以cM（centi-Morgan）为单位。遗传距离的计算采用Kosambi函数，通过多次迭代计算，使遗传图谱的标记顺序和遗传距离更加准确。结合霜霉病抗性鉴定获得的表型数据，运用QTLIciMapping4.2软件进行QTLs定位分析。采用复合区间作图法（CompositeIntervalMapping，CIM）进行QTLs的检测，以LOD（LogarithmofOdds）值大于2.5作为QTLs存在的阈值。在分析过程中，设置控制背景遗传效应的标记数为5，窗口大小为10cM，步长为1cM。通过QTLs定位分析，确定抗霜霉病QTLs在染色体上的位置、效应和遗传贡献率。为了验证QTLs定位结果的可靠性，采用以下方法进行验证：一是利用不同年份、不同环境条件下的霜霉病抗性鉴定数据，对定位到的QTLs进行重复性验证；二是选择与定位到的QTLs紧密连锁的分子标记，在不同的黄瓜-酸黄瓜渐渗系材料中进行验证，观察标记与抗性表型的共分离情况；三是通过构建近等基因系（Near-IsogenicLines，NILs），进一步验证QTLs的效应和功能。2.2.5RNA提取与表达谱测序在霜霉菌接种后的0h（未接种对照）、6h、12h、24h、48h，分别采集抗霜霉病渐渗系和感病对照（‘津优2号’）的叶片样本，每个时间点选取3株植株，每株植株采集第3-4片真叶，混合作为一个生物学重复，每个处理设置3个生物学重复。采用Trizol法提取叶片总RNA，具体步骤如下：取0.1-0.2g叶片组织，置于预冷的研钵中，加入适量液氮，迅速研磨成粉末状。将研磨好的叶片粉末转移至1.5mL离心管中，加入1mLTrizol试剂，剧烈振荡混匀，室温静置5min，使细胞充分裂解。加入200μL氯仿，剧烈振荡15s，室温静置3min。12000rpm离心15min，将上层水相转移至新的1.5mL离心管中。加入等体积的异丙醇，轻轻颠倒混匀，室温静置10min。12000rpm离心10min，弃上清液，用75%乙醇洗涤RNA沉淀2-3次，每次洗涤后12000rpm离心5min，弃去乙醇。将RNA沉淀在室温下晾干，加入适量的RNase-free水溶解RNA，置于-80℃冰箱中保存备用。利用提取的总RNA进行文库构建和测序。文库构建采用IlluminaTruSeqRNASamplePreparationKit试剂盒，按照说明书进行操作。首先，利用Oligo(dT)磁珠富集mRNA，然后将mRNA打断成短片段，以短片段mRNA为模板，合成cDNA第一链和第二链。对cDNA进行末端修复、加A尾、连接测序接头等处理后，通过PCR扩增富集文库片段。文库构建完成后，使用Agilent2100Bioanalyzer对文库的质量和浓度进行检测，确保文库的质量符合测序要求。将合格的文库在IlluminaHiSeq2500测序平台上进行双端测序，测序读长为150bp。测序完成后，对原始数据进行质量控制和过滤，去除低质量reads、接头序列和含N比例过高的reads。将过滤后的高质量reads与黄瓜参考基因组（CucumberGenomeDatabase，/）进行比对，使用Hisat2软件进行比对分析，统计每个基因的reads数和表达量。表达量的计算采用FPKM（FragmentsPerKilobaseofexonperMillionreadsmapped）法，公式为：FPKM=10⁶×C/(N×L/10³)，其中C为比对到某基因的reads数，N为比对到所有基因的总reads数，L为该基因的外显子长度（bp）。通过比较抗霜霉病渐渗系和感病对照在不同时间点的基因表达量，筛选出差异表达基因，筛选标准为：|log2（foldchange）|≥1且FDR（FalseDiscoveryRate）＜0.05。2.2.6差异表达基因的功能注释与富集分析利用多个数据库对筛选出的差异表达基因进行功能注释，包括NCBI（NationalCenterforBiotechnologyInformation）的Nr数据库、Swiss-Prot数据库、KEGG（KyotoEncyclopediaofGenesandGenomes）数据库和GO（GeneOntology）数据库。使用Blast软件将差异表达基因的氨基酸序列与数据库中的序列进行比对，E-value阈值设置为1e-5。根据比对结果，获取差异表达基因的功能注释信息，包括基因的功能描述、所属的代谢途径和生物学过程等。对差异表达基因进行GO功能富集分析和KEGG代谢通路富集分析，以揭示差异表达基因在生物学过程、细胞组成和分子功能等方面的富集情况，以及参与的主要代谢途径和信号转导通路。GO功能富集分析使用clusterProfilerR包，通过超几何检验计算每个GOterm的富集显著性，以P＜0.05作为富集的阈值。KEGG代谢通路富集分析同样使用clusterProfilerR包，通过Fisher精确检验计算每个KEGGpathway的富集显著性，以P＜0.05作为富集的阈值。通过富集分析，找出与黄瓜抗霜霉病相关的关键生物学过程、分子功能和代谢通路，为进一步研究黄瓜抗霜霉病的分子机制提供线索。2.2.7实时荧光定量PCR验证从差异表达基因中挑选10-15个基因进行实时荧光定量PCR（qRT-PCR）验证。挑选验证基因的原则主要包括：一是在表达谱测序结果中差异表达倍数较高，且FDR值较小，表明这些基因在抗霜霉病过程中可能发挥重要作用；二是选择参与已知抗病相关途径或具有潜在抗病功能的基因，如参与植物激素信号转导、防御反应、细胞壁合成等途径的基因；三是考虑基因在不同组织和发育阶段的表达特异性，优先选择在叶片中高表达且与抗霜霉病相关的基因。根据挑选的差异表达基因序列，利用PrimerPremier5.0软件设计qRT-PCR引物。引物设计的基本原则包括：引物长度为18-25bp，GC含量在40%-60%之间，Tm值在58-62℃之间，引物3'端避免出现连续的3个以上相同碱基，且引物之间不能形成二聚体和发夹结构。引物由上海生工生物工程有限公司合成。qRT-PCR反应体系为20μL，包括2×SYBRGreenMasterMix10μL，10μmol/L上下游引物各0.8μL，模板cDNA2μL，ddH₂O6.4μL。反应程序为：95℃预变性30s；95℃变性5s，60℃退火30s，72℃延伸30s，共40个循环；72℃延伸10min。反应结束后，进行熔解曲线分析，以验证扩增产物的特异性。每个样品设置3个技术重复。采用2⁻ΔΔCt法计算基因的相对表达量。首先，以黄瓜的Actin基因作为内参基因，计算目的基因和内参基因的Ct值。然后，计算每个样品中目的基因相对于内参基因三、黄瓜酸黄瓜渐渗系抗霜霉病QTLs定位结果与分析3.1霜霉病抗性表型数据统计与分析对黄瓜-酸黄瓜渐渗系群体及双亲进行霜霉病抗性鉴定，统计分析其病情指数，结果如表1所示。母本‘津优2号’的病情指数高达78.56，表现为高感霜霉病；父本酸黄瓜的病情指数为15.23，表现为高抗霜霉病。渐渗系群体的病情指数呈现连续分布，范围在12.35-85.67之间，平均值为46.58，表现出明显的数量性状遗传特征。通过Shapiro-Wilk检验对病情指数数据进行正态性检验，结果显示P值为0.785，大于0.05，表明病情指数数据符合正态分布，这进一步说明霜霉病抗性受多基因控制，符合数量性状的遗传特点。[此处插入霜霉病抗性表型数据统计分析表]表1黄瓜-酸黄瓜渐渗系群体及双亲霜霉病抗性表型数据统计分析[此处插入霜霉病抗性表型数据统计分析表]表1黄瓜-酸黄瓜渐渗系群体及双亲霜霉病抗性表型数据统计分析表1黄瓜-酸黄瓜渐渗系群体及双亲霜霉病抗性表型数据统计分析材料病情指数平均值±标准差最小值最大值变异系数（%）偏度峰度Shapiro-Wilk检验P值‘津优2号’78.56±5.6772.344-0.560.856酸黄瓜15.23±3.4510.2320.1222.650.450.230.921渐渗系群体46.58±18.5612.3585.6739.85-0.23-0.120.785对渐渗系群体病情指数的频率分布进行分析，结果如图2所示。从图中可以看出，病情指数的分布呈现出近似正态分布的特征，在病情指数为40-50区间的植株数量最多，表明该区间的抗性水平在群体中较为常见。同时，在低病情指数（高抗）和高病情指数（高感）区间也有一定数量的植株分布，这为后续筛选高抗和高感材料提供了可能。群体中病情指数的变异系数为39.85%，说明渐渗系群体在霜霉病抗性上存在丰富的遗传变异，这为QTLs定位提供了良好的材料基础。[此处插入霜霉病抗性表型数据频率分布图]图2黄瓜-酸黄瓜渐渗系群体霜霉病抗性表型数据频率分布图[此处插入霜霉病抗性表型数据频率分布图]图2黄瓜-酸黄瓜渐渗系群体霜霉病抗性表型数据频率分布图图2黄瓜-酸黄瓜渐渗系群体霜霉病抗性表型数据频率分布图3.2遗传图谱的构建利用筛选出的多态性SSR和SNP分子标记，采用JoinMap4.0软件构建黄瓜-酸黄瓜渐渗系群体的遗传图谱。最终构建的遗传图谱包含7个连锁群，对应黄瓜的7条染色体，共包含256个分子标记，其中SSR标记156个，SNP标记100个。遗传图谱的总长度为1056.8cM，平均图距为4.13cM。各连锁群的长度在112.5-205.6cM之间，标记数量在28-45个之间，平均图距在3.98-4.56cM之间。各连锁群的详细信息如表2所示。[此处插入遗传图谱连锁群信息表]表2黄瓜-酸黄瓜渐渗系群体遗传图谱连锁群信息[此处插入遗传图谱连锁群信息表]表2黄瓜-酸黄瓜渐渗系群体遗传图谱连锁群信息表2黄瓜-酸黄瓜渐渗系群体遗传图谱连锁群信息连锁群染色体标记数量连锁群长度（cM）平均图距（cM）LG1Chr145205.64.56LG2Chr238165.44.35LG3Chr332135.84.24LG4Chr428112.53.98LG5Chr540178.24.46LG6Chr635148.74.25LG7Chr738110.62.91遗传图谱的构建质量直接影响QTLs定位的准确性和可靠性。本研究构建的遗传图谱覆盖了黄瓜的整个基因组，标记分布较为均匀，平均图距较小，能够满足QTLs定位的要求。同时，通过与已发表的黄瓜遗传图谱进行比较，验证了本研究构建的遗传图谱的准确性和可靠性。例如，在与某已发表的黄瓜高密度遗传图谱对比时，发现本研究构建的遗传图谱中大部分标记的顺序和位置与该图谱基本一致，仅有少数标记存在差异，这可能是由于不同研究中所用材料和标记筛选方法的差异导致的。3.3QTLs定位结果结合霜霉病抗性鉴定获得的表型数据，运用QTLIciMapping4.2软件进行QTLs定位分析。采用复合区间作图法，以LOD值大于2.5作为QTLs存在的阈值，共检测到5个与黄瓜霜霉病抗性相关的QTLs，分别位于Chr1、Chr3、Chr5和Chr7染色体上。这些QTLs的详细信息如表3所示。[此处插入抗霜霉病QTLs定位结果表]表3黄瓜-酸黄瓜渐渗系群体抗霜霉病QTLs定位结果[此处插入抗霜霉病QTLs定位结果表]表3黄瓜-酸黄瓜渐渗系群体抗霜霉病QTLs定位结果表3黄瓜-酸黄瓜渐渗系群体抗霜霉病QTLs定位结果QTL名称染色体位置（cM）LOD值加性效应显性效应贡献率（%）qDM1.1Chr145.6-52.33.25-6.542.3412.56qDM3.1Chr378.5-85.62.87-5.231.8710.23qDM5.1Chr532.4-40.13.56-7.212.5615.67qDM5.2Chr555.6-62.33.02-6.122.1113.45qDM7.1Chr725.3-32.12.68-4.891.678.98从表3中可以看出，5个QTLs的LOD值在2.68-3.56之间，加性效应均为负值，表明这些QTLs的增效等位基因均来自抗霜霉病的酸黄瓜亲本。贡献率在8.98%-15.67%之间，其中位于Chr5染色体上的qDM5.1贡献率最高，达到15.67%，说明该QTL对黄瓜霜霉病抗性的影响较大，是一个主效QTL。各QTL的加性效应和显性效应表明，黄瓜霜霉病抗性既受到加性基因效应的影响，也受到显性基因效应的影响。为了更直观地展示QTLs在染色体上的位置，绘制了QTLs在染色体上的分布示意图，如图3所示。从图中可以清晰地看到5个QTLs在不同染色体上的分布情况，以及它们所在的区间。这为进一步研究这些QTLs的功能和相互作用提供了便利。[此处插入QTLs在染色体上的分布示意图]图3黄瓜-酸黄瓜渐渗系群体抗霜霉病QTLs在染色体上的分布示意图[此处插入QTLs在染色体上的分布示意图]图3黄瓜-酸黄瓜渐渗系群体抗霜霉病QTLs在染色体上的分布示意图图3黄瓜-酸黄瓜渐渗系群体抗霜霉病QTLs在染色体上的分布示意图3.4QTLs的验证与分析为了验证QTLs定位结果的可靠性，采用不同年份、不同环境条件下的霜霉病抗性鉴定数据，对定位到的QTLs进行重复性验证。在2023年和2024年，分别在温室和露地条件下对黄瓜-酸黄瓜渐渗系群体进行霜霉病抗性鉴定，并利用相同的分子标记和分析方法进行QTLs定位。结果显示，在不同年份和环境条件下，均检测到了与2022年定位结果相同或相近的QTLs，表明这些QTLs具有较好的稳定性。选择与定位到的QTLs紧密连锁的分子标记，在不同的黄瓜-酸黄瓜渐渗系材料中进行验证，观察标记与抗性表型的共分离情况。例如，对于位于Chr5染色体上的qDM5.1，选取了与其紧密连锁的SSR标记SSR5-1和SNP标记SNP5-2，在10份不同的黄瓜-酸黄瓜渐渗系材料中进行验证。结果发现，在抗霜霉病的材料中，均检测到了与抗霜霉病相关的标记基因型，而在感病材料中，则检测到了与感病相关的标记基因型，表明这些标记与霜霉病抗性表型具有良好的共分离关系，进一步验证了QTLs定位结果的可靠性。通过构建近等基因系，进一步验证QTLs的效应和功能。以含有qDM5.1的黄瓜-酸黄瓜渐渗系材料为供体，以感病的栽培黄瓜品种‘津优2号’为轮回亲本，进行多代回交和自交，构建了近等基因系NIL-qDM5.1。对NIL-qDM5.1和‘津优2号’进行霜霉病抗性鉴定，结果显示，NIL-qDM5.1的病情指数显著低于‘津优2号’，表明qDM5.1具有显著的抗霜霉病效应，能够有效提高黄瓜对霜霉病的抗性。将本研究定位到的QTLs与前人的研究结果进行比较，发现部分QTLs的位置和效应与前人报道的结果具有一定的相似性。例如，本研究在Chr5染色体上定位到的qDM5.1和qDM5.2，与前人研究中在该染色体上定位到的抗霜霉病QTLs位置相近，且贡献率也较为接近。这表明这些QTLs可能是黄瓜抗霜霉病的重要位点，在不同的研究中具有一定的保守性。然而，也有一些QTLs是本研究首次报道的，这可能是由于所用材料、定位方法和环境条件等因素的差异导致的。这些新发现的QTLs为黄瓜抗霜霉病遗传研究提供了新的线索，丰富了黄瓜抗霜霉病的遗传资源。四、黄瓜酸黄瓜渐渗系与霜霉菌早期互作的表达谱分析结果4.1表达谱测序数据质量评估对霜霉菌接种后不同时间点的抗霜霉病渐渗系和感病对照的叶片样本进行RNA-Seq测序，共获得18个测序文库（每个处理3个生物学重复，6个时间点）。测序数据质量评估结果显示，各文库的原始数据产量均在5Gb以上，Q30碱基百分比均大于90%，表明测序数据质量较高，能够满足后续分析的要求。通过去除低质量reads、接头序列和含N比例过高的reads，得到高质量的cleanreads。将cleanreads与黄瓜参考基因组进行比对，比对率在85%-92%之间，平均比对率为88.5%，说明大部分reads能够准确比对到参考基因组上，数据可靠性高。各文库的测序数据质量评估详细信息如表4所示。[此处插入表达谱测序数据质量评估表]表4表达谱测序数据质量评估[此处插入表达谱测序数据质量评估表]表4表达谱测序数据质量评估表4表达谱测序数据质量评估文库编号样本时间点原始reads数原始碱基（Gb）Q30（%）GC含量（%）cleanreads数比对率（%）LIB1抗霜霉病渐渗系0h352456785.2991.2545.683256789087.56LIB2抗霜霉病渐渗系0h348765435.2390.8745.453212345686.89LIB3抗霜霉病渐渗系0h350123455.2591.0245.563234567887.23LIB4抗霜霉病渐渗系6h365432105.4892.0146.233365432188.98LIB5抗霜霉病渐渗系6h361234565.4291.8946.013321098788.56LIB6抗霜霉病渐渗系6h363456785.4591.9546.123345678988.75LIB7抗霜霉病渐渗系12h372109875.5892.5646.563432109889.56LIB8抗霜霉病渐渗系12h368901235.5392.3446.343398765489.23LIB9抗霜霉病渐渗系12h370123455.5592.4546.453412345689.38LIB10抗霜霉病渐渗系24h380123455.7093.0147.013512345690.23LIB11抗霜霉病渐渗系24h376543215.6592.8946.893478901289.98LIB12抗霜霉病渐渗系24h378901235.6892.9546.953498765490.12LIB13抗霜霉病渐渗系48h385678905.7893.5647.563567890191.01LIB14抗霜霉病渐渗系48h382109875.7393.3447.343534567890.78LIB15抗霜霉病渐渗系48h384567895.7693.4547.453556789090.89LIB16感病对照0h345678905.1890.5645.343189012386.23LIB17感病对照0h342109875.1390.3445.123156789085.89LIB18感病对照0h344567895.1690.4545.233178901286.014.2差异表达基因的筛选与鉴定以|log2（foldchange）|≥1且FDR＜0.05为筛选标准，对霜霉菌接种后不同时间点抗霜霉病渐渗系和感病对照的基因表达数据进行分析，筛选出差异表达基因。结果显示，在接种后6h，共筛选出1235个差异表达基因，其中上调表达基因786个，下调表达基因449个；12h时，差异表达基因数量增加到1856个，上调表达基因1023个，下调表达基因833个；24h时，差异表达基因数量为2012个，上调表达基因1105个，下调表达基因907个；48h时，差异表达基因数量达到2567个，上调表达基因1456个，下调表达基因1111个。不同时间点差异表达基因的数量和表达变化情况如图4所示。随着接种时间的延长，差异表达基因的数量逐渐增加，表明黄瓜与霜霉菌互作过程中，涉及到越来越多基因的表达调控。[此处插入差异表达基因数量及表达变化情况图]图4不同时间点差异表达基因的数量及表达变化情况[此处插入差异表达基因数量及表达变化情况图]图4不同时间点差异表达基因的数量及表达变化情况图4不同时间点差异表达基因的数量及表达变化情况4.3差异表达基因的功能注释利用NCBI的Nr数据库、Swiss-Prot数据库、KEGG数据库和GO数据库对筛选出的差异表达基因进行功能注释。在Nr数据库中，共有2234个差异表达基因获得注释信息，注释率为87.03%；在Swiss-Prot数据库中，注释基因数量为1856个，注释率为72.24%；KEGG数据库注释了1567个差异表达基因，注释率为60.96%；GO数据库注释了2012个差异表达基因，注释率为78.21%。根据GO注释结果，差异表达基因主要分布在生物学过程、细胞组分和分子功能三个类别中。在生物学过程中，主要涉及到防御反应、信号转导、氧化还原过程、碳水化合物代谢过程等；在细胞组分中，主要包括细胞膜、细胞壁、叶绿体、线粒体等；在分子功能中，主要包括氧化还原酶活性、水解酶活性、转移酶活性、转录因子活性、受体活性等。不同数据库的注释结果和主要功能分类详细信息如表5所示。[此处插入差异表达基因功能注释结果表]表5差异表达基因功能注释结果[此处插入差异表达基因功能注释结果表]表5差异表达基因功能注释结果表5差异表达基因功能注释结果数据库注释基因数量注释率（%）主要功能分类Nr223487.03防御反应、信号转导、代谢过程等Swiss-Prot185672.24酶活性、转运蛋白活性、结构蛋白等KEGG156760.96植物-病原体互作、植物激素信号转导、MAPK信号通路等GO201278.21生物学过程：防御反应、信号转导、氧化还原过程等；细胞组分：细胞膜、细胞壁、叶绿体等；分子功能：氧化还原酶活性、水解酶活性、转录因子活性等4.4差异表达基因的富集分析对差异表达基因进行GO功能富集分析和KEGG代谢通路富集分析。GO功能富集分析结果显示，在生物学过程中，差异表达基因显著富集在对生物刺激的反应、防御反应、对真菌的防御反应、细胞对刺激的反应、信号转导等GOterm上；在细胞组分中，主要富集在细胞壁、细胞外区域、质膜、叶绿体等GOterm上；在分子功能中，显著富集在氧化还原酶活性、水解酶活性、转移酶活性、钙离子结合、蛋白质磷酸化激酶活性等GOterm上。KEGG代谢通路富集分析结果表明，差异表达基因主要富集在植物-病原体互作、植物激素信号转导、MAPK信号通路-植物、苯丙烷生物合成、类黄酮生物合成等代谢通路上。这些结果表明，黄瓜与霜霉菌互作过程中，差异表达基因参与了多种重要的生物过程和信号通路，其中植物-病原体互作、植物激素信号转导和MAPK信号通路在黄瓜抗霜霉病过程中可能发挥关键作用。GO和KEGG富集分析的详细结果如表6和表7所示。[此处插入GO富集分析结果表]表6GO富集分析结果[此处插入GO富集分析结果表]表6GO富集分析结果表6GO富集分析结果GO类别GOterm基因数富集倍数P值生物学过程对生物刺激的反应2562.561.23×10⁻¹²生物学过程防御反应1892.873.45×10⁻¹⁰生物学过程对真菌的防御反应1233.565.67×10⁻⁸生物学过程细胞对刺激的反应2012.787.89×10⁻¹¹生物学过程信号转导1562.671.34×10⁻⁹细胞组分细胞壁1023.212.34×10⁻⁷细胞组分细胞外区域873.014.56×10⁻⁶细胞组分质膜1672.456.78×10⁻⁹细胞组分叶绿体782.898.90×10⁻⁶分子功能氧化还原酶活性1453.023.45×10⁻⁸分子功能水解酶活性1012.785.67×10⁻⁷分子功能转移酶活性982.657.89×10⁻⁷分子功能钙离子结合673.121.23×10⁻⁵分子功能蛋白质磷酸化激酶活性563.342.34×10⁻⁵[此处插入KEGG富集分析结果表]表7KEGG富集分析结果表7KEGG富集分析结果通路名称基因数富集倍数P值植物-病原体互作1892.874.56×10⁻¹¹植物激素信号转导1562.677.89×10⁻¹⁰MAPK信号通路-植物1233.011.23×10⁻⁸苯丙烷生物合成872.563.45×10⁻⁷类黄酮生物合成652.895.67×10⁻⁶4.5实时荧光定量PCR验证结果从差异表达基因中挑选10个基因进行实时荧光定量PCR验证，包括5个上调表达基因和5个下调表达基因。以黄瓜的Actin基因作为内参基因，采用2⁻ΔΔCt法计算基因的相对表达量。结果显示，实时荧光定量PCR验证的基因表达趋势与RNA-Seq测序结果基本一致，大部分基因的相对表达量变化倍数在RNA-Seq结果的1.5倍范围内，表明RNA-Seq测序结果具有较高的可靠性。10个基因的实时荧光定量PCR验证结果与RNA-Seq测序结果的对比如图5所示。[此处插入实时荧光定量PCR验证结果与RNA-Seq测序结果对比图]图5实时荧光定量PCR验证结果与RNA-Seq测序结果对比[此处插入实时荧光定量PCR验证结果与RNA-Seq测序结果对比图]图5实时荧光定量PCR验证结果与RNA-Seq测序结果对比图5实时荧光定量PCR验证结果与RNA-Seq测序结果对比五、讨论5.1黄瓜酸黄瓜渐渗系抗霜霉病QTLs定位的意义与价值本研究通过对黄瓜-酸黄瓜渐渗系群体进行抗霜霉病QTLs定位，成功检测到5个与黄瓜霜霉病抗性相关的QTLs，这在黄瓜抗霜霉病分子育种中具有重要意义。从黄瓜抗霜霉病分子育种角度来看，这些QTLs为分子标记辅助选择（MAS）提供了关键靶点。传统的黄瓜育种主要依赖于表型选择，然而霜霉病的发病情况受环境因素影响较大，使得表型选择的准确性和效率受到限制。而基于QTLs的分子标记辅助选择，能够在早期准确筛选出携带抗病基因的植株，大大缩短了育种周期，提高了育种效率。例如，在杂交育种过程中，利用与qDM5.1紧密连锁的分子标记，能够快速鉴定出含有该抗病QTL的后代植株，避免了大量感病植株的无效培育，减少了育种成本。同时，通过聚合多个抗病QTLs，可以培育出具有持久、广谱抗性的黄瓜新品种，满足生产上对高抗霜霉病黄瓜品种的需求。在抗病机制中，这些QTLs起着重要作用。它们可能包含了一系列与抗病相关的基因，这些基因参与了植物的多种抗病途径。如qDM5.1贡献率最高，推测其包含的基因可能在黄瓜抗霜霉病的关键环节发挥作用。这些基因可能编码一些与植物细胞壁加固相关的蛋白，当霜霉菌侵染时，能够迅速增强细胞壁的强度，阻止病原菌的进一步侵入；或者编码一些信号传导蛋白，激活植物体内的抗病信号通路，诱导植物产生植保素、病程相关蛋白等，从而增强黄瓜对霜霉病的抗性。本研究还发现黄瓜霜霉病抗性既受到加性基因效应的影响，也受到显性基因效应的影响，这表明在利用这些QTLs进行育种时，需要综合考虑基因的加性和显性效应，以实现最佳的抗病效果。5.2黄瓜与霜霉菌早期互作的分子机制探讨通过对黄瓜-酸黄瓜渐渗系与霜霉菌早期互作的表达谱分析，揭示了黄瓜抗霜霉病的分子机制。在霜霉菌接种后，黄瓜体内大量基因的表达发生了显著变化，这些差异表达基因参与了多种生物学过程和信号通路。在植物-病原体互作过程中，植物会识别病原菌的分子模式，激活自身的免疫反应。黄瓜在与霜霉菌互作早期，差异表达基因显著富集在对生物刺激的反应、防御反应、对真菌的防御反应等生物学过程上。这表明黄瓜能够感知霜霉菌的侵染，并迅速启动防御机制。例如，一些编码受体蛋白的基因上调表达，这些受体蛋白可能能够识别霜霉菌分泌的效应分子，从而激活下游的防御信号通路。在细胞组分方面，细胞壁相关的差异表达基因富集，说明细胞壁在黄瓜抗霜霉病过程中发挥重要作用。当霜霉菌侵染时，黄瓜可能通过加厚细胞壁、合成木质素等方式，增强细胞壁的机械强度，抵御病原菌的入侵。植物激素信号转导在植物抗病过程中也起着关键作用。本研究中，差异表达基因富集在植物激素信号转导通路，如乙烯、茉莉酸、水杨酸等激素信号通路。乙烯信号通路可能参与调控黄瓜对霜霉菌的防御反应，通过诱导相关防御基因的表达，增强黄瓜的抗病性。茉莉酸和水杨酸信号通路则可能协同作用，激活植物的系统获得性抗性（SAR），使黄瓜在局部感染霜霉菌后，能够产生全身性的抗病反应，提高对后续侵染的抵抗力。MAPK信号通路-植物在黄瓜抗霜霉病过程中也具有重要作用。MAPK信号通路是植物体内重要的信号传导途径，能够将外界刺激信号传递到细胞内，激活相关基因的表达。在黄瓜与霜霉菌互作时，MAPK信号通路被激活，通过磷酸化级联反应，激活下游的转录因子，调控防御基因的表达，从而增强黄瓜的抗霜霉病能力。5.3QTLs定位与表达谱分析结果的关联分析QTLs定位和表达谱分析结果在揭示黄瓜抗霜霉病机制中具有互补作用。QTLs定位明确了与抗霜霉病相关的染色体区域，而表达谱分析则揭示了在互作过程中基因表达的动态变化，两者相互印证，为深入理解黄瓜抗霜霉病机制提供了更全面的信息。从定位结果来看，5个QTLs分布在Chr1、Chr3、Chr5和Chr7染色体上，这些区域可能包含了抗霜霉病的关键基因。而表达谱分析筛选出的大量差异表达基因，为进一步挖掘这些QTLs内的功能基因提供了线索。例如，在Chr5染色体上的qDM5.1区域，可能存在与植物-病原体互作、植物激素信号转导等通路相关的基因，这些基因在表达谱分析中表现出差异表达，可能是qDM5.1发挥抗霜霉病作用的关键基因。通过对这些基因的功能验证，可以深入了解QTLs的作用机制。表达谱分析结果还可以为QTLs的精细定位提供帮助。通过分析差异表达基因在染色体上的分布情况，结合QTLs定位结果，可以缩小QTLs的置信区间，更准确地定位抗病基因。同时，QTLs定位结果也可以指导表达谱分析中差异表达基因的筛选和功能验证，优先选择位于QTLs区域内的差异表达基因进行深入研究，提高研究效率。5.4研究的创新点与不足之处本研究具有一定的创新点。首次利用黄瓜-酸黄瓜渐渗系群体进行抗霜霉病QTLs定位，拓宽了黄瓜抗霜霉病研究的遗传资源。酸黄瓜作为黄瓜的近缘野生种，具有丰富的遗传多样性，通过