畜牧兽医学报  2022, Vol. 53 Issue (1): 272-281. DOI: 10.11843/j.issn.0366-6964.2022.01.027    PDF    
引用本文
朱静静, 戴政列, 汪涵, 李向臣, 赵阿勇, 周晓龙, 杨松柏. 乙型脑炎病毒感染PK15细胞的lncRNA差异表达谱分析[J]. 畜牧兽医学报, 2022, 53(1): 272-281.  
ZHU Jingjing, DAI Zhenglie, WANG Han, LI Xiangchen, ZHAO Ayong, ZHOU Xiaolong, YANG Songbai. Analysis of Differential Expression Profile of LncRNA in PK15 Cells Infected with Japanese Encephalitis Virus[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(1): 272-281.  
乙型脑炎病毒感染PK15细胞的lncRNA差异表达谱分析
朱静静, 戴政列, 汪涵, 李向臣, 赵阿勇, 周晓龙, 杨松柏     
浙江农林大学动物科技学院·动物医学院, 临安 311300
收稿日期：2021-04-19
基金项目：浙江省自然科学基金（LY19C170003）；浙江省重点研发计划（2021C04034）
作者简介：朱静静(1995-), 安徽安庆人, 硕士生, 主要从事猪分子育种研究, E-mail: 2018120012009@stu.zafu.edu.cn.
通信作者：周晓龙, 主要从事猪遗传育种研究, E-mail: zhouxiaolong@zafu.edu.cn; 杨松柏, 主要从事动物遗传育种研究, E-mail: sbyang@zafu.edu.cn.
摘要：旨在分析乙型脑炎病毒（Japanese encephalitis virus，JEV）感染猪肾上皮细胞PK15后的lncRNA差异表达谱，探究宿主lncRNA在JEV感染和宿主防御中的潜在功能。本研究中利用JEV感染PK15细胞，感染36 h，收集细胞，同时设立未感染对照组，每组3个生物学重复。提取细胞总RNA，构建cDNA文库后进行高通量测序，筛选出差异表达lncRNA。通过对差异表达lncRNA与mRNA的位置关系以及序列相似性预测顺式和反式作用的靶基因，将预测的靶基因进行GO和KEGG分析。最后随机选出9个差异表达lncRNAs，利用实时荧光定量PCR（quantitative reverse transcription PCR，qRT-PCR）技术进行表达验证。结果表明，本研究中共鉴定出856个差异表达lncRNAs，其中452个lncRNAs显著上调表达，404个lncRNAs显著下调表达。GO分析发现，lncRNA的靶基因主要富集于先天性免疫反应；KEGG通路分析发现lncRNA的靶基因显著富集于肿瘤坏死因子、Toll样受体和NF-κB等信号通路。通过qRT-PCR验证发现，定量结果和测序结果基本一致。本研究为进一步探究lncRNA参与宿主免疫反应调控JEV复制奠定了理论基础。
关键词：猪    乙型脑炎病毒    先天免疫反应    PK15细胞    高通量测序    
Analysis of Differential Expression Profile of LncRNA in PK15 Cells Infected with Japanese Encephalitis Virus
ZHU Jingjing, DAI Zhenglie, WANG Han, LI Xiangchen, ZHAO Ayong, ZHOU Xiaolong, YANG Songbai     
College of Animal Science and Technology·College of Veterinary Medicine, Zhejiang A&F University, Lin'an 311300, China
Corresponding author: ZHOU Xiaolong, E-mail: zhouxiaolong@zafu.edu.cn; YANG Songbai, E-mail: sbyang@zafu.edu.cn.
Abstract: The objective of this study was to analyze the differential expression profile of lncRNA in pig kidney epithelial cells PK15 infected by Japanese encephalitis virus (JEV), and to explore the potential function of host lncRNA in JEV infection. In this study, JEV was used to infect PK15 cells for 36 hours, and the cells were collected. Meanwhile, an uninfected control group was established. There were 3 biological replicates in each group. The total RNA was extracted, and the cDNA libraries were constructed and then subjected to high-throughput sequencing to identify differentially expressed lncRNA. The cis- and trans-acting target genes of lncRNA were obtained by calculating the positional relationship and sequence similarity between differentially expressed lncRNA and mRNA. Next, the predicted target genes were analyzed by GO and KEGG. Finally, 9 differentially expressed lncRNAs were randomly selected for expression verification using real-time fluorescent quantitative PCR (quantitative reverse transcription PCR, qRT-PCR) technology. The results showed that a total of 856 differentially expressed lncRNAs were identified, of which 452 lncRNAs were significantly up-regulated and 404 lncRNAs were significantly down-regulated. GO analysis results showed that the target genes of lncRNAs were mainly enriched in innate immune response; KEGG pathway analysis results showed that the target genes of lncRNA were significantly enriched in signal pathways such as tumor necrosis factor, Toll-like receptor and NF-κB. The quantitative results are the same as the sequencing results through qRT-PCR verification. This study lays a theoretical foundation for further exploring the roles of lncRNA in JEV infection.
Key words: swine    Japanese encephalitis virus    innate immune response    PK15 cells    high-throughput sequencing    

乙型脑炎病毒(Japanese encephalitis virus, JEV)感染引起的乙型脑炎是一种严重危害人类健康和养猪业的人畜共患病。乙型脑炎主要流行于亚洲和西太平洋地区，据统计有超过30亿人口生活在这些地区[1]。其中我国是乙型脑炎的高发区，曾一度占世界的80%。尽管一些国家在疫苗接种工作中作出了很大努力，乙脑的发病率也已经有了明显下降，但每年仍有约6.8万人感染，其中75%是15岁以下青少年，并造成超过1万人死亡，而其中一半的幸存者会存在运动障碍、认知功能障碍和癫痫等后遗症[2]。猪是JEV的重要扩增宿主，蚊子是将病毒从猪传播给人类的重要载体[3]。因此，猪在JEV感染周期中起着至关重要的作用。此外，JEV感染也给养猪业造成了巨大的经济损失，母猪感染JEV引起流产和死胎，公猪感染引起无精症[4-5]。而猪肾上皮PK15细胞作为一种良好的感染模型，已被广泛应用于JEV感染研究[6-9]。

在人类基因组中，能够转录并翻译成蛋白质的序列只占2%，而基因组中有大量序列会转录成非编码RNA(ncRNA)，如果ncRNA的长度超过200个核苷酸，将它定义为长ncRNA(lncRNA)[10]。像mRNA一样，大多数lncRNA都具有5′帽结构，并且具有聚腺苷酸化尾巴。长期以来。由于lncRNA不编码蛋白质，因此被认为是“暗物质”或“垃圾RNA”[11]。后续越来越多研究结果表明，lncRNA参与多种生物学过程，包括细胞生长、分化和凋亡[12-13]。利用高通量测序发现宿主lncRNA在多种病毒感染期后差异表达，包括PRRSV[14]、PRV[15]和PEDV[16]，随后通过构建lncRNA-mRNA共表达网络，结果显示，多个差异表达lncRNA与免疫反应相关。

此外，已经证实lncRNA参与免疫反应过程来调控病毒复制。例如，lncRNA-Lnczc3 h7a结合TRIM25蛋白增强RIG-Ⅰ介导的抗病毒先天免疫反应[17]。lncRNA#32通过正向调节干扰素刺激基因表达抑制病毒复制[18]。LncRNA-IVRPIE可以通过促进干扰素β1和干扰素刺激基因的表达来抑制甲型流感病毒的复制[19]。利用微阵列芯片分析发现，在JEV感染小鼠后，在脑组织中鉴定出618个差异表达lncRNAs，富集分析发现，这些lncRNA在先天免疫反应和炎症反应中可能起关键作用[20]。但是，JEV感染后的猪lncRNA表达模式及lncRNA在JEV感染过程中的潜在功能仍不清楚。

本研究中，利用高通量技术对JEV感染后的PK15细胞进行转录组测序，鉴定JEV感染后差异表达lncRNA，预测其靶基因并对靶基因进行富集分析，最后随机选择9个差异表达lncRNAs，利用荧光定量PCR进行验证。本研究结果将为进一步探索lncRNA参与免疫反应调控JEV感染分子机制奠定基础。

1 材料与方法 1.1 细胞培养和病毒感染

将PK15细胞(中国典型培养物保藏中心)接种于6孔细胞培养板中，加入含有10%胎牛血清(FBS，HyClone)和1%非必需氨基酸(Gibco)MEM培养基(HyClone)，置于37 ℃和5%CO2培养箱中培养。第2天，待细胞融合度90%时，感染JEV(SA14-14-2，MOI=1)，吸附1 h后，将细胞用PBS洗涤3次，加入含2% FBS的MEM培养基置于细胞培养箱中继续培养。未感染的细胞用作对照组，每组试验重复3次。感染后36 h收集细胞提取RNA，并命名为JEV-1、JEV-2和JEV-3；对照组分别命名为Mock-1、Mock-2和Mock-3。

1.2 文库构建和测序

利用TRIzol试剂(Thermo Fisher)提取总RNA。随后利用NanoPhotometer分光光度计(IMPLEN)和1%琼脂糖凝胶电泳分别检测RNA的浓度、纯度和完整性。然后每个样品取3 μg RNA，通过Ribo-Zero Magnetic试剂盒(Epicenter)去除核糖体RNA，利用NEBNext Ultra RNA文库制备试剂盒(New England Biolabs)进行文库构建。最后，利用Illumina HiSeq 2500系统(Illumina)对文库进行测序。整个流程委托北京博奥晶典生物技术有限公司完成。

1.3 数据质控与lncRNA鉴定

为了保证数据的质量与可靠性，需要对测序的原始数据进行质控。包括删除带接头(adapters)的reads，去除包含N＞5%的reads和低质量(Q≤20)的reads，从而获得clean reads。

利用TopHat2将clean reads比对到参考基因组上并使用Cufflinks软件进行组装。利用Cuffcompare软件将组装好转录本和参考基因组进行比较，得到lncRNA和mRNA的相互位置关系。选用CPC、CNCI和Pfam软件对lncRNA进行编码潜能预测，3个软件预测都没有编码能力的即为候选lncRNA。

1.4 表达量分析

使用cuffquant和cuffnorm软件进行表达量分析，利用FPKM进行标准化处理样本，确定lncRNA和mRNA的表达水平。然后使用Deseq2软件鉴定差异表达lncRNA和mRNA，筛选标准为log2Fold Change≥1且P value≤0.05。

1.5 靶基因预测

LncRNA可以调控自身和其邻近编码基因的转录。因此根据lncRNA及其邻近编码基因之间的位置关系，将lncRNA分为5种类型：有义、反义、双向、内含子和基因间lncRNA。有义，当它们与蛋白质编码基因的有义链重叠时；反义，当它们与蛋白质编码基因的反义链重叠时；双向，当它们位于与蛋白质编码基因相反的链上时，其转录起始距离少于1 000个碱基对；内含子，当它们完全来自另一个转录本的内含子时；基因间，不在任何其他蛋白质编码基因座的10 000个碱基之内。将这5种类型的lncRNA归类为顺式调控lncRNA。使用blat软件对lncRNA和mRNA(3′UTR)序列进行比对，筛选序列相似lncRNA-mRNA对。因此得到多个顺式作用和反式作用的lncRNA-mRNA对，在这些lncRNA-mRNA对中，只有lncRNA和mRNA都为显著差异表达，才判定为lncRNA的靶基因。

1.6 GO和KEGG富集分析

利用GO(gene ontology)和KEGG(kyoto encyclopedia of genes and genomes)对差异表达lncRNA的顺式和反式靶基因进行功能富集分析。利用超几何分布进行假设检验得到富集结果P值，如果P≤0.05，则认为该富集结果显著。

1.7 荧光定量PCR验证差异表达lncRNA

随机选择9个差异表达lncRNAs进行qRT-PCR验证。在JEV感染后的不同时间点(12、24、36、48 h)收集PK15细胞，提取总RNA，每个样品3个重复。然后利用反转录试剂盒(康为世纪公司)获得cDNA，通过SYBR Green定量试剂盒(TaKaRa)进行荧光定量PCR反应。PCR反应体系：5 μL SYBR Green Mater Mix，上下游引物各0.3 μL, cDNA 1 μL，ddH2O 3.4 μL，总体系为10 μL。反应条件：95 ℃预变性5 min；95 ℃变性15 s, 60 ℃退火30 s, 72℃延伸20 s, 40个循环。依据2－ΔΔCt方法计算lncRNA相对表达量，GAPDH作为内参基因。数据表示为“平均值(x)±标准误差(sx)”，每个试验3个生物学重复，统计学分析采用双尾T检验，*.P＜0.05；**.P＜0.01。引物序列见表 1。

表 1 实时荧光定量PCR引物 Table 1 Primers for qRT-PCR
2 结果 2.1 测序数据评估

对JEV感染组和未感染组共6个样品进行高通量测序，得到70.6~95.3 Mb不等的原始序列。经过质控之后，共得到68.5~92.7 Mb不等的有效数据，有效数据比例达到了96.8 %以上，Q20在97.7 %以上，Q20在93.5 %以上(表 2)。测序数据和猪参考基因组的比对率均达到了70.3%以上，97.9 %以上的序列在参考基因上有唯一的比对位置(表 3)。表明测序数据较好，满足后续的数据分析。

表 2 测序数据质控 Table 2 Quality control of sequencing data
表 3 参考基因组比对情况 Table 3 Reference genome alignment
2.2 lncRNA鉴定与差异表达分析

通过CPC、CNCI和Pfam软件进行筛选后，本研究共鉴定出5 059个新lncRNAs(图 1A)，

A. CPC、CNCI和Pfam预测的lncRNA的维恩图；B. 鉴定的lncRNA和mRNA的数量；C. lncRNA和mRNA的外显子数目；D. lncRNA和mRNA的长度 A. The Venn diagram of lncRNAs predicted by CPC, CNCI, and the Pfam database; B. The number of lncRNAs and mRNAs identified; C. Exon numbers of lncRNAs and mRNAs; D. Transcript lengths of the lncRNAs and mRNAs 图 1 差异lncRNA的鉴定与基因组特征 Fig. 1 Prediction and genomic features of lncRNAs

已知lncRNA共3 435个，已知mRNA共12 818个，新mRNA共23 440个(图 1B)。另外分别对lncRNA和mRNA外显子数量和转录本长度进行比较分析，发现和mRNA相比，lncRNA具有较少的外显子数量和较短的转录本(图 1C、D)。对鉴定出的lncRNA和mRNA进行差异表达分析，结果表明，差异表达的lncRNA共856个，其中452个显著上调表达，404个显著下调表达(图 2A、B)；差异表达的mRNA总数为4 165个，其中显著上调表达的为2 253个，显著下调表达的为1 912个(图 2A)。从热图分析来看，同一个lncRNA或mRNA在同一组的表达模式相似，说明组内的样本重复性较好(图 2C、D)。

A. JEV感染组中上调和下调的差异表达lncRNA和mRNA的数量；B. 差异表达lncRNA的火山图，红点表示上调的基因，绿点表示下调的lncRNAs，黑点表示没有明显变化的基因；C. 差异表达lncRNA的聚类热图分析；D. 差异表达mRNA的聚类热图分析，红色表示相对较高的表达水平，绿色表示相对较低的表达水平 A. The number of upregulated and downregulated differentially expressed lncRNAs and mRNAs in the JEV infected group compared to the mock infected group; B. Volcano plot of differentially expressed lncRNAs. Red dots refer to upregulated and green dots refer to downregulated lncRNAs, and the black dots represent genes with no significant changes; C. Clustering and heatmap analysis of differentially expressed lncRNA; D. Clustering and heatmap analysis of differentially expressed mRNA. Red refers to relative higher expression level, and green refers to relatively lower expression level 图 2 差异表达lncRNA和mRNA鉴定与表达分析 Fig. 2 Identification and expression analysis of differentially expressed lncRNA and mRNA
2.3 差异表达lncRNA靶基因GO和KEGG分析

在差异表达lncRNA和mRNA中，通过比较或者比对lncRNA和mRNA的位置关系以及序列相似性，得到差异表达lncRNA的顺式和反式作用的靶基因。结果显示，共得到潜在靶基因1 041个，其中顺式作用靶基因657个，反式作用靶基因384个。为了进一步探究这些靶基因参与的生物学功能，利用GO数据库，将差异lncRNA靶基因进行功能注释。GO功能注释包括生物学过程(biological process)、细胞组分(cellular component)和分子功能(molecular function)。结果显示，共有821个GO条目显著富集。3个功能分类最显著的条目分别是先天免疫反应(innate immune response)、呼吸链(respiratory chain)和NADH脱氢酶活性(NADH dehydrogenase activity)(图 3)。KEGG富集分析结果显示，lncRNA的靶基因共富集到156个通路中，其中10个通路显著富集，包括肿瘤坏死因子信号通路(TNF signaling pathway)、NF-κB信号通路(NF-κB signaling pathway)和Toll样受体信号通路(Toll-like receptor signaling pathway)(图 4)。

图 3 差异表达lncRNA靶基因GO分析 Fig. 3 GO analysis of target genes of differentially expressed lncRNAs
图 4 差异表达lncRNA靶基因KEGG分析 Fig. 4 KEGG analysis of target genes of differentially expressed lncRNA
2.4 定量PCR验证差异表达lncRNA

本研究利用荧光定量PCR方法对测序结果进行验证。首先将JEV感染PK15细胞，感染不同时间点收集细胞，提取RNA，反转录成cDNA，然后随机选取9个差异表达lncRNAs(6个显著上调表达，3个显著下调表达)，设计引物，进行荧光定量PCR扩增。结果显示，与0 h未感染组相比，6个显著上调表达lncRNAs(TCONS_00091590、TCONS_00061989、TCONS_00057383、TCONS_00035357、TCONS_00022902、TCONS_00071284)在JEV感染36 h显著上调表达，而且在感染的整个过程中呈现上调表达趋势；3个显著下调表达lncRNAs(TCONS_00035689、TCONS_00036059、TCONS_00087848)在JEV感染36 h显著下调表达，并且在JEV感染过程中呈现下调表达趋势(图 5)。表明荧光定量PCR结果和测序结果一致，说明测序结果可靠。

与0 h未感染组相比，*.P＜0.05；**.P＜0.01 Compared with 0 h of uninfected group, *.P < 0.05;**.P < 0.01 图 5 荧光定量PCR验证差异表达lncRNA Fig. 5 Validation of differentially expressed lncRNAs by qRT-PCR
3 讨论

目前，关于JEV感染后引起宿主基因mRNA及miRNA[21]等的表达谱变化已有报道，但JEV感染后引起的猪lncRNA表达谱变化尚未见报道。研究表明，lncRNA在调控病毒宿主基因互作以及病毒增殖中起到关键作用[22]。在一项小鼠的研究中，利用微阵列芯片方法在JEV感染的小鼠脑组织中鉴定出多个差异表达的lncRNAs，通过lncRNA-mRNA共表达网络分析，发现这些lncRNA可能调控JEV诱导的炎症反应[20]。本课题组前期通过荧光定量PCR技术检测到4个lncRNAs在JEV感染PK15细胞后差异表达，并鉴定出lncRNA-SUSAJ1显著抑制JEV增殖[23-24]。本研究中，通过高通量测序方法鉴定JEV感染PK15细胞后的lncRNA表达谱。结果显示，共鉴定出5 059个新lncRNA。和mRNA相比，这些lncRNA具有较短的转录本和较少的外显子，这些lncRNA特征和之前的研究结果一致[25-27]。进一步通过表达分析发现，和对照组相比，感染组有856个差异表达lncRNAs。

在病毒感染细胞的不同阶段，宿主基因以及lncRNA往往会呈现出不同的表达谱。杜程涛等[23]在研究JEV感染PK15细胞动力学发现，在感染后的36 h时，病毒的滴度已经达到较高水平，而且这时细胞并没有发生明显病变。因此本研究选择了JEV感染36 h的细胞进行转录组测序。为了验证测序数据的准确性以及选择JEV感染时间点的普遍代表性，利用JEV感染PK15细胞，在感染不同时间点(0、12、24、36和48 h)收集细胞，提取总RNA，反转录成cDNA。随后随机选取转录组测序筛选到的9个差异表达lncRNAs进行荧光定量PCR验证。结果显示，这些lncRNAs在感染后36 h以及在病毒感染的整个过程中的表达趋势和测序结果一致。说明测序结果的可靠性高以及选择感染时间点具有代表性。

lncRNA能够通过顺式作用调控邻近基因的表达或通过反式作用调控不同染色体上的基因表达[28]。因此，本研究通过判定lncRNA和mRNA的相对位置以及序列相似性预测差异表达lncRNA的顺式和反式调控的靶基因。进一步通过GO分析发现这些靶基因在生物学功能中最显著富集于先天免疫反应。说明这些差异表达lncRNA可以通过调控宿主先天性免疫反应激活模式识别受体信号通路从而诱导干扰素和细胞因子的表达。例如，在一项丙肝病毒(HCV)的研究中发现，干扰素处理或HCV感染能诱导lncRNA-CMPK2的表达，同时在干扰素刺激的肝细胞中，敲低lncRNA-CMPK2表达后正调控干扰素信号通路进而抑制HCV复制[29]。lncRNA-NRAV通过组蛋白修饰抑制干扰素刺激基因的转录进而促进甲型流感病毒的复制[30]。lncRNA-ITPRIP-1通过激活MDA5诱导先天免疫反应，而lncRNA-EPS通过调控核小体的状态抑制免疫反应基因转录[31-32]。lncRNA-LncATV抑制RIG-Ⅰ介导的先天免疫从而促进多种RNA病毒(包括丙型肝炎病毒、寨卡病毒)的复制[33]。因此本研究中筛选到的差异表达lncRNA参与先天性免疫反应调控JEV复制需要进一步验证。

差异表达lncRNA靶基因的KEGG富集结果中，显著富集的通路包括肿瘤坏死因子信号通路、NF-κB信号通路和Toll样受体信号通路。研究表明，Toll样受体(TLR)通过识别病毒复制产物，传递一系列信号激活细胞因子和Ⅰ类干扰素表达启动先天性免疫反应发挥抗病毒作用，其中TLR3、TLR7和TLR8可以识别病毒RNA[34]。Toll样受体在宿主抗JEV感染发挥关键作用。例如，敲低TLR3基因表达后显著促进JEV复制[35-36]。在小鼠的脑组织中特异性敲低TLR7基因的表达后再感染JEV，发现在脑组织中IFN-α和抗病毒蛋白显著降低，同时病毒的载量显著升高[37]。研究表明，lncRNA能调控Toll样受体信号通路基因表达。比如lncRNA MEG3通过抑制TLR4信号通路基因表达促进呼吸道合胞病毒(RSV)的感染[38]。本研究筛选到多个差异表达lncRNA靶向免疫相关信号通路，因此, 推测这些lncRNA介导先天性免疫反应相关通路调控JEV感染。

4 结论

利用高通量测序技术对乙型脑炎病毒(JEV)感染后的PK15细胞进行了转录组测序，筛选到多个差异表达lncRNA，并预测这些差异表达lncRNA的顺式作用和反式作用的靶基因。进一步通过GO和KEGG富集分析发现，这些靶基因主要富集于免疫反应相关通路。推测lncRNA通过参与免疫反应调控JEV复制。本研究为进一步从lncRNA角度探究宿主基因参与调控JEV感染提供了理论依据。

参考文献
[1]
HEFFELFINGER J D, LI X, BATMUNKH N, et al. Japanese encephalitis surveillance and immunization-Asia and Western Pacific Regions, 2016[J]. MMWR Morb Mortal Wkly Rep, 2017, 66(22): 579-583. DOI:10.15585/mmwr.mm6622a3
[2]
NANISHI E, HOSHINA T, SANEFUJI M, et al. A nationwide survey of pediatric-onset Japanese Encephalitis in Japan[J]. Clin Infect Dis, 2019, 68(12): 2099-2104. DOI:10.1093/cid/ciy816
[3]
OLIVEIRA A R S, COHNSTAEDT L W, NORONHA L E, et al. Perspectives regarding the risk of introduction of the japanese encephalitis virus (JEV) in the United States[J]. Front Vet Sci, 2020, 7: 48. DOI:10.3389/fvets.2020.00048
[4]
TAKASHIMA I, WATANABE T, OUCHI N, et al. Ecological studies of Japanese encephalitis virus in Hokkaido: interepidemic outbreaks of swine abortion and evidence for the virus to overwinter locally[J]. Am J Trop Med Hyg, 1988, 38(2): 420-427. DOI:10.4269/ajtmh.1988.38.420
[5]
BURNS K F. Congenital Japanese B encephalitis infection of swine[J]. Proc Soc Exp Biol Med, 1950, 75(2): 621-625. DOI:10.3181/00379727-75-18285
[6]
LIU K, LIAO X W, ZHOU B, et al. Porcine alpha interferon inhibit Japanese encephalitis virus replication by different ISGs in vitro[J]. Res Vet Sci, 2013, 95(3): 950-956. DOI:10.1016/j.rvsc.2013.08.008
[7]
YANG S B, HE M H, LIU X D, et al. Japanese encephalitis virus infects porcine kidney epithelial PK15 cells via clathrin- and cholesterol-dependent endocytosis[J]. Virol J, 2013, 10(1): 258. DOI:10.1186/1743-422X-10-258
[8]
CAI Y H, ZHU L, ZHOU Y C, et al. Identification and analysis of differentially-expressed microRNAs in Japanese encephalitis virus-infected PK-15 cells with deep sequencing[J]. Int J Mol Sci, 2015, 16(1): 2204-2219. DOI:10.3390/ijms16012204
[9]
YANG S B, PEI Y, LI X Y, et al. miR-124 attenuates Japanese encephalitis virus replication by targeting DNM2[J]. Virol J, 2016, 13(1): 105. DOI:10.1186/s12985-016-0562-y
[10]
LAMSISI M, ENNAJI M M. Involvement and roles of long noncoding RNAs in the molecular mechanisms of emerging and reemerging viral infections[J]. Emerg Reemerg Viral Path, 2020, 2: 71-92.
[11]
DOOLITTLE W F. Is junk DNA bunk?A critique of ENCODE[J]. Proc Natl Acad Sci U S A, 2013, 110(14): 5294-5300. DOI:10.1073/pnas.1221376110
[12]
FATICA A, BOZZONI I. Long non-coding RNAs: new players in cell differentiation and development[J]. Nat Rev Genet, 2014, 15(1): 7-21.
[13]
HU W Q, ALVAREZ-DOMINGUEZ J R, LODISH H F. Regulation of mammalian cell differentiation by long non-coding RNAs[J]. Embo Rep, 2012, 13(11): 971-983. DOI:10.1038/embor.2012.145
[14]
ZHEN Y R, WANG F Q, LIANG W, et al. Identification of differentially expressed non-coding RNA in porcine alveolar macrophages from Tongcheng and Large white pigs responded to PRRSV[J]. Sci Rep, 2018, 8(1): 15621. DOI:10.1038/s41598-018-33891-0
[15]
JIN Y, ZHANG K L, HUANG W X, et al. Identification of functional lncRNAs in pseudorabies virus type Ⅱ infected cells[J]. Vet Microbiol, 2020, 242: 108564. DOI:10.1016/j.vetmic.2019.108564
[16]
CHEN J N, ZHANG C Y, ZHANG N, et al. Porcine endemic diarrhea virus infection regulates long noncoding RNA expression[J]. Virology, 2019, 527: 89-97. DOI:10.1016/j.virol.2018.11.007
[17]
LIN H Y, JIANG M H, LIU L, et al. The long noncoding RNA Lnczc3 h7a promotes a TRIM25-mediated RIG-Ⅰ antiviral innate immune response[J]. Nat Immunol, 2019, 20(7): 812-823. DOI:10.1038/s41590-019-0379-0
[18]
NISHITSUJI H, UJINO S, YOSHIO S, et al. Long noncoding RNA #32 contributes to antiviral responses by controlling interferon-stimulated gene expression[J]. Proc Natl Acad Sci U S A, 2016, 113(37): 10388-10393. DOI:10.1073/pnas.1525022113
[19]
ZHAO L N, XIA M, WANG K Y, et al. A long non-coding RNA IVRPIE promotes host antiviral immune responses through regulating interferon β1 and ISG expression[J]. Front Microbiol, 2020, 11: 260. DOI:10.3389/fmicb.2020.00260
[20]
LI Y C, ZHANG H, ZHU B B, et al. Microarray analysis identifies the potential role of long Non-Coding RNA in regulating neuroinflammation during Japanese encephalitis virus infection[J]. Front Immunol, 2017, 8: 1237. DOI:10.3389/fimmu.2017.01237
[21]
杜加茹, 陈浩南, 王方雨, 等. 乙型脑炎病毒感染小鼠原代神经元细胞的miRNA表达谱差异分析[J]. 畜牧兽医学报, 2020, 51(9): 2319-2323.
DU J R, CHEN H N, WANG F Y, et al. analysis of the miRNA expression profiles in the primary neurons of mice infected with Japanese encephalitis virus[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(9): 2319-2323. (in Chinese)
[22]
WANG P. The opening of Pandora's Box: an emerging role of long noncoding RNA in viral infections[J]. Front Immunol, 2018, 9: 3138.
[23]
杜程涛, 汪涵, 杨松柏, 等. lncRNA在乙脑病毒感染PK15细胞过程中的作用研究[J]. 中国畜牧兽医, 2019, 46(7): 2045-2052.
DU C T, WANG H, YANG S B, et al. Study on the role of lncRNA in the process of Japanese encephalitis virus infecting PK15 cells[J]. China Animal Husbandry & Veterinary Medicine, 2019, 46(7): 2045-2052. (in Chinese)
[24]
ZHOU X L, YUAN Q Y, ZHANG C, et al. Inhibition of Japanese encephalitis virus proliferation by long non-coding RNA SUSAJ1 in PK-15 cells[J]. Virol J, 2021, 18(1): 29. DOI:10.1186/s12985-021-01492-5
[25]
MA X L, ZHAO X M, WANG K L, et al. Identification and analysis of long non-coding RNAs that are involved in inflammatory process in response to transmissible gastroenteritis virus infection[J]. BMC Genomics, 2019, 20(1): 806. DOI:10.1186/s12864-019-6156-5
[26]
MAJEWSKA M, LIPKA A, PAUKSZTO L, et al. Preliminary RNA-Seq analysis of long non-coding RNAs expressed in human term placenta[J]. Int J Mol Sci, 2018, 19(7): 1894. DOI:10.3390/ijms19071894
[27]
SHI G L, CHEN L, CHEN G T, et al. Identification and functional prediction of long intergenic Non-coding RNAs related to subcutaneous adipose development in pigs[J]. Front Genet, 2019, 10: 160. DOI:10.3389/fgene.2019.00160
[28]
LIU W W, DING C. Roles of LncRNAs in viral infections[J]. Front Cell Infect Microbiol, 2017, 7: 205. DOI:10.3389/fcimb.2017.00205
[29]
KAMBARA H, NIAZI F, KOSTADINOVA L, et al. Negative regulation of the interferon response by an interferon-induced long non-coding RNA[J]. Nucleic Acids Res, 2014, 42(16): 10668-10680. DOI:10.1093/nar/gku713
[30]
OUYANG J, ZHU X M, CHEN Y H, et al. NRAV, a long noncoding RNA, modulates antiviral responses through suppression of interferon-stimulated gene transcription[J]. Cell Host Microbe, 2014, 16(5): 616-626. DOI:10.1016/j.chom.2014.10.001
[31]
XIE Q Y, CHEN S W, TIAN R Y, et al. Long noncoding RNA ITPRIP-1 positively regulates the innate immune response through promotion of oligomerization and activation of MDA5[J]. J Virol, 2018, 92(17): e00507-18.
[32]
ATIANAND M K, HU W Q, SATPATHY A T, et al. A long noncoding RNA lincRNA-EPS Acts as a transcriptional brake to restrain inflammation[J]. Cell, 2016, 165(7): 1672-1685. DOI:10.1016/j.cell.2016.05.075
[33]
FAN J J, CHENG M, CHI X J, et al. A human long non-coding RNA LncATV promotes virus replication through restricting RIG-Ⅰ-mediated innate immunity[J]. Front Immunol, 2019, 10: 1711. DOI:10.3389/fimmu.2019.01711
[34]
ARPAIA N, BARTON G M. Toll-like receptors: key players in antiviral immunity[J]. Curr Opin Virol, 2011, 1(6): 447-454. DOI:10.1016/j.coviro.2011.10.006
[35]
HAN Y W, CHOI J Y, UYANGAA E, et al. Distinct dictation of Japanese encephalitis virus-induced neuroinflammation and lethality via triggering TLR3 and TLR4 signal pathways[J]. PLoS Pathog, 2014, 10(9): e1004319. DOI:10.1371/journal.ppat.1004319
[36]
JIANG R, YE J, ZHU B B, et al. Roles of TLR3 and RIG-Ⅰ in mediating the inflammatory response in mouse microglia following Japanese encephalitis virus infection[J]. J Immunol Res, 2014, 2014: 787023.
[37]
NAZMI A, MUKHERJEE S, KUNDU K, et al. TLR7 is a key regulator of innate immunity against Japanese encephalitis virus infection[J]. Neurobiol Dis, 2014, 69: 235-247. DOI:10.1016/j.nbd.2014.05.036
[38]
TAO X W, ZENG L K, WANG H Z, et al. LncRNA MEG3 ameliorates respiratory syncytial virus infection by suppressing TLR4 signaling[J]. Mol Med Rep, 2018, 17(3): 4138-4144.

(编辑   白永平)