番茄SMIFH2处理下根系表型鉴定及差异基因表达分析

目录TOC\o"2-3"\h\z\u\t"毕设标题1,1"1前言 22材料与方法 42.1植物材料、生长条件与处理 42.1.1植物材料与种子处理 42.1.2幼苗培养与生长条件 42.2根系表型分析 42.2.1主根长度测量 52.2.2侧根数量与发生时间 52.2.3根毛密度观察与量化 52.3RNA分离与测序 52.3.1样品收集 52.3.2RNA提取、文库构建与测序 52.4转录组数据分析 52.5统计学分析 53.结果 63.1SMIFH2处理对番茄主根生长的影响 63.2SMIFH2处理对番茄侧根的发育的影响 73.3SMIFH2处理对番茄根毛密度的影响 93.440µMSMIFH2处理96h番茄根尖RNA-Seq检测结果分析 93.5差异表达基因的功能富集分析 104.讨论 124.1SMIFH2对番茄根系生长的抑制效应及其与肌动蛋白功能紊乱的关联 124.2高浓度SMIFH2对首个侧根发生的潜在双重效应 124.3转录组分析提示细胞壁代谢和胁迫响应通路参与SMIFH2效应 134.4研究局限性与展望 13参考文献 14番茄SMIFH2处理下根系表型鉴定及差异基因表达分析摘要：成蛋白通过调控肌动蛋白微丝形成与动态重组参与植物发育，但其在番茄根系发育中的作用尚未明确。本研究利用成蛋白抑制剂SMIFH2处理番茄幼苗，分析其对根系表型的影响及潜在分子机制。实验设置0µM（对照）、20µM和40µMSMIFH2处理组，在不同时间点测量主根长度、侧根数量、首个侧根发生时间以及主根特定区域（0-5mm,5-10mm）的根毛密度。对处理96小时后的对照组与40µMSMIFH2处理组的根尖样品进行转录组分析以筛选差异表达基因，随后进行GO和KEGG富集分析。本研究证实SMIFH2对番茄根系生长具有显著的抑制作用，但高浓度下可能加速早期侧根的发生。与对照组相比，SMIFH2处理组共鉴定出64个差异表达基因，显著富集于细胞壁代谢、氧化还原反应和胁迫响应等GO条目，以及戊糖与葡萄糖醛酸转化、苯丙烷类和类黄酮生物合成等KEGG通路。研究结果表明，SMIFH2可能通过干扰细胞壁动态与诱导胁迫相关代谢致使番茄根系结构与功能的异常发育。关键词：番茄；SMIFH2；成蛋白；根系发育；转录组；差异表达基因RootphenotypiccharacterizationanddifferentialgeneexpressionanalysisintomatounderSMIFH2treatmentAbstract:Forminsregulatetheformationanddynamicreorganizationofactinfilamentsandareessentialforplantdevelopment,yettheirspecificrolesintomatorootdevelopmentremainunclear.Inthisstudy,weappliedtheformininhibitorSMIFH2totomatoseedlingstoinvestigateitseffectsonrootmorphologyandunderlyingmolecularmechanisms.Seedlingsweretreatedwith0µM(control),20µM,and40µMSMIFH2,andphenotypicparametersincludingprimaryrootlength,lateralrootnumber,thetimingofthefirstlateralrootemergence,androothairdensityinthe0–5mmand5–10mmregionsoftheprimaryrootweremeasuredatdifferenttimepoints.Transcriptomicanalysiswasconductedonroottipscollected96hoursaftertreatmentfromboththecontroland40µMSMIFH2-treatedgroupstoidentifydifferentiallyexpressedgenes(DEGs),followedbyGOandKEGGenrichmentanalyses.TheresultsdemonstratedthatSMIFH2significantlyinhibitedtomatorootgrowth,whilehighconcentrationsappearedtopromoteearlierinitiationoflateralroots.Atotalof64DEGswereidentifiedintheSMIFH2-treatedgroupcomparedtothecontrol,withsignificantenrichmentinGOtermsrelatedtocellwallmetabolism,redoxprocesses,andstressresponses,andinKEGGpathwaysincludingpentoseandglucuronateinterconversions,phenylpropanoidbiosynthesis,andflavonoidbiosynthesis.ThesefindingssuggestthatSMIFH2maydisruptcellwalldynamicsandactivatestress-relatedmetabolicpathways,leadingtostructuralandfunctionalabnormalitiesintomatorootdevelopment.Keywords:Solanumlycopersicum;SMIFH2;Formin;Rootdevelopment;Transcriptome;Differentiallyexpressedgenes1前言番茄（Solanumlycopersicum），作为茄科植物的重要代表，不仅是全球范围内广泛消费的蔬菜作物，其多样的遗传资源和相对清晰的遗传背景也使其成为研究植物生长发育、果实成熟及胁迫响应机制的重要模式植物。其中AilsaCraig(AC)因其广泛应用和良好的实验特性，常被用于各类生理及分子生物学研究。植物的根系系统是其生命活动的基础，承担着固定植株、吸收土壤水分和矿质营养、合成部分激素以及与土壤微生物互作等多重关键功能。根系的构型，包括主根的生长、侧根的发生与分布、以及根毛的发育状况，共同影响着植物从土壤中获取资源的能力，进而深刻影响地上部分的生长表现、最终产量及对环境变化的适应能力。例如，对番茄根系在盐胁迫下发育可塑性的表型分析揭示了侧根位置等参数对耐盐性的重要性[1]。因此，深入解析番茄根系发育的调控网络，对于指导番茄遗传改良和优化栽培管理措施具有重要的理论与实践意义。植物细胞的形态建成与功能发挥高度依赖于动态变化的细胞骨架系统，该系统主要由微管、肌动蛋白微丝和中间纤维构成。其中，肌动蛋白细胞骨架通过其动态组装和重构调控着细胞分裂、形态建立、极性生长、细胞器运输等关键过程[2]。特别是在根系中，肌动蛋白的作用尤为重要，因为根细胞发育涉及不同阶段的精确控制，如从细胞分裂到伸长、分化、毛状细胞形成等。细胞骨架药物的超敏反应证明了肌动蛋白与微管网络的动态互作（cross-talk）在调控细胞拟南芥根细胞的各向异性伸长中起着重要作用[3]。肌动蛋白骨架的重塑对于番茄细胞响应生物胁迫也显示出不同的模式[4]。成蛋白（Formin）是连接细胞骨架与内膜系统的重要分子枢纽，其功能具有广泛性与复杂性[5]。它能够从头开始成核肌动蛋白，并通过成束蛋白同源域（ForminHomology,FH）1和2启动肌动蛋白单体的聚合，形成无分支的线性肌动蛋白丝[6][7]。植物基因组通常编码大量的成蛋白，暗示了它们在植物生命周期中参与调控多样化的细胞过程。已有的研究表明，植物成蛋白参与了细胞极性生长、细胞分裂过程中细胞板的形成与稳定、胞间连丝的功能调节以及对生物和非生物胁迫的响应等多个方面。例如，在模式植物拟南芥中的研究揭示了不同成蛋白在根毛发育、细胞分裂调控等方面的特定功能。SMIFH2（SmallMoleculeInhibitorofForminHomology2domain）作为一种被专门设计用于抑制成蛋白家族的活性的小分子化合物，其通过特异性结合成蛋白的FH2结构域，有效阻断其介导的肌动蛋白成核与聚合过程，并减少成蛋白与微管的结合[8]。类似地，拟南芥中的AtFH8formin也被证实在低剂量肌动蛋白抑制剂LatrunculinB存在下参与根系发育和细胞分裂过程[9]。在拟南芥表皮细胞中，SMIFH2处理证明了Profilin依赖性的肌动蛋白成核和组装对于细胞伸长的重要性[10]。SMIFH2也被用于解析成蛋白在烟草和拟南芥细胞分裂过程中调控肌动蛋白-微管互作的功能[11]。最后，SMIFH2处理能够显著干扰拟南芥肌动蛋白细胞骨架的动态与组织，并影响其微管结构，进而导致一系列发育缺陷，包括根系生长受到抑制、细胞扩张模式改变以及细胞分裂过程异常[12]。这些在拟南芥中获得的分子和表型层面的证据，凸显了利用SMIFH2作为研究成蛋白功能的有效药理学工具的潜力。尽管SMIFH2在拟南芥等模式植物中的应用揭示了成蛋白在细胞骨架调控和发育中的广泛作用，但其在番茄根系中的具体效应及下游分子机制仍知之甚少。目前关于番茄根系的研究多集中于其他方面，如机械刺激、盐胁迫、生长素信号传导、病原菌侵染等。且现有关于番茄根系的研究虽然积累了一定的转录组数据，但这些研究主要聚焦于激素响应或环境胁迫等其他生物学背景，尚无研究将SMIFH2这类细胞骨架特异性抑制剂的处理效应，与其诱导的番茄根系表型变化，以及全基因组范围内的基因表达响应直接关联起来进行整合分析。因此，本研究旨在探索SMIFH2扰动下番茄根系表型与转录组谱之间联系，通过施用高低浓度（40µM和20µM）的SMIFH2，系统地鉴定并量化其对番茄幼苗主根伸长、侧根数量、首个侧根出现时间以及根毛密度的影响。在此表型分析的基础上，结合转录组测序技术（RNA-Seq），全面解析SMIFH2处理下番茄根系中差异表达的基因谱。通过对这些差异基因进行功能注释和通路富集分析，期望实现表型与转录组数据的整合，揭示SMIFH2影响番茄根系发育的关键下游基因和可能涉及的生物学通路，为深入理解成蛋白和肌动蛋白细胞骨架在调控番茄根系结构建成中的复杂作用网络提供新的实验证据和分子层面的见解。2材料与方法2.1植物材料、生长条件与处理2.1.1植物材料与种子处理本研究选用AC品种的番茄种子。为确保实验材料的均一性，挑选外观无损伤、大小均匀且颗粒饱满的种子备用。为了避免种子表面可能携带的微生物污染影响后续无菌培养，在超净工作台中进行严格的表面消毒程序：首先，将种子完全浸没于75%（v/v）分析纯乙醇溶液中，并伴随持续轻柔摇动1分钟，以确保乙醇与种子表面充分接触；迅速弃去乙醇后，立即用大量无菌去离子水（ddH₂O）清洗三次，每次冲洗持续2分钟并伴随摇动，以彻底去除残留乙醇。接着，将种子转移至新配的1:2体积比的次氯酸钠溶液中，再次浸泡并持续摇动10分钟进行深度消毒。随后，为彻底去除次氯酸钠残留，避免其对种子萌发和幼苗生长的抑制作用，用足量无菌去离子水重复清洗三次，每次清洗时间延长至5分钟。完成消毒后的种子置于含有约10mL无菌去离子水的组培瓶中。为打破种子休眠并促进萌发同步性，将其置于4℃恒温冰箱中进行低温层积处理，处理时间精确控制在48小时。2.1.2幼苗培养与生长条件层积处理结束后，在无菌操作条件下，将种子均匀播种在铺有两层无菌定性滤纸培养皿中。向培养皿中加入约3mL无菌去离子水，以保持滤纸湿润，为种子萌发提供必要水分。将培养皿封口后垂直置于光照培养箱中，设置温度为25°C，黑暗条件下进行萌发，持续2-3天。每日观察种子萌发情况，待幼苗的初生根长度发育至约1-1.5cm，且下胚轴初步伸长时，选取整体长势、根长、下胚轴长度等方面表现一致且健康的幼苗用于后续的接种实验。本实验采用改良的固体培养基对幼苗进行培养。基础培养基成分为1/2浓度的MurashigeandSkoog（MS）基础盐类混合物（2.15g/L），额外添加蔗糖（15g/L）作为碳源，并使用琼脂（8g/L）作为凝固剂。在定容前，使用数字pH计并滴加1MNaOH或1MHCl溶液将培养基pH值精确调节至5.8±0.02。将配制好的培养基装于耐高压玻璃瓶中，于121°C、103kPa条件下进行高压蒸汽灭菌0.5h，以确保培养基彻底无菌。在超净工作台中，将完成灭菌并冷却至约55°C的培养基溶液倒入12cm×12cm见方的无菌方形培养皿中，每皿分装约40-50mL，水平放置待其凝固后，即为对照组培养基。对于处理组，将预先使用DMSO（二甲基亚砜）溶解并通过0.22µm无菌滤膜过滤除菌的SMIFH2储存液，用移液枪精确加入已灭菌并冷却至同样温度的1/2MS基础培养基中，轻轻混匀以避免产生气泡，最终使培养基的SMIFH2浓度分别为20µM和40µM。然后将含有不同浓度SMIFH2的培养基分别倒入10cm×10cm的无菌方形培养皿中，每皿体积约为40mL。待培养基完全凝固后，小心地将前述筛选好的生长状态一致的番茄幼苗，根尖朝下，转移到培养基表面，确保根部与培养基良好接触。每个培养皿均匀放置10-13株幼苗。随后用透气封口膜密封培养皿边缘，并将培养皿垂直放置于可控环境的组培室中。设定培养条件为：恒定温度20°C±1°C，控制光周期为每日16小时光照后8小时黑暗（16h/8hL/D）。光照由顶部排列的冷白荧光灯管提供，光照强度约60μmol·m⁻²·s⁻¹。为保证实验结果的可靠性，每个处理设置3个独立的生物学重复，每个生物学重复包含10-13株重复的幼苗。2.2根系表型分析2.2.1主根长度测量在幼苗接种后的第48小时、72小时和96小时这三个关键时间点，将培养皿从培养架上取出，背面朝上水平放置在带有光源的拍摄平台上。在培养皿旁放置钢制刻度尺作为参照，拍摄包含所有幼苗根系的清晰数字图像。随后，使用图像分析软件ImageJ（v1.8.0）进行主根长度的量化分析。首先，利用图像中刻度尺的已知长度设定软件的比例尺。然后，使用软件中的折线选择工具，测量每株幼苗从主根基部到根尖最末端的长度。记录每个时间点、每个处理、每个重复中所有单株幼苗的主根长度数据。2.2.2侧根数量与发生时间在接种后48h、72h和96h时，记录每株幼苗主根上已萌发突破表皮的侧根数量。自接种后，每隔8h观察并记录每株幼苗第一条侧根出现的时间。必要时，使用光学显微镜观察以确定侧根是否产生。在接种后48小时、72小时和96小时这三个时间点，在测量主根长度的同时，仔细观察并记录每株幼苗主根上已经萌发并突破主根表皮的可见侧根数量。侧根的计数标准为：已形成清晰可见的、突出主根表面的侧根原基或已伸长的侧根。为追踪侧根的起始发生时间，自接种开始后，每隔8小时通过培养皿非侵入式地观察每株幼苗的主根，并记录首次观察到侧根突破表皮的时间点。若在特定时间点肉眼难以明确判断侧根是否已产生，则使用光学显微镜进行确认观察，以确保侧根发生判定的准确性。2.2.3根毛密度观察与量化在接种后96h时，使用光学显微镜对每株幼苗的主根进行观察并拍照（含比例尺）。计数照片中距离根尖0-5mm和5-10mm两个区域内的根毛数量，计算根毛密度（根毛数/mm）。在接种培养至96h时，选取每个处理组每个重复中具有代表性的幼苗10-20株进行根毛密度的细致观察与量化。使用配备有数码摄像头的光学显微镜对每株选定幼苗的主根特定区域进行观察并拍照（含比例尺）。重点关注主根根尖后方的两个连续区域：距离根尖生长点0-5mm（根毛起始区和伸长区）和5-10mm（根毛成熟区）。在获取的显微图像上，计数这两个指定区域内可见的根毛数量，随后计算每个区域的根毛密度，表示为单位长度主根上的根毛数量（根毛数/mm）。2.3RNA分离与测序2.3.1样品收集在接种后96h时，分别收集对照组和40µMSMIFH2处理组的番茄幼苗根尖样品0.01g。每个处理组设置3个独立的生物学重复样本用于RNA测序。样品收集后迅速置于已预冷的无菌冻存管中，并立即投入液氮中进行速冻，以最大限度抑制RNA酶活性，保证RNA完整性。随后将冻存管转移至-80℃超低温冰箱中长期保存，直至进行RNA提取。2.3.2RNA提取、文库构建、测序与测序数据的预处理冷冻的根尖样品交由商业测序服务提供商完成总RNA的提取、使用DNase处理去除基因组DNA、进行RNA质量评估、基于mRNA富集的链特异性RNA-Seq文库构建以及最终的Illumina平台高通量测序。测序数据的处理，包括测序原始图像数据转化为原始测序序列、对原始测序数据进行质量控制和预处理、把清洁读段比对至番茄参考基因组、统计每个基因的读段计数等步骤同样由测序服务商完成。2.4转录组数据分析所有转录组生物信息学分析均使用R语言（v4.4.2）完成，除非另有说明。利用BioConductor中的DESeq2包对基因表达计数矩阵进行基因差异表达分析。该分析旨在比较40µMSMIFH2处理组与对照组之间基因表达水平的统计学差异。筛选差异表达基因（DifferentiallyExpressedGenes,DEGs）的标准设定为：经过Benjamini-Hochberg方法校正后的p值（adjustedp-value,padj）<0.05并且基因表达量变化的绝对值对数（以2为底）（|Log2(FoldChange)|）≥1。使用ggplot2包，以|Log2(FoldChange)|为横坐标，-Log10(padj)为纵坐标，绘制火山图（VolcanoPlot）展示DEGs的整体分布情况，并用不同颜色或标记区分上调、下调及无显著差异的基因。为了探索筛选出的DEGs潜在的生物学功能和参与的代谢或信号通路，使用ClusterProfiler包进行GeneOntology（GO）功能富集分析和KEGG（KyotoEncyclopediaofGenesandGenomes）通路富集分析。GO分析依据GeneOntologyConsortium提供的数据库，将DEGs注释到生物过程（BiologicalProcess,BP）、细胞组分（CellularComponent,CC）和分子功能（MolecularFunction,MF）三个主要类别中，并识别哪些GO条目在DEGs列表中被显著富集。KEGG分析则将DEGs映射到KEGG数据库收录的已知代谢通路和信号转导通路图谱上，展示DEGs显著参与的生物学通路。这两种富集分析的显著性阈值设定为p<0.05。富集分析的结果通过条形图（Barplot）和气泡图（Bubbleplot）形式进行数据可视化。2.5统计学分析所有涉及根系表型测量的实验数据均采用IBMSPSSStatistics19软件进行统计分析。在满足方差齐性的前提下，采用单因素方差分析（One-wayANOVA）检验不同SMIFH2浓度处理组之间在各个测量指标上是否存在显著差异，采用Tukey’s-b法进行多重比较，以确定差异达到统计学显著性水平的具体处理组，并在统计图中以不同字母标注表示处理间的差异。所有统计显著性水平设定为p<0.05。实验数据在图表中表示为平均值±标准误差（Mean±SE）。RNA-Seq数据的统计分析在R中完成，如2.4节所述。3结果3.1SMIFH2处理对番茄主根生长的影响如图1所示，对照组（Control）幼苗的主根长度随时间稳步增长，在处理后48h、72h和96h的平均长度分别达到3.74cm、5.33cm和6.57cm。与对照组相比，两种浓度的SMIFH2处理均对主根的伸长产生了显著的抑制作用。在处理后48小时，20µM和40µMSMIFH2处理组的主根长度均值为2.89cm和3.24cm，已显著低于对照组（图1）。这种抑制效应在72小时和96小时持续存在。至处理后96小时，20µM和40µMSMIFH2处理组的主根长度均值为5.03cm和5.56cm，均显著短于对照组且长度差值更大（图1），表明SMIFH2处理能够有效抑制番茄主根的生长。20µMSMIFH2处理下的主根比40µM下的略短，然而长度差异始终不显著。图2直观地展示了不同处理不同时间的根系表型之间的主根长度差别。图SEQ图\*ARABIC1不同浓度SMIFH2处理主根长度图SEQ图\*ARABIC2不同处理不同时间根系表型对照组不同时间点根系表型（图2-a,图2-b,图2-c）；20µMSMIFH2处理下不同时间点根系表型（图2-d,图2-e,图2-f）；40µMSMIFH2处理不同时间点根系表型（图2-g,图2-h,图2-i）3.2SMIFH2处理对番茄侧根的发育的影响侧根数量和首次出现的时间构成SMIFH2对番茄侧根发育影响的分析。如图3所示，在接种后48h时，各处理组间的侧根数量无显著差异。在接种后72h时，20µMSMIFH2处理组的侧根数显著低于对照组，而40µM处理组与对照组相比无显著差异，但表现出减少的趋势。至接种后96h时，20µM和40µMSMIFH2处理组的侧根数量均显著少于对照组。这说明SMIFH2处理，尤其是在较长时间点，抑制了番茄侧根的生长。值得注意的是，实验结果显示高浓度的SMIFH2处理对首个侧根出现的时间具有促进效应（图4）。对照组和20µMSMIFH2处理组的首个侧根平均分别出现在处理后52.78小时和49.32小时，两者之间无显著差异。然而，40µMSMIFH2处理显著将首个侧根的出现的平均时间提前至31.23小时。这表明尽管高浓度SMIFH2最终减少了侧根总数，但促进了侧根的发生。图2直观地反映了这个结果，在48h时，多个经SMIFH2处理的幼苗已有侧根发生（图2-d，图2-g），而对照组中没有幼苗产生侧根（图2-a）。图SEQ图\*ARABIC3不同浓度SMIFH2处理侧根数图SEQ图\*ARABIC4不同浓度SMIFH2处理首个侧根发生时间3.3SMIFH2处理对番茄根毛密度的影响根毛是植物吸收水分和养分的重要结构。实验中选取定植96小时后不同处理组在主根近根尖区域（0-5mm）和稍成熟区域（5-10mm）的根毛密度进行比较。如图5所示，SMIFH2处理对根毛密度产生了抑制作用，SMIFH2处理组的根毛更为稀疏，且长度更短。在0-5mm区域，40µMSMIFH2处理组的根毛密度显著低于对照组和20µM处理组。在5-10mm区域，40µMSMIFH2处理组的根毛密度同样显著低于对照组，而20µM处理组的根毛密度介于对照组和40µM组之间，与两者均无显著差异（图6）。这些结果表明，SMIFH2处理，特别是40µM的高浓度，显著降低了番茄主根的根毛密度。图SEQ图\*ARABIC5光学显微镜下不同浓度SMIFH2处理的根毛密度对照组96h根毛密度的观察（图5-a）；20µMSMIFH296h根毛密度的观察（图5-b）；40µMSMIFH2根毛密度的观察（图5-c）图SEQ图\*ARABIC6不同浓度SMIFH2处理平均根毛密度3.440µMSMIFH2处理96h番茄根尖RNA-Seq检测结果分析对RNA-seq检测结果的分析体现出SMIFH2处理与对照组在转录水平上的差异。通过设定padj<0.05且|log2(FoldChange)|>1.0作为阈值，共筛选到64个显著差异表达基因DEGs。DEGs火山图（图8）对所有检测到的基因的表达变化倍数和统计显著性进行了可视化。结果显示，与对照组相比，SMIFH2处理后有41个基因显著上调，23个基因显著下调。此外，图8中标注了一些差异表达最为显著的基因编号（GeneID），如gene-LOC101267672（下调）和gene-LOC101249934（上调）。通过这些GeneID获得基因的功能描述，进一步分析分子层面的细节变化。层次聚类热图（图7）直观地展示了DEGs在不同样品间的表达模式。其中，对照组（WT）的三个生物学重复聚为一簇，SMIFH2处理组（S40）的三个生物学重复聚为另一簇，表明两组样品间存在稳定且显著不同的基因表达谱。对基因的聚类分析(行聚类)也显示出不同的表达模式簇，存在一组基因在SMIFH2处理组样品中普遍高表达而在对照组样品中低表达，以及另一组基因在SMIFH2处理组样品中普遍低表达而在对照组样品中高表达。这表明SMIFH2处理与对照组相比存在系统性的基因表达调控变化。3.5差异表达基因的功能富集分析GeneOntology（GO）富集分析体现了这些DEGs潜在的生物学功能（图9）。在生物学过程（BP）类别中，DEGs显著富集于与碳水化合物代谢和细胞壁修饰相关的过程，如碳水化合物运输、果胶分解、半乳糖醛酸代谢和多糖分解。此外，与胁迫响应相关的对氧化应激的反应和防御反应等条目也被显著富集。在细胞组分（CC）方面，DEGs主要富集在膜、膜内在组分、质膜和细胞外区域。分子功能（MF）富集分析则指向了氧化还原酶活性、裂解酶活性、跨膜转运蛋白活性以及与细胞壁降解相关的酶活性，如果胶酯酶活性和多聚半乳糖醛酸酶活性。进一步的KEGG通路富集分析揭示了DEGs主要参与的代谢通路（图10）。显著富集的通路包括核心的代谢途径和次生代谢物的生物合成。与植物特有代谢产物合成相关的通路，如苯丙烷类生物合成、异喹啉生物碱生物合成和类黄酮生物合成也被显著富集。DEGs还富集在与细胞壁组分代谢直接相关的戊糖和葡萄糖醛酸相互转化的通路上，提示SMIFH2处理可能影响细胞壁的动态平衡。此外，DEGs少量富集在ABC转运蛋白和氨酰tRNA生物合成等通路中。DEGs显著富集在果胶、多糖分解代谢的GO条目与戊糖和葡萄糖醛酸相互转化KEGG通路，体现其与细胞壁代谢相关。其中，编码果胶裂解酶（gene-9612,上调）、果胶酯酶（gene-LOC101260941,上调）和多聚半乳糖醛酸酶（gene-LOC101263946,下调）以及参与细胞壁松弛的膨胀素（EXPA5,上调）的基因均在DEGs之列。此外，DEGs富集在氧化还原酶活性、对氧化应激的反应等条目，并且多个参与次生代谢物生物合成的基因表达上调，如苯丙烷类、异喹啉生物碱和类黄酮的编码基因，这暗示了发育情况与胁迫状态的关联。图SEQ图\*ARABIC9差异表达基因的GO功能富集分析图SEQ图\*ARABIC9差异表达基因的GO功能富集分析图SEQ图\*ARABIC8差异表达基因火山图图SEQ图\*ARABIC7差异表达基因的聚类热图图SEQ图\*ARABIC图SEQ图\*ARABIC10差异表达基因的KEGG通路富集分析4讨论实验结果表明，SMIFH2处理显著抑制了番茄的主根伸长、侧根生长和根毛发育，这与成蛋白蛋白在调控肌动蛋白动态和驱动细胞生长中的核心作用相符。同时，转录组学分析揭示了SMIFH2处理主要影响与细胞壁代谢、跨膜运输以及胁迫响应相关的基因表达，为理解表型变化提供了分子层面的线索。4.1SMIFH2对番茄根系生长的抑制效应及其与肌动蛋白功能紊乱的关联在本研究中，我们观察到20µM和40µMSMIFH2处理均显著抑制了番茄幼苗的主根伸长，减少了96小时的侧根总数，并降低了根毛密度。这些表型变化与SMIFH2干扰肌动蛋白功能的结果一致。主根的生长依赖于根尖分生区细胞的持续分裂和伸长区细胞的快速、定向扩张，肌动蛋白丝在胞质环流、驱动囊泡运输以及维持细胞极性中不可或缺。SMIFH2处理可能通过破坏这些过程中肌动蛋白丝的组织与功能，导致细胞分裂受阻、细胞伸长减缓，最终表现为主根长度缩短。类似地，侧根的起始、萌发和后续生长也涉及精密的细胞分裂和扩张调控，肌动蛋白丝的紊乱同样会抑制这些过程，导致侧根数量减少。根毛作为典型的顶端生长细胞，其快速伸长依赖于高度动态的尖端肌动蛋白结构，该结构引导含有细胞壁前体和调节因子的囊泡定向运输至生长点。SMIFH2对成蛋白活性的抑制很可能破坏了根毛尖端的肌动蛋白聚焦和囊泡运输效率，从而降低了根毛的生长速率或起始频率，导致观察到的根毛密度下降。这些结果与先前在拟南芥中使用SMIFH2或其他肌动蛋白干扰药物的研究结果相似，即破坏肌动蛋白动态普遍抑制根系的生长发育[12]。4.2高浓度SMIFH2对首个侧根发生的潜在双重效应一个值得注意的发现是，虽然SMIFH2处理总体上抑制侧根数量，但40µM的高浓度处理却显著提前了首个侧根的出现时间（图4）。这一现象似乎与SMIFH2的抑制功能相悖，暗示了可能存在更复杂的调控机制。一种可能的解释是，高浓度的SMIFH2处理除了干扰正常的细胞骨架功能外，还可能触发了植物的胁迫响应。植物在逆境条件下会通过胁迫诱导的形态建成响应（Stress-InducedMorphogenicResponse,SIMR）来调整其发育程序，比如改变根系结构，通过抑制主根的生长并促进侧根的发生，以此重塑生长模式来适应性地探索新的土壤资源区域或避开胁迫源。[13]在转录组数据中，差异基因富集到了对氧化应激的反应、防御反应等GO条目以及与次生代谢物合成相关的KEGG通路（图9，10），这间接支持了SMIFH2处理可能诱导了胁迫状态。此外，激素信号的改变也可能在其中发挥作用。肌动蛋白的组织状态对生长素的响应具有重要影响。吲哚乙酸（IAA）在减少生长区长度时增加肌动蛋白束。2,4-二氯苯氧乙酸（2,4-D）在降低细胞产生率时去除肌动蛋白[14]。生长素是调控侧根起始的关键信号，在特定的根周鞘细胞中形成生长素的局部积累峰可促使侧根原基的启动[15]。高浓度SMIFH2可能在处理早期短暂地、局部地改变了根内潜伏的侧根起始位点的生长素积累模式，触发了侧根原基的快速启动，尽管后续的生长发育过程仍然受到肌动蛋白功能紊乱的抑制。当然，也不能完全排除在40µM浓度下，SMIFH2存在显著的脱靶效应。例如SMIFH2对肌球蛋白活性的影响[16]可能除成蛋白之外其他参与侧根起始调控的靶点也受到了影响。这一现象需要未来进一步深入研究以确认其原因。4.3转录组分析提示细胞壁代谢和胁迫响应通路参与SMIFH2效应40µMSMIFH2处理96小时后番茄根尖转录组的结果显示，共有64个基因发生显著表达差异（图8），虽然数量不多，但其功能富集分析为理解SMIFH2的作用机制提供了重要线索。其中DEGs显著富集在与细胞壁代谢相关的GO条目和KEGG通路。植物细胞壁既为细胞提供刚性的结构支撑又为其扩张提供可调节的弹性，其中动态的细胞骨架在控制和组织细胞壁的结构中起着重要作用[17]。SMIFH2处理干扰肌动蛋白功能，可能直接导致细胞壁合成与修饰失衡。细胞可能通过调节相关酶的基因表达来适应这种变化，亦或这些表达变化是细胞壁完整性受损的结果。无论如何，细胞壁代谢通路的显著变化与根系生长受到抑制的表型高度吻合。此外，转录组数据也反映了细胞可能处于胁迫状态。DEGs富集在氧化还原酶活性、对氧化应激的反应等条目，这些通路与植物应对氧化胁迫时的反应相关[18]。肌动蛋白骨架的破坏可能导致活性氧（reactiveoxygenspecies,ROS）信号的改变，进而触发下游的防御和代谢调整。值得注意的是，在本次分析的64个DEGs中，并未发现直接编码肌动蛋白亚基或已知番茄成蛋白的基因。这可能表明，在处理96小时这个时间点，SMIFH2对细胞骨架的主要影响可能更多地体现在蛋白质水平的功能干扰和动态变化上，而观察到的转录水平变化是下游的次级效应或适应性调整。未来的研究需要在更早的时间点进行取样，以捕捉可能存在的对细胞骨架组分基因本身的直接调控。4.4研究局限性与展望本研究首次结合表型和转录组学方法，分析了成蛋白抑制剂SMIFH2对番茄根系发育的影响，揭示了其抑制作用及相关的分子通路。然而，研究仍存在一些局限性。首先，SMIFH2作为化学抑制剂，其特异性是相对的，不能完全排除潜在的脱靶效应。因此，观察到的表型和基因表达变化是否完全由成蛋白抑制引起，还需要通过构建和分析番茄成蛋白基因突变体进行验证。其次，RNA-Seq仅在处理后96小时进行，可能无法反映SMIFH2处理早期的、更直接的分子响应。未来的研究应包含多个时间点的转录组分析，以捕捉动态变化。此外，使用整个根尖组织进行RNA提取得到的是混合细胞群的信息，可能掩盖了特定细胞类型的特异性反应，如侧根原基细胞、根毛表皮细胞。基于本研究的发现，未来的研究可以从以下几个方面深入：1）利用遗传学工具CRISPR/Cas9直接编辑番茄中的关键成蛋白基因或本次鉴定出的关键DEGs，验证它们在根系发育中的功能。2）结合活细胞成像技术，使用荧光标记的肌动蛋白探针，在活体番茄根中实时观察SMIFH2处理对肌动蛋白动态和组织结构的影响。3）通过检测内源激素水平和分布的变化，探究高浓度SMIFH2促进首个侧根发生的机制。这些研究将有助于更全面、更深入地理解成蛋白介导的肌动蛋白细胞骨架重组在番茄根系发育调控中的复杂作用网络。参考文献[1]GANDULLOJ,AHMADS,DARWISHE,etal.PhenotypingTomatoRootDevelopmentalPlasticityinResponsetoSalinityinSoilRhizotrons[J].PlantPhenomics.2021,2021.[2]BLANCAFLOREB,WANGYS,MOTESCM.Organizationandfunctionoftheactincytoskeletonindevelopingrootcells[J].IntRevCytol.2006,252:219-264.[3]COLLINGSD,LILLA,HIMMELSPACHR,etal.Hypersensitivitytocytoskeletalantagonistsdemonstratesmicrotubule-microfilamentcross-talkinthecontrolofrootelongationinArabidopsisthaliana[J].TheNewphytologist.2006,170(2):275-290.[4]HUMBERTC,AIMÉS,ALABOUVETTEC,etal.RemodellingofactincytoskeletonintomatocellsinresponsetoinoculationwithabiocontrolstrainofFusariumoxysporumincomparisontoapathogenicstrain[J].PlantPathology.2015,64(6):1366-1374.[5]CVRČKOVÁF,OULEHLOVÁD,ŽÁRSKÝV.Formins:LinkingCytoskeletonandEndomembranesinPlantCells[J].InternationalJournalofMolecularSciences.2014,16(1):1-18.[6]KOVARDR.Moleculardetailsofformin‐mediatedactinassembly[J].CurrOpinCellBiol.2006,18(1):1-7.[7]WALLARBJ,ALBERTSAS.Theformins:Activescaffoldsthatremodelthecytoskeleton[J].Trend