多版本参考基因组与严格统计流程下的番茄转录组再分析：南方番茄病毒（STV）感染仅引起乙烯相关基因轻微下调

《Discover Plants》：Reanalysis of the tomato transcriptome revealed host responses to Southern tomato virus infection

先前研究报道，感染南方番茄病毒（STV）的番茄植株表现出株高增加和果实产量提高。本研究中，研究人员旨在利用更新的生物信息学流程和最新基因组资源重新检验该原始数据集，以验证并扩展上述结果。原始单端RNA测序（RNA-Seq）读段用Rsubread重新比对至当前番茄参考基因组（SL3.0），并使用featureCounts获得基因水平计数。差异表达分析采用DESeq2，在严格阈值（错误发现率（FDR）<0.05，|log₂FC|≥1，对应倍数变化（FC）≥2）下进行。没有基因满足这些标准。与早期报道相反，宿主转录组受STV感染后基本未发生改变。在效应量最大的转录本中，乙烯生物合成基因ACO1（Solyc07g049530）的Wald log₂FC为?1.70（收缩后?1.47，校正P值=0.225），乙烯响应转录因子ERF16（Solyc12g009240）的Wald log₂FC为?1.25（收缩后?1.12，校正P值=0.083），但二者均未通过预设的Benjamini-Hochberg（BH）校正P值<0.05且|log₂FC|≥1标准。当流程针对SL4.0和SL6.0番茄参考组装运行时，结果相同。在任何参考下，零个宿主基因满足BH校正P值<0.05且|log₂FC|≥1。研究人员的RNA-Seq再分析揭示，STV感染引起番茄叶片基因表达的最小变化，最显著的是乙烯相关基因的轻微下调，而无广泛转录组重编程。这些结果缓和了先前关于广泛转录组效应的主张，反而支持STV对宿主影响轻微、可能依赖品种的观点。该研究强调更新参考基因组和严格统计方法的重要性。未来采用更大样本量或不同发育阶段的研究将具有价值。

研究背景与问题：持久性双链RNA（dsRNA）病毒可长期低滴度感染多种作物，常无显著症状，并经种子高效传播，部分研究将其视为互惠或共生体而非典型病原。STV是番茄持久性病毒，2009年由Sabanadzovic et al.首次表征，基因组为单分体dsRNA（约3.5 kb），含两个重叠开放阅读框，归入Amalgavirus属。早期在矮化、丛生和果实变色植株中检测到STV，但其也常见于无症状植株或与Pepino mosaic virus等急性病毒混合感染；其在多种番茄种质中的流行而未引发明显病害暴发，提示STV对宿主影响可能很小甚至有益。为直接评估STV影响，Fukuhara et al.在可控条件下比较M82番茄的STV感染与无STV株系，发现STV感染植株更高、果实和种子更多、种子发芽率更高，并提出STV可能赋予宿主互惠优势；叶片转录组分析还提示STV感染改变宿主基因表达，尤其是乙烯生物合成和信号通路基因下调。由于乙烯可抑制营养生长并促进果实成熟，较低乙烯活性被认为可机制性解释感染M82植株更高产的表型；原研究突出ACO1、NAC1和EREB等乙烯相关基因表达降低。然而，该转录组证据基于相对较小的RNA-Seq数据集和当时标准方法：原分析使用2012年发布的SL2.50参考基因组及TopHat/Cufflinks流程计算每千碱基每百万读段片段数（FPKM），差异表达主要依据FPKM差异并对少数基因作qPCR验证，未报告类似FDR校正的全基因组统计检验。此后番茄参考基因组显著改进（SL3.0具更完整基因模型），RNA-Seq分析转向基于读段计数的更稳健统计流程；对旧RNA-Seq数据再分析既可更有信心地确认原发现，也可能揭示旧方法造成的假阳性或假阴性。鉴于持久病毒效应细微且样本量小，假阴性和假阳性风险均较高；近期RNA-Seq分析指南建议使用旧注释再分析时进行参考敏感性检查。因此，研究人员对GSE137303数据集开展全面再分析，比较STV感染与无STV的M82植株，主要使用最新SL3.0参考，并扩展至SL4.0（Heinz 1706，ITAG4.0）和端粒到端粒SL6.0组装；同时进行无监督聚类和基因本体（GO）富集等探索性分析，以比较再分析与原发现，澄清STV感染是否重编程宿主基因表达（尤其乙烯通路），还是宿主反应极小。该问题对理解持久病毒—植物相互作用范式具有基础意义，并有助于审视分子互惠假说及植物宿主—共生体转录组研究的最佳实践。

