Combining guilt-by-association and guilt-by-profiling to predict Saccharomyces cerevisiae gene function

Fibrinaloid Microclots and Atrial Fibrillation

Atrial fibrillation (AF) is a comorbidity of a variety of other chronic, inflammatory diseases for which fibrinaloid microclots are a known accompaniment (and in some cases, a cause, with a mechanistic basis). Clots are, of course, a well-known consequence of atrial fibrillation. We here ask the question whether the fibrinaloid microclots seen in plasma or serum may in fact also be a cause of (or contributor to) the development of AF. We consider known 'risk factors' for AF, and in particular, exogenous stimuli such as infection and air pollution by particulates, both of which are known to cause AF. The external accompaniments of both bacterial (lipopolysaccharide and lipoteichoic acids) and viral (SARS-CoV-2 spike protein) infections are known to stimulate fibrinaloid microclots when added in vitro, and fibrinaloid microclots, as with other amyloid proteins, can be cytotoxic, both by inducing hypoxia/reperfusion and by other means. Strokes and thromboembolisms are also common consequences of AF. Consequently, taking a systems approach, we review the considerable evidence in detail, which leads us to suggest that it is likely that microclots may well have an aetiological role in the development of AF. This has significant mechanistic and therapeutic implications.

Transcriptomic analysis reveals mitochondrial pathways associated with distinct adolescent behavioral phenotypes and stress response

Adolescent individuals exhibit great variability in cortical dynamics and behavioral outcomes. The developing adolescent brain is highly sensitive to social experiences and environmental insults, influencing how personality traits emerge. A distinct pattern of mitochondrial gene expression in the prefrontal cortex (PFC) during adolescence underscores the essential role of mitochondria in brain maturation and the development of mental illnesses. Mitochondrial features in certain brain regions account for behavioral differences in adulthood. However, it remains unclear whether distinct adolescent behavioral phenotypes and the behavioral consequences of early adolescent stress exposure in rats are accompanied by changes in PFC mitochondria-related genes and mitochondria respiratory chain capacity. We performed a behavioral characterization during late adolescence (postnatal day, PND 47-50), including naïve animals and a group exposed to stress from PND 31-40 (10 days of footshock and 3 restraint sessions) by z-normalized data from three behavioral domains: anxiety (light-dark box tests), sociability (social interaction test) and cognition (novel-object recognition test). Employing principal component analysis, we identified three clusters: naïve with higher-behavioral z-score (HBZ), naïve with lower-behavioral z-score (LBZ), and stressed animals. Genome-wide transcriptional profiling unveiled differences in the expression of mitochondria-related genes in both naïve LBZ and stressed animals compared to naïve HBZ. Genes encoding subunits of oxidative phosphorylation complexes were significantly down-regulated in both naïve LBZ and stressed animals and positively correlated with behavioral z-score of phenotypes. Our network topology analysis of mitochondria-associated genes found Ndufa10 and Cox6a1 genes as central identifiers for naïve LBZ and stressed animals, respectively. Through high-resolution respirometry analysis, we found that both naïve LBZ and stressed animals exhibited a reduced prefrontal phosphorylation capacity and redox dysregulation. Our findings identify an association between mitochondrial features and distinct adolescent behavioral phenotypes while also underscoring the detrimental functional consequences of adolescent stress on the PFC.

Ferroptosis-related genes as diagnostic markers for major depressive disorder and their correlations with immune infiltration

Background

Major depression disorder (MDD) is a devastating neuropsychiatric disease, and one of the leading causes of suicide. Ferroptosis, an iron-dependent form of regulated cell death, plays a pivotal role in numerous diseases. The study aimed to construct and validate a gene signature for diagnosing MDD based on ferroptosis-related genes (FRGs) and further explore the biological functions of these genes in MDD.

Methods

The datasets were downloaded from the Gene Expression Omnibus (GEO) database and FRGs were obtained from the FerrDb database and other literatures. Least absolute shrinkage and selection operator (LASSO) regression and stepwise logistic regression were performed to develop a gene signature. Receiver operating characteristic (ROC) curves were utilized to assess the diagnostic power of the signature. Gene ontology (GO) enrichment analysis was used to explore the biological roles of these diagnostic genes, and single sample gene set enrichment analysis (ssGSEA) algorithm was used to evaluate immune infiltration in MDD. Animal model of depression was constructed to validate the expression of the key genes.

Results

Eleven differentially expressed FRGs were identified in MDD patients compared with healthy controls. A signature of three FRGs (ALOX15B, RPLP0, and HP) was constructed for diagnosis of MDD. Afterwards, ROC analysis confirmed the signature's discriminative capacity (AUC = 0.783, 95% CI = 0.719-0.848). GO enrichment analysis revealed that the differentially expressed genes (DEGs) related to these three FRGs were mainly involved in immune response. Furthermore, spearman correlation analysis demonstrated that these three FRGs were associated with infiltrating immune cells. ALOX15B and HP were significantly upregulated and RPLP0 was significantly downregulated in peripheral blood of the lipopolysaccharide (LPS)-induced depressive model.

Conclusion

Our results suggest that the novel FRG signature had a good diagnostic performance for MDD, and these three FRGs correlated with immune infiltration in MDD.

DapBCH: a disease association prediction model Based on Cross-species and Heterogeneous graph embedding

The study of comorbidity can provide new insights into the pathogenesis of the disease and has important economic significance in the clinical evaluation of treatment difficulty, medical expenses, length of stay, and prognosis of the disease. In this paper, we propose a disease association prediction model DapBCH, which constructs a cross-species biological network and applies heterogeneous graph embedding to predict disease association. First, we combine the human disease-gene network, mouse gene-phenotype network, human-mouse homologous gene network, and human protein-protein interaction network to reconstruct a heterogeneous biological network. Second, we apply heterogeneous graph embedding based on meta-path aggregation to generate the feature vector of disease nodes. Finally, we employ link prediction to obtain the similarity of disease pairs. The experimental results indicate that our model is highly competitive in predicting the disease association and is promising for finding potential disease associations.

Identification of hypoxia-related genes and exploration of their relationship with immune cells in ischemic stroke

Ischemic stroke (IS) is a major threat to human health, and it is the second leading cause of long-term disability and death in the world. Impaired cerebral perfusion leads to acute hypoxia and glucose deficiency, which in turn induces a stroke cascade response that ultimately leads to cell death. Screening and identifying hypoxia-related genes (HRGs) and therapeutic targets is important for neuroprotection before and during brain recanalization to protect against injury and extend the time window to further improve functional outcomes before pharmacological and mechanical thrombolysis. First, we downloaded the GSE16561 and GSE58294 datasets from the NCBI GEO database. Bioinformatics analysis of the GSE16561 dataset using the limma package identified differentially expressed genes (DEGs) in ischemic stroke using adj. p. values < 0.05 and a fold change of 0.5 as thresholds. The Molecular Signature database and Genecards database were pooled to obtain hypoxia-related genes. 19 HRGs associated with ischemic stroke were obtained after taking the intersection. LASSO regression and multivariate logistic regression were applied to identify critical biomarkers with independent diagnostic values. ROC curves were constructed to validate their diagnostic efficacy. We used CIBERSORT to analyze the differences in the immune microenvironment between IS patients and controls. Finally, we investigated the correlation between HRGs and infiltrating immune cells to understand molecular immune mechanisms better. Our study analyzed the role of HRGs in ischemic stroke. Nineteen hypoxia-related genes were obtained. Enrichment analysis showed that 19 HRGs were involved in response to hypoxia, HIF-1 signaling pathway, autophagy, autophagy of mitochondrion, and AMPK signaling pathway. Because of the good diagnostic properties of SLC2A3, we further investigated the function of SLC2A3 and found that it is closely related to immunity. We have also explored the relevance of other critical genes to immune cells. Our findings suggest that hypoxia-related genes play a crucial role in the diversity and complexity of the IS immune microenvironment. Exploring the association between hypoxia-related critical genes and immune cells provides innovative insights into the therapeutic targets for ischemic stroke.

