
==== Front
Biochem Biophys Rep
Biochem Biophys Rep
Biochemistry and Biophysics Reports
2405-5808
Elsevier

S2405-5808(24)00189-4
10.1016/j.bbrep.2024.101825
101825
Research Article
Exogenous interactome analysis of bovine viral diarrhea virus-host using network based-approach and identification of hub genes and important pathways involved in virus pathogenesis
Rezatofighi Seyedeh Elham e.tofighi@scu.ac.ir
e.tofighi@yahoo.com

Biology Department, Faculty of Science, Shahid Chamran University of Ahvaz, Ahvaz, Iran
16 9 2024
12 2024
16 9 2024
40 10182510 2 2024
8 9 2024
11 9 2024
© 2024 The Author. Published by Elsevier B.V.
2024

https://creativecommons.org/licenses/by/4.0/ This is an open access article under the CC BY license (http://creativecommons.org/licenses/by/4.0/).
Bovine viral diarrhea (BVD) is one of the most important diseases in livestock, caused by BVD virus (BVDV). During the pathogenesis of the virus, many interactions occur between host and viral proteins. Studying these interactions can help better understand the pathogenesis of the virus, identify putative functional proteins, and find new treatment and prevention strategies. To this aim, a BVDV-host protein-protein interaction (PPI) network map was constructed using Cytoscape and analyzed with cytoHubba, Kyoto Encyclopedia of Genes and Genomics (KEGG), Gene Ontology (GO), and Protein Analysis Through Evolutionary Relationships (PANTHER). Npro with 125 connections had the greatest number of interactions with host proteins. CD46, EEF-2, and TXN genes were detected as hub genes using different ranking algorithms in cytoHubba. BVDV interactions with its host mainly focus on targeting translation, protein synthesis, and cellular metabolism pathways. Different classes of proteins including translational proteins, nucleic acid metabolism proteins, metabolite interconversion enzymes, and protein-modifying enzymes are affected by BVDV. These findings improve our understanding of the effects of the virus on the cell. Hub genes and key pathways identified in the present study can serve as targets for novel BVDV prevention or treatment strategies.

Highlights

• BVDV-host interactions are mainly focused on targeting the translation, protein synthesis, and cellular metabolism pathways.

• Npro of BVDV has the greatest number of interactions with host proteins.

• CD46, EEF-2, and TXN genes are hob genes.

Keywords

BVDV
Cytoscape
GO
Hub gene
KEGG
Protein-protein interaction
==== Body
pmc1 Introduction

Bovine viral diarrhea (BVD) is a major disease in livestock, causing considerable economic losses to this industry. BVD is associated with gastrointestinal, respiratory, and reproductive disorders in cattle. The etiological agent of the disease is the BVD virus (BVDV), which belongs to the Flaviviridae family, Pestivirus genus. According to latest classification of Pestivirus, this genus has 11 species including pestivirus A (BVDV-1), pestivirus B (BVDV-2), pestivirus C (classical swine fever virus), pestivirus D (border disease virus), and pestiviruses E to K [1]. BVDV is a positive-sense single-stranded RNA virus that replicates in the cytoplasm. The RNA genome is approximately 12.3 kb in length and contains an open reading frame flanked by untranslated regions (UTRs). The ORF encodes a polyprotein precursor that is cleaved into 11 or 12 proteins in noncytopathic (ncp) or cytopathic (cp) virus biotypes, respectively. Mature viral proteins from N-terminus to C-terminus include Npro, capsid, Erns, E1, E2, p7, NS2-3, or NS2 and NS3, NS4A, NS4B, NS5A, and NS5B. BVDV-1 and BVDV-2 cause the same disease, but they are differ in the 5′-UTR sequences of their viral genomes [2].

Based on the cytopathic effects of BVDV on cell cultures, the virus is classified into two biotypes: cp and ncp [3]. The different features of BVDV biotypes result from differences in the processing of BVDV NS2-3. In cpBVDV infections, the NS2-3 protein is continuously cleaved by NS2 autoprotease activity to form NS3. However, in ncpBVDV infections, NS2-3 cleavage occurs during the early stages, and free NS3 is detected only up to 9 h post-infection. Therefore, the effects of NS2-3 and NS3 on cellular signaling appear to be different [4].

Viruses are intracellular parasites that need cells to replicate and survive. During various stages of viral infection, numerous interactions occur between host and viral proteins [5,6]. Investigating these interactions can improve our understanding of viral pathogenesis, help identify putative functional proteins, and facilitate the discovery of new treatment and prevention strategies. Host-viral protein-protein interactions (PPIs) have been studied through experimental approaches such as yeast two-hybrid, co-immunoprecipitation, tandem affinity purification-mass spectrometry, affinity chromatography, protein microarrays, protein-fragment complementation, and/or computational programs. Researchers have applied network biology approaches to construct disease networks that explore host-pathogen relationships and disease pathways [[6], [7], [8], [9], [10]]. Modeling virus-host interactome maps can help identify functionally important proteins [6], sub-networks, and potential hubs involved in the pathogenesis of specific diseases [7]. Researchers can then study these components using various bioinformatics methods. Viral proteins target bottlenecks, hubs, and rich clubs within cellular protein networks, thereby affecting cellular processes such as transcription, translation, cell signaling, cell trafficking, and cell cycle regulation [5,[11], [12], [13]].

