《Journal of Virology》：Expanding the Orchidaceae virome through analysis of public-domain sequencing data
摘要
兰科（Orchidaceae）是最大的植物科之一，包含众多具有重要观赏、经济和生态价值的物种。然而，尽管其多样性极高，与兰科相关的病毒却鲜有深入研究。高通量测序数据的日益丰富为在广泛的分类尺度上探索植物病毒多样性提供了新的机遇。本研究通过对2789个公开可用的RNA-seq数据集进行大规模荟萃分析，这些数据代表了18个亚族和5个亚科的41个属，从而对兰科相关病毒进行了系统性研究。研究人员鉴定出对应于163种病毒的序列，其中包括13种与兰科的新关联。值得注意的是，其中117条序列很可能代表推定的新物种，包括在甜菜坏死黄脉病毒科（Benyviridae）、β-弯曲病毒科（Betaflexiviridae）、芬纽病毒科（Phenuiviridae）和黄化病毒科（Tymoviridae）内的暂定新属。该分析还揭示了兰科病毒组内病毒共发生的模式，并发现了病毒类群与兰科生活型之间的关联。总体而言，这些发现极大地拓宽了已知兰科相关病毒的多样性，并为了解其组成、生态互作和潜在范围提供了一个全面的全球框架。
研究背景与科学问题
兰科（Orchidaceae）是全球最大的植物科之一，拥有超过800个属和30000余个公认物种，广泛分布于世界各地。许多兰科植物因其极高的观赏价值而被广泛栽培，例如卡特兰属（Cattleya）、兰属（Cymbidium）、石斛属（Dendrobium）、兜兰属（Paphiopedilum）和蝴蝶兰属（Phalaenopsis）在全球花卉市场中占据重要份额。此外，香荚兰（Vanilla spp.）作为最具经济价值的兰科产品之一，是最昂贵的香料之一。除了经济价值，兰科植物在维持传粉者多样性、通过与菌根真菌共生贡献土壤微生物多样性以及作为环境变化指示生物方面也发挥着关键的生态作用。
尽管兰科植物种类繁多且生态意义重大，但迄今报道的可感染兰科的病毒数量相对较少。最普遍且造成严重经济损失的病原体是建兰花叶病毒（Cymbidium mosaic virus, CymMV）和齿舌兰环斑病毒（Odontoglossum ringspot virus, ORSV）。其他已报道的病毒还包括兰花花叶病毒（Orchid fleck virus, OFV）和香荚兰花叶病毒（Vanilla mosaic virus, VanMV）等。鉴于兰科植物的巨大多样性，已记录的病毒数量表明兰科病毒组在很大程度上尚未被充分探索。全球病毒多样性调查指出，未被发现的病毒数量可能远超现有描述的分类群。近期在澳大利亚、欧洲和亚洲的小规模测序调查已从少量兰科样本中发现了多种此前未描述的病毒。同时，病毒也是对兰科保护的一个额外威胁，因为约30%的兰科物种正面临灭绝风险。保护计划常依赖超低温保存和离体繁殖，但若缺乏可靠的诊断工具和对兰科相关病毒的全面了解，受感染的繁殖材料可能无意中将病毒引入自然生态系统。因此，获取病毒基因组数据对于开发灵敏的检测方法至关重要。
为解决这一知识空白，研究人员利用日益增长的公共高通量测序数据资源，首次对兰科相关病毒进行了大规模的全球性评估。研究人员整合了已发表的病毒记录，并对公开可用的RNA-seq数据集进行了系统性挖掘，旨在全面表征与兰科植物相关的病毒多样性、分布及其进化关系。
主要技术方法
- 1.
数据获取与预处理：从美国国家生物技术信息中心（NCBI）的序列读取档案库（Sequence Read Archive, SRA）下载了2789个兰科RNA-seq数据集，涵盖203个兰科物种和杂交种。使用FastQC评估原始读段质量，并使用Seqtk进行修剪。每个数据集独立使用rnaSPAdes进行从头组装。
- 2.
病毒序列鉴定：构建了一个包含所有NCBI RefSeq病毒蛋白、兰科注释蛋白及Ensembl非冗余蛋白的参考数据库。使用DIAMOND BLASTX对组装的重叠群（contigs）进行筛选，采用严格的90%最低比对覆盖度阈值以减少假阳性。随后通过BLASTn和BLASTx比对NCBI非冗余（nr）数据库进行验证，并剔除映射到非病毒蛋白或基因组的序列。
- 3.
序列注释与系统发育分析：通过Magic-BLAST回帖读段和IGV手动检查验证病毒样序列的组装质量。利用BLASTn和DIAMOND计算检出的病毒序列与参考序列间的成对核苷酸或氨基酸序列一致性。依据国际病毒分类委员会（International Committee on Taxonomy of Viruses, ICTV）的划界标准进行初步分类。使用IQ-TREE基于核苷酸或氨基酸比对进行最大似然法系统发育分析，并通过ModelFinder自动选择替代模型。
- 4.
关联分析与多样性估计：使用费希尔精确检验（Fisher's exact test）评估病毒属和科在不同兰科属和族中的富集情况。利用NetworkX构建病毒-宿主关联网络，并使用Louvain算法识别网络社群。通过稀疏分析（Rarefaction analysis）拟合负指数和双曲线模型，估算潜在的病毒物种总丰富度。
研究结果
全球结果描述
通过筛选2789个RNA-seq数据集，研究人员鉴定出163种推定病毒，与114个兰科物种或杂交种相关，跨越41个属、18个亚族和5个亚科。其中117种很可能代表新病毒物种，超过了研究当时已知的96种兰科相关病毒。平均每种兰科物种检出1.6种病毒。费希尔精确检验显示，Phalaenopsis（病毒比率为0.42）和Cymbidium（0.86）属的病毒数量低于平均水平。病毒多样性最高的是Paphiopedilum属（21种病毒），其次是Cypripedium属（10种病毒，来自7个科）。病毒科级别的富集分析显示，杓兰亚科（Cypripedioideae）的属在Secoviridae和Partitiviridae中富集；红门兰亚科（Orchidoideae）的属在Partitiviridae中富集；而树兰亚科（Epidendroideae）则在Virgaviridae和Alphaflexiviridae中显著富集。在族水平上，头蕊兰族（Neottieae）在Benyviridae中富集，而沼兰族（Malaxideae）在Potyviridae中富集。
Alphaflexiviridae
鉴定了四个成员，包括CymMV、仙人掌X病毒（Cactus virus X, CVX）、鹤顶兰X病毒（Phaius virus X, PhVX）和香荚兰潜隐病毒（Vanilla latent virus, VLV）。CymMV最为普遍，在217个数据集中检出，涉及15个属的36个兰科物种或杂交种。系统发育分析显示CymMV存在两个分支。CVX在兰科中的检出尚属首次报道。
Amalgaviridae
在Ophrys iricolor中确认了Ophrys fusca amalgavirus 1（OfAV1），并在Paphiopedilum armeniacum中发现了一个可能代表新种的序列，暂命名为Paphiopedilum armeniacum-associated amalgavirus 1（PaaAV1）。
Aspiviridae
确认了数种兰科感染的蛇形病毒（ophioviruses），并鉴定出两个推定新种，分别暂命名为platanthera-associated ophiovirus 1（PgaOV1）和calanthe-associated ophiovirus 1（CalanOV1）。
Benyviridae
在Dendrobium officinale中意外检测到甜菜坏死黄脉病毒（Beet necrotic yellow vein virus, BNYVV）的RNA片段，但因BNYVV寄主范围通常限于双子叶植物，此发现需谨慎对待。此外，鉴定出六个与Dactylorhiza hatagirea beny-like virus（DhBLV）密切相关的beny-like病毒，它们形成一个独立于Benyvirus属的独特分支，可能代表Benyviridae科内的一个新谱系。这些病毒基因组编码一个大的RdRp蛋白和一个假定的膜蛋白，但未检测到运动蛋白和外壳蛋白。
Betaflexiviridae
发现了与兰科的新关联及推定新属。在马铃薯H病毒（Potato virus H, PotVH）中首次报告了与兰科的关联。在Capillovirus属中鉴定出三个推定新种。在Divavirus属中发现了一个新种。更重要的是，发现了一个与Divavirus属姐妹的分支，其基因组特征独特，可能代表Betaflexiviridae科下的一个新属。另一个由四种新病毒构成的分支与Prunevirus和Citrivirus属姐妹，但其基因组更小，可能也代表一个新属。
Bromoviridae
黄瓜花叶病毒（Cucumber mosaic virus, CMV）在Cymbidium sinense和Dendrobium moniliforme中被新检出，其分离株遗传差异较大，暗示可能源自独立的跨物种传播事件。
Caulimoviridae
从Mexipedium xerophyticum、Dendrobium catenatum和Paphiopedilum armeniacum的转录组中恢复了七个候选badnavirus-like序列。由于转录组证据无法区分游离病毒和内源副逆转录病毒元件（EPRVs），这些序列被视为候选关联，有待实验验证。
Endornaviridae
鉴定出六个新的alphaendornavirus-like序列，与多个兰科物种相关。由于其也可能源自组织样本中的真菌或卵菌，其真正宿主有待澄清。
Geminiviridae
在Phalaenopsis equestris中鉴定出一个与鹰嘴豆红叶病毒2（Chickpea redleaf virus 2, CpRLV2）高度相似的序列。在Calanthe appendiculata中检测到一个与番茄黄化曲叶病毒（Tomato yellow leaf curl virus, TYLCV）几乎完全一致的片段，但TYLCV通常不感染单子叶植物，此发现同样需谨慎对待。
Kitaviridae
鉴定出两种新的kitavirus，分别属于Cilevirus属（暂命名vandeae-associated cilevirus 1, VaCV1）和Blunervirus属（暂命名dendrobium officinale-associated blunervirus 1, DoaBV1）。
Partitiviridae
重点关注了Deltapartitivirus属。确认了已知病毒的存在，并鉴定出72个推定新种。根据RdRp氨基酸相似性（90%阈值），这些新种被分为三个主要分支和若干亚支。
Phenuiviridae
首次在兰科中发现该科病毒。鉴定出四个rubodvirus-like序列，它们在系统发育树上形成独立于Rubodvirus属的独特分支，可能代表Phenuiviridae科下此前未描述的属。
Potyviridae
检出10种马铃薯Y病毒（potyviruses），代表六个进化谱系，包括四个推定新种。获得了calanthe mild mosaic virus（CalMMV）的首个近乎完整基因组。发现了一个可能由Dendrobium相关病毒构成的、此前未被识别的兰科相关马铃薯Y病毒谱系。首次报告了东亚西番莲病毒（East Asian passiflora virus, EAPV）与兰科的关联。
Rhabdoviridae
鉴定出分属于Betanucleorhabdovirus、Cytorhabdovirus和Varicosavirus属的序列。发现了数个推定新种，如cypripedium betacytorhabdovirus 2（CypBCRV2）和orchis italica-associated varicosavirus-like 1（OiVL1）。
Secoviridae
在Fabavirus属中，首次报告了葡萄蚕豆萎蔫病毒（Grapevine fabavirus, GFabV）和蚕豆萎蔫病毒2（Broad bean wilt virus 2, BBWV2）与兰科的关联。在Stralavirus属中发现了一个新病毒，暂命名cremastra appendiculata stralavirus（CaSLV）。
Solemoviridae
在Dendrobium williamsonii中发现了一个可能属于Polerovirus属的新种，暂命名dendrobium williamsonii-associated polerovirus 1（DWaPV1）。
Tombusviridae
在Dendrobium nobile中检测到香石竹斑驳病毒（Carnation mottle virus, CarMV），将其与一个新的兰科物种关联起来。
Tospoviridae
仅在Oncidium sp.中检测到辣椒褪绿病毒（Capsicum chlorosis virus, CaCV）。
Tymoviridae
首次在兰科中发现该科病毒。鉴定出两个tymo-like病毒，分别与Epidendrum puniceoluteum和Cremastra appendiculata相关。
Virgaviridae
ORSV是Virgaviridae科中最常检出的病毒。在Gastrodia elata中发现了一个与烟草脆裂病毒（Tobacco rattle virus, TRV）相近的推定新tobravirus。在Calanthe lyroglossa中发现了一个新的tobamovirus。
病毒关联分析
通过结合本研究与NCBI Virus数据库的数据，兰科病毒组共有约213种病毒与70个兰科属相关。大多数病毒（89.9%）仅在一个属中检出，而CymMV和ORSV在最多属中检出。网络分析将共享病毒检出的兰科属划分为四个社群，分别以Paphiopedilum（CymMV主导）、ORSV、Phalaenopsis/Oncidium（OFV主导）和Diuris（BYMV主导）为核心。病毒共发生分析显示，Alphaflexiviridae与Virgaviridae之间存在最强的正关联，而Alphaflexiviridae与Partitiviridae之间存在显著的负关联。聚类分析还揭示了附生兰和地生兰的病毒群落组成存在差异。
兰科相关病毒的预计多样性
稀疏分析预测兰科病毒物种总数可能在1032至1980种之间。但由于目前仅有约40个兰科属有RNA-seq数据，采样严重不足，该预测应视为探索性结果，主要说明了当前未记录的病毒多样性的潜在规模，而非一个稳健的估计值。
讨论与结论
本研究极大地扩展了对兰科相关病毒多样性的认识，并表明兰科病毒组仍远未被完全揭示。研究人员强调，由于研究依赖于公开的转录组数据，检出的病毒序列代表与兰科组织相关的病毒RNA，而非直接证据表明发生了活跃感染，序列可能来源于植物病毒、组织中其他生物的病毒或环境来源。尽管采用了严格过滤策略，但对于与兰科生物学特性不符的病毒关联（如BNYVV和TYLCV）需谨慎解读。
这些发现对兰科保护具有重要意义。研究观察到CymMV和ORSV的高流行率，这些广寄主范围的病毒可通过种质和亲本材料污染带来风险。此外，33种兰科相关病毒（15.5%）也在兰科以外的植物中被检出，表明作物与兰科间可能存在交叉传播。
研究人员指出，病毒与宿主的相互作用是一个连续谱，从致病到潜在有益。这突显了进一步表征新发现病毒生物学特性的必要性，特别是Benyviridae、Betaflexiviridae、Phenuiviridae和Tymoviridae科中的推定新属。例如，beny-like病毒可能与根部寄生性传播媒介有关，这与相关兰科的地生习性一致。而隐症病毒在兰科中的作用也值得关注，它们常建立持续的无症状感染，可能影响宿主防御和胁迫耐受性。附生兰与地生兰病毒组的差异也提示宿主生态学可能塑造病毒群落结构。
最终，植物病毒对农业和自然生态系统构成重大威胁，其出现预计会随着全球贸易、人口流动和气候变化而增加。这项工作展示了大规模数据挖掘作为一种成本效益高的策略在探究植物病毒多样性方面的强大潜力，尤其是在资源有限、难以实施大规模先进技术调查的地区。总之，本研究结果突显了整合基因组分析在拓展我们对植物病毒多样性及其关联理解方面的力量。