SUBATOMIC: a SUbgraph BAsed mulTi-OMIcs clustering framework to analyze integrated multi-edge networks

BACKGROUND

Representing the complex interplay between different types of biomolecules across different omics layers in multi-omics networks bears great potential to gain a deep mechanistic understanding of gene regulation and disease. However, multi-omics networks easily grow into giant hairball structures that hamper biological interpretation. Module detection methods can decompose these networks into smaller interpretable modules. However, these methods are not adapted to deal with multi-omics data nor consider topological features. When deriving very large modules or ignoring the broader network context, interpretability remains limited. To address these issues, we developed a SUbgraph BAsed mulTi-OMIcs Clustering framework (SUBATOMIC), which infers small and interpretable modules with a specific topology while keeping track of connections to other modules and regulators.

RESULTS

SUBATOMIC groups specific molecular interactions in composite network subgraphs of two and three nodes and clusters them into topological modules. These are functionally annotated, visualized and overlaid with expression profiles to go from static to dynamic modules. To preserve the larger network context, SUBATOMIC investigates statistically the connections in between modules as well as between modules and regulators such as miRNAs and transcription factors. We applied SUBATOMIC to analyze a composite Homo sapiens network containing transcription factor-target gene, miRNA-target gene, protein-protein, homologous and co-functional interactions from different databases. We derived and annotated 5586 modules with diverse topological, functional and regulatory properties. We created novel functional hypotheses for unannotated genes. Furthermore, we integrated modules with condition specific expression data to study the influence of hypoxia in three cancer cell lines. We developed two prioritization strategies to identify the most relevant modules in specific biological contexts: one considering GO term enrichments and one calculating an activity score reflecting the degree of differential expression. Both strategies yielded modules specifically reacting to low oxygen levels.

CONCLUSIONS

We developed the SUBATOMIC framework that generates interpretable modules from integrated multi-omics networks and applied it to hypoxia in cancer. SUBATOMIC can infer and contextualize modules, explore condition or disease specific modules, identify regulators and functionally related modules, and derive novel gene functions for uncharacterized genes. The software is available at https://github.com/CBIGR/SUBATOMIC .

Identification and functional annotation of long intergenic non-coding RNAs in Brassicaceae

Long intergenic noncoding RNAs (lincRNAs) are a large yet enigmatic class of eukaryotic transcripts that can have critical biological functions. The wealth of RNA-sequencing (RNA-seq) data available for plants provides the opportunity to implement a harmonized identification and annotation effort for lincRNAs that enables cross-species functional and genomic comparisons as well as prioritization of functional candidates. In this study, we processed >24 Tera base pairs of RNA-seq data from >16,000 experiments to identify ∼130,000 lincRNAs in four Brassicaceae: Arabidopsis thaliana, Camelina sativa, Brassica rapa, and Eutrema salsugineum. We used nanopore RNA-seq, transcriptome-wide structural information, peptide data, and epigenomic data to characterize these lincRNAs and identify conserved motifs. We then used comparative genomic and transcriptomic approaches to highlight lincRNAs in our data set with sequence or transcriptional conservation. Finally, we used guilt-by-association analyses to assign putative functions to lincRNAs within our data set. We tested this approach on a subset of lincRNAs associated with germination and seed development, observing germination defects for Arabidopsis lines harboring T-DNA insertions at these loci. LincRNAs with Brassicaceae-conserved putative miRNA binding motifs, small open reading frames, or abiotic-stress modulated expression are a few of the annotations that will guide functional analyses into this cryptic portion of the transcriptome.

Signatures of 4 autophagy-related genes as diagnostic markers of MDD and their correlation with immune infiltration

BACKGROUND

Major depressive disorder (MDD) is a debilitating mental illness and one of the primary causes of suicide. This study attempted to develop and validate a multigene joint signature for diagnosing MDD based on autophagy-related genes (ARGs) and to explore their biological role in MDD.

METHODS

We downloaded data from the Gene Expression Omnibus (GEO) database and retrieved ARGs from the Human Autophagy Database. The limma package in R software was used to identify differentially expressed genes (DEGs). We used CIBERSORT to analyze differences in the immune microenvironment between MDD patients and controls. Finally, we examined the correlation between diagnostic markers and infiltrating immune cells to better understand the molecular immune mechanism.

RESULTS

In this study, we identified 20 differentially expressed ARGs in MDD compared to controls. A signature of 4 autophagy-related genes (GPR18, PDK4, NRG1 and EPHB2) was obtained. ROC analysis showed that our model has good diagnostic performance (AUC=0.779, 95% CI=0.709-0.848). Bioinformatics analysis validated that GPR18 may represent a new candidate gene for MDD. Correlation analysis revealed that GPR18 was positively correlated with regulatory T cells (Treg), CD8+ T cells, naive B cells, and memory B cells and negatively correlated with M0 macrophages and neutrophils in MDD.

LIMITATIONS

This was a second mining of previously published data sets. Independent studies are warranted to validate and improve the clinical utility of the identified signature.

CONCLUSIONS

We identified a novel four-ARG gene signature that has good diagnostic performance and identified an association between ARG genes and the immune microenvironment in MDD.

IID 2021: towards context-specific protein interaction analyses by increased coverage, enhanced annotation and enrichment analysis

Improved bioassays have significantly increased the rate of identifying new protein-protein interactions (PPIs), and the number of detected human PPIs has greatly exceeded early estimates of human interactome size. These new PPIs provide a more complete view of disease mechanisms but precise understanding of how PPIs affect phenotype remains a challenge. It requires knowledge of PPI context (e.g. tissues, subcellular localizations), and functional roles, especially within pathways and protein complexes. The previous IID release focused on PPI context, providing networks with comprehensive tissue, disease, cellular localization, and druggability annotations. The current update adds developmental stages to the available contexts, and provides a way of assigning context to PPIs that could not be previously annotated due to insufficient data or incompatibility with available context categories (e.g. interactions between membrane and cytoplasmic proteins). This update also annotates PPIs with conservation across species, directionality in pathways, membership in large complexes, interaction stability (i.e. stable or transient), and mutation effects. Enrichment analysis is now available for all annotations, and includes multiple options; for example, context annotations can be analyzed with respect to PPIs or network proteins. In addition to tabular view or download, IID provides online network visualization. This update is available at http://ophid.utoronto.ca/iid.

Network-based methods for gene function prediction