The interaction of some BVDV proteins with host cell proteins has been previously investigated using experimental methods; however, there is no a single comprehensive network that encompasses all the research conducted on the BVDV interactome. This study aimed to combine previous works on BVDV-host interactions and construct a PPI network that encompasses all aspects of this topic. The goal of this study was to identify the proteins that significantly contribute to infection pathways, viral pathogenesis, and the virus's effects on host cells.

2 Material and methods

2.1 Data Collection

To perform this study, in the first step, all data related to the BVDV protein-protein interactions were collected. For this purpose, the advanced search option of PubMed was used. Multiple keywords related to interactions between BVDV and host proteins or virus proteins were applied. The search yielded 2690 research studies in which the virus was studied in different aspects including pathogenesis, epidemiology, and others. After a careful review of the articles, 24 papers contained the PPI data (Table 1) [2,4,[14], [15], [16], [17], [18], [19], [20], [21], [22], [23], [24], [25], [26], [27], [28], [29], [30], [31], [32], [33], [34], [35]]. In these articles, various PPI detection methods such as yeast two-hybrid, immunoprecipitation, co-immunoprecipitation and mass spectrometry were carried out. Besides the data, other PPI database including STRING Viruses [36], VirusMINT [37], IntAct [38], VirHostNet [39], and VirusMentha [40] were investigated. The information related to every single interaction of BVDV-1 or BVDV-2 with host were gathered. Because the BVDV-1 and BVDV-2 genes are identical and differ only in the 5′-UTR region, one network was drawn for both. In the context of virus–host interactions, protein-protein interfaces can be classified as exogenous or endogenous. Host-host or virus-virus PPIs are classified as endogenous interfaces, while interactions between different proteomes including virus–host PPIs are grouped as exogenous interfaces [41,42]. In the present study, host-host endogenous interfaces were excluded. All data gathered were convert to a unique format in order to construct a network in Cytoscape, because all data must be in the same format. Different identifiers including protein name, Uniprot, EMBL/GenBank, and others were converted to the same format, which here was the name gene. This work was done by an online Uniprot ID Mapping (http://www.uniprot.org/mapping/) [43]. Then, the data gathered from all resources were merged and duplicates were removed.Table 1 List of studies which reported the protein-protein interaction data of BVDV-host.

Table 1No	Paper	No. of interactions	References	
1	Weiskircher et al., 2009	2	2	
2	Yamane et al., 2009	2	4	
3	Lazar et al., 2003	1	14	
4	Hilton et al., 2003	1	15	
5	Zahoor et al., 2010	1	16	
6	Gong et al., 2020	1	17	
7	Johnson et al., 2001	1	18	
8	Leveringhaus et al., 2020	1	19	
9	Mu et al., 2021	1	20	
10	Branza-Nichita et al., 2004	1	21	
11	Fu et al., 2014	4	22	
12	Darweesh et al., 2018	1	23	
13	Schaut et al., 2015	1	24	
14	Xu et al., 1997	4	25	
15	Branza-Nichita et al., 2001	1	26	
16	Maurer et al., 2004	1	27	
17	Lackner et al., 2005	1	28	
18	Chen et al., 2007	1	29	
19	Dubrau et al., 2017	2	30	
20	Mohamed et al., 2014	1	31	
21	Wang et al., 2011	7	32	
22	Jefferson et al., 2014	89	33	
23	Ammari et al., 2017	72	34	
24	Gladue et al., 2014	31	35	

2.2 PPI network Construction and analysis

The latest version of Cytoscape (v3.9.0) [44] was used to construct the BVDV-host PPI network. Cytoscape is a user-friendly platform which help to construct, visualize and analyze of the biomolecular interaction networks. This platform has various tools and plugins for understanding and exploration of the various aspects of the network. Cytoscape Analyzer tool was used to define the topological characters of the network including diameter, radius, average number of neighbors, density, clustering coefficient, and so on.

2.3 Hub genes and modules analysis

CytoHubba (v0.1) [45], a plugin in Cytoscape was used to identify the hub genes. Top ten genes were selected with five algorithms in cytoHubba including MCC, Density of Maximum Neighborhood Component (DMNC), Edge Percolated Component (EPC), Radiality centrality, and Stress centrality. Different clustering algorithms including MCODE (v2.0.0) [46], SCODE (v1.0.4) [47], CytoCluster (v2.1.0) [48], and ClusterViz (v1.0.3) [49] were applied to construct modules and sub-networks. These are plugins in Cytoscape software for clustering nodes based on different algorithms. These plugins can find densely connected nodes [50], important protein complexes [2] and valuable functional motifs. Default parameters were applied for all plugins to construct sub-networks.

2.4 Functional enrichment analysis

Gene enrichment analysis was performed to evaluate the biological pathways of proteins that might be involved in the development and pathogenesis of BVDV. KEGG pathway enrichment analysis was performed using free web based search engine DAVID program (v6.8) with P < 0.05 threshold (https://david.ncifcrf.gov/) [51]. GO enrichment analysis was done using AgriGO (v2.0) to evaluate the CC, MF, and BP of involving host proteins. Interacting host proteins were analyzed with Singular Enrichment Analysis (SEA) tool. Gene names were converted to UniProtKB format because AgriGO allowed ID types including Bos taurus NCBI-GI, ENSEMBL ID, Bovine genome locus (Bovine Genome Database), GenBank ID, DDBJ ID, EMBL ID, UniProt ID, RefSeq Peptide ID, PDB ID, and Bovine Affymetrix Genome Array. The AgriGO can be found at (http://bioinfo.cau.edu.cn/agriGO/) [52].

BVDV-interacting host proteins were entered into the PANTHER classification program (v16.0). The PANTHER system classifies proteins and their genes according to family and subfamily, molecular function, biological process, and pathway. The PANTHER is online freely program and can be found at (http://www.pantherdb.org/) [53].

3 Results

3.1 The BVDV-host protein interaction network

After extracting data from literature and databases, the PPI network was constructed using Cytoscape. In total, 224 interactions were confirmed by experiment-derived data from previous studies. These interactions involved viral proteins including Npro (99 interactions), NS3 (77 interactions), E2 (34 interactions), Erns (5 interactions), NS5A (4 interactions), NS4A (one interaction), C (one interaction), E1 (one interaction), NS4B (one interaction), and E1 (one interaction) with host or viral proteins. The total number of interactions derived from previous studies and data bases was 561; however, after removing duplicates, 279 nodes and 436 interactions were identified (Fig. 1). Among the 279 nodes, 12 were viral proteins, while remaining were host proteins. Npro with 125 connections, had the greatest number of interactions with host proteins, followed by NS2-3, NS3, and E2, with 94, 87, and 65 associations, respectively. Fig. 2 shows the number of interactions between each BVDV-protein and host proteins. The PPI network was analyzed using the network analysis tool, and important statistical features are presented in Table 2.Fig. 1 BVDV-Host protein interaction network constructed and visualized using Cytoscape. Red nodes represent host genes, green nodes represent virus genes, and size of nodes represent the number of genes which interacted with the node. Larger nodes have more interactions with the genes. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)

Fig. 1

Fig. 2 A pie chart representing the number of BVDV genes interaction with host genes. In each section, the name of the virus gene and the number of interactions with the host genes are represented.

Fig. 2

Table 2 Topological characters of BVDV-host network.

Table 2Number of nodes	279	
Number of edges	436	
Avg. number of neighbors	3.118	
Network diameter	8	
Characteristics path length	3.506	
Clustering coefficient	0.012	
Network density	0.006	

3.2 Screening of hub genes and modules analysis

Different ranked algorithms in CytoHubba were used to identify host hub genes. The ranks of top ten hub genes selected by each feature are shown in Table 3. CD46 and EEF2 genes were identified by all five algorithms, also TXN2 was detected by four of five ranked methods in cytoHubba.Table 3 Top ten hub genes detected using different topological algorithms in cytoHubba.

Table 3Algorithm	MCC	DMNC	MNC	Radiality	Stress	
1	CD46a	CD46a	CD46a	CD46a	CD46a	
2	EEF2a	EEF2a	EEF2a	TXN2b	TXN2b	
3	TXN2b	ACTB	TXN2b	EEF2a	EEF2a	
4	ACTB	HAX1	TNEM70	RPLPO	RPLPO	
5	HAX1	PCBP1	CNDP2	RPL7A	RPL7A	
6	PCBP1	EZR	PNMA1	RPL24	RPL24	
7	EZR	ACTG1	ACTG1	RPL3	RPL3	
8	ACTG1	RACK1	TMED 10	RPS3	RPS3	
9	RACK1	GALK1	AIP	ACTB	ELOC	
10	GALK1	G1K1R6	G1K1R6	HAX1	UPF1	
MCC: Maximal Clique Centrality; DMNC: Density of Maximum Neighborhood Component; EPC: Edge Percolated Component.

a The hub genes identified by all algorithms in cytoHubba.

b The hub genes identified by four ranked algorithms in cytoHubba.

Different sub-networks were constructed using various clustering algorithms to find important protein complexes and functionally modules in BVDV-host protein interaction network. Fig. 3 shows multiple clusters made for the whole network by Molecular Complex Detection (MCODE), CytoCluster, ClusterViz, and Supervised Complex Detection (SCODE). Two sub-cluster were constructed by MCODE. Sub-cluster 1 had 6 nodes in which 4 nodes belonged to viral proteins and 2 belonged to host (Fig. 3A). Sub-cluster 2 consisted of 4 nodes, two of which belonged to host genes (Fig. 3B). In the cluster made by CytoCluster, there were 3 nodes among which NS2 belonged to viral proteins while CD46 and DNAJC14 were bovine proteins (Fig. 3C). According to cluster analysis performed by ClusterVis, four smaller sub-networks were constructed. First sub-network (Fig. 3D) contained 106 nodes among which NS2, NS3, NS2-3, NS4B, and NS5B were viral proteins and the others belonged to host proteins. Second sub-network consisted of 92 nodes, in which Npro was viral protein (Fig. 3E). Third sub-network had 63 nodes among which E2 was viral protein (Fig. 3F). Forth sub-network was consisted of viral protein NS5A and five host proteins (Fig. 3G). Cluster constructed by SCODE had 3 nodes of Erns and E2 (viral proteins), and ATG14 (host protein).Fig. 3 Key moddules for the whole network by MCODE, CytoCluster, ClusterViz, and SCODE. (A and B) Sub-clusters constructed using the MCODE plugin. (C) Sub-cluster identified using CytoCluster. (D, E, F, and G) Sub-clusters constructed using the ClusterVis plugin. Virus genes are showon in red, while host genes are presented in purple. (H) Modules constructed using SCODE. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.)

Fig. 3

3.3 KEGG and GO enrichment analysis

Kyoto Encyclopedia of Genes and Genomics (KEGG) pathway analysis was performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) database for all host proteins interacted with BVDV and specifically top ten hub genes ranked by CytoHubba based on Maximal Clique Centrality (MCC) algorithm. Table 4 shows the KEGG pathway analysis for all BVDV-interacting host proteins. The interacting proteins were highly enriched in ribosome, spliceosome, Influenza A, protein processing in endoplasmic reticulum, proteoglycans in cancer, measles, RNA transport, and protein export while KEGG pathway analysis of top 10 hub genes revealed that these genes were involved in leukocyte transendothelial migration. The Gene Ontology (GO) function analysis was carried out using AgriGO. The results showed that the proteins were significantly involved in biological processes such as metabolic (98 genes), cellular (122 genes), localization (28 genes), response to stimulus (21 genes), biosynthetic (60 genes), gene expression (63 genes), translation (43 genes), multicellular organismal (23 genes), developmental (19 genes), etc processes (Fig. 4A). For GO cellular component analysis, the proteins were mainly associated with cytoplasm (116 genes), organelle (112 genes), nucleus part (41 genes), ribosome (34 genes), cytoskeleton (20 genes), mitochondrion (14 genes), and endoplasmic reticulum (12 genes) (Fig. 4B). The proteins significantly enriched in molecular functions including binding (125 genes), protein binding (91 genes), structural molecule activity (38 genes), structural constitute of ribosome (32 genes), RNA binding (32 genes), translation factor activity (5 genes), etc (Fig. 4C). The GO cellular component (CC), and molecular function (MF) of the top hub genes ranked by MCC algorithm are shown in Supplementary Fig. 1. The GO biological process (BP) of hub genes was not significant (P > 0.05).Table 4 Kyoto encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of host genes interacted with BVDV genes.

