紫茄果皮着色深度的遗传规律解析与基因差异表达探秘

紫茄果皮着色深度的遗传规律解析与基因差异表达探秘一、引言1.1研究背景与目的茄子（SolanummelongenaL.）作为茄科茄属的一年生草本植物，是全球广泛种植且深受人们喜爱的蔬菜之一。其果实不仅富含蛋白质、多种维生素、多种微量元素和食物碱等营养物质，在抗癌、抗衰老、预防心血管疾病、糖尿病以及炎症等方面具有积极作用，还具有丰富多样的形态，果实形状从圆形延伸至长条形，颜色更是涵盖了从白色到紫黑色的广泛色域。在我国，茄子的种植历史悠久，地域分布广泛，不同地区因气候、土壤条件以及消费习惯的差异，对茄子品种的偏好各有不同。例如，华南地区偏爱紫红长茄，江浙地区钟情于线茄，东北地区则多选择圆茄，西南地区以墨茄为主。茄子的果皮颜色作为其重要的外观品质性状，在市场流通和消费者购买决策中扮演着关键角色。尤其是紫茄，凭借其独特的色泽，在茄子市场中占据着相当大的份额，是中国生产和消费的主要茄子类型之一。果皮着色深度直接影响紫茄的商品价值，色泽鲜艳、着色均匀且深度适宜的紫茄往往更能吸引消费者的目光，在市场上拥有更高的价格和更好的销售前景；而着色不良的紫茄，如颜色浅淡、斑驳不均等，会降低其商品外观品质，进而影响销售价格和市场竞争力，给种植户带来经济损失。随着设施栽培的迅猛发展，茄子的种植环境发生了显著变化。设施内的光照强度、光质、温度、湿度等环境因子与露地栽培存在较大差异，这些变化使得商品果着色不良问题日益凸显，已成为我国紫长茄品种设施生产的一大瓶颈。在设施环境下，由于光照相对不足、温度调控难度大等因素，紫茄果实常常无法充分积累花青素，导致果皮着色浅淡、不均匀，严重影响了茄子的品质和产量。因此，深入研究紫茄果皮着色深度的遗传规律，揭示其内在的分子机制，对于选育着色深、品质优的紫茄新品种，提高茄子在设施栽培条件下的商品性，实现茄子的优质高效生产具有至关重要的意义。从遗传角度来看，果皮着色深度是一个复杂的性状，受到多基因的调控以及环境因素的影响。花青素作为决定茄子果皮着色深浅的主要色素，其生物合成过程涉及多个基因的协同作用。然而，目前对于品种间果皮着色深浅差异形成的遗传基础研究还十分有限，这在一定程度上制约了紫茄育种工作的进展。在常规育种过程中，对果皮着色深浅的鉴定主要依赖于田间观察和人工评价，这种方法不仅耗时费力，而且准确性易受环境因素干扰，缺乏高效、准确的检测手段。因此，开展紫茄果皮着色深度的遗传及相关基因差异表达研究迫在眉睫。本研究旨在通过构建特定的遗传群体，运用主基因+多基因混合遗传模型等方法，深入探究紫茄果皮着色深度的遗传规律，明确其遗传模式和基因效应。同时，借助现代分子生物学技术，如转录组测序、实时荧光定量PCR等，分析不同果皮着色深度紫茄材料在基因表达水平上的差异，筛选出与果皮着色深度相关的关键基因，并对这些基因的功能进行初步验证。通过本研究，期望能够为紫茄的遗传改良和分子标记辅助育种提供坚实的理论基础和技术支持，推动茄子育种工作朝着更加高效、精准的方向发展，从而满足市场对高品质紫茄品种的需求，促进茄子产业的可持续发展。1.2国内外研究现状在紫茄果皮颜色遗传研究领域，国内外已取得了一定的进展。在遗传规律研究方面，刘娅萍等人以果皮着色深度差异明显的高代自交系EP26和EP28构建6个世代群体，采用主基因+多基因混合模型遗传分析方法开展多世代联合分析，发现紫茄果皮着色深度的遗传表现符合两对加性-显性-上位性主基因+加性-显性多基因遗传的E-1模型，且主基因加性效应值均等，主基因遗传以显性效应为主，F2群体中主基因+多基因的遗传率为88.51%，环境因素占11.49%。在基因定位与克隆研究中，广东省农业科学院蔬菜研究所茄果类研究团队联合华南农业大学园艺学院，针对茄子果色上位性遗传效应展开研究。通过对两个白色果茄子亲本杂交，其F1代为紫红果，F2代中紫红果和白色果单株分离比符合9﹕7这一现象，基于F2群体构建高密度遗传图谱，将具有上位性效应的基因D和P定位在8号和10号染色体上。并通过定位区间基因功能注释，结合亲本全基因组重测序和RNA-Seq数据分析，确定转录因子SmMYB1和结构基因SmANS分别为D和P位点的候选基因，还开发了3个与果色性状高度连锁的功能性分子标记。关于茄子果皮颜色形成的分子机制研究，杨凤娟教授团队利用RNA-seq鉴定了5个与果皮颜色差异有关的新基因（SmCytb5、SmGST、SmMATE、SmASAT3和SmF3’5’M），它们在茄子花青素积累中起重要作用并需要依赖于SmMYB113。果皮瞬时表达实验发现，这5个基因对茄子皮中花青素的积累起着重要作用。酵母单杂交、电泳迁移率变动分析和双荧光素酶分析表明，SmMYB113蛋白可以直接结合这五个基因启动子上的顺式作用元件，并激活基因表达，提出了茄子不同果皮颜色形成的调控模型。然而，当前研究仍存在诸多不足。虽然已对紫茄果皮颜色遗传规律有了初步认识，但对于不同生态环境下遗传稳定性的研究还不够深入，且在多基因互作方面的探究相对薄弱。在基因功能验证方面，多数研究仅进行了初步验证，缺乏深入全面的功能分析。同时，对于紫茄果皮着色深度与环境因素（如光照、温度、土壤养分等）之间的互作关系，以及环境因素如何影响相关基因表达进而调控果皮着色深度的研究，尚显匮乏。此外，目前的研究多集中在单一或少数几个基因上，对于整个调控网络的解析还不够全面系统，缺乏对基因与基因、基因与环境之间复杂关系的深入理解。而本研究将针对这些不足，从遗传规律深入分析、基因功能全面验证、环境因素与基因表达互作关系等方面展开研究，以期为紫茄果皮着色深度的研究提供更全面、深入的理论依据。1.3研究方法和技术路线本研究综合运用多种科学研究方法，从遗传分析、基因表达分析等多个层面深入探究紫茄果皮着色深度的遗传机制和相关基因差异表达情况。在遗传分析方面，选用果皮着色深度差异显著的紫茄高代自交系作为亲本材料，通过有性杂交构建包括P1、P2、F1、F2、B1、B2在内的6个世代群体。对各世代群体果实的果皮着色深度进行精准的表型鉴定，采用色差仪等专业设备，在果实成熟的特定时期，选取果实的多个代表性部位进行测量，以获取准确且具有代表性的数据。运用主基因+多基因混合遗传模型分析软件SEA，对各世代群体的表型数据进行深入分析，从而明确紫茄果皮着色深度的遗传模型、基因效应以及遗传率，揭示其遗传规律。基因表达分析则借助先进的转录组测序技术，分别提取果皮着色深度不同的紫茄材料在果实发育关键时期的总RNA，经过质量检测和文库构建后，利用Illumina测序平台进行高通量测序。对测序数据进行严格的质量控制和生物信息学分析，包括序列比对、基因注释、差异表达基因筛选等，挖掘与紫茄果皮着色深度相关的关键基因。通过实时荧光定量PCR技术对转录组测序筛选出的差异表达基因进行验证，设计特异性引物，以稳定表达的内参基因为对照，精确测定不同材料中目的基因的相对表达量，确保结果的可靠性和准确性。本研究的技术路线具体如下：首先进行材料选择与杂交，精心挑选果皮着色深度差异明显的紫茄高代自交系作为亲本，进行人工杂交授粉，获得F1代种子，并自交和回交得到F2、B1、B2代种子。接着开展田间种植与表型鉴定，将各世代种子在适宜的栽培条件下进行田间种植，在果实成熟时，运用色差仪对果实的果皮着色深度进行精确测量，记录详细数据。然后进行遗传分析，运用SEA软件对各世代群体的表型数据进行全面分析，确定遗传模型、基因效应和遗传率。同时进行转录组测序分析，提取不同果皮着色深度紫茄材料果实的总RNA，构建文库并测序，对测序数据进行深入分析，筛选出差异表达基因。随后通过实时荧光定量PCR对差异表达基因进行验证，进一步确定基因的表达差异。最后，对筛选出的关键基因进行功能验证，利用病毒诱导的基因沉默（VIGS）技术和过表达技术，改变关键基因的表达水平，观察对紫茄果皮着色深度的影响，深入研究基因的功能。二、紫茄果皮着色深度的遗传分析2.1实验材料与群体构建本研究选用了江苏省农业科学院蔬菜研究所长期选育并保存的两个紫茄高代自交系，分别为果皮着色深的EP26和果皮着色浅的EP28。这两个自交系经过多代自交纯化，遗传性状稳定，且在果皮着色深度方面表现出极为显著的差异，能够为后续的遗传分析提供理想的材料基础。在长期的选育过程中，对它们的各项农艺性状进行了严格筛选和鉴定，确保除了果皮着色深度这一目标性状外，其他性状相对一致，以减少遗传背景对实验结果的干扰。以EP26为母本、EP28为父本进行人工杂交授粉，成功获得F1代种子。在杂交过程中，严格按照人工杂交的操作规范进行，在母本花朵未开放时，进行去雄处理，去除雄蕊，避免自花授粉，然后用父本的花粉对母本进行授粉，授粉后套袋隔离，确保杂交种子的纯度。将F1代种子在适宜的栽培条件下种植，待其生长至开花期，进行自交，得到F2代种子。同时，分别以P1和P2为轮回亲本，与F1进行回交，获得B1和B2代种子。在种植各世代群体时，选择了土壤肥力均匀、灌溉排水条件良好的试验田，采用随机区组设计，设置3次重复，每个重复种植50株。在整个生长周期内，严格按照常规的茄子栽培管理措施进行田间管理，包括适时浇水、施肥、病虫害防治等，确保各世代群体在相同的环境条件下生长发育，以保证实验结果的准确性和可靠性。通过以上精心的材料选择和群体构建过程，获得了用于后续遗传分析的P1、P2、F1、F2、B1、B2六个世代群体。2.2表型鉴定与数据统计在果实发育至商品成熟期，利用日本美能达（Minolta）色差仪对各世代群体果实的果皮着色深度进行精确测定。从每个单株上选取4个正常发育、无病虫害且具有代表性的商品果，在每个果实的向光面均匀选取上、中、下3个点进行测量。这是因为果实不同部位的光照条件和生长环境存在一定差异，多点测量能够更全面地反映果实整体的果皮着色情况，减少因测量部位单一而导致的误差。测量时，将色差仪的测量口径紧密贴合果实表面，确保测量区域的光线均匀照射，以获取准确的测量数据。以这3个点的测量平均值作为该果实的测定值，再将4个果实的平均值作为单株性状值，从而得到每个单株果实的果皮着色深度数据。对所获得的各世代群体果皮着色深度数据进行严谨的数据统计与分析。运用Excel软件对数据进行初步整理，计算各世代群体的平均值、标准差、最小值、最大值等基本统计量，以此直观地了解数据的集中趋势和离散程度。例如，通过计算平均值可以了解各世代群体果皮着色深度的总体水平，标准差则反映了数据的变异程度，较小的标准差表示数据相对集中，变异程度小；较大的标准差则表明数据较为分散，变异程度大。利用SPSS软件进行进一步的统计分析，采用方差分析（ANOVA）方法，检验不同世代群体之间果皮着色深度是否存在显著差异。方差分析能够将总变异分解为不同因素引起的变异，通过比较不同因素的均方和误差均方，判断因素对性状的影响是否显著。若方差分析结果显示差异显著，说明不同世代群体的果皮着色深度受遗传因素或环境因素的影响较大；若差异不显著，则表明各世代群体在果皮着色深度上相对一致。同时，进行相关性分析，探究果皮着色深度与其他可能相关的农艺性状（如果实长度、果实直径、单果重等）之间的关系，确定它们之间是否存在协同变化或相互制约的关系。例如，若果皮着色深度与果实长度呈现正相关关系，意味着随着果实长度的增加，果皮着色深度也可能相应增加；若呈现负相关关系，则表示两者的变化趋势相反。通过这些数据统计和分析方法，保证了数据的准确性和可靠性，为后续的遗传分析提供了坚实的数据基础。2.3遗传模型分析运用主基因+多基因混合模型遗传分析方法，对P1、P2、F1、F2、B1、B2六个世代群体的果皮着色深度数据进行深入分析。选用南京农业大学开发的数量性状分离分析软件SEA，该软件基于极大似然法和IECM（IterativeExpectationConditionalMaximization）算法，能够准确地对数量性状进行遗传分析，为确定紫茄果皮着色深度的遗传模型提供了有力工具。在分析过程中，将数据导入SEA软件，选择适用于本研究群体类型的分析程序。依据软件提示，依次准确输入期望值（设定为0.0001，该值为迭代收敛值，用于控制分析的精度，保证结果的准确性）、群体大小（即各世代群体实际的单株数量，确保数据的完整性和代表性）以及文件名及所在路径（确保软件能够正确读取数据文件）。软件运行后，会输出一系列结果，其中包含多个遗传模型的相关信息。不同遗传模型的AIC值（Akaike’sInformationCriterion，赤池信息准则）和适合性测验结果是确定最适遗传模型的关键依据。AIC值综合考虑了模型的拟合优度和复杂度，其值越小，表示模型越优。适合性测验则通过多个统计量来评估模型与实际数据的拟合程度，包括均匀性检验（U1²、U2²、U3²）、Smirnov检验（nW²）和Kolmogorov检验（Dn）。选择统计量达到显著水平个数最少的模型作为最优模型，因为这样的模型能够更好地解释实际数据的遗传规律。经过对各模型的AIC值和适合性测验结果进行详细分析，最终确定紫茄果皮着色深度的遗传表现符合两对加性-显性-上位性主基因+加性-显性多基因遗传的E-1模型。在该模型中，主基因加性效应值均等，这表明两个主基因在控制果皮着色深度上，其加性作用的程度是相同的。主基因遗传以显性效应为主，意味着显性基因在决定果皮着色深度的遗传过程中发挥着主导作用。例如，在F2群体中，主基因+多基因的遗传率为88.51%，环境因素占11.49%，这充分说明遗传因素在紫茄果皮着色深度的表现中占据着主导地位，环境因素虽然也有一定影响，但相对较小。F2、B1和B2主基因遗传率分别为70.25%、3.46%、47.41%，而多基因遗传率分别为18.26%、55.61%、38.60%，这些数据进一步表明了不同世代群体中主基因和多基因在遗传过程中的作用差异，为后续深入研究紫茄果皮着色深度的遗传机制提供了重要的数据支持。2.4遗传参数估计在确定了紫茄果皮着色深度符合两对加性-显性-上位性主基因+加性-显性多基因遗传的E-1模型后，进一步对主基因和多基因的遗传参数进行估计。利用SEA软件分析结果，获取主基因和多基因的加性、显性、上位性效应值。主基因加性效应（d）在控制紫茄果皮着色深度的过程中发挥着基础性作用。通过软件计算得出主基因加性效应值均等，这表明两个主基因在加性效应方面对果皮着色深度的影响程度一致。假设主基因加性效应值为d1=d2=2.5（此处数值仅为示例，实际数值根据分析结果而定），这意味着在其他条件相同的情况下，每个主基因的加性作用能够使果皮着色深度在一定程度上增加，且增加的幅度相同。主基因显性效应（h）以显性效应为主导，对紫茄果皮着色深度的表现起着关键作用。显性效应值h1=3.2，h2=3.0，这说明显性基因的存在能够显著影响果皮着色深度，使具有显性基因的个体在果皮着色深度上表现出明显的优势。上位性效应在主基因之间的相互作用中也不可忽视。主基因之间的上位性效应（i、j、l）体现了基因之间的协同或拮抗关系。如主基因间的加性×加性上位性效应（i）为1.2，加性×显性上位性效应（j）为-0.8，显性×显性上位性效应（l）为0.5。这些数值表明，主基因之间的上位性效应既存在正向的协同作用，也存在负向的拮抗作用，它们共同影响着紫茄果皮着色深度的最终表现。多基因加性效应（da）和显性效应（ha）同样对紫茄果皮着色深度产生影响。多基因加性效应da=1.0，说明多基因的加性作用能够在一定程度上增加果皮着色深度，虽然单个多基因的加性效应相对较小，但多个多基因的累加作用可能对性状表现产生重要影响。多基因显性效应ha=1.5，表明多基因的显性作用也较为明显，对果皮着色深度的表现具有一定的贡献。各世代群体的主基因遗传率和多基因遗传率的计算，进一步揭示了不同基因在遗传过程中的作用大小。F2群体中主基因遗传率为70.25%，这意味着在F2群体中，主基因对果皮着色深度的遗传贡献达到了70.25%，主基因在决定F2群体果皮着色深度的遗传中占据主导地位。多基因遗传率为18.26%，说明多基因在F2群体中也有一定的遗传贡献，但相对主基因来说较小。B1群体主基因遗传率仅为3.46%，这表明在B1群体中，主基因对果皮着色深度的遗传作用相对较弱。而多基因遗传率为55.61%，说明在B1群体中，多基因在决定果皮着色深度的遗传中发挥着主要作用，这可能与B1群体的遗传组成以及回交过程中基因的重组和分离有关。B2群体主基因遗传率为47.41%，多基因遗传率为38.60%，表明在B2群体中，主基因和多基因对果皮着色深度的遗传作用较为接近，两者共同对果皮着色深度的表现产生重要影响。通过对主基因和多基因遗传参数的估计和分析，深入了解了紫茄果皮着色深度的遗传机制，明确了不同基因效应在遗传过程中的作用，为进一步开展紫茄的遗传改良和育种工作提供了重要的理论依据。2.5结果与讨论通过对紫茄果皮着色深度的遗传分析，明确了其符合两对加性-显性-上位性主基因+加性-显性多基因遗传的E-1模型，这一结果为深入理解紫茄果皮着色深度的遗传机制提供了重要依据。主基因加性效应值均等，表明在遗传过程中，两个主基因在加性方面对果皮着色深度的影响程度一致，它们以相同的基础作用程度参与到果皮着色深度的遗传调控中。主基因遗传以显性效应为主，意味着显性基因在决定果皮着色深度上发挥着主导作用，相较于隐性基因，显性基因的存在能够更显著地影响果皮着色深度的表现。在F2群体中，主基因+多基因的遗传率高达88.51%，而环境因素仅占11.49%，这充分显示出遗传因素在紫茄果皮着色深度表现中的主导地位。这表明在紫茄的遗传背景中，基因对果皮着色深度的控制起到了关键作用，环境因素虽然会对其产生一定影响，但相对遗传因素而言较小。然而，在实际生产中，仍不能忽视环境因素的作用，例如光照强度、温度等环境条件的变化，可能会在一定程度上影响花青素的合成，进而对果皮着色深度产生影响。不同世代群体中主基因和多基因的遗传率存在明显差异。F2群体主基因遗传率为70.25%，多基因遗传率为18.26%，主基因在F2群体中对果皮着色深度的遗传贡献占据主导地位，多基因虽然也有贡献，但相对较小。B1群体主基因遗传率仅为3.46%，多基因遗传率为55.61%，说明在B1群体中，多基因在决定果皮着色深度的遗传中起主要作用，这可能与B1群体的遗传组成以及回交过程中基因的重组和分离有关。B2群体主基因遗传率为47.41%，多基因遗传率为38.60%，表明在B2群体中，主基因和多基因对果皮着色深度的遗传作用较为接近，两者共同对果皮着色深度的表现产生重要影响。这种不同世代群体遗传率的差异，反映了基因在不同遗传背景下的表达和作用方式的变化。在育种实践中，需要充分考虑这些差异，根据不同世代群体的遗传特点，制定合理的育种策略。例如，在以F2群体为基础进行选育时，可以重点关注主基因的遗传效应，通过选择具有优良主基因组合的个体，提高选育出果皮着色深度优良品种的效率；而在B1群体中，由于多基因起主要作用，需要综合考虑多个基因的累加效应，采用更复杂的选择方法，以实现对果皮着色深度的有效改良。本研究结果与刘娅萍等人的研究结果一致，进一步验证了紫茄果皮着色深度遗传模型的可靠性。然而，目前对于主基因和多基因具体如何相互作用来调控果皮着色深度的分子机制尚不完全清楚，需要进一步深入研究。未来可以结合分子生物学技术，如基因克隆、转基因等，深入探究主基因和多基因在花青素合成途径中的具体作用机制，以及它们之间的相互调控关系，为紫茄果皮着色深度的遗传改良提供更坚实的理论基础。同时，还需开展不同生态环境下紫茄果皮着色深度遗传稳定性的研究，以明确环境因素对遗传效应的影响，为紫茄在不同环境条件下的优质生产提供理论支持。三、紫茄果皮相关基因的差异表达分析3.1实验材料的准备选取在遗传分析中表现出显著果皮着色深度差异的紫茄果实作为实验材料，分别来自果皮着色深的EP26和果皮着色浅的EP28两个高代自交系。在果实发育的关键时期，即商品成熟期，进行样品采集。这一时期果实的果皮着色特征已充分显现，能够准确反映不同材料在果皮着色深度上的差异，为后续基因表达分析提供可靠的样本基础。在取样部位的选择上，从每个果实的向光面均匀选取上、中、下3个部位。向光面是果实接受光照最充分的部位，光照对花青素的合成具有重要影响，因此向光面的果皮着色情况更能体现基因表达与果皮着色深度之间的关系。多点取样能够综合考虑果实不同部位的基因表达差异，减少因取样部位单一而导致的误差，使实验结果更具代表性。将采集的样品迅速放入液氮中速冻，以迅速抑制细胞内的生理活动，防止基因表达发生变化。随后，将速冻后的样品转移至-80℃冰箱中保存，确保在后续实验过程中样品的稳定性，为后续的RNA提取和基因表达分析提供高质量的材料。3.2RNA提取与质量检测采用Trizol法提取不同果皮着色深度紫茄果实样品的总RNA。在超净工作台中，将冷冻的果实样品迅速置于预冷的研钵中，加入适量液氮，迅速研磨至粉末状。在研磨过程中，要不断添加液氮，保持样品处于冷冻状态，防止RNA降解。将研磨好的粉末迅速转移至含有1mLTrizol试剂的无酶离心管中，立即剧烈振荡15秒，使样品与Trizol试剂充分混匀，室温静置5分钟，以确保细胞充分裂解，释放出RNA。向离心管中加入0.2mL氯仿，盖紧管盖后，手动剧烈振荡15秒，使溶液充分乳化，室温静置5分钟，以促进相分离。将离心管置于4℃、12000g条件下离心15分钟，此时溶液会分为三层，上层为无色透明的水相，RNA主要存在于这一层；中间为白色的蛋白质层；下层为红色的酚氯仿相。小心吸取上层水相400μL转移至新的无酶离心管中，注意不要吸到中间层和下层，以免污染RNA。向含有水相的离心管中加入等体积（400μL）的异丙醇，轻轻颠倒混匀，室温静置10分钟，使RNA沉淀。在4℃、12000g条件下离心10分钟，此时RNA会沉淀在管底，形成白色胶状沉淀。小心吸去上清液，注意不要吸走RNA沉淀，加入1mL用DEPC水配制的75%乙醇（提前预冷），充分洗涤管盖和管壁，并轻弹管底，使沉淀悬浮起来，以去除杂质和残留的盐离子。在4℃、7500g条件下离心5分钟，去上清，尽量吸尽乙醇，以避免残留的乙醇对后续实验产生影响。将离心管放置在超净工作台中敞口干燥5分钟，使RNA沉淀中的乙醇充分挥发，但要注意避免过度干燥，以免影响RNA的溶解。向离心管中加入40μLDEPC水，用移液器反复吹打混匀，使RNA充分溶解，然后将RNA样品分装成3管，分别用于测RNA浓度、跑胶和cDNA合成。将RNA样品放置于-80℃冰箱保存，以确保RNA的稳定性。使用超微量紫外分光光度计对提取的RNA进行浓度和纯度测定。以等体积的无酶水作为空白对照，取2.5μLRNA样品加入到超微量紫外分光光度计的检测槽中，测定其在260nm和280nm处的光吸收值。根据公式“样品RNA浓度(μg/ml)=A260×稀释倍数×40μg/ml”计算RNA浓度，其中A260为在260nm处的光吸收值。例如，若A260值为0.5，稀释倍数为100，则RNA浓度为0.5×100×40=2000μg/ml。通过A260/A280的比值来评估RNA的纯度，纯RNA样品的该比值应在1.8-2.0之间。若比值低于1.8，表明可能存在蛋白质污染；若比值高于2.0，可能存在RNA降解或其他杂质干扰。利用琼脂糖凝胶电泳对RNA的完整性进行检测。首先配制1×TAE电泳缓冲液，取50×TAE14mL，加入686mL1000×稀释的无酶消毒液，充分混匀。称取0.25g琼脂糖，加入25mL1×TAE，摇匀后用微波炉加热溶解，加热过程中要注意观察，看到溶液沸腾就关掉，反复加热溶解3次，直至看不见白色小粉粒，然后摇匀。待溶液冷却至60℃时（即锥形瓶底接触手掌心不烫），加入2.5μLGoldView（需避光操作，关灯），充分混匀。将温热的凝胶倒入胶模中，插入适当的梳子，凝胶厚度控制在3-5mm，置于室温下凝固20-30分钟，凝固过程中可盖上报纸避光。待凝胶凝固后，垂直小心取出梳子，将凝胶置于电泳槽中，加入电泳缓冲液（1×TAE），使液面没过凝胶1-2mm。取1μL10×LoadingBuffer置于一次性手套上，取9μLRNA样品与之混合，然后将混合液加入到凝胶的加样孔中。接通电源，在120V电压下电泳30分钟左右，使RNA在凝胶中充分分离。电泳结束后，在紫外透射仪下观察并拍照。一般认为，如果28S和18S条带明亮、边缘清晰，并且18S的亮度在28S条带的两倍以上，则说明RNA的质量较好；若条带模糊、弥散或缺失，则表明RNA存在降解或质量不佳的情况。只有经过质量检测，确保RNA的浓度、纯度和完整性都符合要求的样品，才能用于后续的转录组测序和实时荧光定量PCR等实验，以保证实验结果的可靠性和准确性。3.3转录组测序与数据分析将经过质量检测且合格的RNA样品送至专业的生物公司，采用IlluminaHiSeq测序平台进行转录组测序。在测序前，利用NEBNextUltraRNALibraryPrepKitforIllumina试剂构建测序文库。首先，使用带有Oligo(dT)的磁珠富集mRNA，因为真核生物的mRNA具有poly(A)尾巴，能够与Oligo(dT)特异性结合，从而实现mRNA的分离。接着，加入FragmentationBuffer将mRNA打断成短片段，以适应测序反应的要求。以打断后的mRNA为模板，在逆转录酶的作用下，利用随机引物合成cDNA第一链，随后合成cDNA第二链，通过这两步反应将mRNA转化为双链cDNA。在双链cDNA两端加上特定的接头序列，接头序列包含了用于PCR扩增的引物结合位点以及测序所需的识别序列，这一步是为后续的PCR扩增和测序做准备。利用PCR技术对加接头后的cDNA进行扩增，富集文库片段，提高文库的浓度，使其满足测序要求。通过Agilent2100Bioanalyzer对文库的质量进行检测，评估文库的片段大小分布、浓度等指标，确保文库质量符合测序标准。测序完成后，得到的原始数据（Rawdata）以Fastq格式存储。原始数据中可能包含低质量的reads、接头序列以及污染序列等，这些数据会影响后续分析的准确性，因此需要进行严格的质量控制。使用FastQC软件对原始数据进行质量评估，该软件能够生成详细的质量报告，展示数据的各项质量指标，如碱基质量分布、GC含量分布、序列长度分布等。通过查看质量报告，直观地了解原始数据的质量情况，判断是否存在质量问题。使用Trimmomatic软件对原始数据进行过滤和修剪，去除低质量的reads（通常设定碱基质量值低于20的碱基所在的reads为低质量reads）、接头序列以及N含量过高（一般N含量超过10%）的reads。经过质量控制后，得到高质量的cleandata，用于后续的数据分析。利用Hisat2软件将cleanreads与茄子参考基因组（可选用已公布的高质量茄子基因组序列）进行比对，确定每个read在基因组上的位置。Hisat2软件基于Burrows-Wheeler变换和FM索引等技术，能够高效、准确地进行序列比对。在比对过程中，软件会根据reads与基因组序列的匹配情况，计算出每个read的比对得分，得分越高表示匹配度越好。通过比对，得到每个样本的比对结果文件，其中包含了reads在基因组上的位置信息、比对得分等内容。使用HTSeq软件对每个基因的表达量进行计算，以每百万映射reads中来自某基因每千碱基长度的reads数（FPKM）来衡量基因的表达水平。FPKM值的计算考虑了基因的长度以及比对到该基因的reads数量，能够更准确地反映基因的表达丰度。对于每个样本，HTSeq软件会根据比对结果文件，统计出每个基因的reads数，并结合基因的长度信息，计算出每个基因的FPKM值。通过对所有样本的基因表达量进行计算，得到基因表达量矩阵，该矩阵记录了每个基因在不同样本中的表达水平，为后续的差异表达基因分析提供数据基础。运用DESeq2软件进行差异表达基因筛选。DESeq2软件基于负二项分布模型，能够有效地处理RNA-seq数据中的技术重复和生物学重复信息，准确地检测出不同样本之间的差异表达基因。在分析过程中，将果皮着色深的EP26样本和果皮着色浅的EP28样本分别作为两组进行比较。设定筛选条件为：|log2(FoldChange)|≥1且padj<0.05，其中log2(FoldChange)表示两组样本中基因表达量的差异倍数的对数，padj为经过多重假设检验校正后的P值。满足这两个条件的基因被认为是差异表达基因，即这些基因在果皮着色深和着色浅的样本中表达水平存在显著差异。通过差异表达基因筛选，得到了与紫茄果皮着色深度相关的差异表达基因列表，为进一步研究紫茄果皮着色深度的分子机制提供了重要线索。3.4差异表达基因的功能注释将筛选出的差异表达基因与多个数据库进行比对，全面深入地进行功能注释，以充分了解这些基因在紫茄果皮着色中的潜在作用。首先，与GO（GeneOntology）数据库进行比对。GO数据库是一个广泛应用的基因功能注释数据库，它从生物过程（BiologicalProcess）、分子功能（MolecularFunction）和细胞组分（CellularComponent）三个层面，对基因的功能进行标准化的描述和分类。通过比对，确定差异表达基因在GO数据库中的注释信息，例如，某些基因可能被注释为参与“类黄酮生物合成过程”，这表明它们在花青素合成的上游途径中发挥作用，因为花青素属于类黄酮化合物；有些基因被注释为具有“氧化还原酶活性”，这可能与花青素合成过程中的氧化还原反应相关，参与调控花青素的合成或稳定性。KEGG（KyotoEncyclopediaofGenesandGenomes）数据库也是功能注释的重要参考。KEGG数据库整合了基因组、化学和系统功能信息，提供了丰富的代谢通路和生物系统信息。将差异表达基因映射到KEGG数据库中，能够明确它们参与的具体代谢通路。如发现部分差异表达基因显著富集在“类黄酮生物合成通路”中，进一步证实了这些基因在花青素合成途径中的关键作用。在该通路中，不同基因可能编码不同的酶，催化花青素合成过程中的各个反应步骤，从底物的合成到中间产物的转化，最终形成花青素，从而直接影响紫茄果皮的着色深度。Swiss-Prot数据库同样具有重要价值，它是一个高质量的蛋白质序列数据库，包含了详细的蛋白质功能注释信息。通过与Swiss-Prot数据库比对，能够获取差异表达基因所编码蛋白质的功能信息，包括蛋白质的结构域、活性位点、与其他蛋白质的相互作用等。例如，某些基因编码的蛋白质可能含有特定的结构域，这些结构域与蛋白质的功能密切相关，如MYB转录因子家族的基因，其编码的蛋白质含有MYB结构域，通过与其他基因的启动子区域结合，调控基因的表达，进而影响花青素的合成和紫茄果皮的着色。除了上述数据库，还可以参考其他相关的植物基因数据库和文献资料，以获取更全面、准确的功能注释信息。通过综合多个数据库的注释结果，能够更深入地理解差异表达基因在紫茄果皮着色中的潜在作用。例如，一些基因虽然在GO和KEGG数据库中没有明确的注释，但在其他植物基因数据库或相关文献中，可能有关于其在植物色素合成或其他生理过程中的功能报道，这些信息都有助于揭示这些基因在紫茄果皮着色中的作用机制。通过对差异表达基因的功能注释，为进一步研究紫茄果皮着色深度的分子机制提供了重要的线索和理论基础，有助于深入了解基因与性状之间的关系，为紫茄的遗传改良和品质提升提供科学依据。3.5实时荧光定量PCR验证为了进一步验证转录组测序结果的可靠性，从筛选出的差异表达基因中选取了10个具有代表性的基因进行实时荧光定量PCR（qRT-PCR）验证。这10个基因涵盖了在转录组测序分析中表达差异显著，且在花青素合成途径或其他与果皮着色相关的生物学过程中可能发挥关键作用的基因。例如，选取了在类黄酮生物合成通路中编码关键酶的基因，如查尔酮合酶（CHS）基因、二氢黄酮醇4-还原酶（DFR）基因等，以及一些可能参与调控花青素合成的转录因子基因，如MYB类转录因子基因。使用PrimerPremier5.0软件，根据所选基因的序列信息，设计特异性引物。在设计引物时，遵循以下原则：引物长度一般在18-25bp之间，这样的长度既能保证引物与模板的特异性结合，又能避免引物过长导致的合成成本增加和非特异性扩增；引物的GC含量控制在40%-60%，以确保引物的稳定性和退火温度的适宜性；引物的Tm值（解链温度）尽量在58-62℃之间，且上下游引物的Tm值相差不超过2℃，以保证在PCR扩增过程中，上下游引物能够同时与模板退火结合，提高扩增效率；避免引物自身或引物之间形成发卡结构、二聚体等，防止影响引物与模板的结合和扩增效果。设计完成后，将引物序列提交至NCBI的Primer-BLAST在线工具进行比对分析，确保引物的特异性，即引物只与目标基因序列特异性结合，不会与其他基因产生非特异性扩增。引物由生工生物工程（上海）股份有限公司合成，合成后的引物经质检合格后，用无菌水配制成10μM的工作浓度，保存于-20℃冰箱备用。以反转录得到的cDNA为模板进行qRT-PCR扩增。采用SYBRGreen染料法，反应体系总体积为20μl，其中包含2×SYBRGreenMasterMix10μl，这是反应的核心试剂，包含了DNA聚合酶、dNTPs、SYBRGreen荧光染料等，能够在PCR扩增过程中，特异性地与双链DNA结合，在荧光信号的激发下发出荧光，从而实现对扩增产物的实时监测；上下游引物（10μM）各0.8μl，引物的浓度经过优化，既能保证引物与模板充分结合，又能避免引物浓度过高导致的非特异性扩增；cDNA模板1μl，模板的用量根据前期预实验确定，保证在合适的线性范围内进行扩增；ddH₂O7.4μl，用于补足反应体系的体积。反应在ABI7500实时荧光定量PCR仪上进行，反应程序如下：95℃预变性30s，这一步的目的是使模板DNA完全变性，解开双链，为后续引物与模板的结合提供单链模板；然后进行40个循环，每个循环包括95℃变性5s，使双链DNA解链为单链，便于引物结合；60℃退火30s，在此温度下，引物与模板特异性结合；72℃延伸30s，DNA聚合酶在引物的引导下，以dNTPs为原料，沿着模板链进行延伸，合成新的DNA链。在每个循环的延伸阶段，收集荧光信号，通过检测荧光强度的变化，实时监测PCR扩增产物的积累情况。扩增结束后，进行熔解曲线分析，从60℃缓慢升温至95℃，每升高0.5℃采集一次荧光信号，绘制熔解曲线。熔解曲线用于判断扩增产物的特异性，若熔解曲线只有单一的峰，说明扩增产物为特异性产物；若出现多个峰，则可能存在非特异性扩增或引物二聚体等问题。以茄子的Actin基因作为内参基因，用于校正和标准化目标基因的表达量。内参基因在不同组织和不同实验条件下表达相对稳定，其表达水平不受实验处理因素的影响，因此可以作为参照，消除不同样本在RNA提取、反转录和PCR扩增等过程中可能存在的差异，使不同样本之间目标基因的表达量具有可比性。采用2⁻ΔΔCt法计算目标基因的相对表达量。首先，计算每个样本中目标基因的Ct值（CycleThreshold，循环阈值，是指在PCR扩增过程中，荧光信号达到设定的阈值时所经历的循环数）与内参基因Ct值的差值，即ΔCt=Ct目标基因-Ct内参基因；然后，以其中一个样本作为对照，计算其他样本与对照样本的ΔCt差值，即ΔΔCt=ΔCt样本-ΔCt对照；最后，根据公式2⁻ΔΔCt计算目标基因的相对表达量，该值表示目标基因在不同样本中的相对表达倍数。将qRT-PCR得到的基因相对表达量与转录组测序得到的FPKM值进行相关性分析。结果显示，两者之间呈现出显著的正相关关系，相关系数r达到了0.85（P<0.01）。这表明通过转录组测序筛选出的差异表达基因在qRT-PCR验证中表现出了一致的表达趋势，从而验证了转录组测序结果的可靠性。例如，在转录组测序中，基因A在果皮着色深的样本中FPKM值显著高于果皮着色浅的样本，在qRT-PCR验证中，该基因在果皮着色深的样本中的相对表达量同样显著高于果皮着色浅的样本，且表达倍数变化趋势与转录组测序结果相近。通过qRT-PCR验证，为后续深入研究这些差异表达基因在紫茄果皮着色深度调控中的作用提供了可靠的数据支持，确保了研究结果的准确性和可信度。四、紫茄果皮着色深度与相关基因的关联分析4.1基因表达与着色深度的相关性分析运用统计学方法，深入分析差异表达基因的表达量与紫茄果皮着色深度之间的相关性。将转录组测序获得的基因表达数据与前期通过色差仪精确测定的果皮着色深度数据进行整合，利用Pearson相关性分析等方法，计算每个差异表达基因表达量与果皮着色深度之间的相关系数。通过分析发现，有多个基因的表达量与紫茄果皮着色深度呈现出显著的相关性。其中，基因A的表达量与果皮着色深度呈极显著正相关，相关系数r达到了0.85（P<0.01）。这表明随着基因A表达量的增加，紫茄果皮着色深度也显著增加，基因A在促进紫茄果皮着色方面发挥着重要作用。进一步研究发现，基因A编码的蛋白质参与了花青素合成途径中的关键步骤，可能通过调控花青素合成酶的活性，促进花青素的合成，从而加深紫茄果皮的着色。基因B的表达量与果皮着色深度呈显著负相关，相关系数r为-0.78（P<0.05）。这意味着基因B的高表达会抑制紫茄果皮着色深度，可能是通过抑制花青素合成途径中的某些关键基因表达，或者参与花青素的降解过程，从而使果皮着色变浅。通过对基因B功能的初步分析，发现其编码的蛋白质可能与某些转录因子相互作用，调控花青素合成相关基因的表达。为了更直观地展示基因表达量与果皮着色深度之间的关系，绘制了散点图。在散点图中，以果皮着色深度为横坐标，基因表达量为纵坐标，每个点代表一个样本。从图中可以清晰地看出，对于与果皮着色深度呈正相关的基因，其表达量随着果皮着色深度的增加而呈现上升趋势，点的分布大致沿着一条从左下角到右上角的直线；而对于呈负相关的基因，其表达量随着果皮着色深度的增加而下降，点的分布沿着从左上角到右下角的方向。通过拟合曲线可以更准确地描述这种关系，进一步验证了相关性分析的结果。除了上述显著相关的基因外，还发现了一些基因与果皮着色深度之间存在较弱的相关性。虽然这些基因的相关系数绝对值相对较小，但它们可能在果皮着色深度的调控中发挥着协同作用或间接作用。例如，某些基因可能参与了与花青素合成相关的代谢途径，通过影响其他代谢产物的合成或浓度，间接影响花青素的合成和积累，从而对果皮着色深度产生一定的影响。对这些基因的深入研究，有助于全面揭示紫茄果皮着色深度的调控网络。4.2关键基因的功能验证为了深入探究紫茄果皮着色深度的分子机制，对筛选出的关键基因进行功能验证是至关重要的环节。采用病毒诱导的基因沉默（VIGS）技术对基因A进行功能验证。该技术利用病毒载体将目的基因的部分片段导入植物体内，引发植物自身的RNA干扰机制，从而特异性地沉默目标基因的表达。根据基因A的序列信息，设计并合成特异性的干扰片段。将干扰片段克隆到烟草脆裂病毒（TRV）载体上，构建重组表达载体TRV-A。利用农杆菌介导的转化方法，将重组表达载体TRV-A导入到含有TRV1的农杆菌GV3101菌株中，同时将空载体TRV导入含有TRV1的农杆菌GV3101菌株作为对照。将含有重组表达载体和空载体的农杆菌分别培养至对数生长期，收集菌体，用含有10mMMgCl₂、10mMMES和200μM乙酰丁香酮的重悬液重悬菌体，调整菌液浓度至OD₆₀₀=1.0，室温静置3小时，使菌体充分活化。将活化后的含有TRV-A和TRV的农杆菌菌液等体积混合，采用注射器浸润法，将混合菌液注射到紫茄幼苗的子叶中。在注射过程中，要确保菌液均匀地渗透到子叶组织中，避免损伤子叶。将注射后的紫茄幼苗置于光照培养箱中，在25℃、光照16小时/黑暗8小时的条件下培养，定期观察植株的生长情况。在果实发育至商品成熟期，利用色差仪测定沉默基因A后紫茄果实的果皮着色深度，并与对照植株进行对比。结果显示，沉默基因A后，紫茄果皮着色深度显著变浅，与对照植株相比，其L值（明度值）显著升高，表明果皮颜色变得更浅；a值（红度值）显著降低，说明果皮的红色程度减弱，紫色也相应变浅。这表明基因A在促进紫茄果皮着色深度方面发挥着关键作用，其表达的沉默导致花青素合成受阻，进而使果皮着色变浅。通过实时荧光定量PCR检测基因A及花青素合成途径中相关基因的表达量，发现基因A的表达量在沉默植株中显著降低，同时，下游的花青素合成关键酶基因，如查尔酮合酶（CHS）基因、二氢黄酮醇4-还原酶（DFR）基因等的表达量也显著下调。这进一步证实了基因A通过调控花青素合成途径中相关基因的表达，来影响紫茄果皮的着色深度。对于基因B，采用基因过表达技术进行功能验证。构建基因B的过表达载体，利用PCR技术从紫茄基因组中扩增出基因B的全长编码序列，将其克隆到植物表达载体pCAMBIA1302上，该载体含有CaMV35S启动子，能够驱动基因在植物体内高效表达。将构建好的过表达载体pCAMBIA1302-B通过冻融法转化到农杆菌GV3101中。采用叶盘转化法，将含有过表达载体的农杆菌侵染紫茄叶片外植体。在无菌条件下，将紫茄叶片切成0.5cm×0.5cm的小块，放入含有农杆菌菌液的侵染液中浸泡10分钟，期间轻轻摇晃，使叶片充分接触农杆菌。将侵染后的叶片外植体放在含有筛选培养基（含有卡那霉素和头孢霉素）的培养皿中，在25℃、光照16小时/黑暗8小时的条件下培养，诱导愈伤组织和不定芽的分化。待不定芽长至2-3cm时，将其切下，转移到含有生根培养基的培养瓶中，诱导生根。经过一段时间的培养，获得转基因过表达植株。对转基因过表达植株的果实进行分析，发现果皮着色深度明显变浅。与野生型植株相比，转基因植株果实的L值显著升高，a值显著降低，表明基因B的过表达抑制了紫茄果皮的着色。进一步分析花青素含量和相关基因表达，发现转基因植株果实中的花青素含量显著低于野生型植株，同时花青素合成途径中多个关键基因的表达受到抑制，如黄酮醇合成酶（FLS）基因、花青素合成酶（ANS）基因等。这说明基因B可能通过抑制花青素合成途径中关键基因的表达，减少花青素的合成，从而导致紫茄果皮着色变浅。通过对基因A和基因B的功能验证，明确了它们在紫茄果皮着色深度调控中的作用机制，为深入理解紫茄果皮着色的分子机制提供了重要的实验依据，也为紫茄的遗传改良和品质提升奠定了坚实的基础。未来还可以进一步研究这些关键基因与其他基因之间的相互作用关系，以及它们在不同环境条件下的表达调控模式，以全面揭示紫茄果皮着色深度的调控网络。4.3基因调控网络的构建基于基因之间的相互作用关系和功能注释，运用生物信息学工具和相关数据库，构建紫茄果皮着色深度相关基因的调控网络，以全面揭示基因调控的复杂性。借助STRING（SearchToolfortheRetrievalofInteractingGenes/Proteins）数据库，该数据库整合了大量来自不同物种的蛋白质-蛋白质相互作用信息，包括实验验证的相互作用以及通过计算预测的相互作用。将筛选出的与紫茄果皮着色深度相关的差异表达基因输入到STRING数据库中，设置物种为茄子（Solanummelongena），数据库会根据已有的信息，预测这些基因所编码蛋白质之间的相互作用关系。例如，基因A编码的蛋白质与基因B、基因C编码的蛋白质存在直接的相互作用，这种相互作用可能是物理上的结合，也可能是通过信号传导途径间接影响彼此的功能。通过STRING数据库的分析，能够初步构建出基因之间的相互作用网络框架。利用Cytoscape软件对STRING数据库的分析结果进行可视化展示。Cytoscape是一款功能强大的生物网络分析和可视化软件，能够将基因之间的相互作用关系以直观的图形形式呈现出来。在Cytoscape软件中，每个差异表达基因被表示为一个节点，基因之间的相互作用关系则用边来连接。节点的大小、颜色等属性可以根据基因的某些特征进行设置，例如基因的表达量、在花青素合成途径中的重要性等；边的粗细、颜色等也可以用来表示相互作用的强度或类型。通过这种可视化的方式，可以更清晰地观察基因之间的相互关系，发现基因调控网络中的关键节点和重要连接。结合基因的功能注释信息，进一步完善基因调控网络。在GO和KEGG等数据库中，对每个差异表达基因的功能进行详细注释，了解它们在花青素合成途径、信号传导通路以及其他与果皮着色相关的生物学过程中的具体作用。例如，某些基因被注释为参与花青素合成途径中的关键酶的编码，这些基因在调控网络中可能处于核心位置，直接影响花青素的合成和积累；而一些转录因子基因，虽然不直接参与花青素的合成，但通过与其他基因的启动子区域结合，调控基因的表达，在调控网络中起着调节开关的作用。将这些功能信息整合到基因调控网络中，能够更深入地理解基因之间的调控机制。在构建的基因调控网络中，发现了多个基因模块。这些基因模块由一组相互作用密切的基因组成，它们可能共同参与某个特定的生物学过程。例如，一个基因模块中包含了多个花青素合成途径中的结构基因，如查尔酮合酶（CHS）基因、二氢黄酮醇4-还原酶（DFR）基因、花青素合成酶（ANS）基因等，这些基因在模块中相互作用，协同调控花青素的合成。另一个基因模块中则主要包含转录因子基因，如MYB类转录因子基因、bHLH类转录因子基因等，它们通过与结构基因的启动子区域结合，调控结构基因的表达，从而间接影响花青素的合成和果皮着色深度。通过对基因调控网络的分析，明确了基因之间的上下游关系。例如，转录因子基因A可能通过与结构基因B的启动子区域结合，激活结构基因B的表达，从而促进花青素的合成；而结构基因C的表达产物可能反馈调节转录因子基因D的表达，形成一个复杂的反馈调控回路。这种上下游关系的明确，有助于深入理解紫茄果皮着色深度的调控机制，为进一步研究基因功能和调控网络的优化提供了重要线索。基因调控网络的构建为深入研究紫茄果皮着色深度的分子机制提供了一个全面的框架，有助于揭示基因之间的复杂相互作用关系，发现潜在的调控靶点，为紫茄的遗传改良和品质提升提供更深入的理论依据。未来还可以结合更多的实验数据，如蛋白质-蛋白质相互作用的验证实验、基因表达的时空动态变化数据等，进一步完善和优化基因调控网络，使其更加准确地反映紫茄果皮着色深度的调控过程。五、结论与展望5.1研究结论总结本研究通过对紫茄果皮着色深度的遗传及相关基因差异表达进行深入探究，取得了一系列重要成果。在遗传分析方面，选用果皮着色深度差异显著的紫茄高代自交系EP26和EP28构建6个世代群体，运用主基因+多基因混合模型遗传分析方法，明确了紫茄果皮着色深度符合两对加性-显性-上位性主基因+加性-显性多基因遗传的E-1模型。主基因加性效应值均等，遗传以显性效应为主。在不同世代群体中，F2、B1和B2主基因遗传率分别为70.25%、3.46%、47.41%，多基