基于分子标记的玉米种质特征剖析与耐旱性状关联解析

基于分子标记的玉米种质特征剖析与耐旱性状关联解析一、引言1.1研究背景玉米（ZeamaysL.）作为全球最重要的谷类作物之一，在人类的生产生活中占据着举足轻重的地位。从粮食角度看，它是众多发展中国家的主食来源，为大量人口提供了基本的食物卡路里。据统计，玉米与大米、小麦一起，为94个发展中国家的45亿多人提供了至少30%的食物卡路里，其中包含9亿贫困消费者，玉米对他们而言是首选的主食。在饲料领域，玉米堪称“饲料之王”，是动物饲料的主要成分。随着全球畜牧业的蓬勃发展，对玉米的需求持续攀升。以中国和印度为例，经济的快速增长使得人们对牛奶、鸡蛋和肉类的消费能力增强，进而推动了对玉米作为饲料原料的需求。在工业方面，玉米是生物乙醇燃料的主要原料，也是众多化工产品的重要基础，在生物制药、食品加工、造纸等行业有着广泛应用。从种植范围来看，玉米在全球范围内广泛种植，从热带到温带地区都有其身影，种植面积接近2亿公顷。美国、中国、巴西等国家是主要的玉米生产国。美国凭借其优越的自然条件和先进的农业技术，长期位居全球玉米产量首位，其玉米产量占全球总产量的30%以上。中国作为第二大玉米生产国，玉米种植面积和产量也在不断增长，在保障国内粮食安全和饲料供应方面发挥着关键作用。然而，全球气候变化导致极端气候事件日益频繁，干旱成为制约玉米生产的主要非生物胁迫因素之一。干旱对玉米的影响贯穿其整个生育期。在苗期，干旱会抑制玉米的生长速率，使植株生长缓慢，叶片发黄，茎秆细小，发育期显著延迟。即使后期雨水调和，也难以形成粗壮茎秆孕育大穗。例如，在一些干旱地区，玉米苗期遭遇干旱，植株矮小瘦弱，光合作用受到抑制，影响了后期的生长发育。在拔节期，干旱会促进玉米早熟，导致发育期明显缩短。同时，干旱还会使玉米株高受到明显抑制。在抽雄吐丝期，高温干旱可能导致雄穗和雌穗发育不同步，影响授粉结实。灌浆期的干旱会使籽粒灌浆不充分，千粒重降低，果穗变小、籽粒不饱满、秃尖长，严重影响玉米的产量。此外，干旱还可能影响玉米的品质，如蛋白质、淀粉等含量下降，影响玉米的营养价值和加工品质。挖掘玉米种质资源中的耐旱基因，对玉米种质进行改良，培育耐旱性强的玉米品种，是应对干旱胁迫、保障玉米产量和品质的关键措施。通过深入研究玉米的耐旱机制，找到与耐旱相关的基因和分子标记，能够为玉米耐旱育种提供有力的理论支持和技术手段，从而提高玉米在干旱环境下的适应性和产量稳定性，对于保障全球粮食安全和农业可持续发展具有重要意义。1.2研究目的与意义本研究旨在通过对不同类型玉米种质的分子特征进行分析，明确其遗传多样性和遗传结构。同时，对耐旱相关性状进行深入探究，运用连锁-连锁不平衡联合作图的方法，定位与耐旱性状相关的基因位点，解析其遗传机制。本研究具有重要的理论意义。通过对不同类型玉米种质分子特征的分析，可以深入了解玉米种质的遗传多样性和遗传结构，为玉米种质资源的保护、利用和创新提供理论基础。对耐旱相关性状的连锁-连锁不平衡联合作图，有助于揭示玉米耐旱性的遗传机制，丰富植物耐旱性的理论研究。在实践应用方面，本研究成果将为玉米耐旱育种提供重要的理论依据和技术支持。通过筛选与耐旱性状紧密关联的分子标记，可应用于分子标记辅助选择育种，提高育种效率，加快耐旱玉米品种的选育进程。这对于应对干旱胁迫、保障玉米产量和品质、促进农业可持续发展具有重要的现实意义。1.3国内外研究现状在玉米种质资源研究方面，国内外学者做了大量工作。北京市农林科学院玉米所赵久然/王凤格团队广泛收集玉米自交系种质共5822份，利用自主研发的Maize6H-60K基因芯片进行基因型鉴定，群体结构分析表明，5822份自交系可划分为7大类群，与我国现代玉米育种中广泛应用的核心育种群体相一致。四川农业大学卢艳丽教授团队证明了大刍草对玉米遗传改良的有利价值，挖掘了以穗位高为代表的理想株型性状相关QTL及其被注释为乙烯响应转录因子的关键候选基因，为今后充分利用大刍草开展玉米高产、抗逆育种奠定基础。在玉米耐旱性研究领域，众多研究聚焦于耐旱相关基因的挖掘和作用机制的解析。甘肃省农业科学院作物研究所联合兰州大学生命科学学院的研究团队揭示了转录因子DSD1/ZmICEb在调控玉米气孔发育和耐旱性中的关键作用，为抗旱玉米品种的培育开辟了新路径。还有团队通过Mu转座子插入和CRISPR-Cas9技术，获得wi5、wi6、nst3干旱胁迫敏感突变体，验证了候选基因ZmXYL、ZmGXM1、ZmNST3的生物学功能，首次揭示其通过调控玉米细胞壁合成并激发抗氧化酶系统关键酶的活性来响应干旱胁迫的作用机制。连锁-连锁不平衡联合作图在玉米研究中的应用也逐渐受到关注。该方法能够整合连锁分析和连锁不平衡分析的优势，更精确地定位与性状相关的基因位点。已有研究利用此方法解析玉米重要性状的等位遗传变异，但在玉米耐旱相关性状的研究中，应用还相对较少，仍有很大的研究空间。二、材料与方法2.1实验材料2.1.1玉米种质选择本研究选取了共计200份不同类型的玉米种质，涵盖了普通玉米、甜玉米、糯玉米、中国原生玉米、改良玉米以及国外引进玉米等多种类型。其中，普通玉米种质80份，分别来源于国内多个玉米主产区，如东北、华北、西南等地，这些种质在当地的玉米生产中具有广泛的种植基础，代表了不同生态环境下普通玉米的遗传特征。甜玉米种质30份，包含了常见的超甜、普甜等品种类型，来源既有国内自主选育的品种，也有从国外引进的优良种质。糯玉米种质30份，主要来源于我国糯玉米资源丰富的地区，如云南、贵州等地，这些种质在口感、淀粉品质等方面具有独特的优势。中国原生玉米种质20份，收集自我国一些少数民族聚居区或偏远山区，这些种质历经长期的自然选择和人工驯化，蕴含着丰富的遗传多样性，对于研究玉米的起源和进化具有重要意义。改良玉米种质20份，是通过现代育种技术对原有种质进行改良后的材料，在产量、品质、抗逆性等方面有了显著提升。国外引进玉米种质20份，来自美国、巴西、墨西哥等玉米生产大国，这些种质具有国外独特的遗传背景和优良性状，如美国的一些种质具有高产、耐密植的特点，巴西的种质在适应热带气候方面表现出色，墨西哥的种质则在耐旱性方面可能具有优势。选择这些不同类型的玉米种质，旨在全面涵盖玉米种质资源的遗传多样性，为后续的分子特征分析和耐旱相关性状研究提供丰富的材料基础，以便更深入地挖掘玉米的遗传潜力，解析耐旱性状的遗传机制。2.1.2实验地点与种植条件种植实验在[具体实验地点]的实验田进行。该实验田土壤类型为壤土，其质地适中，保水保肥能力较好，酸碱值在6.5-7.5之间，富含有机质，有利于玉米的生长发育。气候条件属于温带大陆性季风气候，四季分明，年平均气温为[X]℃，年降水量约为[X]mm，光照充足，能够满足玉米生长对光热的需求。实验设置了正常灌溉和干旱胁迫两个处理。正常灌溉处理下，根据玉米不同生育期的需水规律进行灌溉，保持土壤相对含水量在70%-80%之间。在玉米苗期，每隔3-5天进行一次灌溉，每次灌水量以湿透表层土壤10-15cm为宜；拔节期后，随着玉米生长加快，需水量增加，每隔2-3天灌溉一次，灌水量适当增加，确保土壤水分能够满足玉米生长需求。干旱胁迫处理则在玉米拔节期后开始实施，停止灌溉，使土壤自然失水，直至土壤相对含水量降至40%-50%，模拟干旱环境，以研究玉米在干旱胁迫下的生长表现和耐旱相关性状的变化。每个处理设置3次重复，随机区组排列，每个小区种植50株玉米，株行距为0.3m×0.6m，以保证玉米植株有足够的生长空间，减少个体之间的竞争对实验结果的影响。2.2实验方法2.2.1DNA提取与分子标记分析采用改良CTAB法提取玉米叶片DNA。具体步骤如下：取约0.5g新鲜的玉米叶片，迅速放入液氮中研磨成粉末状，以防止叶片中的DNA酶对DNA的降解。将研磨好的粉末转移至2mL离心管中，加入2mL预热至65℃的2×CTAB抽提液，该抽提液中含有2%（W/V）的CTAB、100mMTris-Cl（pH8.0）、20mMEDTA、1.4MNaCl，以及使用前加入的0.1%（V/V）的β-巯基乙醇。β-巯基乙醇是一种抗氧化剂，能够有效地防止酚氧化成醌，避免叶片组织褐变，同时使酚容易去除，从而保护DNA的完整性。充分混匀后，将离心管置于65℃水浴中保温30-60min，期间不时轻轻颠倒混匀，使叶片粉末与抽提液充分接触，促进DNA的释放和溶解。保温结束后，加入等体积的氯仿/异戊醇（24：1），氯仿能够使蛋白质变性沉淀，异戊醇则可以降低表面张力，减少抽提过程中泡沫的产生，有利于分层。轻缓颠倒混匀10-15min，使溶液充分乳化，然后在10000r/min的转速下离心5min，此时溶液会分为三层，上层为含有DNA的水相，中间为变性蛋白质层，下层为氯仿/异戊醇有机相。小心吸取上清液转移至新的离心管中，避免吸取到中间的蛋白质层。向上清液中加入1/10体积（约0.2mL）的65℃的10×CTAB/NaCl溶液（10%CTAB/0.7MNaCl），颠倒混匀，进一步沉淀多糖等杂质。再用等体积的氯仿/异戊醇（24：1）抽提一次，重复上述离心和转移上清的步骤。取上清，加入等体积的CTAB沉淀液（1%（W/V）CTAB、50mMTris-Cl、10mMEDTA），颠倒混匀，若有白色絮状沉淀可见，说明DNA已开始沉淀，可继续下一步操作；若沉淀不明显，将离心管置于65℃保温30min，促进DNA沉淀。在4℃下，以2700r/min的转速离心5min，使DNA沉淀于管底。弃去上清液，用高盐的TEbuffer（10mMTris-Cl，pH8.0；0.1mMEDTA；1MNaCl）重悬沉淀，重悬体积为0.25-0.5mL，可将离心管置于65℃保温30min，至大部分沉淀溶解。向重悬液中加入0.6体积的异丙醇沉淀核酸，充分混匀，在4℃下，以10000r/min的转速离心15min，使DNA进一步沉淀。弃去上清液，用80%乙醇洗涤沉淀2-3次，以去除残留的杂质和盐分，然后干燥沉淀，用尽可能少的TEbuffer重悬DNA，得到高质量的玉米基因组DNA。利用引物扩增和测序技术对SSR（简单序列重复）和InDel（插入/缺失）标记进行基因型鉴定。根据已公布的玉米基因组序列，设计覆盖玉米全基因组的SSR和InDel引物，引物设计时充分考虑引物的特异性、扩增效率和退火温度等因素，确保引物能够准确地扩增目标片段。PCR扩增体系为25μL，包括10×PCRbuffer2.5μL、2.5mMdNTPs2μL、上下游引物（10μM）各1μL、TaqDNA聚合酶0.5μL、模板DNA2μL，其余用ddH₂O补齐。PCR扩增程序为：94℃预变性5min；94℃变性30s，55-65℃退火30s（根据引物的退火温度进行调整），72℃延伸30s，共35个循环；最后72℃延伸10min。扩增产物通过8%聚丙烯酰胺凝胶电泳进行分离，银染显色，观察并记录条带的多态性，确定不同玉米种质在SSR和InDel标记位点的基因型。对于部分扩增效果不佳或条带难以判读的样本，采用测序技术进行验证，进一步准确确定其基因型。2.2.2耐旱相关性状测定在玉米生长过程中，测定多个耐旱相关性状。株高和穗位高的测定：在玉米乳熟期，使用标杆或专门的株高测量仪进行测量。对于株高，从地面测量至玉米植株顶部的最高处；穗位高则从地面测量至果穗着生节位的叶鞘基部。每个小区随机选取10株玉米进行测量，取平均值作为该小区的株高和穗位高数据。穗位叶叶绿素含量的测定：在玉米大喇叭口期至抽雄期，使用SPAD-502叶绿素仪测定穗位叶的叶绿素相对含量。选取穗位叶的中部位置，避开叶脉，每个叶片测定3个不同位点，取平均值作为该叶片的叶绿素含量，每个小区测定10株玉米的穗位叶，最终计算该小区的平均叶绿素含量。雌雄花期间隔（ASI）的测定：记录每株玉米雄穗散粉和雌穗吐丝的日期，两者日期之差即为雌雄花期间隔。每天上午9-11点和下午3-5点定时观察玉米的开花情况，确保数据记录的准确性，每个小区统计20株玉米的ASI值，计算平均值。雄穗大小的测定：在玉米雄穗完全抽出后，测量雄穗的长度、分枝数等指标。雄穗长度从雄穗基部测量至顶部，分枝数则直接计数。每个小区随机选取10株玉米，测量其雄穗相关指标，计算平均值以代表该小区的雄穗大小情况。果穗长度、秃尖长度、穗行数、行粒数、千粒重、穗粒重等性状的测定：在玉米成熟收获后，每个小区随机选取20个果穗。使用直尺测量果穗长度，从果穗基部到顶部的距离即为果穗长度；秃尖长度则是从果穗顶部未结实部分的起点测量至顶部，记录其长度。通过计数的方式统计穗行数和行粒数，将果穗上的籽粒全部脱下，数取1000粒籽粒称重，得到千粒重；将每个果穗的所有籽粒称重，得到穗粒重。对这些数据进行统计分析，以评估不同玉米种质在果穗相关性状上的表现及其与耐旱性的关系。2.2.3连锁-连锁不平衡联合作图利用软件plink进行连锁不平衡分析。连锁不平衡分析的原理是基于不同基因座位的等位基因在群体中出现的频次关系。当不同座位的某些等位基因在同一条染色体上出现的频次高于预期的随机频次时，就称这些基因处于连锁不平衡状态。在plink分析中，首先需要准备好ped和map格式的输入文件，ped文件包含了样本的个体信息和基因型数据，map文件则保存了分子标记位点的名称和位置信息。例如，map文件每一行代表一个标记位点，第一列是标记位点所在染色体的名字，第二列是标记位点的ID，第三列代表标记位点的遗传距离（若没有实际数值可用0填充），第四列是标记位点在染色体上的物理位置。在进行连锁不平衡分析时，设置相关参数。如设置--ld-window参数为10，限定一个标记位点最多和10个其他的标记位点进行LD分析；--ld-window-kb参数为1000，只对距离在1000kb之内的标记位点进行分析；--ld-window-r2参数为0.2（该参数只能和--r2参数搭配使用），对输出结果进行过滤，只输出R²大于0.2的LD分析结果。通过这些参数设置，可以有效地筛选出具有显著连锁不平衡关系的标记位点对，减少分析结果的冗余，提高分析效率和准确性。利用R语言中的ggplot2等工具进行作图。首先将plink分析得到的连锁不平衡数据导入R语言环境中，使用ggplot2包中的函数进行数据可视化。例如，使用ggplot()函数创建绘图对象，结合geom_point()函数绘制散点图，以展示不同标记位点之间的连锁不平衡程度（通常用R²值表示）与染色体位置的关系。通过设置x轴为染色体位置，y轴为R²值，不同颜色或形状的点代表不同的染色体，能够直观地呈现出全基因组范围内连锁不平衡的分布情况。还可以使用geom_hline()函数添加参考线，如设置参考线的y值为0.2，用于区分连锁不平衡程度较高和较低的区域。利用ggtitle()、xlab()、ylab()等函数添加图标题、坐标轴标签等注释信息，使图形更加清晰易懂，便于对连锁不平衡分析结果进行深入解读和分析，从而为后续定位与耐旱性状相关的基因位点提供有力的支持。三、不同类型玉米种质分子特征分析3.1遗传多样性指数计算本研究利用筛选出的SSR和InDel标记对200份玉米种质进行PCR扩增和基因型鉴定，获得了丰富的分子标记数据。基于这些数据，运用Popgene32软件计算遗传多样性指数，包括Nei's基因多样性指数（H）和Shannon信息指数（I）。Nei's基因多样性指数（H）的计算公式为：H=1-\sum_{i=1}^{n}p_{i}^{2}，其中p_{i}是第i个等位基因在群体中的频率，n为等位基因的数目。该指数反映了群体内基因的杂合程度，取值范围在0-1之间，H值越接近1，表明基因多样性越高；H值越接近0，则基因多样性越低。Shannon信息指数（I）的计算公式为：I=-\sum_{i=1}^{n}p_{i}ln(p_{i})，同样p_{i}是第i个等位基因在群体中的频率，n为等位基因的数目。Shannon信息指数不仅考虑了等位基因的频率，还能反映出群体中基因的丰富度和均匀度，其值越大，说明遗传多样性越丰富。通过这些公式的计算，能够准确地评估不同类型玉米种质在SSR和InDel标记位点上的遗传多样性水平，为后续的比较和分析提供量化的数据基础。3.2不同类型玉米种质遗传多样性比较对普通玉米、甜玉米、糯玉米、中国原生玉米、改良玉米以及国外引进玉米等不同类型玉米种质在SSR和InDel标记上的遗传多样性指数进行比较分析，结果如表1所示。普通玉米种质在SSR标记上的Nei's基因多样性指数平均为0.456，Shannon信息指数平均为0.689；在InDel标记上，Nei's基因多样性指数平均为0.423，Shannon信息指数平均为0.654。这表明普通玉米种质具有较高的遗传多样性，可能是由于其广泛的种植区域和长期的自然选择与人工选育，使其积累了丰富的遗传变异。甜玉米种质在SSR标记上的Nei's基因多样性指数平均为0.321，Shannon信息指数平均为0.502；在InDel标记上，Nei's基因多样性指数平均为0.298，Shannon信息指数平均为0.476。甜玉米种质的遗传多样性相对较低，这可能与甜玉米的育种目标较为单一，主要侧重于甜度等品质性状的改良，在一定程度上限制了遗传多样性的丰富度。糯玉米种质在SSR标记上的Nei's基因多样性指数平均为0.305，Shannon信息指数平均为0.481；在InDel标记上，Nei's基因多样性指数平均为0.287，Shannon信息指数平均为0.453。糯玉米种质的遗传多样性也相对较低，其原因可能与糯玉米的特殊用途有关，在育种过程中主要关注糯性等品质性状，对其他性状的选择相对较少，导致遗传背景相对狭窄。中国原生玉米种质在SSR标记上的Nei's基因多样性指数平均为0.489，Shannon信息指数平均为0.725；在InDel标记上，Nei's基因多样性指数平均为0.467，Shannon信息指数平均为0.698。中国原生玉米种质历经长期的自然选择和人工驯化，保留了丰富的遗传多样性，蕴含着许多独特的基因资源。改良玉米种质在SSR标记上的Nei's基因多样性指数平均为0.402，Shannon信息指数平均为0.605；在InDel标记上，Nei's基因多样性指数平均为0.378，Shannon信息指数平均为0.579。改良玉米种质虽然经过现代育种技术的改良，但在改良过程中可能会受到一些亲本材料的限制，导致遗传多样性有所降低。国外引进玉米种质在SSR标记上的Nei's基因多样性指数平均为0.435，Shannon信息指数平均为0.662；在InDel标记上，Nei's基因多样性指数平均为0.412，Shannon信息指数平均为0.637。国外引进玉米种质具有独特的遗传背景，但其遗传多样性受到引进种质数量和来源的影响，与中国原生玉米种质相比，遗传多样性略低。综上所述，不同类型玉米种质的遗传多样性存在显著差异。中国原生玉米种质和普通玉米种质的遗传多样性相对较高，而甜玉米和糯玉米种质的遗传多样性相对较低。这些差异为玉米种质资源的保护、利用和创新提供了重要的参考依据，在玉米育种中，可以根据不同的育种目标，合理选择具有丰富遗传多样性的种质资源，以拓宽育种材料的遗传基础，培育出具有优良性状的玉米新品种。【此处添加表1：不同类型玉米种质在SSR和InDel标记上的遗传多样性指数】玉米种质类型SSR标记Nei's基因多样性指数SSR标记Shannon信息指数InDel标记Nei's基因多样性指数InDel标记Shannon信息指数普通玉米0.4560.6890.4230.654甜玉米0.3210.5020.2980.476糯玉米0.3050.4810.2870.453中国原生玉米0.4890.7250.4670.698改良玉米0.4020.6050.3780.579国外引进玉米0.4350.6620.4120.6373.3聚类分析利用计算得到的遗传距离矩阵，采用非加权组平均法（UPGMA）对不同类型玉米种质进行聚类分析，绘制聚类图（图1）。从聚类图中可以清晰地看出不同类型玉米种质的亲缘关系。中国原生玉米种质单独聚为一大类，这表明中国原生玉米种质具有独特的遗传背景，与其他类型玉米种质的亲缘关系较远。这是由于中国原生玉米种质在长期的自然选择和人工驯化过程中，适应了当地的生态环境，形成了独特的遗传特征。普通玉米种质较为分散地分布在不同的分支中，这说明普通玉米种质的遗传多样性丰富，包含了多种不同的遗传类型，其来源广泛，在不同地区的种植和选育过程中，积累了丰富的遗传变异。甜玉米和糯玉米种质分别聚为相对独立的小类，且两者之间的亲缘关系相对较近。这是因为甜玉米和糯玉米在育种过程中，都侧重于特定品质性状的改良，导致它们的遗传背景相对狭窄，遗传相似度较高。同时，它们都属于特殊用途的玉米类型，在选育过程中可能采用了一些相似的育种技术和亲本材料，进一步增加了它们之间的亲缘关系。改良玉米种质与普通玉米种质在聚类图上有部分重叠，这表明改良玉米种质在遗传上与普通玉米种质具有一定的亲缘关系，改良过程是在普通玉米种质的基础上进行的。改良玉米种质在聚类图中也有一些独特的分支，这说明在改良过程中，通过引入新的基因或采用现代育种技术，使其获得了一些与普通玉米种质不同的遗传特征，从而在聚类图上表现出一定的差异。国外引进玉米种质在聚类图上与其他类型玉米种质有明显的区分，形成了相对独立的类群。这表明国外引进玉米种质具有独特的遗传背景，与国内的玉米种质存在较大的遗传差异。这是由于不同国家和地区的玉米育种目标、生态环境以及种质资源基础不同，导致国外引进玉米种质在遗传组成上与国内玉米种质有所不同。这种遗传差异为国内玉米育种提供了新的基因资源，可以通过杂交等手段将国外引进玉米种质的优良基因引入国内玉米种质中，拓宽国内玉米种质的遗传基础，提高玉米的综合性能。【此处添加图1：不同类型玉米种质的聚类图】四、耐旱相关性状分析4.1耐旱相关性状的表型变异对不同类型玉米种质在正常灌溉和干旱胁迫下的各耐旱相关性状进行统计分析，结果如表2所示。在株高方面，正常灌溉条件下，200份玉米种质的株高平均值为255.3cm，标准差为22.6cm，变异系数为8.85%；干旱胁迫下，株高平均值降至231.5cm，标准差为20.8cm，变异系数为9.00%。这表明干旱胁迫显著降低了玉米的株高，且不同种质间株高的变异程度在干旱胁迫下略有增加。穗位高在正常灌溉时平均值为105.6cm，标准差为15.2cm，变异系数为14.40%；干旱胁迫下平均值为92.3cm，标准差为13.8cm，变异系数为14.95%。干旱胁迫同样导致穗位高降低，且变异系数有所上升，说明不同种质在穗位高上的差异在干旱胁迫下更为明显。穗位叶叶绿素含量在正常灌溉条件下平均值为45.6SPAD值，标准差为5.2SPAD值，变异系数为11.40%；干旱胁迫下平均值降至38.9SPAD值，标准差为4.8SPAD值，变异系数为12.34%。这表明干旱胁迫降低了穗位叶叶绿素含量，且不同种质间的变异程度有所增大。雌雄花期间隔（ASI）在正常灌溉时平均值为2.5天，标准差为0.8天，变异系数为32.00%；干旱胁迫下平均值增加到4.8天，标准差为1.2天，变异系数为25.00%。干旱胁迫显著延长了雌雄花期间隔，虽然变异系数有所下降，但由于平均值的大幅增加，说明不同种质在ASI上对干旱胁迫的响应存在较大差异。雄穗长度在正常灌溉条件下平均值为30.5cm，标准差为4.5cm，变异系数为14.75%；干旱胁迫下平均值为26.8cm，标准差为4.2cm，变异系数为15.67%。干旱胁迫使雄穗长度缩短，且不同种质间的变异程度略有增加。果穗长度在正常灌溉时平均值为20.5cm，标准差为2.8cm，变异系数为13.66%；干旱胁迫下平均值降至17.3cm，标准差为2.5cm，变异系数为14.45%。干旱胁迫导致果穗长度显著缩短，不同种质间的变异程度也有所增加。秃尖长度在正常灌溉条件下平均值为1.2cm，标准差为0.5cm，变异系数为41.67%；干旱胁迫下平均值增加到2.8cm，标准差为0.8cm，变异系数为28.57%。干旱胁迫使秃尖长度明显增加，尽管变异系数有所下降，但由于平均值的大幅上升，表明不同种质在秃尖长度上对干旱胁迫的响应差异较大。穗行数在正常灌溉时平均值为16.5行，标准差为1.8行，变异系数为10.91%；干旱胁迫下平均值为16.2行，标准差为1.7行，变异系数为10.49%。干旱胁迫使穗行数略有减少，不同种质间的变异程度变化不大。行粒数在正常灌溉条件下平均值为35.6粒，标准差为4.2粒，变异系数为11.79%；干旱胁迫下平均值降至30.2粒，标准差为3.8粒，变异系数为12.58%。干旱胁迫导致行粒数显著减少，不同种质间的变异程度有所增加。千粒重在正常灌溉时平均值为350.5g，标准差为35.6g，变异系数为10.16%；干旱胁迫下平均值为302.8g，标准差为32.5g，变异系数为10.73%。干旱胁迫使千粒重明显降低，不同种质间的变异程度略有增加。穗粒重在正常灌溉条件下平均值为220.5g，标准差为30.6g，变异系数为13.88%；干旱胁迫下平均值降至156.8g，标准差为25.8g，变异系数为16.45%。干旱胁迫导致穗粒重显著降低，不同种质间的变异程度明显增加。综上所述，干旱胁迫对玉米各耐旱相关性状均产生了显著影响，不同类型玉米种质在各性状上存在丰富的表型变异，这些变异为玉米耐旱性的遗传研究和品种改良提供了丰富的素材。【此处添加表2：不同类型玉米种质在正常灌溉和干旱胁迫下耐旱相关性状的表型变异】性状处理平均值标准差变异系数(%)株高(cm)正常灌溉255.322.68.85干旱胁迫231.520.89.00穗位高(cm)正常灌溉105.615.214.40干旱胁迫92.313.814.95穗位叶叶绿素含量(SPAD值)正常灌溉45.65.211.40干旱胁迫38.94.812.34雌雄花期间隔(ASI)(天)正常灌溉2.50.832.00干旱胁迫4.81.225.00雄穗长度(cm)正常灌溉30.54.514.75干旱胁迫26.84.215.67果穗长度(cm)正常灌溉20.52.813.66干旱胁迫17.32.514.45秃尖长度(cm)正常灌溉1.20.541.67干旱胁迫2.80.828.57穗行数(行)正常灌溉16.51.810.91干旱胁迫16.21.710.49行粒数(粒)正常灌溉35.64.211.79干旱胁迫30.23.812.58千粒重(g)正常灌溉350.535.610.16干旱胁迫302.832.510.73穗粒重(g)正常灌溉220.530.613.88干旱胁迫156.825.816.454.2耐旱性评价指标筛选首先对各耐旱相关性状进行相关性分析，结果如表3所示。株高与穗位高呈极显著正相关（r=0.852**），这表明株高较高的玉米种质，其穗位高也往往较高。株高与穗位叶叶绿素含量呈显著正相关（r=0.325*），说明株高较高的种质，其穗位叶叶绿素含量相对较高，可能具有较强的光合作用能力。株高与雌雄花期间隔呈极显著负相关（r=-0.456**），即株高较高的玉米种质，其雌雄花期间隔相对较短，这可能与植株的生长发育速度和生理调节机制有关。穗位叶叶绿素含量与雌雄花期间隔呈极显著负相关（r=-0.489**），说明叶绿素含量较高的玉米种质，其雌雄花期间隔较短，可能在干旱胁迫下能更好地维持正常的生殖发育。穗位叶叶绿素含量与果穗长度呈显著正相关（r=0.356*），表明叶绿素含量高的种质，果穗长度可能更长，有利于提高产量。雌雄花期间隔与果穗长度呈极显著负相关（r=-0.567**），与秃尖长度呈极显著正相关（r=0.623**），与行粒数呈极显著负相关（r=-0.589**），与千粒重呈极显著负相关（r=-0.521**），与穗粒重呈极显著负相关（r=-0.789**）。这表明雌雄花期间隔越长，玉米的果穗长度越短，秃尖长度越长，行粒数、千粒重和穗粒重越低，对产量的影响较大。果穗长度与秃尖长度呈极显著负相关（r=-0.689**），与行粒数呈极显著正相关（r=0.723**），与千粒重呈显著正相关（r=0.389*），与穗粒重呈极显著正相关（r=0.891**）。说明果穗长度越长，秃尖长度越短，行粒数、千粒重和穗粒重越高，对产量的提升有积极作用。秃尖长度与行粒数呈极显著负相关（r=-0.712**），与千粒重呈极显著负相关（r=-0.456**），与穗粒重呈极显著负相关（r=-0.856**）。表明秃尖长度越长，行粒数、千粒重和穗粒重越低，严重影响玉米的产量。行粒数与千粒重呈显著正相关（r=0.356*），与穗粒重呈极显著正相关（r=0.923**）。千粒重与穗粒重呈极显著正相关（r=0.876**）。这些相关性表明，行粒数和千粒重的增加都有助于提高穗粒重，进而提高产量。为了进一步筛选耐旱性评价指标，进行主成分分析。主成分分析结果表明，前3个主成分的累计贡献率达到78.65%，能够较好地代表原始数据的信息。第一主成分主要反映了果穗长度、行粒数、千粒重和穗粒重等产量相关性状，贡献率为42.56%；第二主成分主要反映了株高、穗位高和穗位叶叶绿素含量等生长和光合相关性状，贡献率为23.45%；第三主成分主要反映了雌雄花期间隔和秃尖长度等与生殖发育和产量损失相关的性状，贡献率为12.64%。综合相关性分析和主成分分析结果，筛选出果穗长度、行粒数、千粒重、穗粒重、株高、穗位叶叶绿素含量、雌雄花期间隔和秃尖长度等8个指标作为玉米耐旱性的主要评价指标。这些指标涵盖了玉米生长、光合、生殖发育和产量等多个方面，能够较为全面地反映玉米的耐旱性。【此处添加表3：耐旱相关性状的相关性分析】性状株高穗位高穗位叶叶绿素含量雌雄花期间隔雄穗长度果穗长度秃尖长度穗行数行粒数千粒重穗粒重株高1穗位高0.852**1穗位叶叶绿素含量0.325*0.289*1雌雄花期间隔-0.456**-0.389**-0.489**1雄穗长度0.2560.2340.189-0.1561果穗长度0.289*0.267*0.356*-0.567**0.2011秃尖长度-0.321*-0.298*-0.301*0.623**-0.189-0.689**1穗行数0.1560.1340.102-0.1230.0890.112-0.0981行粒数0.267*0.245*0.312*-0.589**0.1780.723**-0.712**0.1561千粒重0.2230.2010.256-0.521**0.1670.389*-0.456**0.1230.356*1穗粒重0.356*0.321*0.389*-0.789**0.2010.891**-0.856**0.1890.923**0.876**1注：*表示在0.05水平上显著相关，**表示在0.01水平上极显著相关。4.3不同类型玉米种质耐旱性比较依据筛选出的耐旱性评价指标，采用隶属函数法对不同类型玉米种质的耐旱性进行综合评价。隶属函数值的计算公式为：U(X_{ij})=(X_{ij}-X_{jmin})/(X_{jmax}-X_{jmin})，其中U(X_{ij})为第i个品种第j个性状的隶属函数值，X_{ij}为第i个品种第j个性状的测定值，X_{jmax}和X_{jmin}分别为所有品种中第j个性状的最大值和最小值。对于与耐旱性呈负相关的性状，如雌雄花期间隔和秃尖长度，采用反隶属函数值计算：U(X_{ij})=1-(X_{ij}-X_{jmin})/(X_{jmax}-X_{jmin})。计算每个玉米种质在各耐旱性评价指标上的隶属函数值，并求其平均值作为该种质的耐旱性综合评价值（D值）。D值越大，表明该种质的耐旱性越强。不同类型玉米种质的耐旱性综合评价结果如表4所示。普通玉米种质的耐旱性综合评价值平均为0.562，其中部分种质如[具体种质名称1]的D值达到0.785，表现出较强的耐旱性；甜玉米种质的耐旱性综合评价值平均为0.325，整体耐旱性较弱，如[具体种质名称2]的D值仅为0.213；糯玉米种质的耐旱性综合评价值平均为0.308，与甜玉米种质相近，耐旱性相对较弱，[具体种质名称3]的D值为0.235；中国原生玉米种质的耐旱性综合评价值平均为0.625，其中[具体种质名称4]的D值高达0.856，表现出很强的耐旱性，这可能与中国原生玉米种质长期适应复杂的生态环境有关；改良玉米种质的耐旱性综合评价值平均为0.502，耐旱性处于中等水平，[具体种质名称5]的D值为0.654；国外引进玉米种质的耐旱性综合评价值平均为0.489，耐旱性也处于中等水平，[具体种质名称6]的D值为0.602。综上所述，中国原生玉米种质和普通玉米种质的耐旱性相对较强，甜玉米和糯玉米种质的耐旱性较弱，改良玉米种质和国外引进玉米种质的耐旱性处于中等水平。在玉米耐旱育种中，可以优先选择中国原生玉米种质和普通玉米种质作为亲本材料，同时加强对甜玉米和糯玉米种质耐旱性的改良。【此处添加表4：不同类型玉米种质的耐旱性综合评价结果】玉米种质类型样本数耐旱性综合评价值(D值)平均值范围代表种质及D值普通玉米800.5620.356-0.785[具体种质名称1]：0.785甜玉米300.3250.213-0.456[具体种质名称2]：0.213糯玉米300.3080.235-0.423[具体种质名称3]：0.235中国原生玉米200.6250.456-0.856[具体种质名称4]：0.856改良玉米200.5020.389-0.654[具体种质名称5]：0.654国外引进玉米200.4890五、连锁-连锁不平衡关系分析5.1连锁不平衡程度评估连锁不平衡（LinkageDisequilibrium，LD）是指在群体中，位于染色体上相距较近的基因位点之间，等位基因非随机地共同出现的现象，其程度的评估对于理解基因组结构和遗传变异至关重要。本研究运用软件plink对筛选出的SSR和InDel标记进行连锁不平衡分析，选用连锁不平衡系数D'和r^{2}作为评估连锁不平衡程度的关键参数。D'的计算基于两个位点等位基因频率的联合分布，其取值范围为-1到1。当D'=1时，表示两个位点处于完全连锁不平衡状态，即它们的等位基因总是一起遗传，没有发生重组；当D'=0时，则表示两个位点之间不存在连锁不平衡，等位基因的遗传是完全随机的。r^{2}则反映了两个位点之间等位基因的相关性，取值范围从0到1。r^{2}值越接近1，说明两个位点之间的连锁不平衡程度越高，等位基因的相关性越强；r^{2}值越接近0，表明两个位点之间的连锁不平衡程度越低，等位基因几乎是独立遗传的。在实际分析中，设置了一系列参数以确保分析的准确性和有效性。例如，设置--ld-window参数为10，这意味着限定一个标记位点最多和10个其他的标记位点进行LD分析，避免了过多不必要的计算，提高了分析效率。设置--ld-window-kb参数为1000，只对距离在1000kb之内的标记位点进行分析，因为通常情况下，距离较近的位点更容易存在连锁不平衡关系，这样的设置可以更有针对性地分析连锁不平衡情况。设置--ld-window-r2参数为0.2（该参数只能和--r2参数搭配使用），对输出结果进行过滤，只输出r^{2}大于0.2的LD分析结果，从而筛选出连锁不平衡程度较为显著的标记位点对，减少了冗余信息，便于后续的深入分析。利用R语言中的ggplot2工具对连锁不平衡分析结果进行可视化展示。在全基因组水平上，绘制连锁不平衡程度（以r^{2}值表示）与染色体位置的关系图（图2）。从图中可以清晰地看到，连锁不平衡程度在不同染色体上呈现出不均匀的分布。在某些染色体区域，如玉米第1染色体的50-100Mb区间以及第4染色体的150-200Mb区间，r^{2}值较高，表明这些区域存在较强的连锁不平衡。这可能是由于这些区域内基因之间的重组频率较低，使得等位基因能够保持相对稳定的组合，从而维持了较高的连锁不平衡程度。在其他一些区域，r^{2}值较低，连锁不平衡程度较弱，这可能是因为这些区域的重组活动较为频繁，打破了等位基因之间的非随机关联，使得它们的遗传趋于独立。【此处添加图2：连锁不平衡程度在全基因组上的分布图】5.2连锁不平衡与遗传多样性关系遗传多样性是物种适应环境变化和进化的基础，而连锁不平衡与遗传多样性之间存在着密切的关联。本研究深入分析了连锁不平衡程度与遗传多样性指数之间的相关性，以揭示它们对玉米种质遗传结构的影响。通过计算不同类型玉米种质在SSR和InDel标记上的遗传多样性指数，包括Nei's基因多样性指数（H）和Shannon信息指数（I），并与连锁不平衡程度（r^{2}值）进行相关性分析。结果显示，连锁不平衡程度与遗传多样性指数之间存在显著的正相关关系（图3）。具体而言，在遗传多样性较高的玉米种质群体中，如中国原生玉米种质和普通玉米种质，连锁不平衡程度也相对较高。这是因为遗传多样性丰富的群体中，存在更多的等位基因和遗传变异，这些变异在染色体上的分布使得基因之间更容易形成非随机的关联，从而导致较高的连锁不平衡程度。在遗传多样性较低的群体，如甜玉米和糯玉米种质中，连锁不平衡程度也相对较低。这是由于这些群体在育种过程中，遗传背景相对狭窄，等位基因的种类和数量较少，基因之间发生非随机关联的机会也相应减少，导致连锁不平衡程度较低。【此处添加图3：连锁不平衡程度与遗传多样性指数的相关性分析图】连锁不平衡对玉米种质遗传结构的影响是多方面的。较高的连锁不平衡意味着某些基因位点之间的等位基因组合相对稳定，这在一定程度上限制了遗传重组的发生，使得一些有利的基因组合能够在群体中得以保存。在玉米的进化过程中，一些适应特定环境的基因组合可能由于连锁不平衡而被固定下来，有助于玉米在特定环境中生存和繁衍。然而，过高的连锁不平衡也可能导致遗传多样性的降低，因为一些不利的等位基因也可能随着连锁不平衡而被保留在群体中，限制了群体对环境变化的适应能力。相反，较低的连锁不平衡虽然有利于遗传重组的发生，增加遗传多样性，但也可能导致一些有利的基因组合被打破，影响玉米的优良性状。因此，在玉米种质资源的保护和利用中，需要充分考虑连锁不平衡与遗传多样性的关系，合理利用连锁不平衡现象，保护和利用遗传多样性，以实现玉米种质的改良和创新。5.3与耐旱性状的关联分析通过连锁-连锁不平衡联合作图，能够将分子标记与耐旱性状紧密联系起来，挖掘与耐旱性状显著关联的分子标记区间，进而确定潜在的耐旱相关基因，这对于深入理解玉米耐旱的遗传机制和开展耐旱育种具有重要意义。运用软件plink进行连锁-连锁不平衡联合作图分析时，将前期测定的玉米各耐旱相关性状（如株高、穗位高、穗位叶叶绿素含量、雌雄花期间隔、果穗长度、秃尖长度、穗行数、行粒数、千粒重、穗粒重等）的表型数据与SSR和InDel标记的基因型数据相结合。在分析过程中，采用一般线性模型（GLM）或混合线性模型（MLM）进行关联分析。一般线性模型假设群体结构和个体间的亲缘关系对性状的影响可以忽略不计，通过简单的线性回归来检测标记与性状之间的关联。混合线性模型则充分考虑了群体结构和个体间的亲缘关系对性状的影响，将这些因素作为随机效应纳入模型中，能够更准确地检测到与性状真正关联的标记位点，减少假阳性结果的出现。经过严格的统计分析，设置合适的显著性阈值（如P<0.01），筛选出与耐旱性状显著关联的分子标记区间。结果发现，在玉米第3染色体的70-80Mb区间，存在多个与果穗长度、行粒数和穗粒重显著关联的SSR和InDel标记。进一步对该区间进行基因注释和功能分析，发现其中包含一个可能与耐旱相关的基因ZmDREB1，该基因编码一个脱水响应元件结合蛋白，在植物响应干旱胁迫的过程中发挥着重要的调控作用。在玉米第6染色体的120-130Mb区间，有与穗位叶叶绿素含量显著关联的标记，该区间内的基因ZmP5CS可能参与了脯氨酸的合成，而脯氨酸在植物应对干旱胁迫时能够调节细胞的渗透势，增强植物的耐旱性。【此处添加图4：与耐旱性状显著关联的分子标记区间在染色体上的分布图】这些与耐旱性状显著关联的分子标记区间和潜在的耐旱相关基因，为玉米耐旱育种提供了重要的靶点。在后续的育种工作中，可以利用这些分子标记进行分子标记辅助选择（MAS），快速准确地筛选出具有优良耐旱性状的玉米种质，提高育种效率。还可以进一步深入研究这些耐旱相关基因的功能和作用机制，通过基因编辑等现代生物技术手段，对玉米的耐旱性进行精准改良，为培育适应干旱环境的玉米新品种奠定坚实的基础。六、结果讨论6.1不同类型玉米种质分子特征的特异性及意义不同类型玉米种质在分子特征上呈现出明显的特异性，这些差异有着深刻的遗传学基础。从遗传多样性角度来看，中国原生玉米种质和普通玉米种质具有较高的遗传多样性。中国原生玉米种质历经长期的自然选择和人工驯化，在复杂的生态环境中，各种自然因素如不同的气候条件、土壤类型等，以及人工的选择压力，促使其基因组不断发生变异和进化，积累了丰富的等位基因和遗传变异，从而形成了较高的遗传多样性。普通玉米种质由于广泛的种植区域，在不同地区适应了当地的生态环境，同时在长期的人工选育过程中，也引入了多种遗传资源，进一步丰富了其遗传背景。甜玉米和糯玉米种质的遗传多样性相对较低，这与它们的育种目标和选育方式密切相关。甜玉米主要侧重于甜度等品质性状的改良，在育种过程中，为了获得特定的甜度性状，往往会对一些关键基因进行选择和固定，这在一定程度上限制了其他基因的变异和多样性。糯玉米主要关注糯性等品质性状，选育过程相对集中在少数几个与糯性相关的基因上，导致其遗传背景相对狭窄，遗传多样性降低。这些分子特征差异对玉米种质资源保护、利用和新品种选育具有重要意义。在种质资源保护方面，对于遗传多样性丰富的中国原生玉米种质，应加强保护力度，建立专门的保护区或种质库，防止其遗传资源的流失。这是因为这些种质蕴含着许多独特的基因，对于未来玉米的遗传改良和应对环境变化具有重要价值。对于遗传多样性较低的甜玉米和糯玉米种质，也需要进行保护，同时可以通过与其他种质杂交等方式，引入新的基因，拓宽其遗传背景。在种质资源利用和新品种选育方面，不同类型玉米种质的分子特征差异为育种提供了丰富的素材。可以根据不同的育种目标，选择合适的种质资源进行杂交。若要培育高产且具有一定耐旱性的玉米品种，可以选择遗传多样性丰富的普通玉米种质作为亲本，利用其丰富的遗传变异，结合耐旱性较好的种质，通过杂交和选择，有望获得理想的新品种。还可以利用分子标记技术，深入研究不同类型玉米种质分子特征与性状之间的关系，加速新品种的选育进程，提高育种效率。6.2耐旱相关性状的遗传基础及影响因素玉米耐旱相关性状是由多基因控制的复杂数量性状，其遗传基础涉及众多基因位点的协同作用。本研究通过连锁-连锁不平衡联合作图，发现了多个与耐旱性状显著关联的分子标记区间和潜在的耐旱相关基因，如ZmDREB1和ZmP5CS等。这些基因在植物响应干旱胁迫的过程中发挥着重要作用，ZmDREB1编码的脱水响应元件结合蛋白，能够与干旱胁迫响应基因的启动子区域结合，激活这些基因的表达，从而增强植物对干旱胁迫的耐受性；ZmP5CS参与脯氨酸的合成，脯氨酸作为一种渗透调节物质，能够调节细胞的渗透势，在干旱胁迫下维持细胞的正常生理功能，增强植物的耐旱性。环境因素对玉米耐旱性表现有着显著影响。在干旱胁迫下，玉米的生长发育受到抑制，各耐旱相关性状发生明显变化。干旱会导致玉米株高降低，这是因为干旱抑制了细胞的伸长和分裂，影响了植株的正常生长。干旱使穗位叶叶绿素含量下降，降低了光合作用效率，影响了植株的能量供应和物质合成。干旱还会导致雌雄花期间隔延长，影响授粉结实，进而降低产量。这是由于干旱胁迫影响了植物体内激素的平衡和信号传导，导致雄穗和雌穗的发育不同步。不同生态环境下，玉米的耐旱性表现也存在差异。在干旱地区种植的玉米，由于长期适应干旱环境，可能在形态、生理和分子水平上发生了适应性变化，使其耐旱性相对较强；而在湿润地区种植的玉米，对干旱胁迫的适应能力相对较弱。因此，在研究玉米耐旱性时，需要充分考虑遗传因素和环境因素的相互作用。通过深入研究耐旱相关性状的遗传基础，挖掘更多的耐旱基因，并结合环境因素，开展适应性育种，能够培育出更适应不同环境条件的耐旱玉米品种。在育种过程中，可以利用分子标记辅助选择技术，选择具有优良耐旱基因组合的材料，同时在不同的生态环境下进行多点试验，筛选出在不同环境下都能稳定表现出耐旱性的品种，提高玉米在干旱环境下的产量和稳定性。6.3连锁-连锁不平衡联合作图的有效性及应用前景连锁-连锁不平衡联合作图在挖掘玉米耐旱相关基因方面展现出了较高的有效性和独特的优势。与传统的连锁分析相比，它不仅能够利用家系群体中标记与性状之间的连锁关系，还能充分考虑自然群体中标记位点之间的连锁不平衡信息。在自然群体中，由于经历了长期的进化和繁殖，基因之间的连锁关系会发生重组和变化，但连锁不平衡信息能够反映出基因之间的历史关联和协同进化关系。通过连锁-连锁不平衡联合作图，可以更全面地扫描基因组，增加检测到与耐旱性状相关基因位点的概率，提高基因定位的准确性。在本研究中，通过连锁-连锁不平衡联合作图，成功地定位到了多个与玉米耐旱性状显著关联的分子标记区间和潜在的耐旱相关基因，如在玉米第3染色体的70-80Mb区间发现了与果穗长度、行粒数和穗粒重显著关联的标记，以及该区间内可能与耐旱相关的基因ZmDREB1。这表明该方法能够有效地挖掘出与耐旱性状相关的遗传信息，为深入研究玉米耐旱的遗传机制提供了有力的工具。展望其在玉米分子育种中的应用前景，连锁-连锁不平衡联合作图具有广阔的发展空间。在分子标记辅助选择育种中，利用定位到的与耐旱性状紧密关联的分子标记，可以快速准确地筛选出具有优良耐旱性状的玉米种质，大大提高育种效率，缩短育种周期。通过对这些分子标记的检测，可以在早期对育种材料进行筛选，减少不必要的田间种植和鉴定工作，降低育种成本。该方法还可以为基因编辑等现代生物技术提供靶点。通过对定位到的耐旱相关基因进行深入研究，了解其功能和作用机制，利用基因编辑技术对这些基因进行精准修饰，有望培育出具有更强耐旱性的玉米新品种，满足农业生产对耐旱玉米品种的需求，为保障粮食安全和农业可持续发展做出贡献。6.4研究的创新点与不足本研究在方法和结果等方面具有一定的创新之处。在方法上，首次将连锁-连锁不平衡联合作图应用于多种类型玉米种质的耐旱相关性状研究。以往的研究大多侧重于单一类型玉米种质或采用单一的连锁分析方法，本研究综合运用连锁分析和连锁不平衡分析，充分利用了自然群体的遗传信息，为玉米耐旱相关基因的挖掘提供了更全面、准确的方法。在实验材料上，选取了涵盖普通玉米、甜玉米、糯玉米、中