Table 4Term	Count (%)	P-Value	Benjamini	
Ribosome	35 (14.5)	8.7E-30	1.4E-27	
Spliceosome	13 (5.4)	5.7E-6	4.7E-4	
Influenza A	10 (4.1)	5.7E-3	2.4E-1	
Protein processing in endoplasmic reticulum	9 (3.7)	1.5E-2	5.1E-1	
Proteoglycans in cancer	8 (3.3)	9.3E-2	1.0E0	
Measles	7 (2.9)	5.1E-2	1.0E0	
RNA transport	7 (2.9)	8.3E-2	1.0E0	
Protein export	5 (2.1)	1.0E-3	5.5E-2	

Fig. 4 The Gene Ontology (GO) enrichment analysis using AgriGO. (A) GO terms biological process (BP). (B) GO terms cellular compenets (CC). (C) GO terms molecular functions (MF). GO term codes are shown in braces.

Fig. 4

3.4 PANTHER classification of proteins

The host proteins that interact with BVDV were classified using the Protein Analysis Through Evolutionary Relationships (PANTHER) Classification System. According to the protein classes of host targeted genes, translational proteins, nucleic acid metabolism proteins, metabolite interconversion enzymes, and protein modifying enzymes were the most likely candidates as targets of viral proteins (Fig. 5). Some of these proteins involve in various biological categories as cellular proteins often play several roles in several pathways.Fig. 5 Panther pie-chart of host interacting proteins classification. The three main GO terms were subclassified to child categories (pie charts at the bottom of the image). GO term codes are shown in braces.