研究开展、结论与意义：研究人员采用更新流程和最新基因组资源重新分析公共RNA-Seq数据，发现STV感染未使任何宿主基因在全基因组严格标准下显著差异表达；乙烯相关基因仅呈方向一致的轻微下调趋势，未达统计显著。该结果不支持STV引起广泛转录组重编程，提示其在M82叶片中对基因表达影响很小，可能为低影响、近中性或轻度互惠的持久病毒；同时表明更新参考基因组和严格统计方法对于重新评估宿主—病毒互作、控制假阳性具有重要意义。

主要关键技术方法（不超过250字）：数据为NCBI GEO GSE137303（Fukuhara et al., 2020）中1月龄温室M82番茄叶片6个mRNA文库（3个STV+、3个STV?；STV?为cured M82亚系，STV+自然携带STV，原研究经dsRNA提取和RT-PCR确认）。原始FASTQ经Trimmomatic质控后用Rsubread比对至SL3.0+STV（GenBank KX949574）联合参考，featureCounts作基因计数；DESeq2负二项广义线性模型、Wald检验、BH校正，阈值FDR<0.05且|log₂FC|≥1，并行SL4.0、SL6.0参考敏感性、留一法、主成分分析（PCA）与GO-BP富集。

研究结果：
3.1 测序与比对概况：六个RNA-Seq文库产生14.3–19.9百万单端读段，质控后保留超过99%；在SL3.0联合宿主+STV索引上，97.5–97.9%读段可比对，78.7–82.8%映射读段被featureCounts分配到注释基因特征；SL4.0比对率90.9–93.5%、分配率84.9–87.8%，SL6.0比对率91.2–93.5%、分配率80.4–83.3%。STV_KX949574特征读段在三个无STV文库中均为零，在三个STV感染文库中分别为386、580和304，且在不同宿主参考下计数一致，独立确认了样本感染状态及无STV株系的清洁分离。
3.2 差异基因表达分析未检出显著宿主基因：在SL3.0参考下，18,121个受检宿主基因中零个满足BH校正P值<0.05且|Wald log₂FC|≥1的预设标准；全测试集最小校正P值为0.0831，见于ERF16（Solyc12g009240）。火山图显示绝大多数基因聚集于log₂FC=0附近，无基因越过显著性阈值。预设候选基因ACO1和ERF16方向均与原报道一致，即在STV感染植株中下调，但均未通过全基因组预设标准：ACO1的Wald log₂FC为?1.70（收缩后?1.47，校正P=0.225），ERF16为?1.25（收缩后?1.12，校正P=0.083）。
3.3 参考组装敏感性确认阴性结果：在SL4.0（Heinz 1706，ITAG4.0）下，18,859个受检宿主基因中零个满足预设标准，最小校正P值为0.0874，仍见于ERF16；在端粒到端粒SL6.0组装下，17,050个受检宿主基因中零个满足标准，最小校正P值为0.0606，见于T2T特异性标识Solyc12T2T000337。由于SL6.0基因标识系统不同，Solyc07g049530/Solyc12g009240候选标识不存在于SL6.0 T2T注释，因此SL6.0仅支持全基因组层面结果。
3.4 留一法影响诊断：对SL3.0对比进行六次留一法再拟合；省略任一STV感染文库均得到零个显著宿主基因，省略STV_neg_rep2、STV_pos_rep1或STV_pos_rep2亦为零；省略STV_neg_rep3时6个基因通过标准，ACO1为最高命中；省略STV_neg_rep1时1个基因显著，为ERF16。说明主要分析的阴性结果受单一对照文库STV_neg_rep3驱动，该对照中ACO1表达异常偏低。
3.5 预设候选基因的表达趋势：虽无基因达全基因组显著，ACO1（Solyc07g049530）和ERF16（Solyc12g009240）在STV感染样本中相对无STV样本均下调，与原报道方向一致；但无STV组内变异较大，ACO1面板中一个无STV重复表达与感染组相当，降低组间分离；ERF16分离略清晰，三个感染重复均低于无STV均值，但校正P=0.083。排名靠前列表还含少量其他乙烯和防御相关转录本，但均未达全基因组显著，且因未预设、每组n=3时个体证据弱，未作为候选基因报告。
3.6 全局转录组模式未按STV状态分离：对变异最大500个基因的方差稳定化表达作PCA，PC1解释49.5%总方差，PC2解释27.8%；PC1将单个无STV样本STV_neg_rep3与其余五个样本分开，无任何轴按感染状态分离，其余无STV样本与STV感染样本聚类。样本距离热图亦未将STV感染与无STV文库分为两个明显簇，最分歧样本仍为STV_neg_rep3；最小校正P值前50个宿主特征热图同样未见STV感染组与无STV组明显分离，说明STV所致转录组差异小于植株个体间残余生物学变异。
3.7 基因本体富集分析：与微小差异表达一致，GO生物过程（BP）富集未获显著结果；研究人员分别尝试名义P<0.05基因集、名义P值最小100个基因及校正P阈值放宽至0.20的基因集，均无GO-BP条目通过FDR<0.05。最高评分条目涉及一般胁迫响应和激素代谢过程，例如“乙烯生物合成过程”（GO:0009693）和“对真菌响应”（GO:0009620）因ACO1和PR蛋白等基因出现，但富集P值约0.1且每个条目仅一两个基因贡献，统计不稳健，未解释为通路水平功能改变。

