Widespread variation in molecular interactions and regulatory properties among transcription factor isoforms

Most human Transcription factors (TFs) genes encode multiple protein isoforms differing in DNA binding domains, effector domains, or other protein regions. The global extent to which this results in functional differences between isoforms remains unknown. Here, we systematically compared 693 isoforms of 246 TF genes, assessing DNA binding, protein binding, transcriptional activation, subcellular localization, and condensate formation. Relative to reference isoforms, two-thirds of alternative TF isoforms exhibit differences in one or more molecular activities, which often could not be predicted from sequence. We observed two primary categories of alternative TF isoforms: "rewirers" and "negative regulators", both of which were associated with differentiation and cancer. Our results support a model wherein the relative expression levels of, and interactions involving, TF isoforms add an understudied layer of complexity to gene regulatory networks, demonstrating the importance of isoform-aware characterization of TF functions and providing a rich resource for further studies.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays or AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold-Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

N-terminal proteoforms may engage in different protein complexes

Alternative translation initiation and alternative splicing may give rise to N-terminal proteoforms, proteins that differ at their N-terminus compared with their canonical counterparts. Such proteoforms can have altered localizations, stabilities, and functions. Although proteoforms generated from splice variants can be engaged in different protein complexes, it remained to be studied to what extent this applies to N-terminal proteoforms. To address this, we mapped the interactomes of several pairs of N-terminal proteoforms and their canonical counterparts. First, we generated a catalogue of N-terminal proteoforms found in the HEK293T cellular cytosol from which 22 pairs were selected for interactome profiling. In addition, we provide evidence for the expression of several N-terminal proteoforms, identified in our catalogue, across different human tissues, as well as tissue-specific expression, highlighting their biological relevance. Protein-protein interaction profiling revealed that the overlap of the interactomes for both proteoforms is generally high, showing their functional relation. We also showed that N-terminal proteoforms can be engaged in new interactions and/or lose several interactions compared with their canonical counterparts, thus further expanding the functional diversity of proteomes.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays and AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

Next-generation large-scale binary protein interaction network for Drosophila melanogaster

Generating reference maps of interactome networks illuminates genetic studies by providing a protein-centric approach to finding new components of existing pathways, complexes, and processes. We apply state-of-the-art methods to identify binary protein-protein interactions (PPIs) for Drosophila melanogaster. Four all-by-all yeast two-hybrid (Y2H) screens of > 10,000 Drosophila proteins result in the 'FlyBi' dataset of 8723 PPIs among 2939 proteins. Testing subsets of data from FlyBi and previous PPI studies using an orthogonal assay allows for normalization of data quality; subsequent integration of FlyBi and previous data results in an expanded binary Drosophila reference interaction network, DroRI, comprising 17,232 interactions among 6511 proteins. We use FlyBi data to generate an autophagy network, then validate in vivo using autophagy-related assays. The deformed wings (dwg) gene encodes a protein that is both a regulator and a target of autophagy. Altogether, these resources provide a foundation for building new hypotheses regarding protein networks and function.

Assessment of community efforts to advance network-based prediction of protein-protein interactions

Comprehensive understanding of the human protein-protein interaction (PPI) network, aka the human interactome, can provide important insights into the molecular mechanisms of complex biological processes and diseases. Despite the remarkable experimental efforts undertaken to date to determine the structure of the human interactome, many PPIs remain unmapped. Computational approaches, especially network-based methods, can facilitate the identification of previously uncharacterized PPIs. Many such methods have been proposed. Yet, a systematic evaluation of existing network-based methods in predicting PPIs is still lacking. Here, we report community efforts initiated by the International Network Medicine Consortium to benchmark the ability of 26 representative network-based methods to predict PPIs across six different interactomes of four different organisms: A. thaliana, C. elegans, S. cerevisiae, and H. sapiens. Through extensive computational and experimental validations, we found that advanced similarity-based methods, which leverage the underlying network characteristics of PPIs, show superior performance over other general link prediction methods in the interactomes we considered.

A proteome-scale map of the SARS-CoV-2-human contactome

Understanding the mechanisms of coronavirus disease 2019 (COVID-19) disease severity to efficiently design therapies for emerging virus variants remains an urgent challenge of the ongoing pandemic. Infection and immune reactions are mediated by direct contacts between viral molecules and the host proteome, and the vast majority of these virus-host contacts (the 'contactome') have not been identified. Here, we present a systematic contactome map of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with the human host encompassing more than 200 binary virus-host and intraviral protein-protein interactions. We find that host proteins genetically associated with comorbidities of severe illness and long COVID are enriched in SARS-CoV-2 targeted network communities. Evaluating contactome-derived hypotheses, we demonstrate that viral NSP14 activates nuclear factor κB (NF-κB)-dependent transcription, even in the presence of cytokine signaling. Moreover, for several tested host proteins, genetic knock-down substantially reduces viral replication. Additionally, we show for USP25 that this effect is phenocopied by the small-molecule inhibitor AZ1. Our results connect viral proteins to human genetic architecture for COVID-19 severity and offer potential therapeutic targets.

