《Frontiers in Cellular and Infection Microbiology》：RNA virome comparison between sylvatic and urban-interface mosquitoes from Southeastern Brazil
编辑推荐：
引言：蚊类（双翅目：蚊科 Culicidae）是高影响力公共卫生虫媒病毒（arbovirus）的主要媒介，但其核心病毒组尤其在新生热带（Neotropical）林区谱系中仍缺乏系统表征。本研究调查巴西南圣保罗州自城市至森林梯度下多种蚊种（含被忽视的林区分类单元
引言：蚊类（双翅目：蚊科 Culicidae）是高影响力公共卫生虫媒病毒（arbovirus）的主要媒介，但其核心病毒组尤其在新生热带（Neotropical）林区谱系中仍缺乏系统表征。本研究调查巴西南圣保罗州自城市至森林梯度下多种蚊种（含被忽视的林区分类单元 Sabethes、Psorophora、Shannoniana 及 Wyeomyia）的RNA病毒组。方法：研究人员在圣保罗州城市至森林梯度采集多蚊种，开展RNA病毒组分析，并进行生态分析以评估宿主分类与环境对病毒组构成的作用，通过网络分析探究病毒-宿主关联及跨生态界面的病毒共享。结果：分析鉴定出919条病毒重叠群（contig），隶属217个病毒种、37个病毒科，在主要以林区蚊种中揭示大量低氨基酸一致性（中位<40%）的“病毒暗物质”（viral dark matter）。虽检出含已知虫媒病毒的 Flaviviridae、Phenuiviridae 与 Peribunyaviridae 等科，但在采样与测序深度敏感限度内未发现高后果人类病原体。生态分析显示病毒组构成受宿主分类与环境强烈结构化（R2=0.570，p=0.001），宿主种解释32.9%独特方差（PERMANOVA，R2=0.329，p=0.001），生境（ecotope）居次（R2=0.029，p=0.001）。网络分析揭示以宿主限制性专化病毒为主的高度模块化病毒-宿主结构，少量桥接种促进跨生态界面病毒共享。讨论：结果表明蚊类内在生物学是病毒群落结构主驱因素，环境梯度居次，并凸显宿主关联过程在新热带森林-城市界面病毒多样性塑造中的重要性。
研究背景与立项依据：蚊科（Culicidae）昆虫维持并传播登革、寨卡、基孔肯雅与黄热等高影响公共卫生虫媒病毒，但医学相关病原仅占蚊类关联病毒多样性一小部分，完整病毒集合即病毒组（virome）还包含昆虫特异性病毒（insect-specific viruses, ISVs）、植物/真菌关联病毒及大量未分类RNA病毒。巴西高强度城市化与高生物多样性并存，圣保罗州集城市、城郊与保存林区于一体，是考察生态梯度如何结构化蚊媒病毒群的天然场景。现有新热带区域研究多忽略 Sabethes、Psorophora、Shannoniana、Wyeomyia 等林区属，病毒组是否被随机组装、宿主演化史与生态上下文如何共同约束病毒群落仍不清楚，故研究人员开展此次跨城市-森林界面多蚊种RNA病毒组与生态网络研究。论文发表于《Frontiers in Cellular and Infection Microbiology》。
主要技术方法概要：研究人员于圣保罗州10个市（Barbosa、Bebedouro、Castilho、Cardoso、Guarulhos、Igaratá、Mat?o、Pindamonhangaba、S?o Paulo、Santa Rita do Passa Quatro）森林（含林缘）与城市-林区界面（动物园内）两生态背景，按地面与林冠两层采集1604头蚊，鉴定为17种、归为45个单种池。总RNA提取后建Illumina MiSeq双端150 bp文库；Raw reads经FastQC质检、Trimmomatic修剪、Bowtie2比对蚊科参考基因组去除宿主读段、SortMeRNA去rRNA，RnaviralSPAdes进行de novo组装（保留>500 bp contig）；Kraken2初分、BLASTn/x/p确证分类并用taxonomizr标准化 lineage，BLASTn筛查并剔除内源性病毒元件（endogenous viral elements, EVEs，identity>85%、cover>50%）；清洁读段回贴量化（均深≥5×、覆盖≥90%留用）。生态上构建样本-病毒矩阵、rarefy后用Bray–Curtis距离做PCoA，adonis2（PERMANOVA）与betadisper（PERMANOVA-based dispersion）分层检验生境/分类/空间/技术协变量，igraph/tidygraph构建蚊-病毒二部网络并以ggraph力导向布局可视化。
研究结果
3.1 宏基因组概览：组装得103714条contig，仅1.0%归病毒、20.51%未定性；筛除反转录转座子/EVE/噬菌体后保留919条病毒contig（217种），87种（40.1%）出现于≥2池。Mika virus与Culex Iflavi-like virus最常检出（各20/19池）。37个病毒科中 Flaviviridae、Phenuiviridae、Peribunyaviridae 等医学相关科存在；Partitiviridae、Iflaviridae、Solemoviridae、Nodaviridae 合计占47%且跨12、12、9、9种蚊，呈广宿主格局；Closteroviridae、Qinviridae 仅见于 Psorophora varipes，Endornaviridae、Solinviviridae 仅 Sabethes aurescens，Polycipiviridae、Artoviridae、Xinmoviridae 仅 Mansonia humeralis，Caliciviridae、Mesoniviridae、Totiviridae 仅 Culex spp.，Lispiviridae 仅 Culex chidesteri，证明强宿主专化。
3.2 病毒组架构的生态决定因素：PCoA（Bray–Curtis）前两轴释23.3%方差，呈宿主驱动结构；Shannon α多样性在林区与城市-林区界面间无差（p=0.89），但 Sabethes spp.与 Shannoniana fluviatilis 等林区专化种病毒丰富度显著更高。全模型PERMANOVA（R2=0.570，F=1.985，p=0.001）中宿主分类与市为主预测因子；宿主种释32.9%（R2=0.329，p=0.001），属级17.1%、族级13.4%，市级在种级模型仅2.9%但族级达14.7%；生境在含市变量全模型中被掩盖，在简约模型中具次要显著效应。多数方差（68.5%）未解释，指向微环境或随机驱动。
3.3 病毒共享的网络枢纽：二部网络呈模块化星型，多数病毒度=1（单宿主），Culex（尤 Cx. chidesteri、Cx. quinquefasciatus）为最高度枢纽；Aedes scapularis 与 Sabethes albiprivus 作桥接种连多分类单元，Aedes aenigmaticus 与 Sabethes identicus 成孤立卫星模块。空间网络显示 S?o Paulo–Guarulhos 共享最高（至多8种），Mantel检验地理距离与病毒共享/相异无关（r=-0.02，p=0.535；r=0.017，p=0.461），表明远距扩散由人为廊道与宿主移动而非连续空间扩散驱动。
讨论浓缩：研究人员指出，虽检出 Phenuiviridae、Flaviviridae、Rhabdoviridae、Peribunyaviridae 等医学相关科，但敏感限度内无高后果虫媒病毒，可能源于常住ISV通过超感染排斥（superinfection exclusion）或固有免疫预激活构成生物屏障降低媒介效能，亦受测序深度限制。宿主种级解释力跃升印证“物种特异屏障”，垂直冠层无差异进一步支持内在生物学压倒生境垂直分区。桥接种（Aedes scapularis、Sabethes albiprivus）与新热带林区“病毒暗物质”（AAI<40%）提示森林界面是新颖RNA病毒储备库；Culex 关联病毒高模块化暗示ISV封闭储库。局限含市内/外采样季节异步（2019-10至2020-03 vs 2020-05至12）及机会性抽样。
结论译文：本研究大幅拓展新热带蚊类病毒多样性认知，揭示尤以林区属为主的未表征病毒储库。病毒组强宿主驱动结构证明蚊类内在生物学为病毒群落构主决因素，环境与空间因素居次。以宿主限制性专化病毒为主的模块化病毒-宿主网络凸显约束跨种病毒共享的生态与系统发育屏障。尽管检出若干医学相关病毒科但未现高后果虫媒病毒，与本土病毒组可能构筑抵御外源感染生物屏障的假说一致。结果强调整合宏基因组与生态框架对理解病毒动态及预测快速变化环境中虫媒病毒涌现的价值。
