CytoGEDEVO-global alignment of biological networks with Cytoscape.

MOTIVATION: In the systems biology era, high-throughput omics technologies have enabled the unraveling of the interplay of some biological entities on a large scale (e.g. genes, proteins, metabolites or RNAs). Huge biological networks have emerged, where nodes correspond to these entities and edges between them model their relations. Protein-protein interaction networks, for instance, show the physical interactions of proteins in an organism. The comparison of such networks promises additional insights into protein and cell function as well as knowledge-transfer across species. Several computational approaches have been developed previously to solve the network alignment (NA) problem, but only a few concentrate on the usability of the implemented tools for the evaluation of protein-protein interactions by the end users (biologists and medical researchers). RESULTS: We have created CytoGEDEVO, a Cytoscape app for visual and user-assisted NA. It extends the previous GEDEVO methodology for global pairwise NAs with new graphical and functional features. Our main focus was on the usability, even by non-programmers and the interpretability of the NA results with Cytoscape. AVAILABILITY AND IMPLEMENTATION: CytoGEDEVO is publicly available from the Cytoscape app store at http://apps.cytoscape.org/apps/cytogedevo In addition, we provide stand-alone command line executables, source code, documentation and step-by-step user instructions at http://cytogedevo.compbio.sdu.dk CONTACT: malek@tugraz.at SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Identification of HNRNPK as regulator of hepatitis C virus particle production.

Hepatitis C virus (HCV) is a major cause of chronic liver disease affecting around 130 million people worldwide. While great progress has been made to define the principle steps of the viral life cycle, detailed knowledge how HCV interacts with its host cells is still limited. To overcome this limitation we conducted a comprehensive whole-virus RNA interference-based screen and identified 40 host dependency and 16 host restriction factors involved in HCV entry/replication or assembly/release. Of these factors, heterogeneous nuclear ribonucleoprotein K (HNRNPK) was found to suppress HCV particle production without affecting viral RNA replication. This suppression of virus production was specific to HCV, independent from assembly competence and genotype, and not found with the related Dengue virus. By using a knock-down rescue approach we identified the domains within HNRNPK required for suppression of HCV particle production. Importantly, HNRNPK was found to interact specifically with HCV RNA and this interaction was impaired by mutations that also reduced the ability to suppress HCV particle production. Finally, we found that in HCV-infected cells, subcellular distribution of HNRNPK was altered; the protein was recruited to sites in close proximity of lipid droplets and colocalized with core protein as well as HCV plus-strand RNA, which was not the case with HNRNPK variants unable to suppress HCV virion formation. These results suggest that HNRNPK might determine efficiency of HCV particle production by limiting the availability of viral RNA for incorporation into virions. This study adds a new function to HNRNPK that acts as central hub in the replication cycle of multiple other viruses.

A genome-wide association scan of nonsynonymous SNPs identifies a susceptibility variant for Crohn disease in ATG16L1.

We performed a genome-wide association study of 19,779 nonsynonymous SNPs in 735 individuals with Crohn disease and 368 controls. A total of 7,159 of these SNPs were informative. We followed up on all 72 SNPs with P 0.4), these data suggest that the underlying biological process may be specific to Crohn disease.

Analyzing and visualizing residue networks of protein structures.

The study of individual amino acid residues and their molecular interactions in protein structures is crucial for understanding structure-function relationships. Recent work has indicated that residue networks derived from 3D protein structures provide additional insights into the structural and functional roles of interacting residues. Here, we present the new software tools RINerator and RINalyzer for the automatized generation, 2D visualization, and interactive analysis of residue interaction networks, and highlight their use in different application scenarios.

Association between variants of PRDM1 and NDP52 and Crohn's disease, based on exome sequencing and functional studies.