Abacavir inhibits but does not cause self-reactivity to HLA-B*57:01-restricted EBV specific T cell receptors

Pre-existing pathogen-specific memory T cell responses can contribute to multiple adverse outcomes including autoimmunity and drug hypersensitivity. How the specificity of the T cell receptor (TCR) is subverted or seconded in many of these diseases remains unclear. Here, we apply abacavir hypersensitivity (AHS) as a model to address this question because the disease is linked to memory T cell responses and the HLA risk allele, HLA-B*57:01, and the initiating insult, abacavir, are known. To investigate the role of pathogen-specific TCR specificity in mediating AHS we performed a genome-wide screen for HLA-B*57:01 restricted T cell responses to Epstein-Barr virus (EBV), one of the most prevalent human pathogens. T cell epitope mapping revealed HLA-B*57:01 restricted responses to 17 EBV open reading frames and identified an epitope encoded by EBNA3C. Using these data, we cloned the dominant TCR for EBNA3C and a previously defined epitope within EBNA3B. TCR specificity to each epitope was confirmed, however, cloned TCRs did not cross-react with abacavir plus self-peptide. Nevertheless, abacavir inhibited TCR interactions with their cognate ligands, demonstrating that TCR specificity may be subverted by a drug molecule. These results provide an experimental road map for future studies addressing the heterologous immune responses of TCRs including T cell mediated adverse drug reactions.

HLA-B*57:01 restricted EBV-specific T-cell receptor specificity is altered by abacavir, suggesting a potentially inhibitory effect of abacavir on HLA-B*57:01 restricted TCR recognition.

EROS is a selective chaperone regulating the phagocyte NADPH oxidase and purinergic signalling

EROS (essential for reactive oxygen species) protein is indispensable for expression of gp91phox, the catalytic core of the phagocyte NADPH oxidase. EROS deficiency in humans is a novel cause of the severe immunodeficiency, chronic granulomatous disease, but its mechanism of action was unknown until now. We elucidate the role of EROS, showing it acts at the earliest stages of gp91phox maturation. It binds the immature 58 kDa gp91phox directly, preventing gp91phox degradation and allowing glycosylation via the oligosaccharyltransferase machinery and the incorporation of the heme prosthetic groups essential for catalysis. EROS also regulates the purine receptors P2X7 and P2X1 through direct interactions, and P2X7 is almost absent in EROS-deficient mouse and human primary cells. Accordingly, lack of murine EROS results in markedly abnormal P2X7 signalling, inflammasome activation, and T cell responses. The loss of both ROS and P2X7 signalling leads to resistance to influenza infection in mice. Our work identifies EROS as a highly selective chaperone for key proteins in innate and adaptive immunity and a rheostat for immunity to infection. It has profound implications for our understanding of immune physiology, ROS dysregulation, and possibly gene therapy.

The HTLV-1 viral oncoproteins Tax and HBZ reprogram the cellular mRNA splicing landscape

Viral infections are known to hijack the transcription and translation of the host cell. However, the extent to which viral proteins coordinate these perturbations remains unclear. Here we used a model system, the human T-cell leukemia virus type 1 (HTLV-1), and systematically analyzed the transcriptome and interactome of key effectors oncoviral proteins Tax and HBZ. We showed that Tax and HBZ target distinct but also common transcription factors. Unexpectedly, we also uncovered a large set of interactions with RNA-binding proteins, including the U2 auxiliary factor large subunit (U2AF2), a key cellular regulator of pre-mRNA splicing. We discovered that Tax and HBZ perturb the splicing landscape by altering cassette exons in opposing manners, with Tax inducing exon inclusion while HBZ induces exon exclusion. Among Tax- and HBZ-dependent splicing changes, we identify events that are also altered in Adult T cell leukemia/lymphoma (ATLL) samples from two independent patient cohorts, and in well-known cancer census genes. Our interactome mapping approach, applicable to other viral oncogenes, has identified spliceosome perturbation as a novel mechanism coordinated by Tax and HBZ to reprogram the transcriptome.

Tax and HBZ are two viral regulatory proteins encoded by the human T-cell leukemia virus type 1 (HTLV-1) via sense and antisense transcripts, respectively. Both proteins are known to drive oncogenic processes that culminate in a T-cell neoplasm, known as Adult T cell leukemia/lymphoma (ATLL). We measured the effects of Tax and HBZ on host gene expression pathway by analyzing the interactome with cellular transcriptional and post-transcriptional regulators, and the transcriptome and mRNA splicing of cell lines expressing either Tax or HBZ. We compared our results with data obtained from independent cohorts of Japanese and Afro-Caribbean patients, and identified common splicing changes that might represent clinically useful biomarkers for ATLL. Finally, we provide evidence that the viral protein Tax can reprogram initial steps of the T-cell transcriptome diversification by hijacking the U2AF complex, a key cellular regulator of pre-mRNA splicing.

