《International Microbiology》：Identification of a novel mitovirus in Colocasia esculenta through high?throughput sequencing and in-silico miRNA prediction for sequence-specific gene silencing
编辑推荐：
为解决植物相关有丝分裂病毒的发现与多样性认知不足的问题，研究人员对呈现病毒症状的芋头（Colocasia esculenta）样本，开展了高通量测序分析和生物信息学预测。他们鉴定出一个暂命名为Colocasia esculenta associated mitovirus 2 (CeaMV2)的新型有丝分裂病毒，其基因组编码的RNA依赖的RNA聚合酶（RdRp）与已知同类病毒相似性低，支持其归类为Mitoviridae科新物种。此外，通过计算分析预测了植物成熟miRNA（如sbi-miR6231-3p）可靶向CeaMV2的RdRp基因，提示宿主可能通过RNA干扰（RNAi）途径参与病毒互作。该研究扩展了我们对有丝分裂病毒多样性的认识，并为利用HTS揭示隐性病毒感染和潜在RNAi防御机制提供了范例。
提到芋头，你想到的可能是软糯可口的芋泥，但你可能不知道，这种重要的热带根茎作物，也正遭受着病毒的隐秘威胁。在病毒的世界里，有一类特殊的成员——有丝分裂病毒（mitoviruses），长期以来，科学家们认为它们只在真菌的“地盘”上生活，即寄生于真菌的线粒体中，是结构最简单、功能最“精简”的RNA病毒之一。然而，随着高通量测序（High-throughput Sequencing, HTS）这把“万能钥匙”的普及，科学家们意外地打开了新世界的大门：这类病毒的“脚印”也频频出现在植物的转录组中。它们真的能跨界感染植物吗？还是仅仅作为历史的“遗迹”——内源性病毒元件（EVE）——存在于植物的基因组里？这引发了学术界对病毒宿主范围和进化历程的重新思考。
芋头作为一种重要的经济作物，其病毒感染情况尚未被完全探明。尽管此前有来自中国的研究报道了芋头中首个相关有丝分裂病毒CeaMV1，但从印度等地获得的数据依然匮乏。为了填补这一空白，揭开芋头中潜藏的病毒多样性，一个来自印度的研究团队开启了一项探索之旅。他们从印度比哈尔邦采集了出现叶缘褪绿和羽状褪绿斑驳等病毒症状的芋头叶片，试图用最前沿的技术，去“聆听”这些病叶中可能存在的病毒“低语”。这项研究最终发表在《International Microbiology》期刊上，揭示了芋头中一个全新的病毒成员，并探讨了宿主可能如何利用自身的“分子剪刀”来对抗它。
为了回答上述问题，研究人员采用了几个关键的技术方法。首先，他们对采集自印度的、表现出病毒样症状的芋头叶片样本进行总RNA提取，并利用Illumina NovaSeq 6000平台进行核糖体去除后的高通量测序。其次，他们通过生物信息学流程处理原始测序数据，去除低质量序列、与芋头参考基因组比对后，对未比对的序列进行从头组装，并通过BLASTX比对和保守结构域分析鉴定潜在的病毒序列。接着，利用RT-PCR和Sanger测序对计算分析发现的病毒序列进行了实验验证。最后，他们通过系统发育分析（使用MEGA-X软件）和序列同源性比对（使用SDT 5.4工具）来评估新病毒的进化地位，并运用三种生物信息学算法（psRNATarget, miRanda, RNA22）对已公开的植物成熟miRNA进行靶点预测，以评估宿主RNA干扰（RNAi）途径靶向该病毒的可能性。
研究结果
新型有丝分裂病毒的鉴定与特征
对病叶RNA-seq数据的分析，组装出了一条2,644个核苷酸长的序列。该序列经BLASTX比对，显示与有丝分裂病毒具有同源性，但其编码的假定RNA依赖的RNA聚合酶（RNA-dependent RNA polymerase, RdRp）蛋白与已知有丝分裂病毒的氨基酸同源性最高仅为60.76%，远低于国际病毒分类委员会（ICTV）设定的70%的物种划分阈值。因此，研究人员将其暂命名为Colocasia esculenta associated mitovirus 2 (CeaMV2)。该基因组编码一个包含544个氨基酸的单开放阅读框。Pfam数据库分析证实，在氨基酸第298-518位存在有丝分裂病毒特异的RdRp结构域（PF05919）。值得注意的是，该RdRp蛋白的长度相对较短，这与之前报道的植物相关有丝分裂病毒RdRp长度（约500个氨基酸或更短）的特征相符，而与真菌有丝分裂病毒的RdRp（约700个氨基酸或更长）不同。该序列可同时用植物标准密码子和真菌线粒体密码子进行翻译，暗示了其可能适应植物线粒体环境。
系统发育分析与序列同源性
基于RdRp氨基酸序列构建的系统发育树显示，CeaMV2与先前报道的其他植物相关有丝分裂病毒（如葡萄有丝分裂病毒1、槐叶萍有丝分裂病毒1、大麻有丝分裂病毒1等）聚为一支，共同归属于有丝分裂病毒科（Mitoviridae）。同源性分析进一步表明，CeaMV2与最接近的植物有丝分裂病毒（大麻有丝分裂病毒1）的RdRp氨基酸序列一致性也仅为40.74%，与已知的物种均存在显著差异。此外，CeaMV2与此前唯一报道的芋头有丝分裂病毒CeaMV1的序列相似性也仅为38.14%，表明两者属于不同的进化谱系。
通过RT-PCR进行验证
研究人员设计了特异性引物，通过RT-PCR对HTS的发现进行了实验验证，成功扩增出约525bp的预期条带。Sanger测序确认了该扩增产物的序列来源于CeaMV2，从而在分子水平上证实了该病毒在芋头样品中的存在。
miRNA的鉴定及其在CeaMV2基因组中的预测靶点
为探索宿主植物可能存在的抗病毒RNAi机制，研究人员进行了计算分析。他们利用三种预测算法，从公共数据库的植物成熟miRNA库中筛选可能靶向CeaMV2 RdRp基因的序列。分析发现，多个植物来源的成熟miRNA可能靶向RdRp基因，其中包括sbi-miR6231-3p、mtr-miR5759、hbr-miR6171、gra-miR7484c和bdi-miR7746-3p等家族成员。其中，sbi-miR6231-3p展现出最高的结合亲和力（解链目标位点所需的能量UPE为13.912）。这一预测结果暗示，宿主植物的miRNA可能通过RNAi途径，靶向并沉默病毒的RdRp基因，从而构成了潜在的抗病毒防御层。
结论与讨论
本研究在印度的芋头中鉴定并描述了一个新型的有丝分裂病毒CeaMV2。其基因组特征、显著低的RdRp序列同源性（22.55%-40.74%，远低于物种阈值）以及系统发育树中与植物相关有丝分裂病毒聚类的结果，共同支持将其归类为有丝分裂病毒科内的一个新物种。CeaMV2的发现，与来自中国的CeaMV1的报告一起，凸显了有丝分裂病毒在重要经济作物芋头中的地理多样性，也进一步证实了植物可以作为这类病毒的宿主，拓展了我们对有丝分裂病毒宿主范围和多样性的认知。其相对较短的RdRp长度、可被植物和真菌两种线粒体密码子翻译的特性，都指向了其对植物线粒体环境的潜在适应性，为理解这类病毒在动植物界间的进化与跨宿主适应提供了新的素材。
研究的一个重要延伸在于，通过计算预测揭示了植物内源性miRNA（尤其是sbi-miR6231-3p）靶向CeaMV2 RdRp基因的可能性。这提示了在芋头与CeaMV2的相互作用中，可能存在由宿主miRNA介导的RNA干扰防御机制。这不仅是理解植物-病毒互作分子机制的一个切入点，也为未来开发基于人工miRNA（amiRNA）的基因工程策略，培育抗此类病毒作物品种提供了潜在的靶点和理论依据。然而，该研究也指出了当前的局限性，例如，检测到的病毒序列可能来源于整合到宿主基因组中的内源性病毒元件，而非活动性病毒感染；观察到的叶片症状也可能由其他生物胁迫因素引起，而非CeaMV2直接导致。因此，未来需要通过PCR排除内源化、亚细胞定位、真菌微生物组分析以及5‘和3’RACE（cDNA末端快速扩增）等技术，来进一步明确CeaMV2的宿主特异性、基因组完整性及其确切的生物学作用。
总而言之，这项研究不仅利用高通量测序成功“捕捞”到一个隐藏的病毒新成员，丰富了我们对病毒多样性的认知图谱，还通过计算预测，为理解宿主如何利用“沉默武器”对抗病毒入侵提供了线索。它再次证明了高通量测序技术在揭示隐性病毒感染方面的强大威力，并为从基础研究通向作物病毒病可持续防控的应用探索，架起了一座充满潜力的桥梁。
