《Veterinary Research》：Natural polymorphisms in the bovine leukemia virus microRNA cluster modulate miRNA expression and host regulatory pathways
编辑推荐：
研究人员对自然感染的53头牛的外周血白细胞进行分析，扩增并测序了长度为554核苷酸（nt）的牛白血病病毒（Bovine Leukemia Virus, BLV）miRNA编码区。选择含有RNA聚合酶III（Pol III）启动子基序、种子区及终止信号突变的田间
研究人员对自然感染的53头牛的外周血白细胞进行分析，扩增并测序了长度为554核苷酸（nt）的牛白血病病毒（Bovine Leukemia Virus, BLV）miRNA编码区。选择含有RNA聚合酶III（Pol III）启动子基序、种子区及终止信号突变的田间分离株进行功能评估。参考序列与各变异株在受控条件下与缺失miRNA的BLV感染性克隆（BLVΔ-miRNA）共表达于HEK293T细胞。成熟miRNA水平通过茎环逆转录定量PCR（stem-loop RT–qPCR）定量，宿主转录组反应采用寡核苷酸微阵列分析，并利用生物信息学预测miRNA靶基因。序列分析鉴定出84个多态位点，其中大量位于Pol III调控元件和种子区。变异株显示特定成熟miRNA积累改变及预测的依赖种子的靶基因谱变化。转录组分析表明，变异株与先天免疫信号组分、干扰素应答基因、抗原呈递通路、肿瘤抑制网络及细胞外基质相关过程的调控相关。富集分析显示预测的miRNA靶基因与下调转录本存在显著重叠。结果表明，BLV miRNA簇的自然多态性可调节miRNA表达，并在受控实验系统中与不同的宿主调控特征相关联，提示病毒miRNA区的序列变异可能参与差异性的宿主–病毒相互作用并影响支持病毒持久性的机制。
本研究由研究人员发表于《Veterinary Research》，围绕牛白血病病毒（BLV）miRNA编码区的自然单核苷酸多态性（SNP）展开。BLV属于δ逆转录病毒属，可在B淋巴细胞中建立终身潜伏感染，其miRNA簇在潜伏期仍高水平表达，但其自然序列变异对miRNA生成及宿主基因调控的影响尚不清楚。此前的研究已证实BLV编码的miRNA可通过调控宿主信号网络促进病毒持久性和致癌潜能，但在田间病毒群体中的序列多样性及其功能后果缺乏系统性解析。研究人员旨在揭示这些自然多态性是否通过改变miRNA的生物发生、稳定性或靶基因特异性，进而重塑宿主调控网络。
在技术方法上，研究纳入波兰地方性牛白血病的监测项目中的53份自然感染BLV的牛全血样本，分离外周血白细胞并提取基因组DNA。通过巢式PCR扩增并测序BLV miRNA编码区，筛选含Pol III启动子、种子区及终止信号突变的田间分离株。将参考序列及变异株克隆至启动子缺失载体，并与BLVΔ-miRNA共转染HEK293T细胞，采用茎环RT–qPCR检测成熟miRNA水平，利用寡核苷酸微阵列分析宿主转录组变化，结合生物信息学预测miRNA靶基因并进行功能富集分析。统计学分析采用非参数Kruskal–Wallis检验及超几何分布检验。
SNP identification in the miRNA-coding locus
对53条554 nt长的BLV miRNA序列与参考株JC613347比对，发现84个多态位点，包括80个SNP和4个插入/缺失。约72.6%的多态位于调控区，其中隐秘A盒多态最多（28个），其次为B盒（19个）、A盒（10个）、种子区（9个）、隐秘B盒（3个）及终止信号（4个）。
Detailed analysis of SNPs within miRNA functional motifs of selected BLV isolates
选取7个代表性分离株，携带跨越种子区、Pol III启动子元件及终止信号的33个突变。V29分离株的变异最为广泛，包括miR-B1-5p、miR-B5种子区替换，多个Pol III启动子基序变异及插入/缺失，提示其对miRNA靶特异性及转录效率可能产生显著影响。其余分离株则集中在热点区域，如B1/A-盒70位、B3/A簇211位及终止区427位的突变。
Expression analysis of BLV miRNA variants
在HEK293T细胞中，与BLVΔ-miRNA共表达的变异株显示，V29显著降低miR-B1-5p、miR-B2-3p及miR-B5-3p水平，V90降低miR-B2-3p水平，其余变异株与参考序列无显著差异。各组的病毒转录本（tax/rex、gag、env）表达无统计学差异，提示转染效率一致。
Microarray-based analysis of DEGs across experimental groups
转录组分析显示，V29、V90、参考组及空载体对照组分别有1747、1954、2062、2034个差异表达基因（DEGs），均以蛋白编码基因和长链非编码RNA（lncRNA）为主，且下调基因数量多于上调基因。所有组共享885个核心DEGs，但各组亦有独特基因子集，表明BLV miRNA多态性并未引发全新的转录表型，而是调节已有转录反应的幅度。
Distinct pathway enrichment signatures across experimental groups
通路富集分析显示，所有组均富集通用转录通路、整合素细胞表面相互作用、特发性肺纤维化信号、病原体诱导的细胞因子风暴信号、rRNA加工及破骨细胞在类风湿性关节炎中的作用通路。V29偏向细胞周期调控及RNA代谢，V90则在细胞外基质重塑、免疫及癌症信号、细胞周期及应激反应等方面富集更广泛。
Differential regulation of immune signaling pathways
与参考组相比，V29和V90在病原体识别受体通路、干扰素α/β信号、抗原呈递及巨噬细胞经典激活信号等免疫相关通路中表现出基因表达差异，提示miRNA变异可能调节抗病毒及炎症信号，从而影响病毒持久性与宿主防御。
Altered expression of cancer-associated and tissue remodeling pathways
V29和V90在癌症相关通路（如组织因子在癌症中的作用、DNA损伤修复蛋白SUMO化、细胞周期检查点、分子致癌机制、RHO GTP酶循环及VEGF信号）及细胞外基质重塑、胶原合成与降解、伤口愈合信号等过程中表现出差异表达，提示这些变异可能影响肿瘤相关微环境。
RT–qPCR validation of microarray data
RT–qPCR验证了微阵列结果，多数基因的表达趋势一致，如EGR1和TNFAIP6上调，ATP7A、RPL23AP32及KLF12下调。
Predicted structural and accessibility differences of miRNA variants
热力学预测显示，V29和V90的miRNA序列变异改变了成熟miRNA的二级结构和可及性，尤其是miR-B5-5p（V29）和miR-B2-3p（V90）的种子区结构稳定性下降，提示靶位点结合效率可能受影响。
Impact of seed-region polymorphisms in BLV miRNAs on host gene regulation
种子区多态导致预测靶基因谱的改变，约半数下调基因被预测为BLV miRNA靶基因。功能富集显示参考组偏向发育程序，V29偏向细胞周期与RNA代谢，V90偏向免疫及细胞死亡通路。多个预测的靶基因属于已知肿瘤抑制基因，包括PTEN、TSC2、NF1、TGFBR2等。
Bovine orthologs of human genes
在牛基因组的同源基因中，约45%–60%保留了预测的BLV miRNA结合位点，分布于3′UTR和编码区，提示在人类细胞中的调控模式可能在自然宿主中部分保守。
在讨论部分，研究人员指出BLV miRNA区的自然多态性并非中性变异，而是功能性调控多样性的来源。Pol III启动子的冗余结构可能是进化上确保miRNA稳定转录的策略，而种子区多态可重塑靶基因网络，影响免疫信号、细胞周期及肿瘤抑制通路。这些变化可能通过削弱先天免疫识别和改变细胞外基质微环境，促进病毒感染细胞的存活与克隆扩增。尽管研究在HEK293T细胞中进行，但通过牛基因同源分析及保守通路预测，为后续在天然宿主中的验证提供了基础。
结论
BLV miRNA簇的自然多态性可调节miRNA表达及宿主基因调控网络，尤其在免疫调节、细胞应激及肿瘤抑制相关通路中发挥潜在作用。miRNA序列变异是病毒遗传景观的一部分，可能与病毒持久性及致病性差异相关。进一步在自然宿主及体内模型中的研究将有助于明确这些多态性在病毒流行病学与临床进展中的意义。