ORF Capture-Seq as a versatile method for targeted identification of full-length isoforms

Most human protein-coding genes are expressed as multiple isoforms, which greatly expands the functional repertoire of the encoded proteome. While at least one reliable open reading frame (ORF) model has been assigned for every coding gene, the majority of alternative isoforms remains uncharacterized due to (i) vast differences of overall levels between different isoforms expressed from common genes, and (ii) the difficulty of obtaining full-length transcript sequences. Here, we present ORF Capture-Seq (OCS), a flexible method that addresses both challenges for targeted full-length isoform sequencing applications using collections of cloned ORFs as probes. As a proof-of-concept, we show that an OCS pipeline focused on genes coding for transcription factors increases isoform detection by an order of magnitude when compared to unenriched samples. In short, OCS enables rapid discovery of isoforms from custom-selected genes and will accelerate mapping of the human transcriptome.

A reference map of the human binary protein interactome

Global insights into cellular organization and genome function require comprehensive understanding of the interactome networks that mediate genotype-phenotype relationships^1,2. Here, we present a human “all-by-all” reference interactome map of human binary protein interactions, or “HuRI”. With ~53,000 high-quality protein-protein interactions (PPIs), HuRI has approximately four times more such interactions than high-quality curated interactions from small-scale studies. Integrating HuRI with genome³, transcriptome⁴, and proteome⁵ data enables the study of cellular function within most physiological or pathological cellular contexts. We demonstrate the utility of HuRI in identifying specific subcellular roles of PPIs. Inferred tissue-specific networks reveal general principles for the formation of cellular context-specific functions and elucidate potential molecular mechanisms underlying tissue-specific phenotypes of Mendelian diseases. HuRI represents a systematic proteome-wide reference linking genomic variation to phenotypic outcomes.

Maximizing binary interactome mapping with a minimal number of assays

Complementary assays are required to comprehensively map complex biological entities such as genomes, proteomes and interactome networks. However, how various assays can be optimally combined to approach completeness while maintaining high precision often remains unclear. Here, we propose a framework for binary protein-protein interaction (PPI) mapping based on optimally combining assays and/or assay versions to maximize detection of true positive interactions, while avoiding detection of random protein pairs. We have engineered a novel NanoLuc two-hybrid (N2H) system that integrates 12 different versions, differing by protein expression systems and tagging configurations. The resulting union of N2H versions recovers as many PPIs as 10 distinct assays combined. Thus, to further improve PPI mapping, developing alternative versions of existing assays might be as productive as designing completely new assays. Our findings should be applicable to systematic mapping of other biological landscapes.

APID database: redefining protein-protein interaction experimental evidences and binary interactomes