The rapid development of high-throughput technology has generated a large number of biological networks. Network-based methods are able to provide rich information for inferring gene function. This is composed of analyzing the topological characteristics of genes in related networks, integrating biological information, and considering data from different data sources. To promote network biology and related biotechnology research, this article provides a survey for the state of the art of advanced methods of network-based gene function prediction and discusses the potential challenges.

The STRING database in 2021: customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets

Cellular life depends on a complex web of functional associations between biomolecules. Among these associations, protein–protein interactions are particularly important due to their versatility, specificity and adaptability. The STRING database aims to integrate all known and predicted associations between proteins, including both physical interactions as well as functional associations. To achieve this, STRING collects and scores evidence from a number of sources: (i) automated text mining of the scientific literature, (ii) databases of interaction experiments and annotated complexes/pathways, (iii) computational interaction predictions from co-expression and from conserved genomic context and (iv) systematic transfers of interaction evidence from one organism to another. STRING aims for wide coverage; the upcoming version 11.5 of the resource will contain more than 14 000 organisms. In this update paper, we describe changes to the text-mining system, a new scoring-mode for physical interactions, as well as extensive user interface features for customizing, extending and sharing protein networks. In addition, we describe how to query STRING with genome-wide, experimental data, including the automated detection of enriched functionalities and potential biases in the user's query data. The STRING resource is available online, at https://string-db.org/.

Identifying functionally relevant candidate genes for inflexible ethanol intake in mice and humans using a guilt-by-association approach

Gene prioritization approaches are useful tools to explore and select candidate genes in transcriptome studies. Knowing the importance of processes such as neuronal activity, intracellular signal transduction, and synapse plasticity to the development and maintenance of compulsive ethanol drinking, the aim of the present study was to explore and identify functional candidate genes associated with these processes in an animal model of inflexible pattern of ethanol intake. To do this, we applied a guilt-by-association approach, using the GUILDify and ToppGene software, in our previously published microarray data from the prefrontal cortex (PFC) and striatum of inflexible drinker mice. We then tested some of the prioritized genes that showed a tissue-specific pattern in postmortem brain tissue (PFC and nucleus accumbens (NAc)) from humans with alcohol use disorder (AUD). In the mouse brain, we prioritized 44 genes in PFC and 26 in striatum, which showed opposite regulation patterns in PFC and striatum. The most prioritized of them (i.e., Plcb1 and Prkcb in PFC, and Dnm2 and Lrrk2 in striatum) were associated with synaptic neuroplasticity, a neuroadaptation associated with excessive ethanol drinking. The identification of transcription factors among the prioritized genes suggests a crucial role for Irf4 in the pattern of regulation observed between PFC and striatum. Lastly, the differential transcription of IRF4 and LRRK2 in PFC and nucleus accumbens in postmortem brains from AUD compared to control highlights their involvement in compulsive ethanol drinking in humans and mice.

Natural display of nuclear-encoded RNA on the cell surface and its impact on cell interaction

BACKGROUND

Compared to proteins, glycans, and lipids, much less is known about RNAs on the cell surface. We develop a series of technologies to test for any nuclear-encoded RNAs that are stably attached to the cell surface and exposed to the extracellular space, hereafter called membrane-associated extracellular RNAs (maxRNAs).

RESULTS

We develop a technique called Surface-seq to selectively sequence maxRNAs and validate two Surface-seq identified maxRNAs by RNA fluorescence in situ hybridization. To test for cell-type specificity of maxRNA, we use antisense oligos to hybridize to single-stranded transcripts exposed on the surface of human peripheral blood mononuclear cells (PBMCs). Combining this strategy with imaging flow cytometry, single-cell RNA sequencing, and maxRNA sequencing, we identify monocytes as the major type of maxRNA+ PBMCs and prioritize 11 candidate maxRNAs for functional tests. Extracellular application of antisense oligos of FNDC3B and CTSS transcripts inhibits monocyte adhesion to vascular endothelial cells.

CONCLUSIONS

Collectively, these data highlight maxRNAs as functional components of the cell surface, suggesting an expanded role for RNA in cell-cell and cell-environment interactions.

Hierarchy and levels: analysing networks to study mechanisms in molecular biology

Network representations are flat while mechanisms are organized into a hierarchy of levels, suggesting that the two are fundamentally opposed. I challenge this opposition by focusing on two aspects of the ways in which large-scale networks constructed from high-throughput data are analysed in systems biology: identifying clusters of nodes that operate as modules or mechanisms and using bio-ontologies such as gene ontology (GO) to annotate nodes with information about where entities appear in cells and the biological functions in which they participate. Of particular importance, GO organizes biological knowledge about cell components and functions hierarchically. I illustrate how this supports mechanistic interpretation of networks with two examples of network studies, one using epistatic interactions among genes to identify mechanisms and their parts and the other using deep learning to predict phenotypes. As illustrated in these examples, when network research draws upon hierarchical information such as provided by GO, the results not only can be interpreted mechanistically but provide new mechanistic knowledge. This article is part of the theme issue 'Unifying the essential concepts of biological networks: biological insights and philosophical foundations'.

ISOGO: Functional annotation of protein-coding splice variants

