首页|基于转录组测序技术分析猪德尔塔冠状病毒感染ST细胞机制的初步研究

基于转录组测序技术分析猪德尔塔冠状病毒感染ST细胞机制的初步研究

A preliminary study to analyse the mechanism of ST cell infection by porcine deltacoronavirus based on transcriptome sequencing technology

扫码查看
猪德尔塔冠状病毒(PDCoV)感染可导致仔猪出现水样腹泻和死亡,严重危害猪的健康.为研究PDCoV感染引起宿主细胞基因转录水平的变化及初步研究其感染机制,本研究提取PDCoV感染组和不感染该病毒的阴性对照组ST细胞样品总RNA,构建二者的cDNA文库后采用Illumina HiSeq 6000高通量测序平台进行转录组测序(RNA-Seq).结果显示,所有样品的测序错误率均≤0.03%;Q20/%和Q30/%分别达98%和94%以上;GC含量在51.09%~55.15%;相关性分析结果显示,阴性对照组及感染组细胞样品各组内的R值均>0.9,上述结果表明测序数据可靠,可用于后续转录组学数据分析.采用ggplot2软件分析两组细胞组内测序结果的相关性,并以相关系数R>0.9判定各组内测序数据的准确性;采用DESeq2软件并以P<0.05和log2(Fold change)>1为条件筛选两组细胞中转录显著差异基因;采用ClusterProfiler R软件对转录显著差异基因进行GO功能和KEGG信号通路的富集分析;采用Pheatmap软件对各组ST细胞中与免疫反应相关的转录显著差异基因进行聚类分析;随机选择9个转录显著差异基因经RT-qPCR验证RNA-Seq的结果.筛选结果显示,与阴性对照组细胞相比,共筛选出5 719个转录显著差异基因,其中3 573个转录显著上调基因,2 146个转录显著下调基因.GO功能和KEGG信号通路富集分析结果显示,转录显著上调的基因与ST细胞的免疫应答、防御反应等功能有关;转录显著下调的基因与细胞的氧化还原和新陈代谢等功能有关;转录显著上调的基因参与细胞因子受体互作、天然免疫反应相关、肿瘤坏死因子等的信号通路;转录显著下调的基因与细胞物质代谢信号通路有关,包括氨基酸代谢、脂代谢、乙醛酸和二羧酸代谢等信号通路.与免疫反应相关的转录显著差异基因的聚类分析结果显示,与阴性对照组细胞相比,PDCoV感染的ST细胞中白细胞介素基因(IL-1R1、1L-6、IL-10RB、IL-15和IL-12A等)、抗病毒基因(MX1、OAS)和干扰素基因(INFAR1、IFN-ALPHAOMEGA)的转录水平均显著上调.上述结果表明PDCoV感染后,ST细胞的天然免疫反应被激活,但细胞的代谢反应降低,这更有利于减弱病毒在ST细胞中的复制,起到一定的抗病毒效果.9个转录显著差异基因的RT-qPCR结果与RNA-seq结果基本一致,表明转录组测序数据可靠.本研究为进一步解析PDCoV的感染机制提供参考依据.
Porcine deltacoronavirus(PDCoV)infection can cause watery diarrhoea and death in piglets,which is a serious health hazard for pigs.In order to study the changes of gene transcription level in host cell caused by PDCoV infection,total RNA was extracted from ST cells that were infected with PDCoV and negative control group not infected with the virus,and their cDNA libraries were constructed and then subjected to transcriptome sequencing(RNA-Seq)by using the Illumina HiSeq 6000 high-throughput sequencing platform.The results showed that the sequencing error rate of all samples was≤0.03%;Q20%and Q30%were above 98%and 94%,respectively;the GC content ranged from 51.09%to 55.15%;the correlation analysis results showed that the R of each group was>0.9,which indicated that the sequencing data were reliable and could be used for the subsequent analysis of transcriptomics data.The ggplot2 software was used to analyse the correlation of the sequencing results within each group,and the accuracy of the sequencing data was determined by the correlation coefficient R>0.9;DESeq2 software was used to screen for significantly differentially transcribed genes that met the requirement of log2(Fold change)>1 and P<0.05;ClusterProfiler R software was used to perform GO function and KEGG signaling pathway enrichment analysis on genes with significant differences in transcription;pheatmap software was used to analyse the clustering of transcriptionally significantly different genes that related to immune response in each group;nine transcriptionally significanfly different genes were randomly selected to validate the results of RNA-Seq by RT-qPCR.The screening results showed that compared with the cells in the negative control group,a total of 5719 transcriptionally significantly different genes were screened,including 3573 transcriptionally significantly up-regulated genes and 2146 transcriptionally significantly down-regulated genes.The results of GO function and KEGG signalling pathway enrichment showed that the transcriptionally significantly up-regulated genes were related to the immune response and defence response of ST cells,and involved in the signalling pathways of cytokine receptor interactions,natural immune response,and tumour necrosis factor,etc,while the transcriptionally significantly down-regulated genes were related to the redox and metabolism of cells,and were related to cellular metabolism signalling pathways,including amino acid metabolism,lipid metabolism,glyoxylate and dicarboxylic acid metabolism signalling pathways.The clustering analysis of transcriptionally significantly different genes associated with immune response showed that compared with the negative control group,the transcriptional levels of interleukin genes(IL1R1,1L-6,IL-10RB,IL-15 and IL-12A,etc.),antiviral genes(MX1,OAS)and interferon genes(INFAR1,IFN-ALPHAOMEGA)were significantly up-regulated in the PDCoV-infected group.These results indicated that the natural immune response of ST cells was activated after PDCoV infection,but the metabolic response of the cells was reduced,which was more conducive to weakening the replication of the virus in ST cells and exerting a certain antiviral effect.The RT-qPCR results of the nine transcriptionally significantly different genes were largely consistent with the RNA-seq results,indicating the high reliability of the transcriptome sequencing data.This study provides a theoretical basis for further analysis of the infection mechanism of PDCoV.

porcine deltacoronavirusRNA-seqtranscription significantly different genes

谢永生、张昆丽、张建峰、贺东生

展开 >

华南农业大学兽医学院/广东省动物源性人兽共患病预防与控制重点实验室,广东广州 510642

宜春学院生命科学与资源环境学院,江西宜春 336000

广东省农业科学院动物卫生研究所/广东省畜禽疫病防治研究重点实验室,广东广州 510640

猪德尔塔冠状病毒 转录组测序 转录显著差异基因

国家重点研发技术项目(十四五)广东省农业科技创新十大主攻方向"揭榜挂帅"项目(十四五)广东省省级科技计划猪禽种业全国重点实验室项目猪禽种业全国重点实验室项目岭南现代农业科学与技术广东省实验室肇庆分中心项目宜春学院博士科研启动基金

2022YFD18008012022SDZG022023B12120600402023QZ-NK13ZQQZ-55P20211154-03023360122008

2024

中国预防兽医学报
中国农业科学院哈尔滨兽医研究所

中国预防兽医学报

CSTPCD北大核心
影响因子:0.674
ISSN:1008-0589
年,卷(期):2024.46(3)
  • 21