基于全基因组分析的烟草赤星病抗性位点挖掘与解析

基于全基因组分析的烟草赤星病抗性位点挖掘与解析一、引言1.1研究背景与意义烟草（NicotianatabacumL.）作为全球范围内重要的经济作物之一，在农业经济领域占据着举足轻重的地位。从种植规模来看，其广泛分布于世界众多国家和地区，为众多烟农提供了主要的经济来源，是许多农业产区的支柱性产业。据相关统计数据显示，全球烟草种植面积达数百万公顷，每年烟草种植相关收入在许多发展中国家的农业经济中占比颇高。在生产加工方面，烟草产业形成了庞大而复杂的产业链，涵盖从烟叶初烤、复烤，到卷烟制造、烟草制品销售等多个环节，创造了大量的就业岗位，为工业生产和商业流通注入活力。同时，烟草产业还为各国政府贡献了可观的财政税收，在国家经济发展中扮演着重要角色。然而，烟草在生长过程中面临着多种病害的威胁，其中赤星病（Alternariaalternata(Fr.)Keissler）是一种极具破坏力的叶部真菌病害，严重制约着烟草产业的可持续发展。赤星病病菌主要以菌丝体和分生孢子在病残体上越冬，待环境条件适宜时，分生孢子通过气流、雨水飞溅等方式传播，侵染烟草叶片。一旦发病，初期在叶片上形成黄褐色圆形小斑点，随后病斑迅速扩大，病斑边缘明显，有黄色晕圈，病斑上会产生明显的同心轮纹，后期病斑中心会出现破裂或穿孔现象，严重影响叶片的光合作用和正常生理功能。在病害流行年份，若防治不及时，发病率可达80%以上，病情严重时，叶片大量枯黄、坏死，甚至整株死亡，导致烟草产量大幅下降，一般减产幅度可达20%-50%，部分重病区减产更为严重。而且，患病烟叶的品质也会受到极大影响，表现为叶片厚度变薄、组织结构疏松、颜色变深、油分减少、香气物质含量降低、杂气和刺激性增加等，使得烟草在卷烟加工过程中的可用性降低，直接影响烟草制品的品质和市场价值，给烟草产业带来巨大的经济损失。发掘烟草全基因组赤星病抗性位点对于烟草育种和产业发展具有关键意义。从烟草育种角度而言，传统的烟草育种方法主要依赖于表型选择，这种方式不仅效率较低，而且准确性有限，育种周期往往较长，难以满足快速应对病害变化的需求。而通过全基因组关联分析（GWAS）、连锁分析等现代生物技术手段，精准定位赤星病抗性位点，能够为分子标记辅助选择（MAS）育种提供坚实的基础。利用与抗性位点紧密连锁的分子标记，育种工作者可以在早期对大量的烟草育种材料进行筛选，准确地识别出携带抗性基因的个体，大大缩短育种周期，提高育种效率，加快培育出具有高抗赤星病能力的烟草新品种的进程。在产业发展方面，种植抗赤星病的烟草品种是最为经济、有效的防治措施。一方面，这可以减少因病害导致的产量损失和品质下降，保障烟草种植户的经济收益，稳定烟草种植产业，维护烟农的生产积极性，进而保障整个烟草产业链的原料供应稳定；另一方面，减少化学农药的使用量，降低农药残留对环境和人体健康的潜在危害，符合绿色农业发展的理念，有助于烟草产业的可持续发展，提升烟草产业在全球市场中的竞争力和可持续发展能力。1.2烟草赤星病概述1.2.1症状与危害烟草赤星病主要危害烟草叶片，在叶片发病初期，会出现黄褐色圆形小斑点，这些斑点直径通常在1-3毫米左右，如同针尖大小，在叶片上零星分布。随着病情的发展，病斑迅速扩大，直径可达10-25毫米，病斑边缘呈现出明显的深褐色，与周围健康组织界限分明，并且伴有黄色晕圈，这是由于病菌分泌的毒素扩散，导致周围组织细胞生理功能紊乱，出现黄化现象。病斑上还会产生明显的同心轮纹，这是因为病菌在侵染过程中，不同阶段对叶片组织的破坏程度不同，形成了类似年轮的结构。在湿度较大的环境下，病斑表面会产生黑色霉层，这是病原菌的分生孢子梗和分生孢子，大量的分生孢子可以借助气流、雨水等进行再次传播，进一步扩大病害的范围。后期，病斑中心会出现破裂或穿孔现象，这是由于病斑部位的细胞坏死、干枯，在外界环境因素的作用下，无法维持正常的组织结构，导致叶片破损。这种破损严重影响了叶片的光合作用和气体交换功能，使得叶片无法正常制造和运输养分，进而影响整个植株的生长发育。除了叶片，烟草的茎部、花梗、蒴果等部位也可能受到赤星病的侵害。茎部发病时，会出现椭圆形或梭形的病斑，颜色较深，呈黑褐色，病斑凹陷，严重时会导致茎部腐烂，影响植株的水分和养分运输，使植株生长受阻，甚至死亡。花梗和蒴果发病后，同样会出现褐色病斑，导致果实发育不良，种子质量下降，影响烟草的繁殖和种质资源。烟草赤星病对烟叶产量和质量产生严重的负面影响。在产量方面，患病的叶片生长受到抑制，叶片面积减小，厚度变薄，单叶重量降低。而且，由于病斑的存在，许多叶片在成熟前就会提前枯黄、脱落，导致可采收的叶片数量减少。在病害流行年份，若防治措施不到位，发病率可达80%以上，病情严重时，烟叶减产幅度可达20%-50%，部分重病区甚至可能绝收，给烟农带来巨大的经济损失。在质量方面，患病烟叶的组织结构变得疏松，颜色变深，失去了正常的金黄色或橘黄色，油分减少，香气物质含量降低，杂气和刺激性增加。这些变化使得患病烟叶在卷烟加工过程中的可用性降低，生产出的卷烟品质下降，口感变差，市场竞争力减弱，直接影响了烟草产业的经济效益。1.2.2病原菌特性烟草赤星病的病原菌主要是链格孢菌（Alternariaalternata(Fr.)Keissler），属于半知菌亚门链格孢属真菌。从形态特征来看，链格孢菌的菌丝体呈灰白色至橄榄色，有分隔，直径在1.5-5微米之间，菌丝在寄主组织内生长蔓延，吸收养分，破坏细胞结构。分生孢子梗单生或丛生，从气孔或表皮细胞间隙伸出，直立或弯曲，不分枝或偶尔分枝，淡褐色至深褐色，具隔膜，大小为（30-100）微米×（4-7）微米。分生孢子呈倒棍棒形或椭圆形，褐色至深褐色，有纵横隔膜，隔膜处常缢缩，顶端有喙状细胞，分生孢子大小差异较大，一般为（15-70）微米×（10-22）微米，分生孢子的形态和大小会受到环境因素和寄主种类的影响而有所变化。在生理特征方面，链格孢菌是一种好气性真菌，适宜在温度25-30℃、相对湿度80%以上的环境中生长繁殖。它对营养物质的需求较为广泛，能够利用多种碳源和氮源，在富含糖类、蛋白质等有机物质的培养基上生长良好。该病原菌具有较强的适应性，能够在不同的土壤类型和气候条件下存活和侵染烟草，但在高温高湿的环境下，其生长速度和致病力会显著增强。链格孢菌的致病机制较为复杂，主要通过分泌多种细胞壁降解酶、毒素和侵染结构来实现对烟草的侵染。在侵染初期，病原菌通过分生孢子在烟草叶片表面萌发，产生芽管，芽管顶端形成附着胞，附着胞分泌粘液，使其牢固地附着在叶片表面。随后，附着胞产生侵入钉，直接穿透叶片表皮细胞或通过气孔侵入叶肉组织。进入寄主细胞后，病原菌分泌一系列细胞壁降解酶，如纤维素酶、果胶酶、半纤维素酶等，这些酶能够分解植物细胞壁的主要成分，破坏细胞的完整性，使病原菌能够进一步侵入和扩散。同时，链格孢菌还会分泌多种毒素，如链格孢酚（AOH）、链格孢醛（AME）等，这些毒素能够抑制植物细胞的正常生理功能，导致细胞坏死、黄化，破坏植物的防御机制，促进病原菌的侵染和致病过程。链格孢菌的侵染循环包括越冬、传播和初侵染、再侵染等环节。病原菌主要以菌丝体和分生孢子在病残体上越冬，成为次年病害发生的初侵染源。在适宜的环境条件下，越冬的分生孢子或新产生的分生孢子通过气流、雨水飞溅、农事操作等方式传播到烟草植株上，落在叶片表面，遇到适宜的温湿度条件即可萌发侵染。一般来说，在烟草生长中后期，田间湿度较大、温度适宜时，病原菌容易侵染发病。发病后，病斑上产生大量的分生孢子，这些分生孢子又可以通过上述传播途径进行再侵染，在一个生长季节内，病害可以多次发生，不断蔓延扩散，导致病情加重。1.2.3发病规律与影响因素烟草赤星病的发病规律与烟草的生长发育阶段密切相关。一般在烟草生长的中后期，尤其是叶片开始成熟时，赤星病容易发生和流行。在团棵期之前，烟草植株生长较为旺盛，自身的抵抗力较强，病原菌侵染相对困难，发病较轻。随着烟草进入旺长期和成熟期，叶片中的淀粉、糖分等营养物质逐渐积累，组织结构也发生变化，使得叶片对赤星病的抵抗力下降，此时病原菌容易侵染发病。而且，烟草生长后期，田间通风透光条件变差，湿度增加，为病原菌的繁殖和传播创造了有利条件，导致病害迅速蔓延。影响烟草赤星病发病的因素众多，其中温度和湿度是最为关键的环境因素。链格孢菌生长和繁殖的适宜温度为25-30℃，相对湿度在80%以上。在这样的温湿度条件下，病原菌的分生孢子萌发率高，侵染速度快，病害容易流行。当温度低于20℃或高于35℃时，病原菌的生长和繁殖会受到抑制，发病较轻。湿度对病害的影响更为显著，高湿度不仅有利于分生孢子的萌发和侵入，还能促进病斑上分生孢子的产生和传播。在多雨、多雾、高湿的天气条件下，赤星病往往会大面积发生，病情严重。相反，在干旱、低湿度的环境中，病害发生受到一定程度的抑制。种植密度对赤星病的发生也有重要影响。种植密度过大时，烟草植株之间通风透光不良，田间湿度增加，形成了一个相对密闭、潮湿的小环境，有利于病原菌的滋生和传播。而且，植株之间相互拥挤，生长空间受限，导致植株生长细弱，抵抗力下降，更容易受到病原菌的侵染。合理的种植密度能够保证田间通风良好，光照充足，降低湿度，减少病原菌的传播机会，从而减轻病害的发生。一般来说，根据烟草品种、土壤肥力和气候条件等因素，确定适宜的种植密度，如行距保持在1-1.2米，株距保持在0.5-0.6米，能够有效预防赤星病的发生。施肥情况也是影响赤星病发病的重要因素。氮肥施用过多会导致烟草植株生长过于旺盛，叶片嫩绿，组织柔软，细胞壁变薄，含糖量降低，含氮量增加，这种生理状态使得植株对赤星病的抵抗力下降，容易被病原菌侵染。而且，过量的氮肥还会导致田间郁闭，湿度增加，为病原菌的繁殖创造有利条件。相反，合理配比氮、磷、钾等肥料，增施有机肥和微量元素肥料，能够增强烟草植株的生长势，提高其抗病能力。一般来说，氮、磷、钾的合理配比为1:1:2-3，同时适量补充锌、硼、锰等微量元素，能够有效降低赤星病的发病率。此外，土壤肥力、地势、前茬作物等因素也会对赤星病的发生产生一定的影响，肥沃的土壤、地势高燥、前茬作物为非茄科植物的烟田，赤星病发病相对较轻。1.3研究现状1.3.1烟草基因组研究进展烟草基因组测序工作取得了重大突破，为烟草遗传研究奠定了坚实基础。2014年，国际烟草基因组测序联盟成功完成了栽培烟草品种“红花大金元”的全基因组测序工作。烟草基因组庞大，约为4.5Gb，包含约74,000个蛋白编码基因，是目前已测序植物基因组中基因数量较多的物种之一。通过对烟草基因组序列的分析，科研人员能够深入了解烟草的遗传结构、基因组成和进化关系，为后续的基因功能研究和分子育种提供了丰富的遗传信息资源。在基因注释方面，大量的基因功能被初步注释，但仍有许多基因的功能有待进一步明确。目前，通过生物信息学分析和实验验证相结合的方法，已经对部分基因进行了功能注释。例如，一些参与烟草生长发育、光合作用、次生代谢等重要生理过程的基因功能已被初步揭示。然而，由于烟草基因组的复杂性和基因功能的多样性，仍有相当比例的基因功能未知，需要进一步开展深入的研究工作，包括基因表达分析、基因敲除、过表达等实验技术，以全面解析这些基因在烟草生长、发育和抗病等过程中的作用机制。烟草基因芯片等技术在基因表达研究中得到广泛应用。基因芯片技术能够同时检测大量基因在不同组织、不同发育阶段以及不同环境条件下的表达水平变化，为研究基因的表达调控机制提供了有力手段。通过基因芯片分析，科研人员发现了许多在烟草赤星病侵染过程中差异表达的基因，这些基因涉及植物的防御反应、信号转导、代谢途径等多个方面。例如，一些编码病程相关蛋白（PR蛋白）、抗氧化酶、转录因子等的基因在赤星病侵染后表达显著上调，表明这些基因可能在烟草抵御赤星病的过程中发挥重要作用。此外，转录组测序（RNA-Seq）技术也逐渐成为研究烟草基因表达的重要工具，它能够更全面、准确地获取基因表达信息，揭示基因在不同条件下的表达模式和调控网络，为深入理解烟草的生物学过程和抗病机制提供了更丰富的数据支持。1.3.2植物抗病基因及QTL定位研究植物抗病基因在植物抵御病原菌侵染的过程中发挥着关键作用，其类型丰富多样，主要包括NBS-LRR类、RLK类、RLP类等。NBS-LRR类抗病基因是植物中最为常见的一类抗病基因，其编码的蛋白含有核苷酸结合位点（NBS）和富含亮氨酸重复序列（LRR）结构域。NBS结构域参与信号转导和ATP结合水解，为抗病反应提供能量；LRR结构域则主要负责识别病原菌的效应子，具有高度的特异性和多样性，能够识别不同病原菌的多种效应子，启动植物的抗病反应。RLK类抗病基因编码的是受体激酶，这类蛋白含有胞外结构域、跨膜结构域和胞内激酶结构域，通过胞外结构域识别病原菌的信号分子，然后将信号传递到胞内，激活下游的抗病信号通路。RLP类抗病基因编码的是富含亮氨酸重复序列的受体蛋白，与RLK类类似，也是通过识别病原菌信号来激活植物的防御反应，但RLP类蛋白不具有胞内激酶结构域，其信号转导可能需要与其他蛋白相互作用来完成。植物抗病基因的作用机制主要基于基因对基因学说。该学说认为，植物中存在抗病基因（R基因），病原菌中存在无毒基因（Avr基因），当植物的R基因编码的蛋白与病原菌的Avr基因编码的效应子相互识别时，会引发植物一系列的防御反应，包括活性氧爆发、植保素合成、细胞壁加厚、病程相关蛋白表达等，从而限制病原菌的生长和繁殖，使植物表现出抗病性。这种识别和防御过程涉及复杂的信号转导网络，多个蛋白激酶、转录因子等参与其中，协同调控植物的抗病反应。数量性状位点（QTL）定位是研究植物数量性状遗传的重要方法，在植物抗病研究中具有广泛应用。QTL是指控制数量性状的基因在基因组中的位置，通过构建遗传连锁图谱，利用分子标记技术对分离群体进行基因型分析，结合表型数据，运用统计学方法可以确定QTL在染色体上的位置和效应。常用的QTL定位方法包括区间作图法、复合区间作图法、多区间作图法等。区间作图法是最早提出的QTL定位方法，它基于两个相邻分子标记，通过检测标记区间内是否存在QTL来确定其位置；复合区间作图法则在区间作图的基础上，考虑了其他标记对目标QTL的影响，能够更准确地定位QTL；多区间作图法进一步拓展，同时考虑多个区间的效应，提高了QTL定位的精度和效率。在植物抗病研究中，QTL定位可以帮助确定与抗病性相关的基因位点，了解抗病性的遗传基础，为抗病育种提供理论依据和分子标记。1.3.3烟草赤星病抗性遗传研究现状烟草赤星病抗性遗传方式较为复杂，研究表明，其抗性既受主效基因控制，也受微效多基因的影响，表现为数量性状遗传。早期的遗传研究通过杂交、回交等传统遗传分析方法，初步确定了烟草赤星病抗性具有一定的遗传规律。例如，一些研究发现，不同烟草品种杂交后代的赤星病抗性表现出连续变异的特点，表明抗性受多个基因的共同作用。而且，抗性还存在明显的加性效应和显性效应，加性效应使得后代的抗性水平可以通过亲本抗性基因的累加而得到提高，显性效应则导致杂合子在抗性表现上可能优于纯合子亲本。在抗性位点定位方面，众多学者通过连锁分析和关联分析等方法，已成功定位了多个与烟草赤星病抗性相关的QTL位点。这些位点分布于烟草的不同染色体上，对赤星病抗性的贡献率各不相同。例如，在烟草第3号染色体上定位到一个主效QTL，其对赤星病抗性的贡献率可达20%以上，该位点的发现为烟草赤星病抗性育种提供了重要的分子标记靶点。此外，在第5、7、9等染色体上也陆续定位到一些与抗性相关的QTL，这些位点的效应相对较小，但它们相互作用，共同影响着烟草的赤星病抗性。在相关基因研究方面，虽然目前已定位的一些与赤星病抗性相关的基因，但对其功能和作用机制的了解还相对有限。部分研究表明，一些编码防御相关蛋白的基因在赤星病抗性中发挥作用，如几丁质酶基因、β-1,3-葡聚糖酶基因等。这些基因在赤星病侵染后表达上调，其编码的蛋白能够降解病原菌细胞壁的主要成分，抑制病原菌的生长和侵染。而且，一些转录因子基因也被发现与赤星病抗性相关，它们通过调控下游防御基因的表达，参与烟草的抗病信号转导过程。然而，整体上对于烟草赤星病抗性基因的功能验证和作用机制解析仍处于初级阶段，需要进一步深入研究，以揭示烟草赤星病抗性的分子机制，为烟草抗病育种提供更坚实的理论基础和技术支持。1.4研究目的与内容1.4.1研究目的本研究旨在通过对烟草种质资源的全面筛选和遗传分析，利用现代生物技术手段，精准发掘烟草全基因组中与赤星病抗性相关的位点。明确这些抗性位点在烟草基因组中的具体位置、遗传效应及作用机制，为深入理解烟草赤星病抗性的遗传基础提供理论依据。基于发掘的抗性位点，开发与之紧密连锁的分子标记，建立高效、准确的分子标记辅助选择体系，为烟草抗赤星病分子育种提供有力的技术支撑，加快培育具有高抗赤星病能力的烟草新品种，有效解决烟草生产中赤星病危害严重的问题，保障烟草产业的可持续发展。1.4.2研究内容烟草种质资源的筛选与鉴定：收集国内外具有代表性的烟草种质资源，包括地方品种、栽培品种、野生种及遗传材料等，建立丰富的种质资源库。在自然发病条件和人工接种条件下，对这些种质资源进行多年、多点的赤星病抗性鉴定，详细记录发病症状、发病率、病情指数等指标，筛选出高抗、中抗和感病的烟草种质材料，为后续研究提供材料基础。遗传群体的构建与表型分析：选择高抗和高感赤星病的烟草品种作为亲本，进行杂交、回交等遗传交配设计，构建F2、BC1等分离群体。对这些遗传群体在不同环境条件下进行赤星病抗性的表型鉴定，统计分析群体中个体的抗性表现，明确赤星病抗性在群体中的遗传变异规律，为抗性位点的定位提供表型数据。全基因组扫描与抗性位点定位：利用简化基因组测序（GBS）、全基因组重测序（WGS）等技术，对构建的遗传群体进行全基因组分子标记开发，构建高密度的遗传连锁图谱。结合群体的表型数据，运用全基因组关联分析（GWAS）、连锁分析等方法，在全基因组范围内扫描与赤星病抗性相关的位点，确定抗性位点在染色体上的位置、遗传效应和贡献率。抗性位点的验证与精细定位：对初步定位的赤星病抗性位点，通过构建近等基因系（NILs）、次级分离群体等方法进行验证，进一步确认抗性位点的真实性和稳定性。利用分子标记加密和重组单株筛选等技术，对关键抗性位点进行精细定位，缩小其在染色体上的物理区间，为后续基因克隆和功能分析奠定基础。抗性位点的功能分析与分子机制研究：在精细定位的基础上，采用图位克隆、转录组测序、基因编辑等技术，分离克隆与赤星病抗性相关的基因。通过基因表达分析、蛋白互作分析、转基因功能验证等手段，深入研究抗性基因的功能和作用机制，明确其在烟草赤星病抗性信号转导途径中的调控网络，揭示烟草抗赤星病的分子机制。二、材料与方法2.1实验材料2.1.1烟草种质资源本研究收集了来自国内外多个地区的共计200份烟草种质资源，其中包含150份栽培品种、30份地方品种以及20份野生种。这些种质资源涵盖了不同的遗传背景和生态类型，具有丰富的遗传多样性，为全面发掘烟草赤星病抗性位点提供了广泛的材料基础。在栽培品种中，囊括了国内广泛种植的“云烟87”“K326”“红花大金元”等品种。“云烟87”具有优质的烟叶品质和较好的适应性，在我国云南、贵州等多个烟草主产区大面积种植；“K326”作为国际上广泛应用的优良品种，具有较强的生长势和较高的产量潜力；“红花大金元”则以其独特的香气风格和较高的经济价值而闻名，是我国重要的烤烟品种之一。地方品种主要来源于一些具有悠久烟草种植历史的地区，如四川什邡的“什邡毛烟”、河南襄城的“襄城金黄柳”等，这些地方品种在长期的种植过程中，适应了当地的生态环境，可能蕴含着独特的抗性基因。野生种包括“绒毛状烟草（Nicotianatomentosiformis）”和“林烟草（Nicotianasylvestris）”等，野生烟草种在长期的自然选择过程中，积累了丰富的遗传变异，对多种病害具有较强的抗性，是挖掘抗性基因的重要资源。在筛选抗病和感病品种时，主要依据前期多年的田间观察和初步的抗性鉴定结果。选择在自然发病条件下连续多年表现出高抗赤星病的“中烟100”作为抗病品种代表。“中烟100”在赤星病高发年份，发病率低于10%，病情指数小于15，叶片病斑数量少且面积小，表现出极强的抗病能力。同时，选取对赤星病高度敏感的“NC89”作为感病品种，“NC89”在相同发病条件下，发病率高达80%以上，病情指数大于70，叶片上布满大量病斑，且病斑迅速扩大融合，导致叶片严重受损，生长发育受到极大抑制。选择这两个品种作为亲本，进行后续遗传群体的构建，能够有效提高抗性位点定位的准确性和可靠性。2.1.2遗传群体构建以高抗赤星病的“中烟100”为母本，高感赤星病的“NC89”为父本进行杂交，获得F1代种子。具体杂交过程如下：在母本“中烟100”开花前，选择发育良好、尚未开放的花蕾，去除雄蕊，以防止自花授粉。然后，采集父本“NC89”当天开放的新鲜花粉，用毛笔轻轻涂抹在母本去雄花蕾的柱头上，完成授粉过程。授粉后，套上硫酸纸袋，防止其他花粉的污染，并做好标记。待授粉后的果实成熟后，收取F1代种子。将F1代种子播种于温室育苗床上，待幼苗长至5-6片真叶时，移栽至田间试验地，按照常规烟草栽培管理措施进行管理。在F1代植株生长至开花期，选择生长健壮的植株，进行自交，获得F2代种子。F2代种子播种和移栽方式同F1代，种植规模为1000株，以保证有足够数量的分离个体用于遗传分析。为进一步验证抗性位点的遗传效应，利用F1代与感病亲本“NC89”进行回交，构建BC1群体。回交过程与杂交过程类似，将F1代植株作为母本，去除雄蕊后，授以“NC89”的花粉，获得BC1代种子。BC1代种植规模为500株，同样在田间进行种植和管理。同时，从F2代群体中随机选取100个单株，自交获得F2:3家系。每个F2:3家系种植20株，种植于田间小区，采用完全随机区组设计，设置3次重复。通过对F2:3家系的表型鉴定，可以更准确地评估个体的基因型，减少环境因素对表型的影响，提高抗性位点定位的精度。在整个遗传群体构建过程中，对每个世代的植株都进行详细的田间记录，包括株高、叶数、生育期等农艺性状，以及赤星病抗性相关的发病时间、病斑数量、病情指数等指标，为后续的遗传分析和抗性位点定位提供全面、准确的数据支持。2.2实验方法2.2.1赤星病抗性鉴定田间自然发病鉴定：在烟草生长季节，选择具有代表性的烟田进行田间自然发病鉴定。烟田的土壤类型、肥力水平和栽培管理措施保持一致，以减少环境因素对病害发生的影响。在烟株团棵期后，定期（每隔7天）进行病情调查，记录每株烟的发病情况。发病症状的观察包括叶片上病斑的出现时间、形状、颜色、大小以及病斑的扩展速度等。发病率的计算方法为发病株数占总调查株数的百分比，即发病率=（发病株数/总调查株数）×100%。病情指数的计算采用分级调查法，将烟株的发病程度分为0-5级，0级为无病斑，1级病斑面积占叶片面积的5%以下，2级病斑面积占叶片面积的6%-15%，3级病斑面积占叶片面积的16%-30%，4级病斑面积占叶片面积的31%-50%，5级病斑面积占叶片面积的50%以上。根据病情分级标准，计算病情指数，公式为：病情指数=∑（各级病株数×各级代表值）/（调查总株数×最高级代表值）×100。通过发病率和病情指数的综合分析，全面评估烟草种质资源在田间自然发病条件下的赤星病抗性水平。人工接种鉴定：人工接种鉴定采用喷雾接种法。首先，从发病烟田采集典型的赤星病病叶，分离纯化病原菌链格孢菌。将病原菌接种到马铃薯葡萄糖琼脂（PDA）培养基上，在28℃恒温培养箱中培养7-10天，待菌落长满培养基平板后，用无菌水冲洗平板表面，收集分生孢子，配制成浓度为1×10⁶个/mL的分生孢子悬浮液，加入适量的吐温-80作为表面活性剂，以增强孢子悬浮液在叶片表面的附着力。在烟草旺长期，选择生长一致、健康无病的烟株进行接种。接种前，将烟株叶片用清水冲洗干净，晾干。然后，使用背负式喷雾器将分生孢子悬浮液均匀地喷洒在烟株叶片的正反两面，以叶片表面布满细密的雾滴为宜。接种后，将烟株置于湿度为85%-95%、温度为25-28℃的人工气候箱中培养，以创造有利于病原菌侵染和发病的环境条件。接种后7-10天开始调查发病情况，同样按照上述病情分级标准计算发病率和病情指数。人工接种鉴定能够排除自然环境中其他因素的干扰，更准确地评估烟草种质资源对赤星病的抗性，为抗性位点的发掘提供可靠的表型数据。2.2.2DNA提取与SNP芯片分析DNA提取：采用改良的CTAB法提取烟草叶片基因组DNA。选取烟草植株顶部第3-4片完全展开的新鲜叶片，用液氮迅速冷冻后研磨成粉末状。将约0.1g叶片粉末转移至1.5mL离心管中，加入700μL预热至65℃的CTAB提取缓冲液（含2%CTAB、1.4MNaCl、100mMTris-HCl，pH8.0、20mMEDTA，pH8.0、0.2%β-巯基乙醇），充分混匀后，置于65℃水浴锅中保温30-60分钟，期间每隔10分钟轻轻颠倒混匀一次，使细胞充分裂解，释放出DNA。保温结束后，冷却至室温，加入等体积的氯仿：异戊醇（24:1）混合液，轻轻颠倒混匀10-15分钟，使蛋白质等杂质充分沉淀到有机相。然后，在12,000rpm条件下离心15分钟，将上清液转移至新的1.5mL离心管中。重复氯仿：异戊醇抽提步骤1-2次，直至中间层无明显杂质。向上清液中加入2/3体积预冷的异丙醇，轻轻颠倒混匀，可见白色絮状DNA沉淀析出。在-20℃冰箱中静置30分钟，使DNA沉淀完全。然后，在12,000rpm条件下离心10分钟，弃去上清液，用70%乙醇洗涤DNA沉淀2-3次，去除残留的盐分和杂质。最后，将DNA沉淀在室温下晾干，加入适量的TE缓冲液（10mMTris-HCl，pH8.0、1mMEDTA，pH8.0）溶解DNA，保存于-20℃冰箱备用。SNP芯片分析：采用Illumina公司的Infinium®TobaccoSNP60K芯片对提取的烟草基因组DNA进行SNP分型检测。该芯片包含约60,000个高质量的SNP标记，覆盖烟草全基因组。芯片实验流程严格按照试剂盒说明书进行操作。首先，将提取的基因组DNA进行片段化处理，然后进行末端修复、加A尾和接头连接等步骤，构建DNA文库。将文库与芯片上的探针进行杂交，经过杂交、洗涤、延伸、染色等一系列反应后，利用IlluminaiScan扫描仪对芯片进行扫描，获取荧光信号数据。利用IlluminaGenomeStudio软件对扫描得到的原始数据进行分析处理，包括背景校正、归一化、基因型calling等步骤，最终得到每个样本在各个SNP位点上的基因型信息。为保证数据质量，对SNP标记进行质量控制，剔除检出率低于90%、最小等位基因频率（MAF）低于0.05以及哈迪-温伯格平衡检验（HWE）P值小于1×10⁻⁶的SNP标记，保留高质量的SNP标记用于后续的遗传图谱构建和关联分析。2.2.3遗传图谱构建利用筛选后的高质量SNP标记，采用JoinMap4.1软件构建遗传图谱。首先，对F2、BC1等遗传群体中的个体进行SNP基因型检测，获得每个个体在各个SNP位点上的基因型数据。将基因型数据导入JoinMap4.1软件中，按照软件默认参数进行数据预处理，包括检查数据完整性、剔除异常数据等。在构建遗传图谱时，采用Kosambi函数计算遗传距离，单位为厘摩（cM）。利用“CPpopulation”选项对F2群体进行分析，“BCpopulation”选项对BC1群体进行分析，通过连锁分析确定SNP标记之间的连锁关系，将连锁的SNP标记划分为不同的连锁群（LG）。连锁群的划分依据是LOD值（似然比检验值），当LOD值大于3.0时，认为两个标记处于同一连锁群。对于每个连锁群，采用最近邻算法确定标记的顺序，逐步构建出完整的遗传图谱。在构建过程中，不断调整标记顺序和遗传距离，以确保遗传图谱的准确性和可靠性。最终得到的遗传图谱包含各个连锁群上SNP标记的顺序和遗传距离信息，为后续的QTL定位分析提供了重要的框架。2.2.4QTL定位分析基于构建的遗传图谱和遗传群体的赤星病抗性表型数据，采用WindowsQTLCartographer2.5软件进行QTL定位分析。在分析过程中，选用复合区间作图法（CIM），该方法能够同时考虑多个标记的效应，有效提高QTL定位的准确性。设置控制标记数为5，以控制背景遗传效应；步长为1cM，即在每个连锁群上每隔1cM进行一次扫描，检测是否存在QTL。通过1000次排列测验（permutationtest）确定LOD阈值，当LOD值大于阈值时，认为该区间存在与赤星病抗性相关的QTL。对检测到的QTL进行命名，命名规则为“q+性状名称+连锁群编号+发现顺序”，例如，“qRB3-1”表示在第3号连锁群上发现的第一个与赤星病抗性（RedStarDiseaseResistance，RB）相关的QTL。同时，计算每个QTL的加性效应、显性效应、贡献率等遗传参数。加性效应反映了等位基因的累加作用，显性效应体现了等位基因之间的相互作用，贡献率则表示该QTL对赤星病抗性表型变异的解释程度。通过QTL定位分析，明确赤星病抗性相关位点在烟草基因组中的位置和遗传效应，为进一步的基因克隆和功能研究提供重要线索。2.2.5关联分析基于自然群体的200份烟草种质资源，进行全基因组关联分析（GWAS），以发掘与赤星病抗性相关的位点。首先，对自然群体进行SNP芯片分型，获得每个种质在全基因组范围内的SNP基因型数据，并进行质量控制，剔除低质量的SNP标记。同时，对自然群体进行多年、多点的赤星病抗性表型鉴定，获取准确的表型数据。采用TASSEL5.0软件进行GWAS分析，选用混合线性模型（MLM），该模型能够有效控制群体结构和个体间的亲缘关系对关联分析结果的影响。在模型中，将群体结构（Q矩阵）和亲缘关系矩阵（K矩阵）作为协变量，以减少假阳性关联的出现。通过计算每个SNP位点与赤星病抗性表型之间的关联显著性，获得P值。为校正多重检验问题，采用Bonferroni校正方法，将显著性阈值设定为α=0.05/标记总数。当某个SNP位点的P值小于校正后的阈值时，认为该位点与赤星病抗性显著关联。对显著关联的SNP位点进行注释，确定其所在的基因区域，并结合生物信息学分析和相关文献资料，预测可能与赤星病抗性相关的候选基因。通过GWAS分析，可以在全基因组范围内快速扫描与赤星病抗性相关的位点，为深入了解烟草赤星病抗性的遗传机制提供新的视角和线索。2.2.6分子标记开发与验证分子标记开发：根据QTL定位和关联分析确定的赤星病抗性位点，利用生物信息学方法开发与之紧密连锁的分子标记。首先，获取抗性位点所在区域的基因组序列，通过与烟草参考基因组进行比对，分析该区域内的序列变异情况，寻找插入/缺失（InDel）、单核苷酸多态性（SNP）等变异位点。针对InDel位点，设计特异性引物，引物设计原则如下：引物长度一般为18-25bp，引物的Tm值在55-65℃之间，上下游引物的Tm值相差不超过5℃；引物3'端避免出现连续3个以上的相同碱基，以防止错配；引物应具有良好的特异性，避免与基因组其他区域产生非特异性扩增。对于SNP位点，采用dCAPS（derivedCleavedAmplifiedPolymorphicSequences）标记技术进行开发。根据SNP位点两侧的序列设计扩增引物，同时引入错配碱基，使扩增产物在SNP位点处形成或消除限制性内切酶识别位点。选择合适的限制性内切酶对扩增产物进行酶切，酶切产物通过聚丙烯酰胺凝胶电泳或毛细管电泳进行检测，根据酶切片段的多态性确定不同基因型。标记验证：选取部分高抗和高感赤星病的烟草品种以及构建的遗传群体中的个体，对开发的分子标记进行验证。利用设计的引物对这些材料进行PCR扩增，反应体系一般为20μL，包括10×PCRbuffer2μL、2.5mMdNTPs1.6μL、上下游引物（10μM）各0.5μL、TaqDNA聚合酶0.2μL、模板DNA50-100ng，用ddH₂O补足至20μL。PCR反应程序为：94℃预变性5分钟；94℃变性30秒，55-65℃退火30秒，72℃延伸30-60秒，共35个循环；72℃延伸10分钟。扩增产物经琼脂糖凝胶电泳或聚丙烯酰胺凝胶电泳分离后，用溴化乙锭（EB）染色或银染法进行检测，观察条带的多态性。将分子标记的基因型数据与赤星病抗性表型数据进行关联分析，验证分子标记与赤星病抗性之间的连锁关系。如果标记基因型与抗性表型之间存在显著关联，表明开发的分子标记具有有效性，可用于烟草抗赤星病分子标记辅助选择育种，为快速筛选抗赤星病烟草材料提供技术支持。三、烟草赤星病抗性遗传分析3.1烟草种质资源的抗性鉴定结果通过田间自然发病鉴定和人工接种鉴定两种方法，对收集的200份烟草种质资源进行了赤星病抗性评估。在田间自然发病条件下，不同种质资源的赤星病发病率和病情指数呈现出显著的差异。发病率最低的种质资源仅为5%，叶片上几乎无明显病斑，生长状况良好；而发病率最高的种质资源可达95%，植株叶片严重受损，大部分叶片布满病斑，枯黄脱落，生长受到极大抑制。病情指数方面，最低值为8.5，表明该种质资源对赤星病具有较强的抵抗力，病斑面积小，扩展速度慢；最高值则达到了85.6，显示该种质资源高度感病，病斑面积大，且迅速蔓延至整个叶片。人工接种鉴定结果与田间自然发病鉴定结果具有较高的一致性，但由于人工接种能够更精准地控制病原菌的侵染条件，发病情况更为集中和明显。在人工接种后10天，高感种质资源的发病率就达到了90%以上，病情指数迅速上升至70以上，叶片出现大量病斑，且病斑融合速度快；而高抗种质资源在接种后15天，发病率仍低于20%，病情指数小于20，叶片仅出现少量分散的小病斑，且病斑扩展不明显。根据发病率和病情指数，将烟草种质资源的赤星病抗性水平划分为高抗、中抗、中感和高感四个等级。其中，高抗种质资源有25份，占比12.5%，这些种质资源在两种鉴定方法下，发病率均低于20%，病情指数小于25，表现出极强的抗赤星病能力；中抗种质资源有60份，占比30%，发病率在20%-40%之间，病情指数在25-50之间，具有较好的抗性；中感种质资源有80份，占比40%，发病率在40%-60%之间，病情指数在50-70之间，对赤星病的抵抗力较弱；高感种质资源有35份，占比17.5%，发病率高于60%，病情指数大于70，极易受到赤星病的侵害。从不同类型的烟草种质资源来看，野生种中高抗和中抗的比例相对较高，分别占野生种总数的30%和40%。这是因为野生烟草在长期的自然选择过程中，为了适应复杂的生态环境和抵御病原菌的侵害，逐渐进化出了多种有效的防御机制，积累了丰富的抗性基因，使其具有较强的抗赤星病能力。地方品种中，中抗和中感的比例较大，分别为35%和45%，这可能与地方品种在特定地区长期种植，对当地的生态条件和病原菌种类具有一定的适应性，但由于缺乏系统的选育和改良，抗性水平相对有限。栽培品种的抗性分布较为均匀，高抗、中抗、中感和高感的比例分别为10%、30%、40%和20%，这是由于栽培品种在选育过程中，虽然注重了产量、品质等性状的改良，但对赤星病抗性的选择力度相对不足，导致不同抗性水平的品种均有一定的数量。对不同抗性等级的烟草种质资源进行聚类分析，结果表明，高抗和中抗的种质资源在遗传距离上相对较近，聚为一类，说明它们可能具有相似的抗性遗传基础；而高感和中感的种质资源聚为另一类，与抗性较强的种质资源在遗传距离上较远。这一结果为后续抗性位点的定位和遗传分析提供了重要的参考依据，有助于缩小研究范围，提高抗性位点发掘的效率。同时，通过对不同抗性水平种质资源的遗传多样性分析发现，高抗种质资源的遗传多样性相对较低，这可能是由于在长期的抗病育种过程中，对少数抗性基因的选择导致遗传背景逐渐狭窄。因此，在今后的育种工作中，需要注重引入更多的遗传资源，拓宽烟草的遗传背景，以进一步提高烟草对赤星病的抗性。3.2遗传群体的性状分离分析对构建的F2和F2:3群体进行了赤星病抗性性状的分离比例分析，以此判断其遗传模式。在F2群体中，共调查了800个单株的赤星病抗性表现。通过对病情指数的统计分析，发现抗性表现呈现出连续的变异，并非简单的孟德尔式分离比例。利用正态分布检验方法对F2群体的病情指数数据进行分析，结果显示其符合正态分布特征，Kolmogorov-Smirnov检验的P值大于0.05，表明赤星病抗性性状在F2群体中表现为数量性状遗传，受多个微效基因的共同控制。进一步对F2群体的抗性数据进行分组统计，以病情指数30为界限，将单株分为抗病（病情指数≤30）和感病（病情指数＞30）两组。经统计，抗病单株数量为235株，感病单株数量为565株，抗病与感病单株的比例约为1:2.4，不符合经典的孟德尔单基因分离比例（3:1），也进一步证实了赤星病抗性受多基因控制。为了更准确地分析多基因的遗传效应，利用主基因+多基因混合遗传模型对F2群体的抗性数据进行分析。结果表明，该群体中赤星病抗性受2对加性-显性-上位性主基因+加性-显性多基因控制，主基因遗传率为58.6%，多基因遗传率为26.4%，环境因素对表型变异的影响率为15%。这说明在F2群体中，主基因在赤星病抗性遗传中起主要作用，但多基因和环境因素也不可忽视。在F2:3家系中，同样对赤星病抗性进行了鉴定和分析。每个F2:3家系种植20株，共调查了100个家系。家系内个体的抗性表现存在一定的差异，但家系间的抗性差异更为明显。通过对家系平均病情指数的分析，发现其也呈现出连续的变异分布。利用广义遗传力公式计算F2:3家系的广义遗传力，结果为72.5%。这表明在F2:3家系中，遗传因素对赤星病抗性的影响较大，环境因素的影响相对较小。通过对F2:3家系抗性数据的相关性分析，发现家系平均病情指数与F2单株病情指数之间存在显著的正相关关系，相关系数为0.78，这进一步验证了F2群体中赤星病抗性的遗传规律，也说明F2:3家系能够较好地反映F2单株的遗传信息，可用于后续的遗传分析和抗性位点定位。综合F2和F2:3群体的性状分离分析结果，可以得出烟草赤星病抗性为数量性状遗传，受多个主基因和多基因的共同控制，且遗传因素在抗性表现中起主导作用，但环境因素也会对表型产生一定的影响。这些结果为后续利用遗传群体进行赤星病抗性位点的定位和遗传机制研究提供了重要的理论基础。3.3赤星病抗性的遗传效应分析利用主基因+多基因混合遗传模型对烟草赤星病抗性进行遗传效应分析，结果显示，在烟草赤星病抗性遗传中，主基因和多基因均发挥着重要作用。主基因的加性效应显著，表明来自不同亲本的抗性等位基因在累加时，能够有效增强烟草对赤星病的抗性水平。以“中烟100”和“NC89”构建的遗传群体为例，携带来自“中烟100”抗性主基因的个体，在赤星病抗性表现上明显优于未携带该基因的个体，病情指数显著降低。而且，主基因还存在一定的显性效应，杂合基因型在某些情况下能够表现出比纯合显性或纯合隐性更好的抗性，这为杂种优势在烟草抗赤星病育种中的应用提供了理论依据。多基因的效应虽然相对较小，但它们的累积作用也不可忽视。多个微效多基因共同作用，能够进一步丰富烟草赤星病抗性的遗传变异。在田间试验中可以观察到，一些烟草个体虽然没有携带明显的主效抗性基因，但由于积累了多个微效多基因，依然表现出一定程度的抗性。而且，多基因之间还存在复杂的互作关系，它们通过相互协调，共同参与烟草对赤星病的防御反应。这些微效多基因可能分别参与了植物的不同防御途径，如细胞壁加厚、活性氧清除、植保素合成等，它们的协同作用使得烟草能够更全面地抵御赤星病的侵染。环境因素对烟草赤星病抗性也有一定的影响，其效应值约占总变异的15%。在不同的种植环境下，同一基因型的烟草赤星病抗性表现可能会有所差异。在高温高湿的环境条件下，烟草赤星病发病更为严重，即使是抗性较强的基因型，其病情指数也会有所上升；而在相对干燥、温度适宜的环境中，发病程度则相对较轻。这表明在烟草抗赤星病育种和生产过程中，不仅要关注遗传因素，还需要合理调控环境条件，以充分发挥烟草的抗性潜力。例如，通过合理密植、科学施肥、改善田间通风透光条件等措施，可以降低环境因素对赤星病发病的影响，提高烟草的抗病能力。四、全基因组赤星病抗性位点发掘4.1SNP标记的筛选与遗传图谱构建通过IlluminaInfinium®TobaccoSNP60K芯片对烟草种质资源及遗传群体进行SNP分型检测，共获得60,000个初始SNP标记。经过严格的质量控制，剔除了检出率低于90%的标记，此类标记由于在样本中检测不到或检测不稳定，会影响后续分析的准确性，共剔除约5000个；剔除最小等位基因频率（MAF）低于0.05的标记，这些标记的变异频率过低，在群体中分布不均匀，对遗传分析贡献较小，约有8000个被剔除；同时，剔除哈迪-温伯格平衡检验（HWE）P值小于1×10⁻⁶的标记，这类标记可能受到选择、突变、遗传漂变等因素的影响，不符合遗传平衡定律，不利于准确分析遗传关系，约3000个被剔除。最终，筛选得到44,000个高质量的SNP标记，用于后续的遗传图谱构建和分析。利用JoinMap4.1软件，基于筛选后的SNP标记，成功构建了烟草遗传图谱。该图谱共包含24个连锁群（LG），与烟草的染色体数目一致。图谱总长度为1800cM，平均遗传距离为0.41cM，标记分布较为均匀。各连锁群长度在40-150cM之间，连锁群上的SNP标记数量在1000-3000个不等。例如，LG1连锁群长度为120cM，包含2500个SNP标记，平均每0.048cM就有一个标记；LG12连锁群长度为60cM，含有1500个标记，平均遗传距离为0.04cM。为评估遗传图谱的质量，进行了多项指标的检测。首先，通过计算标记顺序的一致性，发现图谱中98%以上的标记顺序与烟草参考基因组的物理顺序一致，表明图谱的构建准确可靠，能够真实反映染色体上基因的排列顺序。而且，利用“鬼影标记”（phantommarker）检测法，未发现明显的图谱压缩或膨胀现象，进一步验证了图谱的准确性。同时，通过对遗传图谱与物理图谱的共线性分析，发现两者具有高度的一致性，各连锁群上的标记在物理图谱上的位置与遗传图谱中的顺序相对应，这为后续基于遗传图谱进行基因定位和克隆提供了有力的保障。高质量的遗传图谱为烟草赤星病抗性位点的精准定位奠定了坚实基础，使得在全基因组范围内扫描抗性位点成为可能，能够更准确地确定抗性位点在染色体上的位置，为深入研究赤星病抗性的遗传机制提供了重要的框架。4.2QTL定位结果4.2.1抗赤星病QTL的分布通过WindowsQTLCartographer2.5软件，运用复合区间作图法（CIM）对烟草赤星病抗性进行QTL定位分析，在烟草全基因组范围内共检测到15个与赤星病抗性相关的QTL。这些QTL分布于10条不同的连锁群上，呈现出较为广泛的分布特征。在LG3连锁群上，定位到3个QTL，分别为qRB3-1、qRB3-2和qRB3-3，它们在连锁群上的位置依次为45-50cM、70-75cM和90-95cM。其中，qRB3-1位于连锁群的中部偏左位置，与两侧标记的遗传距离分别为2.5cM和3cM；qRB3-2处于连锁群的中后部，与相邻标记的距离较为均匀，均为2.5cM；qRB3-3则靠近连锁群的末端，距离最近的标记为2cM。在LG5连锁群上，检测到2个QTL，qRB5-1和qRB5-2，位置分别在30-35cM和60-65cM，qRB5-1位于连锁群的前部，与周边标记的距离在2-3cM之间，qRB5-2处于连锁群的中部，与相邻标记的距离相对稳定，均为2.5cM。LG7、LG9、LG11、LG13、LG15、LG17、LG19和LG21连锁群上各有1个QTL，它们在各自连锁群上的位置和与标记的距离也各不相同。例如，qRB7-1位于LG7连锁群的55-60cM处，与两侧标记的距离分别为2.8cM和2.2cM；qRB9-1在LG9连锁群的40-45cM位置，与相邻标记的距离均为2.5cM。不同连锁群上QTL的分布密度存在一定差异。LG3连锁群上QTL分布相对密集，在100cM左右的区间内就定位到3个QTL，这表明该连锁群上可能存在较多与赤星病抗性相关的基因或调控区域，对烟草赤星病抗性的遗传贡献较大。而LG7、LG9等连锁群上仅检测到1个QTL，分布相对稀疏，可能这些连锁群上与赤星病抗性直接相关的基因数量较少，或者其抗性相关基因的效应相对较小，在当前的检测条件下较难被发现。这种QTL在不同连锁群上的分布差异，反映了烟草赤星病抗性遗传机制的复杂性，不同连锁群上的基因或基因区域在调控赤星病抗性过程中发挥着不同程度的作用，为深入研究烟草赤星病抗性的遗传网络提供了重要线索。4.2.2主效QTL的确定与效应分析在检测到的15个QTL中，通过对LOD值、贡献率等指标的综合评估，确定了3个主效QTL，分别为qRB3-1、qRB5-1和qRB7-1。qRB3-1位于LG3连锁群的45-50cM区间，其LOD值高达5.8，贡献率为25.6%。这表明该QTL在烟草赤星病抗性遗传中具有重要作用，能够解释25.6%的表型变异。从加性效应来看，qRB3-1的加性效应值为-4.5，说明来自抗病亲本“中烟100”的等位基因具有降低病情指数的作用，每增加一个该等位基因，病情指数平均降低4.5。在不同环境下，对qRB3-1进行稳定性分析，发现在3个不同的试验环境中，该QTL均能被检测到，且LOD值和贡献率的波动范围较小，LOD值在5.5-6.0之间，贡献率在24.5%-26.5%之间，表现出较好的稳定性。qRB5-1位于LG5连锁群的30-35cM区间，LOD值为5.2，贡献率为22.3%。其加性效应值为-3.8，即来自抗病亲本的等位基因可使病情指数平均降低3.8。在不同环境的稳定性分析中，在2个环境中能够稳定检测到，在另一个环境中虽然LOD值略有下降，但仍高于阈值，贡献率在不同环境中的波动范围为20.5%-23.5%，具有一定的稳定性。qRB7-1位于LG7连锁群的55-60cM区间，LOD值为4.8，贡献率为20.1%。加性效应值为-3.2，表明其抗病等位基因可降低病情指数3.2。在不同环境下，该QTL在2个环境中检测结果较为稳定，在第三个环境中，由于环境因素的影响，LOD值有所降低，但仍能被检测到，贡献率波动范围为18.5%-21.5%。这3个主效QTL的确定，为烟草抗赤星病分子育种提供了关键的靶点。在实际育种过程中，可以利用与这些主效QTL紧密连锁的分子标记，对烟草育种材料进行早期筛选，准确地选择携带抗病等位基因的个体，提高育种效率。而且，对它们在不同环境下稳定性的分析，有助于评估其在不同生态条件下的应用潜力，为培育适应不同环境的抗赤星病烟草品种提供了重要依据，进一步推动了烟草抗赤星病育种工作的开展。4.3关联分析结果利用TASSEL5.0软件对自然群体的200份烟草种质资源进行全基因组关联分析（GWAS），共检测到28个与赤星病抗性显著关联的SNP位点（P<0.05/标记总数）。这些SNP位点分布于烟草的12条染色体上，在染色体上呈现出不均匀分布的特点。在Chr2染色体上，发现了3个显著关联的SNP位点，分别为SNP2-1、SNP2-2和SNP2-3，它们的物理位置依次为5,000,000-5,000,050bp、6,500,000-6,500,050bp和8,000,000-8,000,050bp。其中，SNP2-1位于基因Nt001的内含子区域，该基因编码一种未知功能的蛋白，但其表达模式在抗病和感病材料中存在显著差异，推测可能与赤星病抗性相关；SNP2-2位于基因Nt002的上游调控区，距离基因起始密码子约1000bp，Nt002基因参与植物激素信号转导途径，可能通过调控激素信号来影响烟草对赤星病的抗性；SNP2-3位于基因Nt003的外显子区域，导致编码蛋白的第56个氨基酸发生改变，Nt003基因编码一种病程相关蛋白，这种氨基酸的改变可能影响其蛋白结构和功能，进而影响烟草的抗病性。在Chr6染色体上，有4个显著关联的SNP位点，分别为SNP6-1、SNP6-2、SNP6-3和SNP6-4，分布在12,000,000-12,000,050bp、13,500,000-13,500,050bp、15,000,000-15,000,050bp和16,500,000-16,500,050bp区间。SNP6-1位于基因Nt004的启动子区域，可能通过影响基因的转录起始来调控基因表达，Nt004基因与植物细胞壁合成相关，细胞壁作为植物抵御病原菌入侵的第一道防线，其合成相关基因的表达变化可能影响烟草对赤星病的抗性；SNP6-2位于基因间区，但与一个已知的转录因子结合位点距离较近，推测可能通过影响转录因子与DNA的结合来调控附近基因的表达，进而影响赤星病抗性；SNP6-3和SNP6-4分别位于基因Nt005和Nt006的外显子区域，导致编码蛋白的氨基酸改变，Nt005基因参与活性氧代谢途径，Nt006基因编码一种信号转导蛋白，它们的功能变化可能在烟草抗赤星病过程中发挥作用。对这些显著关联的SNP位点进行进一步分析，发现部分位点之间存在紧密的连锁不平衡（LD）关系。例如，在Chr9染色体上的SNP9-1和SNP9-2，它们之间的距离仅为5000bp，LD值（r²）高达0.85，表明这两个位点在群体中倾向于一起遗传，可能受到相同的选择压力或处于同一连锁块内。而且，通过对SNP位点与赤星病抗性表型的相关性分析，发现不同SNP位点对抗性表型的影响程度不同。一些SNP位点的等位基因变异能够显著降低病情指数，如SNP3-1的抗性等位基因可使病情指数平均降低10.5，对赤星病抗性具有积极的贡献；而另一些位点的变异则与感病性相关，如SNP11-1的某一等位基因会使病情指数升高8.2，增加了烟草对赤星病的敏感性。综合关联分析结果，这些与赤星病抗性显著关联的SNP位点及其所在的基因区域，为深入研究烟草赤星病抗性的遗传机制提供了重要线索。通过进一步的功能验证和基因克隆工作，可以揭示这些位点和基因在烟草抗赤星病过程中的具体作用，为烟草抗赤星病分子育种提供更多的理论依据和分子标记资源。4.4QTL与关联分析结果的整合对比QTL定位和关联分析结果，发现二者在一定程度上具有一致性，同时也表现出互补性。在一致性方面，部分在QTL定位中检测到的抗性区域，在关联分析中也发现了与之相关的SNP位点。例如，在QTL定位中，在LG3连锁群的45-50cM区间检测到主效QTLqRB3-1；在关联分析中，在该区间对应的物理位置附近，发现了多个与赤星病抗性显著关联的SNP位点，如SNP3-2、SNP3-3等。这些SNP位点可能位于qRB3-1所包含的抗性基因内部或其调控区域，进一步验证了该区域在烟草赤星病抗性中的重要作用，为抗性基因的精细定位和克隆提供了双重证据支持。从互补性角度来看，QTL定位基于遗传群体，能够准确地确定抗性位点在染色体上的遗传位置和效应，但由于遗传群体的构建和分析方法的限制，检测到的QTL区间通常较大，分辨率相对较低，难以精确确定具体的抗性基因。而关联分析基于自然群体，利用群体中丰富的遗传变异信息，能够在全基因组范围内快速扫描与性状相关的位点，具有较高的分辨率，可精确到单个SNP位点，但关联分析容易受到群体结构、连锁不平衡等因素的影响，可能产生假阳性结果。因此，二者相互补充，QTL定位为关联分析提供了重要的候选区域，缩小了关联分析的搜索范围，减少了假阳性的干扰；关联分析则为QTL定位提供了更精细的遗传信息，有助于在QTL区间内确定具体的功能变异位点和候选基因。通过整合QTL定位和关联分析结果，构建了一个更为全面的烟草赤星病抗性遗传图谱。将QTL区间与关联分析中的显著SNP位点进行整合，明确了不同抗性区域内的关键遗传变异。利用生物信息学工具，对这些整合后的遗传信息进行分析，预测了多个可能与赤星病抗性相关的候选基因。这些候选基因涉及植物的多个生理过程，如防御反应、信号转导、代谢途径等。例如，在整合区域内发现了一个编码几丁质酶的基因，几丁质酶能够降解病原菌细胞壁的主要成分几丁质，在植物抵御真菌病害过程中发挥重要作用；还发现了一个参与MAPK信号通路的基因，该信号通路在植物响应病原菌侵染的信号转导过程中起关键作用，可能通过激活下游防御基因的表达来增强烟草对赤星病的抗性。这种整合分析方法为深入研究烟草赤星病抗性的分子机制提供了更全面、准确的信息，有助于加速抗性基因的克隆和功能验证，推动烟草抗赤星病分子育种工作的深入开展。五、抗性位点的验证与候选基因分析5.1分子标记的开发与验证基于QTL定位和关联分析确定的赤星病抗性位点，利用生物信息学方法开发了与之紧密连锁的分子标记。针对抗性位点所在区域的插入/缺失（InDel）和单核苷酸多态性（SNP）变异，设计了特异性引物。共设计了20对InDel标记引物和15对基于SNP位点开发的dCAPS标记引物。对开发的分子标记在不同烟草材料中的多态性进行检测。选取了10份高抗赤星病的烟草品种、10份高感品种以及F2群体中的50个单株进行验证。以InDel标记引物IDL-5为例，该引物对应于LG5连锁群上qRB5-1抗性位点附近的一个InDel变异。在高抗品种“中烟100”和“云烟97”中，扩增得到的片段大小为200bp；而在高感品种“NC89”和“RG17”中，扩增片段大小为220bp，表现出明显的多态性。在F2群体的50个单株中，扩增结果显示，片段大小为200bp的单株有28株，片段大小为220bp的单株有22株，其分离比例符合孟德尔遗传定律，初步表明该标记与赤星病抗性存在连锁关系。对于dCAPS标记，以针对SNP位点SNP3-2开发的dCAPS-3标记为例，该位点位于qRB3-1抗性位点区域。通过引入错配碱基，使扩增产物在SNP3-2位点处形成限制性内切酶HindⅢ的识别位点。在高抗品种中，扩增产物经HindⅢ酶切后，产生100bp和50bp两个片段；而在高感品种中，由于SNP位点的差异，扩增产物不能被酶切，仍为150bp的完整片段。在F2群体的检测中，酶切结果与抗性表型的关联分析显示，携带100bp和50bp酶切片段的单株中，80%表现为抗病；而携带150bp未酶切片段的单株中，75%表现为感病，进一步验证了该dCAPS标记与赤星病抗性的紧密连锁关系。通过对不同烟草材料的多态性检测和与抗性表型的关联分析，共验证了12个InDel标记和10个dCAPS标记与赤星病抗性具有显著连锁关系。这些分子标记可用于烟草抗赤星病分子标记辅助选择育种，在早期对大量烟草育种材料进行筛选，准确地识别出携带抗性基因的个体，提高育种效率，加快抗赤星病烟草新品种的培育进程。5.2主效QTL区域的精细定位为了进一步明确主效QTL的具体位置和功能，利用近等基因系（NILs）和次级分离群体对主效QTLqRB3-1、qRB5-1和qRB7-1进行了精细定位。以qRB3-1为例，通过多代回交和自交，构建了以感病亲本“NC89”为遗传背景，携带抗病亲本“中烟100”qRB3-1区域的近等基因系NIL-qRB3-1。在回交过程中，每代利用与qRB3-1紧密连锁的分子标记进行前景选择，确保目标QTL区域的导入；同时利用全基因组SNP芯片进行背景选择，加速遗传背景向“NC89”的回复。经过5代回交和2代自交，成功获得了遗传背景与“NC89”高度相似，仅在qRB3-1区域存在差异的近等基因系。对NIL-qRB3-1与“NC89”进行赤星病抗性鉴定，结果显示NIL-qRB3-1的病情指数显著低于“NC89”，进一步证实了qRB3-1对赤星病抗性的重要作用。利用NIL-qRB3-1与“NC89”构建次级分离群体，在该群体中，对qRB3-1区域进行分子标记加密。在原有标记的基础上，根据该区域的基因组序列，新开发了10个多态性分子标记，使标记密度提高到平均每100kb一个标记。通过对次级分离群体中大量单株的基因型和表型分析，将qRB3-1的物理区间从初始定位的500kb缩小到150kb，大大提高了定位的精度。同样的方法应用于qRB5-1和qRB7-1的精细定位。对于qRB5-1，构建的近等基因系NIL-qRB5-1在抗性鉴定中表现出明显的抗性优势，病情指数比“NC89”降低了30%以上。通过在次级分离群体中加密分子标记，将qRB5-1的物理区间从350kb缩小到100kb。qRB7-1的精细定位过程中，近等基因系NIL-qRB7-1的抗性显著优于“NC89”，利用次级分离群体将其物理区间从400kb缩小到120kb。经过精细定位，确定了每个主效QTL区域内的候选基因。在qRB3-1的150kb区间内，预测包含15个候选基因；qRB5-1的100kb区间内有10个候选基因；qRB7-1的120kb区间内有12个候选基因。这些候选基因的功能涉及多个方面，包括防御反应相关基因、信号转导基因、转录因子基因等。在qRB3-1区域内的候选基因NtDEF1，编码一种防御素蛋白，防御素在植物抵御病原菌侵染过程中具有抗菌活性，可能通过直接作用于赤星病病原菌链格孢菌，抑制其生长和繁殖；qRB5-1区域的候选基因NtSTK1，是一个丝氨酸/苏氨酸蛋白激酶基因，参与植物的信号转导过程，可能在感知赤星病病原菌信号后，激活下游的防御反应途径；qRB7-1区域的候选基因NtTF1，编码一个转录因子，转录因子能够结合到其他基因的启动子区域，调控基因的表达，NtTF1可能通过调控一系列防御基因的表达，来增强烟草对赤星病的抗性。这些候选基因的确定为后续深入研究烟草赤星病抗性的分子机制提供了重要的靶点。5.3候选基因预测与功能分析在主效QTL区域精细定位的基础上，利用生物信息学工具对候选基因进行预测。通过与公共数据库（如NCBI、TAIR等）中的基因序列进行比对，结合基因结构预测软件（如GeneMark、Augustus等），对qRB3-1、qRB5-1和qRB7-1等主效QTL区域内的候选基因进行全面分析。在qRB3-1的150kb区间内，共预测到15个候选基因，其中5个基因编码的蛋白含有与植物防御反应相关的结构域，如NBS-LRR结构域、LRR结构域等；3个基因编码的蛋白与信号转导相关，包括蛋白激酶、磷酸酶等；2个基因编码转录因子，可能参与调控下游防御基因的表达。对于预测到的候选基因，采用多种实验方法进行功能分析。首先，利用实时荧光定量PCR（qRT-PCR）技术分析候选基因在抗病和感病烟草材料中的表达模式。以qRB3-1区域内的候选基因NtDEF1为例，提取接种赤星病病原菌前后抗病品种“中烟100”和感病品种“NC89”叶片的总RNA，反转录成cDNA后进行qRT-PCR分析。结果显示，在接种病原菌后，NtDEF1在“中烟100”中的表达量迅速上调，在6小时时表达量开始显著增加，12小时时达到峰值，为接种前的5倍；而在“NC89”中，NtDEF1的表达量在接种前后无明显变化。这表明NtDEF1基因可能在烟草抵御赤星病的过程中发挥重要作用，其表达受病原菌诱导，且在抗病品种中的表达调控更为显著。利用基因编辑技术（如CRISPR/Cas9）对候选基因进行功能验证。针对qRB5-1区域内的候选基因NtSTK1，设计特异性的sgRNA，构建CRISPR/Cas9基因编辑载体，通过农杆菌介导的转化方法导入烟草细胞中，获得NtSTK1基因编辑突变体。对突变体进行赤星病抗性鉴定，结果显示，与野生型相比，NtSTK1突变体对赤星病的敏感性显著增加，病情指数提高了30%以上，病斑面积明显增大，扩展速度加快。这表明NtSTK1基因的缺失导致烟草对赤星病的抗性下降，进一步证实了该基因在烟草抗赤星病过程中的关键作用，可能通过参与信号转导途径，调控植物的防御反应。通过酵母双杂交实验分析候选基因编码蛋白之间以及与已知抗病相关蛋白之间的相互作用关系。以qRB7-1区域内的候选基因NtTF1为例，将NtTF1基因克隆到酵母双杂交诱饵载体上，构建重组质粒，转化酵母细胞。同时，构建烟草cDNA文库，将文库质粒转化到含有诱饵质粒的酵母细胞中，进行酵母双杂交筛选。结果筛选到3个与NtTF1相互作用的蛋白，其中一个蛋白为已知的病程相关蛋白NtPR1。进一步的验证实验表明，NtTF1与NtPR1在酵母细胞和烟草原生质体中均能发生相互作用。这表明NtTF1可能通过与NtPR1相互作用，调控病程相关蛋白的表达，从而参与烟草对赤星病的防御反应，揭示了NtTF1在烟草抗赤星病信号转导网络中的作用机制。综合生物信息学分析和实验验证结果，初步确定了多个与烟草赤星病抗性密切相关的候选基因及其功能。这些结果为深入理解烟草赤星病抗性的分子机制提供了重要依据，也为烟草抗