牛白血病病毒（BLV）感染临床分期的先天免疫转录组重塑：无症状感染中的受控抗病毒激活与持续性淋巴细胞增多症（PL）中的炎症重编程

《Virus Genes》：Innate immune transcriptomic remodeling across clinical stages of bovine leukemia virus infection

牛白血病病毒（bovine leukemia virus, BLV）是一种致癌性δ逆转录病毒（deltaretrovirus），可在牛中建立终身感染，并可能从无症状携带进展为持续性淋巴细胞增多症（persistent lymphocytosis, PL）和淋巴瘤。尽管适应性免疫改变已被广泛研究，但BLV感染中先天免疫的阶段特异性转录景观仍未完全明确。在本研究中，研究人员利用RNA-seq数据开展靶向先天免疫转录组分析，数据来源于墨西哥克雷塔罗（Querétaro）筛查的130头动物群体中选出的18头荷斯坦奶牛。根据血清学、前病毒载量和血液学标准将动物分为三组：BLV血清阴性对照（n=5）、无症状BLV感染牛（asymptomatic BLV-infected cattle, AL）和伴PL的BLV感染牛。对外周血白细胞进行RNA提取，并在Illumina NovaSeq平台测序。使用DESeq2进行差异基因表达分析，显著性定义为校正P值（adjusted P-value, padj）<0.05且|log₂ fold change|≥1。分析了一个精选的先天免疫相关基因panel，涵盖模式识别受体（pattern recognition receptors, PRRs）、干扰素信号、炎症介质、炎症小体（inflammasome）组分和细胞毒性效应分子。对照与无症状牛比较鉴定出8个显著改变基因（padj<0.05）。无症状动物显示抗病毒和免疫调节基因上调，包括EIF2AK2（log₂FC=1.10, padj=0.011）、KLRK1（log₂FC=0.98, padj=0.009）和IL18（log₂FC=0.77, padj=0.015），表明干扰素相关和炎症通路激活。TNFRSF13C（log₂FC=3.25, padj=0.046）和TNFAIP3（log₂FC=2.03, padj=0.049）显著上调，提示B细胞存活信号和NF-κB调控调节。相反，NCR1（log₂FC=?1.04, padj=0.044）、LMX1B（log₂FC=?0.80, padj=0.030）和TNFSF12（log₂FC=?0.78, padj=0.044）下调，反映细胞毒性和TNF超家族通路的选择性调整。PL动物表现出更广泛的先天免疫转录改变，与持续炎症激活一致。这些发现表明，无症状BLV感染以抗病毒和调节性先天免疫通路的受控激活为特征，而向PL疾病进展与扩大的炎症重编程相关。阶段特异性先天免疫特征可能因此促进BLV感染牛的病毒持续和临床结局。

研究背景：BLV（bovine leukemia virus）是致癌性δ逆转录病毒，引起地方性牛白血病，主要感染CD5⁺ IgM⁺ B淋巴细胞并整合前病毒基因组，建立终身持续感染。多数感染牛无症状，约三分之一发生PL，1–5%经长潜伏期进展为恶性B细胞淋巴瘤；感染还导致免疫功能改变、继发感染易感、生产力下降和经济损失。BLV与人类T淋巴细胞病毒1型（HTLV-1）相关，是研究逆转录病毒持续和白血病发生的大动物模型。病毒持续主要依赖感染B细胞克隆扩增和受控病毒基因表达。先天免疫是抗病毒第一道防线，通过PRRs如Toll样受体（TLRs）、RIG-I样受体（RLRs）和胞质核酸感受器识别病毒，诱导I型干扰素（type I interferons, IFN-I）、干扰素刺激基因（interferon-stimulated genes, ISGs）和炎症介质。但既往BLV研究多关注适应性免疫，先天免疫转录程序在不同临床阶段仍定义不足，且常依赖靶向基因或细胞因子，难以捕捉网络协调。RNA-seq可无偏全局分析。研究人员因此对已发表RNA-seq数据集进行假设驱动的二次分析，聚焦先天免疫相关基因在无症状、PL等阶段的协调重塑；目的包括刻画BLV感染先天免疫通路转录改变、比较不同临床阶段特征、识别可能参与病毒持续或疾病进展的干扰素调节基因、炎症介质和抗病毒效应分子；意义在于深化宿主–病毒互作理解，为疾病控制、遗传抗性策略和生物标志物提供分子见解。

关键技术方法（≤250字）：本研究为已发表全转录组RNA-seq数据集的聚焦二次分析。样本来源：墨西哥Querétaro一个商业牛群130头荷斯坦奶牛筛查群体中选18头临床健康泌乳中期（第三或第四泌乳期）牛；按血清学、PCR和qPCR前病毒载量（proviral load, PVL）及血液学分为BLV血清阴性对照（n=5）、无症状BLV感染（AL）和BLV感染伴PL。取总外周血白细胞（peripheral blood leukocytes, PBLs）RNA，Illumina NovaSeq 6000平台150-bp双端测序，平均约29 million read pairs。FastQC质控，HISAT2比对Bos taurus ARS-UCD1.2，featureCounts计数，DESeq2差异表达，阈值padj<0.05且|log₂FC|≥1。预先构建101个先天免疫相关基因panel，分PRRs、IFN-I信号/ISGs、促炎细胞因子/NF-κB、炎症小体/程序性细胞死亡、NK细胞细胞毒性五类，并进行差异表达、聚类和通路解释。