BACKGROUND & AIMS: Genome-wide association studies (GWAS) have identified 140 Crohn's disease (CD) susceptibility loci. For most loci, the variants that cause disease are not known and the genes affected by these variants have not been identified. We aimed to identify variants that cause CD through detailed sequencing, genetic association, expression, and functional studies. METHODS: We sequenced whole exomes of 42 unrelated subjects with CD and 5 healthy subjects (controls) and then filtered single nucleotide variants by incorporating association results from meta-analyses of CD GWAS and in silico mutation effect prediction algorithms. We then genotyped 9348 subjects with CD, 2868 subjects with ulcerative colitis, and 14,567 control subjects and associated variants analyzed in functional studies using materials from subjects and controls and in vitro model systems. RESULTS: We identified rare missense mutations in PR domain-containing 1 (PRDM1) and associated these with CD. These mutations increased proliferation of T cells and secretion of cytokines on activation and increased expression of the adhesion molecule L-selectin. A common CD risk allele, identified in GWAS, correlated with reduced expression of PRDM1 in ileal biopsy specimens and peripheral blood mononuclear cells (combined P = 1.6 × 10(-8)). We identified an association between CD and a common missense variant, Val248Ala, in nuclear domain 10 protein 52 (NDP52) (P = 4.83 × 10(-9)). We found that this variant impairs the regulatory functions of NDP52 to inhibit nuclear factor κB activation of genes that regulate inflammation and affect the stability of proteins in Toll-like receptor pathways. CONCLUSIONS: We have extended the results of GWAS and provide evidence that variants in PRDM1 and NDP52 determine susceptibility to CD. PRDM1 maps adjacent to a CD interval identified in GWAS and encodes a transcription factor expressed by T and B cells. NDP52 is an adaptor protein that functions in selective autophagy of intracellular bacteria and signaling molecules, supporting the role of autophagy in the pathogenesis of CD.

NetworkPrioritizer: a versatile tool for network-based prioritization of candidate disease genes or other molecules.

SUMMARY: The prioritization of candidate disease genes is often based on integrated datasets and their network representation with genes as nodes connected by edges for biological relationships. However, the majority of prioritization methods does not allow for a straightforward integration of the user's own input data. Therefore, we developed the Cytoscape plugin NetworkPrioritizer that particularly supports the integrative network-based prioritization of candidate disease genes or other molecules. Our versatile software tool computes a number of important centrality measures to rank nodes based on their relevance for network connectivity and provides different methods to aggregate and compare rankings. AVAILABILITY: NetworkPrioritizer and the online documentation are freely available at http://www.networkprioritizer.de

Sarcoidosis is associated with a truncating splice site mutation in BTNL2.

Sarcoidosis is a polygenic immune disorder with predominant manifestation in the lung. Genome-wide linkage analysis previously indicated that the extended major histocompatibility locus on chromosome 6p was linked to susceptibility to sarcoidosis. Here, we carried out a systematic three-stage SNP scan of 16.4 Mb on chromosome 6p21 in as many as 947 independent cases of familial and sporadic sarcoidosis and found that a 15-kb segment of the gene butyrophilin-like 2 (BTNL2) was associated with the disease. The primary disease-associated variant (rs2076530; P(TDT) = 3 x 10(-6), P(case-control) = 1.1 x 10(-8); replication P(TDT) = 0.0018, P(case-control) = 1.8 x 10(-6)) represents a risk factor that is independent of variation in HLA-DRB1. BTNL2 is a member of the immunoglobulin superfamily and has been implicated as a costimulatory molecule involved in T-cell activation on the basis of its homology to B7-1. The G --> A transition constituting rs2076530 leads to the use of a cryptic splice site located 4 bp upstream of the affected wild-type donor site. Transcripts of the risk-associated allele have a premature stop in the spliced mRNA. The resulting protein lacks the C-terminal IgC domain and transmembrane helix, thereby disrupting the membrane localization of the protein, as shown in experiments using green fluorescent protein and V5 fusion proteins.

Characterization of the EGFR interactome reveals associated protein complex networks and intracellular receptor dynamics.

Growth factor receptor mediated signaling is meanwhile recognized as a complex signaling network, which is initiated by recruiting specific patterns of adaptor proteins to the intracellular domain of epidermal growth factor receptor (EGFR). Approaches to globally identify EGFR-binding proteins are required to elucidate this network. We affinity-purified EGFR with its interacting proteins by coprecipitation from lysates of A431 cells. A total of 183 proteins were repeatedly detected in high-resolution MS measurements. For 15 of these, direct interactions with EGFR were listed in the iRefIndex interaction database, including Grb2, shc-1, SOS1 and 2, STAT 1 and 3, AP2, UBS3B, and ERRFI. The newly developed Cytoscape plugin ModuleGraph allowed retrieving and visualizing 93 well-described protein complexes that contained at least one of the proteins found to interact with EGFR in our experiments. Abundances of 14 proteins were modulated more than twofold upon EGFR activation whereof clathrin-associated adaptor complex AP-2 showed 4.6-fold enrichment. These proteins were further annotated with different cellular compartments. Finally, interactions of AP-2 proteins and the newly discovered interaction of CIP2A could be verified. In conclusion, a powerful technique is presented that allowed identification and quantitative assessment of the EGFR interactome to provide further insight into EGFR signaling.