Fig. 5

4 Discussion

We constructed a comprehensive BVDV-host protein interaction map to describe the interactome between BVDV and host proteins, aiming to better understand the pathogenicity of this virus and identify potential drug targets for the treating of BVDV infection. Although BVDV, like other RNA viruses, produces few proteins, these proteins often interact with multiple host proteins to create favorable conditions for the virus's survival and replication [14]. Identifying of host proteins targeted by viral proteins and investigating under both of their functions in normal and infectious conditions can help us to uncover the key processes required for infection [14]. In this context, network biology can be a useful tool for studying host-pathogen relationships at the molecular level. Various researchers have used network biology to discover key viral proteins and their targets in infection processes [6,7,9]. In the present study, all previous small-scale, large-scale, and computational investigations on BVDV-host interactions were integrated to detect highly BVDV-associated host proteins and their roles in viral pathogenesis.

In cellular network analysis, proteins with a wide range of connecting partners are called “hubs” [54]. After constructing the PPI network, 10 hub genes were identified using five different ranking methods. The CD46, EEF-2, and TXN genes were detected by at least four out of the five ranking methods in cytoHubba. CD46, the most prominent identified hub gene, is a ubiquitously expressed “multitasker,” found on all nucleated cells, protecting them from autologous damage by the complement system [55,56]. CD46 regulates both adaptive immunity and the complement system [55]. It has been identified as interacting with other pathogens, such as like the measles virus [55], human herpesvirus 6 [56], group B adenoviruses [57], and pathogenic Neisseria [58]. TXN-2, also known as Trx-2, is involved in various cellular processes, including the regulation of ROS production, apoptosis, cell growth, and the modulation of transcription factors [59]. Researchers have reported that TXN-2 plays a role in the life cycle of the classical swine fever virus (CSFV) and BVDV [59,60]. TXN-2 inhibits CSFV replication, and in the case of BVDV, it is downregulated by the virus during early infection. EEF-2 encodes eukaryotic elongation factor 2, which is essential for protein synthesis [61].

