引言
猫传染性腹膜炎(FIP)是猫最具毁灭性的病毒性疾病之一,由猫传染性腹膜炎病毒(FIPV)引起,FIPV是猫冠状病毒的一种高毒力生物型[1]。尽管经过数十年的研究,FIP仍然几乎普遍致命,有效治疗手段有限,对猫的健康和兽医诊疗构成了重大负担。临床上,FIP表现为脓肉芽肿性炎症、渗出性腹膜炎和多器官衰竭,反映了广泛的内皮损伤和免疫失调[2,3]。鉴于其高致死率和全身性特征,FIPV不仅是一个紧迫的临床靶标,也是剖析冠状病毒-宿主相互作用的理想模型,尤其是那些决定病毒RNA积累和宿主限制逃逸的相互作用。
冠状病毒已进化出复杂的分子策略来操控宿主细胞并确保高效复制。这场病毒-宿主军备竞赛的一个关键前沿在于细胞内在RNA限制系统,该系统能够识别并降解外源RNA转录本以阻止感染[4,5]。然而,病毒对抗这一防御的全谱反击机制仍未完全阐明。因此,理解致病性冠状病毒(如FIPV)如何颠覆这些RNA监视通路,对于解释其如何维持病毒RNA积累、高效复制及感染相关致病机制具有重要意义。
在冠状病毒结构蛋白中,核衣壳(N)蛋白已成为一种多功能调控因子,通过与病毒RNA及复制-转录复合体(RTC)组分结合来促进病毒RNA合成[6]。通过其高亲和力RNA结合活性以及与nsp3和nsp12等非结构蛋白的相互作用,N蛋白稳定病毒RNA模板、促进RTC组装并增强基因组复制和亚基因组转录的效率[7-9]。除这一经典功能外,N蛋白还越来越多地被发现参与调控宿主RNA代谢[10,11]和抗病毒信号传导[12,13],提示其可能还影响病毒RNA对宿主RNA监视通路的敏感性。然而,这一推测的RNA保护功能是否涉及对抗宿主细胞质RNA水平限制——一种作为病毒复制早期细胞内屏障的关键抗病毒通路[14]——目前仍知之甚少。阐明这一联系将为理解冠状病毒N蛋白如何将病毒RNA合成与逃逸宿主RNA限制系统相整合提供关键见解。
为解决这一问题,我们在FIPV感染细胞中进行了非偏倚的免疫共沉淀联合质谱(Co-IP/MS)筛选,以鉴定与N蛋白结合的宿主蛋白。分析揭示了与宿主RNA结合因子的广泛相互作用,其中限制因子ZC3HAV1(也称为ZAP或PARP13)作为显著候选因子出现。ZC3HAV1是一种干扰素刺激基因产物,可结合病毒RNA,并在多种病毒系统中通过招募RNA监视和降解相关因子来抑制复制[15,16]。越来越多的证据已确立ZC3HAV1作为针对多种病毒(包括多种正链RNA病毒和冠状病毒)的抗病毒限制因子。在冠状病毒中,已有报道SARS-CoV-2和猪流行性腹泻病毒(PEDV)受ZC3HAV1介导的限制[17,18],表明尽管许多冠状病毒基因组普遍具有CpG抑制特征,冠状病毒RNA仍可对ZC3HAV1相关的RNA水平限制敏感。然而,低CpG冠状病毒RNA仍对ZC3HAV1相关限制敏感的机制基础仍未完全阐明,特别是残余CpG分布、其他RNA序列背景特征或病毒特异性RNA结合蛋白是否影响这一过程。此外,既往冠状病毒研究主要聚焦于ZC3HAV1本身的抗病毒活性,而冠状病毒如何对抗ZC3HAV1介导的RNA限制则了解甚少。特别是,冠状病毒结构蛋白是否能干扰ZC3HAV1-病毒RNA轴尚未完全明确。值得注意的是,FIPV具有低CpG基因组,而ZC3HAV1经典定义为CpG偏向性RNA结合抗病毒因子[17]。这一明显悖论使FIPV成为研究CpG耗竭的冠状病毒RNA如何仍对ZC3HAV1相关RNA水平限制敏感,以及冠状病毒N蛋白是否通过拮抗ZC3HAV1-病毒RNA轴来维持病毒RNA积累和持续存在的有用系统。
我们的初步数据表明,ZC3HAV1与FIPV N蛋白相互作用,并作为病毒复制的负调控因子发挥作用。我们假设FIPV N通过降低ZC3HAV1与病毒RNA的结合来拮抗ZC3HAV1相关的RNA水平限制,从而维持病毒RNA的积累和持续存在。通过聚焦于一种临床致死性猫冠状病毒,本研究不仅揭示了一种RNA水平的反向限制机制,还为冠状病毒如何维持转录本稳定性以维持感染提供了更广泛的见解。我们确定了ZC3HAV1含RBD区域是参与FIPV N结合的关键界面,证明FIPV N损害ZC3HAV1与病毒RNA的结合,并表明ZC3HAV1对FIPV RNA的识别不能完全由整体CpG丰度解释。总之,这些发现通过揭示冠状病毒结构蛋白如何颠覆内在RNA限制,桥接了兽医病毒学和分子病毒学,从而拓展了我们对冠状病毒-宿主在转录后水平相互作用的理解。
材料与方法
病毒、载体、细胞和动物
SeV、FIPV毒株SH2021(GenBank登录号OR295209.1)和日本脑炎病毒(JEV)由上海兽医研究所(SHVRI)伴侣动物疾病团队提供。质粒pCold-SUMO-FIPV-N、pET-28a、pLKO.1、psPAX2、pMD2.G、pCAGGS、pCMV14-Flag、pCMV14-血凝素(HA)、pCMV14-Myc、pCMV14-绿色荧光蛋白(GFP)、pGAS-Luc-IFN-β(人)、pGL3-IFN-β(猫)和pRL-TK也由SHVRI提供。CRFK和HEK293T细胞购自中国科学院细胞库。BHK-21细胞由SHVRI伴侣动物疾病团队提供。
质粒构建
质粒通过同源重组使用One Step Cloning Kit(Vazyme Biotech Co., Ltd.,南京,中国)构建。N基因从pCold-SUMO-FIPV-N质粒中扩增,宿主基因从CRFK细胞RNA合成的cDNA中扩增。PCR使用2 × Phanta Max Master Mix(Vazyme)按照制造商说明进行。质粒通过双酶切线性化,插入片段通过同源重组组装。原核表达质粒转化至大肠杆菌Rosetta(DE3),真核质粒转化至大肠杆菌DH5α。引物序列列于S1表。
蛋白质印迹(WB)分析
细胞接种并以Lipofectamine 3000(Thermo Fisher Scientific,Waltham,MA,USA)按照制造商说明在铺板后16–24小时达到80–90%汇合度时转染质粒。12小时后,培养基更换为维持培养基(293T细胞除外)。在指示时间点收集细胞或病毒感染样本,于冰上用RIPA缓冲液裂解30分钟,随后在4°C下以12,000 rpm离心5分钟。上清液与4 × SDS上样缓冲液混合,煮沸10分钟,保存于−40°C。蛋白质经SDS-PAGE分离,转移至0.22 μm PVDF膜(Merck Millipore,Burlington,MA,USA),用5%脱脂奶粉封闭,并与一抗(按制造商推荐)在4°C过夜或25°C孵育1小时。与二抗孵育后,使用ECL底物检测信号,并通过成像系统(Cat#4600;Biotanon,上海,中国)可视化。
50%组织培养感染剂量(TCID50)的测定
为测定FIPV的TCID50,收集全细胞裂解物并进行三次冻融循环,随后在4°C下以3,000 × g离心10分钟。收集上清液并在Eagle最低必需培养基中进行10倍系列稀释(10-1至10-8)。在96孔板中接种的CRFK细胞洗涤以去除培养基。每孔加入100微升各病毒稀释液,每个稀释度设八个重复。包括未接种病毒的阴性对照孔。平板在37°C、5% CO2气氛中孵育5天。记录细胞病变效应(CPE),并使用Reed–Muench法计算TCID50值。
RT-qPCR
为研究宿主对FIPV感染的转录反应,CRFK细胞感染FIPV(感染复数[MOI] = 0.01)或模拟处理。在感染后12和24小时(hpi)使用TRIzol试剂(Accurate Biotechnology Co., Ltd.,上海,中国)提取总RNA,用Evo M-MLV(Accurate Biotechnology Co., Ltd.)逆转录为cDNA,并使用SYBR Green Pro Taq HS Mix(Accurate Biotechnology Co., Ltd.)和基因特异性引物(S2表)进行qPCR。相对mRNA水平以甘油醛3-磷酸脱氢酶(GAPDH)归一化,并使用2−ΔΔCt方法计算。
对于绝对定量,FIPV RNA使用靶向N基因或非翻译区(UTR)的引物组进行定量。基于N基因的定量,引物为正向5′-TAAGGAACTCGCTGAGAGGTGG-3′和反向5′-CTGAAACTGTGGCGGTATCTTA-3′。生成范围为10²至10⁸拷贝/μL的标准曲线(y = −3.3351x + 37.367,R² = 0.9958)。基于UTR的定量,引物为正向5′-TCCGAGGAATTACTGGTCATC-3′和反向5′-CACTAGATCCAGACGTTAGCTC-3′。生成范围为10²至10⁸拷贝/μL的标准曲线(y = −3.3141x + 36.214,R² = 0.9947)。病毒RNA载量根据相应标准曲线计算,并以每微克总RNA的拷贝数表示。
免疫沉淀(IP)与免疫共沉淀(Co-IP)
IP采用抗N单克隆抗体或基于标签的磁珠进行。对于直接IP,将感染FIPV(MOI = 0.01)的CRFK细胞于感染后16小时裂解,细胞裂解液与抗FIPV N单克隆抗体(克隆10E7)孵育,随后按照制造商说明加入Protein A/G琼脂糖珠(货号36403ES25;Yeasen Biotech Co., Ltd.,上海,中国)。洗脱的蛋白质通过WB和MS进行分析。
对于Co-IP,使用Lipofectamine 3000将标记蛋白构建体转染至CRFK或293T细胞。在指示情况下,转染后CRFK细胞感染FIPV。转染后24小时,裂解细胞,按照制造商说明使用抗标签磁珠(ShareBio,上海,中国)免疫沉淀标记蛋白。洗涤后,洗脱结合的蛋白质进行SDS-PAGE,随后进行WB或MS。
液相色谱(LC)-MS
切取目标凝胶条带并进行胶内胰蛋白酶消化。肽段由上海APT Biotechnology Co., Ltd.(上海,中国)通过LC-串联MS(LC-MS/MS)分析进行蛋白质鉴定。
宿主蛋白的筛选
为鉴定与FIPV N蛋白相互作用的宿主蛋白,分析了N蛋白复合物的IP MS数据。选择高丰度宿主蛋白进行进一步研究。构建宿主蛋白的过表达质粒并转染至CRFK细胞。将CRFK细胞培养至80%–90%汇合度后,使用Lipofectamine 3000转染表达不同宿主蛋白的质粒或空pCAGGS对照质粒。转染12小时后,将培养基更换为维持培养基。24小时后,洗涤细胞并以MOI为0.01感染FIPV。感染16小时后,细胞经历冻融循环,离心收集上清液。使用TCID50方法测定上清液中的FIPV滴度。
为分析FIPV感染对宿主基因表达的影响,进行qPCR以测量感染过程中宿主基因的转录变化。将CRFK细胞接种于6孔板中,以MOI为0.01感染FIPV,未感染细胞作为对照。在感染后12和24小时(hpi)收集样品。使用TRIZOL提取总RNA,并使用Evo M-MLV逆转录试剂盒合成cDNA。随后使用SYBR Green Pro Taq HS qPCR试剂盒进行qPCR,以评估可能参与FIPV复制的宿主基因的转录水平。
ZC3HAV1过表达对FIPV复制的影响
为研究ZC3HAV1在FIPV复制中的作用,将pCAGGS-ZC3HAV1或对照质粒pCAGGS-Flag转染至CRFK细胞。转染后12小时,更换培养基,随后细胞感染FIPV(MOI = 0.01)。在感染后指示时间点收集细胞,通过Western blotting分析FIPV N蛋白,通过绝对RT-qPCR分析病毒RNA水平,并通过TCID50测定分析病毒滴度。此外,通过转染递增量的ZC3HAV1表达质粒评估ZC3HAV1的剂量依赖性效应,随后在感染后16小时分析FIPV N蛋白。
为评估FIPV感染是否影响内源性ZC3HAV1表达,将CRFK细胞感染FIPV(MOI = 0.01)并在指示时间点收集,通过RT-qPCR和Western blotting进行分析。此外,为评估FIPV N蛋白是否直接调控ZC3HAV1表达,将Flag-ZC3HAV1与递增量的HA标记N质粒共转染至CRFK细胞,并通过Western blotting测定ZC3HAV1蛋白水平。
敲低内源性ZC3HAV1对FIPV复制的影响
设计并合成了三条靶向ZC3HAV1的小干扰RNA(siRNA)序列。其序列列于S3表。将每种siRNA或非靶向对照(siNC)转染至CRFK细胞,在转染后48小时通过qPCR评估敲低效率,并在60小时通过WB评估,以鉴定最有效的siRNA。为评估其对FIPV复制的影响,将最佳siRNA或siNC转染至CRFK细胞,48小时后感染FIPV(MOI = 0.01)。在感染后12、16和24小时收集样品,用于WB分析N蛋白和内源性ZC3HAV1表达。病毒RNA水平和滴度分别通过qPCR和TCID50测定定量。
双荧光素酶报告基因实验
ZC3HAV1和FIPV N对IFN-β转录的影响:将293T细胞接种于24孔板中,共转染pGAS-Luc-IFN-β(60 ng)、pRL-TK(10 ng)以及pCAGGS-Flag(400 ng)、pCAGGS-FIPV N(400 ng)或pCAGGS-ZC3HAV1(400 ng)。18小时后,用SeV刺激细胞12小时,按照制造商说明使用双荧光素酶报告基因检测试剂盒(Easen Biotechnology Co., Ltd.,上海,中国)测定荧光素酶活性。
ZC3HAV1对FIPV N介导的IFN-β调控的拮抗效应:将293T细胞接种于24孔板,待融合度达80–90%时转染pGAS-Luc-IFN-β(人源,60 ng)、pRL-TK(10 ng)以及以下质粒组合之一:pCAGGS-Flag(400 ng);pCAGGS-FIPV N(200 ng、400 ng或800 ng);pCAGGS-ZC3HAV1(200 ng)+ pCAGGS-Flag(200 ng);或pCAGGS-FIPV N(200 ng)+ pCAGGS-ZC3HAV1(200 ng)。每种条件设三个复孔。转染16 h后,用SeV刺激细胞额外12 h。按照试剂盒说明书,使用双荧光素酶报告基因检测试剂盒测定荧光素酶活性。对于CRFK细胞,采用相同操作,但将pGAS-Luc-IFN-β(人源)替换为pGL3-IFN-β(猫源)以进行物种特异性验证。
ZC3HAV1与FIPV N相互作用的IP和Co-IP验证
通过IP和Co-IP实验验证了ZC3HAV1与FIPV N蛋白之间的相互作用。为排除RNA介导该相互作用的可能性,在Co-IP实验的裂解过程中向细胞裂解液中加入RNase(终浓度:50 µg/mL)。
共聚焦显微镜分析
为检测ZC3HAV1与FIPV N的亚细胞定位及共定位,将CRFK细胞接种于玻璃盖玻片上,并转染以下表达质粒组合:pCMV14-Myc-FIPV N + pCMV14-HA(对照)、pCMV14-HA-ZC3HAV1 + pCMV14-Myc(对照)或pCMV14-Myc-FIPV N + pCMV14-HA-ZC3HAV1(实验组)。转染后24 h,用4%多聚甲醛固定细胞,用2% Triton X-100透化,并使用抗血凝素(HA)和抗Myc抗体孵育,随后用AF594和AF488偶联的二抗孵育。细胞核用DAPI复染。封片后,使用激光扫描共聚焦显微镜观察细胞。
同时,为确认生理感染条件下内源性蛋白的亚细胞定位,将感染FIPV(MOI = 0.01)的CRFK细胞用4%多聚甲醛固定,并用2% Triton X-100透化。以模拟感染细胞作为阴性对照。用5% BSA封闭后,将细胞与针对内源性ZC3HAV1(小鼠多克隆抗体,由上海兽医研究所伴侣动物研究组惠赠)和FIPV N蛋白的一抗在4°C孵育过夜,随后与Alexa Fluor偶联的二抗孵育。细胞核用DAPI染色,使用激光扫描共聚焦显微镜采集图像。
谷胱甘肽S-转移酶(GST)pull-down实验
为确定ZC3HAV1与FIPV N之间的相互作用是否为直接作用,进行了体外GST pull-down实验。使用pET-28a-ZC3HAV1、pEGX-4T-1-FIPV N和空载pEGX-4T-1质粒在大肠杆菌Rosetta(DE3)细胞中表达重组ZC3HAV1和FIPV N蛋白。为提高原核表达过程中ZC3HAV1在上清中的产量,培养物在16°C下用1 mM IPTG诱导16 h。通过超声处理收获菌体沉淀,收集上清用于SDS-PAGE和WB分析。对于pull-down实验,将GST磁珠(货号C650031;BBI Life Sciences,上海,中国)与表达GST或GST-FIPV N的细菌裂解液孵育1 h,洗涤两次,然后与His-ZC3HAV1裂解液在25°C下再孵育1 h。洗涤五次后,通过在SDS上样缓冲液中煮沸洗脱结合蛋白,并进行WB分析。
FIPV N蛋白与ZC3HAV1相互作用结构域的鉴定
基于直接相互作用的证据,进一步进行实验以确定ZC3HAV1与FIPV N之间的结合区域。在一组实验中,将pCMV14-Myc-FIPV N与一系列ZC3HAV1截短突变体(Z1–Z6)共转染HEK293T细胞。在平行的一组中,将pCMV14-HA-ZC3HAV1与一组FIPV N截短构建体(N1–N5)共转染细胞。所有构建体均克隆至pCAGGS。两种实验中均纳入表达GFP的细胞(pCAGGS-GFP)作为阴性对照。转染后24 h,使用抗Myc或抗HA偶联磁珠对细胞裂解液进行免疫沉淀,随后进行Western blot分析。
此外,使用AlphaFold3在线服务器(https://alphafoldserver.com/)预测ZC3HAV1与FIPV N之间的蛋白–蛋白相互作用模型。使用预测局部距离差异检验(pLDDT)评分评估模型置信度。所得模型进一步使用PyMOL软件可视化,以评估结合结构域的空间定位。
稳定ZC3HAV1敲低细胞的构建及回补实验
使用慢病毒介导的shRNA方法构建稳定敲低ZC3HAV1的CRFK细胞。简言之,将pLKO.1-ZC3HAV1与psPAX2和pMD2.G共转染HEK293T细胞以产生慢病毒颗粒。收集含慢病毒颗粒的上清并用于转导CRFK细胞。转导后的细胞用嘌呤霉素(6 μg/mL)(货号ST551;Biotanon,上海,中国)筛选,并在嘌呤霉素选择压力下连续传代十次以构建稳定敲低细胞系。shRNA序列列于补充S1表。通过Western blot和RT-qPCR确认敲低效率。
为进行回补实验,将表达野生型ZC3HAV1或ZC3HAV1截短突变体Z5(该突变体缺失参与FIPV N结合的含RBD区域)的质粒转染至CRFK-shZC3HAV1细胞,并共表达或不共表达FIPV N蛋白。转染后16小时,以FIPV(MOI=0.01)感染细胞。在指定时间点收集样本,通过Western blotting、RT-qPCR和TCID50试验评估病毒复制。
为进一步评估FIPV N中与ZC3HAV1相关的区域是否有助于拮抗ZC3HAV1介导的抗病毒活性,将野生型ZC3HAV1与全长FIPV N或指定的N蛋白截短突变体(N1–N5)共转染至CRFK-shZC3HAV1细胞。转染相应空载体的细胞作为对照。转染后16小时,以FIPV(MOI=0.01)感染细胞。在指定时间点收集样本。通过Western blotting分析蛋白表达,同时分别使用FIPV特异性引物和TCID50试验通过RT-qPCR测定病毒RNA水平和感染性病毒滴度。
RNA结合缺陷型FIPV N突变体的构建与验证
为构建FIPV N蛋白的RNA结合缺陷型突变体,使用特异性设计的诱变引物通过重叠PCR将N端RNA结合结构域内的三个保守残基(Y93、R117和K123)替换为丙氨酸。将所得突变体(Mut-N;Y93A/R117A/K123A)克隆至pCMV-Myc表达载体,并通过DNA测序确认。引物序列列于补充S1表。
为验证RNA结合活性,在HEK293T细胞中表达野生型N和Mut-N蛋白。首先通过Co-IP试验检测Mut-N与ZC3HAV1之间的相互作用,以确定引入的突变是否影响蛋白–蛋白相互作用。随后按下文所述通过RNA免疫沉淀(RIP)试验评估RNA结合能力。通过RT-qPCR定量病毒RNA富集程度,并以input RNA进行归一化。
FIPV gRNA和sgRNA的qPCR检测方法的建立
为研究ZC3HAV1是否调控FIPV基因组RNA(gRNA)复制并影响sgRNA转录,基于冠状病毒不连续转录机制以及gRNA和所有sgRNA共有的5¢非翻译区(UTR)前导序列,建立了qPCR检测方法。首先以FIPV感染CRFK细胞的总RNA合成的cDNA为模板,使用覆盖整个5¢ UTR的多个正向引物(引物序列列于S4表)和靶向N基因开放阅读框(ORF)5¢近端区域的反向引物,通过PCR鉴定FIPV的5¢ UTR前导序列。扩增结果表明上游序列属于前导区域的一部分。
确认前导序列后,设计特异性sgRNA检测引物,正向引物位于5¢ UTR前导序列内,反向引物靶向目标ORF的5¢区域,确保仅扩增目标sgRNA。对于gRNA特异性检测,引物设计于pp1a/pp1b ORF的3¢区域(正向)和S基因ORF的5¢区域(反向),以避免扩增sgRNA。详细引物序列见S4表。
FIPV感染过程中sgRNA转录水平的鉴定
冠状病毒复制过程中sgRNA的合成遵循时间顺序,不同基因的sgRNA在病毒复制的不同阶段(包括早期、中期和晚期)可能表现出不同的表达模式。在后续分析ZC3HAV1对FIPV sgRNA与gRNA相对表达的影响时,FIPV复制过程中不同时间点sgRNA表达的差异可能对结果产生影响。为探究FIPV sgRNA的转录差异,我们进一步研究了CRFK细胞在感染后12和16小时(hpi)的sgRNA表达。将生长状态良好的CRFK细胞接种于12孔板中。待细胞达到100%汇合后,加入FIPV悬液(MOI=0.01),并设置未感染对照组。在12和16 hpi收集样本,提取RNA并逆转录。使用qPCR检测FIPV gRNA,并以sgRNA评估表达水平。
ZC3HAV1对FIPV亚基因组RNA表达影响的定量评估
为研究ZC3HAV1是否调控FIPV sgRNA的转录,以FIPV(MOI=0.01)感染CRFK细胞。在12和16 hpi收集总RNA。使用建立的区分gRNA和sgRNA的qPCR方法分析sgRNA的表达水平。为确定ZC3HAV1的具体调控作用,分别使用pCAGGS-ZC3HAV1和siRNA进行过表达和敲低实验。sgRNA转录水平在分别以GAPDH和病毒gRNA归一化后进行定量。具体而言,ΔCt值相对于GAPDH计算。使用ΔΔCt法将sgRNA水平以gRNA归一化:ΔΔCtsgRNA/gRNA = (ΔCtsgRNA − ΔCtgRNA)处理组−(ΔCtsgRNA − ΔCtgRNA)对照组。sgRNA的相对表达量以2−ΔΔCt sgRNA/gRNA计算,从而可独立于整体病毒RNA复制评估ZC3HAV1介导的调控。
基于GS-441524的FIPV RNA持续性检测
为确定ZC3HAV1是否独立于正在进行的病毒RNA合成来调控FIPV RNA的丰度,使用RNA依赖性RNA聚合酶抑制剂GS-441524阻断新生的病毒RNA合成。CRFK细胞感染FIPV(MOI = 0.01),并在感染后4小时加入终浓度为200 μM的GS-441524。未处理的感染组作为阳性对照。在加药后0、4、8和12小时收集细胞,提取总RNA,通过RT-qPCR定量病毒基因组RNA(gRNA)和亚基因组RNA(sgRNA),以此确定后续实验的最佳时间窗口。为研究ZC3HAV1的作用,平行进行了功能获得和功能缺失两种方法。对于过表达,CRFK细胞在感染前12小时转染pCMV14-HA-ZC3HAV1或空载体。对于敲低,细胞在感染前60小时转染siZC3HAV1或阴性对照siRNA(siNC)。随后细胞感染FIPV(MOI = 0.01),并在感染后4小时加入GS-441524(200 μM)以抑制新生的病毒RNA合成。在加药后10小时收集细胞,提取总RNA并通过RT-qPCR定量病毒gRNA和sgRNA。
FIPV中CpG分布分析及与其他冠状病毒的比较基因组分析
为探究ZC3HAV1是否通过CpG识别优先靶向FIPV RNA,我们分析了FIPV基因组中的CpG分布,并与其他冠状病毒进行了比较分析。从GenBank获取了FIPV和五种代表性冠状病毒的全长基因组序列——HCoV-229E(GenBank:KU291448.1)、CCoV(KP981644.1)、SARS-CoV-2(MZ054889.1)、IBV(OR824987.1)和PDCoV(MN942260.2)。使用R以100 bp非重叠窗口计算全基因组的CpG频率,并进行可视化以评估种间差异。对于FIPV,进一步使用滑动窗口分析计算了特定基因组区域(如5′ UTR、ORF1a/1ab、S、N)的CpG观察值/期望值(O/E)比值。区域分布以小提琴箱线图可视化,5′ UTR以10 bp窗口进行分析,以捕捉前导序列附近的CpG变异。
ZC3HAV1与FIPV基因组RNA和亚基因组RNA结合的RNA免疫沉淀(RIP)实验
为验证ZC3HAV1是否靶向冠状病毒gRNA和sgRNA,使用BeyoRIP RIP Assay Kit(P0107;Beyotime)按照制造商说明进行RIP。将CRFK细胞接种于6孔板,在90%融合时转染2 μg pCAGGS-ZC3HAV1、pCAGGS-ZC3HAV1-Z2(RBD结构域)、pCAGGS-ZC3HAV1-Z4(WWE+PARP样结构域)、pCMV14-Flag-FIPV N或pCAGGS-Flag,以未转染细胞作为对照。20小时后,细胞感染FIPV(MOI = 0.01)16小时。在补充有二硫苏糖醇(DTT)、RNase抑制剂和苯甲基磺酰氟(PMSF)的RIP裂解缓冲液中制备裂解物,随后使用预结合到Protein A/G琼脂糖珠(P2012;Beyotime)上的抗Flag抗体(F1804;Sigma-Aldrich)或小鼠IgG同型对照(M07001;Abmart,Berkeley Heights,NJ,USA)进行免疫沉淀。RNA在55°C洗脱,并使用RNAeasy Animal RNA Extraction Kit(R0027;Beyotime)纯化,保留20 μL用于WB输入验证。共沉淀的RNA进行逆转录并通过qPCR分析以检测FIPV N sgRNA。RIP样品中靶RNA的相对富集使用%Input和Fold Change方法计算:ΔCtRIP = CtRIP−(CtInput−log2(Input稀释倍数));%Input = 2^−ΔCtRIP;ΔΔCtRIP/Flag = CtRIP − CtFlag;Fold Change = 2^−ΔΔCtRIP/Flag。
竞争性RIP实验评估FIPV N蛋白对ZC3HAV1–病毒RNA相互作用的干扰
为探究FIPV N蛋白是否能拮抗ZC3HAV1与病毒sgRNA之间的相互作用,在体外进行了竞争性RIP实验。由于N蛋白结合病毒RNA且ZC3HAV1直接与N蛋白相互作用,对实验流程进行了优化以避免N蛋白结合FIPV RNA的干扰。简言之,将CRFK细胞转染pCAGGS-ZC3HAV1、pCMV14-Myc和pCMV14-Myc-N质粒,孵育24小时后收集。使用含DTT、RNase抑制剂和PMSF的裂解缓冲液裂解细胞,收集细胞裂解物。对于RNA制备,CRFK细胞感染FIPV(MOI = 0.01)16小时,使用RNAeasy Animal RNA Extraction Kit提取总RNA。Protein A/G琼脂糖预先与HA抗体孵育以捕获ZC3HAV1。然后将蛋白复合物与表达ZC3HAV1的细胞裂解物在25°C孵育1小时,随后洗涤。接着,实验组加入Myc-N蛋白,对照组加入Myc蛋白,混合物再孵育1小时。洗涤后,两组均加入FIPV RNA(5 µg),反应在4°C孵育4小时。使用ΔΔCt方法分析RIP RNA,计算实验组和对照组中FIPV基因组RNA(gRNA)和sgRNA的富集倍数变化。该实验设计能够评估FIPV N蛋白是否与ZC3HAV1竞争结合病毒RNA。
为确定FIPV N蛋白是否能够从预先形成的ZC3HAV1–RNA复合物中置换病毒RNA,使用预结合至Protein A/G琼脂糖珠的抗HA抗体,从转染CRFK细胞裂解液中免疫沉淀HA-ZC3HAV1蛋白复合物。转染空HA载体的细胞(HA对照)作为阴性对照。将从感染FIPV(MOI = 0.01)16 h的CRFK细胞中提取的总量5 µg RNA与珠结合的HA-ZC3HAV1蛋白在4°C孵育4 h,以允许形成ZC3HAV1–RNA复合物。洗涤后,加入等量的Myc-N、Myc-Mut-N或Myc蛋白。样品在25°C再孵育1 h。充分洗涤后,使用RNAeasy Animal RNA Extraction Kit从免疫沉淀珠中提取RNA。通过RT-qPCR定量病毒RNA,并将FIPV gRNA和sgRNA的丰度对输入RNA进行归一化。
为评估FIPV N蛋白是否能破坏ZC3HAV1与宿主mRNA的结合,使用HSPG2 mRNA(先前报道的ZC3HAV1相关宿主转录本)作为代表性宿主靶标。除用于复合物形成的RNA为从未感染CRFK细胞中提取的5 µg总RNA外,其余步骤相同。所有后续步骤如上所述进行。通过RT-qPCR定量HSPG2 mRNA水平。
FIPV N对ZC3HAV1介导的病毒抑制的拮抗效应分析
为评估FIPV N蛋白是否能拮抗ZC3HAV1的抗病毒作用,我们使用JEV作为模型。ZC3HAV1已被证明靶向JEV RNA的3¢ UTR以抑制复制。将pCAGGS-Flag、pCMV14-Flag-FIPV N、pCAGGS-ZC3HAV1或这些质粒的组合转染BHK-21细胞。12 h后,细胞感染JEV(MOI = 0.1),并在感染后12、24和36 h收集。进行WB分析以检测JEV NS1蛋白表达,指示病毒复制水平。该设置允许检查FIPV N蛋白是否能抵消ZC3HAV1的抗病毒活性。
AlphaFold3结构预测
使用AlphaFold3在线服务器预测FIPV N与ZC3HAV1的相互作用模型。预测结构的置信度使用AlphaFold3提供的预测局部距离差异测试(pLDDT)评分评估,pLDDT > 90表示非常高的置信度,70–90表示可信预测,50–70表示低置信度,< 50表示非常低置信度。pLDDT评分通过预测模型中的颜色编码可视化。所得结构使用PyMOL可视化。由于预测复合物未通过残基水平诱变或结构实验验证,AlphaFold3模型用于指导域映射结果的解释,而非定义实验确认的氨基酸接触。
统计分析
所有实验数据使用GraphPad Prism版本9(GraphPad Software,San Diego,CA,USA)进行统计分析。组间差异使用多重t检验评估,p < 0.05被认为具有统计学显著性。
结果
蛋白质组学分析鉴定ZC3HAV1为关键N相关宿主限制因子
为系统映射与FIPV核衣壳(N)蛋白相互作用的宿主因子,我们采用两种互补的免疫共沉淀(Co-IP)策略结合质谱(MS)(图1A)。在第一种方法中,使用抗N单克隆抗体10E7从感染FIPV的CRFK细胞中共沉淀内源性N相关宿主蛋白。同时,使用抗Flag抗体从转染细胞中免疫沉淀外源表达的Flag标签N,以在受控表达条件下捕获宿主相互作用因子。Western blot分析确认两种抗体均高效且特异性地拉下N及相关宿主蛋白,对照组中非特异性背景最小(图1B),验证了两种流程用于N相互作用宿主因子蛋白质组学筛选的可靠性。
[图:与FIPV N蛋白相关的宿主蛋白的鉴定与表征。(A) c的实验流程示意图]

无标记LC–MS分析分别在CRFK、CRFK+FIPV、pCMV14-FIPV N和pCMV14-GFP组中鉴定出2,197、1,975、3,589和3,630个蛋白条目(S5表)。值得注意的是,FIPV感染细胞中检测到的蛋白条目数比未感染细胞低约10%。这种减少可能归因于多种因素,包括病毒诱导的宿主蛋白丰度、溶解度或亚细胞定位的改变,免疫沉淀过程中的竞争或表位遮蔽效应,以及感染相关的样本复杂度或细胞应激变化。为确定候选N相关宿主蛋白,我们计算了每种蛋白相对于其相应对照的富集比。在FIPV感染的抗N IP样本相对于未感染CRFK对照IP样本中,或在Flag-FIPV N IP样本相对于Flag-GFP对照IP样本中,富集比>2的蛋白被保留为候选互作蛋白。仅在对照样本中检测到或相对于相应背景对照缺乏富集的蛋白作为非特异性背景蛋白被排除。然后对来自感染和转染数据集的过滤后蛋白列表进行整合和人工整理以去除冗余条目,最终获得98个高置信度N相关宿主蛋白用于后续功能分析。
为探索与FIPV N蛋白相互作用的宿主因子的生物学背景,我们使用STRING数据库构建了蛋白–蛋白相互作用(PPI)网络(图1C)。该网络揭示了一个高度互联的结构,以明确标注的RPL和RPS家族核糖体蛋白以及代表性翻译相关因子如EIF3D和EIF3M为主,表明N与宿主翻译机器广泛互作。若干RNA相关因子,包括ZC3HAV1、MOV10、HNRNPF、DDX18和RALY,也占据网络中的中心节点,提示病毒RNA处理与宿主RNA结合过程之间存在潜在联系。
基因本体(GO)富集分析进一步从三个功能维度表征了这些互作蛋白的分子特征。在细胞组分类别中,蛋白主要富集于胞质核糖体和核糖核蛋白复合物(图1D),表明许多N相关宿主因子是翻译组装体的组成部分。在分子功能类别中,富集以核糖体结构成分、RNA结合和核酸结合活性为主(图1E),反映了与RNA代谢和翻译调控的强关联。在生物学过程类别中,翻译、肽和酰胺代谢过程、胞质翻译及核糖体生物发生显著富集(图1F),揭示N互作蛋白在功能上汇聚于蛋白质合成和RNA加工通路。
总之,这些整合的蛋白质组学和生物信息学分析表明,FIPV N蛋白与参与翻译和RNA代谢的宿主因子密集网络发生相互作用——这两条通路常与内在转录后限制和抗病毒RNA监视相重叠。这提示N可能占据病毒RNA合成与宿主RNA调控系统之间的关键功能界面,潜在地使病毒能够将高效基因组复制与逃避细胞RNA水平限制相耦合。从MS分析鉴定的98个N相关宿主蛋白中,我们接下来基于其在RNA代谢、翻译调控或抗病毒防御中的注释功能优先选择了18个候选蛋白(S6表)。FIPV感染期间的转录谱分析揭示了特征性的双相模式——大多数基因在感染后12小时短暂上调,随后在24小时下调(图1G)。这种时间模式可能代表宿主初始稳定RNA代谢和限制病毒复制的努力,随后在病毒获得对宿主翻译机器的控制后被克服。
为评估其功能相关性,将18个候选基因分别克隆至pCAGGS中在CRFK细胞中过表达,随后进行FIPV攻毒和病毒滴度定量。六种RNA相关蛋白——DDX18、MOV10、RPL36、ZC3HAV1、HNRNPF和RALY——显著降低病毒复制(图1H),将其鉴定为N互作组中的推定限制因子。值得注意的是,其中若干因子,如MOV10和ZC3HAV1,已在其他RNA病毒中被认为参与转录后抗病毒限制[19,20],进一步支持FIPV利用宿主RNA监视网络组分的可能性。
其中,ZC3HAV1因其在多种病毒家族(包括冠状病毒)中公认的RNA靶向限制作用而尤为突出[18,21,22]。其作为N相关宿主蛋白和FIPV复制负调控因子的双重身份促使我们研究ZC3HAV1如何限制FIPV以及FIPV N如何拮抗这种限制。
ZC3HAV1通过不依赖IFN的RNA水平机制限制FIPV复制
在确定ZC3HAV1是一种强效的N结合限制因子后,我们接下来旨在验证其抗病毒活性并阐明其作用机制。为确定宿主限制因子ZC3HAV1是否影响FIPV复制,将ZC3HAV1表达质粒转染CRFK细胞,随后以FIPV感染(MOI = 0.01)。在指定时间点收集样本进行Western blotting、TCID₅₀测定和绝对RT–qPCR定量。Western blot分析显示,ZC3HAV1过表达显著抑制了FIPV复制,表现为感染后12、16和24小时N蛋白丰度降低(图2A–B)。随着ZC3HAV1质粒量的增加,还观察到N表达呈剂量依赖性降低(图2C–D)。绝对RT–qPCR证实病毒N转录本在8至32小时显著减少(图2E),而TCID₅₀测定进一步证明病毒滴度相应降低(图2F)。总之,这些数据表明ZC3HAV1是FIPV复制的强效抑制因子。
[图:ZC3HAV1通过不依赖IFN的抗病毒机制抑制FIPV复制。(A) FIPV N的Western blot分析]

FIPV感染后内源性ZC3HAV1 mRNA表达增加,在感染后12小时观察到最高诱导,且在感染后24小时仍可检测到持续升高(图2G)。这些数据表明ZC3HAV1对FIPV感染具有转录响应性。由于外源性表达的ZC3HAV1在感染后期蛋白丰度降低(图2A),我们进一步检查这种降低是否可归因于FIPV N介导的直接降解。将HA-ZC3HAV1与递增量的Flag-FIPV N共转染并未明显降低HA-ZC3HAV1蛋白水平,光密度分析显示HA-ZC3HAV1丰度无显著降低(S1A-B图)。这些结果不支持FIPV N单独直接诱导ZC3HAV1降解的模型。
我们还检测了FIPV感染期间内源性ZC3HAV1蛋白水平。与外源性表达ZC3HAV1观察到的后期降低相反,在感染后12、16和24小时,FIPV感染细胞中内源性ZC3HAV1蛋白较模拟感染细胞增加,同时作为生产性感染标志物的FIPV N表达也增加(S1C-D图)。结合图2G所示的ZC3HAV1 mRNA诱导,这些数据表明FIPV感染诱导而非抑制内源性ZC3HAV1表达。因此,外源性表达ZC3HAV1的后期降低可能反映转染后的晚期蛋白周转或感染相关的细胞蛋白稳态变化,而非FIPV N的直接降解。
为进一步验证其抑制角色,设计了三条siRNA以敲低ZC3HAV1表达。其中,siZC3HAV1–1在mRNA和蛋白水平均达到最高沉默效率(图2H–J)。敲低ZC3HAV1显著增强了FIPV复制,表现为N蛋白表达增加(图2K–L)、病毒RNA水平升高(图2M)和病毒滴度更高(图2N)。这些互补的过表达和敲低结果一致支持ZC3HAV1是FIPV复制的负调控因子。
由于ZC3HAV1是一种干扰素刺激基因(ISG),其抗病毒活性可能通过I型干扰素(IFN)通路激活间接产生,或通过内在的RNA靶向限制直接产生。为区分这些可能性,我们接下来检查ZC3HAV1的抗病毒效应是依赖于IFN信号激活,还是反映内在的转录后限制机制。双荧光素酶报告基因实验显示,ZC3HAV1和FIPV N蛋白均独立抑制293T细胞中仙台病毒(SeV)诱导的IFN-β启动子活性(图2O)。然而,在293T和CRFK细胞中的共表达实验显示,两种蛋白在调控IFN-β转录方面无协同或拮抗相互作用(图2P–Q)。这些结果表明,在该报告基因实验条件下,ZC3HAV1介导的FIPV复制抑制不太可能由I型干扰素(IFN)通路激活来解释。因此,我们接下来检查ZC3HAV1是否通过与病毒RNA直接结合来限制FIPV。
ZC3HAV1减少FIPV基因组和亚基因组RNA的积累与持续存在
鉴于其强效抗病毒活性,我们接下来研究ZC3HAV1是否通过直接干扰病毒转录本积累在RNA水平限制FIPV复制。为此,我们首先建立了一种定量检测系统,可区分基因组RNA(gRNA)和通过冠状病毒不连续转录机制产生的各亚基因组RNA(sgRNA)。如图3A所示,5’ UTR内的一系列正向引物(Lp1–Lp12)与位于N基因5’端的反向引物配对。以病毒cDNA为模板进行PCR扩增,使用Lp1–Lp5产生了明显产物(图3B),表明FIPV的前导序列位于该区域内。与代表性冠状病毒的序列比对显示,前导序列中的FIPV转录调控序列(TRS-L)为ACGAAA,并确定了完整的前导序列(图3C)。
[图:ZC3HAV1减少FIPV基因组RNA和亚基因组RNA的积累与持续存在。(A) 引物设计策略示意图]

基于该映射关系,我们针对每个sgRNA(S、3a、3c、E、M、N、7a和7b)设计了特异性引物组,并设计了一组用于gRNA检测的引物。每对引物均扩增出单一预期条带,证实了引物特异性,并建立了一种可靠的FIPV gRNA和sgRNA定量检测方法。
我们接下来通过将各sgRNA归一化至gRNA水平,分析了FIPV sgRNA的相对累积模式。由于gRNA反映病毒基因组的总体丰度并作为sgRNA合成的模板,sgRNA/gRNA比值提供了在校正总病毒RNA水平差异后sgRNA累积的归一化读数。在检测到的sgRNA中,N sgRNA的sgRNA/gRNA比值最高,其次为M和S,而若干辅助sgRNA的相对丰度较低(图3D)。该模式与经典冠状病毒不连续转录机制一致,即结构基因sgRNA的累积通常较辅助转录本更为显著。大多数转录本的sgRNA/gRNA比值在感染后12小时和16小时之间基本保持可比,仅存在适度的转录本特异性波动(图3D)。因此,图3D建立了FIPV感染期间sgRNA与gRNA关系的参考图谱,使我们能够评估ZC3HAV1是改变病毒RNA的总体累积还是选择性重塑sgRNA的相对谱。
为评估ZC3HAV1是否影响这些病毒RNA,在感染前将ZC3HAV1表达质粒或siRNA转染至CRFK细胞。ZC3HAV1过表达显著降低了12小时和16小时的FIPV gRNA水平(图3E),而其敲低则产生相反效应(图3F)。sgRNA也观察到类似模式:ZC3HAV1过表达显著降低、敲低则增加多个结构和辅助sgRNA的转录本丰度(图3G–H)。这些数据表明,ZC3HAV1降低基因组和亚基因组FIPV RNA的稳态累积。
由于gRNA作为sgRNA合成的模板,我们进一步将sgRNA丰度归一化至gRNA水平,以评估ZC3HAV1是否独立于病毒基因组总体丰度的变化而改变sgRNA的相对谱。尽管ZC3HAV1过表达降低、敲低增加gRNA和多个sgRNA的绝对水平(图3E–H),但大多数sgRNA种类的sgRNA/gRNA比值仅显示有限变化(图3I–J)。这些结果表明,ZC3HAV1主要降低FIPV RNA的总体累积,而非显著重塑各sgRNA的相对丰度。因此,sgRNA水平的下降很大程度上与病毒RNA丰度的广泛降低一致,而非sgRNA相对组成的的选择性改变。
然而,仅凭稳态RT-qPCR分析无法确定病毒RNA丰度的降低是源于病毒RNA合成受损还是已合成病毒RNA的持续性降低。为确定ZC3HAV1是否影响已合成病毒RNA的持续性,我们在FIPV RNA累积后进行了基于GS-441524的RNA持续性检测。GS-441524处理有效抑制了FIPV N RNA的持续累积,并随时间降低归一化sgRNA/gRNA表达,证实该条件下持续进行的RdRp依赖性RNA合成被强烈抑制(图S2A–B)。在这种复制被抑制的条件下,ZC3HAV1过表达降低了FIPV gRNA和多个sgRNA的水平,而ZC3HAV1敲低则增加了其丰度(图S2C–E)。归一化至gRNA后,ZC3HAV1过表达同样降低、而敲低则增加若干sgRNA的相对丰度(图S2F–G)。这些结果表明,ZC3HAV1介导的限制不能仅用病毒RNA合成减少或sgRNA转录改变来解释,而是涉及已合成FIPV基因组和亚基因组RNA持续性的降低。
有趣的是,ZC3HAV1对病毒RNA累积的抑制效应在感染后期时间点减弱,提示FIPV可能随着感染进展拮抗ZC3HAV1相关的RNA水平限制。总之,这些发现确定ZC3HAV1是一种限制因子,至少部分通过降低RNA合成发生后病毒RNA的持续性来限制FIPV基因组和亚基因组RNA的累积。尽管这些数据与ZC3HAV1相关病毒RNA的去稳定化或降解一致,但并未明确所涉及的具体下游RNA衰减机制。总之,这些发现确定ZC3HAV1是一种限制因子,主要通过降低RNA合成发生后病毒RNA的持续性来限制FIPV基因组和亚基因组RNA的累积。由于我们的蛋白质组学分析鉴定FIPV N为主要的ZC3HAV1相关病毒蛋白,该RNA水平限制表型为检验FIPV N是否拮抗ZC3HAV1介导的病毒RNA持续性调控提供了依据。
FIPV N与ZC3HAV1结合并拮抗ZC3HAV1介导的抗病毒限制
为验证我们关于FIPV N干扰ZC3HAV1相关RNA水平限制的假说,我们接下来检测了N是否与ZC3HAV1相互作用。在感染FIPV(MOI = 0.01)的CRFK细胞中,免疫共沉淀(Co-IP)实验显示,Flag标签的ZC3HAV1能有效沉淀内源性N蛋白,而IgG或载体对照中未检测到相互作用(图4A)。该相互作用在共表达Myc标签N和Flag-ZC3HAV1的293T细胞中进一步得到证实,表明它们的结合不局限于CRFK细胞背景(图4B)。重要的是,RNase处理并未消除该相互作用,提示该结合并非主要由RNA桥接介导(图4B,右)。
[图:ZC3HAV1通过其含RBD的结构域直接与FIPV N蛋白相互作用。(A) 免疫共沉淀分析]

为进一步检测ZC3HAV1与FIPV N能否在体外结合,将His-ZC3HAV1与含GST-N或GST对照蛋白的澄清可溶性细菌裂解液共孵育。GST pull-down后行免疫印迹显示,His-ZC3HAV1与GST-N共沉淀,而不与单独GST共沉淀,支持在本实验条件下ZC3HAV1与FIPV N之间存在体外结合(图4C)。共聚焦免疫荧光分析显示,过表达的ZC3HAV1与FIPV N在细胞质中存在部分重叠(图4E),这与该相互作用发生于细胞质(病毒RNA复制和RNA水平限制发生的场所)的可能性一致。为进一步评估在感染条件下是否也能观察到这种空间关联,我们检测了CRFK细胞中的内源性ZC3HAV1和感染来源的FIPV N。在未感染的CRFK细胞中,内源性ZC3HAV1呈细胞质染色,而未检测到FIPV N染色。FIPV感染后,N蛋白主要在细胞质中蓄积,并与内源性ZC3HAV1阳性信号显示部分重叠(S3图)。这些数据支持在过表达系统中观察到的ZC3HAV1–N空间关联在感染背景下的相关性。
为定位该相互作用所需的区域,我们利用两种蛋白的截短突变体进行了交互结构域定位。对于FIPV N,构建了六个涵盖其N端结构域(NTD)和C端结构域(CTD)不同组合的构建体(N1–N6)(图4D)。为排除过表达蛋白的非特异性沉淀,在Co-IP实验中纳入GFP作为无关蛋白阴性对照。Co-IP实验显示,ZC3HAV1与全长N及同时含NTD和CTD的N4片段高效相互作用,而GFP及缺失任一区域的N片段则显示弱或检测不到的沉淀(图4F)。因此,仅当NTD和CTD同时存在时才观察到与ZC3HAV1的高效结合。这些数据提示,NTD和CTD的结构完整性是与ZC3HAV1高效相互作用所必需的。由于任一结构域的缺失也可能影响N的折叠、寡聚化或RNA结合特性,这些截短结果被解释为支持性结构域定位证据,而非对最小结合界面的确定性鉴定。
类似地,设计了ZC3HAV1截短构建体以精细定位N相互作用区域。ZC3HAV1的RBD包含五个CCCH型锌指基序,命名为ZnF1–ZnF5,而ZC3HAV1截短构建体命名为Z1–Z6(图4G)。纳入GFP作为无关蛋白阴性对照以评估Co-IP实验中的非特异性沉淀。FIPV N能有效共沉淀全长ZC3HAV1及含ZnF1–ZnF5的RBD构建体,而其他截短构建体未检测到沉淀(图4H)。这些结果表明,ZC3HAV1中包含ZnF1–ZnF5的RBD区域是与FIPV N高效结合所必需的。
为进一步可视化推测的分子界面,生成了AlphaFold3预测的结构模型并在PyMOL中进行分析(图4I–J)。该模型提示,FIPV N蛋白的NTD和CTD可能形成一个推测的复合表面,位于ZC3HAV1 N端区域(涵盖RBD,残基227–634)附近。连接区L1和第五个CCCH型锌指基序位于预测接触区域附近,提供了与结构域定位结果一致的结构假说。由于该模型为计算预测,其用于指导对结构域定位数据的解释,而非作为确定性结构证据。
总之,这些发现确立了FIPV N与ZC3HAV1结合,且高效相互作用需要ZC3HAV1的含RBD区域和FIPV N的完整结构组织。我们接下来研究了该相互作用是否有助于FIPV N对抗ZC3HAV1介导的抗病毒限制的能力。
为进一步评估N–ZC3HAV1相互作用的功能相关性,我们在ZC3HAV1敲低的CRFK-shZC3HAV1细胞中进行了回补实验。在这些细胞中重建了全长HA-ZC3HAV1或HA-Z5(一种缺失参与FIPV N相互作用的RBD包含区的截短体),并在指定条件下将Myc-FIPV N与全长HA-ZC3HAV1共表达,随后进行FIPV感染。用全长ZC3HAV1回补恢复了抗病毒限制作用,表现为感染后12和24小时FIPV-M蛋白表达降低以及病毒滴度随时间下降(S4A–C图)。FIPV N的共表达部分逆转了全长ZC3HAV1的抑制作用,导致FIPV M蛋白丰度增加、病毒滴度高于单独回补ZC3HAV1组(S4A–C图)。这些结果支持FIPV N对ZC3HAV1介导的限制具有功能性拮抗作用。
相比之下,与全长ZC3HAV1相比,HA-Z5的抗病毒活性明显受损,表现为更高的FIPV M蛋白水平和更高的病毒滴度(S4A–C图)。这一发现表明,RBD包含区不仅是FIPV N结合所必需的,也是ZC3HAV1内在抗病毒功能所必需的。因此,HA-Z5应被解释为一种功能受损的截短体,而非严格意义上相互作用缺陷但仍具备限制能力的突变体。
为进一步确定FIPV N中与ZC3HAV1相关的区域是否参与其拮抗活性,我们利用相互作用定位实验中鉴定的N截短突变体进行了额外的功能分析(图4D和F)。在ZC3HAV1回补系统中评估了全长N和各N截短突变体。表达分析证实,在实验条件下可检测到N片段,从而允许对不同构建体进行功能比较(S4D图)。与全长N相比,各N截短突变体的拮抗活性因其保留的结构区域不同而有所差异。在截短突变体中,N4片段在定位实验中显示出高效的ZC3HAV1结合,表现出相对较强的缓解ZC3HAV1介导限制的能力;而ZC3HAV1结合弱或检测不到的片段则拮抗活性降低,以病毒蛋白蓄积和病毒RNA水平衡量(S4E和F图)。一致地,感染性病毒产生分析显示N截短突变体之间呈现相似的模式(S4G图)。这些结果提示,FIPV N中与ZC3HAV1相关的区域有助于有效拮抗ZC3HAV1介导的抗病毒限制。
总之,这些回补分析支持一个模型,即N–ZC3HAV1相互作用有助于FIPV N介导的对宿主RNA水平限制的拮抗,而N破坏ZC3HAV1–病毒RNA结合的机制将在下文进一步研究。
ZC3HAV1通过其RBD包含区靶向FIPV基因组和亚基因组RNA,尽管CpG频率较低
ZC3HAV1经典上被定义为一种CpG偏向性抗病毒因子,优先与富含CpG的病毒RNA结合并限制其蓄积。为确定这一经典CpG偏向性识别模型是否适用于FIPV,我们接下来检查了FIPV及相关冠状病毒的全基因组CpG分布模式。全基因组分析显示,FIPV的整体CpG频率明显较低,富集仅限于5′区域(图5A–B)。相比之下,PDCoV和PEDV等猪冠状病毒显示出更高且更均匀的CpG密度,而HCoV-229E和FIPV具有相似的模式,特征为单个5′端CpG簇且其他区域分布稀疏。
[图:尽管FIPV整体CpG频率较低,ZC3HAV1仍通过其含RBD的区域结合FIPV基因组RNA和亚基因组RNA]

对FIPV基因组更精细的窗口分析(图5C)证实,在5′非翻译区(UTR)附近存在局部CpG富集,而基因组其余部分显示明显的CpG耗竭。由于所有冠状病毒sgRNA共享源自该UTR的共同5′前导序列,我们推测ZC3HAV1可能通过富含CpG的前导序列识别sgRNA。然而,分析FIPV 5′UTR内的CpG组成(图5D)发现,CpG富集发生在经典前导–主体连接区之外,提示FIPV中ZC3HAV1–RNA识别并不能完全由含前导区域内的局部CpG富集来解释。
为通过实验验证ZC3HAV1是否结合FIPV RNA,我们在过表达不同ZC3HAV1结构域构建体——全长ZC3HAV1、分离的RNA结合结构域(RBD)或WWE+PARP区域——的FIPV感染CRFK细胞中进行了RNA免疫沉淀(RIP)实验,并设置FIPV N和Flag对照。Western blotting证实各Flag标签构建体均被有效沉淀(图5E)。随后对免疫沉淀RNA的RT-qPCR分析显示,全长ZC3HAV1及其RBD显著富集多种FIPV RNA,包括gRNA以及编码结构蛋白(S、M、N、E)和辅助蛋白(3a、3c、7a、7b)的sgRNA(图5F)。各sgRNA的富集信号存在差异,部分辅助sgRNA(如3c)回收率相对较低。这种差异可能反映细胞内转录本丰度不同、低拷贝RNA种类的检测灵敏度降低,或与含ZC3HAV1的核糖核蛋白复合体结合较弱。相比之下,WWE+PARP片段富集极少,与IgG和Flag对照相当。这些RIP数据支持ZC3HAV1在感染细胞中与FIPV RNA结合,但本身并不能确定具体的CpG基序、RNA结构或ZC3HAV1结合病毒RNA的下游命运。
与既往报道一致,即冠状病毒N蛋白在感染细胞中可与基因组和亚基因组病毒RNA结合,FLAG-FIPV N在RIP实验中也富集了多种FIPV RNA种类,包括gRNA和若干sgRNA(图5F)。由于RIP实验使用的是感染细胞裂解物,这些信号反映的是感染期间细胞内与N结合的病毒RNA。相对较强的sgRNA信号与不连续转录过程中产生的冠状病毒sgRNA在细胞内丰度较高相一致。
这些发现表明,尽管FIPV基因组整体CpG频率较低,ZC3HAV1仍通过其含RBD的区域与FIPV基因组和亚基因组RNA结合。鉴于FIPV基因组大部分区域CpG密度低,ZC3HAV1的识别可能涉及其他RNA特征,如局部RNA结构、亚基因组RNA丰度、局部序列背景、局部CpG富集区域或病毒核糖核蛋白组织,而非仅由整体CpG频率解释。因此,我们的数据支持ZC3HAV1以非经典或CpG受限模式结合,但并未确立完全独立于CpG。
FIPV N破坏ZC3HAV1与病毒RNA的结合以减弱RNA水平限制
在确定FIPV N与ZC3HAV1相互作用且ZC3HAV1与FIPV基因组和亚基因组RNA结合后,我们接下来研究了N是否干扰ZC3HAV1与病毒RNA的结合。我们首先进行了竞争性RIP实验,其中HA-ZC3HAV1预先结合到Protein A/G偶联的抗HA磁珠上,然后在RNA捕获前与Myc标签FIPV N或Myc对照孵育(图6A)。在这些条件下,N易于与磁珠结合的HA-ZC3HAV1结合(图6B),但其存在显著降低了ZC3HAV1 RIP中多种FIPV转录本的回收率,包括gRNA和若干sgRNA(图6C)。这些结果表明,预先结合的FIPV N可损害ZC3HAV1捕获FIPV基因组和亚基因组RNA的能力。然而,由于该实验主要检验N在RNA暴露前与ZC3HAV1结合是否限制后续RNA回收,因此未能确定N是否能破坏已形成的ZC3HAV1-病毒RNA复合体。
[图:FIPV N破坏ZC3HAV1与病毒RNA的结合,从而减弱ZC3HAV1介导的RNA水平限制。(A) 示意图]

为解决该问题,我们接下来使用WT FIPV N和RNA结合受损的N突变体进行了置换RIP实验。Mut-N通过将FIPV N中三个保守RNA结合残基Y93、R117和K123替换为丙氨酸(Y93A/R117A/K123A)而生成(S5A图)。在RNase处理条件下,Mut-N与HA-ZC3HAV1保持可检测的相互作用,与WT N相当,表明这些突变并未废除N与ZC3HAV1的结合(S5B图)。然而,RIP-qPCR分析显示,与WT N相比,Mut-N对FIPV gRNA和多种sgRNA的富集显著降低,证实该突变体病毒RNA结合受损但仍保留ZC3HAV1结合能力(S5C图)。
我们随后设计了置换RIP实验,以检验FIPV N在ZC3HAV1-RNA复合体已形成后是否能减少病毒RNA结合。在该实验中,磁珠结合的HA-ZC3HAV1首先与FIPV RNA或宿主RNA孵育以允许RNA结合,洗涤去除未结合RNA,随后在RIP-qPCR分析前加入Myc标签WT N或Mut-N(S5D图)。WT N显著降低了ZC3HAV1结合的FIPV gRNA和多种sgRNA水平,而Mut-N的置换活性明显较弱(S5E图)。这些结果表明,FIPV N可在ZC3HAV1-病毒RNA复合体形成后减少ZC3HAV1结合的病毒RNA,且有效破坏需要N的病毒RNA结合能力。
为确定该效应是否反映了对ZC3HAV1结合RNA的非特异性破坏,我们将HSPG2 mRNA作为代表性ZC3HAV1相关宿主RNA对照纳入实验。与病毒RNA结果不同,WT N和Mut-N均未显著降低ZC3HAV1与HSPG2 mRNA的结合(S5F图)。尽管对单一宿主转录本的分析不能排除对所有细胞RNA的影响,但该结果不支持在我们的实验条件下ZC3HAV1–RNA复合物发生普遍的非特异性解离。综合来看,竞争性RIP和置换RIP实验支持一种偏嗜病毒RNA的拮抗模型,即FIPV N通过协同的ZC3HAV1相互作用和病毒RNA结合活性干扰ZC3HAV1–病毒RNA结合。
为确定FIPV N介导的ZC3HAV1相关限制性减弱是否严格依赖于FIPV复制背景,我们使用日本脑炎病毒(JEV)作为异源正链RNA病毒模型。选择JEV是因为其在细胞质中复制,并提供了一个实验上易于操作的RNA病毒系统,可在独立于FIPV复制复合物的情况下评估ZC3HAV1相关抗病毒活性。在感染JEV的BHK-21细胞中,ZC3HAV1表达降低了JEV NS1蛋白积累,而FIPV N表达则与NS1蛋白水平升高相关(相对于表达ZC3HAV1的条件)(图6D–E)。当FIPV N与ZC3HAV1共表达时,与单独表达ZC3HAV1相比,JEV NS1蛋白表达部分恢复(图6F–G)。这些结果表明,FIPV N可在异源RNA病毒系统中在病毒蛋白积累水平上减弱ZC3HAV1相关限制性,提示这种反限制活性并非严格局限于FIPV复制背景。
总体而言,这些结果支持一种协同拮抗模型,即FIPV N利用蛋白质–蛋白质和蛋白质–RNA相互作用来对抗ZC3HAV1相关的RNA水平限制。N–ZC3HAV1相互作用可能将N定位于ZC3HAV1–病毒RNA界面,而N固有的病毒RNA结合活性有助于减少ZC3HAV1相关的病毒RNA。因此,FIPV N似乎并非仅通过空间位阻阻断ZC3HAV1含RBD区域来拮抗ZC3HAV1;相反,其可以破坏预先形成的ZC3HAV1–病毒RNA复合物,且该方式至少部分依赖于N介导的病毒RNA结合。HSPG2对照进一步支持该效应偏嗜病毒RNA,而非反映ZC3HAV1相关RNA的全面破坏。异源JEV实验提供了额外证据,表明FIPV N介导的ZC3HAV1相关抗病毒活性减弱可在FIPV感染系统之外观察到。
讨论
冠状病毒已进化出复杂的分子策略来逃避宿主RNA监视并维持病毒RNA的积累和持续存在。这些RNA监视和降解系统构成了细胞内防御的重要第一道防线,能够在干扰素信号激活之前识别并降解非自身RNA[23,24]。在冠状病毒结构蛋白中,核衣壳(N)蛋白在该过程中被越来越多地认为发挥重要作用,不仅参与RNA衣壳化和复制,还调节宿主抗病毒通路[25]。对于猫传染性腹膜炎病毒(FIPV)——一种导致猫严重且常为致死性全身性疾病的冠状病毒——了解N如何规避宿主RNA限制尤为重要,因为有效控制病毒RNA合成和周转与猫传染性腹膜炎(FIP)的病毒载量动态和免疫病理结局密切相关。因此,本研究以FIPV作为致病模型,剖析冠状病毒结构蛋白如何拮抗宿主RNA水平防御。通过蛋白质组学筛选和功能验证,我们鉴定出宿主抗病毒限制因子ZC3HAV1是FIPV N靶向的主要宿主因子,揭示了一种此前未被认识到的转录后反限制机制。
蛋白质组学分析显示,FIPV N蛋白与广泛的宿主RNA结合和翻译调控因子网络相互作用,形成的相互作用图谱与SARS-CoV-2[26]和小鼠肝炎病毒(MHV)[27]中描述的图谱部分重叠。在这些病毒中,N已被证明与核糖体蛋白和应激颗粒组分相互作用,重塑宿主RNA调控环境以利于病毒翻译和RNA稳定性[28,29]。功能验证证实,这些N相关宿主蛋白中的若干成员——包括ZC3HAV1和MOV10,二者均以转录后抗病毒作用著称[30,31]——限制了FIPV复制。在这些候选因子中,ZC3HAV1尤其值得注意,因为其已被确立为RNA结合限制因子,促使我们进一步明确ZC3HAV1如何限制FIPV以及FIPV N如何对抗这种限制。
鉴于ZC3HAV1作为RNA结合抗病毒因子的已知作用,我们接下来检测了ZC3HAV1限制FIPV复制的水平。ZC3HAV1同时降低了基因组和亚基因组病毒RNA的积累,这与其先前报道的对SARS-CoV-2的限制作用一致[32]。重要的是,仅通过RT-qPCR测定的稳态gRNA和sgRNA水平降低,无法确定ZC3HAV1是通过抑制病毒RNA合成还是通过降低已合成病毒RNA的稳定性来发挥作用。为区分这两种可能性,我们进行了基于GS-441524的RNA持久性实验。由于GS-441524强烈抑制冠状病毒RNA依赖性RNA聚合酶活性[33,34],在病毒RNA已积累后进行处理可限制持续的病毒RNA合成,从而能够评估预先存在的FIPV RNA的持久性。在此条件下,ZC3HAV1加速了FIPV基因组和亚基因组RNA的丢失,表明ZC3HAV1介导的限制主要并非归因于病毒RNA合成减少或sgRNA转录降低。相反,这些数据支持一种合成后RNA去稳定化机制,即ZC3HAV1在FIPV RNA产生后降低其持久性。尽管这一发现与ZC3HAV1结合的病毒RNA的直接或辅因子辅助的去稳定化一致,但参与FIPV RNA丢失的具体RNA降解途径和宿主辅因子仍有待鉴定。这一机制上的细化很重要,因为它将ZC3HAV1的限制作用定位在病毒RNA持久性层面,并为FIPV N蛋白通过破坏ZC3HAV1与病毒RNA的关联从而在感染期间维持病毒RNA积累提供了直接的理论依据。
在确定ZC3HAV1在RNA水平上限制FIPV复制后,我们接下来试图阐明其如何识别病毒转录本。值得注意的是,尽管FIPV基因组CpG含量较低——这是包括SARS-CoV-2在内的许多冠状病毒共有的特征[17,35]——ZC3HAV1仍能有效限制FIPV。在我们的CpG谱分析中,FIPV在大部分基因组中表现出广泛的CpG缺失,仅在5′区域附近有局部富集。对5′ UTR的进一步定位显示,这种CpG富集与冠状病毒sgRNA共有的经典前导-主体连接处并不重合。这些观察结果表明,全局CpG频率、甚至含前导序列区域内局部CpG分布,都不足以完全解释ZC3HAV1相关的FIPV RNA限制。这一解释与先前研究一致,即ZC3HAV1介导的限制并不总是能由全局CpG丰度严格预测,且ZC3HAV1可限制CpG频率相对较低的病毒或通过情境依赖的RNA识别模式发挥作用[36,37]。因此,我们的数据并非证明完全独立于CpG,而是支持在CpG受抑制的冠状病毒RNA背景下ZC3HAV1的一种非经典结合模式。
我们的RIP实验进一步显示,全长ZC3HAV1及其含RBD的区域富集了FIPV gRNA和多种sgRNA,而WWE + PARP区域则显示极少的富集。这些数据表明,含RBD的区域是ZC3HAV1在感染细胞中与FIPV RNA有效结合所必需的。值得注意的是,编码结构蛋白(包括N和S)的几种sgRNA被强烈富集,这可能反映了它们较高的细胞内丰度和/或促进ZC3HAV1结合的情境特征。由于RIP不能确定ZC3HAV1识别的精确RNA基序或结构,其他特征——包括局部CpG富集区域、RNA丰度、二级结构、局部序列背景或病毒核糖核蛋白组织——可能有助于ZC3HAV1与FIPV RNA的结合[15,38]。这为理解FIPV尽管整体CpG缺失但仍对ZC3HAV1相关的RNA水平限制敏感提供了机制框架。
FIPV 5′区域附近的局部CpG富集也可能有助于ZC3HAV1的结合。尽管该富集区域与经典前导-主体连接处并不精确重叠,但它可能通过局部RNA结构、共享的5′ RNA架构或病毒核糖核蛋白组织影响ZC3HAV1的识别。因此,我们的数据并不排除5′ CpG富集区域的作用,而是表明FIPV RNA的识别不能仅由全局CpG频率或简单的前导连接处CpG模型来解释。对5′ CpG富集区域进行靶向突变定位、含前导序列的报告RNA或FIPV微型基因组/复制子系统将有助于确定该区域是否直接参与ZC3HAV1的结合和限制。
这些观察结果还表明,FIPV感染期间ZC3HAV1的识别不太可能局限于单一病毒RNA种类,也不大可能仅由简单的整体CpG阈值决定。FIPV N能够在异源JEV系统中减弱ZC3HAV1介导的限制,结合此前关于PEDV N在流感病毒系统中的观察结果[18],支持ZC3HAV1能够在不同实验背景下与病毒RNA结合的观点。在此背景下,我们的RIP-qPCR数据确定FIPV gRNA和sgRNAs是FIPV感染期间与ZC3HAV1相关的相关病毒RNA种类。重要的是,在置换RIP实验中,HSPG2被纳入作为代表性的ZC3HAV1相关非病毒RNA对照。FIPV N有效减少ZC3HAV1相关的FIPV RNA,但未显著降低ZC3HAV1与HSPG2的结合,这一发现表明N介导的置换并非对所有ZC3HAV1–RNA复合物的普遍非特异性破坏。相反,该结果支持一种偏嗜病毒RNA的拮抗模型,即FIPV N优先干扰ZC3HAV1与病毒RNA的结合,可能通过协调的N–ZC3HAV1相互作用和N介导的病毒RNA结合实现。然而,仍需采用ZC3HAV1 RIP-seq、CLIP-seq或相关RNA相互作用组分析等转录组范围的方法,以确定冠状病毒感染期间ZC3HAV1识别的病毒RNA底物及结合特征的全貌[38]。
我们的数据进一步细化了FIPV N拮抗ZC3HAV1介导的限制的机制。简单的竞争模型——即N仅遮蔽ZC3HAV1的RNA结合区域并阻止病毒RNA接近——不足以解释置换RIP结果。在重新设计的实验中,野生型N在ZC3HAV1–病毒RNA复合物已经形成后减少了ZC3HAV1相关的FIPV RNA,而RNA结合受损的Mut-N保留了可检测的ZC3HAV1相互作用,但置换活性明显减弱。这些结果表明,N介导的拮抗不仅仅需要与ZC3HAV1的蛋白质–蛋白质结合。相反,FIPV N可能通过一种协调机制发挥作用,即其与ZC3HAV1的相互作用将N定位到ZC3HAV1–病毒RNA界面,而其内在的病毒RNA结合活性有助于减少ZC3HAV1相关的病毒RNA。该模型也解释了为何破坏ZC3HAV1–病毒RNA结合而非降解ZC3HAV1本身,是FIPV N介导的主要拮抗机制。在机制上,这种拮抗扩展了已知的病毒对抗ZC3HAV1的策略谱。我们证明FIPV N蛋白与ZC3HAV1的含RBD区域结合,并损害其与病毒RNA的结合。结合置换RIP和Mut-N分析,这些发现支持一种协调的拮抗机制,即FIPV N同时利用ZC3HAV1结合和内在病毒RNA结合活性来减少ZC3HAV1相关病毒RNA的回收。该机制不同于先前描述的拮抗剂,如甲型流感病毒NS1,后者通过破坏辅因子依赖途径间接抑制ZC3HAV1/ZAP介导的RNA水平限制[39];也不同于痘苗病毒C16,后者将ZC3HAV1/ZAP螯合至胞质斑点中以减弱其抗病毒活性[40]。据我们所知,这些发现提供了证据表明,冠状病毒结构蛋白可通过靶向ZC3HAV1–病毒RNA轴来对抗ZC3HAV1相关的RNA限制,而非通过蛋白水解或转录抑制ZC3HAV1。通过这些偶联活性,N可能不仅作为RNA伴侣发挥作用,还作为病毒拮抗剂帮助在ZC3HAV1介导的限制下维持病毒RNA的积累和持续存在。该机制揭示了一种此前未被充分认识的RNA水平反限制策略,FIPV N借此在感染期间帮助维持病毒RNA积累(图7)。
[图:FIPV N介导拮抗ZC3HAV1相关RNA水平限制的推测模型。ZC3HAV1与FIP]

一个尚待解决的问题是,哪些下游宿主辅因子执行ZC3HAV1相关的FIPV RNA限制作用。在其他病毒系统中,据报道ZC3HAV1与RNA监控或降解相关因子协同作用,包括TRIM25、KHNYN、OAS3及相关宿主蛋白,以促进抗病毒RNA限制[41–43]。尽管这些辅因子为未来研究提供了重要候选,但本研究聚焦于界定使ZC3HAV1介导的限制和FIPV N介导的反限制得以实现的上游病毒–宿主界面。具体而言,ZC3HAV1与FIPV基因组和亚基因组RNA结合,减少其积累和持续存在,而FIPV N通过破坏ZC3HAV1–病毒RNA结合来对抗这一过程。通过结合RNA持续性试验、RIP-qPCR和置换RIP分析,我们的研究超越了稳态病毒RNA定量,将ZC3HAV1相关的RNA限制与特定的病毒拮抗机制联系起来。进一步的机制解析将需要靶向敲低或回补候选辅因子、更广泛的相互作用组分析,以及以RNA为中心的方法,如FIPV复制子系统、转录缺陷突变体、代谢RNA标记或Northern印迹。这些方法将有助于确定是经典ZC3HAV1辅因子还是FIPV特异性RNA降解途径执行病毒RNA丢失,并将进一步阐明病毒RNA合成、加工和降解对ZC3HAV1介导限制的相对贡献。
为进一步检验FIPV N介导的ZC3HAV1相关限制减弱是否严格局限于FIPV复制背景,我们使用JEV作为异源正链RNA病毒模型。FIPV N部分缓解ZC3HAV1介导的JEV限制的能力表明,这种反限制活性可在FIPV感染系统之外观察到。这一发现支持FIPV N介导的ZC3HAV1拮抗的功能可转移性,并表明该活性并非完全依赖于FIPV特异性复制机器。然而,还需要更多冠状病毒和不相关RNA病毒来确定FIPV N可对抗ZC3HAV1相关限制的病毒RNA背景广度。
总体而言,这些发现表明FIPV N蛋白除了在病毒RNA包装、复制支持和免疫调节中的已知功能外[44,45],还具有额外的RNA保护作用。这项工作将冠状病毒N蛋白的经典RNA结合功能与细胞RNA水平限制的逃逸联系起来,并鉴定N介导的ZC3HAV1–病毒RNA界面破坏为一种此前未被充分认识的冠状病毒反限制模式。
结论
总之,本研究以临床致死性猫冠状病毒为模型,揭示了一种RNA水平反限制机制,其中核衣壳蛋白通过其与含RNA结合结构域区域的相互作用结合ZC3HAV1,并在功能上对抗ZC3HAV1介导的抗病毒限制。N截短突变体的功能分析进一步表明,FIPV N的ZC3HAV1相关区域有助于有效拮抗ZC3HAV1介导的限制。此外,FIPV N干扰ZC3HAV1–病毒RNA结合,从而帮助维持病毒RNA的积累和持续存在。这些发现为FIPV如何在宿主RNA结合限制因子存在的情况下帮助维持丰富的基因组和亚基因组RNA提供了机制解释,并将冠状病毒N蛋白的功能作用从病毒RNA包装和复制支持扩展到对宿主RNA水平限制的主动调节。重要的是,我们的数据支持一种模型,即ZC3HAV1介导的限制至少部分作用于病毒RNA生命周期的合成后阶段,限制已合成的FIPV RNA的持续存在,而非仅由病毒RNA合成减少或sgRNA转录改变来解释。尽管确切的下游RNA降解机器仍有待界定,FIPV N对ZC3HAV1–病毒RNA限制轴的协同靶向代表了一个特定的宿主–病毒界面,可能为未来旨在恢复内在RNA限制和限制冠状病毒RNA持续存在的抗病毒策略提供概念基础。
支持信息
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
致谢
我们诚挚感谢中国农业科学院上海兽医研究所公共实验室提供实验设施。我们也衷心感谢国家自然科学基金会的财务支持。
数据可用性
所有相关数据均在手稿及其支持信息文件中。
资助声明
本研究得到国家自然科学基金(32172832,负责人GQL)、国家重点研发计划(2023YFD1800705,负责人GQL)、中央级公益性科研院所基本科研业务费专项资金(2026JB09,负责人YQZ)、中国博士后科学基金(2024M753586,负责人YQZ)、中央级公益性科研院所基本科研业务费专项资金(2025JB01 AXT)以及上海市科技重大专项(ZD2021CY001,负责人GQL)的资助。资助方在研究设计、数据收集与分析、发表决定或稿件撰写中均未发挥作用。
参考文献
[内容请查看原文]
关联数据
补充材料
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
数据可用性声明
所有相关数据均包含在稿件及其补充信息文件中。
查看英文原文▾
Introduction
Feline infectious peritonitis (FIP) is one of the most devastating viral diseases in cats, caused by the feline infectious peritonitis virus (FIPV), a highly virulent biotype of feline coronavirus [1]. Despite decades of research, FIP remains almost uniformly fatal, with limited effective treatments available, posing a significant burden on feline health and veterinary care. Clinically, FIP manifests as pyogranulomatous inflammation, effusive peritonitis, and multiorgan failure, reflecting widespread endothelial damage and immune dysregulation [2,3]. Given its high lethality and systemic nature, FIPV provides not only an urgent clinical target but also an excellent model for dissecting coronavirus–host interactions, particularly those that determine viral RNA accumulation and host restriction evasion.
Coronaviruses have evolved intricate molecular strategies to manipulate host cells and ensure efficient replication. A critical frontier in this virus–host arms race lies in the intrinsic cellular RNA restriction system, which detects and degrades foreign RNA transcripts to prevent infection [4,5]. However, the full spectrum of viral countermeasures against this defense remains incompletely defined. Understanding how pathogenic coronaviruses, such as FIPV, subvert these RNA surveillance pathways is therefore important for explaining how they maintain viral RNA accumulation, efficient replication, and infection-associated pathogenesis.
Among coronavirus structural proteins, the nucleocapsid (N) protein has emerged as a multifunctional regulator that facilitates viral RNA synthesis by associating with both viral RNA and components of the replication–transcription complex (RTC) [6]. Through its high-affinity RNA-binding activity and interactions with nonstructural proteins such as nsp3 and nsp12, N stabilizes viral RNA templates, promotes RTC assembly, and enhances the efficiency of genome replication and subgenomic transcription [7–9]. Beyond this canonical role, the N protein has been increasingly implicated in modulating host RNA metabolism [10,11] and antiviral signaling [12,13], suggesting that it may also influence the susceptibility of viral RNAs to host RNA surveillance pathways. However, whether this putative RNA-protective function involves counteraction of host cytoplasmic RNA-level restriction, a crucial antiviral pathway that acts as an early intracellular barrier to viral replication [14], remains poorly understood. Elucidating this link would provide key insight into how coronavirus N proteins integrate viral RNA synthesis with the evasion of host RNA restriction systems.
To address this question, we performed an unbiased co-immunoprecipitation coupled with mass spectrometry (Co-IP/MS) screen in FIPV-infected cells to identify host proteins that associate with N. The analysis revealed extensive interactions with host RNA-binding factors, among which the restriction factor ZC3HAV1 (also known as ZAP or PARP13) emerged as a prominent candidate. ZC3HAV1 is an interferon-stimulated gene product that binds viral RNA and, in several viral systems, suppresses replication by recruiting RNA-surveillance and decay-associated factors [15,16]. Accumulating evidence has established ZC3HAV1 as an antiviral restriction factor against diverse viruses, including several positive-strand RNA viruses and coronaviruses. In coronaviruses, ZC3HAV1-mediated restriction has been reported for SARS-CoV-2 and porcine epidemic diarrhea virus (PEDV) [17,18], indicating that coronavirus RNAs can be susceptible to ZC3HAV1-associated RNA-level restriction despite the generally CpG-suppressed nature of many coronavirus genomes. However, the basis by which low-CpG coronavirus RNAs remain susceptible to ZC3HAV1-associated restriction remains incompletely understood, particularly with regard to whether residual CpG distribution, other RNA-context features, or virus-specific RNA-binding proteins influence this process. Moreover, previous coronavirus studies have mainly focused on the antiviral activity of ZC3HAV1 itself, whereas how coronaviruses counteract ZC3HAV1-mediated RNA restriction remains much less understood. In particular, whether coronavirus structural proteins can interfere with the ZC3HAV1–viral RNA axis has not been fully defined. Notably, FIPV possesses a low-CpG genome, yet ZC3HAV1 is classically defined as a CpG-biased RNA-binding antiviral factor [17]. This apparent paradox makes FIPV a useful system to investigate how a CpG-depleted coronavirus RNA remains sensitive to ZC3HAV1-associated RNA-level restriction and whether the coronavirus N protein preserves viral RNA accumulation and persistence by antagonizing the ZC3HAV1–viral RNA axis.
Our preliminary data suggested that ZC3HAV1 interacts with the FIPV N protein and acts as a negative regulator of viral replication. We hypothesized that FIPV N antagonizes ZC3HAV1-associated RNA-level restriction by reducing ZC3HAV1 association with viral RNA, thereby preserving viral RNA accumulation and persistence. By focusing on a clinically fatal feline coronavirus, this study not only investigates an RNA-level counter-restriction mechanism but also offers broader insight into how coronaviruses preserve transcript stability to sustain infection. We identify the ZC3HAV1 RBD-containing region as a key interface involved in FIPV N association, demonstrate that FIPV N impairs ZC3HAV1–viral RNA association, and show that ZC3HAV1 recognition of FIPV RNA is not fully explained by overall CpG abundance. Collectively, these findings bridge veterinary and molecular virology by revealing how coronavirus structural proteins subvert intrinsic RNA restriction, thereby expanding our understanding of coronavirus–host interplay at the post-transcriptional level.
Materials and methods
Virus, vectors, cells, and animals
SeV, FIPV strain SH2021 (GenBank accession number OR295209.1), and Japanese encephalitis virus (JEV) were provided by the Companion Animal Diseases Team, Shanghai Veterinary Research Institute (SHVRI). Plasmids pCold-SUMO-FIPV-N, pET-28a, pLKO.1, psPAX2, pMD2.G, pCAGGS, pCMV14-Flag, pCMV14-hemmagultinin (HA), pCMV14-Myc, pCMV14-green fluorescent protein (GFP), pGAS-Luc-IFN-β (human), pGL3-IFN-β (feline), and pRL-TK were also provided by SHVRI. CRFK and HEK293T cells were purchased from the Chinese Academy of Sciences Cell Bank. BHK-21 cells were provided by the Companion Animal Diseases Team at SHVRI.
Plasmid construction
Plasmids were constructed by homologous recombination using the One Step Cloning Kit (Vazyme Biotech Co., Ltd., Nanjing, China). The N gene was amplified from the pCold-SUMO-FIPV-N plasmid, and host genes from cDNA synthesized from CRFK cell RNA. PCR was performed with 2 × Phanta Max Master Mix (Vazyme) following the manufacturer’s instructions. Plasmids were linearized by double digestion, and inserts were assembled by homologous recombination. Prokaryotic expression plasmids were transformed into Escherichia coli Rosetta (DE3), and eukaryotic plasmids into E. coli DH5α. Primer sequences are listed in S1 Table.
Western blot (WB) analysis
Cells were seeded and transfected with plasmids using Lipofectamine 3000 (Thermo Fisher Scientific, Waltham, MA, USA) according to the manufacturer’s instructions when they reached 80–90% confluency at 16–24 h after plating. After 12 h, the medium was replaced with maintenance medium (except for 293T cells). Cells or virus-infected samples were collected at the indicated time points and lysed in RIPA buffer on ice for 30 min, followed by centrifugation at 12,000 rpm for 5 min at 4°C. Supernatants were mixed with 4 × SDS loading buffer, boiled for 10 min, and stored at −40°C. Proteins were separated by SDS-PAGE, transferred to 0.22 μm PVDF membranes (Merck Millipore, Burlington, MA, USA), blocked with 5% skim milk, and incubated with primary antibodies (as recommended by the manufacturer) overnight at 4°C or for 1 h at 25°C. After incubation with secondary antibodies, signals were detected using ECL substrate and visualized with an imaging system (Cat#4600; Biotanon, Shanghai, China).
Determination of 50% tissue culture infectious dose (TCID50)
To determine the TCID50 of FIPV, whole-cell lysates were collected and subjected to three freeze-thaw cycles, followed by centrifugation at 3,000 × g for 10 min at 4°C. The supernatants were collected and serially diluted 10-fold (10-1 to 10-8) in Eagle’s minimum essential medium. CRFK cells seeded in 96-well plates were washed to remove the culture medium. One hundred microliters of each virus dilution was added to each well using eight replicates per dilution. Negative control wells were included without virus inoculation. Plates were incubated at 37°C in a 5% CO2 atmosphere for 5 days. The cytopathic effect (CPE) was recorded, and the TCID50 value was calculated using the Reed–Muench method.
RT-qPCR
To investigate host transcriptional responses to FIPV infection, CRFK cells were infected with FIPV (multiplicity of infection [MOI] = 0.01) or mock-treated. Total RNA was extracted at 12 and 24 h post-infection (hpi) using TRIzol reagent (Accurate Biotechnology Co., Ltd., Shanghai, China), reverse transcribed into cDNA with Evo M-MLV (Accurate Biotechnology Co., Ltd.), and subjected to qPCR using SYBR Green Pro Taq HS Mix (Accurate Biotechnology Co., Ltd.) and gene-specific primers (S2 Table). Relative mRNA levels were normalized to glyceraldehyde 3-phosphate dehydrogenase (GAPDH) and calculated using the 2−ΔΔCt method.
For absolute quantification, FIPV RNA was quantified using primer sets targeting either the N gene or the untranslated region (UTR). For N gene-based quantification, the primers were forward 5′-TAAGGAACTCGCTGAGAGGTGG-3′ and reverse 5′-CTGAAACTGTGGCGGTATCTTA-3′. A standard curve ranging from 10² to 10⁸ copies/μL was generated (y = −3.3351x + 37.367, R² = 0.9958). For UTR-based quantification, the primers were forward 5′-TCCGAGGAATTACTGGTCATC-3′ and reverse 5′-CACTAGATCCAGACGTTAGCTC-3′. A standard curve ranging from 10² to 10⁸ copies/μL was generated (y = −3.3141x + 36.214, R² = 0.9947). Viral RNA loads were calculated from the corresponding standard curves and expressed as copies per μg of total RNA.
Immunoprecipitation (IP) and co-immunoprecipitation (Co-IP)
IP was performed using either anti-N monoclonal antibodies or tag-based magnetic beads. For direct IP, CRFK cells infected with FIPV (MOI = 0.01) were lysed 16 hpi, and cell lysates were incubated with anti-FIPV N monoclonal antibody (clone 10E7), followed by Protein A/G agarose beads (Cat# 36403ES25; Yeasen Biotech Co., Ltd., Shanghai, China), according to the manufacturer’s instructions. Eluted proteins were analyzed by WB and MS.
For Co-IP, CRFK or 293T cells were transfected with tagged protein constructs using Lipofectamine 3000. Where indicated, CRFK cells were infected with FIPV after transfection. At 24 h post-transfection, cells were lysed, and tagged proteins were immunoprecipitated using anti-tag magnetic beads (ShareBio, Shanghai, China) following the manufacturer’s instructions. After washing, bound proteins were eluted for SDS-PAGE, followed by WB or MS.
Liquid chromatography (LC)-MS
The target gel bands were excised and subjected to in-gel tryptic digestion. Peptides were analyzed by LC-tandem MS (LC-MS/MS) at Shanghai APT Biotechnology Co., Ltd. (Shanghai, China) for protein identification.
Screening of host proteins
To identify host proteins interacting with FIPV N protein, IP MS data of N protein complexes were analyzed. Highly abundant host proteins were selected for further study. Overexpression plasmids for the host proteins were constructed and transfected into CRFK cells. After culturing CRFK cells to 80%–90% confluence, plasmids expressing different host proteins or an empty pCAGGS control plasmid were transfected using Lipofectamine 3000. Following a 12-h transfection, the culture medium was replaced with maintenance medium. After 24 h, the cells were washed and infected with FIPV at an MOI of 0.01. After 16 h of infection, cells underwent freeze-thaw cycles, and the supernatant was collected by centrifugation. FIPV titers in the supernatant were determined using the TCID50 method.
To analyze the effect of FIPV infection on host gene expression, qPCR was performed to measure transcriptional changes in host genes during the infection process. CRFK cells were seeded in 6-well plates and infected with FIPV at an MOI of 0.01, with uninfected cells serving as controls. Samples were collected at 12 and 24 h post-infection (hpi). Total RNA was extracted using TRIZOL, and cDNA was synthesized using the Evo M-MLV reverse transcription kit. qPCR was then conducted using the SYBR Green Pro Taq HS qPCR kit to assess the transcriptional levels of host genes potentially involved in FIPV replication.
Effect of ZC3HAV1 overexpression on FIPV replication
To investigate the role of ZC3HAV1 in FIPV replication, CRFK cells were transfected with pCAGGS-ZC3HAV1 or control plasmid pCAGGS-Flag. At 12 h post-transfection, the medium was replaced, and cells were subsequently infected with FIPV (MOI = 0.01). At the indicated time points post-infection, cells were harvested for analysis of FIPV N protein by Western blotting, viral RNA levels by absolute RT-qPCR, and viral titers by TCID50 assay. In addition, dose-dependent effects of ZC3HAV1 were assessed by transfecting increasing amounts of ZC3HAV1 expression plasmid, followed by analysis of FIPV N protein at 16 h post-infection.
To evaluate whether FIPV infection affects endogenous ZC3HAV1 expression, CRFK cells were infected with FIPV (MOI = 0.01) and collected at the indicated time points for analysis by RT-qPCR and Western blotting. Furthermore, to assess whether FIPV N protein directly regulates ZC3HAV1 expression, CRFK cells were co-transfected with Flag-ZC3HAV1 and increasing amounts of HA-tagged N plasmids, and ZC3HAV1 protein levels were determined by Western blotting.
Effect of knockdown of endogenous ZC3HAV1 on FIPV replication
Three small interfering RNA (siRNA) sequences targeting ZC3HAV1 were designed and synthesized. Their sequences are listed in S3 Table. CRFK cells were transfected with each siRNA or a non-targeting control (siNC), and knockdown efficiency was evaluated by qPCR at 48 h post-transfection and by WB at 60 h to identify the most effective siRNA. To assess its effect on FIPV replication, CRFK cells were transfected with the optimal siRNA or siNC, and infected with FIPV (MOI = 0.01) 48 h later. Samples were collected at 12, 16, and 24 hpi for WB analysis of N protein and endogenous ZC3HAV1 expression. Viral RNA levels and titers were quantified by qPCR and TCID50 assay, respectively.
Dual-luciferase reporter gene assay
Effect of ZC3HAV1 and FIPV N on IFN-β transcription: 293T cells were seeded in 24-well plates and co-transfected with pGAS-Luc-IFN-β (60 ng), pRL-TK (10 ng), and either pCAGGS-Flag (400 ng), pCAGGS-FIPV N (400 ng), or pCAGGS-ZC3HAV1 (400 ng). After 18 h, cells were stimulated with SeV for 12 h, and luciferase activity was measured using the Dual-Luciferase Reporter Gene Assay Kit (Easen Biotechnology Co., Ltd., Shanghai, China) according to the manufacturer’s instructions.
Antagonistic effect of ZC3HAV1 on FIPV N-mediated IFN-β regulation: 293T cells were seeded in 24-well plates and transfected at 80–90% confluency with pGAS-Luc-IFN-β (human, 60 ng), pRL-TK (10 ng), and one of the following plasmid combinations: pCAGGS-Flag (400 ng); pCAGGS-FIPV N (200 ng, 400 ng, or 800 ng); pCAGGS-ZC3HAV1 (200 ng) + pCAGGS-Flag (200 ng); or pCAGGS-FIPV N (200 ng) + pCAGGS-ZC3HAV1 (200 ng). Each condition was performed in triplicate. After 16 h of transfection, cells were stimulated with SeV for an additional 12 h. Luciferase activity was measured using the Dual-Luciferase Reporter Gene Assay Kit according to the manufacturer’s instructions. For CRFK cells, the same procedure was applied, except that pGAS-Luc-IFN-β (human) was replaced with pGL3-IFN-β (feline) for species-specific validation.
IP and Co-IP verification of the interaction between ZC3HAV1 and FIPV N
The interaction between ZC3HAV1 and the FIPV N protein was validated through IP and Co-IP assays. To exclude the possibility that RNA mediates the interaction, RNase was added to the cell lysate (final concentration: 50 µg/mL) during lysis in the Co-IP experiments.
Confocal microscopy analysis
To examine the subcellular localization and colocalization of ZC3HAV1 and FIPV N, CRFK cells were seeded on glass coverslips and transfected with the following combinations of expression plasmids: pCMV14-Myc-FIPV N + pCMV14-HA (control), pCMV14-HA-ZC3HAV1 + pCMV14-Myc (control), or pCMV14-Myc-FIPV N + pCMV14-HA-ZC3HAV1 (experimental group). At 24 h post-transfection, cells were fixed with 4% paraformaldehyde, permeabilized with 2% Triton X-100, and stained using anti-hemagglutinin (HA) and anti-Myc antibodies followed by AF594- and AF488-conjugated secondary antibodies. Nuclei were counterstained with DAPI. After mounting, cells were visualized by confocal laser scanning microscopy.
In parallel, to confirm the subcellular localization of endogenous proteins under physiological infection conditions, CRFK cells infected with FIPV (MOI = 0.01) were fixed with 4% paraformaldehyde and permeabilized with 2% Triton X-100. Mock-infected cells were included as negative controls. After blocking with 5% BSA, cells were incubated with primary antibodies against endogenous ZC3HAV1 (mouse polyclonal antibody, kindly provided by the Companion Animal Research Group, Shanghai Veterinary Research Institute) and FIPV N protein overnight at 4°C, followed by incubation with Alexa Fluor–conjugated secondary antibodies. Nuclei were stained with DAPI, and images were acquired using a laser-scanning confocal microscope.
Glutathione S-transferase (GST) pull-down assay
To determine whether the interaction between ZC3HAV1 and FIPV N is direct, an in vitro GST pull-down assay was performed. Recombinant ZC3HAV1 and FIPV N proteins were expressed in E. coli Rosetta (DE3) cells using the pET-28a-ZC3HAV1, pEGX-4T-1-FIPV N, and empty pEGX-4T-1 plasmids. To enhance the yield of ZC3HAV1 in the supernatant during prokaryotic expression, cultures were induced with 1 mM IPTG at 16°C for 16 h. Bacterial pellets were harvested by sonication, and the supernatants were collected for SDS-PAGE and WB analysis. For the pull-down assay, GST magnetic beads (Cat#C650031; BBI Life Sciences, Shanghai, China) were incubated with bacterial lysates expressing GST or GST-FIPV N for 1 h, washed twice, and then incubated with His-ZC3HAV1 lysates for an additional 1 h at 25°C. After five washes, the bound proteins were eluted by boiling in SDS loading buffer and analyzed by WB.
Identification of interaction domains between FIPV N protein and ZC3HAV1
Based on the evidence of direct interaction, further experiments were performed to map the binding regions between ZC3HAV1 and FIPV N. In one set of experiments, HEK293T cells were co-transfected with pCMV14-Myc-FIPV N and a series of ZC3HAV1 truncation mutants (Z1–Z6). In a parallel set, cells were co-transfected with pCMV14-HA-ZC3HAV1 and a panel of FIPV N truncation constructs (N1–N5). All constructs were cloned into pCAGGS. GFP-expressing cells (pCAGGS-GFP) were included as negative controls in both assays. At 24 h post-transfection, cell lysates were subjected to immunoprecipitation using anti-Myc or anti-HA conjugated magnetic beads, followed by Western blot analysis.
In addition, protein–protein interaction models between ZC3HAV1 and FIPV N were predicted using the AlphaFold3 online server (https://alphafoldserver.com/). Model confidence was assessed using predicted Local Distance Difference Test (pLDDT) scores. The resulting models were further visualized using PyMOL software to assess the spatial localization of binding domains.
Generation of stable ZC3HAV1-knockdown cells and rescue experiments
CRFK cells with stable ZC3HAV1 knockdown were generated using a lentivirus-mediated shRNA approach. Briefly, pLKO.1-ZC3HAV1 was co-transfected with psPAX2 and pMD2.G into HEK293T cells to produce lentiviral particles. Supernatants containing lentiviral particles were collected and used to transduce CRFK cells. Transduced cells were selected with puromycin (6 μg/mL) (Cat# ST551; Biotanon, Shanghai, China) and maintained under puromycin selection for ten consecutive passages to generate stable knockdown cell lines. shRNA sequences are listed in Supplementary S1 Table. Knockdown efficiency was confirmed by Western blotting and RT-qPCR.
For rescue experiments, CRFK-shZC3HAV1 cells were transfected with plasmids expressing wild-type ZC3HAV1 or the ZC3HAV1 truncation mutant Z5, which lacks the RBD-containing region involved in FIPV N association, with or without co-expression of FIPV N protein. At 16 h post-transfection, cells were infected with FIPV (MOI = 0.01). Samples were collected at the indicated time points, and viral replication was evaluated by Western blotting, RT-qPCR, and TCID50 assays.
To further evaluate whether the ZC3HAV1-associated region of FIPV N contributes to counteracting ZC3HAV1-mediated antiviral activity, CRFK-shZC3HAV1 cells were transfected with wild-type ZC3HAV1 together with full-length FIPV N or the indicated N-protein truncation mutants (N1–N5). Cells transfected with the corresponding empty vectors served as controls. At 16 h post-transfection, cells were infected with FIPV (MOI = 0.01). Samples were collected at the indicated time points. Protein expression was analyzed by Western blotting, while viral RNA levels and infectious virus titers were determined by RT-qPCR using FIPV-specific primers and TCID50 assays, respectively.
Construction and validation of an RNA-binding-deficient FIPV N mutant
To generate an RNA-binding-deficient mutant of the FIPV N protein, three conserved residues (Y93, R117, and K123) within the N-terminal RNA-binding domain were substituted with alanine using overlap PCR with specifically designed mutagenic primers. The resulting mutant (Mut-N; Y93A/R117A/K123A) was cloned into the pCMV-Myc expression vector and confirmed by DNA sequencing. Primer sequences are listed in Supplementary S1 Table.
To validate RNA-binding activity, wild-type N and Mut-N proteins were expressed in HEK293T cells. The interaction between Mut-N and ZC3HAV1 was first examined by Co-IP assays to determine whether the introduced mutations affected protein–protein interaction. RNA-binding capacity was then assessed by RNA immunoprecipitation (RIP) assays as described below. Viral RNA enrichment was quantified by RT-qPCR and normalized to input RNA.
Establishment of qPCR detection methods for FIPV gRNA and sgRNA
To investigate whether ZC3HAV1 regulates FIPV genomic RNA (gRNA) replication and affects sgRNA transcription, qPCR assays were developed based on the coronavirus discontinuous transcription mechanism and the shared 5¢ untranslated region (UTR) leader sequence in gRNA and all sgRNAs. The 5¢ UTR leader sequence of FIPV was first identified by PCR using multiple forward primers spanning the entire 5¢ UTR (primer sequences are listed in S4 Table) and a reverse primer targeting the 5¢ proximal region of the N gene open reading frame (ORF), with cDNA synthesized from FIPV-infected CRFK total RNA as the template. Amplification indicated that the upstream sequences were part of the leader region.
Once the leader sequence was confirmed, specific sgRNA detection primers were designed with the forward primer located within the 5¢ UTR leader sequence and the reverse primer targeting the 5¢ region of the desired ORF, ensuring amplification of only the sgRNA of interest. For gRNA-specific detection, primers were designed within the 3¢ region of the pp1a/pp1b ORF (forward) and the 5¢ region of the S gene ORF (reverse), avoiding amplification of sgRNAs. Detailed primer sequences are provided in S4 Table.
Identification of sgRNA transcription levels during FIPV infection
The synthesis of sgRNA during coronavirus replication follows a temporal order, with sgRNAs from different genes potentially exhibiting distinct expression patterns at various stages of viral replication, including early, middle, and late phases. In subsequent analyses of the impact of ZC3HAV1 on the relative expression of FIPV sgRNA to gRNA, the differences in sgRNA expression at various time points during FIPV replication may have an impact on the results. To explore the transcriptional differences of FIPV sgRNAs, we further investigated sgRNA expression at 12 and 16 hpi in CRFK cells. CRFK cells in good growth condition were seeded in 12-well plates. Upon reaching 100% confluence, FIPV suspension (MOI = 0.01) was added, with an uninfected control group included. Samples were collected at 12 and 16 hpi, and RNA was extracted and reverse-transcribed. qPCR was used to detect FIPV gRNA, and sgRNA was used to assess expression levels.
Quantitative assessment of ZC3HAV1 effects on FIPV subgenomic RNA expression
To investigate whether ZC3HAV1 regulates the transcription of FIPV sgRNAs, CRFK cells were infected with FIPV (MOI = 0.01). Total RNA was collected at 12 and 16 hpi. The expression levels of sgRNAs were analyzed using a qPCR method established to distinguish between gRNA and sgRNAs. To determine the specific regulatory effect of ZC3HAV1, overexpression and knockdown experiments were performed using pCAGGS-ZC3HAV1 and siRNAs, respectively. sgRNA transcription levels were quantified following normalization to both GAPDH and viral gRNA. Specifically, ΔCt values were calculated relative to GAPDH. The sgRNA levels were normalized to gRNA using the ΔΔCt method: ΔΔCtsgRNA/gRNA = (ΔCtsgRNA − ΔCtgRNA)treatment−(ΔCtsgRNA − ΔCtgRNA)control. The relative expression of sgRNA was calculated as 2−ΔΔCt sgRNA/gRNA, allowing assessment of ZC3HAV1-mediated regulation independent of overall viral RNA replication.
GS-441524-based assay for FIPV RNA persistence
To determine whether ZC3HAV1 regulates the abundance of FIPV RNAs independently of ongoing viral RNA synthesis, the RNA-dependent RNA polymerase inhibitor GS-441524 was used to block de novo viral RNA synthesis. CRFK cells were infected with FIPV (MOI = 0.01), and GS-441524 was added at a final concentration of 200 μM at 4 hpi. An untreated infected group served as a positive control. Cells were harvested at 0, 4, 8, and 12 h post-treatment, and total RNA was extracted for RT-qPCR-based quantification of viral genomic RNA (gRNA) and subgenomic RNAs (sgRNAs), which was used to determine the optimal time window for subsequent experiments.To investigate the role of ZC3HAV1, both gain- and loss-of-function approaches were performed in parallel. For overexpression, CRFK cells were transfected with pCMV14-HA-ZC3HAV1 or empty vector 12 h prior to infection. For knockdown, cells were transfected with siZC3HAV1 or negative control siRNA (siNC) 60 h prior to infection. Cells were then infected with FIPV (MOI = 0.01), and GS-441524 (200 μM) was added at 4 hpi to inhibit de novo viral RNA synthesis. Cells were harvested 10 h after drug addition for total RNA extraction and RT-qPCR-based quantification of viral gRNA and sgRNAs.
Analysis of CpG distribution in FIPV and comparative genomic analysis with other coronaviruses
To investigate whether ZC3HAV1 preferentially targets FIPV RNA via CpG recognition, we analyzed CpG distribution in the FIPV genome and performed a comparative analysis with other coronaviruses. Full-length genomic sequences of FIPV and five representative coronaviruses—HCoV-229E (GenBank: KU291448.1), CCoV (KP981644.1), SARS-CoV-2 (MZ054889.1), IBV (OR824987.1), and PDCoV (MN942260.2)—were retrieved from GenBank. Using R, CpG frequency was calculated in 100 bp non-overlapping windows across genomes, and visualized to assess interspecies differences. For FIPV, CpG observed/expected (O/E) ratios were further calculated in defined genomic regions (e.g., 5′ UTR, ORF1a/1ab, S, N) using sliding window analysis. Regional distributions were visualized with violin-box plots, and the 5′ UTR was analyzed in 10 bp windows to capture CpG variation near the leader sequence.
RNA IP (RIP) assay of ZC3HAV1 binding to genomic and subgenomic RNAs of FIPV
To verify whether ZC3HAV1 targets coronavirus gRNA and sgRNA, RIP was performed using the BeyoRIP RIP Assay Kit (P0107; Beyotime) following the manufacturer’s instructions. CRFK cells were seeded in 6-well plates and transfected at 90% confluence with 2 μg of pCAGGS-ZC3HAV1, pCAGGS-ZC3HAV1-Z2 (RBD domain), pCAGGS-ZC3HAV1-Z4 (WWE+PARP-like domains), pCMV14-Flag-FIPV N, or pCAGGS-Flag, with non-transfected cells as controls. After 20 h, cells were infected with FIPV (MOI = 0.01) for 16 h. Lysates were prepared in RIP lysis buffer supplemented with dithiothreitol (DTT), RNase inhibitor, and phenylmethanesulfonyl fluoride (PMSF), followed by IP using anti-Flag antibody (F1804; Sigma-Aldrich) or mouse IgG isotype control (M07001; Abmart, Berkeley Heights, NJ, USA) prebound to Protein A/G agarose beads (P2012; Beyotime). RNA was eluted at 55°C and purified using RNAeasy Animal RNA Extraction Kit (R0027; Beyotime), with 20 μL reserved for input validation by WB. Co-precipitated RNA was reverse transcribed and analyzed by qPCR to detect FIPV N sgRNA. The relative enrichment of target RNAs in RIP samples was calculated using the %Input and Fold Change methods: ΔCtRIP = CtRIP−(CtInput−log2(Input dilution factor)); %Input = 2^−ΔCtRIP; ΔΔCtRIP/Flag = CtRIP − CtFlag; Fold Change = 2^−ΔΔCtRIP/Flag.
Competitive RIP assay to evaluate interference of FIPV N protein with ZC3HAV1–viral RNA interaction
To investigate whether FIPV N protein can antagonize the interaction between ZC3HAV1 and viral sgRNA, a competitive RIP assay was performed in vitro. Since N protein binds to viral RNA and ZC3HAV1 directly interacts with N protein, the experimental procedures were optimized to avoid interference from N protein binding to FIPV RNA. Briefly, CRFK cells were transfected with pCAGGS-ZC3HAV1, pCMV14-Myc, and pCMV14-Myc-N plasmids and incubated for 24 h before harvesting. The cells were lysed using lysis buffer containing DTT, RNase inhibitor, and PMSF, and the cell lysates were collected. For RNA preparation, CRFK cells were infected with FIPV (MOI = 0.01) for 16 h, and total RNA was extracted using the RNAeasy Animal RNA Extraction Kit. Protein A/G Agarose was pre-incubated with HA antibody to capture ZC3HAV1. The protein complexes were then incubated with cell lysates expressing ZC3HAV1 for 1 h at 25°C, followed by washing. Subsequently, Myc-N protein was added to the experimental group and Myc protein to the control group, and the mixtures were incubated for an additional 1 h. After washing, FIPV RNA (5 µg) was added to both groups, and the reactions were incubated for 4 h at 4°C. RIP RNA was analyzed using the ΔΔCt method to calculate the enrichment fold change of FIPV genomic RNA (gRNA) and sgRNA in both the experimental and control groups. This experimental setup allowed for the assessment of whether FIPV N protein competes with ZC3HAV1 for binding to viral RNA.
To determine whether FIPV N protein could displace viral RNA from preformed ZC3HAV1–RNA complexes, HA-ZC3HAV1 protein complexes were immunoprecipitated from lysates of transfected CRFK cells using anti-HA antibody pre-bound to Protein A/G agarose beads. Cells transfected with empty HA vector (HA control) were included as negative controls. A total of 5 µg RNA, extracted from CRFK cells infected with FIPV (MOI = 0.01) for 16 h, was incubated with bead-bound HA-ZC3HAV1 protein for 4 h at 4°C to allow formation of ZC3HAV1–RNA complexes. After washing, equivalent amounts of Myc-N, Myc-Mut-N, or Myc proteins were added. Samples were incubated for an additional 1 h at 25°C. After extensive washing, RNA was extracted from the immunoprecipitated beads using the RNAeasy Animal RNA Extraction Kit. Viral RNAs were quantified by RT-qPCR, and the abundance of FIPV gRNA and sgRNAs was normalized to input RNA.
To evaluate whether FIPV N protein can disrupt ZC3HAV1 binding to host mRNA, HSPG2 mRNA, a previously reported ZC3HAV1-associated host transcript, was used as a representative host target. The same procedure was performed, except that the RNA used for complex formation was 5 µg total RNA extracted from uninfected CRFK cells. All subsequent steps were performed as described above. HSPG2 mRNA levels were quantified by RT-qPCR.
Analysis of the antagonistic effect of FIPV N on ZC3HAV1-mediated viral suppression
To assess whether FIPV N protein can antagonize the antiviral effect of ZC3HAV1, we used JEV as a model. ZC3HAV1 has been shown to target the 3¢ UTR of JEV RNA to inhibit replication. BHK-21 cells were transfected with pCAGGS-Flag, pCMV14-Flag-FIPV N, pCAGGS-ZC3HAV1, or combinations of these plasmids. After 12 h, cells were infected with JEV (MOI = 0.1) and collected at 12, 24, and 36 hpi. WB analysis was performed to measure JEV NS1 protein expression, indicating viral replication levels. This setup allowed an examination of whether FIPV N protein can counteract the antiviral activity of ZC3HAV1.
AlphaFold3 structural prediction
The interaction model of FIPV N and ZC3HAV1 was predicted using the AlphaFold3 online server. The confidence of the predicted structure was evaluated using the predicted local distance difference test (pLDDT) score provided by AlphaFold3, with pLDDT > 90 indicating very high confidence, 70–90 indicating confident prediction, 50–70 indicating low confidence, and < 50 indicating very low confidence. The pLDDT score was visualized by color coding in the predicted model. The resulting structures were visualized using PyMOL. Because the predicted complex was not validated by residue-level mutagenesis or structural experiments, the AlphaFold3 model was used to guide interpretation of the domain-mapping results rather than to define experimentally confirmed amino acid contacts.
Statistical analysis
All experimental data were statistically analyzed using GraphPad Prism version 9 (GraphPad Software, San Diego, CA, USA). Differences between groups were evaluated using multiple t-tests, with p < 0.05 considered statistically significant.
Results
Proteomic profiling identifies ZC3HAV1 as a key N-associated host restriction factor
To systematically map host factors interacting with the FIPV nucleocapsid (N) protein, we employed two complementary co-immunoprecipitation (Co-IP) strategies coupled with mass spectrometry (MS) (Fig 1A). In the first approach, endogenous N-associated host proteins were co-precipitated from FIPV-infected CRFK cells using the anti-N monoclonal antibody 10E7. In parallel, ectopically expressed Flag-tagged N was immunoprecipitated from transfected cells with an anti-Flag antibody to capture host interactors under controlled expression conditions. Western blot analysis confirmed that both antibodies efficiently and specifically pulled down N along with associated host proteins, with minimal nonspecific background in the controls (Fig 1B), validating the reliability of both pipelines for proteomic screening of N-interacting host factors.
[图:Identification and characterization of host proteins associated with the FIPV N protein. (A) Schematic workflow of the c]

Label-free LC–MS analysis identified 2,197, 1,975, 3,589, and 3,630 protein entries in the CRFK, CRFK + FIPV, pCMV14-FIPV N, and pCMV14-GFP groups, respectively (S5 Table). Notably, the number of protein entries detected in FIPV-infected cells was approximately 10% lower than that in uninfected cells. This reduction may be attributed to multiple factors, including virus-induced alterations in host protein abundance, solubility, or subcellular localization, competition or epitope-masking effects during immunoprecipitation, and infection-associated changes in sample complexity or cellular stress. To define candidate N-associated host proteins, we calculated the enrichment ratio for each protein relative to its corresponding control. Proteins with an enrichment ratio >2 in the FIPV-infected anti-N IP sample relative to the uninfected CRFK control IP sample, or in the Flag-FIPV N IP sample relative to the Flag-GFP control IP sample, were retained as candidate interactors. Proteins detected exclusively in the control samples or lacking enrichment over the corresponding background control were excluded as nonspecific background proteins. The filtered protein lists from the infection- and transfection-based datasets were then integrated and manually curated to remove redundant entries, resulting in a final set of 98 high-confidence N-associated host proteins for subsequent functional analyses.
To explore the biological context of host factors interacting with the FIPV N protein, we constructed a protein–protein interaction (PPI) network using the STRING database (Fig 1C). The network revealed a highly interconnected architecture dominated by explicitly labeled ribosomal proteins of the RPL and RPS families, together with representative translation-related factors such as EIF3D and EIF3M, indicating that N interacts extensively with the host translational machinery. Several RNA-associated factors, including ZC3HAV1, MOV10, HNRNPF, DDX18, and RALY, also occupied central nodes within the network, suggesting potential links between viral RNA handling and host RNA-binding processes.
Gene Ontology (GO) enrichment analysis further characterized the molecular features of these interactors across three functional dimensions. In the cellular component category, proteins were mainly enriched in the cytosolic ribosome and ribonucleoprotein complexes (Fig 1D), indicating that many N-associated host factors are integral parts of translation assemblies. In the molecular function category, enrichment was dominated by structural constituents of ribosomes, RNA-binding, and nucleic acid-binding activities (Fig 1E), reflecting a strong association with RNA metabolism and translational regulation. In the biological process category, significant enrichment was observed for translation, peptide and amide metabolic processes, cytoplasmic translation, and ribosome biogenesis (Fig 1F), revealing that N-interacting proteins converge functionally on protein synthesis and RNA processing pathways.
Together, these integrated proteomic and bioinformatic analyses demonstrate that the FIPV N protein engages a dense network of host factors involved in translation and RNA metabolism—two pathways that often overlap with intrinsic post-transcriptional restriction and antiviral RNA surveillance. This suggests that N may occupy a pivotal functional interface between viral RNA synthesis and host RNA-regulatory systems, potentially enabling the virus to couple efficient genome replication with evasion of cellular RNA-level restriction. From the 98 N-associated host proteins identified by MS analysis, we next prioritized 18 candidates based on their annotated functions in RNA metabolism, translation regulation, or antiviral defense S6 Table. Transcriptional profiling during FIPV infection revealed a characteristic biphasic pattern—most genes were transiently upregulated at 12 h post-infection and subsequently downregulated by 24 h (Fig 1G). This temporal pattern may represent an initial host effort to stabilize RNA metabolism and restrict viral replication, which is subsequently overridden as the virus gains control over the host translational machinery.
To evaluate their functional relevance, each of the 18 candidates was cloned into pCAGGS for overexpression in CRFK cells, followed by FIPV challenge and viral titer quantification. Six RNA-associated proteins—DDX18, MOV10, RPL36, ZC3HAV1, HNRNPF, and RALY—significantly reduced viral replication (Fig 1H), identifying them as putative restriction factors within the N interactome. Notably, several of these factors, such as MOV10 and ZC3HAV1, have been implicated in post-transcriptional antiviral restriction in other RNA viruses [19,20], reinforcing the possibility that FIPV engages components of the host RNA surveillance network.
Among these, ZC3HAV1 stood out for its well-established role in RNA-targeted restriction across diverse viral families, including coronaviruses [18,21,22]. Its dual identification as an N-associated host protein and as a negative regulator of FIPV replication prompted us to investigate how ZC3HAV1 restricts FIPV and how FIPV N counteracts this restriction.
ZC3HAV1 restricts FIPV replication through an IFN-independent RNA-level mechanism
Having identified ZC3HAV1 as a potent N-binding restriction factor, we next sought to validate its antiviral activity and elucidate its mechanism of action. To determine whether the host restriction factor ZC3HAV1 affects FIPV replication, CRFK cells were transfected with a ZC3HAV1 expression plasmid and subsequently infected with FIPV (MOI = 0.01). Samples were collected at indicated time points for Western blotting, TCID₅₀ assays, and absolute RT–qPCR quantification. Western blot analysis showed that ZC3HAV1 overexpression markedly suppressed FIPV replication, as reflected by the decreased N protein abundance at 12, 16, and 24 hpi (Fig 2A–B). A dose-dependent reduction in N expression was also observed upon increasing amounts of ZC3HAV1 plasmid (Fig 2C–D). Absolute RT–qPCR confirmed a significant decrease in viral N transcripts from 8 to 32 h (Fig 2E), while TCID₅₀ assays further demonstrated a corresponding reduction in viral titers (Fig 2F). Together, these data indicate that ZC3HAV1 acts as a potent suppressor of FIPV replication.
[图:ZC3HAV1 suppresses FIPV replication through an IFN-independent antiviral mechanism. (A) Western blot analysis of FIPV N ]

Endogenous ZC3HAV1 mRNA expression increased following FIPV infection, with the highest induction observed at 12 hpi and sustained elevation still detectable at 24 hpi (Fig 2G). These data indicate that ZC3HAV1 is transcriptionally responsive to FIPV infection. Because ectopically expressed ZC3HAV1 showed reduced protein abundance at later infection time points (Fig 2A), we further examined whether this decrease could be attributed to direct degradation mediated by FIPV N. Co-transfection of HA-ZC3HAV1 with increasing amounts of Flag-FIPV N did not noticeably reduce HA-ZC3HAV1 protein levels, and densitometric analysis showed no significant decrease in HA-ZC3HAV1 abundance (S1A-B Fig). These results do not support a model in which FIPV N alone directly induces ZC3HAV1 degradation.
We also examined endogenous ZC3HAV1 protein levels during FIPV infection. In contrast to the late decrease observed for ectopically expressed ZC3HAV1, endogenous ZC3HAV1 protein was increased in FIPV-infected cells compared with mock-infected cells at 12, 16, and 24 hpi, accompanied by increasing FIPV N expression as a marker of productive infection (S1C-D Fig). Together with the ZC3HAV1 mRNA induction shown in Fig 2G, these data indicate that FIPV infection induces endogenous ZC3HAV1 expression rather than suppressing it. Therefore, the late decrease in ectopically expressed ZC3HAV1 likely reflects late post-transfection protein turnover or infection-associated changes in cellular protein homeostasis, rather than direct degradation by FIPV N.
To further validate its inhibitory role, three siRNAs were designed to knock down ZC3HAV1 expression. Among them, siZC3HAV1–1 achieved the highest silencing efficiency at both mRNA and protein levels (Fig 2H–J). Knockdown of ZC3HAV1 significantly enhanced FIPV replication, as evidenced by increased N protein expression (Fig 2K–L), elevated viral RNA levels (Fig 2M), and higher viral titers (Fig 2N). These complementary overexpression and knockdown results consistently support ZC3HAV1 as a negative regulator of FIPV replication.
Because ZC3HAV1 is an interferon-stimulated gene (ISG), its antiviral activity could either arise indirectly through type I interferon (IFN) pathway activation or directly through intrinsic RNA-targeted restriction. To distinguish between these possibilities, we next examined whether the antiviral effect of ZC3HAV1 depends on IFN signaling activation or instead reflects an intrinsic post-transcriptional restriction mechanism. Dual-luciferase reporter assays revealed that both ZC3HAV1 and the FIPV N protein independently suppressed Sendai virus (SeV)–induced IFN-β promoter activity in 293T cells (Fig 2O). However, co-expression experiments in both 293T and CRFK cells showed no synergistic or antagonistic interaction between the two proteins in regulating IFN-β transcription (Fig 2P–Q). These results indicate that, under the conditions of this reporter assay, ZC3HAV1-mediated suppression of FIPV replication is unlikely to be explained by type I interferon (IFN) pathway activation. We therefore next examined whether ZC3HAV1 restricts FIPV through direct association with viral RNA.
ZC3HAV1 reduces the accumulation and persistence of FIPV genomic and subgenomic RNAs
Given its potent antiviral activity, we next investigated whether ZC3HAV1 restricts FIPV replication at the RNA level through direct interference with viral transcript accumulation. To this end, we first established a quantitative detection system that distinguishes between genomic RNA (gRNA) and individual subgenomic RNAs (sgRNAs) produced through the coronavirus discontinuous transcription mechanism. As illustrated in Fig 3A, a series of forward primers (Lp1–Lp12) within the 5’ UTR were paired with a reverse primer located at the 5’ end of the N gene. PCR amplification using viral cDNA as the template yielded distinct products with Lp1–Lp5 (Fig 3B), indicating that the leader sequence of FIPV resides within this region. Sequence alignment with representative coronaviruses revealed that the FIPV transcription-regulatory sequence in the leader (TRS-L) is ACGAAA, and the complete leader sequence was determined (Fig 3C).
[图:ZC3HAV1 reduces the accumulation and persistence of FIPV genomic and subgenomic RNAs. (A) Schematic of primer design str]

Based on this mapping, we designed specific primer sets for each sgRNA (S, 3a, 3c, E, M, N, 7a, and 7b) and one for gRNA detection. Each pair amplified a single, expected band, confirming primer specificity and establishing a robust assay for quantifying FIPV gRNA and sgRNAs.
We next analyzed the relative accumulation pattern of FIPV sgRNAs by normalizing each sgRNA to gRNA levels. Because gRNA reflects overall viral genome abundance and serves as the template for sgRNA synthesis, the sgRNA/gRNA ratio provides a normalized readout of sgRNA accumulation after accounting for differences in total viral RNA levels. Among the detected sgRNAs, the N sgRNA showed the highest sgRNA/gRNA ratio, followed by M and S, whereas several accessory sgRNAs displayed lower relative abundance (Fig 3D). This pattern is consistent with canonical coronavirus discontinuous transcription, in which structural gene sgRNAs generally accumulate more prominently than accessory transcripts. The sgRNA/gRNA ratios of most transcripts remained broadly comparable between 12 h and 16 h post-infection, with only modest transcript-specific fluctuations (Fig 3D). Thus, Fig 3D establishes a reference profile of sgRNA-to-gRNA relationships during FIPV infection, allowing us to assess whether ZC3HAV1 alters overall viral RNA accumulation or selectively reshapes the relative sgRNA profile.
To assess whether ZC3HAV1 influences these viral RNAs, CRFK cells were transfected with ZC3HAV1 expression plasmids or siRNAs prior to infection. Overexpression of ZC3HAV1 markedly reduced FIPV gRNA levels at 12 h and 16 h (Fig 3E), whereas its knockdown had the opposite effect (Fig 3F). A similar pattern was observed for sgRNAs: ZC3HAV1 overexpression significantly decreased, and knockdown increased, the transcript abundance of multiple structural and accessory sgRNAs (Fig 3G–H). These data indicate that ZC3HAV1 reduces the steady-state accumulation of both genomic and subgenomic FIPV RNAs.
Because gRNA serves as the template for sgRNA synthesis, we further normalized sgRNA abundance to gRNA levels to evaluate whether ZC3HAV1 altered the relative sgRNA profile independently of changes in overall viral genome abundance. Although ZC3HAV1 overexpression reduced, and ZC3HAV1 knockdown increased, the absolute levels of both gRNA and multiple sgRNAs (Fig 3E–H), the sgRNA/gRNA ratios showed only limited changes for most sgRNA species (Fig 3I–J). These results suggest that ZC3HAV1 primarily reduces overall FIPV RNA accumulation rather than markedly reshaping the relative abundance of individual sgRNAs. Thus, the decrease in sgRNA levels is largely consistent with a broad reduction in viral RNA abundance, rather than a selective alteration in the relative sgRNA composition.
However, steady-state RT-qPCR analysis alone cannot determine whether reduced viral RNA abundance results from impaired viral RNA synthesis or from reduced persistence of already synthesized viral RNAs. To determine whether ZC3HAV1 affects the persistence of already synthesized viral RNAs, we performed a GS-441524-based RNA persistence assay after FIPV RNAs had accumulated. GS-441524 treatment effectively suppressed the continued accumulation of FIPV N RNA and reduced normalized sgRNA/gRNA expression over time, confirming that ongoing RdRp-dependent RNA synthesis was strongly inhibited under these conditions (S2A–B Fig). In this replication-suppressed setting, ZC3HAV1 overexpression reduced the levels of FIPV gRNA and multiple sgRNAs, whereas ZC3HAV1 knockdown increased their abundance (S2C–E Fig). After normalization to gRNA, ZC3HAV1 overexpression also decreased, whereas ZC3HAV1 knockdown increased, the relative abundance of several sgRNAs (S2F–G Fig). These results indicate that ZC3HAV1-mediated restriction is not explained solely by reduced viral RNA synthesis or altered sgRNA transcription, but involves reduced persistence of already synthesized FIPV genomic and subgenomic RNAs.
Interestingly, the suppressive effect of ZC3HAV1 on viral RNA accumulation weakened at later infection time points, suggesting that FIPV may counteract ZC3HAV1-associated RNA-level restriction as infection progresses. Together, these findings identify ZC3HAV1 as a restriction factor that limits FIPV genomic and subgenomic RNA accumulation, at least in part, by reducing viral RNA persistence after RNA synthesis has occurred. Although these data are consistent with destabilization or degradation of ZC3HAV1-associated viral RNAs, they do not define the precise downstream RNA-decay machinery involved. Together, these findings identify ZC3HAV1 as a restriction factor that limits FIPV genomic and subgenomic RNA accumulation primarily by reducing viral RNA persistence after RNA synthesis has occurred. Because our proteomic analysis identified FIPV N as a major ZC3HAV1-associated viral protein, this RNA-level restriction phenotype provided the rationale for examining whether FIPV N counteracts ZC3HAV1-mediated control of viral RNA persistence.
FIPV N engages ZC3HAV1 and counteracts ZC3HAV1-mediated antiviral restriction
To test our hypothesis that FIPV N interferes with ZC3HAV1-associated RNA-level restriction, we next examined whether N interacts with ZC3HAV1. In CRFK cells infected with FIPV (MOI = 0.01), co-immunoprecipitation (Co-IP) assays revealed that Flag-tagged ZC3HAV1 efficiently pulled down the endogenous N protein, whereas no interaction was detected in IgG or vector controls (Fig 4A). This interaction was further confirmed in 293T cells co-expressing Myc-tagged N and Flag-ZC3HAV1, indicating that their association is not limited to the CRFK cell background (Fig 4B). Importantly, treatment with RNase did not abolish the interaction, suggesting that the association is not primarily mediated by RNA bridging (Fig 4B, right).
[图:ZC3HAV1 directly interacts with the FIPV N protein via its RBD-containing domain. (A) Co-immunoprecipitation analysis of]

To further examine whether ZC3HAV1 and FIPV N can associate in vitro, His-ZC3HAV1 was incubated with clarified soluble bacterial lysates containing GST-N or GST control protein. GST pull-down followed by immunoblotting showed that His-ZC3HAV1 was recovered with GST-N but not with GST alone, supporting an in vitro association between ZC3HAV1 and FIPV N under the conditions of this assay (Fig 4C). Confocal immunofluorescence analysis showed partial cytoplasmic overlap between overexpressed ZC3HAV1 and FIPV N (Fig 4E), consistent with the possibility that the interaction occurs in the cytoplasm, where viral RNA replication and RNA-level restriction take place. To further evaluate whether this spatial association could also be observed under infection conditions, we examined endogenous ZC3HAV1 and infection-derived FIPV N in CRFK cells. In uninfected CRFK cells, endogenous ZC3HAV1 showed cytoplasmic staining, whereas FIPV N staining was not detected. Upon FIPV infection, N protein accumulated mainly in the cytoplasm and showed partial overlap with endogenous ZC3HAV1-positive signals (S3 Fig). These data support the infection-context relevance of the ZC3HAV1–N spatial association observed in the overexpression system.
To map the regions required for this interaction, we performed reciprocal domain mapping using truncation mutants of both proteins. For FIPV N, six constructs were generated (N1–N6) encompassing different combinations of its N-terminal domain (NTD) and C-terminal domain (CTD) (Fig 4D). To exclude nonspecific recovery of overexpressed proteins, GFP was included as an unrelated protein negative control in the Co-IP assay. Co-IP assays revealed that ZC3HAV1 interacted efficiently with full-length N and the N4 fragment containing both NTD and CTD, whereas GFP and N fragments lacking either region showed weak or undetectable recovery (Fig 4F). Thus, efficient association with ZC3HAV1 was observed only when both the NTD and CTD were present. These data suggest that the structural integrity of both the NTD and CTD is required for efficient interaction with ZC3HAV1. Because deletion of either domain may also affect N folding, oligomerization, or RNA-binding properties, these truncation results were interpreted as supportive domain-mapping evidence rather than definitive identification of a minimal binding interface.
Similarly, ZC3HAV1 truncation constructs were designed to refine the N-interacting region. The RBD of ZC3HAV1 contains five CCCH-type zinc finger motifs, designated ZnF1–ZnF5, whereas the ZC3HAV1 truncation constructs are designated Z1–Z6 (Fig 4G). GFP was included as an unrelated protein negative control to evaluate nonspecific recovery in the Co-IP assay. FIPV N efficiently co-immunoprecipitated full-length ZC3HAV1 and the RBD-containing construct encompassing ZnF1–ZnF5, whereas the other truncation constructs were not detectably recovered (Fig 4H). These results indicate that the RBD-containing region of ZC3HAV1, which encompasses ZnF1–ZnF5, is required for efficient association with FIPV N.
To further visualize the putative molecular interface, an AlphaFold3-predicted structural model was generated and analyzed in PyMOL (Fig 4I–J). The model suggested that the NTD and CTD of the FIPV N protein may form a putative composite surface positioned near the N-terminal region of ZC3HAV1 encompassing the RBD (residues 227–634). The linker L1 and the fifth CCCH-type zinc finger motif were positioned near the predicted contact region, providing a structural hypothesis consistent with the domain-mapping results. Because this model is computationally predicted, it was used to guide interpretation of the domain-mapping data rather than as definitive structural evidence.
Together, these findings establish that FIPV N associates with ZC3HAV1 and that efficient interaction requires the RBD-containing region of ZC3HAV1 and the intact structural organization of FIPV N. We next investigated whether this interaction contributes to the ability of FIPV N to counteract ZC3HAV1-mediated antiviral restriction.
To further assess the functional relevance of the N–ZC3HAV1 interaction, we performed rescue experiments in ZC3HAV1-depleted CRFK-shZC3HAV1 cells. Full-length HA-ZC3HAV1 or HA-Z5, a truncation lacking the RBD-containing region involved in FIPV N interaction, was reconstituted in these cells, and Myc-FIPV N was co-expressed with full-length HA-ZC3HAV1 where indicated, followed by FIPV infection. Reconstitution with full-length ZC3HAV1 restored antiviral restriction, as indicated by reduced FIPV-M protein expression at 12 and 24 hpi and decreased viral titers over time (S4A–C Fig). Co-expression of FIPV N partially reversed the inhibitory effect of full-length ZC3HAV1, resulting in increased FIPV M protein abundance and higher viral titers compared with ZC3HAV1 reconstitution alone (S4A–C Fig). These results support a functional antagonistic effect of FIPV N on ZC3HAV1-mediated restriction.
In contrast, HA-Z5 showed markedly impaired antiviral activity compared with full-length ZC3HAV1, as reflected by higher FIPV M protein levels and increased viral titers (S4A–C Fig). This finding indicates that the RBD-containing region is required not only for FIPV N association but also for the intrinsic antiviral function of ZC3HAV1. Therefore, HA-Z5 was interpreted as a functionally impaired truncation rather than a strict interaction-deficient but restriction-competent mutant.
To further determine whether the ZC3HAV1-associated region of FIPV N contributes to its antagonistic activity, we performed additional functional analyses using the N truncation mutants identified in the interaction mapping assay (Fig 4D and F). Full-length N and individual N truncation mutants were evaluated in the ZC3HAV1 rescue system. Expression analysis confirmed that the N fragments were detectable under the experimental conditions, allowing functional comparison among different constructs (S4D Fig). Compared with full-length N, the antagonistic activity of individual N truncation mutants varied according to their retained structural regions. Among the truncation mutants, the N4 fragment, which showed efficient ZC3HAV1 association in the mapping assay, displayed relatively stronger capacity to relieve ZC3HAV1-mediated restriction, whereas fragments with weak or undetectable ZC3HAV1 association showed reduced antagonistic activity, as measured by viral protein accumulation and viral RNA levels (S4E and F Fig.). Consistently, analysis of infectious virus production showed a similar pattern among N truncation mutants (S4G Fig). These results suggest that the ZC3HAV1-associated region of FIPV N contributes to efficient antagonism of ZC3HAV1-mediated antiviral restriction.
Together, these rescue analyses support a model in which the N–ZC3HAV1 interaction contributes to FIPV N-mediated antagonism of host RNA-level restriction, while the mechanism by which N disrupts ZC3HAV1–viral RNA association was further investigated below.
ZC3HAV1 targets FIPV genomic and subgenomic RNAs through its RBD-containing region despite low CpG frequency
ZC3HAV1 has been classically defined as a CpG-biased antiviral factor that preferentially associates with CpG-enriched viral RNAs and restricts their accumulation. To determine whether this canonical CpG-biased recognition model applies to FIPV, we next examined genome-wide CpG distribution patterns across FIPV and related coronaviruses. Genome-wide profiling revealed that the overall CpG frequency of FIPV was markedly low, with enrichment confined to the 5′ region (Fig 5A–B). In contrast, porcine coronaviruses such as PDCoV and PEDV displayed higher and more uniform CpG densities, while HCoV-229E and FIPV shared a similar pattern characterized by a single 5′-end CpG cluster and sparse distribution elsewhere.
[图:ZC3HAV1 binds FIPV genomic and subgenomic RNAs through its RBD-containing region despite the overall low CpG frequency o]

A more refined window analysis of the FIPV genome (Fig 5C) confirmed a localized CpG enrichment near the 5′ untranslated region (UTR), with the remainder of the genome showing substantial CpG depletion. Because all coronavirus sgRNAs share a common 5′ leader derived from this UTR, we hypothesized that ZC3HAV1 might recognize sgRNAs through CpG-rich leader sequences. However, mapping the CpG composition within the FIPV 5′UTR (Fig 5D) revealed that CpG enrichment occurs outside the canonical leader–body junction, suggesting that ZC3HAV1–RNA recognition in FIPV is not fully explained by local CpG enrichment within the leader-containing region.
To experimentally verify whether ZC3HAV1 binds FIPV RNA, we performed RNA immunoprecipitation (RIP) assays in FIPV-infected CRFK cells overexpressing different ZC3HAV1 domain constructs—full-length ZC3HAV1, the isolated RNA-binding domain (RBD), or the WWE + PARP region—along with FIPV N and Flag controls. Western blotting confirmed efficient precipitation of each Flag-tagged construct (Fig 5E). Subsequent RT-qPCR analysis of immunoprecipitated RNA revealed that full-length ZC3HAV1 and its RBD markedly enriched multiple FIPV RNAs, including gRNA and sgRNAs encoding structural (S, M, N, E) and accessory proteins (3a, 3c, 7a, 7b) (Fig 5F). The enrichment signal varied among individual sgRNAs, with some accessory sgRNAs, such as 3c, showing relatively low recovery. This variation may reflect differences in intracellular transcript abundance, reduced assay sensitivity for low-copy RNA species, or weaker association with ZC3HAV1-containing ribonucleoprotein complexes. In contrast, the WWE + PARP fragment showed minimal enrichment, comparable to the IgG and Flag controls. These RIP data support ZC3HAV1 association with FIPV RNAs in infected cells, but do not by themselves define the precise CpG motif, RNA structure, or downstream fate of ZC3HAV1-bound viral RNAs.
Consistent with previous reports that coronavirus N proteins can associate with both genomic and subgenomic viral RNAs in infected cells, FLAG-FIPV N also enriched multiple FIPV RNA species in the RIP assay, including gRNA and several sgRNAs (Fig 5F). Because the RIP assay was performed using infected cell lysates, these signals reflect intracellular N-associated viral RNAs during infection. The relatively strong sgRNA signals are consistent with the high intracellular abundance of coronavirus sgRNAs generated during discontinuous transcription.
These findings show that ZC3HAV1 associates with FIPV genomic and subgenomic RNAs through its RBD-containing region despite the overall low CpG frequency of the FIPV genome. Given the low CpG density across most of the FIPV genome, ZC3HAV1 recognition may involve additional RNA features, such as local RNA structure, subgenomic RNA abundance, local sequence context, localized CpG-enriched regions, or viral ribonucleoprotein organization, rather than being explained solely by global CpG frequency. Thus, our data support a noncanonical or CpG-limited mode of ZC3HAV1 engagement, but do not establish complete CpG independence.
FIPV N disrupts ZC3HAV1–viral RNA association to attenuate RNA-level restriction
Having established that FIPV N interacts with ZC3HAV1 and that ZC3HAV1 associates with FIPV genomic and subgenomic RNAs, we next investigated whether N interferes with ZC3HAV1–viral RNA association. We first performed a competitive RIP assay in which HA-ZC3HAV1 was pre-bound to Protein A/G–coupled anti-HA beads and then incubated with Myc-tagged FIPV N or Myc control before RNA capture (Fig 6A). Under these conditions, N readily associated with bead-bound HA-ZC3HAV1 (Fig 6B), but its presence significantly reduced the recovery of multiple FIPV transcripts in the ZC3HAV1 RIP, including gRNA and several sgRNAs (Fig 6C). These results indicate that pre-associated FIPV N can impair the ability of ZC3HAV1 to capture FIPV genomic and subgenomic RNAs. However, because this assay mainly tested whether N bound to ZC3HAV1 before RNA exposure limits subsequent RNA recovery, it did not determine whether N can disrupt pre-formed ZC3HAV1–viral RNA complexes.
[图:FIPV N disrupts ZC3HAV1–viral RNA association to attenuate ZC3HAV1-mediated RNA-level restriction. (A) Schematic diagram]

To address this issue, we next performed a displacement RIP assay using WT FIPV N and an RNA-binding-impaired N mutant. Mut-N was generated by substituting three conserved RNA-binding residues in FIPV N, Y93, R117, and K123, with alanine (Y93A/R117A/K123A) (S5A Fig.). Under RNase-treated conditions, Mut-N retained detectable interaction with HA-ZC3HAV1, comparable to WT N, indicating that these mutations did not abolish N–ZC3HAV1 association (S5B Fig). However, RIP-qPCR analysis showed that Mut-N displayed markedly reduced enrichment of FIPV gRNA and multiple sgRNAs compared with WT N, confirming that the mutant was impaired in viral RNA association while retaining ZC3HAV1-binding capacity (S5C Fig).
We then designed a displacement RIP assay to test whether FIPV N can reduce viral RNA association after ZC3HAV1–RNA complexes have already formed. In this assay, bead-bound HA-ZC3HAV1 was first incubated with FIPV RNA or host RNA to allow RNA binding, unbound RNA was removed by washing, and Myc-tagged WT N or Mut-N was subsequently added before RIP-qPCR analysis (S5D Fig). WT N markedly reduced the levels of ZC3HAV1-associated FIPV gRNA and multiple sgRNAs, whereas Mut-N showed substantially weaker displacement activity (S5E Fig). These results demonstrate that FIPV N can reduce ZC3HAV1-associated viral RNA after ZC3HAV1–viral RNA complexes have formed, and that efficient disruption requires the viral RNA-binding capacity of N.
To determine whether this effect reflected nonspecific disruption of ZC3HAV1-bound RNAs, we included HSPG2 mRNA as a representative ZC3HAV1-associated host RNA control. In contrast to the viral RNA results, neither WT N nor Mut-N significantly reduced ZC3HAV1 association with HSPG2 mRNA (S5F Fig). Although analysis of a single host transcript cannot exclude effects on all cellular RNAs, this result argues against a general nonspecific dissociation of ZC3HAV1–RNA complexes under our assay conditions. Together, the competitive RIP and displacement RIP experiments support a viral RNA-biased antagonism model in which FIPV N interferes with ZC3HAV1–viral RNA association through coordinated ZC3HAV1 interaction and viral RNA-binding activity.
To determine whether FIPV N-mediated attenuation of ZC3HAV1-associated restriction was strictly dependent on the FIPV replication context, we used Japanese encephalitis virus (JEV) as a heterologous positive-strand RNA virus model. JEV was selected because it replicates in the cytoplasm and provides an experimentally tractable RNA virus system in which ZC3HAV1-associated antiviral activity can be evaluated independently of FIPV replication complexes. In BHK-21 cells infected with JEV, ZC3HAV1 expression reduced JEV NS1 protein accumulation, whereas FIPV N expression was associated with increased NS1 protein levels relative to the ZC3HAV1-expressing condition (Fig 6D–E). When FIPV N was co-expressed with ZC3HAV1, JEV NS1 protein expression was partially restored compared with ZC3HAV1 expression alone (Fig 6F–G). These results indicate that FIPV N can attenuate ZC3HAV1-associated restriction in a heterologous RNA virus system at the level of viral protein accumulation, suggesting that this counter-restriction activity is not strictly limited to the FIPV replication context.
Overall, these results support a coordinated antagonism model in which FIPV N uses both protein–protein and protein–RNA interactions to counteract ZC3HAV1-associated RNA-level restriction. The N–ZC3HAV1 interaction likely positions N at the ZC3HAV1–viral RNA interface, whereas the intrinsic viral RNA-binding activity of N contributes to reducing ZC3HAV1-associated viral RNA. Thus, FIPV N does not appear to antagonize ZC3HAV1 solely by sterically blocking the ZC3HAV1 RBD-containing region; rather, it can disrupt pre-formed ZC3HAV1–viral RNA complexes in a manner that depends, at least in part, on N-mediated viral RNA binding. The HSPG2 control further supports that this effect is biased toward viral RNA rather than reflecting global disruption of ZC3HAV1-associated RNAs. The heterologous JEV assay provides additional evidence that FIPV N-mediated attenuation of ZC3HAV1-associated antiviral activity can be observed outside the FIPV infection system.
Discussion
Coronaviruses have evolved intricate molecular strategies to evade host RNA surveillance and preserve viral RNA accumulation and persistence. These RNA surveillance and decay systems represent an essential first line of intracellular defense, capable of recognizing and degrading non-self RNAs before the activation of interferon signaling [23,24]. Among coronavirus structural proteins, the nucleocapsid (N) protein has been increasingly implicated in this process, functioning not only in RNA encapsidation and replication but also in modulating host antiviral pathways [25]. In the case of feline infectious peritonitis virus (FIPV), a coronavirus responsible for severe and often fatal systemic disease in cats, understanding how N circumvents host RNA restriction is particularly important, as efficient control of viral RNA synthesis and turnover is closely linked to the viral load dynamics and immunopathological outcome of feline infectious peritonitis (FIP). This study therefore used FIPV as a pathogenic model to dissect how coronavirus structural proteins antagonize host RNA-level defenses. Through proteomic screening and functional validation, we identified the host antiviral restriction factor ZC3HAV1 as a major host factor targeted by FIPV N, revealing a previously unrecognized post-transcriptional counter-restriction mechanism.
Proteomic profiling revealed that the FIPV N protein engages an extensive network of host RNA-binding and translational regulators, forming an interaction landscape partially overlapping with those described for SARS-CoV-2 [26] and mouse hepatitis virus (MHV) [27]. In these viruses, N has been shown to interact with ribosomal proteins and stress-granule components, reshaping host RNA regulatory environments to favor viral translation and RNA stability [28,29]. Functional validation confirmed that several of these N-associated host proteins—including ZC3HAV1 and MOV10, both known for their post-transcriptional antiviral roles [30,31]—restricted FIPV replication. Among these candidates, ZC3HAV1 was particularly notable because of its established role as an RNA-binding restriction factor, prompting us to further define how ZC3HAV1 restricts FIPV and how FIPV N counteracts this restriction.
Given its established role as an RNA-binding antiviral factor, we next examined the level at which ZC3HAV1 restricts FIPV replication. ZC3HAV1 reduced both genomic and subgenomic viral RNA accumulation, consistent with its previously reported restriction of SARS-CoV-2 [32]. Importantly, reduced steady-state gRNA and sgRNA levels measured by RT-qPCR alone cannot determine whether ZC3HAV1 acts by inhibiting viral RNA synthesis or by reducing the stability of already synthesized viral RNAs. To distinguish these possibilities, we performed a GS-441524-based RNA persistence assay. Because GS-441524 strongly inhibits coronavirus RNA-dependent RNA polymerase activity [33,34], treatment after viral RNAs had accumulated limited ongoing viral RNA synthesis, thereby enabling assessment of pre-existing FIPV RNA persistence. Under these conditions, ZC3HAV1 accelerated the loss of FIPV genomic and subgenomic RNAs, indicating that ZC3HAV1-mediated restriction is not primarily attributable to reduced viral RNA synthesis or sgRNA transcription. Rather, these data support a post-synthetic RNA destabilization mechanism in which ZC3HAV1 decreases the persistence of FIPV RNAs after they are produced. Although this finding is consistent with direct or cofactor-assisted destabilization of ZC3HAV1-bound viral RNAs, the specific RNA-decay pathway and host cofactors involved in FIPV RNA loss remain to be identified. This mechanistic refinement is important because it places ZC3HAV1 restriction at the level of viral RNA persistence and provides a direct rationale for why FIPV N, by disrupting ZC3HAV1–viral RNA association, can preserve viral RNA accumulation during infection.
Having established that ZC3HAV1 restricts FIPV replication at the RNA level, we next sought to elucidate how it recognizes viral transcripts. Remarkably, ZC3HAV1 efficiently restricted FIPV despite the virus’s low genomic CpG content, a feature shared by many coronaviruses, including SARS-CoV-2 [17,35]. In our CpG profiling analysis, FIPV showed broad CpG depletion across most of the genome, with only localized enrichment near the 5′ region. Further mapping of the 5′ UTR showed that this CpG enrichment did not coincide with the canonical leader–body junction shared by coronavirus sgRNAs. These observations suggest that global CpG frequency, and even local CpG distribution within the leader-containing region, is insufficient to fully explain ZC3HAV1-associated restriction of FIPV RNA. This interpretation is consistent with previous studies showing that ZC3HAV1-mediated restriction is not always strictly predicted by global CpG abundance and that ZC3HAV1 can restrict viruses with relatively low CpG frequencies or through context-dependent RNA recognition modes [36,37]. Therefore, rather than demonstrating complete CpG independence, our data support a noncanonical mode of ZC3HAV1 engagement in a CpG-suppressed coronavirus RNA context.
Our RIP assays further showed that full-length ZC3HAV1 and its RBD-containing region enriched both FIPV gRNA and multiple sgRNAs, whereas the WWE + PARP region showed minimal enrichment. These data indicate that the RBD-containing region is required for efficient association of ZC3HAV1 with FIPV RNAs in infected cells. Notably, several sgRNAs encoding structural proteins, including N and S, were strongly enriched, which may reflect their high intracellular abundance and/or RNA-context features that promote ZC3HAV1 engagement. Because RIP does not define the precise RNA motif or structure recognized by ZC3HAV1, additional features—including localized CpG-enriched regions, RNA abundance, secondary structure, local sequence context, or viral ribonucleoprotein organization—may contribute to ZC3HAV1 engagement of FIPV RNAs [15,38]. This provides a mechanistic framework for understanding how FIPV, despite overall CpG depletion, remains sensitive to ZC3HAV1-associated RNA-level restriction.
The localized CpG enrichment near the FIPV 5′ region may also contribute to ZC3HAV1 engagement. Although this enriched region does not overlap precisely with the canonical leader–body junction, it could influence ZC3HAV1 recognition through local RNA structure, shared 5′ RNA architecture, or viral ribonucleoprotein organization. Thus, our data do not exclude a role for the 5′ CpG-enriched region, but rather indicate that FIPV RNA recognition is not explained by global CpG frequency alone or by a simple leader-junction CpG model. Targeted mutational mapping of the 5′ CpG-enriched region, leader-containing reporter RNAs, or FIPV minigenome/replicon systems will be useful to determine whether this region directly contributes to ZC3HAV1 binding and restriction.
These observations also suggest that ZC3HAV1 recognition during FIPV infection is unlikely to be restricted to a single viral RNA species or dictated solely by a simple global CpG threshold. The ability of FIPV N to attenuate ZC3HAV1-mediated restriction in a heterologous JEV system, together with previous observations for PEDV N in an influenza virus system [18], supports the view that ZC3HAV1 is capable of engaging viral RNAs in different experimental contexts. In this context, our RIP-qPCR data identify FIPV gRNA and sgRNAs as relevant ZC3HAV1-associated viral RNA species during FIPV infection. Importantly, HSPG2 was included as a representative ZC3HAV1-associated non-viral RNA control in the displacement RIP assay. The finding that FIPV N efficiently reduced ZC3HAV1-associated FIPV RNAs but did not significantly reduce ZC3HAV1 association with HSPG2 suggests that N-mediated displacement is not a general nonspecific disruption of all ZC3HAV1–RNA complexes. Instead, this result supports a viral RNA-biased antagonism model in which FIPV N preferentially interferes with ZC3HAV1 association with viral RNAs, likely through coordinated N–ZC3HAV1 interaction and N-mediated viral RNA binding. Nevertheless, transcriptome-wide approaches such as ZC3HAV1 RIP-seq, CLIP-seq, or related RNA interactome analyses will be required to define the full range of viral RNA substrates and binding features recognized by ZC3HAV1 during coronavirus infection [38].
Our data further refine the mechanism by which FIPV N counteracts ZC3HAV1-mediated restriction. A simple competition model, in which N simply masks the ZC3HAV1 RNA-binding region and prevents viral RNA access, is insufficient to explain the displacement RIP results. In the redesigned assay, WT N reduced ZC3HAV1-associated FIPV RNAs after ZC3HAV1–viral RNA complexes had already formed, whereas the RNA-binding-impaired Mut-N retained detectable ZC3HAV1 interaction but showed markedly weaker displacement activity. These results indicate that N-mediated antagonism requires more than protein–protein association with ZC3HAV1. Instead, FIPV N likely acts through a coordinated mechanism in which its interaction with ZC3HAV1 positions N at the ZC3HAV1–viral RNA interface, while its intrinsic viral RNA-binding activity contributes to reducing ZC3HAV1-associated viral RNA. This model also explains why disruption of ZC3HAV1–viral RNA association, rather than degradation of ZC3HAV1 itself, represents a major antagonistic mechanism mediated by FIPV N. Mechanistically, this antagonism expands the known spectrum of viral countermeasures against ZC3HAV1. We show that the FIPV N protein associates with the RBD-containing region of ZC3HAV1 and impairs its association with viral RNA. Together with the displacement RIP and Mut-N analyses, these findings support a coordinated antagonistic mechanism in which FIPV N uses both ZC3HAV1 association and intrinsic viral RNA-binding activity to reduce ZC3HAV1-associated viral RNA recovery. This mechanism differs from previously described antagonists such as influenza A virus NS1, which suppresses ZC3HAV1/ZAP-mediated RNA-level restriction indirectly by disrupting cofactor-dependent pathways [39], and vaccinia virus C16, which sequesters ZC3HAV1/ZAP into cytoplasmic puncta to attenuate its antiviral activity [40]. To our knowledge, these findings provide evidence that a coronavirus structural protein can counteract ZC3HAV1-associated RNA restriction by targeting the ZC3HAV1–viral RNA axis, rather than through proteolytic or transcriptional suppression of ZC3HAV1. Through these coupled activities, N may function not only as an RNA chaperone but also as a viral antagonist that helps preserve viral RNA accumulation and persistence under ZC3HAV1-mediated restriction. This mechanism identifies a previously underappreciated RNA-level counter-restriction strategy by which FIPV N helps sustain viral RNA accumulation during infection (Fig 7).
[图:Proposed model of FIPV N-mediated counteraction of ZC3HAV1-associated RNA-level restriction. ZC3HAV1 associates with FIP]

A remaining question is which downstream host cofactors execute ZC3HAV1-associated restriction of FIPV RNAs. In other viral systems, ZC3HAV1 has been reported to cooperate with RNA-surveillance or decay-associated factors, including TRIM25, KHNYN, OAS3, and related host proteins, to promote antiviral RNA restriction [41–43]. Although these cofactors provide important candidates for future investigation, the present study focuses on defining the upstream virus–host interface that enables ZC3HAV1-mediated restriction and FIPV N-mediated counter-restriction. Specifically, ZC3HAV1 associates with FIPV genomic and subgenomic RNAs, reduces their accumulation and persistence, and FIPV N counteracts this process by disrupting ZC3HAV1–viral RNA association. By combining RNA persistence assays, RIP-qPCR, and displacement RIP analyses, our study moves beyond steady-state viral RNA quantification and links ZC3HAV1-associated RNA restriction to a specific viral antagonism mechanism. Further mechanistic resolution will require targeted depletion or rescue of candidate cofactors, broader interactome analyses, and RNA-centered approaches such as FIPV replicon systems, transcription-defective mutants, metabolic RNA labeling, or Northern blotting. These approaches would help determine whether canonical ZC3HAV1 cofactors or FIPV-specific RNA-decay pathways execute viral RNA loss and would further clarify the relative contributions of viral RNA synthesis, processing, and decay to ZC3HAV1-mediated restriction.
To further test whether FIPV N-mediated attenuation of ZC3HAV1-associated restriction was strictly confined to the FIPV replication context, we used JEV as a heterologous positive-strand RNA virus model. The ability of FIPV N to partially alleviate ZC3HAV1-mediated restriction of JEV suggests that this counter-restriction activity can be observed outside the FIPV infection system. This finding supports the functional transferability of FIPV N-mediated ZC3HAV1 antagonism and indicates that the activity is not entirely dependent on FIPV-specific replication machinery. However, additional coronaviruses and unrelated RNA viruses will be needed to define the breadth of viral RNA contexts in which FIPV N can counteract ZC3HAV1-associated restriction.
Collectively, these findings suggest that the FIPV N protein has an additional RNA-protective role beyond its established functions in viral RNA packaging, replication support, and immune modulation [44,45]. This work links the classical RNA-binding function of coronavirus N proteins to the evasion of cellular RNA-level restriction and identifies N-mediated disruption of the ZC3HAV1–viral RNA interface as a previously underappreciated mode of coronavirus counter-restriction.
Conclusions
In summary, this study uses a clinically fatal feline coronavirus as a model to reveal an RNA-level counter-restriction mechanism in which the nucleocapsid protein engages ZC3HAV1 through its interaction with the RNA-binding domain-containing region and functionally counteracts ZC3HAV1-mediated antiviral restriction. Functional analysis of N truncation mutants further indicates that the ZC3HAV1-associated region of FIPV N contributes to efficient antagonism of ZC3HAV1-mediated restriction. In addition, FIPV N interferes with ZC3HAV1–viral RNA association, thereby helping preserve viral RNA accumulation and persistence. These findings provide a mechanistic explanation for how FIPV helps sustain abundant genomic and subgenomic RNAs despite the presence of a host RNA-binding restriction factor, and extend the functional role of coronavirus N protein from viral RNA packaging and replication support to active modulation of host RNA-level restriction. Importantly, our data support a model in which ZC3HAV1-mediated restriction acts, at least in part, at a post-synthetic stage of the viral RNA life cycle, limiting the persistence of already synthesized FIPV RNAs rather than being explained solely by reduced viral RNA synthesis or altered sgRNA transcription. Although the precise downstream RNA-decay machinery remains to be defined, the coordinated targeting of the ZC3HAV1–viral RNA restriction axis by FIPV N represents a specific host–virus interface that may provide a conceptual basis for future antiviral strategies aimed at restoring intrinsic RNA restriction and limiting coronavirus RNA persistence.
Supporting information
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
Acknowledgments
We sincerely acknowledge the Public Laboratory of the Shanghai Veterinary Research Institute, Chinese Academy of Agricultural Sciences, for providing experimental facilities. We also gratefully acknowledge the financial support from the National Natural Science Foundation of China.
Data Availability
All relevant data are within the manuscript and its Supporting Information files.
Funding Statement
This study was financially supported by the National Natural Science Foundation of China (32172832 to GQL), the National Key R&D Program of China (2023YFD1800705 to GQL), the Central Public-interest Scientific Institution Basal Research Fund (2026JB09 to YQZ), the China Postdoctoral Science Foundation (2024M753586 to YQZ), the Central Public-interest Scientific Institution Basal Research Fund (2025JB01 AXT) and the Shanghai Municipal Science and Technology Major Project (ZD2021CY001 to GQL). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
References
[内容请查看原文]
Associated Data
Supplementary Materials
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
[内容请查看原文]
Data Availability Statement
All relevant data are within the manuscript and its Supporting Information files.
兽研