Identifying Pathway Proteins in Networks using Convergence

One of the key goals of systems biology concerns the analysis of experimental biological data available to the scientific public. New technologies are rapidly developed to observe and report whole-scale biological phenomena; however, few methods exist with the ability to produce specific, testable hypotheses from this noisy ‘big’ data. In this work, we propose an approach that combines the power of data-driven network theory along with knowledge-based ontology to tackle this problem. Network models are especially powerful due to their ability to display elements of interest and their relationships as internetwork structures. Additionally, ontological data actually supplements the confidence of relationships within the model without clouding critical structure identification. As such, we postulate that given a (gene/protein) marker set of interest, we can systematically identify the core of their interactions (if they are indeed working together toward a biological function), via elimination of original markers and addition of additional necessary markers. This concept, which we refer to as “convergence,” harnesses the idea of “guilt-by-association” and recursion to identify whether a core of relationships exists between markers. In this study, we test graph theoretic concepts such as shortest-path, k-Nearest- Neighbor and clustering) to identify cores iteratively in data- and knowledge-based networks in the canonical yeast Pheromone Mating Response pathway. Additionally, we provide results for convergence application in virus infection, hearing loss, and Parkinson’s disease. Our results indicate that if a marker set has common discrete function, this approach is able to identify that function, its interacting markers, and any new elements necessary to complete the structural core of that function. The result below find that the shortest path function is the best approach of those used, finding small target sets that contain a majority or all of the markers in the gold standard pathway. The power of this approach lies in its ability to be used in investigative studies to inform decisions concerning target selection

Structure-based prediction of insertion-site preferences of transposons into chromosomes

Mobile genetic elements with the ability to integrate genetic information into chromosomes can cause disease over short periods of time and shape genomes over eons. These elements can be used for functional genomics, gene transfer and human gene therapy. However, their integration-site preferences, which are critically important for these uses, are poorly understood. We analyzed the insertion sites of several transposons and retroviruses to detect patterns of integration that might be useful for prediction of preferred integration sites. Initially we found that a mathematical description of DNA-deformability, called V(step), could be used to distinguish preferential integration sites for Sleeping Beauty (SB) transposons into a particular 100 bp region of a plasmid [G. Liu, A. M. Geurts, K. Yae, A. R. Srinivassan, S. C. Fahrenkrug, D. A. Largaespada,J. Takeda, K. Horie, W. K. Olson and P. B. Hackett (2005) J. Mol. Biol., 346, 161–173 ]. Based on these findings, we extended our examination of integration of SB transposons into whole plasmids and chromosomal DNA. To accommodate sequences up to 3 Mb for these analyses, we developed an automated method, ProTIS(©), that can generate profiles of predicted integration events. However, a similar approach did not reveal any structural pattern of DNA that could be used to predict favored integration sites for other transposons as well as retroviruses and lentiviruses due to a limitation of available data sets. Nonetheless, ProTIS(©) has the utility for predicting likely SB transposon integration sites in investigator-selected regions of genomes and our general strategy may be useful for other mobile elements once a sufficiently high density of sites in a single region are obtained. ProTIS analysis can be useful for functional genomic, gene transfer and human gene therapy applications using the SB system

Functional Characteristics of a Highly Specific Integrase Encoded by an LTR-Retrotransposon

Background: The retroviral Integrase protein catalyzes the insertion of linear viral DNA into host cell DNA. Although different retroviruses have been shown to target distinctive chromosomal regions, few of them display a site-specific integration. ZAM, a retroelement from Drosophila melanogaster very similar in structure and replication cycle to mammalian retroviruses is highly site-specific. Indeed, ZAM copies target the genomic 59-CGCGCg-39 consensus-sequences. To enlighten the determinants of this high integration specificity, we investigated the functional properties of its integrase protein denoted ZAM-IN. Principal Findings: Here we show that ZAM-IN displays the property to nick DNA molecules in vitro. This endonuclease activity targets specific sequences that are present in a 388 bp fragment taken from the white locus and known to be a genomic ZAM integration site in vivo. Furthermore, ZAM-IN displays the unusual property to directly bind specific genomic DNA sequences. Two specific and independent sites are recognized within the 388 bp fragment of the white locus: the CGCGCg sequence and a closely apposed site different in sequence. Conclusion: This study strongly argues that the intrinsic properties of ZAM-IN, ie its binding properties and its endonuclease activity, play an important part in ZAM integration specificity. Its ability to select two binding sites and to nick the DNA molecule reminds the strategy used by some site-specific recombination enzymes and forms the basis for site-specifi