结果：
Differential innate immune gene expression across clinical stages of BLV infection（BLV感染临床分期中的差异先天免疫基因表达）：使用预设101基因panel，进行Control vs AL、Control vs PL、PL vs AL比较，以错误发现率（false discovery rate, FDR）控制，padj<0.05为显著，|log₂FC|≥1评估生物相关性。
Control vs. asymptomatic：Control与无症状比较鉴定8个padj<0.05基因。无症状中EIF2AK2、KLRK1、IL18、TNFRSF13C、TNFAIP3上调，NCR1、LMX1B、TNFSF12下调。结论：无症状BLV感染呈混合的先天免疫激活和调节调整，干扰素相关信号、炎症介质、免疫调节和B细胞存活信号升高，部分细胞毒性受体降低，支持受控免疫激活而非广泛炎症失调。
Differential gene expression between control and persistent lymphocytosis：Control与PL组中101基因panel无基因达到padj<0.05且|log₂FC|≥1；最小padj=0.0678，最大上调log₂FC=3.12，最强下调log₂FC=?1.12，但均不显著。结论：在该panel和统计阈值下未检测到稳健差异，可能为细微转录改变或统计功效不足。
Persistent lymphocytosis vs. asymptomatic：PL与无症状比较鉴定5个显著失调基因：TNFAIP3、TNFRSF13C、IL6、EIF2AK2上调，C1QTNF9下调；NCR1降低但padj=0.052为趋势。结论：PL中炎症和免疫调节信号增强，C1QTNF9降低。
Hierarchical clustering of differentially expressed genes：基于Z-score的聚类热图显示CT、AS、PL按临床组分离，阶段依赖免疫重塑。AS中STAT1、STAT2、IRF7、IFNAR1、IFNAR2、MX1、OAS家族等I型干扰素信号基因相对高表达，形成抗病毒簇；PL中许多ISG相对降低，CASP1、NLRP3、GSDMD等炎症小体和焦亡信号较AS低，但TLR3、TLR7、TLR9、IL1B、IL18、CCL2、NOD2等炎症和PRR基因升高且与干扰素模块分开。结论：从AS广谱抗病毒谱转向PL重构炎症景观，干扰素相关表达降低并选择性富集TLR和细胞因子转录物。

讨论总结：研究人员指出这是假设驱动二次分析，不是原基因组范围研究的重复；新颖性在先天免疫系统免疫学解读。先天免疫分析揭示阶段依赖抗病毒和炎症程序重塑：无症状感染相关协调抗病毒特征，PL相关选择性炎症激活和干扰素网络失调。无症状中受控但协调的干扰素相关基因激活，EIF2AK2等上调，无广泛炎症升级；KLRK1（NKG2D）提示NK细胞监视，IL18促炎并桥接先天适应，TNFRSF13C（B细胞活化因子受体，BAFFR）和TNFAIP3（A20）提示B细胞存活信号和NF-κB负反馈，NCR1和TNFSF12下调提示效应机制微调。该平衡允许病毒持续而少免疫病理。Control vs PL无显著差异需谨慎解释，可能因细微改变、PL异质性、PVL和克隆扩增、总白细胞组成变化掩盖细胞类型特异信号；Polat et al.、Shrestha et al.、Trueblood et al.等支持背景依赖和克隆增殖。PL vs AL中TNFAIP3、TNFRSF13C、IL6上调，提示NF-κB反馈、BAFF介导B细胞存活扩增和促炎微环境；IL6与慢性炎症和肿瘤过程相关，可能促进淋巴细胞增殖。聚类显示AS中经典ISG协调升高，PL中减弱并富集TLR、IL1B、IL18和趋化因子，提示从抗病毒协调转向炎症重编程；Ivashkiv and Donlin、Deng et al.、Iwan et al.、Florins et al.等描述慢性感染中干扰素信号功能重编程和阶段特异细胞因子改变。总体模型：无症状为受控先天抗病毒信号，PL为炎症激活、调节反馈和干扰素网络选择性减弱，可能降低抗病毒压力、维持促B细胞扩增的慢性炎症，利于病毒持续。与HTLV-1持续携带的免疫调节相似。限制：二次分析、样本量继承、无先验功效、总PBL非分选、无qPCR/流式/蛋白/功能验证，候选标志需独立大队列和功能研究验证。未来方向包括单细胞转录组、PVL与转录组整合、纵向采样、TNFAIP3/TNFRSF13C/IL6功能验证、宿主遗传和表观调控、与HTLV-1比较。

结论部分翻译：本研究提供了基于RNA-seq的聚焦表征，揭示BLV感染不同临床阶段先天免疫转录动态。通过检测一个精选的涵盖抗病毒感知、干扰素信号、炎症、炎症小体和细胞毒性相关基因的panel，研究人员证明BLV感染与阶段依赖性先天免疫重塑相关，而非统一的免疫激活。无症状牛表现出受限但协调的抗病毒转录谱，以干扰素相关基因中度诱导和免疫平衡维持为标志。相反，PL以选择性失调为特征，包括TNFAIP3、TNFRSF13C和IL6等炎症和免疫调节介质表达增强，同时经典干扰素刺激基因聚类模式减弱或重新分布。值得注意的是，在严格统计阈值下，对照与PL组之间缺乏稳健差异表达，提示BLV感染中的免疫改变可能细微、背景依赖，并受白细胞组成和PVL动态影响。总体而言，这些发现支持一个模型：早期BLV感染与受控抗病毒参与相关，而向PL进展涉及炎症重编程和调节反馈机制，可能有利于病毒持续和B细胞扩增。该工作促进了对BLV发病机制中先天免疫应答更整合的理解，并强调慢性逆转录病毒感染中宿主–病毒平衡的复杂性。