The advent of RNA-seq technologies has switched the paradigm of genetic analysis from a genome to a transcriptome-based perspective. Alternative splicing generates functional diversity in genes, but the precise functions of many individual isoforms are yet to be elucidated. Gene Ontology was developed to annotate gene products according to their biological processes, molecular functions and cellular components. Despite a single gene may have several gene products, most annotations are not isoform-specific and do not distinguish the functions of the different proteins originated from a single gene. Several approaches have tried to automatically annotate ontologies at the isoform level, but this has shown to be a daunting task. We have developed ISOGO (ISOform + GO function imputation), a novel algorithm to predict the function of coding isoforms based on their protein domains and their correlation of expression along 11,373 cancer patients. Combining these two sources of information outperforms previous approaches: it provides an area under precision-recall curve (AUPRC) five times larger than previous attempts and the median AUROC of assigned functions to genes is 0.82. We tested ISOGO predictions on some genes with isoform-specific functions (BRCA1, MADD,VAMP7 and ITSN1) and they were coherent with the literature. Besides, we examined whether the main isoform of each gene -as predicted by APPRIS- was the most likely to have the annotated gene functions and it occurs in 99.4% of the genes. We also evaluated the predictions for isoform-specific functions provided by the CAFA3 challenge and results were also convincing. To make these results available to the scientific community, we have deployed a web application to consult ISOGO predictions (https://biotecnun.unav.es/app/isogo). Initial data, website link, isoform-specific GO function predictions and R code is available at https://gitlab.com/icassol/isogo.

The Genetic Perspective of Familial Glucocorticoid Deficiency: In Silico Analysis of Two Novel Variants

Familial glucocorticoid deficiency is a rare autosomal recessive genetic disorder which belongs to a group of primary adrenal insufficiency (PAI) and is mainly caused by mutations in the MC2R and MRAP genes. A comprehensive search was conducted to find the reported variants of MC2R and MRAP genes. In silico pathogenic analysis was performed for the reported variants. PCR amplification and sequencing were performed for three patients. Structural analysis, modeling, and interactome analysis were applied to characterize novel MC2R variants and their proteins. About 80% of MC2R-related cases showed the clinical symptoms which were diagnosed at <2 years old. 107 patients had MC2R mutations (85 homozygotes, 21 compound heterozygotes, and 1 simple heterozygote). 59 variants were found in the MC2R gene. Four mutations were responsible for half of patients. 39 homozygous patients had MRAP mutations; 14 variants were determined in the MRAP gene. Nine proteins were predicted by STRING to associate with the studied proteins. Two novel MC2R variants, c.128T > G (p.Leu43Arg) and c.251T > A (p.Ile84Asn), were found in two patients at the age of above and below 2 years, respectively. Mutations in MC2R and MRAP genes are the main cause of FGD. Genetic studies and in silico analysis will help to confirm the diagnosis.

Unified feature association networks through integration of transcriptomic and proteomic data

High-throughput multi-omics studies and corresponding network analyses of multi-omic data have rapidly expanded their impact over the last 10 years. As biological features of different types (e.g. transcripts, proteins, metabolites) interact within cellular systems, the greatest amount of knowledge can be gained from networks that incorporate multiple types of -omic data. However, biological and technical sources of variation diminish the ability to detect cross-type associations, yielding networks dominated by communities comprised of nodes of the same type. We describe here network building methods that can maximize edges between nodes of different data types leading to integrated networks, networks that have a large number of edges that link nodes of different-omic types (transcripts, proteins, lipids etc). We systematically rank several network inference methods and demonstrate that, in many cases, using a random forest method, GENIE3, produces the most integrated networks. This increase in integration does not come at the cost of accuracy as GENIE3 produces networks of approximately the same quality as the other network inference methods tested here. Using GENIE3, we also infer networks representing antibody-mediated Dengue virus cell invasion and receptor-mediated Dengue virus invasion. A number of functional pathways showed centrality differences between the two networks including genes responding to both GM-CSF and IL-4, which had a higher centrality value in an antibody-mediated vs. receptor-mediated Dengue network. Because a biological system involves the interplay of many different types of molecules, incorporating multiple data types into networks will improve their use as models of biological systems. The methods explored here are some of the first to specifically highlight and address the challenges associated with how such multi-omic networks can be assembled and how the greatest number of interactions can be inferred from different data types. The resulting networks can lead to the discovery of new host response patterns and interactions during viral infection, generate new hypotheses of pathogenic mechanisms and confirm mechanisms of disease.

The EXPANDER Integrated Platform for Transcriptome Analysis

Genome-wide analysis of cellular transcriptomes using RNA-seq or expression arrays is a major mainstay of current biological and biomedical research. EXPANDER (EXPression ANalyzer and DisplayER) is a comprehensive software package for analysis of expression data, with built-in support for 18 different organisms. It is designed as a "one-stop shop" platform for transcriptomic analysis, allowing for execution of all analysis steps starting with gene expression data matrix. Analyses offered include low-level preprocessing and normalization, differential expression analysis, clustering, bi-clustering, supervised grouping, high-level functional and pathway enrichment tests, and networks and motif analyses. A variety of options is offered for each step, using established algorithms, including many developed and published by our laboratory. EXPANDER has been continuously developed since 2003, having to date over 18,000 downloads and 540 citations. One of the innovations in the recent version is support for combined analysis of gene expression and ChIP-seq data to enhance the inference of transcriptional networks and their functional interpretation. EXPANDER implements cutting-edge algorithms and makes them accessible to users through user-friendly interface and intuitive visualizations. It is freely available to users at http://acgt.cs.tau.ac.il/expander/.

IID 2018 update: context-specific physical protein-protein interactions in human, model organisms and domesticated species

Knowing the set of physical protein-protein interactions (PPIs) that occur in a particular context-a tissue, disease, or other condition-can provide valuable insights into key research questions. However, while the number of identified human PPIs is expanding rapidly, context information remains limited, and for most non-human species context-specific networks are completely unavailable. The Integrated Interactions Database (IID) provides one of the most comprehensive sets of context-specific human PPI networks, including networks for 133 tissues, 91 disease conditions, and many other contexts. Importantly, it also provides context-specific networks for 17 non-human species including model organisms and domesticated animals. These species are vitally important for drug discovery and agriculture. IID integrates interactions from multiple databases and datasets. It comprises over 4.8 million PPIs annotated with several types of context: tissues, subcellular localizations, diseases, and druggability information (the latter three are new annotations not available in the previous version). This update increases the number of species from 6 to 18, the number of PPIs from ∼1.5 million to ∼4.8 million, and the number of tissues from 30 to 133. IID also now supports topology and enrichment analyses of returned networks. IID is available at http://ophid.utoronto.ca/iid.

The evolutionary signal in metagenome phyletic profiles predicts many gene functions

BACKGROUND

The function of many genes is still not known even in model organisms. An increasing availability of microbiome DNA sequencing data provides an opportunity to infer gene function in a systematic manner.

RESULTS

We evaluated if the evolutionary signal contained in metagenome phyletic profiles (MPP) is predictive of a broad array of gene functions. The MPPs are an encoding of environmental DNA sequencing data that consists of relative abundances of gene families across metagenomes. We find that such MPPs can accurately predict 826 Gene Ontology functional categories, while drawing on human gut microbiomes, ocean metagenomes, and DNA sequences from various other engineered and natural environments. Overall, in this task, the MPPs are highly accurate, and moreover they provide coverage for a set of Gene Ontology terms largely complementary to standard phylogenetic profiles, derived from fully sequenced genomes. We also find that metagenomes approximated from taxon relative abundance obtained via 16S rRNA gene sequencing may provide surprisingly useful predictive models. Crucially, the MPPs derived from different types of environments can infer distinct, non-overlapping sets of gene functions and therefore complement each other. Consistently, simulations on > 5000 metagenomes indicate that the amount of data is not in itself critical for maximizing predictive accuracy, while the diversity of sampled environments appears to be the critical factor for obtaining robust models.

CONCLUSIONS

In past work, metagenomics has provided invaluable insight into ecology of various habitats, into diversity of microbial life and also into human health and disease mechanisms. We propose that environmental DNA sequencing additionally constitutes a useful tool to predict biological roles of genes, yielding inferences out of reach for existing comparative genomics approaches.

Gene co-opening network deciphers gene functional relationships

Genome sequencing technology has generated a vast amount of genomic and epigenomic data, and has provided us a great opportunity to study gene functions on a global scale from an epigenomic view. In the last decade, network-based studies, such as those based on PPI networks and co-expression networks, have shown good performance in capturing functional relationships between genes. However, the functions of a gene and the mechanism of interaction of genes with each other to elucidate their functions are still not entirely clear. Here, we construct a gene co-opening network based on chromatin accessibility of genes. We show that genes related to a specific biological process or the same disease tend to be clustered in the co-opening network. This understanding allows us to detect functional clusters from the network and to predict new functions for genes. We further apply the network to prioritize disease genes for Psoriasis, and demonstrate the power of the joint analysis of the co-opening network and GWAS data in identifying disease genes. Taken together, the co-opening network provides a new viewpoint for the elucidation of gene associations and the interpretation of disease mechanisms.

Gene function prediction based on Gene Ontology Hierarchy Preserving Hashing

Gene Ontology (GO) uses structured vocabularies (or terms) to describe the molecular functions, biological roles, and cellular locations of gene products in a hierarchical ontology. GO annotations associate genes with GO terms and indicate the given gene products carrying out the biological functions described by the relevant terms. However, predicting correct GO annotations for genes from a massive set of GO terms as defined by GO is a difficult challenge. To combat with this challenge, we introduce a Gene Ontology Hierarchy Preserving Hashing (HPHash) based semantic method for gene function prediction. HPHash firstly measures the taxonomic similarity between GO terms. It then uses a hierarchy preserving hashing technique to keep the hierarchical order between GO terms, and to optimize a series of hashing functions to encode massive GO terms via compact binary codes. After that, HPHash utilizes these hashing functions to project the gene-term association matrix into a low-dimensional one and performs semantic similarity based gene function prediction in the low-dimensional space. Experimental results on three model species (Homo sapiens, Mus musculus and Rattus norvegicus) for interspecies gene function prediction show that HPHash performs better than other related approaches and it is robust to the number of hash functions. In addition, we also take HPHash as a plugin for BLAST based gene function prediction. From the experimental results, HPHash again significantly improves the prediction performance. The codes of HPHash are available at: http://mlda.swu.edu.cn/codes.php?name=HPHash.

Proteomic Profiling of Leishmania donovani Promastigote Subcellular Organelles

To facilitate a greater understanding of the biological processes in the medically important Leishmania donovani parasite, a combination of differential and density-gradient ultracentrifugation techniques were used to achieve a comprehensive subcellular fractionation of the promastigote stage. An in-depth label-free proteomic LC-MS/MS analysis of the density gradients resulted in the identification of ∼50% of the Leishmania proteome (3883 proteins detected), which included ∼645 integral membrane proteins and 1737 uncharacterized proteins. Clustering and subcellular localization of proteins was based on a subset of training Leishmania proteins with known subcellular localizations that had been determined using biochemical, confocal microscopy, or immunoelectron microscopy approaches. This subcellular map will be a valuable resource that will help dissect the cell biology and metabolic processes associated with specific organelles of Leishmania and related kinetoplastids.

Predicting Protein Functions by Using Unbalanced Random Walk Algorithm on Three Biological Networks

With the gap between the sequence data and their functional annotations becomes increasing wider, many computational methods have been proposed to annotate functions for unknown proteins. However, designing effective methods to make good use of various biological resources is still a big challenge for researchers due to function diversity of proteins. In this work, we propose a new method named ThrRW, which takes several steps of random walking on three different biological networks: protein interaction network (PIN), domain co-occurrence network (DCN), and functional interrelationship network (FIN), respectively, so as to infer functional information from neighbors in the corresponding networks. With respect to the topological and structural differences of the three networks, the number of walking steps in the three networks will be different. In the course of working, the functional information will be transferred from one network to another according to the associations between the nodes in different networks. The results of experiment on S. cerevisiae data show that our method achieves better prediction performance not only than the methods that consider both PIN data and GO term similarities, but also than the methods using both PIN data and protein domain information, which verifies the effectiveness of our method on integrating multiple biological data sources.

The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible

A system-wide understanding of cellular function requires knowledge of all functional interactions between the expressed proteins. The STRING database aims to collect and integrate this information, by consolidating known and predicted protein–protein association data for a large number of organisms. The associations in STRING include direct (physical) interactions, as well as indirect (functional) interactions, as long as both are specific and biologically meaningful. Apart from collecting and reassessing available experimental data on protein–protein interactions, and importing known pathways and protein complexes from curated databases, interaction predictions are derived from the following sources: (i) systematic co-expression analysis, (ii) detection of shared selective signals across genomes, (iii) automated text-mining of the scientific literature and (iv) computational transfer of interaction knowledge between organisms based on gene orthology. In the latest version 10.5 of STRING, the biggest changes are concerned with data dissemination: the web frontend has been completely redesigned to reduce dependency on outdated browser technologies, and the database can now also be queried from inside the popular Cytoscape software framework. Further improvements include automated background analysis of user inputs for functional enrichments, and streamlined download options. The STRING resource is available online, at http://string-db.org/.

Extensive complementarity between gene function prediction methods

MOTIVATION

The number of sequenced genomes rises steadily but we still lack the knowledge about the biological roles of many genes. Automated function prediction (AFP) is thus a necessity. We hypothesized that AFP approaches that draw on distinct genome features may be useful for predicting different types of gene functions, motivating a systematic analysis of the benefits gained by obtaining and integrating such predictions.

RESULTS

Our pipeline amalgamates 5 133 543 genes from 2071 genomes in a single massive analysis that evaluates five established genomic AFP methodologies. While 1227 Gene Ontology (GO) terms yielded reliable predictions, the majority of these functions were accessible to only one or two of the methods. Moreover, different methods tend to assign a GO term to non-overlapping sets of genes. Thus, inferences made by diverse genomic AFP methods display a striking complementary, both gene-wise and function-wise. Because of this, a viable integration strategy is to rely on a single most-confident prediction per gene/function, rather than enforcing agreement across multiple AFP methods. Using an information-theoretic approach, we estimate that current databases contain 29.2 bits/gene of known Escherichia coli gene functions. This can be increased by up to 5.5 bits/gene using individual AFP methods or by 11 additional bits/gene upon integration, thereby providing a highly-ranking predictor on the Critical Assessment of Function Annotation 2 community benchmark. Availability of more sequenced genomes boosts the predictive accuracy of AFP approaches and also the benefit from integrating them.

AVAILABILITY AND IMPLEMENTATION

The individual and integrated GO predictions for the complete set of genes are available from http://gorbi.irb.hr/ CONTACT: fran.supek@irb.hrSupplementary information: Supplementary materials are available at Bioinformatics online.

UNIPred: Unbalance-Aware Network Integration and Prediction of Protein Functions

The proper integration of multiple sources of data and the unbalance between annotated and unannotated proteins represent two of the main issues of the automated function prediction (AFP) problem. Most of supervised and semisupervised learning algorithms for AFP proposed in literature do not jointly consider these items, with a negative impact on both sensitivity and precision performances, due to the unbalance between annotated and unannotated proteins that characterize the majority of functional classes and to the specific and complementary information content embedded in each available source of data. We propose UNIPred (unbalance-aware network integration and prediction of protein functions), an algorithm that properly combines different biomolecular networks and predicts protein functions using parametric semisupervised neural models. The algorithm explicitly takes into account the unbalance between unannotated and annotated proteins both to construct the integrated network and to predict protein annotations for each functional class. Full-genome and ontology-wide experiments with three eukaryotic model organisms show that the proposed method compares favorably with state-of-the-art learning algorithms for AFP.

Solute Carrier Family 26 Member a2 (slc26a2) Regulates Otic Development and Hair Cell Survival in Zebrafish

Hearing loss is one of the most prevalent human birth defects. Genetic factors contribute to the pathogenesis of deafness. It is estimated that one-third of deafness genes have already been identified. The current work is an attempt to find novel genes relevant to hearing loss using guilt-by-profiling and guilt-by-association bioinformatics analyses of approximately 80 known non-syndromic hereditary hearing loss (NSHL) genes. Among the 300 newly identified candidate deafness genes, slc26a2 were selected for functional studies in zebrafish. The slc26a2 gene was knocked down using an antisense morpholino (MO), and significant defects were observed in otolith patterns, semicircular canal morphology, and lateral neuromast distributions in morphants. Loss-of-function defects are caused primarily by apoptosis, and morphants are insensitive to sound stimulation and imbalanced swimming behaviours. Morphant defects were found to be partially rescued by co-injection of human SLC26A2 mRNA. All the results suggest that bioinformatics is capable of predicting new deafness genes and this showed slc26a2 is to be a critical otic gene whose dysfunction may induce hearing impairment.

Selecting causal genes from genome-wide association studies via functionally coherent subnetworks

Genome-wide association (GWA) studies have linked thousands of loci to human diseases, but the causal genes and variants at these loci generally remain unknown. Although investigators typically focus on genes closest to the associated polymorphisms, the causal gene is often more distal. Reliance on published work to prioritize candidates is biased toward well-characterized genes. We describe a 'prix fixe' strategy and software that uses genome-scale shared-function networks to identify sets of mutually functionally related genes spanning multiple GWA loci. Using associations from ∼100 GWA studies covering ten cancer types, our approach outperformed the common alternative strategy in ranking known cancer genes. As more GWA loci are discovered, the strategy will have increased power to elucidate the causes of human disease.

Prediction of disease genes using tissue-specified gene-gene network

BACKGROUND

Tissue specificity is an important aspect of many genetic diseases in the context of genetic disorders as the disorder affects only few tissues. Therefore tissue specificity is important in identifying disease-gene associations. Hence this paper seeks to discuss the impact of using tissue specificity in predicting new disease-gene associations and how to use tissue specificity along with phenotype information for a particular disease.

METHODS

In order to find out the impact of using tissue specificity for predicting new disease-gene associations, this study proposes a novel method called tissue-specified genes to construct tissues-specific gene-gene networks for different tissue samples. Subsequently, these networks are used with phenotype details to predict disease genes by using Katz method. The proposed method was compared with three other tissue-specific network construction methods in order to check its effectiveness. Furthermore, to check the possibility of using tissue-specific gene-gene network instead of generic protein-protein network at all time, the results are compared with three other methods.

RESULTS

In terms of leave-one-out cross validation, calculation of the mean enrichment and ROC curves indicate that the proposed approach outperforms existing network construction methods. Furthermore tissues-specific gene-gene networks make a more positive impact on predicting disease-gene associations than generic protein-protein interaction networks.

CONCLUSIONS

In conclusion by integrating tissue-specific data it enabled prediction of known and unknown disease-gene associations for a particular disease more effectively. Hence it is better to use tissue-specific gene-gene network whenever possible. In addition the proposed method is a better way of constructing tissue-specific gene-gene networks.

Hierarchical ensemble methods for protein function prediction

Protein function prediction is a complex multiclass multilabel classification problem, characterized by multiple issues such as the incompleteness of the available annotations, the integration of multiple sources of high dimensional biomolecular data, the unbalance of several functional classes, and the difficulty of univocally determining negative examples. Moreover, the hierarchical relationships between functional classes that characterize both the Gene Ontology and FunCat taxonomies motivate the development of hierarchy-aware prediction methods that showed significantly better performances than hierarchical-unaware "flat" prediction methods. In this paper, we provide a comprehensive review of hierarchical methods for protein function prediction based on ensembles of learning machines. According to this general approach, a separate learning machine is trained to learn a specific functional term and then the resulting predictions are assembled in a "consensus" ensemble decision, taking into account the hierarchical relationships between classes. The main hierarchical ensemble methods proposed in the literature are discussed in the context of existing computational methods for protein function prediction, highlighting their characteristics, advantages, and limitations. Open problems of this exciting research area of computational biology are finally considered, outlining novel perspectives for future research.

Novel cardiovascular gene functions revealed via systematic phenotype prediction in zebrafish

Comprehensive functional annotation of vertebrate genomes is fundamental to biological discovery. Reverse genetic screening has been highly useful for determination of gene function, but is untenable as a systematic approach in vertebrate model organisms given the number of surveyable genes and observable phenotypes. Unbiased prediction of gene-phenotype relationships offers a strategy to direct finite experimental resources towards likely phenotypes, thus maximizing de novo discovery of gene functions. Here we prioritized genes for phenotypic assay in zebrafish through machine learning, predicting the effect of loss of function of each of 15,106 zebrafish genes on 338 distinct embryonic anatomical processes. Focusing on cardiovascular phenotypes, the learning procedure predicted known knockdown and mutant phenotypes with high precision. In proof-of-concept studies we validated 16 high-confidence cardiac predictions using targeted morpholino knockdown and initial blinded phenotyping in embryonic zebrafish, confirming a significant enrichment for cardiac phenotypes as compared with morpholino controls. Subsequent detailed analyses of cardiac function confirmed these results, identifying novel physiological defects for 11 tested genes. Among these we identified tmem88a, a recently described attenuator of Wnt signaling, as a discrete regulator of the patterning of intercellular coupling in the zebrafish cardiac epithelium. Thus, we show that systematic prioritization in zebrafish can accelerate the pace of developmental gene function discovery.

Predicting gene function using similarity learning

BACKGROUND

Computational methods that make use of heterogeneous biological datasets to predict gene function provide a cost-effective and rapid way for annotating genomes. A common framework shared by many such methods is to construct a combined functional association network from multiple networks representing different sources of data, and use this combined network as input to network-based or kernel-based learning algorithms. In these methods, a key factor contributing to the prediction accuracy is the network quality, which is the ability of the network to reflect the functional relatedness of gene pairs. To improve the network quality, a large effort has been spent on developing methods for network integration. These methods, however, produce networks, which then remain unchanged, and nearly no effort has been made to optimize the networks after their construction.

RESULTS

Here, we propose an alternative method to improve the network quality. The proposed method takes as input a combined network produced by an existing network integration algorithm, and reconstructs this network to better represent the co-functionality relationships between gene pairs. At the core of the method is a learning algorithm that can learn a measure of functional similarity between genes, which we then use to reconstruct the input network. In experiments with yeast and human, the proposed method produced improved networks and achieved more accurate results than two other leading gene function prediction approaches.

CONCLUSIONS

The results show that it is possible to improve the accuracy of network-based gene function prediction methods by optimizing combined networks with appropriate similarity measures learned from data. The proposed learning procedure can handle noisy training data and scales well to large genomes.

Prediction and validation of gene-disease associations using methods inspired by social network analyses

Correctly identifying associations of genes with diseases has long been a goal in biology. With the emergence of large-scale gene-phenotype association datasets in biology, we can leverage statistical and machine learning methods to help us achieve this goal. In this paper, we present two methods for predicting gene-disease associations based on functional gene associations and gene-phenotype associations in model organisms. The first method, the Katz measure, is motivated from its success in social network link prediction, and is very closely related to some of the recent methods proposed for gene-disease association inference. The second method, called Catapult (Combining dATa Across species using Positive-Unlabeled Learning Techniques), is a supervised machine learning method that uses a biased support vector machine where the features are derived from walks in a heterogeneous gene-trait network. We study the performance of the proposed methods and related state-of-the-art methods using two different evaluation strategies, on two distinct data sets, namely OMIM phenotypes and drug-target interactions. Finally, by measuring the performance of the methods using two different evaluation strategies, we show that even though both methods perform very well, the Katz measure is better at identifying associations between traits and poorly studied genes, whereas Catapult is better suited to correctly identifying gene-trait associations overall.

The authors want to thank Jon Laurent and Kris McGary for some of the data used, and Li and Patra for making their code available. Most of Ambuj Tewari's contribution to this work happened while he was a postdoctoral fellow at the University of Texas at Austin.

RIDDLE: reflective diffusion and local extension reveal functional associations for unannotated gene sets via proximity in a gene network

The growing availability of large-scale functional networks has promoted the development of many successful techniques for predicting functions of genes. Here we extend these network-based principles and techniques to functionally characterize whole sets of genes. We present RIDDLE (Reflective Diffusion and Local Extension), which uses well developed guilt-by-association principles upon a human gene network to identify associations of gene sets. RIDDLE is particularly adept at characterizing sets with no annotations, a major challenge where most traditional set analyses fail. Notably, RIDDLE found microRNA-450a to be strongly implicated in ocular diseases and development. A web application is available at http://www.functionalnet.org/RIDDLE.

CvManGO, a method for leveraging computational predictions to improve literature-based Gene Ontology annotations

The set of annotations at the Saccharomyces Genome Database (SGD) that classifies the cellular function of S. cerevisiae gene products using Gene Ontology (GO) terms has become an important resource for facilitating experimental analysis. In addition to capturing and summarizing experimental results, the structured nature of GO annotations allows for functional comparison across organisms as well as propagation of functional predictions between related gene products. Due to their relevance to many areas of research, ensuring the accuracy and quality of these annotations is a priority at SGD. GO annotations are assigned either manually, by biocurators extracting experimental evidence from the scientific literature, or through automated methods that leverage computational algorithms to predict functional information. Here, we discuss the relationship between literature-based and computationally predicted GO annotations in SGD and extend a strategy whereby comparison of these two types of annotation identifies genes whose annotations need review. Our method, CvManGO (Computational versus Manual GO annotations), pairs literature-based GO annotations with computational GO predictions and evaluates the relationship of the two terms within GO, looking for instances of discrepancy. We found that this method will identify genes that require annotation updates, taking an important step towards finding ways to prioritize literature review. Additionally, we explored factors that may influence the effectiveness of CvManGO in identifying relevant gene targets to find in particular those genes that are missing literature-supported annotations, but our survey found that there are no immediately identifiable criteria by which one could enrich for these under-annotated genes. Finally, we discuss possible ways to improve this strategy, and the applicability of this method to other projects that use the GO for curation.

Database URL:http://www.yeastgenome.org

Systematic identification of gene annotation errors in the widely used yeast mutation collections

The baker's yeast mutation collections are extensively used genetic resources that are the basis for many genome-wide screens and new technologies. Anecdotal evidence has previously pointed to the putative existence of a neighboring gene effect (NGE) in these collections. NGE occurs when the phenotype of a strain carrying a particular perturbed gene is due to the lack of proper function of its adjacent gene. Here we performed a large-scale study of NGEs, presenting a network-based algorithm for detecting NGEs and validating software predictions using complementation experiments. We applied our approach to four datasets uncovering a similar magnitude of NGE in each (7-15%). These results have important consequences for systems biology, as the mutation collections are extensively used in almost every aspect of the field, from genetic network analysis to functional gene annotation.

A Resource of Quantitative Functional Annotation for Homo sapiens Genes

The body of human genomic and proteomic evidence continues to grow at ever-increasing rates, while annotation efforts struggle to keep pace. A surprisingly small fraction of human genes have clear, documented associations with specific functions, and new functions continue to be found for characterized genes. Here we assembled an integrated collection of diverse genomic and proteomic data for 21,341 human genes and make quantitative associations of each to 4333 Gene Ontology terms. We combined guilt-by-profiling and guilt-by-association approaches to exploit features unique to the data types. Performance was evaluated by cross-validation, prospective validation, and by manual evaluation with the biological literature. Functional-linkage networks were also constructed, and their utility was demonstrated by identifying candidate genes related to a glioma FLN using a seed network from genome-wide association studies. Our annotations are presented—alongside existing validated annotations—in a publicly accessible and searchable web interface.

Integrative approaches for microarray data analysis

Microarrays were one of the first technologies of the genomic revolution to gain widespread adoption, rapidly expanding from a cottage industry to the source of thousands of experimental results. They were one of the first assays for which data repositories and metadata were standardized and researchers were required by many journals to make published data publicly available. Microarrays provide high-throughput insights into the biological functions of genes and gene products; however, they also present a "curse of dimensionality," whereby the availability of many gene expression measurements in few samples make it challenging to distinguish noise from true biological signal. All of these factors argue for integrative approaches to microarray data analysis, which combine data from multiple experiments to increase sample size, avoid laboratory-specific bias, and enable new biological insights not possible from a single experiment. Here, we discuss several approaches to integrative microarray analysis for a diverse range of applications, including biomarker discovery, gene function and interaction prediction, and regulatory network inference. We also show how, by integrating large microarray compendia with diverse genomic data types, more nuanced biological hypotheses can be explored computationally. This chapter provides overviews and brief descriptions of each of these approaches to microarray integration.

Removing bias against membrane proteins in interaction networks

Background

Cellular interaction networks can be used to analyze the effects on cell signaling and other functional consequences of perturbations to cellular physiology. Thus, several methods have been used to reconstitute interaction networks from multiple published datasets. However, the structure and performance of these networks depends on both the quality and the unbiased nature of the original data. Due to the inherent bias against membrane proteins in protein-protein interaction (PPI) data, interaction networks can be compromised particularly if they are to be used in conjunction with drug screening efforts, since most drug-targets are membrane proteins.

Results

To overcome the experimental bias against PPIs involving membrane-associated proteins we used a probabilistic approach based on a hypergeometric distribution followed by logistic regression to simultaneously optimize the weights of different sources of interaction data. The resulting less biased genome-scale network constructed for the budding yeast Saccharomyces cerevisiae revealed that the starvation pathway is a distinct subnetwork of autophagy and retrieved a more integrated network of unfolded protein response genes. We also observed that the centrality-lethality rule depends on the content of membrane proteins in networks.

Conclusions

We show here that the bias against membrane proteins can and should be corrected in order to have a better representation of the interactions and topological properties of protein interaction networks.

Shift in chicken intestinal gene association networks after infection with Salmonella

A primary infection of Salmonella enteritidis causes a spatial-temporal dependent change in the gene expression patterns in the intestine of chickens (Gallus gallus). This is the result of a dynamic intestinal response to adapt to the altered environment and to optimize its 'health' and functionality under the new circumstances. By inferring gene association networks (GANs), the complexities of and changes in biological networks can be uncovered. Within such GANs highly interacting (hub) genes can be identified, which are supposed to be high-level regulators connected to multiple processes. By exploring the intestinal expression of genes differing between control and Salmonella infected chicken in a time-dependent manner differences in GANs were found. In control chickens more developmental processes were observed, whereas in infected chickens relatively more processes were associated to 'defense/pathogen response'. Moreover the conserved protein domains of the identified hub genes in controls were nuclear-associated, whereas hub genes in infected chickens were involved in 'cellular communication'. The shift in topology and functionality of the intestinal GANs in control and Salmonella infected animals and the identification of GAN-specific hubs is a first step to understand the complexity of biological networks and processes regulating intestinal health and functionality under normal and disturbed conditions.

A protein network-guided screen for cell cycle regulators in Drosophila

Background

Large-scale RNAi-based screens are playing a critical role in defining sets of genes that regulate specific cellular processes. Numerous screens have been completed and in some cases more than one screen has examined the same cellular process, enabling a direct comparison of the genes identified in separate screens. Surprisingly, the overlap observed between the results of similar screens is low, suggesting that RNAi screens have relatively high levels of false positives, false negatives, or both.

Results

We re-examined genes that were identified in two previous RNAi-based cell cycle screens to identify potential false positives and false negatives. We were able to confirm many of the originally observed phenotypes and to reveal many likely false positives. To identify potential false negatives from the previous screens, we used protein interaction networks to select genes for re-screening. We demonstrate cell cycle phenotypes for a significant number of these genes and show that the protein interaction network is an efficient predictor of new cell cycle regulators. Combining our results with the results of the previous screens identified a group of validated, high-confidence cell cycle/cell survival regulators. Examination of the subset of genes from this group that regulate the G1/S cell cycle transition revealed the presence of multiple members of three structurally related protein complexes: the eukaryotic translation initiation factor 3 (eIF3) complex, the COP9 signalosome, and the proteasome lid. Using a combinatorial RNAi approach, we show that while all three of these complexes are required for Cdk2/Cyclin E activity, the eIF3 complex is specifically required for some other step that limits the G1/S cell cycle transition.

Conclusions

Our results show that false positives and false negatives each play a significant role in the lack of overlap that is observed between similar large-scale RNAi-based screens. Our results also show that protein network data can be used to minimize false negatives and false positives and to more efficiently identify comprehensive sets of regulators for a process. Finally, our data provides a high confidence set of genes that are likely to play key roles in regulating the cell cycle or cell survival.

True path rule hierarchical ensembles for genome-wide gene function prediction

Gene function prediction is a complex computational problem, characterized by several items: the number of functional classes is large, and a gene may belong to multiple classes; functional classes are structured according to a hierarchy; classes are usually unbalanced, with more negative than positive examples; class labels can be uncertain and the annotations largely incomplete; to improve the predictions, multiple sources of data need to be properly integrated. In this contribution, we focus on the first three items, and, in particular, on the development of a new method for the hierarchical genome-wide and ontology-wide gene function prediction. The proposed algorithm is inspired by the “true path rule” (TPR) that governs both the Gene Ontology and FunCat taxonomies. According to this rule, the proposed TPR ensemble method is characterized by a two-way asymmetric flow of information that traverses the graph-structured ensemble: positive predictions for a node influence in a recursive way its ancestors, while negative predictions influence its offsprings. Cross-validated results with the model organism S. Crevisiae, using seven different sources of biomolecular data, and a theoretical analysis of the the TPR algorithm show the effectiveness and the drawbacks of the proposed approach.

Using computational predictions to improve literature-based Gene Ontology annotations: a feasibility study

Annotation using Gene Ontology (GO) terms is one of the most important ways in which biological information about specific gene products can be expressed in a searchable, computable form that may be compared across genomes and organisms. Because literature-based GO annotations are often used to propagate functional predictions between related proteins, their accuracy is critically important. We present a strategy that employs a comparison of literature-based annotations with computational predictions to identify and prioritize genes whose annotations need review. Using this method, we show that comparison of manually assigned ‘unknown’ annotations in the Saccharomyces Genome Database (SGD) with InterPro-based predictions can identify annotations that need to be updated. A survey of literature-based annotations and computational predictions made by the Gene Ontology Annotation (GOA) project at the European Bioinformatics Institute (EBI) across several other databases shows that this comparison strategy could be used to maintain and improve the quality of GO annotations for other organisms besides yeast. The survey also shows that although GOA-assigned predictions are the most comprehensive source of functional information for many genomes, a large proportion of genes in a variety of different organisms entirely lack these predictions but do have manual annotations. This underscores the critical need for manually performed, literature-based curation to provide functional information about genes that are outside the scope of widely used computational methods. Thus, the combination of manual and computational methods is essential to provide the most accurate and complete functional annotation of a genome.

Database URL:http://www.yeastgenome.org

It's the machine that matters: Predicting gene function and phenotype from protein networks

Increasing knowledge about the organization of proteins into complexes, systems, and pathways has led to a flowering of theoretical approaches for exploiting this knowledge in order to better learn the functions of proteins and their roles underlying phenotypic traits and diseases. Much of this body of theory has been developed and tested in model organisms, relying on their relative simplicity and genetic and biochemical tractability to accelerate the research. In this review, we discuss several of the major approaches for computationally integrating proteomics and genomics observations into integrated protein networks, then applying guilt-by-association in these networks in order to identify genes underlying traits. Recent trends in this field include a rising appreciation of the modular network organization of proteins underlying traits or mutational phenotypes, and how to exploit such protein modularity using computational approaches related to the internet search algorithm PageRank. Many protein network-based predictions have recently been experimentally confirmed in yeast, worms, plants, and mice, and several successful approaches in model organisms have been directly translated to analyze human disease, with notable recent applications to glioma and breast cancer prognosis.

Transcriptional-metabolic networks in beta-carotene-enriched potato tubers: the long and winding road to the Golden phenotype

Vitamin A deficiency is a public health problem in a large number of countries. Biofortification of major staple crops (wheat [Triticum aestivum], rice [Oryza sativa], maize [Zea mays], and potato [Solanum tuberosum]) with β-carotene has the potential to alleviate this nutritional problem. Previously, we engineered transgenic "Golden" potato tubers overexpressing three bacterial genes for β-carotene synthesis (CrtB, CrtI, and CrtY, encoding phytoene synthase, phytoene desaturase, and lycopene β-cyclase, respectively) and accumulating the highest amount of β-carotene in the four aforementioned crops. Here, we report the systematic quantitation of carotenoid metabolites and transcripts in 24 lines carrying six different transgene combinations under the control of the 35S and Patatin (Pat) promoters. Low levels of B-I expression are sufficient for interfering with leaf carotenogenesis, but not for β-carotene accumulation in tubers and calli, which requires high expression levels of all three genes under the control of the Pat promoter. Tubers expressing the B-I transgenes show large perturbations in the transcription of endogenous carotenoid genes, with only minor changes in carotenoid content, while the opposite phenotype (low levels of transcriptional perturbation and high carotenoid levels) is observed in Golden (Y-B-I) tubers. We used hierarchical clustering and pairwise correlation analysis, together with a new method for network correlation analysis, developed for this purpose, to assess the perturbations in transcript and metabolite levels in transgenic leaves and tubers. Through a "guilt-by-profiling" approach, we identified several endogenous genes for carotenoid biosynthesis likely to play a key regulatory role in Golden tubers, which are candidates for manipulations aimed at the further optimization of tuber carotenoid content.