Novel search method for the discovery of functional relationships.

MOTIVATION: Numerous annotations are available that functionally characterize genes and proteins with regard to molecular process, cellular localization, tissue expression, protein domain composition, protein interaction, disease association and other properties. Searching this steadily growing amount of information can lead to the discovery of new biological relationships between genes and proteins. To facilitate the searches, methods are required that measure the annotation similarity of genes and proteins. However, most current similarity methods are focused only on annotations from the Gene Ontology (GO) and do not take other annotation sources into account. RESULTS: We introduce the new method BioSim that incorporates multiple sources of annotations to quantify the functional similarity of genes and proteins. We compared the performance of our method with four other well-known methods adapted to use multiple annotation sources. We evaluated the methods by searching for known functional relationships using annotations based only on GO or on our large data warehouse BioMyn. This warehouse integrates many diverse annotation sources of human genes and proteins. We observed that the search performance improved substantially for almost all methods when multiple annotation sources were included. In particular, our method outperformed the other methods in terms of recall and average precision.

Inferring physical protein contacts from large-scale purification data of protein complexes.

Recent large-scale data sets of protein complex purifications have provided unprecedented insights into the organization of cellular protein complexes. Several computational methods have been developed to detect co-complexed proteins in these data sets. Their common aim is the identification of biologically relevant protein complexes. However, much less is known about the network of direct physical protein contacts within the detected protein complexes. Therefore, our work investigates whether direct physical contacts can be computationally derived by combining raw data of large-scale protein complex purifications. We assess four established scoring schemes and introduce a new scoring approach that is specifically devised to infer direct physical protein contacts from protein complex purifications. The physical contacts identified by the five methods are comprehensively benchmarked against different reference sets that provide evidence for true physical contacts. Our results show that raw purification data can indeed be exploited to determine high-confidence physical protein contacts within protein complexes. In particular, our new method outperforms competing approaches at discovering physical contacts involving proteins that have been screened multiple times in purification experiments. It also excels in the analysis of recent protein purification screens of molecular chaperones and protein kinases. In contrast to previous findings, we observe that physical contacts inferred from purification experiments of protein complexes can be qualitatively comparable to binary protein interactions measured by experimental high-throughput assays such as yeast two-hybrid. This suggests that computationally derived physical contacts might complement binary protein interaction assays and guide large-scale interactome mapping projects by prioritizing putative physical contacts for further experimental screens.

Genetic variation in DLG5 is associated with inflammatory bowel disease.

Crohn disease and ulcerative colitis are two subphenotypes of inflammatory bowel disease (IBD), a complex disorder resulting from gene-environment interaction. We refined our previously defined linkage region for IBD on chromosome 10q23 and used positional cloning to identify genetic variants in DLG5 associated with IBD. DLG5 encodes a scaffolding protein involved in the maintenance of epithelial integrity. We identified two distinct haplotypes with a replicable distortion in transmission (P = 0.000023 and P = 0.004 for association with IBD, P = 0.00012 and P = 0.04 for association with Crohn disease). One of the risk-associated DLG5 haplotypes is distinguished from the common haplotype by a nonsynonymous single-nucleotide polymorphism 113G-->A, resulting in the amino acid substitution R30Q in the DUF622 domain of DLG5. This mutation probably impedes scaffolding of DLG5. We stratified the study sample according to the presence of risk-associated CARD15 variants to study potential gene-gene interaction. We found a significant difference in association of the 113A DLG5 variant with Crohn disease in affected individuals carrying the risk-associated CARD15 alleles versus those carrying non-risk-associated CARD15 alleles. This is suggestive of a complex pattern of gene-gene interaction between DLG5 and CARD15, reflecting the complex nature of polygenic diseases. Further functional studies will evaluate the biological significance of DLG5 variants.

FunSimMat update: new features for exploring functional similarity.