Integrating transposable elements in the 3D genome

Chromosome organisation is increasingly recognised as an essential component of genome regulation, cell fate and cell health. Within the realm of transposable elements (TEs) however, the spatial information of how genomes are folded is still only rarely integrated in experimental studies or accounted for in modelling. Whilst polymer physics is recognised as an important tool to understand the mechanisms of genome folding, in this commentary we discuss its potential applicability to aspects of TE biology. Based on recent works on the relationship between genome organisation and TE integration, we argue that existing polymer models may be extended to create a predictive framework for the study of TE integration patterns. We suggest that these models may offer orthogonal and generic insights into the integration profiles (or "topography") of TEs across organisms. In addition, we provide simple polymer physics arguments and preliminary molecular dynamics simulations of TEs inserting into heterogeneously flexible polymers. By considering this simple model, we show how polymer folding and local flexibility may generically affect TE integration patterns. The preliminary discussion reported in this commentary is aimed to lay the foundations for a large-scale analysis of TE integration dynamics and topography as a function of the three-dimensional host genome

Retroviral insertions in the VISION database identify molecular pathways in mouse lymphoid leukemia and lymphoma

AKXD recombinant inbred (RI) strains develop a variety of leukemias and lymphomas due to somatically acquired insertions of retroviral DNA into the genome of hematopoetic cells that can mutate cellular proto-oncogenes and tumor suppressor genes. We generated a new set of tumors from nine AKXD RI strains selected for their propensity to develop B-cell tumors, the most common type of human hematopoietic cancers. We employed a PCR technique called viral insertion site amplification (VISA) to rapidly isolate genomic sequence at the site of provirus insertion. Here we describe 550 VISA sequence tags (VSTs) that identify 74 common insertion sites (CISs), of which 21 have not been identified previously. Several suspected proto-oncogenes and tumor suppressor genes lie near CISs, providing supportive evidence for their roles in cancer. Furthermore, numerous previously uncharacterized genes lie near CISs, providing a pool of candidate disease genes for future research. Pathway analysis of candidate genes identified several signaling pathways as common and powerful routes to blood cancer, including Notch, E-protein, NFκB, and Ras signaling. Misregulation of several Notch signaling genes was confirmed by quantitative RT-PCR. Our data suggest that analyses of insertional mutagenesis on a single genetic background are biased toward the identification of cooperating mutations. This tumor collection represents the most comprehensive study of the genetics of B-cell leukemia and lymphoma development in mice. We have deposited the VST sequences, CISs in a genome viewer, histopathology, and molecular tumor typing data in a public web database called VISION (Viral Insertion Sites Identifying Oncogenes), which is located at http://www.mouse-genome.bcm.tmc.edu/vision

Novel Escape Mutants Suggest an Extensive TRIM5α Binding Site Spanning the Entire Outer Surface of the Murine Leukemia Virus Capsid Protein

After entry into target cells, retroviruses encounter the host restriction factors such as Fv1 and TRIM5α. While it is clear that these factors target retrovirus capsid proteins (CA), recognition remains poorly defined in the absence of structural information. To better understand the binding interaction between TRIM5α and CA, we selected a panel of novel N-tropic murine leukaemia virus (N-MLV) escape mutants by a serial passage of replication competent N-MLV in rhesus macaque TRIM5α (rhTRIM5α)-positive cells using a small percentage of unrestricted cells to allow multiple rounds of virus replication. The newly identified mutations, many of which involve changes in charge, are distributed over the outer ‘top’ surface of N-MLV CA, including the N-terminal β-hairpin, and map up to 29 Ao apart. Biological characterisation with a number of restriction factors revealed that only one of the new mutations affects restriction by human TRIM5α, indicating significant differences in the binding interaction between N-MLV and the two TRIM5αs, whereas three of the mutations result in dual sensitivity to Fv1n and Fv1b. Structural studies of two mutants show that no major changes in the overall CA conformation are associated with escape from restriction. We conclude that interactions involving much, if not all, of the surface of CA are vital for TRIM5α binding

