Person:

Padi, Megha

Many human diseases, arising from mutations of disease susceptibility genes (genetic diseases), are also associated with viral infections (virally implicated diseases), either in a directly causal manner or by indirect associations. Here we examine whether viral perturbations of host interactome may underlie such virally implicated disease relationships. Using as models two different human viruses, Epstein-Barr virus (EBV) and human papillomavirus (HPV), we find that host targets of viral proteins reside in network proximity to products of disease susceptibility genes. Expression changes in virally implicated disease tissues and comorbidity patterns cluster significantly in the network vicinity of viral targets. The topological proximity found between cellular targets of viral proteins and disease genes was exploited to uncover a novel pathway linking HPV to Fanconi anemia.

Genotypic differences greatly influence susceptibility and resistance to disease. Understanding genotype-phenotype relationships requires that phenotypes be viewed as manifestations of network properties, rather than simply as the result of individual genomic variations. Genome sequencing efforts have identified numerous germline mutations associated with cancer predisposition and large numbers of somatic genomic alterations. However, it remains challenging to distinguish between background, or “passenger” and causal, or “driver” cancer mutations in these datasets. Human viruses intrinsically depend on their host cell during the course of infection and can elicit pathological phenotypes similar to those arising from mutations. To test the hypothesis that genomic variations and tumour viruses may cause cancer via related mechanisms, we systematically examined host interactome and transcriptome network perturbations caused by DNA tumour virus proteins. The resulting integrated viral perturbation data reflects rewiring of the host cell networks, and highlights pathways that go awry in cancer, such as Notch signalling and apoptosis. We show that systematic analyses of host targets of viral proteins can identify cancer genes with a success rate on par with their identification through functional genomics and large-scale cataloguing of tumour mutations. Together, these complementary approaches result in increased specificity for cancer gene identification. Combining systems-level studies of pathogen-encoded gene products with genomic approaches will facilitate prioritization of cancer-causing driver genes so as to advance understanding of the genetic basis of human cancer.

Background: Genome-wide libraries of yeast deletion strains have been used to screen for genes that drive phenotypes such as stress response. A surprising observation emerging from these studies is that the genes with the largest changes in mRNA expression during a state transition are not those that drive that transition. Here, we show that integrating gene expression data with context-independent protein interaction networks can help prioritize master regulators that drive biological phenotypes. Results: Genes essential for survival had previously been shown to exhibit high centrality in protein interaction networks. However, the set of genes that drive growth in any specific condition is highly context-dependent. We inferred regulatory networks from gene expression data and transcription factor binding motifs in Saccharomyces cerevisiae, and found that high-degree nodes in regulatory networks are enriched for transcription factors that drive the corresponding phenotypes. We then found that using a metric combining protein interaction and transcriptional networks improved the enrichment for drivers in many of the contexts we examined. We applied this principle to a dataset of gene expression in normal human fibroblasts expressing a panel of viral oncogenes. We integrated regulatory interactions inferred from this data with a database of yeast two-hybrid protein interactions and ranked 571 human transcription factors by their combined network score. The ranked list was significantly enriched in known cancer genes that could not be found by standard differential expression or enrichment analyses. Conclusions: There has been increasing recognition that network-based approaches can provide insight into critical cellular elements that help define phenotypic state. Our analysis suggests that no one network, based on a single data type, captures the full spectrum of interactions. Greater insight can instead be gained by exploring multiple independent networks and by choosing an appropriate metric on each network. Moreover we can improve our ability to rank phenotypic drivers by combining the information from individual networks. We propose that such integrative network analysis could be used to combine clinical gene expression data with interaction databases to prioritize patient- and disease-specific therapeutic targets. Electronic supplementary material The online version of this article (doi:10.1186/s12918-015-0228-1) contains supplementary material, which is available to authorized users.

Merkel cell polyomavirus (MCPyV) is an etiological agent of Merkel cell carcinoma (MCC), a highly aggressive skin cancer. The MCPyV small tumor antigen (ST) is required for maintenance of MCC and can transform normal cells. To gain insight into cellular perturbations induced by MCPyV ST, we performed transcriptome analysis of normal human fibroblasts with inducible expression of ST. MCPyV ST dynamically alters the cellular transcriptome with increased levels of glycolytic genes, including the monocarboxylate lactate transporter SLC16A1 (MCT1). Extracellular flux analysis revealed increased lactate export reflecting elevated aerobic glycolysis in ST expressing cells. Inhibition of MCT1 activity suppressed the growth of MCC cell lines and impaired MCPyV-dependent transformation of IMR90 cells. Both NF-κB and MYC have been shown to regulate MCT1 expression. While MYC was required for MCT1 induction, MCPyV-induced MCT1 levels decreased following knockdown of the NF-κB subunit RelA, supporting a synergistic activity between MCPyV and MYC in regulating MCT1 levels. Several MCC lines had high levels of MYCL and MYCN but not MYC. Increased levels of MYCL was more effective than MYC or MYCN in increasing extracellular acidification in MCC cells. Our results demonstrate the effects of MCPyV ST on the cellular transcriptome and reveal that transformation is dependent, at least in part, on elevated aerobic glycolysis.