Topological analysis of protein co-abundance networks identifies novel host targets important for HCV infection and pathogenesis
© McDermott et al; licensee BioMed Central Ltd. 2012
Received: 13 March 2012
Accepted: 30 April 2012
Published: 30 April 2012
High-throughput methods for obtaining global measurements of transcript and protein levels in biological samples has provided a large amount of data for identification of 'target' genes and proteins of interest. These targets may be mediators of functional processes involved in disease and therefore represent key points of control for viruses and bacterial pathogens. Genes and proteins that are the most highly differentially regulated are generally considered to be the most important. We present topological analysis of co-abundance networks as an alternative to differential regulation for confident identification of target proteins from two related global proteomics studies of hepatitis C virus (HCV) infection.
We analyzed global proteomics data sets from a cell culture study of HCV infection and from a clinical study of liver biopsies from HCV-positive patients. Using lists of proteins known to be interaction partners with pathogen proteins we show that the most differentially regulated proteins in both data sets are indeed enriched in pathogen interactors. We then use these data sets to generate co-abundance networks that link proteins based on similar abundance patterns in time or across patients. Analysis of these co-abundance networks using a variety of network topology measures revealed that both degree and betweenness could be used to identify pathogen interactors with better accuracy than differential regulation alone, though betweenness provides the best discrimination. We found that though overall differential regulation was not correlated between the cell culture and liver biopsy data, network topology was conserved to an extent. Finally, we identified a set of proteins that has high betweenness topology in both networks including a protein that we have recently shown to be essential for HCV replication in cell culture.
The results presented show that the network topology of protein co-abundance networks can be used to identify proteins important for viral replication. These proteins represent targets for further experimental investigation that will provide biological insight and potentially could be exploited for novel therapeutic approaches to combat HCV infection.
Recent advances in high-throughput methods for taking global measurements of transcript or protein levels from biological samples have driven the field of systems biology. A common application of such methods is to identify genes or proteins that are likely to be involved in the disease process being studied to direct further experimental investigation. These 'targets' are potential mediators of important aspects of the disease, or may be downstream responses to the disease process. Targets are generally identified from the most highly differentially expressed genes or proteins. However, this approach can overlook genes or proteins that are important, but may not be the most highly differentially regulated, such as transcription factors or other upstream mediators of critical processes . In this study we extend our previous work showing that targets can be identified using network approaches based on global proteomics measurements . We show that the differential regulation of a protein is an important factor in predicting biological significance, but that treating the data as a network and using topological measures allows for better prediction of biologically significant targets, provides better ranking of proteins, and allows extension by integrating other kinds of relationships, for example protein-protein interactions. Additionally we show that network topology of proteins is more conserved between experiments than is differential regulation. Our work provides a framework for network analysis of global proteomics data, and shows that this approach can identify biologically interesting targets.
Hepatitis C virus (HCV), a single-stranded positive RNA virus in the Flaviviridae family, is a major cause of liver disease in chronically infected individuals. Chronic infection causes inflammation and fibrosis of the liver and increases the chance of developing more serious hepatocellular carcinoma or cirrhosis in approximately 30% of infected individuals . Current therapies have limited efficacy and numerous side effects  and a major challenge in translational hepatology research is the development of new approaches that target critical processes in the HCV life cycle and progression to disease state. Currently, study of HCV infection has been carried out in cell culture , on liver biopsy samples from infected patients , and in limited animal models , however similarities and differences between these different systems have not been extensively studied.
Previously we used global proteomics and lipidomics to show that HCV can reprogram cellular metabolism and bioenergetics in cell culture . In order to identify possible targets through which HCV regulates metabolic reprogramming we constructed a correlation network based on global proteomics measurements of human hepatoma Huh7.5 cells responding to a time course of HCV infection . We used the topology of the network, specifically proteins with high betweenness or bottlenecks, to identify biologically important proteins. Subsequently we showed that genetic silencing and pharmacological inhibition of one of these predicted targets, DCI, significantly inhibited processes critical for HCV infection . These results showed the utility of network approaches to identify key components and interactions associated with HCV infection in cell culture experiments, but did not delineate how the approaches could be applied to provide the best results, nor if the approach would generalize to other proteomics data sets with very different experimental designs.
While our previous studies used network analysis to identify targets for further experimental investigation, they did not explore the generality and robustness of the approach. Though promising, the approach requires analyses of the parameters used for network generation and target identification, analysis of topological measures beyond betweenness centrality, and application to other similar data sets. Only by exploring these aspects can the significance and applicability of the approach be established. The current study had two principal aims. The first was to evaluate these factors for network-based target identification from proteomics data and to compare this approach with an existing method for identification of important proteins from global proteomics data, differential regulation. This is particularly important work because proteomics technology has recently reached a point where it is possible to generate studies with multiple global proteomics datasets of the system being studied under different conditions and there have been very few reports describing use of proteomics data in network inference approaches. The second aim was to compare network-based analysis of proteomics from HCV infection in cell culture experiments with similar networks generated from liver biopsy samples to identify common targets that have potential translational impact. These aims represent an important and significant advance over our previous work because we systematically compare our network topology approach with traditional approaches to target identification, characterize the impact of network inference parameters on our results, and compare the results obtained in our cell culture studies with those obtained from clinically relevant patient-derived samples.
In order to further explore the identification of novel, translationally relevant pathways and important proteins involved in HCV infection and liver disease progression, we first analyzed the topology of networks inferred from the time course study of HCV infection in cultured Huh7.5 cells with an emphasis on now evaluating the ability of various topological measures to predict proteins known to be targeted by pathogens in general and HCV proteins specifically. As described above, we further integrated protein-protein interaction data in the networks and showed that the integrated networks provide improved discrimination of important proteins using network topology. An important observation from this analysis was that network topology provided better discrimination of important proteins than differential regulation. We then reanalyzed proteomics data from a previous cross-sectional study of HCV infected patients  using the same approaches. We obtained similar trends in this analysis for identification of important proteins in vivo. In addition we found a number of proteins that share important topological roles in networks inferred from both the in vitro system and the in vivo clinical samples. We conclude that considering proteomic data as networks highlight important in vivo proteins from examination of in vitro systems thus, providing valuable insight into translationally relevant disease processes.
We used two datasets in this study that have been previously described. The first is from the Huh7.5 human hepatoma cell line infected with a chimeric HCV genotype 2a virus, J6/JFH-1 . Cells were inoculated with HCV or UV-inactivated virus and samples taken at 24, 48, 72, and 96 hours post-infection. The samples were analyzed by liquid chromatography-mass spectrometry (LC-MS) using the accurate mass-and-time tag (AMT) approach in combination with trypsin-catalyzed 16O/18O labeling for quantitation [8, 9]. Briefly, peptides from time-matched mocks were individually labeled with 18O and spiked at equal amounts into the appropriate HCV or UV-HCV-inoculated sample. The corresponding 18O/16O intensity data from multiple observations of the same protein were then rolled up to compute a final protein abundance ratio for all proteins identified in a given sample and, to identify those proteins exhibiting statistically significant (p < 0.05) changes in abundance compared to the control sample .
The second dataset used was from HCV-infected liver tissues from 15 patients at different stages of fibrosis . This study also employed stable isotope 16O/18O trypsin catalyzed labeling in combination with the AMT tag approach for protein quantitation. In this case, proteins exhibiting statistically significant (by ANOVA on fibrosis stage groups ; p < 0.05) changes in abundance were determined relative to a control sample consisting of peptides generated from a pool of 8 HCV-positive patients with minimal liver disease as previously described .
The list of proteins known to be physically targeted by pathogens (interactors) was taken from supplemental material in . A list of proteins identified in a two-hybrid screen as interacting with HCV proteins was obtained from supplemental material in . Mouse homologs were obtained from the Mouse Genome Database (MGI) . A list of human genes that exhibit positive selection was obtained from the Human PAML Browser  available at http://mendel.gene.cwru.edu/adamslab/cgi-bin/paml/pbrowser.py using a significance threshold of p < 0.01. Protein-protein interactions were obtained from http://cytoscape.wodaklab.org/wiki/Data_Sets.
Proteomics data filtering and network construction
Proteomics data was filtered for significance.
Data was converted to a ratio versus control conditions.
A filter was applied to remove differential abundance ratios below a threshold (abundance filter).
Correlation values were calculated between present values for all pairs of proteins.
A filter was applied to remove correlation values with a number of comparisons below a threshold (correspondence filter).
A filter was applied to remove correlation values below a threshold (correlation filter).
Significance filtering and ratio calculation are described above and in the original papers [2, 5]. The abundance filter (step 3) replaces all values with ratios below the threshold with missing values in the vector of abundance ratios for each protein. Correlation is calculated as the Pearson correlation coefficient for all pairwise complete observations (steps 4). Correlation values with a low number of comparisons are removed (set to 0) according to the correspondence filter, where a single comparison is counted if abundance ratios are observed for the same condition for the pair of proteins being considered. Finally, the correlation filter is used to generate a final adjacency matrix, which is then treated as a network for topological analysis. Previously, the impact of the choice of similarity threshold on construction of coexpression networks has been investigated [14, 15]. However, in this study we have chosen reasonable values for these parameters by evaluating the topological enrichment of the resulting network in proteins known to be targeted by pathogens (see Results). For topological analysis (below) we varied parameters for the three filters listed here to generate multiple different networks.
Topological analysis of networks was performed using in-house scripts in the statistical language R http://www.r-project.org/ that utilize the igraph R library http://igraph.sourceforge.net/. We provide our code in (Additional file 1). Advanced topological analysis was performed using the network analysis software UCINET 6.0 http://www.analytictech.com/ucinet/. Examples of advanced topology metrics are reachability , Katz influence , and Bonacich power centrality . The clustering coefficient, and degree, closeness and betweenness centrality metrics are defined as below [19, 20].
Generally, degree centrality is the fraction of edges for a particular protein out of all possible interactions for that protein in the network.
Generally, closeness is the mean distance between a protein and all other proteins in the network.
Generally, betweenness is the number of shortest paths between all pairs of proteins in the network that pass through a specific node.
then normalize the centrality by the largest value and repeat until the values converge. Generally, a protein with high eigenvector centrality is connected to other proteins who themselves are connected to many other proteins.
Generally, page rank is the importance of a protein in the network.
Generally, the clustering coefficient is defined as the percentage of neighboring proteins that interact with each other
Functional enrichment analysis
Enrichment of a population (for example, the top 20% of proteins in terms of betweenness) for a particular functional label was calculated using the hypergeometric test. Functional labels were defined by the list of pathogen or HCV targets, or positively selected genes. In all cases the background for significance was the total list of proteins determined to be significantly differentially regulated not including the population in question. Significance levels are indicated in the text but in general a p-value of 0.05 or below was considered to be significant. Where indicated, multiple hypothesis correction was applied to p-values using the Bonferroni correction.
Results and discussion
Highly abundant proteins are more likely to be targeted by pathogens
Analysis of patient samples should provide a more direct assessment of the validity of potential clinical targets, compared to in vitro experimental models. We examined data from a previous study analyzing liver biopsy samples from 15 patients at five stages of fibrosis  to compare our observations in cell culture. We compared proteins from the 15 infected patients against a pool of 8 HCV-positive control patient samples. Significance was assessed using ANOVA resulting in 210 significantly changing proteins, and 193 of these proteins were associated with a gene symbol and could be used for enrichment calculations with the pathogen interactor lists. Our results are presented in Figure 1B and show that the most differentially regulated proteins (top 20%) are enriched in proteins that are pathogen interactors in fibrosis stages 1 and 4, however, these differences are not statistically significant due to the smaller number of significant proteins identified.
Network topology identification of HCV targeted proteins
Networks inferred from proteomics show significant enrichment in pathogen targets across many network inference parameters
CoA + PPI
Following our previous approach we combined each network with experimentally determined protein-protein interactions (PPIs) between observed proteins. In this process known PPIs between proteins already in the co-abundance network are added as new edges to the network. The results of this analysis are shown in Figure 3B. Open symbols show the enrichment in the PPI network alone. These results show that the inferred protein association relationships can improve target discrimination using each topological measure except clustering coefficient, but that the best discrimination is provided by degree followed by betweenness. The specificity and sensitivity of both the degree and betweenness approaches are significantly better than either differential regulation or the inferred networks without PPIs. Table 1 provides a summary enrichment across all networks (full results in Additional file 3: Table S2). These results show that all four topological measures provide highly significant (p-value < 1e-16) enrichment in pathogen targets, with betweenness and clustering coefficient displaying highest enrichment, thereby demonstrating the added value of incorporating PPI data into inferred networks for a generalizable approach to identify target regulatory nodes within networks.
It is well known that some topological measures display varying degrees of overlap; for example, proteins may have both high betweenness and high degree and thus be bottleneck-hubs . We were interested in assessing the relationships between topological measures in our integrated network using Spearman rank correlation (to account for differences in the distributions of these measures). The results are presented in Additional file 4: Figure S1 and show that degree and closeness are highly correlated in all of the networks examined, while betweenness was slightly less correlated with these two and clustering coefficient was the least correlated with the other measures. Eigenvector centrality was highly correlated with degree and closeness, but pagerank was not as correlated with the other measures.
Examining degree and betweenness, the most used metrics for biological networks, we found that of 343 bottlenecks and hubs (the top 20% of proteins as ranked by betweenness and degree, respectively), 184 (53%) were shared, reflecting the moderate correlation between degree and betweenness and the fact that they aren't capturing the same characteristics of the networks. To examine this overlap further, we assessed whether the enrichment of bottlenecks in pathogen targets was dependent on their hub status within the integrated cell network. In Figure 4B we show the results from a topological subgroup analysis (similar to ) showing enrichment in pathogen targets for various overlapping groups. Interestingly, these results show that betweenness alone contributes more to importance than does degree, since the hub and hub-nonbottleneck groups are less enriched than the bottleneck, hub-bottleneck, or nonhub-bottleneck groups. Similar results were observed in the previous study by Yu, et al.  for regulatory networks, but not for PPI networks, indicating that our inferred networks combined with PPIs maintain the properties of regulatory networks and are less similar to PPI networks.
Functional characterization of topologically-defined targets
Given their enrichment in pathogen interactors, we hypothesized that proteins with high betweenness might be involved in similar functions. We therefore investigated this in the network with the best enrichment from the analysis above. This was a network that was given by an abundance filter of 0, a correspondence filter of 4, and a correlation filter of 0.9 (see Methods). We then assessed the top 20% of the proteins in this network ranked by betweenness for functional enrichment in gene ontology categories. Despite the fact that the bottleneck proteins from this network were the most enriched for pathogen interactors as well as for HCV-specific interactors, we found no significant enrichment in any functional categories, relative to the other proteins in the network. This indicates that these proteins are united by their importance to the replication of HCV, but span diverse functional categories.
Application of network analysis to clinical proteomics data
Network topology is more conserved than differential regulation
The Huh7.5 cell culture system has been extensively used as a model for HCV infection [2, 24–26]. However, it is unclear at the molecular level how common patterns of expression and regulation might be in terms of HCV pathogenesis and disease progression. To examine this we used two approaches: examining the correspondence in differential regulation between the two datasets and examining the correspondence of topological characteristics in the networks described. For network comparison we chose to compare the networks with the best topological enrichment of pathogen targets, as described for the cell culture network above. For the liver biopsy derived network we chose a correspondence filter of 5, abundance filter of 2, correlation filter of 0.8, with integrated PPIs (see Additional file 3: Table S2). The filters used were different than those used for the Huh7.5-derived network because the structure of the underlying data sets were different. The liver biopsy data contains data from 15 patients and thus the number of corresponding data points used is more (5 versus 4) and the abundance filter is related to the overall range of differential abundance so the difference here (2 versus 0) reflects the larger variance observed in the patient samples. These differences highlight the fact that some care must be used when applying these methods to different kinds of datasets.
We first compared differential regulation in the 148 proteins that were observed in both the cell culture samples and the liver biopsy samples. Differential regulation ratios for infected samples were averaged per protein across all time points or patients. This process provides a reasonable estimate of the overall differential regulation for a protein in each experiment. We then compared the abundance ratios for proteins identified in both datasets using Spearman rank correlation. The two lists displayed no correlation (R = -0.05) indicating that the overall level of differential regulation in the cell culture system is not a good indicator of differential regulation in the liver biopsy samples. To ensure that this result did not reflect the use of averaged differential abundance across patients and time points that could mask true correlation between the two groups, we also calculated the Spearman rank correlation between the differential abundance ratios for all fibrosis stages from patients versus all time points. These results confirmed our findings; the mean correlation in differential abundance between different stages of fibrosis was 0.63, and between different time points was 0.24, whereas the mean correlation between the two sets was 0.03. The maximum correlation between any two fibrosis stages was 0.73, and between any two time points was 0.74, whereas the maximum correlation between any fibrosis stage and time point was 0.23. These results show that differential abundance, in general, is not well conserved between proteins in the cell system and patient samples. This analysis is consistent with our previous analysis  showing that there was a subset of proteins displaying similar temporal progression in the cell culture and liver biopsies since the current analysis compared the regulation of all proteins, and did not incorporate the temporal information. We discuss these findings further in the Conclusion section.
We next examined the agreement between topological measures between the cell-culture derived network and the network derived from liver biopsy samples. The Spearman rank correlation comparing betweenness measures in proteins in both networks was 0.4. Though not perfect agreement, this correlation is much better than the correlation obtained comparing differential regulation. To examine this in a slightly different way, we examined the distribution of betweenness values from the clinical network in topological bottlenecks, using a two-sided t-test. We found that bottlenecks in the cell culture network have a significantly higher mean betweenness (from clinical network) than non-bottleneck proteins (138 versus 42, p-value 0.02). These results indicate that the general topology of the networks is more conserved than is differential regulation of the individual proteins in each dataset, supporting our notion that network topology provides information not provided by differential regulation in some cases.
Network topology provides better target identification than differential regulation
We postulated that network topology could provide better discrimination of interesting proteins than differential regulation. To address this we compared the enrichment of known pathogen interactors, which we consider to be interesting targets, in the most highly differentially regulated proteins from each dataset (see Figure 1) and the proteins with highest betweenness in the network derived from each dataset. This analysis revealed that the top differentially regulated proteins (top 20%) from any time point in the cell culture data set were comprised of 20% known pathogen interactors (relative to 14% in the remaining portion of proteins). The proteins with the top betweenness (bottlenecks; top 20%) from the cell culture-derived network were comprised of 32% known pathogen targets (relative to 10% background). This observation is consistent with results in the liver biopsy data where the maximum enrichment in the top differentially regulated proteins was 17% versus an enrichment of 33% for top ranked bottlenecks. These numbers were also reflected in HCV-specific targets (data not shown). Collectively, these results show that network topology provides better identification of target proteins than does simply ranking by differential regulation.
Conserved bottlenecks between cell culture and clinical samples
glutathione S-transferase kappa 1
inner membrane protein, mitochondrial (mitofilin)
dodecenoyl-Coenzyme A delta isomerase
ATP synthase, H + transporting
tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation
heat shock 70 kDa protein 8
P, H, E
ribosomal protein, large, P1
cytochrome c, somatic
heat shock 70 kDa protein 9 (mortalin)
glutamic-oxaloacetic transaminase 2
To further investigate potential functional roles of these conserved bottleneck proteins we performed functional enrichment on the first-order networks of each protein. The first-order networks for several conserved bottlenecks are shown in (Additional file 8: Figure S2) and the functional categories significantly enriched in each neighborhood are listed in (Additional file 9: Table S6). Though many of the individual conserved bottleneck proteins were linked to mitochondria (Table 2), the functions of their neighborhoods are fairly diverse. However, two neighborhood networks were significantly enriched in processes related to fatty acid metabolism and its regulation (DCI and YWHAQ). Additionally, we provide the topology of the conserved targets in both networks in (Additional file 10: Table S7). These results show that many of the bottlenecks are also hubs (highly connected proteins) in both networks including CALR, ETFA, IMMT, and RPLP1. Interestingly, all of the targets have low clustering coefficients. The clustering coefficient is the percentage of a node's neighbors that are also linked to each other and reflects the density of edges in that portion of the network. Given that betweenness is a primary driver of importance in the network (Figure 4B) this is not a surprising observation. That is, even hubs having many neighbors may be playing connecting roles in the network because they are also bottlenecks, and a high density of edges in their neighborhoods would decrease their betweenness since this would provide multiple routes through their neighborhood.
We previously described network analyses of cell culture data to define interactions between host and pathogen and identified mitochondrial fatty acid oxidation enzymes that are predicted to function as central points for connecting and controlling metabolic pathways and as such, key targets in HCV-associated metabolic reprogramming . In fact, dys-regulations in mitochondrial function are evidenced by wide-spread perturbation of related proteins across every HCV model system we have studied [2, 5, 7, 26, 27]. Thus, the modeling efforts reported here leveraged these data to further investigate whether the parameters used in our prior in vitro modeling activities (abundance filter of 0, correspondence filter of 4, and correlation filter of 0.9) were the best to help identify new targets. Indeed, our previous study did not examine whether the use of a network topology approach could identify important targets any better than a standard approach such as considering highly differentially regulated proteins. Additionally, though we previously showed that this approach was valuable in cell culture studies it also remained unclear how it would perform on data from very different kinds of samples, such as those from liver biopsies of HCV-positive patients.
In the current study we build upon our previous findings to determine if there are other topological metrics (for example clustering coefficient and closeness) that identify targets in these networks, to define the parameters for network construction that provide the best target identification, and to characterize the relationship between networks derived from the cell culture data and those derived from data from patient biopsies. Our results indicate that the method of network construction has a significant impact on the results obtained. We found that betweenness was the most effective metric for defining important targets in our network but that other topological metrics (degree, clustering coefficient and closeness) could also discriminate targets to a statistically significant extent. From previous work examining the properties of topological bottlenecks in networks inferred from global transcriptomics data we have postulated that bottlenecks may represent mediators of transitions between states of the system [1, 28–30], and therefore represent critical points of control for the disease process. We have speculated that this is because bottlenecks link functional modules that represent groups of genes or proteins coexpressed under similar conditions. The transition between modules may represent state changes in the system, and the position of bottlenecks makes them candidates for regulators of these transitions. Our finding that degree was also a good predictor of importance in the system reiterates previous findings in other undirected biological networks , though the primary contribution to importance we found to be betweenness, similar to findings in regulatory networks. Our findings here are consistent with the idea that bottlenecks in coabundance networks represent transitions between functional modules, and show that bottlenecks and hubs from proteomics-based networks may have similar properties as those from transcriptomics-based networks.
We note that modeling activities involving integrated genomic-proteomic analyses is an important area of research aimed at understanding the differences between co-expression at the transcript and protein level. However, our initial modeling efforts centered on the utilization of proteomic and metabolomic data indicating a temporal regulation of cellular metabolic homeostasis that was not detected by the accompanying gene expression profiles. Indeed our prior in vitro studies were unique in part because they described a previously un-identified role for post-transcriptional regulatory mechanisms in the metabolic rerouting that was observed . For this reason, the scope of the current manuscript has focused on extending our analyses specifically to comparison with in vivo protein co-expression networks.
Upon optimization of network construction, subsequent comparative analyses revealed that topologically-defined bottleneck proteins in the cell culture-derived network were generally more differentially regulated in patients with advanced fibrosis than their non-bottleneck counterparts. Interestingly, this was not observed when comparing differential abundance alone between the two datasets, indicating that topological analysis may identify more clinically relevant targets from cell culture studies than relative expression. It is important to note that we previously identified a subset of proteins that showed strongly conserved patterns of differential abundance  between the cell culture and liver biopsy samples. In the current analysis we show that as an overall measure, differential abundance does not correlate well between the two data sets. Additionally, bottlenecks in the cell culture network were more likely to be bottlenecks in the clinical network. This shows that our approach can identify proteins of interest based on cell culture studies that are important in human disease and that these proteins would not be identified by examining differential abundance alone. Importantly, these findings point to the limitations of identifying/prioritizing pathogen-host targets based solely on highly differential regulation, a common approach to the identification of targets for further investigation.
Throughout this study we refer to target proteins as proteins that are important for HCV replication and/or fibrosis development. Some of these proteins have been defined using two-hybrid screens [11, 21], but our working hypothesis is that there are proteins that are important for replication that have not been previously defined. These are proteins that may or may not be direct interaction partners with HCV proteins but could contribute to metabolic or signaling pathways necessary for HCV replication and/or liver disease progression. We previously proposed an important role for temporal regulation of mitochondrial fatty acid oxidation and energy production in HCV infection and liver disease progression. Briefly, we described early increases in mitochondrial fatty acid oxidation that contribute to the creation of a "pro-viral" environment immediately preceding the subsequent increase in viral replication observed in vitro . This was eventually followed by a decline in fatty acid oxidation that accompanied the appearance of a cytopathic effect in vitro and liver disease progression in vivo . The down-regulation of mitochondrial fatty acid oxidation would favor an increase in hepatocellular lipid content (for example, steatosis), a common occurrence in HCV, and histological feature observed among 4 of 6 patients with advanced fibrosis in our in vivo studies . The conservation of protein abundance changes associated with pathogenesis in vitro (e.g. cytopathic effect) and liver disease progression in vivo, and the corresponding mitochondrial bottlenecks reported here, including DCI, raises the interesting prospect that these proteins play an important role in the viral life cycle and pathogenesis.
Our previous findings and those described in the current study prompted us to further explore the predicted influence of HCV-associated disruptions in mitochondrial fatty acid oxidation, including consideration of whether these perturbations would be reflected by disease-related patterns detected in blood. From a clinical perspective, biomarker discovery efforts in body fluids represent an attractive alternative to tissue samples owing to the relative ease and less invasive nature of collection and the large volumes that normally can be obtained. We have observed the accumulation of both substrates for enoyl-CoA isomerase activity (e.g. DCI) as well as dicarboxylic acids well known to reflect alternative fatty acid catabolism through ω-oxidation pathways, findings consistent with our predictions regarding an important role for DCI, the essential link between saturated and unsaturated β-oxidation, in the impaired mitochondrial fatty acid catabolism occurring during HCV-associated liver disease progression . Thus, the identification of disease-related fatty acid patterns in the blood of patients with HCV-associated liver disease progression provides a potentially useful noninvasive diagnostic link to the previously described alterations in hepatic mitochondrial fatty acid oxidation occurring during HCV infection and pathogenesis. Importantly, we have unequivocally validated a biologically relevant role for DCI in the HCV life cycle using a combination of gene silencing and pharmacologic inhibition approaches [2, 5, 7]. In summary, our data from multiple model systems and clinically relevant physiologic compartments provide evidence confirming our original modeling predictions regarding a requirement for DCI in the HCV life cycle  and demonstrate a physiologically relevant association of temporal declines in fatty acid oxidation that coincide with pathogenesis in vitro and in vivo. Taken together, we believe these data provide proof of principle for the utility of integrated in vitro/in vivo modeling efforts to identify key host targets of HCV infection and pathogenesis.
The biological interpretation of the remaining top 10% bottlenecks, 4 out of 5 of which are mitochondrial proteins with links to fatty acid oxidation and energy production, was predicated on the wealth of data described for the representative example DCI as highlighted above together with the growing literature on the important role of altered mitochondrial function in HCV infection and pathogenesis (for an excellent review on the interactions between HCV and mitochondria we recommend ). Among the additional bottlenecks identified was glutathione-S-transferase kappa 1 (GSTK1), a protein that localizes to the mitochondria and peroxisome and has pleiotropic functions including glutathione conjugation, peroxidase, and disulphide-bond-forming oxidoreductase activities . Interestingly, GSTK1 has recently been shown to play an important role in the oligomeric assembly and secretion of adiponectin, a cytokine that stimulates fatty acid oxidation through interaction with the hepatic receptor AdipoR2 and subsequent activation of peroxisome proliferator-activated receptor (PPAR)-alpha [31, 32]. HCV-associated targeting of GSTK1 and DCI may serve to provide multiple control points for modulating catabolic flux of fatty acids during metabolic reprogramming. GSTK1 may promote further cross-talk between metabolic signaling and biochemical pathways by modulating the folding and assembly of oligomeric proteins directly involved in lipid synthesis and/or catabolism, including the trimeric DCI protein. A similar role in the folding of lipid metabolism enzymes has been suggested in Caenorhabditis elegans where GSTK1 silencing was associated with a decline in the biosynthesis of the monounsaturated fatty acid cis-vaccenic acid . It is worth noting that the differential abundance of cis-vaccenic acid was observed to impact lipid droplet remodeling under pathogenic conditions of defective peroxisomal β-oxidation in C. elegans . Taken together, these findings suggest interesting new avenues of research aimed at exploring the interplay between GSTK1 and DCI during metabolic reprogramming and the lipid remodeling events predicted to provide important constituents in the various structural entities supporting the HCV life cycle, including the lipid droplet and membranous replicase compartments.
Among the other bottlenecks detected in our analyses was mitofilin, also known as mitochondrial inner membrane protein (IMMT). Mitofilin is a protein localized to the inner mitochondrial membrane whose presence is essential for tubular cristae formation and the increased surface-to-volume ratio of the inner membrane that occurs during increased metabolic output . While the molecular basis for these alterations in mitochondrial cristae morphology are not well understood, mitofilin depletion has been shown to induce aberrant structural changes in the inner membrane that are associated with abrogation of ATP production despite increased flux of fatty acid substrates through the β-oxidation pathway thus, suggesting an adverse impact on the oxidative phosphorylation machinery that resides in the inner membrane . We suspect that the putative HCV targeting of mitofilin reflects a coordinated effort to maximize energy production in support of the significant macromolecular biosynthesis necessary for viral growth . Consistent with this idea we further identified ATP5B, the major catalytic subunit of F1 ATP synthase, as a conserved bottleneck in our studies. A similarly important pro-viral role for ATP5B has recently been reported for herpes simplex virus-1 (HSV-1) . In a series of elegant experiments aimed at exploring the effect of host microRNAs on HSV-1 replication, Zheng et al, identified a point of cross talk between host cell and virus that results in the progressive induction of host cell miR-101 levels that is accompanied by concomitant declines in ATP5B expression and HSV-1 replication . The interplay between virus and the miR-101/ATP5B regulatory network suggests a potential link between modulation of this host defense mechanism and the establishment of long-term HSV-1 latency . This latter point is of particular interest as we and others have proposed a similar role for modulation of fatty acid oxidation and energy production in the establishment of persistent HCV and measles virus infection [2, 37].
It is important to note that our intent is not to provide a network representation that is faithful to the underlying true network of interactions in the cell, but rather to use topology in these simply defined association networks to identify target proteins for further experimental investigation. The networks generated using this approach are based on correlation of protein abundance over many different observations (time points in the cell culture data and patients in the clinical data). As such they represent the information flow in the system. For example, closely coordinated proteins are close together in the networks, while those with little or no coordination are far apart. It is likely that this organization allows use of topology to query the network for more important proteins, since bottlenecks in particular represent points constriction in information flow in the system . In a fashion analogous to that for DCI, additional conserved bottleneck proteins represent particularly attractive targets for further investigation of their functional significance during HCV infection and liver disease progression. In this regard, recent efforts to link these findings with clinical protein profiling studies of serial liver biopsies obtained from HCV-positive liver transplant recipients revealed a statistically significant up-regulation of the protein bottleneck GSTK1 in patients who developed severe liver injury . Importantly, the increased abundance of GSTK1 occurred prior to histologic evidence of fibrosis. Collectively, these findings merit further investigation to understand the functional, regulatory and/or prognostic significance of this protein bottleneck during HCV-associated liver disease progression.
In summary, the results presented in this study show that a network approach to consideration of global proteomics data is a powerful way to identify important target proteins and to elucidate potential mechanisms of pathogenesis. Previous results in yeast [23, 38], fruit fly and worm , pathogenic bacteria [40, 41], cyanobacteria , mouse macrophages , mouse blood  and human cell culture [2, 7] support the notion that our approach is generally applicable, though these have been focused on analysis of coexpression networks from transcriptomics. We have recently published on the network analysis of proteomics data from Salmonella under infectious-like conditions, and have found that these networks show a similar kind of enrichment of bottlenecks in proteins important to the system . In the current work we fully characterize the application of this approach to protein co-abundance networks showing that it works very well to identify important nodes in the network. In this study we show that topological betweenness provides the best identification of important target proteins, but that other topological measures can also be used to identify targets. Importantly, we show that this approach can be applied successfully to global proteomic data derived from liver biopsies of HCV-positive fibrosis patients. Key findings of the study were validated in a patient cohort by metabolic profiling in serum . Interestingly, the topology of cell culture networks provides better insight into important proteins in the liver biopsy data than does differential regulation, showing that it is a viable alternative or complement to standard analysis methods. Our approach represents a generally applicable method for using global proteomics data as a systems biology tool that goes beyond differential abundance of individual proteins. The finding that other metrics could also identify targets suggests that combining network metrics in some way may provide improved discrimination over the individual measures. Our initial results using a simple mean, geometric mean, or minimum of protein rank from each of the four metrics revealed that the results were not improved (data not shown). We are currently investigating more sophisticated methods for integrating multiple topological measures to improve our results.
This work was supported by the National Institute on Drug Abuse grant 1P30DA01562501 to M.G.K. Portions of this research were also supported by the NIH National Center for Research Resources (RR18522 to RDS). Portions of the research were performed at the W.R. Wiley Environmental Molecular Sciences Laboratory, a national scientific user facility sponsored by US Department of Energy's Office of Biological and Environmental Research (BER) program located at PNNL. PNNL is operated for the US Department of Energy by Battelle under contract DE-AC05-76RLO-1830.
- McDermott JE, Costa M, Janszen D, Singhal M, Tilton SC: Separating the drivers from the driven: Integrative network and pathway approaches aid identification of disease biomarkers from high-throughput data. Dis Markers. 2010, 28 (4): 253-266.View ArticleGoogle Scholar
- Diamond DL, Syder AJ, Jacobs JM, Sorensen CM, Walters KA, Proll SC, McDermott JE, Gritsenko MA, Zhang Q, Zhao R, et al: Temporal proteome and lipidome profiles reveal hepatitis C virus-associated reprogramming of hepatocellular metabolism and bioenergetics. PLoS pathogens. 2010, 6 (1): e1000719-10.1371/journal.ppat.1000719.View ArticleGoogle Scholar
- Alter MJ, Margolis HS, Krawczynski K, Judson FN, Mares A, Alexander WJ, Hu PY, Miller JK, Gerber MA, Sampliner RE, et al: The natural history of community-acquired hepatitis C in the United States. The Sentinel Counties Chronic non-A, non-B Hepatitis Study Team. N Engl J Med. 1992, 327 (27): 1899-1905. 10.1056/NEJM199212313272702.View ArticleGoogle Scholar
- Ikeda M, Kato N: Modulation of host metabolism as a target of new antivirals. Adv Drug Deliv Rev. 2007, 59 (12): 1277-1289. 10.1016/j.addr.2007.03.021.View ArticleGoogle Scholar
- Diamond DL, Jacobs JM, Paeper B, Proll SC, Gritsenko MA, Carithers RL, Larson AM, Yeh MM, Camp DG, Smith RD, et al: Proteomic profiling of human liver biopsies: hepatitis C virus-induced fibrosis and mitochondrial dysfunction. Hepatology. 2007, 46 (3): 649-657. 10.1002/hep.21751.View ArticleGoogle Scholar
- Dorner M, Horwitz JA, Robbins JB, Barry WT, Feng Q, Mu K, Jones CT, Schoggins JW, Catanese MT, Burton DR, et al: A genetically humanized mouse model for hepatitis C virus infection. Nature. 2011, 474 (7350): 208-211. 10.1038/nature10168.View ArticleGoogle Scholar
- Rasmussen A, Diamond D, McDermott J, Metz T, Gao X, Matzke M, Carter V, Belisle S, Korth M, Waters K, et al: Systems Virology Identifies a Mitochondrial Fatty Acid Oxidation Enyzme, Dodecenoyl-CoA Delta Isomerase (DCI), Required for HCV Replication and Pathogen. J Virol. 2011Google Scholar
- Smith RD, Anderson GA, Lipton MS, Pasa-Tolic L, Shen Y, Conrads TP, Veenstra TD, Udseth HR: An accurate mass tag strategy for quantitative and high-throughput proteome measurements. Proteomics. 2002, 2 (5): 513-523. 10.1002/1615-9861(200205)2:5<513::AID-PROT513>3.0.CO;2-W.View ArticleGoogle Scholar
- Qian WJ, Monroe ME, Liu T, Jacobs JM, Anderson GA, Shen Y, Moore RJ, Anderson DJ, Zhang R, Calvano SE, et al: Quantitative proteome analysis of human plasma following in vivo lipopolysaccharide administration using 16O/18O labeling and the accurate mass and time tag approach. Mol Cell Proteomics. 2005, 4 (5): 700-709. 10.1074/mcp.M500045-MCP200.View ArticleGoogle Scholar
- Dyer MD, Murali TM, Sobral BW: The landscape of human proteins interacting with viruses and other pathogens. PLoS Pathog. 2008, 4 (2): e32-10.1371/journal.ppat.0040032.View ArticleGoogle Scholar
- de Chassey B, Navratil V, Tafforeau L, Hiet MS, Aublin-Gex A, Agaugue S, Meiffren G, Pradezynski F, Faria BF, Chantier T, et al: Hepatitis C virus infection protein network. Mol Syst Biol. 2008, 4: 230-View ArticleGoogle Scholar
- Eppig JT, Bult CJ, Kadin JA, Richardson JE, Blake JA, Anagnostopoulos A, Baldarelli RM, Baya M, Beal JS, Bello SM, et al: The Mouse Genome Database (MGD): from genes to mice--a community resource for mouse biology. Nucleic Acids Res. 2005, D471-D475. 33 DatabaseGoogle Scholar
- Nickel GC, Tefft D, Adams MD: Human PAML browser: a database of positive selection on human genes using phylogenetic methods. Nucleic Acids Res. 2008, D800-D808. 36 DatabaseGoogle Scholar
- Butte AJ, Kohane IS: Mutual information relevance networks: functional genomic clustering using pairwise entropy measurements. Pac Symp Biocomput. 2000, 5: 418-429.Google Scholar
- Borate BR, Chesler EJ, Langston MA, Saxton AM, Voy BH: Comparison of threshold selection methods for microarray gene co-expression matrices. BMC Res Notes. 2009, 2: 240-10.1186/1756-0500-2-240.View ArticleGoogle Scholar
- Higley J, Hoffman-Lange U, Kadushin C, Moore G: Elite integration in stable democracies: a reconsideration. European Sociological Review. 1991, 7: 35-53.Google Scholar
- Katz L: A new index derived from sociometric data analysis. Psychometrika. 1953, 18: 39-43. 10.1007/BF02289026.View ArticleGoogle Scholar
- Bonacich P: Power and centrality: a family of measures. American Journal of Sociology. 1987, 92: 1170-1182. 10.1086/228631.View ArticleGoogle Scholar
- Wasserman S, Faust K: Social Network Analysis. 1994, Cambridge: Cambridge University PressView ArticleGoogle Scholar
- Watts DJ, Strogatz SH: Collective dynamics of 'small-world' networks. Nature. 1998, 393 (6684): 440-442. 10.1038/30918.View ArticleGoogle Scholar
- Tripathi LP, Kataoka C, Taguwa S, Moriishi K, Mori Y, Matsuura Y, Mizuguchi K: Network based analysis of hepatitis C virus core and NS4B protein interactions. Mol Biosyst. 2010, 6 (12): 2539-2553. 10.1039/c0mb00103a.View ArticleGoogle Scholar
- Brin S, Page L: The Anatomy of a Large-Scale Hypertextual Web Search Engine. 7th World-Wide Web Conference: 1998; Brisbane, Australia. 1998Google Scholar
- Yu H, Kim PM, Sprecher E, Trifonov V, Gerstein M: The importance of bottlenecks in protein networks: correlation with gene essentiality and expression dynamics. PLoS Comput Biol. 2007, 3 (4): e59-10.1371/journal.pcbi.0030059.View ArticleGoogle Scholar
- Cai Z, Zhang C, Chang KS, Jiang J, Ahn BC, Wakita T, Liang TJ, Luo G: Robust production of infectious hepatitis C virus (HCV) from stably HCV cDNA-transfected human hepatoma cells. Journal of virology. 2005, 79 (22): 13963-13973. 10.1128/JVI.79.22.13963-13973.2005.View ArticleGoogle Scholar
- Lindenbach BD, Evans MJ, Syder AJ, Wolk B, Tellinghuisen TL, Liu CC, Maruyama T, Hynes RO, Burton DR, McKeating JA, et al: Complete replication of hepatitis C virus in cell culture. Science. 2005, 309 (5734): 623-626. 10.1126/science.1114016.View ArticleGoogle Scholar
- Walters KA, Syder AJ, Lederer SL, Diamond DL, Paeper B, Rice CM, Katze MG: Genomic analysis reveals a potential role for cell cycle perturbation in HCV-mediated apoptosis of cultured hepatocytes. PLoS Pathog. 2009, 5 (1): e1000269-10.1371/journal.ppat.1000269.View ArticleGoogle Scholar
- Jacobs JM, Diamond DL, Chan EY, Gritsenko MA, Qian W, Stastna M, Baas T, Camp DG, Carithers RL, Smith RD, et al: Proteome analysis of liver cells expressing a full-length hepatitis C virus (HCV) replicon and biopsy specimens of posttransplantation liver from HCV-infected patients. Journal of virology. 2005, 79 (12): 7558-7569. 10.1128/JVI.79.12.7558-7569.2005.View ArticleGoogle Scholar
- Diamond DL, Krasnoselsky AL, Burnum KE, Monroe ME, Webb-Robertson BJ, McDermott JE, Yeh MM, Golib Dzib JF, Susnow N, Strom S, et al: Proteome and computational analyses reveal new insights into the mechanisms of hepatitis C virus mediated liver disease post-transplantation. Hepatology. 2012, doi:10.1002/hep.25649Google Scholar
- Piccoli C, Quarato G, Ripoli M, D'Aprile A, Scrima R, Cela O, Boffoli D, Moradpour D, Capitanio N: HCV infection induces mitochondrial bioenergetic unbalance: causes and effects. Biochim Biophys Acta. 2009, 1787 (5): 539-546. 10.1016/j.bbabio.2008.11.008.View ArticleGoogle Scholar
- Shield AJ, Murray TP, Cappello JY, Coggan M, Board PG: Polymorphisms in the human glutathione transferase Kappa (GSTK1) promoter alter gene expression. Genomics. 2010, 95 (5): 299-305. 10.1016/j.ygeno.2010.02.007.View ArticleGoogle Scholar
- Bertolani C, Marra F: Role of adipocytokines in hepatic fibrosis. Curr Pharm Des. 2010, 16 (17): 1929-1940. 10.2174/138161210791208857.View ArticleGoogle Scholar
- Liu M, Zhou L, Xu A, Lam KS, Wetzel MD, Xiang R, Zhang J, Xin X, Dong LQ, Liu F: A disulfide-bond A oxidoreductase-like protein (DsbA-L) regulates adiponectin multimerization. Proceedings of the National Academy of Sciences of the United States of America. 2008, 105 (47): 18302-18307. 10.1073/pnas.0806341105.View ArticleGoogle Scholar
- Petit E, Michelet X, Rauch C, Bertrand-Michel J, Terce F, Legouis R, Morel F: Glutathione transferases kappa 1 and kappa 2 localize in peroxisomes and mitochondria, respectively, and are involved in lipid metabolism and respiration in Caenorhabditis elegans. FEBS J. 2009, 276 (18): 5030-5040. 10.1111/j.1742-4658.2009.07200.x.View ArticleGoogle Scholar
- Zhang SO, Box AC, Xu N, Le Men J, Yu J, Guo F, Trimble R, Mak HY: Genetic and dietary regulation of lipid droplet expansion in Caenorhabditis elegans. Proceedings of the National Academy of Sciences of the United States of America. 2010, 107 (10): 4640-4645. 10.1073/pnas.0912308107.View ArticleGoogle Scholar
- John GB, Shang Y, Li L, Renken C, Mannella CA, Selker JM, Rangell L, Bennett MJ, Zha J: The mitochondrial inner membrane protein mitofilin controls cristae morphology. Mol Biol Cell. 2005, 16 (3): 1543-1554. 10.1091/mbc.E04-08-0697.View ArticleGoogle Scholar
- Zheng SQ, Li YX, Zhang Y, Li X, Tang H: MiR-101 regulates HSV-1 replication by targeting ATP5B. Antiviral Res. 2011, 89 (3): 219-226. 10.1016/j.antiviral.2011.01.008.View ArticleGoogle Scholar
- Takahashi M, Watari E, Shinya E, Shimizu T, Takahashi H: Suppression of virus replication via down-modulation of mitochondrial short chain enoyl-CoA hydratase in human glioblastoma cells. Antiviral Res. 2007, 75 (2): 152-158. 10.1016/j.antiviral.2007.02.002.View ArticleGoogle Scholar
- Jeong H, Mason SP, Barabasi AL, Oltvai ZN: Lethality and centrality in protein networks. Nature. 2001, 411 (6833): 41-42. 10.1038/35075138.View ArticleGoogle Scholar
- Hahn MW, Kern AD: Comparative genomics of centrality and essentiality in three eukaryotic protein-interaction networks. Mol Biol Evol. 2005, 22 (4): 803-806. 10.1093/molbev/msi072.View ArticleGoogle Scholar
- McDermott JE, Taylor RC, Yoon H, Heffron F: Bottlenecks and hubs in inferred networks are important for virulence in Salmonella typhimurium. J Comput Biol. 2009, 16 (2): 169-180. 10.1089/cmb.2008.04TT.View ArticleGoogle Scholar
- Yoon H, Ansong C, McDermott JE, Gritsenko M, Smith RD, Heffron F, Adkins JN: Systems analysis of multiple regulator perturbations allows discovery of virulence factors in Salmonella. BMC systems biology. 2011, 5: 100-10.1186/1752-0509-5-100.View ArticleGoogle Scholar
- McDermott JE, Oehmen CS, McCue LA, Hill E, Choi DM, Stockel J, Liberton M, Pakrasi HB, Sherman LA: A model of cyclic transcriptomic behavior in the cyanobacterium Cyanothece sp. ATCC 51142. Mol Biosyst. 2011, 7: 2407-2418. 10.1039/c1mb05006k.View ArticleGoogle Scholar
- McDermott JE, Archuleta M, Thrall BD, Adkins JN, Waters KM: Controlling the response: predictive modeling of a highly central, pathogen-targeted core response module in macrophage activation. PLoS ONE. 2011, 6 (2): e14673-10.1371/journal.pone.0014673.View ArticleGoogle Scholar
- McDermott JE, Archuleta M, Stevens SL, Stenzel-Poore MP, Sanfilippo A: Defining the players in higher-order networks: predictive modeling for reverse engineering functional influence networks. Pac Symp Biocomput. 2011, 16: 314-325.Google Scholar
- Niemann GS, Brown RN, Gustin JK, Stufkens A, Shaikh-Kidwai AS, Li J, McDermott JE, Brewer HM, Schepmoes A, Smith RD, et al: Discovery of novel secreted virulence factors from Salmonella enterica serovar Typhimurium by proteomic analysis of culture supernatants. Infect Immun. 2011, 79 (1): 33-43. 10.1128/IAI.00771-10.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.