The Ankyrin Repeats and DHHC S-acyl Transferase Domain of AKR1 Act Independently to Regulate Switching from Vegetative to Mating States in Yeast

Signal transduction from G-protein coupled receptors to MAPK cascades through heterotrimeric G-proteins has been described for many eukaryotic systems. One of the best-characterised examples is the yeast pheromone response pathway, which is negatively regulated by AKR1. AKR1-like proteins are present in all eukaryotes and contain a DHHC domain and six ankyrin repeats. Whilst the DHHC domain dependant S-acyl transferase (palmitoyl transferase) function of AKR1 is well documented it is not known whether the ankyrin repeats are also required for this activity. Here we show that the ankyrin repeats of AKR1 are required for full suppression of the yeast pheromone response pathway, by sequestration of the Gβγ dimer, and act independently of AKR1 S-acylation function. Importantly, the functions provided by the AKR1 ankyrin repeats and DHHC domain are not required on the same molecule to fully restore WT phenotypes and function. We also show that AKR1 molecules are S-acylated at locations other than the DHHC cysteine, increasing the abundance of AKR1 in the cell. Our results have important consequences for studies of AKR1 function, including recent attempts to characterise S-acylation enzymology and kinetics. Proteins similar to AKR1 are found in all eukaryotes and our results have broad implications for future work on these proteins and the control of switching between Gβγ regulated pathways

Deciphering the Code for Retroviral Integration Target Site Selection

Upon cell invasion, retroviruses generate a DNA copy of their RNA genome and integrate retroviral cDNA within host chromosomal DNA. Integration occurs throughout the host cell genome, but target site selection is not random. Each subgroup of retrovirus is distinguished from the others by attraction to particular features on chromosomes. Despite extensive efforts to identify host factors that interact with retrovirion components or chromosome features predictive of integration, little is known about how integration sites are selected. We attempted to identify markers predictive of retroviral integration by exploiting Precision-Recall methods for extracting information from highly skewed datasets to derive robust and discriminating measures of association. ChIPSeq datasets for more than 60 factors were compared with 14 retroviral integration datasets. When compared with MLV, PERV or XMRV integration sites, strong association was observed with STAT1, acetylation of H3 and H4 at several positions, and methylation of H2AZ, H3K4, and K9. By combining peaks from ChIPSeq datasets, a supermarker was identified that localized within 2 kB of 75% of MLV proviruses and detected differences in integration preferences among different cell types. The supermarker predicted the likelihood of integration within specific chromosomal regions in a cell-type specific manner, yielding probabilities for integration into proto-oncogene LMO2 identical to experimentally determined values. The supermarker thus identifies chromosomal features highly favored for retroviral integration, provides clues to the mechanism by which retrovirus integration sites are selected, and offers a tool for predicting cell-type specific proto-oncogene activation by retroviruses

Proteasomal Degradation of TRIM5α during Retrovirus Restriction

The host protein TRIM5α inhibits retroviral infection at an early post-penetration stage by targeting the incoming viral capsid. While the detailed mechanism of restriction remains unclear, recent studies have implicated the activity of cellular proteasomes in the restriction of retroviral reverse transcription imposed by TRIM5α. Here, we show that TRIM5α is rapidly degraded upon encounter of a restriction-susceptible retroviral core. Inoculation of TRIM5α-expressing human 293T cells with a saturating level of HIV-1 particles resulted in accelerated degradation of the HIV-1-restrictive rhesus macaque TRIM5α protein but not the nonrestrictive human TRIM5α protein. Exposure of cells to HIV-1 also destabilized the owl monkey restriction factor TRIMCyp; this was prevented by addition of the inhibitor cyclosporin A and was not observed with an HIV-1 virus containing a mutation in the capsid protein that relieves restriction by TRIMCyp IVHIV. Likewise, human TRIM5α was rapidly degraded upon encounter of the restriction-sensitive N-tropic murine leukemia virus (N-MLV) but not the unrestricted B-MLV. Pretreatment of cells with proteasome inhibitors prevented the HIV-1-induced loss of both rhesus macaque TRIM5α and TRIMCyp proteins. We also detected degradation of endogenous TRIM5α in rhesus macaque cells following HIV-1 infection. We conclude that engagement of a restriction-sensitive retrovirus core results in TRIM5α degradation by a proteasome-dependent mechanism