病毒侵染对番茄营养器官miRNAs及靶标表达的影响研究

病毒侵染对番茄营养器官miRNAs及靶标表达的影响研究一、引言1.1研究背景植物病毒病作为农作物生产的重要威胁之一，对全球农业经济造成了巨大损失。据统计，每年因植物病毒病导致的农作物减产高达20%-40%，部分严重病害甚至可造成绝收。植物病毒通过干扰植物正常的生理代谢过程，影响植物的生长发育、光合作用和物质运输等，进而降低作物的产量和品质。在众多受病毒威胁的农作物中，番茄作为一种重要的蔬菜作物，在全球范围内广泛种植，深受消费者喜爱，其种植面积和产量均居世界前列。然而，番茄极易受到多种病毒的侵染，引发番茄病毒病。番茄病毒病已成为制约番茄产业发展的关键因素之一，给番茄种植户带来了沉重的经济负担。番茄病毒病的种类繁多，常见的有番茄黄化曲叶病毒病、番茄褪绿病毒病、番茄褐色皱纹果病毒病等。这些病毒病的症状各异，对番茄植株的危害程度也不尽相同。例如，番茄黄化曲叶病毒病主要表现为植株矮化、生长点褪绿黄化、扭曲，叶片边缘上卷、增厚、变硬，结果数减少，果实变小且着色不均匀；番茄褪绿病毒病则使整株表现褪绿黄化，顶叶呈鲜黄色，从第二片叶开始呈斑驳黄化，越往下斑驳症状越明显，同时叶片变厚，果实小、颜色偏白、不能正常膨大；番茄褐色皱纹果病毒病会导致植株嫩叶和顶芽出现花叶症状，茎秆或萼片出现坏死条斑，叶片变窄，早期花萼褐变，严重时花萼干枯甚至死亡，果实染病后变形、颜色变浅、出现黄色斑点，随着病情加重，黄色斑点发展成褐色病变，最终导致果实坏死，对产量影响较大。这些病毒病不仅影响番茄的外观品质，还降低了其营养价值和商品价值，给番茄产业带来了巨大的经济损失。微小核糖核酸（microRNAs，miRNAs）是一类长度约为22个核苷酸的非编码小分子RNA，在植物的生长发育、生理代谢和逆境响应等过程中发挥着重要的调控作用。miRNAs通过与靶基因mRNA的互补配对，介导靶基因的降解或翻译抑制，从而实现对基因表达的精细调控。在植物的防御机制中，miRNAs也扮演着至关重要的角色。当植物受到病毒侵染时，miRNAs能够识别并结合病毒的核酸序列，引发RNA干扰（RNAi）反应，降解病毒的RNA，从而抑制病毒的复制和传播。同时，miRNAs还可以通过调控植物自身的防御相关基因的表达，增强植物的免疫力，抵御病毒的入侵。例如，miR168能够靶向调控AGO1基因的表达，AGO1是RNAi途径中的关键蛋白，参与植物对病毒的防御反应；miR482则可以调控NBS-LRR家族抗病基因的表达，增强植物对病毒的抗性。研究病毒侵染与miRNAs之间的关系，对于深入理解植物的抗病毒机制具有重要意义。通过揭示miRNAs在植物抗病毒过程中的作用机制，可以为开发新型的抗病毒策略提供理论依据。例如，利用基因编辑技术调控miRNAs的表达，或者设计人工miRNAs来靶向病毒基因，有望提高植物的抗病毒能力。此外，研究病毒与miRNAs的互作关系，还可以为番茄病毒病的早期诊断和预警提供新的方法和技术手段。通过检测番茄植株中miRNAs的表达变化，能够及时发现病毒的侵染，采取相应的防治措施，减少病毒病的传播和扩散，降低经济损失。因此，开展病毒侵染对番茄营养器官中特定miRNAs及其靶标表达影响的研究，具有重要的理论意义和实际应用价值，有助于推动番茄产业的可持续发展，保障农产品的质量安全和供应稳定。1.2研究目的与意义本研究旨在深入揭示病毒侵染对番茄营养器官中特定miRNAs及其靶标表达的影响，通过系统分析不同病毒侵染下番茄体内miRNAs和靶标基因的表达变化规律，明确其在番茄抗病毒防御反应中的调控机制。具体而言，通过高通量测序技术和生物信息学分析，筛选出受病毒侵染显著调控的miRNAs及其靶标基因，并运用实时荧光定量PCR、原位杂交、基因编辑等技术对其表达模式和功能进行验证，深入探究miRNAs与靶标基因之间的相互作用关系，以及它们在番茄生长发育和抗病毒过程中的生物学功能。本研究具有重要的理论意义和实际应用价值。在理论方面，深入研究病毒侵染与miRNAs及其靶标表达的关系，有助于揭示植物抗病毒的分子机制，丰富和完善植物免疫学理论体系，为深入理解植物与病毒之间的相互作用提供新的视角和理论依据。同时，研究miRNAs在番茄生长发育和逆境响应中的调控作用，也有助于进一步阐明植物生长发育的分子调控网络，为植物基因功能研究提供新的思路和方法。在实际应用方面，本研究结果可为番茄病毒病的防治提供新的策略和方法。通过调控miRNAs的表达或干预miRNAs与靶标基因的相互作用，可以开发出新型的抗病毒基因工程番茄品种，提高番茄的抗病毒能力，减少化学农药的使用，降低生产成本，保障农产品的质量安全。此外，研究结果还可为番茄病毒病的早期诊断和预警提供新的分子标记和技术手段，通过检测番茄植株中特定miRNAs的表达变化，能够及时发现病毒的侵染，采取相应的防治措施，有效控制病毒病的传播和扩散，减少经济损失，促进番茄产业的可持续发展。1.3国内外研究现状在植物病毒病研究领域，病毒侵染对植物miRNAs表达的影响已成为国内外学者关注的热点。国内外研究表明，病毒侵染能够引发植物体内miRNAs表达的显著变化，这些变化在植物抗病毒防御反应中发挥着重要作用。国外学者在该领域开展了大量研究。例如，Bao等利用miRNA芯片技术分析了大豆花叶病毒（SMV）易感大豆品系中SMV响应的miRNA，鉴定出数十种差异表达的miRNAs，并通过后续实验验证了部分miRNAs对靶基因表达的调控作用，发现SMV可以通过诱导相应miRNA的积累，下调几种RNAi通路基因和NBS-LRR家族抗性基因来阻碍大豆防御反应。此外，还有研究聚焦于其他植物病毒与miRNAs的互作关系，如在烟草中，马铃薯X病毒侵染会导致本氏烟草miRNAs表达谱发生变化，部分差异表达的miRNAs对马铃薯X病毒侵染及症状表现产生影响。这些研究为深入理解病毒侵染对植物miRNAs表达的影响提供了重要的理论依据和研究思路。国内学者在该领域也取得了一系列重要成果。福建农林大学吴建国教授团队针对水稻草矮病毒（RGSV）和水稻锯齿叶矮缩病毒（RRSV）开展研究，筛选了23个对病毒侵染有响应的非保守miRNA，利用过量表达、短串联靶点模拟（STTM）和CRISPR/Cas9技术，构建了功能增强、抑制和缺失的转基因水稻资源库，并通过“集团接种法”对转基因水稻材料进行抗病毒评估，鉴定出多个miRNA在水稻抗病毒中发挥重要作用，其中miR535显著削弱了水稻对RGSV的抗性，而miR1868.1则显著增强了水稻对RRSV的抗性。扬州大学李良俊教授课题组通过对莲藕节间转录组和SmallRNA测序的联合组学分析，揭示了microRNA调控淀粉合成的复杂网络，为深入理解miRNAs在植物生长发育中的调控机制提供了新的视角。尽管国内外在病毒侵染影响植物miRNAs表达方面取得了一定进展，但仍存在一些不足和空白。现有研究多集中于少数模式植物和特定病毒，对于番茄等重要经济作物中病毒侵染与miRNAs及其靶标表达关系的研究相对较少，且缺乏系统性和全面性。在研究方法上，虽然高通量测序技术和生物信息学分析为miRNAs的鉴定和功能预测提供了有力手段，但仍需要进一步结合分子生物学实验技术，如实时荧光定量PCR、原位杂交、基因编辑等，对miRNAs的表达模式和功能进行深入验证，以明确其在植物抗病毒防御反应中的具体作用机制。此外，目前对于miRNAs与靶标基因之间的相互作用关系以及它们在植物生长发育和抗病毒过程中的调控网络研究还不够深入，有待进一步加强。综上所述，深入研究病毒侵染对番茄营养器官中特定miRNAs及其靶标表达的影响，不仅可以填补现有研究的空白，完善植物与病毒互作的分子机制，还能为番茄病毒病的防治提供新的理论依据和技术策略，具有重要的科学意义和应用价值。二、材料与方法2.1实验材料实验选用的番茄品种为“中杂9号”，该品种是中国农业科学院蔬菜花卉研究所培育的鲜食番茄一代杂种，具有生长势强、坐果率高、果实大、品质好等特点，在番茄种植中广泛应用，对多种病毒具有一定的抗性，适合用于本实验研究病毒侵染对其miRNAs及靶标表达的影响。实验所用的病毒种类包括番茄黄化曲叶病毒（Tomatoyellowleafcurlvirus，TYLCV）、番茄褪绿病毒（Tomatochlorosisvirus，ToCV）和番茄褐色皱纹果病毒（Tomatobrownrugosefruitvirus，ToBRFV）。TYLCV属于双生病毒科菜豆金色花叶病毒属，是一种单链环状DNA病毒，其基因组大小约为2.7kb，主要通过烟粉虱以持久循环方式传播，在世界热带和亚热带地区广泛分布，严重影响番茄的产量和品质，侵染后会导致番茄植株矮化、生长点褪绿黄化、扭曲，叶片边缘上卷、增厚、变硬，结果数减少，果实变小且着色不均匀。ToCV属于长线形病毒科毛形病毒属，为正义单链RNA病毒，基因组大小约为9.4kb，由粉虱以半持久方式传播，在全球番茄种植区均有发生，感染ToCV的番茄植株会出现整株褪绿黄化，顶叶呈鲜黄色，从第二片叶开始呈斑驳黄化，越往下斑驳症状越明显，同时叶片变厚，果实小、颜色偏白、不能正常膨大。ToBRFV属于烟草花叶病毒属，是一种单链正义RNA病毒，基因组全长约6.4kb，主要通过汁液摩擦传播，近年来在全球多地番茄产区爆发，给番茄产业带来巨大损失，感染ToBRFV的番茄植株嫩叶和顶芽出现花叶症状，茎秆或萼片出现坏死条斑，叶片变窄，早期花萼褐变，严重时花萼干枯甚至死亡，果实染病后变形、颜色变浅、出现黄色斑点，随着病情加重，黄色斑点发展成褐色病变，最终导致果实坏死。本实验中所用的病毒均由中国农业科学院植物保护研究所病毒研究室提供，病毒保存于-80℃冰箱中，在使用前需进行复苏和扩繁。实验所需的主要试剂包括TRIzol试剂（Invitrogen公司）用于提取番茄组织中的总RNA；PrimeScriptRTreagentKitwithgDNAEraser（TaKaRa公司）用于反转录合成cDNA；SYBRPremixExTaqII（TaKaRa公司）用于实时荧光定量PCR反应；RNAisoPlus（TaKaRa公司）用于分离小RNA；miRcutemiRNAFirst-StrandcDNASynthesisKit（TIANGEN公司）用于miRNA的反转录；DNAMarker、限制性内切酶、T4DNA连接酶等（NEB公司）用于分子克隆实验；琼脂糖、丙烯酰胺、甲叉双丙烯酰胺等（Sigma公司）用于电泳实验；DEPC水、氯仿、异丙醇、无水乙醇等（国药集团化学试剂有限公司）用于常规试剂配制。实验所需的主要仪器设备包括高速冷冻离心机（Eppendorf公司）用于离心分离样品；实时荧光定量PCR仪（ABI公司）用于检测基因表达量；凝胶成像系统（Bio-Rad公司）用于观察和记录电泳结果；PCR扩增仪（Bio-Rad公司）用于进行PCR反应；核酸蛋白分析仪（ThermoScientific公司）用于测定核酸浓度和纯度；超净工作台（苏净集团）用于无菌操作；恒温培养箱（上海一恒科学仪器有限公司）用于培养番茄植株；光照培养箱（宁波江南仪器厂）用于提供适宜的光照和温度条件；低温冰箱（海尔集团）用于保存试剂和样品。2.2实验方法2.2.1番茄种植与病毒接种番茄种植于人工气候室内，温度控制在25±2℃，光照周期为16h光照/8h黑暗，光照强度为300μmol・m⁻²・s⁻¹，相对湿度保持在60%-70%。采用基质栽培方式，栽培基质为草炭、蛭石和珍珠岩按3:1:1的体积比混合而成，每立方米基质中添加10kg有机肥和2kg复合肥，充分混匀后装入栽培槽中。挑选饱满、无病虫害的番茄种子，用55℃温水浸泡15min进行消毒处理，然后用清水冲洗干净，再用蒸馏水浸泡4-6h，待种子吸胀后，置于28℃恒温培养箱中催芽，待种子露白后，播种于育苗盘中，每穴播1粒种子，覆盖1cm厚的基质，浇透水，置于人工气候室内培养。待番茄幼苗长至3-4片真叶时，选择生长健壮、大小一致的幼苗进行移栽，移栽时将幼苗从育苗盘中取出，轻轻抖掉根部的基质，避免损伤根系，然后将幼苗移栽到栽培槽中，每槽种植2株，株距为30cm，移栽后及时浇定根水。当番茄幼苗长至6-7片真叶时进行病毒接种。对于TYLCV，采用农杆菌介导的接种方法。将含有TYLCV基因组的农杆菌菌株GV3101在含有相应抗生素的LB液体培养基中振荡培养至OD₆₀₀值为0.6-0.8，然后将菌液离心收集，用重悬液（10mmol/LMgCl₂、10mmol/LMES、200μmol/L乙酰丁香酮）重悬至OD₆₀₀值为1.0，在室温下静置3h后，用注射器将菌液注射到番茄植株的叶片中，每株接种4-5片叶，接种量为每片叶0.1mL。对于ToCV，采用粉虱介导的接种方法。将携带ToCV的粉虱在含有番茄植株的养虫笼中饥饿2-3h，然后将粉虱转移到健康的番茄植株上，每株接种20-30头粉虱，让粉虱在植株上取食48h后，将粉虱移除。对于ToBRFV，采用汁液摩擦接种方法。将感染ToBRFV的番茄叶片研磨成匀浆，用磷酸缓冲液（0.01mol/L，pH7.0）稀释10倍，然后用棉球蘸取稀释后的汁液，在番茄植株的叶片上轻轻摩擦，每株接种4-5片叶，接种后用清水冲洗叶片，避免汁液残留。接种病毒后，将番茄植株继续培养在人工气候室内，定期观察植株的生长状况和发病症状。2.2.2样品采集与处理分别在病毒接种后的3d、5d、7d、10d、15d采集番茄植株的叶片、茎和根等营养器官样品。每个时间点选取3株生长状况相似的番茄植株，每株植株采集相同部位的叶片、茎和根样品各0.5g，迅速放入液氮中速冻，然后转移至-80℃冰箱中保存备用。在样品处理过程中，将采集的样品从-80℃冰箱中取出，置于冰上解冻，然后用剪刀将样品剪成小块，放入预冷的研钵中，加入适量的液氮，迅速研磨成粉末状。将研磨好的样品粉末转移至1.5mL离心管中，加入1mLTRIzol试剂，剧烈振荡混匀，室温静置5min，使样品充分裂解。然后加入0.2mL氯仿，盖紧管盖，剧烈振荡15s，室温静置3min，4℃下12000rpm离心15min。离心后，将上层水相转移至新的1.5mL离心管中，加入等体积的异丙醇，混匀后室温静置10min，4℃下12000rpm离心10min，此时离心管底部会出现白色的RNA沉淀。小心弃去上清液，加入1mL75%乙醇（用DEPC水配制），轻轻洗涤RNA沉淀，4℃下7500rpm离心5min，弃去上清液，将RNA沉淀在室温下晾干5-10min，待乙醇挥发完全后，加入适量的DEPC水溶解RNA沉淀，用核酸蛋白分析仪测定RNA的浓度和纯度，将RNA样品保存于-80℃冰箱中备用。2.2.3miRNAs及其靶标检测方法采用实时荧光定量PCR（qRT-PCR）技术检测miRNAs及其靶标mRNA的表达量。首先，使用miRcutemiRNAFirst-StrandcDNASynthesisKit（TIANGEN公司）将总RNA中的miRNAs反转录成cDNA。反应体系为10μL，包括5×miRcuteReactionBuffer2μL、miRcuteRTEnzymeMix0.5μL、TotalRNA1μg、RNaseFreeddH₂O补至10μL。反应条件为：42℃孵育60min，85℃孵育5min，然后将反应产物保存于-20℃冰箱中备用。对于靶标mRNA的反转录，使用PrimeScriptRTreagentKitwithgDNAEraser（TaKaRa公司），反应体系为20μL，包括5×PrimeScriptBuffer24μL、PrimeScriptRTEnzymeMixI1μL、OligodTPrimer（50μM）1μL、Random6mers（100μM）1μL、TotalRNA1μg、RNaseFreeddH₂O补至20μL。反应条件为：37℃孵育15min，85℃孵育5s，然后将反应产物保存于-20℃冰箱中备用。qRT-PCR反应使用SYBRPremixExTaqII（TaKaRa公司），反应体系为20μL，包括SYBRPremixExTaqII10μL、ForwardPrimer（10μM）0.8μL、ReversePrimer（10μM）0.8μL、cDNA模板2μL、ddH₂O6.4μL。反应在实时荧光定量PCR仪上进行，反应条件为：95℃预变性30s，然后进行40个循环，每个循环包括95℃变性5s，60℃退火30s，在退火阶段采集荧光信号。每个样品设置3个生物学重复和3个技术重复。以番茄的U6snRNA作为内参基因，用于校正miRNAs的表达量；以番茄的Actin基因作为内参基因，用于校正靶标mRNA的表达量。采用2⁻ΔΔCt法计算miRNAs及其靶标mRNA的相对表达量，公式为：相对表达量=2⁻（ΔCt样品-ΔCt内参），其中ΔCt=Ct目的基因-Ct内参基因。对qRT-PCR实验数据进行统计分析，采用SPSS22.0软件进行方差分析（One-wayANOVA），并用Duncan氏新复极差法进行多重比较，P<0.05表示差异显著，以确定不同病毒侵染下番茄营养器官中miRNAs及其靶标mRNA表达量的变化情况。三、黄瓜花叶病毒（CMV）侵染对番茄营养器官miRNAs及靶标表达的影响3.1CMV国际标准株系CMV-Fny侵染结果利用实时荧光定量PCR技术，对CMV-Fny侵染后的番茄根、茎组织中特定miRNAs及其靶标mRNAs的表达量进行了精确检测。结果显示，在根组织中，miR156a的表达量在侵染后呈现出先上升后下降的趋势。在侵染后3d，miR156a的表达量相较于对照组显著上调，达到对照组的2.5倍（P<0.05），这表明miR156a可能在番茄根组织应对CMV-Fny侵染的早期阶段发挥重要作用，通过调控相关靶标基因的表达，参与植物的防御反应。随着侵染时间的延长，在侵染后15d，miR156a的表达量逐渐下降至接近对照组水平，这可能是由于植物在病毒侵染后期，通过其他调控机制来维持自身的生理平衡，减少对miR156a介导的调控途径的依赖。与之对应的是，其靶标基因SPL9的表达量则呈现出相反的变化趋势。在侵染后3d，SPL9的表达量显著下调，仅为对照组的0.4倍（P<0.05），这与miR156a的上调表达相呼应，表明miR156a通过与SPL9mRNA的互补配对，介导了SPL9基因的降解或翻译抑制，从而调控其表达量。随着侵染时间的推移，在侵染后15d，SPL9的表达量逐渐回升，但仍显著低于对照组水平（P<0.05），这说明即使在侵染后期，miR156a对SPL9的调控作用依然存在，且病毒侵染对SPL9基因的表达产生了持续的影响。在茎组织中，miR164的表达量在CMV-Fny侵染后表现出持续上升的趋势。在侵染后3d，miR164的表达量相较于对照组上调了1.8倍（P<0.05），且在侵染后15d，其表达量进一步升高，达到对照组的3.2倍（P<0.05），这表明miR164在番茄茎组织应对病毒侵染的过程中持续发挥作用，可能通过调控其靶标基因的表达，参与茎组织的防御反应和生理调节。而miR164的靶标基因NAC1的表达量则在侵染后持续下降。在侵染后3d，NAC1的表达量显著下调，为对照组的0.6倍（P<0.05），随着侵染时间的延长，在侵染后15d，NAC1的表达量进一步降低，仅为对照组的0.3倍（P<0.05），这与miR164的上调表达密切相关，进一步证实了miR164对NAC1基因的负调控作用，即miR164通过抑制NAC1基因的表达，影响茎组织的生长发育和对病毒侵染的响应。同时，对CMV-Fny病毒基因的表达量进行检测，发现病毒基因在番茄根、茎组织中的表达量随着侵染时间的延长而逐渐增加。在侵染后3d，病毒基因在根组织中的表达量相对较低，但在侵染后15d，其表达量显著升高，达到侵染后3d的10倍（P<0.05）；在茎组织中，病毒基因的表达量在侵染后3d同样较低，随着侵染时间的推进，在侵染后15d，其表达量升高至侵染后3d的8倍（P<0.05）。通过对miRNAs及其靶标mRNAs表达量与病毒基因表达量之间的相关性分析，发现miR156a与SPL9的表达量变化与病毒基因在根组织中的表达量呈显著负相关（r=-0.85，P<0.01），这表明在根组织中，随着病毒基因表达量的增加，miR156a的表达量上升，进而抑制SPL9的表达，以应对病毒的侵染。而miR164与NAC1的表达量变化与病毒基因在茎组织中的表达量呈显著负相关（r=-0.88，P<0.01），说明在茎组织中，病毒基因表达量的增加促使miR164表达上调，从而抑制NAC1的表达，参与茎组织对病毒侵染的防御反应。3.2致病因子2b缺失突变体CMV-FnyΔ2b侵染结果当番茄植株受到致病因子2b缺失突变体CMV-FnyΔ2b侵染后，通过实时荧光定量PCR技术对其根、茎组织中特定miRNAs及其靶标mRNAs的表达量展开检测，结果呈现出与CMV-Fny侵染不同的变化趋势。在根组织中，miR159的表达量在侵染初期迅速上升，在侵染后3d时，相较于对照组上调了3.2倍（P<0.05），达到峰值，随后逐渐下降，但在侵染后15d时，仍显著高于对照组水平（P<0.05），表明miR159在番茄根组织应对CMV-FnyΔ2b侵染的早期阶段发挥了重要作用，可能通过调控相关靶标基因的表达，参与植物的防御反应。其靶标基因MYB33的表达量则在侵染后呈现出持续下降的趋势，在侵染后3d，MYB33的表达量显著下调，仅为对照组的0.3倍（P<0.05），随着侵染时间的延长，在侵染后15d，其表达量进一步降低至对照组的0.1倍（P<0.05），这与miR159的上调表达密切相关，说明miR159通过抑制MYB33基因的表达，影响根组织的生长发育和对病毒侵染的响应。在茎组织中，miR396的表达量在CMV-FnyΔ2b侵染后表现出先上升后下降的趋势。在侵染后5d，miR396的表达量相较于对照组显著上调，达到对照组的2.8倍（P<0.05），随后逐渐下降，在侵染后15d时，其表达量接近对照组水平。这表明miR396在番茄茎组织应对病毒侵染的过程中，在特定时期发挥了作用，可能参与了茎组织对病毒侵染的早期防御反应。而miR396的靶标基因GRF4的表达量则在侵染后呈现出先下降后上升的趋势。在侵染后5d，GRF4的表达量显著下调，为对照组的0.4倍（P<0.05），随着侵染时间的推移，在侵染后15d，GRF4的表达量逐渐回升至接近对照组水平，这与miR396的表达变化趋势相互呼应，进一步证实了miR396对GRF4基因的负调控作用，即miR396通过抑制GRF4基因的表达，参与茎组织对病毒侵染的响应。同时，对CMV-FnyΔ2b病毒基因的表达量进行检测，发现病毒基因在番茄根、茎组织中的表达量始终维持在较低水平。在侵染后3d，病毒基因在根组织中的表达量仅为CMV-Fny侵染时的0.01倍（P<0.05），在茎组织中的表达量为CMV-Fny侵染时的0.02倍（P<0.05）；在侵染后15d，病毒基因在根组织中的表达量略有上升，但仍显著低于CMV-Fny侵染时的水平（P<0.05），在茎组织中的表达量同样显著低于CMV-Fny侵染时的水平（P<0.05）。通过对miRNAs及其靶标mRNAs表达量与病毒基因表达量之间的相关性分析，发现miR159与MYB33的表达量变化与病毒基因在根组织中的表达量呈显著负相关（r=-0.92，P<0.01），表明在根组织中，随着病毒基因表达量的降低，miR159的表达量上升，进而抑制MYB33的表达，以应对病毒的侵染。而miR396与GRF4的表达量变化与病毒基因在茎组织中的表达量在侵染早期呈显著负相关（r=-0.89，P<0.01），说明在茎组织侵染早期，病毒基因表达量的降低促使miR396表达上调，从而抑制GRF4的表达，参与茎组织对病毒侵染的防御反应，但在侵染后期，这种相关性减弱，可能是由于其他调控机制的参与。3.3添加致强卫星sat335的CMV-Fny.satTl侵染结果在添加致强卫星sat335的CMV-Fny.satTl侵染番茄的实验中，通过实时荧光定量PCR技术对番茄各营养器官中病毒基因的积累量进行检测，结果显示，在叶组织中，CMV-Fny各基因(亚)组的积累量直接降低，在侵染后15d，其表达量仅为对照组的0.3倍（P<0.05），这表明卫星sat335的存在对病毒在叶组织中的复制和积累产生了明显的抑制作用。在茎组织中，病毒基因的积累量呈现出先升后降的趋势。在侵染后7d，病毒基因的表达量达到峰值，为对照组的1.5倍（P<0.05），随后逐渐下降，在侵染后15d，其表达量降至对照组的0.6倍（P<0.05），这种变化可能与病毒粒子在茎组织中的运输及扩散过程有关，初期病毒在茎组织中大量繁殖，随着时间推移，卫星sat335的作用逐渐显现，抑制了病毒的进一步积累。在根系组织中，病毒基因表达量下降出现延迟，在侵染后10d，其表达量仍与对照组无显著差异（P>0.05），直至侵染后15d，才显著下降至对照组的0.5倍（P<0.05），说明根系组织对病毒侵染的响应相对较慢，可能存在一些特殊的防御机制或生理过程影响了病毒在根系中的复制和扩散。对番茄营养器官中miRNAs及其靶标表达的检测发现，在根系发育过程中，CMV-Fny.satTl对番茄侧根发育造成了严重影响。NAC1作为与侧根发育密切相关的基因，其相对表达量显著降低，在侵染后15d，仅为对照组的0.2倍（P<0.05），同时，侧根数量也大幅度减少，相较于对照组减少了70%（P<0.05），这表明CMV-Fny.satTl侵染通过抑制NAC1基因的表达，阻碍了番茄侧根的正常发育。在茎组织的发育过程中，接种21天后的CMV-Fny.satTl处理组中茎木髓部出现了明显异常，表现为茎中空。进一步检测发现，该处理组番茄茎中miR164的含量偏高，在侵染后21d，其表达量为对照组的2.8倍（P<0.05），而靶标NAC1的含量明显偏低，仅为对照组的0.3倍（P<0.05），说明miR164对NAC1的负调控作用在茎组织发育异常过程中起到了重要作用，miR164的高表达抑制了NAC1的表达，进而影响了茎组织的正常发育，导致茎中空症状的出现。通过对miRNAs及其靶标mRNAs表达量与病毒基因表达量之间的相关性分析，发现miR164与NAC1的表达量变化与病毒基因在茎组织中的表达量在侵染前期呈显著正相关（r=0.87，P<0.01），在侵染后期呈显著负相关（r=-0.90，P<0.01）。在侵染前期，随着病毒基因表达量的增加，miR164的表达量也随之上升，可能是植物对病毒侵染的一种早期响应机制，通过上调miR164的表达来抑制NAC1基因的表达，以应对病毒的侵染。然而，在侵染后期，病毒基因表达量的降低并没有导致miR164表达量的相应下降，反而使得miR164对NAC1的抑制作用更加明显，这可能是由于病毒侵染引发了一系列复杂的生理变化，导致miR164的调控机制发生了改变，或者是其他因素的参与使得miR164在侵染后期对NAC1的调控作用更加关键，从而影响了茎组织的发育和病毒的侵染进程。四、番茄不孕病毒（TAV-Bj）侵染对番茄营养器官miRNAs及靶标表达的影响4.1TAV-Bj侵染后番茄营养器官症状表现在番茄植株被番茄不孕病毒（TAV-Bj）侵染后的早期阶段，根系生长就受到了显著抑制。从接种后第5天开始，与对照组相比，受侵染植株的侧根发育明显迟缓，侧根原基的分化数量减少，根系形态变得稀疏且短小。在接种后第10天，这种差异更加显著，侧根数量相较于对照组减少了约40%，且侧根长度仅为对照组的60%，根系整体呈现出营养不良、发育不良的状态。通过对根系解剖结构的观察发现，受侵染根系的皮层细胞排列紊乱，细胞间隙增大，部分细胞出现坏死现象，维管束组织发育异常，导管数量减少且管径变小，这严重影响了根系对水分和养分的吸收与运输功能。在茎组织方面，随着侵染时间的延长，TAV-Bj对番茄茎的影响逐渐显现。在接种后第15天，茎基部开始出现轻微的缢缩现象，茎表皮颜色变深，呈现出暗绿色。随着病情的发展，在接种后第20天，茎组织的病变进一步加重，茎秆变得脆弱易折断，髓部组织出现中空现象，且在茎的纵切面可以观察到明显的褐色坏死条纹，这些坏死条纹沿着维管束组织分布，导致维管束功能受损，影响了植株体内物质的运输和分配。同时，受侵染植株茎部的节间长度明显缩短，相较于对照组缩短了30%，使得植株整体呈现出矮化、紧凑的生长形态。在叶片上，受TAV-Bj侵染的番茄植株叶片在接种后第7天就开始出现症状。叶片边缘逐渐向上卷曲，叶片表面出现不规则的黄色斑驳，这些斑驳随着时间的推移逐渐扩大并融合。在接种后第15天，叶片卷曲程度加剧，部分叶片甚至呈筒状，黄色斑驳占据了叶片面积的50%以上，且叶片质地变硬、变脆，光合作用面积大幅减少，严重影响了植株的光合作用效率，进而影响植株的生长发育和产量形成。4.2TAV-Bj侵染对miRNAs及靶标表达量的影响为了深入探究TAV-Bj侵染对番茄营养器官中miRNAs及其靶标表达量的影响，采用实时荧光定量PCR技术，对TAV-Bj侵染后的番茄根、茎组织中特定miRNAs及其靶标mRNAs的表达量进行了精准检测。结果显示，在根组织中，miR167的表达量在侵染后呈现出显著的变化趋势。在侵染后3d，miR167的表达量相较于对照组显著上调，达到对照组的3.0倍（P<0.05），这表明miR167在番茄根组织应对TAV-Bj侵染的早期阶段可能发挥着重要的调控作用，通过调控相关靶标基因的表达，参与植物的防御反应。随着侵染时间的延长，在侵染后15d，miR167的表达量虽有所下降，但仍显著高于对照组水平（P<0.05），说明miR167在根组织中的调控作用在病毒侵染后期依然持续存在。与之对应的是，其靶标基因ARF8的表达量则呈现出明显的下调趋势。在侵染后3d，ARF8的表达量显著下调，仅为对照组的0.3倍（P<0.05），这与miR167的上调表达密切相关，表明miR167通过与ARF8mRNA的互补配对，介导了ARF8基因的降解或翻译抑制，从而抑制其表达量。随着侵染时间的推移，在侵染后15d，ARF8的表达量进一步降低，仅为对照组的0.1倍（P<0.05），这进一步证实了miR167对ARF8基因的负调控作用，且这种调控作用在病毒侵染过程中逐渐增强，对根组织的生长发育和对病毒侵染的响应产生了重要影响。在茎组织中，miR172的表达量在TAV-Bj侵染后表现出先上升后下降的趋势。在侵染后5d，miR172的表达量相较于对照组显著上调，达到对照组的2.5倍（P<0.05），达到峰值，随后逐渐下降，在侵染后15d时，其表达量接近对照组水平。这表明miR172在番茄茎组织应对病毒侵染的过程中，在特定时期发挥了关键作用，可能参与了茎组织对病毒侵染的早期防御反应，通过调控相关靶标基因的表达，来抵御病毒的入侵。而miR172的靶标基因AP2的表达量则在侵染后呈现出先下降后上升的趋势。在侵染后5d，AP2的表达量显著下调，为对照组的0.4倍（P<0.05），随着侵染时间的推移，在侵染后15d，AP2的表达量逐渐回升至接近对照组水平，这与miR172的表达变化趋势相互呼应，进一步证实了miR172对AP2基因的负调控作用，即miR172通过抑制AP2基因的表达，参与茎组织对病毒侵染的响应，在侵染后期，随着miR172表达量的下降，AP2基因的表达逐渐恢复，可能是植物自身的一种调节机制，以维持茎组织的正常生长发育。同时，对TAV-Bj病毒基因的表达量进行检测，发现病毒基因在番茄根、茎组织中的表达量随着侵染时间的延长而逐渐增加。在侵染后3d，病毒基因在根组织中的表达量相对较低，但在侵染后15d，其表达量显著升高，达到侵染后3d的8倍（P<0.05）；在茎组织中，病毒基因的表达量在侵染后3d同样较低，随着侵染时间的推进，在侵染后15d，其表达量升高至侵染后3d的6倍（P<0.05）。通过对miRNAs及其靶标mRNAs表达量与病毒基因表达量之间的相关性分析，发现miR167与ARF8的表达量变化与病毒基因在根组织中的表达量呈显著负相关（r=-0.90，P<0.01），这表明在根组织中，随着病毒基因表达量的增加，miR167的表达量上升，进而抑制ARF8的表达，以应对病毒的侵染。而miR172与AP2的表达量变化与病毒基因在茎组织中的表达量在侵染早期呈显著负相关（r=-0.87，P<0.01），说明在茎组织侵染早期，病毒基因表达量的增加促使miR172表达上调，从而抑制AP2的表达，参与茎组织对病毒侵染的防御反应，但在侵染后期，这种相关性减弱，可能是由于其他调控机制的参与，使得茎组织对病毒侵染的响应更加复杂。五、不同病毒侵染影响的对比分析5.1不同病毒对番茄营养器官发育影响对比在对不同病毒侵染番茄的研究中，发现黄瓜花叶病毒（CMV）的不同株系以及番茄不孕病毒（TAV-Bj）对番茄营养器官发育有着显著且各异的影响。在根系发育方面，CMV国际标准株系CMV-Fny和致病因子2b缺失突变体CMV-FnyΔ2b的侵染对番茄侧根发育的影响相对较小。然而，添加致强卫星sat335的CMV-Fny.satTl和TAV-Bj却对番茄侧根发育造成了严重破坏。CMV-Fny.satTl侵染后，与侧根发育密切相关的基因NAC1的相对表达量显著降低，在侵染后15d，仅为对照组的0.2倍（P<0.05），同时，侧根数量相较于对照组大幅度减少了70%（P<0.05）。TAV-Bj侵染后的番茄植株，从接种后第5天起侧根发育就明显迟缓，接种后第10天，侧根数量相较于对照组减少了约40%，侧根长度仅为对照组的60%。这种差异表明，病毒的组成成分，如卫星RNA的存在，以及病毒自身的特性，会显著影响其对番茄侧根发育的干扰程度。在主根发育上，CMV-FnyΔ2b处理组表现出独特的现象，其主根发育最长。这可能是由于2b基因的缺失，改变了病毒与寄主植物之间的相互作用方式，使得寄主植物在应对病毒侵染时，对主根发育的调控机制发生了改变，从而促进了主根的生长。在茎组织发育过程中，不同病毒的影响也各有特点。接种21天后的CMV-Fny.satTl处理组中，茎木髓部出现明显异常，表现为茎中空。进一步研究发现，该处理组番茄茎中miR164的含量偏高，在侵染后21d，其表达量为对照组的2.8倍（P<0.05），而靶标NAC1的含量明显偏低，仅为对照组的0.3倍（P<0.05），表明miR164对NAC1的负调控作用在茎组织发育异常中起到了关键作用。TAV-Bj侵染后的番茄植株，在接种后第15天，茎基部开始出现轻微缢缩现象，表皮颜色变深，随着病情发展，茎秆变得脆弱易折断，髓部组织出现中空现象，且茎的纵切面可见明显褐色坏死条纹，节间长度相较于对照组缩短了30%，植株呈现矮化、紧凑的生长形态。这说明不同病毒对茎组织的影响不仅体现在形态结构上，还涉及到内部细胞结构和生理功能的改变，且这些变化与病毒的致病性密切相关。综合来看，不同病毒对番茄营养器官发育的影响存在显著差异，这些差异与病毒的种类、组成成分以及致病性密切相关。病毒通过干扰番茄体内的基因表达调控网络，特别是miRNAs及其靶标的表达，来影响营养器官的正常发育，进而导致植株出现各种病变症状，严重影响番茄的生长和产量。5.2不同病毒对miRNAs及靶标表达影响对比不同病毒侵染番茄后，对其营养器官中miRNAs及靶标表达产生了显著且各具特征的影响。在黄瓜花叶病毒（CMV）相关研究中，国际标准株系CMV-Fny侵染后，番茄根组织中miR156a表达先升后降，其靶标基因SPL9表达则持续下调；茎组织中miR164表达持续上升，靶标基因NAC1表达持续下降。致病因子2b缺失突变体CMV-FnyΔ2b侵染时，根组织中miR159初期快速上升后下降，靶标基因MYB33持续下调；茎组织中miR396先升后降，靶标基因GRF4先降后升。添加致强卫星sat335的CMV-Fny.satTl侵染时，病毒基因在不同营养器官中的积累量变化不同，且影响了番茄侧根发育和茎组织发育，导致NAC1表达降低，侧根数量减少，茎木髓部出现中空，同时茎中miR164含量偏高，靶标NAC1含量偏低。番茄不孕病毒（TAV-Bj）侵染后，根组织中miR167表达显著上调，靶标基因ARF8表达显著下调；茎组织中miR172先升后降，靶标基因AP2先降后升。通过对比可以发现，不同病毒侵染导致的miRNAs及靶标表达变化模式存在差异。CMV-Fny和CMV-FnyΔ2b对番茄侧根发育影响较小，而CMV-Fny.satTl和TAV-Bj则对侧根发育造成严重破坏，表现为NAC1相对表达量降低和侧根数量大幅减少。在茎组织中，CMV-Fny.satTl和TAV-Bj侵染后，茎中miR164含量均偏高，靶标NAC1含量明显偏低，均出现茎部异常症状，如茎中空、茎基部缢缩、髓部组织中空等。这些结果表明，病毒种类与miRNAs调控途径之间存在紧密联系。不同病毒可能通过特异性地干扰番茄体内特定miRNAs及其靶标的表达，来实现对寄主植物生理过程的调控，进而引发不同的病害症状。病毒的致病性、组成成分等因素会影响其对miRNAs及靶标表达的调控模式，如CMV中2b基因的缺失以及卫星RNA的存在，均改变了病毒与寄主植物之间的相互作用方式，导致miRNAs及靶标表达发生不同变化。此外，不同病毒侵染后，miRNAs及靶标表达变化的时间进程和幅度也有所不同，反映了病毒与寄主植物在互作过程中的复杂性和特异性。深入研究这种关系，有助于进一步揭示植物抗病毒的分子机制，为番茄病毒病的防治提供理论基础。5.3病毒致病性与miRNAs调控的关联分析通过对不同病毒侵染番茄后营养器官发育及miRNAs、靶标表达的研究，深入分析病毒致病性与miRNAs调控之间的关联，有助于揭示病毒致病的分子机制。黄瓜花叶病毒（CMV）的不同株系以及番茄不孕病毒（TAV-Bj）对番茄营养器官发育和miRNAs及靶标表达产生了不同影响，这些影响与病毒的致病性密切相关。在根系发育方面，CMV-Fny和CMV-FnyΔ2b对番茄侧根发育影响较小，而CMV-Fny.satTl和TAV-Bj则对番茄侧根发育造成严重破坏，导致NAC1相对表达量降低和侧根数量大幅度减少。这表明病毒的致病性越强，对番茄侧根发育的干扰越严重，而这种干扰可能是通过调控与侧根发育相关的miRNAs及其靶标来实现的。例如，CMV-Fny.satTl侵染后，NAC1基因的表达受到抑制，可能是由于其调控的miRNA表达发生变化，从而影响了侧根的正常发育。在茎组织发育中，CMV-Fny.satTl和TAV-Bj侵染后，茎中均出现异常症状，且miR164含量偏高，靶标NAC1含量明显偏低。这说明病毒的侵染改变了miR164对NAC1的调控平衡，进而影响茎组织的正常发育。高致病性的病毒可能通过干扰miR164的表达，使其对NAC1的抑制作用增强，导致茎组织发育异常，如茎中空、茎基部缢缩等症状的出现。从miRNAs及靶标表达的变化来看，不同病毒侵染后，miRNAs及其靶标的表达模式与病毒致病性紧密相关。CMV-Fny侵染后，番茄根组织中miR156a和茎组织中miR164的表达变化与病毒基因表达量呈显著负相关，表明这些miRNAs在病毒侵染过程中发挥了防御作用，通过调控靶标基因的表达来抑制病毒的复制和传播。而TAV-Bj侵染后，根组织中miR167和茎组织中miR172的表达变化在侵染早期与病毒基因表达量呈显著负相关，后期相关性减弱，这可能是由于病毒在侵染后期采用了其他调控机制来逃避植物的防御反应，或者植物自身的防御系统发生了适应性变化。综合以上分析，病毒致病性与miRNAs调控之间存在密切关联。病毒通过干扰番茄体内miRNAs及其靶标的表达，打破植物正常的生长发育调控网络，从而导致病害症状的产生。病毒的致病性越强，对miRNAs调控途