Comparative Proteomic Profiling: Cellular Metabolisms Are Mainly Affected in Senecavirus A-Inoculated Cells at an Early Stage of Infection

Comparative transcriptomics analysis on Senecavirus A-infected and non-infected cells

Senecavirus A (SVA) is an emerging virus that causes the vesicular disease in pigs, clinically indistinguishable from other high consequence vesicular diseases. This virus belongs to the genus Senecavirus in the family Picornaviridae. Its genome is a positive-sense, single-stranded RNA, approximately 7,300 nt in length, with a 3' poly(A) tail but without 5'-end capped structure. SVA can efficiently propagate in different cells, including some non-pig-derived cell lines. A wild-type SVA was previously rescued from its cDNA clone using reverse genetics in our laboratory. In the present study, the BSR-T7/5 cell line was inoculated with the passage-5 SVA. At 12 h post-inoculation, SVA-infected and non-infected cells were independently collected for the analysis on comparative transcriptomics. The results totally showed 628 differentially expressed genes, including 565 upregulated and 63 downregulated ones, suggesting that SVA infection significantly stimulated the transcription initiation in cells. GO and KEGG enrichment analyses demonstrated that SVA exerted multiple effects on immunity-related pathways in cells. Furthermore, the RNA sequencing data were subjected to other in-depth analyses, such as the single-nucleotide polymorphism, transcription factors, and protein-protein interactions. The present study, along with our previous proteomics and metabolomics researches, provides a multi-omics insight into the interaction between SVA and its hosts.

Cells at early and late stages of infection with Senecavirus A: Comparative analysis of N6-methyladenosine modification on mRNAs

Infection with Senecavirus A (SVA) causes differential phenotypes in cells. In this study, cells were inoculated with SVA for culture. At 12 and 72 h post infection, cells were independently harvested for high-throughput RNA sequencing, and further methylated RNA immunoprecipitation sequencing. The resultant data were comprehensively analyzed for mapping N⁶-methyladenosine (m⁶A)-modified profiles of SVA-infected cells. More importantly, m⁶A-modified regions were identified in the SVA genome. A dataset of m⁶A-modified mRNAs was generated for screening out differentially m⁶A-modified mRNAs, further subjected to a series of in-depth analyses. This study not only showed statistical differentiation of m⁶A-modified sites between two SVA-infected groups, but also demonstrated that SVA genome, as a positive-sense, single-stranded mRNA, itself could be modified through the m⁶A pattern. Out of the six samples of SVA mRNAs, only three were identified to be m⁶A-modified, implying that the epigenetic effect might not be a crucial driving force for SVA evolution.

The porcine piRNA transcriptome response to Senecavirus a infection

Introduction

Senecavirus A (SVA) belongs to the genus Senecavirus in the family Picornaviridae. PIWI-interacting RNAs (piRNAs) are a class of small Ribonucleic Acids (RNAs) that have been found in mammalian cells in recent years. However, the expression profile of piRNAs in the host during SVA infection and their roles are poorly understood.

Methods

Here, we found the significant differential expression of 173 piRNAs in SVA-infected porcine kidney (PK-15) cells using RNA-seq and 10 significant differentially expressed (DE) piRNAs were further verified by qRT-PCR.

Results

GO annotation analysis showed that metabolism, proliferation, and differentiation were significantly activated after SVA infection. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that significant DE piRNAs were mainly enriched in AMPK pathway, Rap1 pathway, circadian rhythm and VEGF pathway. It was suggested that piRNAs may regulated antiviral immunity, intracellular homeostasis, and tumor activities during SVA infection. In addition, we found that the expression levels of the major piRNA-generating genes BMAL1 and CRY1 were significantly downregulated after SVA infection.

Discussion

This suggests that SVA may affect circadian rhythm and promote apoptosis by inhibiting the major piRNA-generating genes BMAL1 and CRY1. The piRNA transcriptome in PK-15 cells has never been reported before, and this study will further the understanding of the piRNA regulatory mechanisms underlying SVA infections.

Comparative Proteomic Analysis Reveals Mx1 Inhibits Senecavirus A Replication in PK-15 Cells by Interacting with the Capsid Proteins VP1, VP2 and VP3

As an emergent picornavirus pathogenic to pigs, Senecavirus A (SVA) can replicate in pig kidneys and proliferates well in porcine kidney epithelial PK-15 cells. Here, tandem mass tags (TMT) labeling coupled with liquid chromatography–tandem mass spectrometry (LC-MS/MS) was used to analyze the proteome dynamic changes in PK-15 cells during SVA infection. In total, 314, 697 and 426 upregulated differentially expressed proteins (DEPs) and 131, 263 and 342 downregulated DEPs were identified at 12, 24 and 36 hpi, respectively. After ensuring reliability of the proteomic data by quantitative PCR and Western blot testing of five randomly selected DEPs, Mx1, eIF4E, G6PD, TOP1 and PGAM1, all the DEPs were subjected to multiple bioinformatics analyses, including GO, COG, KEGG and STRING. The results reveal that the DEPs were mainly involved in host innate and adaptive immune responses in the early and middle stages of SVA infection, while the DEPs mainly participated in various metabolic processes in the late stage of infection. Finally, we demonstrated that Mx1 protein exerts antiviral activity against SVA by interacting with VP1 and VP2 proteins dependent on its GTPase, oligomerization and interaction activities, while Mx1 interacts with VP3 only depending on its oligomerization activity. Collectively, our study provides valuable clues for further investigation of SVA pathogenesis.

Senecavirus A- and Non-Infected Cells at Early Stage of Infection: Comparative Metabolomic Profiles

Senecavirus A (SVA), classified into the genus Senecavirus in the family Picornaviridae, causes an infectious disease in pigs. This virus can efficiently replicate in some non-pig-derived cells, such as the BHK cell line and its derivative (BSR-T7/5 cell line). We had recovered a wild-type SVA from its cDNA clone previously, and then uncovered the proteomic profile of SVA-infected BSR-T7/5 cells at 12 h post inoculation (hpi). In order to explore the cellular metabolomics further, the SVA-inoculated BSR-T7/5 cell monolayer was collected at 12 hpi for assay via liquid chromatography-tandem mass spectrometry (LC-MS/MS). The resultant data set was comprehensively analyzed using bioinformatics tools. A total of 451 metabolites were identified using in-house and public databases. Out of these metabolites, sixty-one showed significantly differential values (p value < 0.05). The Kyoto Encyclopedia of Genes and Genomes (KEGG) database was used to analyze metabolic pathways of the significantly differential metabolites. There were eighty-one identified KEGG pathways, out of which twenty-seven showed their p values < 0.05. The pyrimidine metabolism revealed the minimum p value and the maximum number of significantly differential metabolites, implying the pyrimidine played a key role in cellular metabolism after SVA infection. SVA replication must rely on the cellular metabolism. The present study on metabolomics would shed light on impacts of SVA-induced multiple interactions among metabolites on cells or even on natural hosts.