KEGG pathway analysis revealed that BVDV proteins primarily target biological pathways involved in translation and protein synthesis, including ribosome, spliceosome, protein processing in the endoplasmic reticulum, and RNA transport. BVDV interacts with ribosomal complexes to ensure that viral RNAs are translated and can produce substantial amounts of viral proteins even under the strict conditions imposed by the host shutoff mechanism, where cellular metabolism is severely limited [14,62]. As found in the present study and reported by other researchers, BVDV manipulates biological processes such as cellular metabolism, gene expression, translation, cell development and differentiation [[62], [63], [64]]. The cellular proteins and pathways targeted by BVDV proteins, including RNA binding, protein binding, spliceosome, translation, and biosynthesis, are also involved in the pathogenesis of other viruses, indicating common strategies among these viruses [54,63,65]. According to KEGG analysis, the enriched gene set is found to be involved in other infectious pathways such as Influenza A virus and measles infections.

GO enrichment analysis revealed that most interacting proteins function in various biological processes, including cellular protein metabolism, cellular macromolecular biosynthesis, gene expression, and translation. These proteins are primarily found in the nucleus, ribosome, cytoskeleton, mitochondrion, and endoplasmic reticulum. Most of the host-interacting proteins have binding molecular functions, such as protein binding, RNA binding, translation factor activity-nucleic acid binding, and cytoskeletal protein binding. Some of these proteins also serve as structural constituents of the ribosome. The PANTHER, GO, and KEGG evidence highlights the widespread effects of BVDV on cellular processes, particularly in the translation and synthesis of proteins, as reported in other studies [62,63].

5 Conclusion

Although BVDV encodes few proteins (like other RNA viruses), it can effectively affect multiple host pathways. The BVDV-host PPI network map, constructed based on small-scale, large-scale, and computational studies, revealed that viral proteins target various protein classes, including translational proteins, nucleic acid metabolism proteins, metabolite interconversion enzymes, and protein-modifying enzymes. CD46, EEF-2, and TXN were identified as hub genes using different ranking algorithms. The present study shows that BVDV interactions with its host primarily focus on targeting translation, protein synthesis, and cellular metabolic pathways. Further studies are needed to explore the effects of these proteins on cellular pathways and the molecular mechanisms of BVDV infection.

Funding

No funding.

Compliance with ethical standards

Not applicable.

Availability of data and materials

The datasets used and/or analyzed during the current study are available through the paper and additional files.

CRediT authorship contribution statement

Seyedeh Elham Rezatofighi: Writing – original draft, Visualization, Validation, Supervision, Project administration, Methodology, Investigation, Formal analysis, Data curation.

Declaration of competing interest

The authors declare no conflict of interests.

Appendix A Supplementary data

The following are the Supplementary data to this article:Multimedia component 1

Multimedia component 1

Multimedia component 2

Multimedia component 2

figs1 figs1

figs2 figs2

Acknowledgements

The authors are very thankful to 10.13039/501100005412 Shahid Chamran University of Ahvaz for the facilities to accomplish the present research project (Grant Number: SCU.SB98.658 ). We also thank Dr. Mahdi Askari Badouei for editing the article.

Appendix A Supplementary data to this article can be found online at https://doi.org/10.1016/j.bbrep.2024.101825.
==== Refs
References