Quantifying the functional similarity of genes and their products based on Gene Ontology annotation is an important tool for diverse applications like the analysis of gene expression data, the prediction and validation of protein functions and interactions, and the prioritization of disease genes. The Functional Similarity Matrix (FunSimMat, http://www.funsimmat.de) is a comprehensive database providing various precomputed functional similarity values for proteins in UniProtKB and for protein families in Pfam and SMART. With this update, we significantly increase the coverage of FunSimMat by adding data from the Gene Ontology Annotation project as well as new functional similarity measures. The applicability of the database is greatly extended by the implementation of a new Gene Ontology-based method for disease gene prioritization. Two new visualization tools allow an interactive analysis of the functional relationships between proteins or protein families. This is enhanced further by the introduction of an automatically derived hierarchy of annotation classes. Additional changes include a revised user front-end and a new RESTlike interface for improving the user-friendliness and online accessibility of FunSimMat.

AltAnalyze and DomainGraph: analyzing and visualizing exon expression data.

Alternative splicing is an important mechanism for increasing protein diversity. However, its functional effects are largely unknown. Here, we present our new software workflow composed of the open-source application AltAnalyze and the Cytoscape plugin DomainGraph. Both programs provide an intuitive and comprehensive end-to-end solution for the analysis and visualization of alternative splicing data from Affymetrix Exon and Gene Arrays at the level of proteins, domains, microRNA binding sites, molecular interactions and pathways. Our software tools include easy-to-use graphical user interfaces, rigorous statistical methods (FIRMA, MiDAS and DABG filtering) and do not require prior knowledge of exon array analysis or programming. They provide new methods for automatic interpretation and visualization of the effects of alternative exon inclusion on protein domain composition and microRNA binding sites. These data can be visualized together with affected pathways and gene or protein interaction networks, allowing a straightforward identification of potential biological effects due to alternative splicing at different levels of granularity. Our programs are available at http://www.altanalyze.org and http://www.domaingraph.de. These websites also include extensive documentation, tutorials and sample data.

Tissue-specific proteins and functional implications.

Tissue-specific gene expression can result in the presence or absence of certain protein interactions and complexes, leading to profound functional differences of biological processes between the tissues. In this study, we integrate human gene expression data based on RNA-sequencing with protein interactions, domains and complexes to analyze the functional implications of their tissue specificity. This reveals that tissue specificity is characterized by much fewer proteins, domains, interactions, and complexes than previously thought. In contrast to previous microarray studies, our analysis based on RNA-sequencing suggests that tissue-specific protein interactions are less common and mainly involved with transmembrane transport and receptor activation. Additionally, tissue-specific protein domains show enrichments in DNA-related functions. This confirms that receptor-activated signaling processes and transcriptional regulation are two key factors for tissue specificity. Furthermore, many protein complexes are widely expressed regardless of their size, and their formation is frequently controlled by very few tissue-specific proteins. Interestingly, the number of alternative transcripts is increased for widely expressed genes. This suggests that alternative splicing plays a prominent role in generating specific functional characteristics of tissues.

Intragenic allele pyramiding combines different specificities of wheat Pm3 resistance alleles.

Some plant resistance genes occur as allelic series, with each member conferring specific resistance against a subset of pathogen races. In wheat, there are 17 alleles of the Pm3 gene. They encode nucleotide-binding (NB-ARC) and leucine-rich-repeat (LRR) domain proteins, which mediate resistance to distinct race spectra of powdery mildew. It is not known if specificities from different alleles can be combined to create resistance genes with broader specificity. Here, we used an approach based on avirulence analysis of pathogen populations to characterize the molecular basis of Pm3 recognition spectra. A large survey of mildew races for avirulence on the Pm3 alleles revealed that Pm3a has a resistance spectrum that completely contains that of Pm3f, but also extends towards additional races. The same is true for the Pm3b and Pm3c gene pair. The molecular analysis of these allelic pairs revealed a role of the NB-ARC protein domain in the efficiency of effector-dependent resistance. Analysis of the wild-type and chimeric Pm3 alleles identified single residues in the C-terminal LRR motifs as the main determinant of allele specificity. Variable residues of the N-terminal LRRs are necessary, but not sufficient, to confer resistance specificity. Based on these data, we constructed a chimeric Pm3 gene by intragenic allele pyramiding of Pm3d and Pm3e that showed the combined resistance specificity and, thus, a broader recognition spectrum compared with the parental alleles. Our findings support a model of stepwise evolution of Pm3 recognition specificities.

Improving disease gene prioritization using the semantic similarity of Gene Ontology terms.

MOTIVATION: Many hereditary human diseases are polygenic, resulting from sequence alterations in multiple genes. Genomic linkage and association studies are commonly performed for identifying disease-related genes. Such studies often yield lists of up to several hundred candidate genes, which have to be prioritized and validated further. Recent studies discovered that genes involved in phenotypically similar diseases are often functionally related on the molecular level. RESULTS: Here, we introduce MedSim, a novel approach for ranking candidate genes for a particular disease based on functional comparisons involving the Gene Ontology. MedSim uses functional annotations of known disease genes for assessing the similarity of diseases as well as the disease relevance of candidate genes. We benchmarked our approach with genes known to be involved in 99 diseases taken from the OMIM database. Using artificial quantitative trait loci, MedSim achieved excellent performance with an area under the ROC curve of up to 0.90 and a sensitivity of over 70% at 90% specificity when classifying gene products according to their disease relatedness. This performance is comparable or even superior to related methods in the field, albeit using less and thus more easily accessible information. AVAILABILITY: MedSim is offered as part of our FunSimMat web service (http://www.funsimmat.de).

DASMIweb: online integration, analysis and assessment of distributed protein interaction data.

In recent years, we have witnessed a substantial increase of the amount of available protein interaction data. However, most data are currently not readily accessible to the biologist at a single site, but scattered over multiple online repositories. Therefore, we have developed the DASMIweb server that affords the integration, analysis and qualitative assessment of distributed sources of interaction data in a dynamic fashion. Since DASMIweb allows for querying many different resources of protein and domain interactions simultaneously, it serves as an important starting point for interactome studies and assists the user in finding publicly accessible interaction data with minimal effort. The pool of queried resources is fully configurable and supports the inclusion of own interaction data or confidence scores. In particular, DASMIweb integrates confidence measures like functional similarity scores to assess individual interactions. The retrieved results can be exported in different file formats like MITAB or SIF. DASMIweb is freely available at http://www.dasmiweb.de.

DASMI: exchanging, annotating and assessing molecular interaction data.

MOTIVATION: Ever increasing amounts of biological interaction data are being accumulated worldwide, but they are currently not readily accessible to the biologist at a single site. New techniques are required for retrieving, sharing and presenting data spread over the Internet. RESULTS: We introduce the DASMI system for the dynamic exchange, annotation and assessment of molecular interaction data. DASMI is based on the widely used Distributed Annotation System (DAS) and consists of a data exchange specification, web servers for providing the interaction data and clients for data integration and visualization. The decentralized architecture of DASMI affords the online retrieval of the most recent data from distributed sources and databases. DASMI can also be extended easily by adding new data sources and clients. We describe all DASMI components and demonstrate their use for protein and domain interactions. AVAILABILITY: The DASMI tools are available at http://www.dasmi.de/ and http://ipfam.sanger.ac.uk/graph. The DAS registry and the DAS 1.53E specification is found at http://www.dasregistry.org/.

Variants in a novel epidermal collagen gene (COL29A1) are associated with atopic dermatitis.

Atopic dermatitis (AD) is a common chronic inflammatory skin disorder and a major manifestation of allergic disease. AD typically presents in early childhood often preceding the onset of an allergic airway disease, such as asthma or hay fever. We previously mapped a susceptibility locus for AD on Chromosome 3q21. To identify the underlying disease gene, we used a dense map of microsatellite markers and single nucleotide polymorphisms, and we detected association with AD. In concordance with the linkage results, we found a maternal transmission pattern. Furthermore, we demonstrated that the same families contribute to linkage and association. We replicated the association and the maternal effect in a large independent family cohort. A common haplotype showed strong association with AD (p = 0.000059). The associated region contained a single gene, COL29A1, which encodes a novel epidermal collagen. COL29A1 shows a specific gene expression pattern with the highest transcript levels in skin, lung, and the gastrointestinal tract, which are the major sites of allergic disease manifestation. Lack of COL29A1 expression in the outer epidermis of AD patients points to a role of collagen XXIX in epidermal integrity and function, the breakdown of which is a clinical hallmark of AD.

FunSimMat: a comprehensive functional similarity database.

Functional similarity based on Gene Ontology (GO) annotation is used in diverse applications like gene clustering, gene expression data analysis, protein interaction prediction and evaluation. However, there exists no comprehensive resource of functional similarity values although such a database would facilitate the use of functional similarity measures in different applications. Here, we describe FunSimMat (Functional Similarity Matrix, http://funsimmat.bioinf.mpi-inf.mpg.de/), a large new database that provides several different semantic similarity measures for GO terms. It offers various precomputed functional similarity values for proteins contained in UniProtKB and for protein families in Pfam and SMART. The web interface allows users to efficiently perform both semantic similarity searches with GO terms and functional similarity searches with proteins or protein families. All results can be downloaded in tab-delimited files for use with other tools. An additional XML-RPC interface gives automatic online access to FunSimMat for programs and remote services.