The collection and integration of all the known protein–protein physical interactions within a proteome framework are critical to allow proper exploration of the protein interaction networks that drive biological processes in cells at molecular level. APID Interactomes is a public resource of biological data (http://apid.dep.usal.es) that provides a comprehensive and curated collection of `protein interactomes’ for more than 1100 organisms, including 30 species with more than 500 interactions, derived from the integration of experimentally detected protein-to-protein physical interactions (PPIs). We have performed an update of APID database including a redefinition of several key properties of the PPIs to provide a more precise data integration and to avoid false duplicated records. This includes the unification of all the PPIs from five primary databases of molecular interactions (BioGRID, DIP, HPRD, IntAct and MINT), plus the information from two original systematic sources of human data and from experimentally resolved 3D structures (i.e. PDBs, Protein Data Bank files, where more than two distinct proteins have been identified). Thus, APID provides PPIs reported in published research articles (with traceable PMIDs) and detected by valid experimental interaction methods that give evidences about such protein interactions (following the `ontology and controlled vocabulary’: www.ebi.ac.uk/ols/ontologies/mi; developed by `HUPO PSI-MI’). Within this data mining framework, all interaction detection methods have been grouped into two main types: (i) `binary’ physical direct detection methods and (ii) `indirect’ methods. As a result of these redefinitions, APID provides unified protein interactomes including the specific `experimental evidences’ that support each PPI, indicating whether the interactions can be considered `binary’ (i.e. supported by at least one binary detection method) or not.

Pooled-matrix protein interaction screens using Barcode Fusion Genetics

High‐throughput binary protein interaction mapping is continuing to extend our understanding of cellular function and disease mechanisms. However, we remain one or two orders of magnitude away from a complete interaction map for humans and other major model organisms. Completion will require screening at substantially larger scales with many complementary assays, requiring further efficiency gains in proteome‐scale interaction mapping. Here, we report Barcode Fusion Genetics‐Yeast Two‐Hybrid (BFG‐Y2H), by which a full matrix of protein pairs can be screened in a single multiplexed strain pool. BFG‐Y2H uses Cre recombination to fuse DNA barcodes from distinct plasmids, generating chimeric protein‐pair barcodes that can be quantified via next‐generation sequencing. We applied BFG‐Y2H to four different matrices ranging in scale from ~25 K to 2.5 M protein pairs. The results show that BFG‐Y2H increases the efficiency of protein matrix screening, with quality that is on par with state‐of‐the‐art Y2H methods.

An inter-species protein-protein interaction network across vast evolutionary distance

In cellular systems, biophysical interactions between macromolecules underlie a complex web of functional interactions. How biophysical and functional networks are coordinated, whether all biophysical interactions correspond to functional interactions, and how such biophysical‐versus‐functional network coordination is shaped by evolutionary forces are all largely unanswered questions. Here, we investigate these questions using an “inter‐interactome” approach. We systematically probed the yeast and human proteomes for interactions between proteins from these two species and functionally characterized the resulting inter‐interactome network. After a billion years of evolutionary divergence, the yeast and human proteomes are still capable of forming a biophysical network with properties that resemble those of intra‐species networks. Although substantially reduced relative to intra‐species networks, the levels of functional overlap in the yeast–human inter‐interactome network uncover significant remnants of co‐functionality widely preserved in the two proteomes beyond human–yeast homologs. Our data support evolutionary selection against biophysical interactions between proteins with little or no co‐functionality. Such non‐functional interactions, however, represent a reservoir from which nascent functional interactions may arise.

Widespread Expansion of Protein Interaction Capabilities by Alternative Splicing

While alternative splicing is known to diversify the functional characteristics of some genes, the extent to which protein isoforms globally contribute to functional complexity on a proteomic scale remains unknown. To address this systematically, we cloned full-length open reading frames of alternatively spliced transcripts for a large number of human genes and used protein-protein interaction profiling to functionally compare hundreds of protein isoform pairs. The majority of isoform pairs share less than 50% of their interactions. In the global context of interactome network maps, alternative isoforms tend to behave like distinct proteins rather than minor variants of each other. Interaction partners specific to alternative isoforms tend to be expressed in a highly tissue-specific manner and belong to distinct functional modules. Our strategy, applicable to other functional characteristics, reveals a widespread expansion of protein interaction capabilities through alternative splicing and suggests that many alternative "isoforms" are functionally divergent (i.e., "functional alloforms").

Widespread macromolecular interaction perturbations in human genetic disorders

How disease-associated mutations impair protein activities in the context of biological networks remains mostly undetermined. Although a few renowned alleles are well characterized, functional information is missing for over 100,000 disease-associated variants. Here we functionally profile several thousand missense mutations across a spectrum of Mendelian disorders using various interaction assays. The majority of disease-associated alleles exhibit wild-type chaperone binding profiles, suggesting they preserve protein folding or stability. While common variants from healthy individuals rarely affect interactions, two-thirds of disease-associated alleles perturb protein-protein interactions, with half corresponding to "edgetic" alleles affecting only a subset of interactions while leaving most other interactions unperturbed. With transcription factors, many alleles that leave protein-protein interactions intact affect DNA binding. Different mutations in the same gene leading to different interaction profiles often result in distinct disease phenotypes. Thus disease-associated alleles that perturb distinct protein activities rather than grossly affecting folding and stability are relatively widespread.

The EBNA3 family of Epstein-Barr virus nuclear proteins associates with the USP46/USP12 deubiquitination complexes to regulate lymphoblastoid cell line growth

The Epstein-Barr virus (EBV) nuclear proteins EBNA3A, EBNA3B, and EBNA3C interact with the cell DNA binding protein RBPJ and regulate cell and viral genes. Repression of the CDKN2A tumor suppressor gene products p16^INK4A and p14^ARF by EBNA3A and EBNA3C is critical for EBV mediated transformation of resting B lymphocytes into immortalized lymphoblastoid cell lines (LCLs). To define the composition of endogenous EBNA3 protein complexes, we generated lymphoblastoid cell lines (LCLs) expressing flag-HA tagged EBNA3A, EBNA3B, or EBNA3C and used tandem affinity purification to isolate each EBNA3 complex. Our results demonstrated that each EBNA3 protein forms a distinct complex with RBPJ. Mass-spectrometry revealed that the EBNA3A and EBNA3B complexes also contained the deubquitylation complex consisting of WDR48, WDR20, and USP46 (or its paralog USP12) and that EBNA3C complexes contained WDR48. Immunoprecipitation confirmed that EBNA3A, EBNA3B, and EBNA3C association with the USP46 complex. Using chromatin immunoprecipitation, we demonstrate that WDR48 and USP46 are recruited to the p14^ARF promoter in an EBNA3C dependent manner. Mapping studies were consistent with WDR48 being the primary mediator of EBNA3 association with the DUB complex. By ChIP assay, WDR48 was recruited to the p14^ARF promoter in an EBNA3C dependent manner. Importantly, WDR48 associated with EBNA3A and EBNA3C domains that are critical for LCL growth, suggesting a role for USP46/USP12 in EBV induced growth transformation.

Epstein-Barr virus (EBV) is a gammaherpesvirus implicated in the pathogenesis of multiple malignancies, including Burkitt lymphoma, Hodgkin lymphoma, post-transplant lymphoproliferative disease (PTLD), nasopharyngeal carcinoma, and gastric carcinoma. EBV infection of resting B-lymphocytes drives them to proliferate as lymphoblastoid cell lines (LCLs), an in vitro model of PTLD. LCLs express a limited EBV gene repertoire, including six nuclear proteins (EBNA1, 2, 3A, 3B, 3C, and LP), three integral membrane proteins (LMP1, 2A, and 2B), and more than 30 micro RNAs. EBNA2 and the EBNA3 proteins are transcription factors that regulate viral and cell gene expression through the cell DNA binding protein RBPJ. In this study, we established LCLs transformed by recombinant EBV genomes in which a Flag-HA epitope tag is fused in-frame to the C-terminus of EBNA3A, EBNA3B or EBNA3C. Using these LCLs, we purified endogenous EBNA3 complexes and identified the USP46 deubiquitinating enzyme (DUB) and its associated chaperones WDR48 and WDR20 as EBNA3 binding proteins. We find that EBNA3s interact primarily with the WDR48 protein and that loss of WDR48 interaction with EBNA3A or EBNA3C impairs LCL growth. This study represents the first characterization of EBNA3 complexes from LCLs and implicates the USP46 DUB complex in EBNA3 mediated gene regulation.

A proteome-scale map of the human interactome network

Just as reference genome sequences revolutionized human genetics, reference maps of interactome networks will be critical to fully understand genotype-phenotype relationships. Here, we describe a systematic map of ?14,000 high-quality human binary protein-protein interactions. At equal quality, this map is ?30% larger than what is available from small-scale studies published in the literature in the last few decades. While currently available information is highly biased and only covers a relatively small portion of the proteome, our systematic map appears strikingly more homogeneous, revealing a "broader" human interactome network than currently appreciated. The map also uncovers significant interconnectivity between known and candidate cancer gene products, providing unbiased evidence for an expanded functional cancer landscape, while demonstrating how high-quality interactome models will help "connect the dots" of the genomic revolution.

Protein interaction network of alternatively spliced isoforms from brain links genetic risk factors for autism

Increased risk for autism spectrum disorders (ASD) is attributed to hundreds of genetic loci. The convergence of ASD variants have been investigated using various approaches, including protein interactions extracted from the published literature. However, these datasets are frequently incomplete, carry biases and are limited to interactions of a single splicing isoform, which may not be expressed in the disease-relevant tissue. Here we introduce a new interactome mapping approach by experimentally identifying interactions between brain-expressed alternatively spliced variants of ASD risk factors. The Autism Spliceform Interaction Network reveals that almost half of the detected interactions and about 30% of the newly identified interacting partners represent contribution from splicing variants, emphasizing the importance of isoform networks. Isoform interactions greatly contribute to establishing direct physical connections between proteins from the de novo autism CNVs. Our findings demonstrate the critical role of spliceform networks for translating genetic knowledge into a better understanding of human diseases.

An RS motif within the Epstein-Barr virus BLRF2 tegument protein is phosphorylated by SRPK2 and is important for viral replication

Epstein-Barr virus (EBV) is a gammaherpesvirus that causes infectious mononucleosis, B cell lymphomas, and nasopharyngeal carcinoma. Many of the genes required for EBV virion morphogenesis are found in all herpesviruses, but some are specific to gammaherpesviruses. One of these gamma-specific genes, BLRF2, encodes a tegument protein that has been shown to be essential for replication in other gammaherpesviruses. In this study, we identify BLRF2 interacting proteins using binary and co-complex protein assays. Serine/Arginine-rich Protein Kinase 2 (SRPK2) was identified by both assays and was further shown to phosphorylate an RS motif in the BLRF2 C-terminus. Mutation of this RS motif (S148A+S150A) abrogated the ability of BLRF2 to support replication of a murine gammaherpesvirus 68 genome lacking the BLRF2 homolog (ORF52). We conclude that the BLRF2 RS motif is phosphorylated by SRPK2 and is important for viral replication.

Identification of FAM111A as an SV40 host range restriction and adenovirus helper factor

The small genome of polyomaviruses encodes a limited number of proteins that are highly dependent on interactions with host cell proteins for efficient viral replication. The SV40 large T antigen (LT) contains several discrete functional domains including the LXCXE or RB-binding motif, the DNA binding and helicase domains that contribute to the viral life cycle. In addition, the LT C-terminal region contains the host range and adenovirus helper functions required for lytic infection in certain restrictive cell types. To understand how LT affects the host cell to facilitate viral replication, we expressed full-length or functional domains of LT in cells, identified interacting host proteins and carried out expression profiling. LT perturbed the expression of p53 target genes and subsets of cell-cycle dependent genes regulated by the DREAM and the B-Myb-MuvB complexes. Affinity purification of LT followed by mass spectrometry revealed a specific interaction between the LT C-terminal region and FAM111A, a previously uncharacterized protein. Depletion of FAM111A recapitulated the effects of heterologous expression of the LT C-terminal region, including increased viral gene expression and lytic infection of SV40 host range mutants and adenovirus replication in restrictive cells. FAM111A functions as a host range restriction factor that is specifically targeted by SV40 LT.

Viruses have evolved numerous mechanisms to counteract host cell defenses to facilitate productive infection. Simian Virus 40 (SV40) replication depends on specific interactions between large T antigen (LT) and a wide variety of host cell proteins. Although the LT C-terminal region has no evident enzymatic activity, mutations or deletions of this region significantly reduce the ability of the virus to replicate in restrictive cell types. Here, we identified host proteins that bind to LT and determined that the LT C-terminal region binds specifically to FAM111A. This physical interaction was required for efficient viral replication and sustained viral gene expression in restrictive cell types. In addition, RNAi-mediated knockdown of FAM111A levels in restrictive cells restored lytic infection of SV40 host range mutants and human adenovirus. These results indicate that FAM111A plays an important role in viral host range restriction. Our study provides insights into the viral-host perturbations caused by SV40 LT and the interaction of viruses with host restriction factors.

Interpreting cancer genomes using systematic host network perturbations by tumour virus proteins

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, and large numbers of somatic genomic alterations, associated with a predisposition to cancer. However, it remains difficult to distinguish background, or 'passenger', cancer mutations from causal, or 'driver', mutations in these data sets. Human viruses intrinsically depend on their host cell during the course of infection and can elicit pathological phenotypes similar to those arising from mutations. Here we test the hypothesis that genomic variations and tumour viruses may cause cancer through related mechanisms, by systematically examining 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, such as Notch signalling and apoptosis, that go awry in cancer. We show that systematic analyses of host targets of viral proteins can identify cancer genes with a success rate on a par with their identification through functional genomics and large-scale cataloguing of tumour mutations. Together, these complementary approaches increase the specificity of cancer gene identification. Combining systems-level studies of pathogen-encoded gene products with genomic approaches will facilitate the prioritization of cancer-causing driver genes to advance the understanding of the genetic basis of human cancer.

Viral perturbations of host networks reflect disease etiology

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.

Epstein-Barr virus nuclear protein 3C binds to the N-terminal (NTD) and beta trefoil domains (BTD) of RBP/CSL; only the NTD interaction is essential for lymphoblastoid cell growth

Association of EBV nuclear proteins EBNA2, EBNA3A and EBNA3C with RBP/CSL, is essential for lymphoblastoid cell line (LCL) proliferation. Conserved residues in the EBNA3 homology domain, required for RBP/CSL interaction, lack the WΦP motif that mediates EBNA2 and Notch binding to the RBP/CSL beta-trefoil domain (BTD). We map RBP/CSL interacting residues within EBNA3A(aa128-204) and EBNA3C(aa211-233). The EBNA3A results are consistent with an earlier report (aa125-222), but the EBNA3C domain is unexpectedly small and includes a "WTP" sequence. This EBNA3C WTP motif confers RBP/CSL binding in vitro, in yeast, and in mammalian cells. Further, an EBNA3C WTP→STP(W227S) mutation impaired BTD binding whereas EBNA3 homology domain mutations disrupted RBP/CSL N-terminal domain (NTD) binding. WTP was not essential for EBNA3C repression of EBNA2 in reporter assays or for maintenance of LCL growth. Our results indicate that EBNA3 proteins interact with multiple RBP/CSL domains, but only NTD interactions are required for LCL growth.

Epstein-Barr virus and virus human protein interaction maps

A comprehensive mapping of interactions among Epstein-Barr virus (EBV) proteins and interactions of EBV proteins with human proteins should provide specific hypotheses and a broad perspective on EBV strategies for replication and persistence. Interactions of EBV proteins with each other and with human proteins were assessed by using a stringent high-throughput yeast two-hybrid system. Overall, 43 interactions between EBV proteins and 173 interactions between EBV and human proteins were identified. EBV-EBV and EBV-human protein interaction, or "interactome" maps provided a framework for hypotheses of protein function. For example, LF2, an EBV protein of unknown function interacted with the EBV immediate early R transactivator (Rta) and was found to inhibit Rta transactivation. From a broader perspective, EBV genes can be divided into two evolutionary classes, "core" genes, which are conserved across all herpesviruses and subfamily specific, or "noncore" genes. Our EBV-EBV interactome map is enriched for interactions among proteins in the same evolutionary class. Furthermore, human proteins targeted by EBV proteins were enriched for highly connected or "hub" proteins and for proteins with relatively short paths to all other proteins in the human interactome network. Targeting of hubs might be an efficient mechanism for EBV reorganization of cellular processes.

Development of herpesvirus-based episomally maintained gene delivery vectors

Successful gene therapy aims to deliver and express therapeutic genes to cure or slow the progression of disease. However, a major obstacle in the application of gene therapy has been the development of the vectors used to deliver heterologous DNA to the cell or tissue of choice. A number of viral- and non-viral-based vector systems have undergone clinical trials with varying success. However, at present, no vector system possesses the full complement of properties that are generally believed necessary in an ideal gene delivery system. Therefore, alongside attempts to improve current gene delivery vectors, the identification and evaluation of new viral vectors is crucial for the long-term success of gene therapy. Herpesviruses are large DNA viruses which possess a number of advantages as gene delivery vectors. These relate to an ability to package large DNA insertions and establish lifelong latent infections in which the genomic material exists as a stable episome. This review aims to highlight the potential of herpesvirus vectors, in particular an alternative vector system based on herpesvirus saimiri (HVS). HVS is capable of infecting a range of human cell lines with high efficiencies, and the viral genome persists as high copy number, circular, non-integrated episomes which segregate to progeny following cell division. This allows the virus-based vector to stably transduce a dividing cell population and provide sustained transgene expression for an extended period of time both in vitro and in vivo. Moreover, the insertion of a bacterial artificial chromosome cassette into the HVS genome simplifies the incorporation of large amounts of heterologous DNA for gene delivery. These properties offer characteristics similar to an artificial chromosome combined with an efficient delivery system and merit its continual development as a possible gene delivery vector for the future.

The herpesvirus saimiri ORF73 gene product interacts with host-cell mitotic chromosomes and self-associates via its C terminus

The herpesvirus saimiri (HVS) ORF73 gene product shares limited homology with the ORF73 protein of Kaposi's sarcoma-associated herpesvirus (KSHV). ORF73 is expressed in an in vitro model of HVS latency, where the genome persists as a non-integrated circular episome. This suggests it may have a similar role to KSHV ORF73 in episomal maintenance, by tethering viral genomes to host-cell chromosomes. Here, the association of ORF73 with host mitotic chromosomes is described. Deletion analysis demonstrates that the distal 123 aa of the ORF73 protein are required for mitotic chromosomal localization and for self-association. Moreover, deletion of the extreme C terminus disrupts both self-association and host mitotic chromosome colocalization. This suggests that HVS ORF73 has a similar role to KSHV ORF73 in episomal maintenance and that association of ORF73 to host mitotic chromosomes is dependent on its ability to form multimers.

Generation and precise modification of a herpesvirus saimiri bacterial artificial chromosome demonstrates that the terminal repeats are required for both virus production and episomal persistence

Herpesvirus saimiri (HVS) is the prototype gamma-2 herpesvirus, and shares considerable homology with the human gammaherpesviruses Kaposi's sarcoma-associated herpesvirus and Epstein–Barr virus. The generation of herpesvirus mutants is a key facet in the study of virus biology. The use of F-factor-based bacterial artificial chromosomes (BACs) to clone and modify the genomes of herpesviruses has enhanced the variety, precision and simplicity of mutant production. Here we describe the cloning of the genome of HVS non-transforming strain A11-S4 into a BAC. The cloning of the BAC elements disrupts open reading frame (ORF) 15 but the HVS-BAC can still replicate at levels similar to wild-type virus, and can persistently infect fibroblasts. The HVS-BAC was modified by RecA-mediated recombination initially to substitute reporter genes and also to delete the terminal repeats (TR). After deletion of the TR, the HVS-BAC fails to enter a productive virus lytic cycle, and cannot establish a persistent episomal infection when transfected into fibroblast cell lines. This shows that while ORF 15 is dispensable for virus function in vitro, the TR is required for both virus latency and lytic virus production. In addition, the HVS-BAC promises to be a valuable tool that can be used for the routine and precise production and analysis of viral mutants to further explore gammaherpesvirus biology.

The Herpesvirus Saimiri open reading frame 73 gene product interacts with the cellular protein p32

The role of the gamma-2 herpesvirus open reading frame (ORF) 73 gene product has become the focus of considerable interest. It has recently been shown that the Kaposi's sarcoma-associated herpesvirus (KSHV) latency-associated nuclear antigen (LANA) is expressed during a latent infection and can modulate both viral and cellular gene expression. The herpesvirus saimiri (HVS) ORF 73 gene product has some sequence homology to LANA; however, the role of HVS ORF 73 is unknown. We have previously demonstrated that HVS ORF73 is expressed in a stably transduced human carcinoma cell line, where HVS genomes persist as nonintegrated circular episomes. This implies that there may be some functional homology between these proteins. To further investigate the role of the HVS ORF 73 protein, the yeast two-hybrid system was employed to identify interacting cellular proteins. We demonstrate that ORF 73 interacts with the cellular protein p32 and triggers the accumulation of p32 in the nucleus. Using reporter gene-based transient-transfection assays, we demonstrate that ORF 73 can transactivate a number of heterologous promoter constructs and also upregulate its own promoter. Moreover, ORF 73 and p32 act synergistically to transactivate these promoters. The binding of ORF 73 to p32 is mediated by an amino-terminal arginine-rich domain, which contains two functionally distinct nuclear localization signals. The p32 binding domains are required for ORF 73 transactivating abilities and for ORF 73 to induce nuclear accumulation of p32. These results suggest that ORF 73 can function as a regulator of gene expression and that p32 is involved in ORF 73-dependent transcriptional activation.

Characterization of the herpesvirus saimiri ORF73 gene product

The herpesvirus saimiri (HVS) gene product encoded by ORF73 shares a limited homology with the ORF73 encoded protein of Kaposi's sarcoma-associated herpesvirus (KSHV). It has recently been shown that the KSHV ORF73 protein is expressed during a latent infection and co-localizes with host cell chromosomes, suggesting that it plays a role in episomal maintenance by tethering viral genomes to host cell chromosomes. At present the role of the HVS ORF73 gene product is unknown. However, the expression of HVS ORF73 in a stably transduced human carcinoma cell line, where the HVS genome persists as a non-integrated circular episome, has recently been shown. In this report, the characterization of the HVS ORF73 protein and the mapping of its functional domains are described. The results suggest that the HVS ORF73 gene encodes a 64 kDa nuclear protein. Moreover, the amino terminus contains two functional nuclear localization signals, whereas the carboxy terminus is required for the distinctive speckled nuclear distribution pattern as observed with both the HVS and KSHV ORF73 proteins.

The carboxy terminus of the herpesvirus saimiri ORF 57 gene contains domains that are required for transactivation and transrepression

Herpesvirus saimiri (HVS) ORF 57 is homologous to genes identified in all classes of herpesviruses. We have previously shown that ORF 57 encodes a multifunctional protein, responsible for both transactivation and repression of viral gene expression at a post-transcriptional level. This suggests that the ORF 57 protein shares some functional similarities with the herpes simplex virus IE63/ICP27 and Epstein-Barr virus Mta proteins. However, little is known about the functional domains responsible for the properties of ORF 57 due to the limited homology shared between these proteins. In this report, we have identified the functional domains responsible for transactivation and repression by the ORF 57 protein. We demonstrate that the carboxy terminus is required for ORF 57 transactivation, repression and an intense SC-35 nuclear spotting. This region contains two highly conserved motifs amongst its homologues, a zinc finger-like motif and a highly hydrophobic domain. We further show that the hydrophobic domain is required for transactivation and is also involved in nuclear localization of the ORF 57 protein, whereas the zinc finger-like domain is required for transactivation, repression and the intense SC-35 nuclear spotting.

Analysis of gene expression in a human cell line stably transduced with herpesvirus saimiri

Herpesvirus saimiri (HVS) is the prototype gamma-2 herpesvirus; it has significant homology to the human gammaherpesviruses Kaposi's sarcoma-associated virus and Epstein-Barr virus and the murine gammaherpesvirus murine herpesvirus 68. HVS causes a persistent asymptomatic infection in its natural host, the squirrel monkey. Both subgroups A and C possess the ability to immortalize common marmoset T lymphocytes to interleukin-2-independent proliferation. However, only subgroup C is capable of transforming human, rabbit, and rhesus monkey lymphocytes in vitro. In addition, HVS can stably transduce a variety of human cell lines where the virus persists as a nonintegrating circular episome. In this study, we have developed a system in which the HVS DNA is stably maintained as a nonintegrated circular episome in the human lung carcinoma cell line A549. Virus production can be reactivated using chemical inducing agents, including tetradecanoyl phorbol acetate and n-butyrate, suggesting that the infection in human A549 cells is latent. To analyze virus gene expression in these stably transduced cells, Northern blot analysis was performed using a series of probes produced from restriction fragments spanning the entire coding region of the HVS genome. This demonstrated that an adjacent set of genes containing open reading frames (ORFs) 71 to 73 are expressed in this stably transduced cell line. Moreover, these genes are transcribed as a polycistronic mRNA species produced from a common promoter upstream of ORF 73. This model may serve as a useful tool in the further analysis of the role of ORFs 71 to 73 in gamma-2 herpesvirus latency.

Specific oncolytic activity of herpesvirus saimiri in pancreatic cancer cells

The potential use of oncolytic viruses in the treatment of cancer has been investigated for some time. A variety of agents have been studied, including some which appear to be selectively replication-competent in cancer cell lines. In this study, we have investigated the ability of herpesvirus saimiri to specifically lyse selected human cancer cell lines. Upon infection with a replication-competent virus carrying the EGFP reporter gene and a neomycin resistance marker, the pancreatic cancer lines MIAPACA and PANC-1 exhibited definite cytopathic effects. In contrast, the colonic carcinoma cell lines SW480 and HCT116 were phenotypically unaltered. In addition, stable cell lines could not be generated from PANC-1 infected cultures, in marked contrast to cultures of cells from other human tissues. Virus recovery assays demonstrated that all of the cell lines produced a small amount of virus post-infection, but that virus replication was minimal after 1 week in culture. In addition, treatment with acyclovir inhibited virus replication but paradoxically increased cytopathic effect. These data suggest that herpesvirus saimiri may have potential as an oncolytic agent for the treatment of pancreatic cancer.