《Genomics & Informatics》：Exploratory strain-associated patterns of antiviral transcriptional responses to Zika virus exposure in developing human neural tissue
编辑推荐：
摘要：寨病毒（Zika virus, ZIKV）感染影响早期人脑发育，但不同发育细胞谱系如何塑造抗病毒转录响应尚不完全清楚。本研究通过对公开的人胎儿脑来源神经组织经巴西ZIKV-BR株、柬埔寨2010年亚洲谱系FSS13025分离株（简称ZIKV-FSS）或干
摘要：寨病毒（Zika virus, ZIKV）感染影响早期人脑发育，但不同发育细胞谱系如何塑造抗病毒转录响应尚不完全清楚。本研究通过对公开的人胎儿脑来源神经组织经巴西ZIKV-BR株、柬埔寨2010年亚洲谱系FSS13025分离株（简称ZIKV-FSS）或干扰素β（interferon-beta, IFNβ）体外（ex vivo）暴露的单细胞转录组数据集进行二次描述性分析来探讨此问题。在各神经谱系中，富集祖细胞和增殖神经祖细胞（neural progenitor cells, NPCs）的放射状胶质样细胞显示出强烈的转录水平I型干扰素/干扰素刺激基因（interferon-stimulated genes, ISGs）响应性，且较其他数种细胞类型满足探索性差异表达标准的基因数目更多。探索性比较提示存在毒株关联的响应结构差异——在选定的祖细胞聚焦及ISG相关汇总中，ZIKV-FSS关联的模式与IFNβ相关谱更相似，而ZIKV-BR暴露则关联更广泛和多样的基因表达变化。敏感性分析显示ZIKV与IFNβ间的转录相似性依赖于所用基因全集（gene universe），在经典ISGs中最强，在更广泛的非ISG及所有表达基因集合中更为保守。对独立人脑类器官（cerebral organoid）bulk RNA-seq数据集及Vero细胞/单核细胞来源树突状细胞（moDC）单细胞数据集的背景分析支持ZIKV暴露系统中普遍存在IFN/ISG相关转录激活，但未将其解释为胎儿神经谱系特异性响应的验证。总体而言，结果提供了发育中人神经组织中抗病毒转录状态的探索性、谱系分辨描述性参考。由于主要胎儿神经数据集每条件仅含一个生物学重复，这些发现应被解读为探索性毒株关联模式而非毒株依赖机制的确定性证据。
论文解读：发育中人神经组织ZIKV暴露后抗病毒转录响应的毒株关联谱系解析二次分析
研究背景与立题依据
寨病毒（Zika virus, ZIKV）感染与先天性寨病毒综合征（congenital Zika syndrome, CZS）密切相关，可导致小头畸形及皮质发育异常，但其为何特异地影响发育中人脑及其前体细胞群体尚未完全阐明。已有研究表明放射状胶质（radial glia, RG）样细胞及增殖神经祖细胞（neural progenitor cells, NPCs）对ZIKV易感，且可强烈诱导干扰素刺激基因（interferon-stimulated genes, ISGs）程序，但这种ISG诱导是代表细胞自主抗病毒启动还是旁分泌I型干扰素（type I interferon, IFN-I）信号尚存争议。此外，不同ZIKV毒株在神经细胞感染性、病毒产量及先天免疫基因诱导上存在表型差异，而单细胞分辨率下谱系特异性宿主响应与IFNβ刺激转录网络的收敛程度仍不清楚。因此研究人员基于已发表的单细胞转录组（single-cell RNA sequencing, scRNA-seq）数据，开展谱系分辨的二次描述性再分析，系统描述ZIKV暴露与IFNβ刺激在人发育神经组织中的转录模式及毒株间探索性差异。
本研究系对公开数据的二次分析，原文未标注新期刊投稿信息，用户提供发表期刊为《Genomics》。
主要技术方法
研究人员使用来自GEO数据库的三个公开数据集：主要分析对象为人胎儿脑来源神经组织scRNA-seq数据集GSE238140（Stokes et al., Cell Rep Med 2025），包含mock、巴西ZIKV-BR株、柬埔寨亚洲谱系FSS13025株（ZIKV-FSS）及重组IFNβ四种ex vivo处理条件，每条件各一个生物学样本；背景对照采用人脑类器官bulk RNA-seq数据集GSE97919（Watanabe et al., Cell Rep 2017）及含ZIKV感染Vero细胞和人单核细胞来源树突状细胞（moDC）的scRNA-seq数据集GSE230571。分析流程包括：质控过滤（线粒体基因比例阈值分别<12%和<10%）、SCTransform归一化、PCA及UMAP降维聚类（分辨率0.3–0.4）、依据经典标记基因（SOX2/PAX6/NES定义放射状胶质样细胞；MKI67/TOP2A/HMGB2定义cyclin NPCs；SLC17A7/NEUROD6为兴奋性神经元；GAD1/GAD2/DLX1为抑制性神经元；GFAP/AQP4为星形胶质细胞；C1QA/C1QB/TYROBP为小胶质细胞；PDGFRA/OLIG1/OLIG2为少突胶质祖细胞OPCs）进行细胞类型注释、按样本×细胞类型构建假 bulk（pseudobulk）矩阵、细胞水平Wilcoxon秩和检验进行探索性差异表达（|log2FC|≥0.25，FDR<0.05，表达于任一组≥10%细胞，因无生物学重复仅作描述性排序）、经典ISG特征（IFITM1、MX1、OAS1/2/3、IFIT1、IFI6、ISG15、STAT1、IRF7、RSAD2）log2FC热图展示及模块评分（AddModuleScore）、ZIKV–IFNβ转录相似性Pearson/Spearman相关性分析（多基因全集敏感性检验：DE表基因、仅经典ISGs、非ISG基因、剔除ISG的全表达基因）、独立bulk RNA-seq用DESeq2进行差异表达及VST（variance-stabilizing transformation）可视化。所有分析用R 4.4.1及Seurat v5.3.1完成。
研究结果
3.1 Global transcriptional structure and dataset integrity（全局转录结构与数据集完整性）
经质控后保留27,267个人胎儿神经细胞及63,889个Vero/moDC细胞。UMAP嵌入显示胎儿脑样本各处理条件及无监督聚类均具清晰转录分区；Vero细胞与moDC分处不同转录空间。肘图支持选前30个主成分。结果表明数据集结构清晰，适合后续谱系水平分析。
3.2 Cellular architecture of fetal neural tissue under ZIKV exposure（ZIKV暴露下胎儿神经组织的细胞架构）
无监督聚类鉴定出7大类群：祖细胞富集放射状胶质样细胞（10,501个）、cyclin NPCs（5459个）、兴奋性神经元（5670个）、抑制性神经元（3668个）、星形胶质细胞（1407个）、小胶质细胞（506个）及少量OPCs（56个）。各处理条件下细胞群转录可辨，条件间细胞比例有差异。按样本×细胞类型聚合为假bulk矩阵（23,183基因×28组合）供后续分析。
3.3 Lineage-specific transcriptional remodeling in response to viral and cytokine stimulation（病毒与细胞因子刺激下谱系特异性转录重塑）
探索性差异表达显示放射状胶质样细胞、cyclin NPCs、星形胶质细胞及小胶质细胞较神经元群体有更多基因达差异表达标准；ISGs在各谱系一致上调。在祖细胞聚焦及ISG相关汇总中，ZIKV-FSS转录模式与IFNβ更相似，ZIKV-BR则引发更广泛多变的变化。表明发育身份与病毒暴露条件关联于不同响应幅度及结构。
3.4 Exploratory strain-associated transcriptional patterns in progenitor populations（祖细胞群体中毒株关联探索性转录模式）
放射状胶质样及cyclin NPC火山图与经典ISG热图显示ZIKV-FSS与IFNβ的ISG相关谱视觉上重叠更强，ZIKV-BR响应更广。因每条件单样本，差异解读为毒株关联探索模式而非确定毒株依赖效应。ISG诱导反映转录IFN/ISG响应性，不能推定细胞自主抗病毒启动或感染率。
3.5 Comparative alignment between ZIKV-associated and IFNβ-associated transcriptional programs（ZIKV关联与IFNβ关联转录程序的比较对齐）
以DE表基因计算ZIKV–IFNβ log2FC相关性较高；但敏感性分析显示相关性受基因全集影响大——经典ISGs间高度一致，非ISG或剔除ISG的全表达基因集Spearman相关明显降低。提示ZIKV与IFNβ共享IFN/ISG响应但有更宽泛保守估计，不等同功能等价。
3.6 Transcriptional response magnitude across developmental identities（跨发育身份的转录响应幅度）
定量显示放射状胶质样细胞差异表达基因数最多，cyclin NPC有明显ISG响应及部分细胞周期基因变化。ZIKV暴露可在无外源细胞因子时诱导祖细胞抗病毒转录本。说明发育身份关联于观测抗病毒转录响应规模差异。
3.7 Independent bulk RNA-seq contextual analysis（独立bulk RNA-seq背景分析）
人脑类器官ZIKV vs mock的bulk RNA-seq中，经典ISGs协调上调，PCA部分分离。支持ZIKV暴露脑类器官中存在IFN/ISG关联转录激活，但因无单细胞分辨率不验证谱系特异性模式。
3.8 Contextual comparison across immune and epithelial cell models（跨免疫与上pi上皮细胞模型背景比较）
GSE230571中moDC显示强ISG模块评分，Vero细胞极弱。符合两系统先天免疫能力差异，为ZIKV暴露伴IFN/ISG程序提供跨体系背景支持，但不验证胎儿神经谱系特异性响应。
讨论与结论翻译
讨论指出本研究是对Stokes et al.数据的谱系聚焦二次再分析，提供标准化细胞类型水平描述性总结、明确无重复的限制并检验ZIKV–IFNβ相似性对基因集选择的敏感性。放射状胶质样及cyclin NPC显示强IFN/ISG转录响应性，ISG诱导可能源于旁分泌信号而非细胞自主启动。ZIKV-BR与ZIKV-FSS在差异表达基因数目与构成上有探索性差别，但因单样本仅为毒株关联模式。ZIKV与IFNβ转录响应在经典ISGs高度一致，相关性受基因全集影响，二者不等同功能等价。跨数据集比较仅作背景不支持验证胎儿神经谱系结论。局限性含ex vivo条件不完全模拟体内、无生物学重复、放射状胶质样群未细分亚型、跨体系不直接可比。
结论（Concluding remarks原文凝练翻译）：
研究人员通过对ZIKV及IFNβ暴露后发育中人神经组织公开scRNA-seq数据的探索性谱系分辨二次分析发现：在各注释神经谱系中，祖细胞富集放射状胶质样细胞及cyclin NPCs表现出强IFN/ISG关联转录响应性及更多达探索性差异表达标准的基因；ZIKV-BR与ZIKV-FSS间存在探索性毒株关联转录响应模式差异，ZIKV关联与IFNβ关联基因表达谱有部分重叠——经典ISGs各条件一致诱导，但病毒暴露附加的基因变化尤在广谱非ISG或全表达基因集中不能被IFNβ单独复现。结果提供了发育中人神经组织抗病毒转录状态的描述性谱系分辨参考，并显示发育身份与病毒暴露条件如何关联于宿主转录响应变异。因主要胎儿神经数据集每条件仅一个生物学样本，结果应被解读为探索性毒株关联模式而非毒株依赖机制的确证证据。
