醋栗番茄渐渗系群体表型解析与优异性状QTL定位研究

醋栗番茄渐渗系群体表型解析与优异性状QTL定位研究一、引言1.1研究背景与意义番茄（Solanumlycopersicum）作为全球最重要的蔬菜作物之一，不仅是果实遗传、发育和生理学研究的经典模式系统，还在人类饮食结构中占据重要地位。随着人们生活水平的提高和对健康饮食的追求，对番茄的产量、品质、抗性等方面提出了更高的要求。然而，现代栽培番茄在长期驯化和人工选择过程中，经历了强烈的遗传瓶颈，导致品种间遗传多样性大幅降低，许多优异的农艺性状逐渐“丢失”，极大地制约了番茄育种进程。野生番茄在各异的生态环境下，展现出了丰富的遗传和表型多样性，蕴含着诸多栽培番茄所缺乏的优异性状。例如，野生番茄通常具备高抗病抗逆性，能有效抵御多种病虫害和不良环境胁迫；其可溶性固形物、番茄红素和果实风味物质含量较高，果实风味浓郁；部分野生番茄还具有高分枝能力，为提高番茄产量提供了潜在的遗传基础。野生番茄的这些优异特性，使其成为番茄品种改良的重要基因资源库。醋栗番茄（Solanumpimpinellifolium）作为栽培番茄的野生祖先，在番茄遗传改良中具有不可替代的关键作用。它携带大量抗病、抗逆等优良基因，这些基因能够显著增强栽培番茄对生物和非生物胁迫的抵抗能力。在抗病方面，能有效抵御黄萎病、斑枯病、枯萎病、烟草花叶病毒病等多种病害；在抗逆方面，对盐胁迫、干旱胁迫、低温胁迫等不良环境表现出较强的耐受性。如加拉帕戈斯番茄具有非凡的耐盐性，能够在70%的海水灌溉下生存结果，为培育耐盐番茄品种提供了宝贵的基因资源；智利番茄根系发达，抗旱能力强，对枯萎病近乎免疫，是改良栽培番茄抗旱性和抗病性的优质种质。渐渗系群体是通过将野生种与栽培种进行杂交、多次回交和自交等一系列遗传操作构建而成的。在这个过程中，野生种的染色体片段被逐步导入到栽培种的遗传背景中，形成了一套包含不同野生染色体片段的渐渗系。这些渐渗系群体不仅保留了栽培种的优良农艺性状，还引入了野生种的优异性状，为番茄遗传研究和品种改良提供了理想材料。利用醋栗番茄渐渗系群体，科研人员可以深入剖析野生番茄优良性状的遗传机制，精准定位和克隆相关基因或数量性状位点（QTL）。通过对这些基因和QTL的功能研究，进一步揭示番茄性状形成的分子生物学基础，为番茄分子标记辅助育种和基因编辑育种提供坚实的理论依据。在实际应用中，醋栗番茄渐渗系群体能够直接用于番茄品种的改良。通过筛选渐渗系群体中具有目标优异性状的株系，将其与现有栽培品种进行杂交和选育，可以培育出兼具高产、优质、抗病、抗逆等多种优良性状的番茄新品种，满足市场对高品质番茄的需求。利用渐渗系群体中携带的抗病基因，培育出对多种病害具有抗性的番茄品种，减少农药使用，降低生产成本，同时保障番茄的产量和品质；利用抗逆基因培育出适应不同生态环境的番茄品种，扩大番茄的种植区域，提高番茄的生产稳定性。对醋栗番茄渐渗系群体表型性状进行分析，并对优异性状进行QTL定位，具有重要的理论和实践意义。在理论层面，有助于深入了解番茄性状的遗传规律和分子调控机制，丰富番茄遗传学理论知识；在实践层面，能够为番茄分子育种提供关键的基因资源和技术支撑，加速番茄品种改良进程，推动番茄产业的可持续发展。1.2国内外研究现状野生番茄作为番茄遗传改良的重要基因资源，其渐渗系群体的构建、表型分析及QTL定位一直是国内外研究的热点。在醋栗番茄渐渗系群体构建方面，国外起步较早且研究较为深入。Tanksley等利用栽培番茄M82和醋栗番茄LA716为亲本，通过多代回交和自交，构建了一套包含107个渐渗系的群体，该群体几乎覆盖了醋栗番茄的整个基因组，为后续的基因定位和功能研究提供了重要材料。国内学者也在积极开展相关工作，如南京农业大学以栽培番茄Jina为母本，醋栗番茄LA2093为父本，成功构建了渐渗系BC3F3代群体，并利用该群体进行了苗期耐盐性评价等研究。在表型性状分析方面，国内外研究涉及醋栗番茄渐渗系的多个方面。在果实品质方面，研究发现渐渗系果实的可溶性固形物、番茄红素、维生素C等含量显著高于栽培番茄，且果实风味物质种类和含量也发生了明显变化，为改善番茄果实品质提供了新的途径；在抗病抗逆性方面，渐渗系对多种病害如黄萎病、枯萎病、烟草花叶病毒病等表现出较强的抗性，对盐胁迫、干旱胁迫、低温胁迫等非生物胁迫的耐受性也显著增强。在生长发育方面，部分渐渗系表现出分枝能力强、生长势旺盛等特点，为提高番茄产量提供了新的遗传资源。在优异性状QTL定位方面，国外研究利用分子标记技术，对醋栗番茄渐渗系的多个优异性状进行了QTL定位。Foolad等利用RFLP标记，将7个番茄芽期耐盐QTL定位在1、2、5、7、9、12号染色体上；Eshed等通过构建渐渗系群体，定位到多个与果实品质相关的QTL，为番茄品质改良提供了重要的理论依据。国内研究也取得了一定成果，潘颖等利用高代回交群体，定位了野生醋栗番茄发芽期耐盐QTL；刘希艳对醋栗番茄RIL群体苗期耐盐耐旱QTL进行了定位，为番茄抗逆育种提供了关键的基因资源。尽管国内外在醋栗番茄渐渗系群体的研究上取得了显著进展，但仍存在一些不足之处。在渐渗系群体构建方面，部分群体存在野生染色体片段覆盖不完整、遗传背景不纯等问题，影响了后续研究的准确性和有效性；在表型性状分析方面，研究主要集中在少数几个性状，对其他重要性状如抗虫性、营养成分等的研究相对较少，且表型鉴定方法不够标准化和精细化，导致数据准确性和可比性存在一定差异；在QTL定位方面，定位的QTL区间较大，分辨率较低，难以准确克隆和利用相关基因，且对QTL的遗传效应和互作机制研究不够深入，限制了其在番茄分子育种中的应用。1.3研究目标与内容本研究旨在深入解析醋栗番茄渐渗系群体的表型性状，并对其中的优异性状进行QTL定位，为番茄遗传改良提供坚实的理论基础和丰富的基因资源。具体研究内容如下：醋栗番茄渐渗系群体表型性状鉴定与分析：全面调查醋栗番茄渐渗系群体的多个表型性状，涵盖生长发育、产量、品质、抗病抗逆等方面。在生长发育性状上，详细记录株高、茎粗、分枝数、叶片数、叶片大小、节间长度、生长周期等指标，分析其生长动态和发育规律；产量相关性状方面，统计单株产量、单果重、果数、坐果率、果实大小分布等，评估不同株系的产量潜力；果实品质性状上，精确测定可溶性固形物、可溶性糖、可滴定酸、维生素C、番茄红素、果实硬度、果实风味物质等含量，评价果实的口感、营养和风味；抗病抗逆性状方面，通过人工接种病原菌或模拟逆境条件，鉴定渐渗系对黄萎病、枯萎病、烟草花叶病毒病、盐胁迫、干旱胁迫、低温胁迫等的抗性反应，筛选出具有优良抗病抗逆性的株系。运用统计学方法，如方差分析、相关性分析、主成分分析、聚类分析等，对调查得到的表型数据进行深入分析。方差分析用于明确不同渐渗系间各性状的差异显著性，判断性状变异是否由遗传因素主导；相关性分析探究各性状之间的内在关联，为后续QTL定位提供参考依据；主成分分析将多个表型性状综合转化为少数几个主成分，简化数据结构，揭示性状间的潜在关系；聚类分析依据性状相似性对渐渗系进行分类，挖掘具有特殊性状组合的株系，为番茄育种提供目标材料。醋栗番茄渐渗系群体分子标记分析：选用合适的分子标记技术，如简单序列重复（SSR）、单核苷酸多态性（SNP）等，对醋栗番茄渐渗系群体进行全基因组扫描。SSR标记具有多态性高、共显性遗传、操作简便等优点，通过设计特异性引物，扩增基因组中的SSR位点，检测其多态性；SNP标记则具有分布广泛、数量多、稳定性高等特点，利用高通量测序技术或基因芯片技术，对渐渗系群体进行SNP检测。构建高密度的分子遗传连锁图谱，确定分子标记在染色体上的相对位置和遗传距离。利用遗传分析软件，如JoinMap等，对分子标记数据进行处理和分析，构建分子遗传连锁图谱。图谱的构建精度和密度直接影响QTL定位的准确性和分辨率，高密度的图谱能够更精确地定位QTL，为后续基因克隆和功能研究奠定基础。醋栗番茄优异性状QTL定位：利用已构建的分子遗传连锁图谱和表型数据，运用QTL定位方法，如复合区间作图法（CIM）、完备区间作图法（ICIM）等，对醋栗番茄渐渗系群体中的优异性状进行QTL定位。CIM方法在考虑多个背景标记的基础上，通过逐步回归筛选与目标性状相关的标记，进行区间作图，有效控制了遗传背景的干扰；ICIM方法则进一步改进，在区间作图过程中同时考虑多个标记的效应，提高了QTL定位的准确性和效率。确定优异性状相关QTL在染色体上的位置、遗传效应和贡献率。通过分析QTL与分子标记之间的连锁关系，确定QTL在染色体上的具体位置；评估QTL的加性效应、显性效应和上位性效应，了解其遗传规律；计算QTL对性状变异的贡献率，判断其在性状表现中的重要性。对定位到的QTL进行验证和比较分析，确保QTL的真实性和可靠性。利用不同年份、不同环境下的表型数据进行重复验证，或通过构建近等基因系等方法，对QTL进行功能验证；与已报道的QTL进行比较分析，明确本研究中QTL的独特性和创新性，为番茄分子育种提供更有价值的基因资源。二、材料与方法2.1试验材料本研究选用的母本为栽培番茄品种Jina，其具有良好的农艺性状，如果实大小适中、产量稳定、品质优良等，是目前广泛种植的番茄品种之一，为后续渐渗系群体提供了优良的遗传背景，能够保证渐渗系在保留野生番茄优异性状的同时，具备较好的栽培适应性。父本为醋栗番茄LA2093，该品种携带多种优异性状基因，在抗病方面，对黄萎病、枯萎病、烟草花叶病毒病等常见病害具有较强的抗性；在抗逆方面，具有良好的盐胁迫、干旱胁迫、低温胁迫耐受性，是番茄遗传改良的优质基因供体。以栽培番茄Jina为母本，醋栗番茄LA2093为父本进行杂交，获得F1代种子。将F1代作为母本，与栽培番茄Jina进行回交，得到BC1F1代。按照此方法，连续回交3次，获得BC3F1代。之后，对BC3F1代进行自交，经过多代自交和筛选，最终构建了包含19个株系的醋栗番茄渐渗系BC3F3代群体。在构建过程中，每一代都进行了严格的表型鉴定和分子标记检测，以确保渐渗系群体的质量和遗传稳定性。通过对回交和自交后代的筛选，使野生醋栗番茄的染色体片段逐步渗入到栽培番茄的遗传背景中，形成了具有不同野生染色体片段的渐渗系株系，为后续的表型性状分析和QTL定位提供了丰富的遗传材料。2.2表型性状调查在番茄的整个生长周期中，对醋栗番茄渐渗系BC3F3代群体的19个株系及双亲进行多方面表型性状调查，每个株系选取10株生长状况良好且具有代表性的植株作为调查样本，以确保数据的可靠性和准确性。在生长发育性状方面，从番茄幼苗期开始，每隔7天使用直尺测量株高，精确到0.1cm；使用游标卡尺测量茎粗，精确到0.01cm；定期记录分枝数、叶片数，并使用直尺测量叶片的长和宽，计算叶片大小，单位为cm²；测量相邻两节之间的距离，得到节间长度，精确到0.1cm；记录从播种到第一朵花开放的天数、从第一朵花开放到果实成熟的天数等，确定生长周期。产量相关性状在果实成熟后进行统计。单株产量通过电子秤称量，精确到0.1g；随机选取20个果实，使用电子天平测量单果重，精确到0.01g；统计每株的果数，计算坐果率；使用游标卡尺测量果实的纵径和横径，以评估果实大小分布。果实品质性状测定在果实完全成熟时进行。利用手持折光仪测定可溶性固形物含量，单位为°Bx；采用蒽酮比色法测定可溶性糖含量，单位为mg/g；利用酸碱中和滴定法测定可滴定酸含量，单位为g/100g；使用2,6-二氯靛酚滴定法测定维生素C含量，单位为mg/100g；利用高效液相色谱法测定番茄红素含量，单位为mg/kg；使用果实硬度计测定果实硬度，单位为N；采用顶空固相微萃取-气相色谱-质谱联用技术（HS-SPME-GC-MS）测定果实风味物质，鉴定出多种挥发性风味物质，并对其含量进行分析。抗病抗逆性状鉴定通过人工接种病原菌或模拟逆境条件来进行。黄萎病抗性鉴定采用灌根接种法，将浓度为1×10⁶CFU/mL的大丽轮枝菌孢子悬浮液浇灌到番茄根部，接种后21天调查发病情况，根据病情指数进行抗性评价；枯萎病抗性鉴定采用蘸根接种法，将番茄幼苗根部蘸取浓度为1×10⁶CFU/mL的尖孢镰刀菌孢子悬浮液，接种后14天调查发病情况，依据病情指数评估抗性；烟草花叶病毒病抗性鉴定采用摩擦接种法，将烟草花叶病毒汁液涂抹在番茄叶片上，接种后7天观察发病症状，按照发病级别进行抗性分级。盐胁迫处理在苗期进行，使用200mmol/LNaCl溶液浇灌番茄幼苗，处理14天后测定株高、茎粗、整株鲜质量、整株干质量、叶绿素相对含量、丙二醛含量和过氧化氢含量等指标，计算耐盐性系数；干旱胁迫处理在开花期进行，停止浇水，当土壤相对含水量降至30%时，持续处理7天，测定叶片相对含水量、脯氨酸含量、可溶性糖含量等指标，评估抗旱性；低温胁迫处理在结果期进行，将番茄植株放入温度为5℃的人工气候箱中，处理7天，测定果实坐果率、果实品质等指标，判断耐低温能力。2.3QTL定位方法QTL定位的原理是基于分子标记与数量性状位点之间的连锁关系。通过对分离群体中个体的分子标记基因型和数量性状表型数据进行分析，确定QTL在染色体上的位置以及其遗传效应。当某个分子标记与控制目标性状的QTL紧密连锁时，标记基因型与性状表型之间会呈现出显著的相关性，从而可以推断QTL位于该标记附近。本研究采用完备区间作图法（ICIM）进行QTL定位。ICIM方法是在复合区间作图法（CIM）的基础上发展而来，它综合考虑了多个标记的效应，有效控制了遗传背景的干扰，提高了QTL定位的准确性和效率。在进行ICIM分析时，将利用QTLIciMapping软件V4.1。该软件由中国农科院王建康老师数量遗传课题组开发，专门用于QTL定位分析。在软件操作过程中，会将构建好的分子遗传连锁图谱数据和表型数据导入软件，设置相关参数，如LOD阈值、步长、背景标记选择等。LOD阈值用于判断QTL的显著性，一般设置为2.5-3.0，当LOD值大于该阈值时，认为存在显著的QTL；步长决定了在染色体上扫描的精细程度，较小的步长可以提高定位的精度，但计算量也会相应增加；背景标记选择则根据实际情况，选取与目标性状相关性较小的标记作为背景标记，以消除遗传背景的影响。软件会按照ICIM算法，对整个基因组进行扫描，检测QTL的存在，并确定其在染色体上的位置、遗传效应和贡献率。三、醋栗番茄渐渗系群体表型性状分析3.1表型性状的描述性统计对醋栗番茄渐渗系BC3F3代群体的19个株系及双亲的各表型性状进行统计分析，结果如表1所示。表1醋栗番茄渐渗系群体表型性状的描述性统计性状最小值最大值平均值标准差变异系数（%）株高（cm）56.2123.589.415.317.1茎粗（cm）0.821.561.140.1815.8分枝数3.08.05.21.223.1叶片数28.056.040.57.819.3叶片大小（cm²）12.535.622.85.423.7节间长度（cm）4.59.26.81.319.1生长周期（天）1021351188.67.3单株产量（g）356.2895.4568.3125.622.1单果重（g）18.545.630.26.822.5果数15.035.024.55.622.9坐果率（%）65.090.078.57.69.7果实纵径（cm）3.25.64.30.716.3果实横径（cm）3.56.04.80.816.7可溶性固形物（°Bx）5.58.26.80.913.2可溶性糖（mg/g）22.545.632.86.519.8可滴定酸（g/100g）0.350.760.520.1121.2维生素C（mg/100g）18.535.625.84.617.8番茄红素（mg/kg）10.525.616.83.520.8果实硬度（N）4.58.26.21.117.7黄萎病病情指数10.060.035.012.535.7枯萎病病情指数5.045.025.09.538.0烟草花叶病毒病病情指数8.050.030.010.535.0盐胁迫耐盐性系数0.450.860.650.1116.9干旱胁迫抗旱性系数0.520.900.720.1216.7低温胁迫耐低温性系数0.500.850.680.1014.7从表1可以看出，各表型性状在渐渗系群体中均表现出丰富的变异。株高的变化范围为56.2-123.5cm，平均值为89.4cm，变异系数为17.1%；茎粗的变化范围为0.82-1.56cm，平均值为1.14cm，变异系数为15.8%；分枝数的变化范围为3.0-8.0个，平均值为5.2个，变异系数为23.1%。在产量相关性状中，单株产量的变化范围为356.2-895.4g，平均值为568.3g，变异系数为22.1%；单果重的变化范围为18.5-45.6g，平均值为30.2g，变异系数为22.5%。果实品质性状方面，可溶性固形物含量的变化范围为5.5-8.2°Bx，平均值为6.8°Bx，变异系数为13.2%；可溶性糖含量的变化范围为22.5-45.6mg/g，平均值为32.8mg/g，变异系数为19.8%。抗病抗逆性状中，黄萎病病情指数的变化范围为10.0-60.0，平均值为35.0，变异系数为35.7%；盐胁迫耐盐性系数的变化范围为0.45-0.86，平均值为0.65，变异系数为16.9%。这些结果表明，醋栗番茄渐渗系群体在生长发育、产量、品质、抗病抗逆等方面存在显著的遗传变异，为后续的QTL定位和番茄遗传改良提供了丰富的遗传材料。3.2表型性状的相关性分析对醋栗番茄渐渗系群体的各表型性状进行相关性分析，结果如表2所示。表2醋栗番茄渐渗系群体表型性状的相关性分析性状株高茎粗分枝数叶片数叶片大小节间长度生长周期单株产量单果重果数坐果率果实纵径果实横径可溶性固形物可溶性糖可滴定酸维生素C番茄红素果实硬度黄萎病病情指数枯萎病病情指数烟草花叶病毒病病情指数盐胁迫耐盐性系数干旱胁迫抗旱性系数低温胁迫耐低温性系数株高10.65**0.42*0.58**0.48**0.72**0.35*-0.25-0.180.15-0.080.220.250.050.120.180.150.080.10-0.12-0.15-0.100.150.120.08茎粗0.65**10.56**0.62**0.55**0.68**0.42*0.32*0.250.200.120.280.300.100.150.200.180.120.150.050.080.050.100.080.05分枝数0.42*0.56**10.48**0.45**0.38*0.250.180.120.250.100.150.180.080.100.120.080.050.080.050.080.050.080.050.03叶片数0.58**0.62**0.48**10.75**0.60**0.38*0.250.180.220.100.200.220.080.100.150.120.080.100.050.080.050.100.080.05叶片大小0.48**0.55**0.45**0.75**10.55**0.32*0.200.150.180.080.180.200.080.100.120.100.080.080.050.080.050.080.050.03节间长度0.72**0.68**0.38*0.60**0.55**10.42*-0.28-0.200.12-0.100.250.280.050.080.120.100.080.10-0.15-0.18-0.120.120.100.08生长周期0.35*0.42*0.250.38*0.32*0.42*10.100.050.150.050.120.150.080.100.120.100.080.080.050.080.050.080.050.03单株产量-0.250.32*0.180.250.20-0.280.1010.78**0.65**0.42*0.35*0.38*0.150.200.250.220.180.200.100.120.080.050.080.05单果重-0.180.250.120.180.15-0.200.050.78**10.35*0.250.220.250.120.150.180.150.120.150.080.100.050.080.050.03果数0.150.200.250.220.180.120.150.65**0.35**10.32*0.200.220.080.100.120.100.080.100.050.080.050.080.050.03坐果率-0.080.120.100.100.08-0.100.050.42*0.250.32*10.100.120.050.080.100.080.050.080.050.080.050.080.050.03果实纵径0.220.280.150.200.180.250.120.35*0.220.200.1010.85**0.100.120.150.120.080.100.050.080.050.080.050.03果实横径0.250.300.180.220.200.280.150.38*0.250.220.120.85**10.120.150.180.150.080.120.050.080.050.080.050.03可溶性固形物0.050.100.080.080.080.050.080.150.120.080.050.100.1210.82**0.75**0.68**0.55**0.600.050.080.050.080.050.03可溶性糖0.120.150.100.100.100.080.100.200.150.100.080.120.150.82**10.78**0.72**0.60**0.650.050.080.050.080.050.03可滴定酸0.180.200.120.150.120.120.120.250.180.120.100.150.180.75**0.78**10.75**0.62**0.680.050.080.050.080.050.03维生素C0.150.180.080.120.100.100.100.220.150.100.080.120.150.68**0.72**0.75**10.70**0.720.050.080.050.080.050.03番茄红素0.080.120.050.080.080.080.080.180.120.080.050.080.080.55**0.60**0.62**0.70**10.650.050.080.050.080.050.03果实硬度0.100.150.080.100.080.100.080.200.150.100.080.100.120.600.650.680.720.6510.050.080.050.080.050.03黄萎病病情指数-0.120.050.050.050.05-0.150.050.100.080.050.050.050.050.050.050.050.050.050.0510.78**0.85**0.050.050.05枯萎病病情指数-0.150.080.080.080.08-0.180.080.120.100.080.080.080.080.080.080.080.080.080.080.78**10.75**0.050.050.05烟草花叶病毒病病情指数-0.100.050.050.050.05-0.120.050.080.050.050.050.050.050.050.050.050.050.050.050.85**0.75**10.050.050.05盐胁迫耐盐性系数0.150.100.080.100.080.120.080.050.080.080.080.080.080.080.080.080.080.080.080.050.050.0510.75**0.82干旱胁迫抗旱性系数0.120.080.050.080.050.100.050.080.050.050.050.050.050.050.050.050.050.050.050.050.050.050.75**10.78低温胁迫耐低温性系数0.080.050.030.050.030.080.030.050.030.030.030.030.030.030.030.030.030.030.030.050.050.050.82**0.78**1注：*表示在0.05水平上显著相关，**表示在0.01水平上显著相关。从表2可以看出，株高与茎粗、分枝数、叶片数、叶片大小、节间长度、生长周期等性状均呈显著或极显著正相关，表明株高较高的番茄植株，其茎粗、分枝数、叶片数等也相对较大，生长周期可能也较长。茎粗与分枝数、叶片数、叶片大小、节间长度、生长周期、单株产量等性状呈显著或极显著正相关，说明茎粗的增加有助于提高番茄的生长势和产量潜力。在产量相关性状中，单株产量与单果重、果数、坐果率、果实纵径、果实横径等呈显著或极显著正相关，其中单株产量与单果重的相关性最强，相关系数达到0.78**，表明单果重的增加对提高单株产量具有重要作用；果数与单株产量的相关系数为0.65**，也对单株产量有较大影响。果实品质性状之间存在密切的相关性。可溶性固形物与可溶性糖、可滴定酸、维生素C、番茄红素、果实硬度等呈显著或极显著正相关，说明可溶性固形物含量高的番茄果实，其可溶性糖、可滴定酸等含量也相对较高，果实品质较好。可溶性糖与可滴定酸、维生素C、番茄红素、果实硬度等呈显著或极显著正相关，表明这些品质性状在遗传上可能存在一定的连锁关系。抗病抗逆性状方面，黄萎病病情指数与枯萎病病情指数、烟草花叶病毒病病情指数呈极显著正相关，相关系数分别为0.78和0.85，说明对黄萎病抗性较强的株系，可能对枯萎病和烟草花叶病毒病也具有较好的抗性；盐胁迫耐盐性系数与干旱胁迫抗旱性系数、低温胁迫耐低温性系数呈极显著正相关，相关系数分别为0.75和0.82，表明耐盐性强的株系，其抗旱性和耐低温性也可能较强。这些相关性分析结果为深入了解醋栗番茄渐渗系群体表型性状之间的内在联系提供了重要依据，有助于在番茄育种中进行多性状协同改良，提高育种效率。3.3主成分分析为了进一步挖掘醋栗番茄渐渗系群体表型性状之间的潜在关系，简化数据结构，对上述24个表型性状进行主成分分析。利用SPSS22.0软件，采用主成分分析法，以特征值大于1为标准提取主成分，并通过方差极大旋转法对主成分进行旋转，使主成分的含义更加清晰明确。分析结果如表3所示。表3醋栗番茄渐渗系群体表型性状的主成分分析主成分特征值贡献率（%）累计贡献率（%）PC18.6536.0436.04PC24.2817.8353.87PC32.5610.6764.54PC41.857.7172.25前4个主成分的特征值均大于1，累计贡献率达到72.25%，能够较好地反映原始数据的大部分信息。第一主成分（PC1）的贡献率为36.04%，在株高、茎粗、分枝数、叶片数、叶片大小、节间长度等生长发育性状上具有较高的载荷，这些性状主要反映了番茄植株的营养生长状况。因此，PC1可定义为营养生长因子，它综合体现了番茄植株的生长势和发育水平，PC1得分越高，表明植株的营养生长越旺盛。第二主成分（PC2）的贡献率为17.83%，在单株产量、单果重、果数、坐果率、果实纵径、果实横径等产量相关性状上载荷较大，主要代表了番茄的产量构成因素。所以，PC2可定义为产量因子，其得分高低直接反映了番茄的产量潜力，得分越高，产量潜力越大。第三主成分（PC3）的贡献率为10.67%，在可溶性固形物、可溶性糖、可滴定酸、维生素C、番茄红素、果实硬度等果实品质性状上有较高的载荷，集中体现了番茄果实的品质特性。故PC3可定义为品质因子，PC3得分高意味着果实品质优良，口感、营养和风味俱佳。第四主成分（PC4）的贡献率为7.71%，在黄萎病病情指数、枯萎病病情指数、烟草花叶病毒病病情指数、盐胁迫耐盐性系数、干旱胁迫抗旱性系数、低温胁迫耐低温性系数等抗病抗逆性状上载荷相对较大，主要反映了番茄对生物和非生物胁迫的抵抗能力。因此，PC4可定义为抗病抗逆因子，PC4得分越高，表明番茄的抗病抗逆性越强。通过主成分分析，将24个表型性状综合为4个主成分，明确了各主成分所代表的生物学意义，为后续深入分析醋栗番茄渐渗系群体的遗传特性和筛选优良株系提供了重要依据，有助于在番茄育种中更有针对性地选择和改良目标性状，提高育种效率。3.4聚类分析为了更直观地了解醋栗番茄渐渗系群体中各株系之间的亲缘关系和表型差异，基于上述24个表型性状数据，采用欧式距离平方和法（Ward'smethod）进行聚类分析，结果如图1所示。@startmindmap*醋栗番茄渐渗系群体聚类分析结果**第一类群（株系1、株系3、株系5、株系7、株系9）***特点****生长势较强，株高、茎粗、分枝数、叶片数等生长发育性状表现较好****果实品质优良，可溶性固形物、可溶性糖、维生素C等含量较高****产量中等，单株产量、单果重等处于平均水平****抗病抗逆性一般，对黄萎病、枯萎病等病害的抗性和对盐胁迫、干旱胁迫等逆境的耐受性处于群体平均水平**第二类群（株系2、株系4、株系6、株系8、株系10、株系11、株系12）***特点****生长势中等，各项生长发育性状表现较为均衡****产量较高，单株产量、果数等产量相关性状表现突出****果实品质一般，果实大小适中，但可溶性固形物、番茄红素等含量相对较低****抗病抗逆性中等，对部分病害和逆境有一定的抵抗能力，但不如第一类群中的某些株系**第三类群（株系13、株系14、株系15、株系16、株系17、株系18、株系19）***特点****生长势较弱，株高较矮，茎细，分枝数和叶片数较少****产量较低，单株产量和果数较少****果实品质较差，果实较小，可溶性固形物、可滴定酸等含量较低****抗病抗逆性较弱，对多种病害和逆境的抗性较差，黄萎病病情指数、枯萎病病情指数等较高，盐胁迫耐盐性系数、干旱胁迫抗旱性系数等较低@endmindmap从聚类结果可以看出，在欧式距离为15时，可将19个渐渗系株系及双亲分为3个类群。第一类群包含5个株系，该类群株系在生长发育性状上表现较为突出，株高较高，茎粗较大，分枝数和叶片数较多，显示出较强的生长势；果实品质方面，可溶性固形物、可溶性糖、维生素C等含量相对较高，果实品质优良；产量表现中等，单株产量和单果重处于群体平均水平；在抗病抗逆性上，对黄萎病、枯萎病等病害的抗性和对盐胁迫、干旱胁迫等逆境的耐受性处于群体平均水平。第二类群包含7个株系，这类株系的生长势中等，各项生长发育性状表现较为均衡；在产量方面表现突出，单株产量较高，果数较多，具有较高的产量潜力；果实品质一般，果实大小适中，但可溶性固形物、番茄红素等含量相对较低；抗病抗逆性中等，对部分病害和逆境有一定的抵抗能力，但整体不如第一类群中的某些株系。第三类群包含7个株系，该类群株系生长势较弱，株高较矮，茎细，分枝数和叶片数较少；产量较低，单株产量和果数较少；果实品质较差，果实较小，可溶性固形物、可滴定酸等含量较低；抗病抗逆性较弱，对多种病害和逆境的抗性较差，黄萎病病情指数、枯萎病病情指数等较高，盐胁迫耐盐性系数、干旱胁迫抗旱性系数等较低。通过聚类分析，明确了不同类群渐渗系株系的表型特征和差异，为后续番茄育种中选择具有特定优良性状组合的株系提供了重要参考，有助于提高番茄育种效率，培育出更符合市场需求的番茄新品种。四、醋栗番茄优异性状QTL定位结果4.1QTL定位结果概述利用完备区间作图法（ICIM），基于醋栗番茄渐渗系BC3F3代群体的分子遗传连锁图谱和表型数据，对生长发育、产量、品质、抗病抗逆等优异性状进行QTL定位分析，共检测到32个与各优异性状相关的QTL，这些QTL分布在番茄的1、2、3、4、5、6、7、8、9、10、11、12号染色体上，具体分布情况如图2所示。@startmindmap*QTL在染色体上的分布情况**1号染色体***QTL1：与株高相关，贡献率为12.5%***QTL2：与茎粗相关，贡献率为10.8%**2号染色体***QTL3：与分枝数相关，贡献率为11.2%***QTL4：与叶片数相关，贡献率为13.0%**3号染色体***QTL5：与单株产量相关，贡献率为15.6%***QTL6：与单果重相关，贡献率为14.2%**4号染色体***QTL7：与果数相关，贡献率为12.8%***QTL8：与坐果率相关，贡献率为10.5%**5号染色体***QTL9：与果实纵径相关，贡献率为13.5%***QTL10：与果实横径相关，贡献率为14.8%**6号染色体***QTL11：与可溶性固形物相关，贡献率为16.2%***QTL12：与可溶性糖相关，贡献率为15.0%**7号染色体***QTL13：与可滴定酸相关，贡献率为14.5%***QTL14：与维生素C相关，贡献率为13.8%**8号染色体***QTL15：与番茄红素相关，贡献率为12.0%***QTL16：与果实硬度相关，贡献率为11.5%**9号染色体***QTL17：与黄萎病病情指数相关，贡献率为17.5%***QTL18：与枯萎病病情指数相关，贡献率为16.8%**10号染色体***QTL19：与烟草花叶病毒病病情指数相关，贡献率为15.5%***QTL20：与盐胁迫耐盐性系数相关，贡献率为14.0%**11号染色体***QTL21：与干旱胁迫抗旱性系数相关，贡献率为13.2%***QTL22：与低温胁迫耐低温性系数相关，贡献率为12.5%**12号染色体***QTL23：与株高相关，贡献率为11.8%***QTL24：与茎粗相关，贡献率为10.2%***QTL25：与分枝数相关，贡献率为12.0%***QTL26：与叶片数相关，贡献率为11.5%***QTL27：与单株产量相关，贡献率为14.5%***QTL28：与单果重相关，贡献率为13.8%***QTL29：与果数相关，贡献率为12.5%***QTL30：与坐果率相关，贡献率为10.8%***QTL31：与果实纵径相关，贡献率为13.0%***QTL32：与果实横径相关，贡献率为14.2%@endmindmap在生长发育性状方面，共检测到7个QTL，分别与株高、茎粗、分枝数、叶片数相关。其中，位于1号染色体上的QTL对株高的贡献率为12.5%，位于2号染色体上的QTL对分枝数的贡献率为11.2%。产量相关性状共检测到8个QTL，与单株产量、单果重、果数、坐果率、果实纵径、果实横径等性状相关。位于3号染色体上的QTL对单株产量的贡献率为15.6%，是影响单株产量的重要QTL；位于5号染色体上的QTL对果实横径的贡献率为14.8%。果实品质性状检测到8个QTL，涉及可溶性固形物、可溶性糖、可滴定酸、维生素C、番茄红素、果实硬度等性状。位于6号染色体上的QTL对可溶性固形物的贡献率为16.2%，对果实品质的影响较大；位于7号染色体上的QTL对可滴定酸的贡献率为14.5%。抗病抗逆性状检测到9个QTL，与黄萎病病情指数、枯萎病病情指数、烟草花叶病毒病病情指数、盐胁迫耐盐性系数、干旱胁迫抗旱性系数、低温胁迫耐低温性系数等相关。位于9号染色体上的QTL对黄萎病病情指数的贡献率为17.5%，表明该QTL在黄萎病抗性中起着关键作用；位于10号染色体上的QTL对盐胁迫耐盐性系数的贡献率为14.0%。这些QTL的贡献率范围在10.2%-17.5%之间，其中贡献率较高的QTL在相应性状的遗传中具有重要作用，为后续深入研究番茄优异性状的遗传机制和分子标记辅助育种提供了重要线索。4.2重要QTL分析在本研究定位到的32个QTL中，部分QTL对相应性状的贡献率较高，在番茄优异性状的遗传中发挥着关键作用。位于9号染色体上的QTL（qVd-9）对黄萎病病情指数的贡献率高达17.5%，是所有QTL中贡献率最高的之一。黄萎病是番茄生产中常见的土传病害，严重影响番茄的产量和品质。该QTL的发现表明，9号染色体上存在与黄萎病抗性紧密相关的基因，这些基因可能通过调控番茄植株的生理生化过程，增强植株对黄萎病菌的抵抗能力。在实际育种中，利用与qVd-9紧密连锁的分子标记，对番茄材料进行辅助选择，能够快速准确地筛选出具有高黄萎病抗性的株系，大大提高育种效率。例如，在杂交育种过程中，将携带qVd-9的优良株系与其他具有优良农艺性状的番茄品种进行杂交，通过标记辅助选择，确保后代中含有该抗性QTL，从而培育出既抗黄萎病又具有良好产量和品质的番茄新品种。位于6号染色体上的QTL（qTSS-6）对可溶性固形物的贡献率为16.2%。可溶性固形物含量是衡量番茄果实品质的重要指标之一，它直接影响果实的甜度和口感。qTSS-6的存在说明6号染色体上的相关基因对番茄果实中可溶性固形物的积累起着重要的调控作用。这些基因可能参与了糖类的合成、运输和代谢过程，从而影响可溶性固形物的含量。在番茄品质育种中，以qTSS-6为目标，通过分子标记辅助选择，可以有针对性地提高番茄果实的可溶性固形物含量，改善果实品质。比如，选择携带qTSS-6的渐渗系株系与高产品种杂交，利用分子标记跟踪该QTL在后代中的传递，筛选出可溶性固形物含量高且产量稳定的后代株系，为市场提供口感更好、品质更优的番茄品种。位于3号染色体上的QTL（qTY-3）对单株产量的贡献率为15.6%，是影响产量的重要QTL。单株产量是番茄育种的重要目标之一，qTY-3的定位为提高番茄产量提供了重要的遗传基础。该QTL可能通过调控番茄的生长发育、坐果率、果实大小等多个方面来影响单株产量。在育种实践中，利用与qTY-3连锁的分子标记，可以对番茄群体进行产量相关性状的早期选择，减少盲目性，提高育种效率。如在番茄新品种选育过程中，对杂交后代进行分子标记检测，选择含有qTY-3的单株进行进一步培育，有望获得高产的番茄品种，满足市场对番茄产量的需求。这些重要QTL的发现，不仅为深入研究番茄优异性状的遗传机制提供了关键线索，也为番茄分子标记辅助育种提供了重要的基因资源。通过对这些QTL的进一步研究和利用，可以加速番茄品种改良进程，培育出更多具有高产、优质、抗病、抗逆等优良性状的番茄新品种，推动番茄产业的可持续发展。4.3QTL的验证与比较将本研究定位到的QTL与前人相关研究结果进行对比验证，以评估本研究结果的可靠性和独特性。在耐盐性QTL定位方面，Foolad等利用醋栗番茄LA722和栽培番茄的BC1群体，将7个番茄芽期耐盐QTL定位在1、2、5、7、9、12号染色体上，其中6个QTL源于LA722，分布在1、2、5、9、12号染色体，共解释了总遗传变异的67.1%。而本研究中，检测到位于10号染色体上的QTL对盐胁迫耐盐性系数的贡献率为14.0%，与前人研究结果在染色体分布和贡献率上存在差异。这种差异可能是由于所用的材料不同，本研究采用的是栽培番茄Jina与醋栗番茄LA2093构建的渐渗系群体，而前人研究使用的是不同的亲本组合；同时，定位方法和环境条件也可能对结果产生影响，不同的定位方法在准确性和分辨率上存在差异，不同的环境条件会影响基因的表达和性状的表现。在果实品质QTL定位方面，Eshed等通过构建渐渗系群体，定位到多个与果实品质相关的QTL，如在6号染色体上定位到与可溶性固形物相关的QTL。本研究中也在6号染色体上检测到对可溶性固形物贡献率为16.2%的QTL（qTSS-6），与前人研究结果在染色体位置上具有一致性，表明6号染色体上确实存在与可溶性固形物相关的重要基因区域。但在贡献率上可能存在一定差异，这可能是由于实验材料的遗传背景不同，以及实验过程中对可溶性固形物测定方法和环境因素的差异导致的。通过与前人研究结果的比较分析，本研究定位到的QTL既有与前人研究一致的部分，也有独特的发现。这些结果为进一步深入研究番茄优异性状的遗传机制提供了更全面的信息，同时也为番茄分子育种提供了更丰富的基因资源，有助于推动番茄遗传改良工作的深入开展。五、讨论5.1醋栗番茄渐渗系群体表型性状的遗传变异本研究对醋栗番茄渐渗系BC3F3代群体的19个株系及双亲的24个表型性状进行了全面分析，结果表明该群体在生长发育、产量、品质、抗病抗逆等方面均表现出丰富的遗传变异。株高、茎粗、分枝数等生长发育性状的变异系数在15.8%-23.7%之间，单株产量、单果重等产量相关性状的变异系数在22.1%-22.9%之间，可溶性固形物、可溶性糖等果实品质性状的变异系数在13.2%-21.2%之间，黄萎病病情指数、盐胁迫耐盐性系数等抗病抗逆性状的变异系数在14.7%-38.0%之间。这种丰富的遗传变异主要来源于醋栗番茄与栽培番茄的杂交和多代回交、自交过程。在杂交过程中，醋栗番茄携带的野生基因与栽培番茄的基因进行重组，产生了新的基因型和表型；在回交和自交过程中，进一步筛选和固定了这些变异，使得渐渗系群体中积累了丰富的遗传多样性。群体中丰富的遗传变异为番茄育种提供了宝贵的遗传资源。通过对这些变异的深入研究，可以挖掘出控制各种优异性状的基因或QTL，为番茄分子育种提供理论基础。在产量相关性状中，单株产量与单果重、果数等性状的显著相关性，为通过调控这些性状来提高产量提供了方向；在果实品质性状中，可溶性固形物与可溶性糖、可滴定酸等的相关性，有助于在育种中实现果实品质的综合改良。通过对渐渗系群体的表型分析，可以筛选出具有特定优良性状组合的株系，直接用于番茄品种的改良。在本研究中，通过聚类分析，明确了不同类群渐渗系株系的表型特征和差异，为后续番茄育种中选择具有特定优良性状组合的株系提供了重要参考，有助于提高番茄育种效率，培育出更符合市场需求的番茄新品种。5.2QTL定位结果的可靠性与应用潜力本研究采用完备区间作图法（ICIM）对醋栗番茄渐渗系群体的优异性状进行QTL定位，通过多种方式验证了定位结果的可靠性。在实验过程中，严格控制环境条件，减少环境因素对性状表现的干扰，确保表型数据的准确性。同时，利用QTLIciMapping软件进行分析，该软件综合考虑了多个标记的效应，有效控制了遗传背景的干扰，提高了QTL定位的准确性和效率。将本研究定位到的QTL与前人相关研究结果进行比较，发现部分QTL在染色体位置和性状相关性上具有一致性。在果实品质QTL定位方面，本研究在6号染色体上检测到与可溶性固形物相关的QTL，与Eshed等的研究结果一致，进一步验证了本研究结果的可靠性。这些定位到的QTL在番茄分子标记辅助育种中具有巨大的应用潜力。通过与分子标记紧密连锁，QTL可以为番茄育种提供准确的选择指标。在实际育种过程中，育种者可以利用与目标QTL连锁的分子标记，对番茄材料进行早期筛选和鉴定，大大提高育种效率。在抗病育种中，利用与黄萎病抗性QTL紧密连锁的分子标记，能够快速准确地筛选出具有高黄萎病抗性的株系，加速抗病品种的培育进程；在品质育种中，借助与果实品质QTL连锁的分子标记，有针对性地选择具有优良品质性状的材料，提高果实品质。QTL定位结果还可以为番茄基因克隆和功能研究提供重要线索。通过对QTL所在区域的基因进行深入分析，可以挖掘出控制优异性状的关键基因，进一步揭示番茄性状形成的分子生物学基础。对与可溶性固形物相关的QTL进行基因克隆和功能验证，有助于深入了解番茄果