讨论部分总结：该再分析表明，STV在M82番茄叶片中引起的基因表达变化极小，符合持久病毒“ stealthy”感染、通常不触发强防御或发育扰动的特征，并强化STV在分子层面为低影响共生体的观点。结果既确认也缓和原GEO研究结论：原研究报道STV感染植株乙烯生物合成/信号基因表达降低，并以此解释更高产表型；再分析检测到ACO1、NAC1及乙烯响应转录因子等乙烯通路组分在STV+植株中较低的相同趋势，但在考虑方差和多检验校正后不显著，说明变化方向可重复而幅度有限且不确定性大。ACO1转录水平在感染植株中平均低2–3倍，与原倍数变化估计一致，但因一个无STV样本ACO1异常低至与感染植株相当，未达显著；这显示小样本可导致不同解释，留一法亦表明省略STV_neg_rep3可使6个基因跨过FDR阈值并以ACO1信号最强，说明边缘性乙烯通路结果对单一对照文库敏感，当前重复水平下轻微宿主反应难以与采样噪声区分。原分析未报告FDR或多重检验校正，其差异表达基因识别实际基于名义P值或倍数变化阈值；再分析采用更严格统计后，提示在全基因组5% FDR下无表达变化稳健。原文认为ACO1等下调可能真实但细微，需更多重复以区分噪声。若STV确实降低ACO1（ACC氧化酶）表达及可能的其他乙烯相关因子，因ACO1催化ACC转化为乙烯，ACO1降低可导致乙烯产生减少；原研究推断STV感染M82植株内源乙烯可能较低，从而减轻乙烯对生长的抑制并允许植株更高。再分析数据既不能确证也不能否定该机制，但与其一致：STV感染植株未显示乙烯通路基因升高，若有趋势则为降低；乙烯相关转录本位于最大倍数变化之列，提示乙烯信号确为STV影响的主要通路。另观察到PR蛋白和WRKY6等防御相关基因在STV感染植株中轻微下调。原研究报道少数非编码RNA和其他基因上调，但再分析中多数未出现或因低计数被过滤；例如原数据中LOC101264228和LOC101260727显示显著倍数增加，而新流程中这些转录本在所有样本中计数近零，可能被排除；注释差异（Ensembl与原ITAG2.4）可能解释部分差异。无论如何，无高表达基因或主要代谢酶因STV而差异表达，支持STV效应狭窄细微，主要触及少量调控基因而非大规模重编程。该结果与许多持久病毒近乎共生、诱导宿主胁迫可忽略的观点一致；植物持久内源RNA病毒通常无明显症状，转录组研究仅发现有限变化。普通豆内源RNA病毒感染研究报告约132个差异表达基因，富集于氧化还原过程，数量虽大于零但远少于急性病毒感染；STV—番茄体系零显著DEGs提示互作更温和或统计功效较低。STV感染植株未诱导水杨酸（SA）或茉莉酸（JA）防御标志基因，亦未激活急性病毒感染常见的RNA干扰（RNAi）基因，这与持久病毒长期进化、可能避免触发宿主防御或宿主已适应耐受/整合病毒存在的观点相符。乙烯信号下调并非所有持久植物病毒共有特征；前述豆内源RNA病毒研究未报告乙烯通路变化，而报告氧化还原和胁迫相关基因变化。原STV研究作者亦指出乙烯效应似为STV在番茄中的独特表现。乙烯广泛影响植物发育，STV轻度下调乙烯可能间接促进有利生长的激素平衡；在M82特定遗传背景中，乙烯略降可能转化为株高和果实增加。由于仅采样叶片，数据不能直接说明果实特异过程。再分析强调更新基因组资源和统计方法的价值：映射至SL3.0可减少比对错误，并获得超过80%读段分配至基因，表明转录本捕获较全面；从FPKM转向基于计数的DESeq2分析可提供误差估计并控制FDR，这对全基因组研究至关重要。“无基因满足FDR 0.05”并不否定原始观察，而提示变化处于边缘且可能来自噪声。未来增加样本量将显著提高检测细微变化的功效。对三个独立番茄参考（含端粒到端粒SL6.0）运行同一流程均产生一致阴性结果，表明STV引起的极小宿主转录组反应是生物系统属性而非特定注释所致；STV效应可能在特定阶段或组织更明显，当前数据仅为1月龄叶片，可能早于主要乙烯介导过程。研究局限在于：本研究为公共RNA-Seq数据集GSE137303的生物信息学再分析，未包含新的独立验证实验；原研究仅对少量乙烯通路基因作RT-qPCR，未覆盖此处全基因组尺度涉及的全部基因，因此将转录组结果报告为具有明确全基因组阴性结局的探索性再分析；在更多生物学重复中开展独立qRT-PCR或靶向RNA定量仍是边缘性ACO1效应的确定性检验，为最重要的后续实验。另需注意M82为常用实验室/参考品种，不同番茄基因型对STV反应可能不同，其他栽培种或野生近缘种可能反应更强或激素平衡不同；鉴于STV可经种子传播，若其确能提高产量且无负面影响，育种者或可有意维持携带STV的品系。

研究结论（译文）：对RNA-Seq数据集GSE137303的再分析揭示，STV感染未引起番茄叶片重大转录重编程，可能的例外是少数乙烯相关基因表达受抑。与许多致病病毒不同，STV未在宿主中触发广泛防御反应或胁迫通路，符合其无症状性质；STV感染植株中乙烯生物合成和信号组分轻微下调，为先前观察到的生长和产量益处提供了合理解释。STV似乎是番茄植株的近中性或轻度互惠栖居者，对宿主细胞过程负担极小。该工作证明更新基因组资源和方法在重新评估宿主—病毒相互作用中的价值，并强调生物学重复和适当统计控制对于区分真实生物学效应与噪声至关重要。未来研究建议在不同条件下考察STV影响，以确定其互惠效应是否依赖情境；直接测定有无STV时的激素水平和适合度结局以巩固因果联系亦值得开展。