1 Walker P.J. Siddell S.G. Lefkowitz E.J. Mushegian A.R. Adriaenssens E.M. Dempsey D.M. Dutilh B.E. Harrach B. Harrison R.L. Hendrickson R.C. Changes to virus taxonomy and the statutes ratified by the international committee on taxonomy of viruses Arch. Virol. 165 2020 2737 2748 10.1007/s00705-023-05797-4 32816125
2 Weiskircher E. Aligo J. Ning G. Konan K.V. Bovine viral diarrhea virus NS4B protein is an integral membrane protein associated with Golgi markers and rearranged host membranes Virol. J. 6 2009 185 10.1186/1743-422X-6-185 19887001
3 Ridpath J. Bolin S. Dubovi E. Segregation of bovine viral diarrhea virus into genotypes Virology 205 1994 66 74 10.1006/viro.1994.1620 7975238
4 Yamane D. Zahoor M.A. Mohamed Y.M. Azab W. Kato K. Tohya Y. Akashi H. Inhibition of sphingosine kinase by bovine viral diarrhea virus NS3 is crucial for efficient viral replication and cytopathogenesis J. Biol. Chem. 284 2009 13648 13659 10.1074/jbc.M807498200 19293152
5 Halehalli R.R. Nagarajaram H.A. Molecular principles of human virus protein-protein interactions Bioinformatics 31 2012 1025 1033 10.1093/bioinformatics/btu763
6 Farooq Q.U.A. Shaukat Z. Zhou T. Aiman S. Gong W. Li C. Inferring Virus-Host relationship between HPV and its host Homo sapiens using protein interaction network Sci. Rep. 10 2020 8719 10.1038/s41598-020-65837-w 32457456
7 Qua Farooq Khan F.F. Construction and analysis of a comprehensive protein interaction network of HCV with its host Homo sapiens BMC Infect. Dis. 19 2019 367 10.1186/s12879-019-4000-9 31039741
8 Wei W. Kong W. Identification of key genes and signaling pathways during Sendai virus infection in vitro Braz. J. Microbiol. 50 2019 13 22 10.1007/s42770-018-0021-6 30637656
9 Farooq Q.U.A. Shaukat Z. Aiman S. Zhou T. Li C. A systems biology-driven approach to construct a comprehensive protein interaction network of influenza A virus with its host BMC Infect. Dis. 20 2020 1 10 10.1186/s12879-020-05214-0
10 Kumar S. Lata K.S. Sharma P. Bhairappanavar S.B. Soni S. Das J. Inferring pathogen-host interactions between Leptospira interrogans and Homo sapiens using network theory Sci. Rep. 9 2019 1 17 10.1038/s41598-018-38329-1 30626917
11 Garamszegi S. Franzosa E.A. Xia Y. Signatures of pleiotropy, economy and convergent evolution in a domain-resolved map of human-virus protein-protein interaction networks PLoS Pathog. 9 2013 e1003778 10.1371/journal.ppat.1003778
12 Pichlmair A. Kandasamy K. Alvisi G. Mulhern O. Sacco R. Habjan M. Binder M. Stefanovic A. Eberle C.A. Goncalves A. Viral immune modulators perturb the human molecular network by common and unique strategies Nature 487 2012 486 490 10.1038/nature11289 22810585
13 Meyniel-Schicklin L. de Chassey B. André P. Lotteau V. Viruses and interactomes in translation Mol. Cell. Proteomics 11 2012 014738 10.1074/mcp.M111.014738 M111
14 Lazar C. Zitzmann N. Dwek R.A. Branza-Nichita N. The pestivirus E(rns) glycoprotein interacts with E2 in both infected cells and mature virions Virology 314 2003 696 705 10.1016/s0042-6822(03)00510-5 14554096
15 Hilton L. Moganeradj K. Zhang G. Chen Y.H. Randall R.E. McCauley J.W. Goodbourn S. The NPro product of bovine viral diarrhea virus inhibits DNA binding by interferon regulatory factor 3 and targets it for proteasomal degradation J. Virol. 80 2006 11723 11732 10.1128/JVI.01145-06 16971436
16 Zahoor M.A. Yamane D. Mohamed Y.M. Suda Y. Kobayashi K. Kato K. Tohya Y. Akashi H. Bovine viral diarrhea virus non-structural protein 5A interacts with NIK- and IKKbeta-binding protein J. Gen. Virol. 91 2010 1939 1948 10.1099/vir.0.020990-0 20444997
17 Gong X. Chen Q. Zheng F. Identification of protein inhibitor of activated STAT 4, a novel host interacting partner that involved in bovine viral diarrhea virus growth Virol. J. 17 2020 59 10.1186/s12985-020-01330-0 32321515
18 Johnson C.M. Perez D.R. French R. Merrick W.C. Donis R.O. The NS5A protein of bovine viral diarrhoea virus interacts with the alpha subunit of translation elongation factor-1 J. Gen. Virol. 82 2001 2935 2943 10.1099/0022-1317-82-12-2935 11714969
19 Leveringhaus E. Cagatay G.N. Hardt J. Becher P. Postel A. Different impact of bovine complement regulatory protein 46 (CD46(bov)) as a cellular receptor for members of the species Pestivirus H and Pestivirus G Emerg. Microb. Infect. 2012 1 38 10.1080/22221751.2021.2011620
20 Mu Y. Tews B.A. Luttermann C. Meyers G. Interaction of pestiviral E1 and E2 sequences in dimer formation and intracellular retention Int. J. Mol. Sci. 22 2021 2 22 10.3390/ijms22147285
21 Branza-Nichita N. Lazar C. Dwek R.A. Zitzmann N. Role of N-glycan trimming in the folding and secretion of the pestivirus protein E(rns) Biochem. Biophys. Res. Commun. 319 2004 655 662 10.1016/j.bbrc.2004.05.039 15178456
22 Fu Q. Shi H. Shi M. Meng L. Bao H. Zhang G. Ren Y. Zhang H. Guo F. Qiao J. Roles of bovine viral diarrhea virus envelope glycoproteins in inducing autophagy in MDBK cells Microb. Pathog. 76 2014 61 66 10.1016/j.micpath.2014.09.011 25245007
23 Darweesh M.F. Rajput M.K. Braun L.J. Rohila J.S. Chase C.C. BVDV Npro protein mediates the BVDV induced immunosuppression through interaction with cellular S100A9 protein Microb. Pathog. 121 2018 341 349 10.1016/j.micpath.2018.05.047 29859294
24 Schaut R.G. McGill J.L. Neill J.D. Ridpath J.F. Sacco R.E. Bovine viral diarrhea virus type 2 in vivo infection modulates TLR4 responsiveness in differentiated myeloid cells which is associated with decreased MyD88 expression Virus Res. 208 2015 44 55 10.1016/j.virusres.2015.05.017 26043978
25 Xu J. Mendez E. Caron P.R. Lin C. Murcko M.A. Collett M.S. Rice C.M. Bovine viral diarrhea virus NS3 serine proteinase: polyprotein cleavage sites, cofactor requirements, and molecular model of an enzyme essential for pestivirus replication J. Virol. 71 1997 5312 5322 10.1128/JVI.71.7.5312-5322.1997 9188600
26 Branza-Nichita N. Durantel D. Carrouée-Durantel S. Dwek R.A. Zitzmann N. Antiviral effect of N-butyldeoxynojirimycin against bovine viral diarrhea virus correlates with misfolding of E2 envelope proteins and impairment of their association into E1-E2 heterodimers J. Virol. 75 2001 3527 3536 10.1128/JVI.75.8.3527-3536.2001 11264342
27 Maurer K. Krey T. Moennig V. Thiel H-Jr Rümenapf T. CD46 is a cellular receptor for bovine viral diarrhea virus J. Virol. 78 2004 1792 1799 10.1128/jvi.78.4.1792-1799.2004 14747544
28 Lackner T. Muller A. Konig M. Thiel H.-J. Tautz N. Persistence of bovine viral diarrhea virus is determined by a cellular cofactor of a viral autoprotease J. Virol. 79 2005 9746 9755 10.1128/JVI.79.15.9746-9755.2005 16014936
29 Chen Z. Rijnbrand R. Jangra R.K. Devaraj S.G. Qu L. Ma Y. Lemon S.M. Li K. Ubiquitination and proteasomal degradation of interferon regulatory factor-3 induced by Npro from a cytopathic bovine viral diarrhea virus Virology 366 2007 277 292 10.1016/j.virol.2007.04.023 17531282
30 Dubrau D. Tortorici M.A. Rey F.A. Tautz N. A positive-strand RNA virus uses alternative protein-protein interactions within a viral protease/cofactor complex to switch between RNA replication and virion morphogenesis PLoS Pathog. 13 2017 e1006134 10.1371/journal.ppat.1006134
31 Mohamed Y.M. Bangphoomi N. Yamane D. Suda Y. Kato K. Horimoto T. Akashi H. Physical interaction between bovine viral diarrhea virus nonstructural protein 4A and adenosine deaminase acting on RNA (ADAR) Arch. Virol. 159 2014 1735 1741 10.1007/s00705-014-1997-3 24500065
32 Wang J. Liu B. Wang N. Lee Y.M. Liu C. Li K. TRIM56 is a virus- and interferon-inducible E3 ubiquitin ligase that restricts pestivirus infection J. Virol. 85 2011 3733 3745 10.1128/JVI.02546-10 21289118
33 Jefferson M. Donaszi-Ivanov A. Pollen S. Dalmay T. Saalbach G. Powell P.P. Host factors that interact with the pestivirus N-terminal protease, Npro, are components of the ribonucleoprotein complex J. Virol. 88 2014 10340 10353 10.1128/JVI.00984-14 24965446
34 Ammari M.G. Pharr G.T. Pendarvis K. Pinchuk L.M. McCarthy F.M. Evaluation of bovine viral diarrhea virus interactions with its host proteome during the course of infection Biomed Res Clin Prac 2 2017 1 7 10.15761/BRCP.1000140
35 Gladue D.P. Baker-Bransetter R. Holinka L.G. Fernandez-Sainz I.J. O'Donnell V. Fletcher P. Lu Z. Borca M.V. Interaction of CSFV E2 protein with swine host factors as detected by yeast two-hybrid system PLoS One 9 2014 e85324 10.1371/journal.pone.0085324
36 Cook H.V. Doncheva N.T. Szklarczyk D. Von Mering C. Jensen L.J. Viruses. STRING: a virus-host protein-protein interaction database Viruses 10 2018 519 10.3390/v10100519 30249048
37 Chatr-Aryamontri A. Ceol A. Peluso D. Nardozza A. Panni S. Sacco F. Tinti M. Smolyar A. Castagnoli L. Vidal M. VirusMINT: a viral protein interaction database Nucleic Acids Res. 37 2009 D669 D673 10.1093/nar/gkn739 18974184
38 Hermjakob H. Montecchi‐Palazzi L. Lewington C. Mudali S. Kerrien S. Orchard S. Vingron M. Roechert B. Roepstorff P. Valencia A. IntAct: an open source molecular interaction database Nucleic Acids Res. 32 2004 D452 D455 10.1093/nar/gkh052 14681455
39 Guirimand T. Delmotte S. Navratil V. VirHostNet 2.0: surfing on the web of virus/host molecular interactions data Nucleic Acids Res. 43 2015 D583 D587 10.1093/nar/gku1121 25392406
40 Calderone A. Licata L. Cesareni G. VirusMentha: a new resource for virus-host protein interactions Nucleic Acids Res. 43 2015 D588 D592 10.1093/nar/gku830 25217587
41 Brito A.F. Pinney J.W. Protein–protein interactions in virus–host systems Front. Microbiol. 8 2017 1557 10.3389/fmicb.2017.01557 28861068
42 Franzosa E.A. Xia Y. Structural principles within the human-virus protein-protein interaction network Proc. Natl. Acad. Sci. USA 108 2011 10538 10543 10.1073/pnas.1101440108 21680884
43 Consortium T.U. UniProt: the universal protein knowledgebase in 2021 Nucleic Acids Res. 49 2020 D480 D489 10.1093/nar/gkaa1100
44 Otasek D. Morris J.H. Bouças J. Pico A.R. Demchak B. Cytoscape Automation: empowering workflow-based network analysis Genome Biol. 20 2019 185 10.1186/s13059-019-1758-4 31477170
45 Chin C.-H. Chen S.-H. Wu H.-H. Ho C.-W. Ko M.-T. Lin C.-Y. cytoHubba: identifying hub objects and sub-networks from complex interactome BMC Syst. Biol. 8 2014 1 7 10.1186/1752-0509-8-S4-S11 24393148
46 Bader G.D. Hogue C.W. An automated method for finding molecular complexes in large protein interaction networks BMC Bioinf. 4 2003 1 27 10.1186/1471-2105-4-2
47 Mohamed S. Janus N. Qi Y. SCODE: a Cytoscape app for supervised complex detection in protein-protein interaction graphs F1000Research 5 2016 1699 10.12688/f1000research.9184.1
48 Li M. Li D. Tang Y. Wu F. Wang J. CytoCluster: a cytoscape plugin for cluster analysis and visualization of biological networks Int. J. Mol. Sci. 18 2017 1880 10.3390/ijms18091880 28858211
49 Wang J. Zhong J. Chen G. Li M. Wu F-x Pan Y. ClusterViz: a cytoscape APP for cluster analysis of biological network IEEE ACM Trans. Comput. Biol. Bioinf 12 2014 815 822 10.1109/TCBB.2014.2361348
50 Hu L. Zhao T. Sun Y. Chen Y. Bai K. Xu F. Bioinformatic identification of hub genes and key pathways in neutrophils of patients with acute respiratory distress syndrome Medicine 99 2020 e19820 10.1097/MD.0000000000019820
51 Dennis G. Sherman B.T. Hosack D.A. Yang J. Gao W. Lane H.C. Lempicki R.A. DAVID: database for annotation, visualization, and integrated discovery Genome Biol. 4 2003 R60 PMID: 12734009
52 Tian T. Liu Y. Yan H. You Q. Yi X. Du Z. Xu W. Su Z. agriGO v2.0: a GO analysis toolkit for the agricultural community, 2017 update Nucleic Acids Res. 45 2017 W122 W129 10.1093/nar/gkx382 28472432
53 Mi H. Ebert D. Muruganujan A. Mills C. Albou L.-P. Mushayamaha T. Thomas P.D. PANTHER version 16: a revised family classification, tree-based classification tool, enhancer regions and extensive API Nucleic Acids Res. 49 2020 D394 D403 10.1093/nar/gkaa1106
54 De Chassey B. Navratil V. Tafforeau L. Hiet M. Aublin‐Gex A. Agaugue S. Meiffren G. Pradezynski F. Faria B. Chantier T. Hepatitis C virus infection protein network Mol. Syst. Biol. 4 2008 230 10.1038/msb.2008.66 18985028
55 Naniche D. Varior-Krishnan G. Cervoni F. Wild T. Rossi B. Rabourdin-Combe Ca Gerlier D. Human membrane cofactor protein (CD46) acts as a cellular receptor for measles virus J. Virol. 67 1993 6025 6032 10.1128/JVI.67.10.6025-6032.1993 8371352
56 Santoro F. Kennedy P.E. Locatelli G. Malnati M.S. Berger E.A. Lusso P. CD46 is a cellular receptor for human herpesvirus 6 Cell 99 1999 817 827 10.1016/s0092-8674(00)81678-5 10619434
57 Gaggar A. Shayakhmetov D.M. Lieber A. CD46 is a cellular receptor for group B adenoviruses Nat. Med. 9 2003 1408 1412 10.1038/nm952 14566335
58 Cattaneo R. Four viruses, two bacteria, and one receptor: membrane cofactor protein (CD46) as pathogens' magnet J. Virol. 78 2004 4385 4388 10.1128/jvi.78.9.4385-4388.2004 15078919
59 Li S. Wang J. He W.R. Feng S. Li Y. Wang X. Liao Y. Qin H.Y. Li L.F. Dong H. Thioredoxin 2 is a novel E2-interacting protein that inhibits the replication of classical swine fever virus J. Virol. 89 2015 8510 8524 10.1128/JVI.00429-15 26041303
60 Villalba M. Fredericksen F. Otth C. Olavarría V. Transcriptomic analysis of responses to cytopathic bovine viral diarrhea virus-1 (BVDV-1) infection in MDBK cells Mol. Immunol. 71 2016 192 202 10.1016/j.molimm.2016.01.009 26919728
61 Jørgensen R. Merrill A.R. Andersen G.R. The life and death of translation elongation factor 2 Biochem. Soc. Trans. 34 2006 1 6 10.1042/BST20060001 16246167
62 Zhang A. He L. Wang Y. Prediction of GCRV virus-host protein interactome based on structural motif-domain interactions BMC Bioinf. 18 2017 145 10.1186/s12859-017-1500-8
63 Goodwin C.M. Xu S. Munger J. Stealing the keys to the kitchen: viral manipulation of the host cell metabolic network Trends Microbiol. 23 2015 789 798 10.1016/j.tim.2015.08.007 26439298
64 Sanchez E.L. Lagunoff M. Viral activation of cellular metabolism Virology 479 2015 609 618 10.1016/j.virol.2015.02.038 25812764
65 Wang L. Zhou L. Zhang H. Li Y. Ge X. Guo X. Yu K. Yang H. Interactome profile of the host cellular proteins and the nonstructural protein 2 of porcine reproductive and respiratory syndrome virus PLoS One 9 2014 e99176 10.1371/journal.pone.0099176
