GUILDify: a web server for phenotypic characterization of genes through biological data integration and network-based prioritization algorithms

Investigating the genomic background of calving-related traits in Canadian Jersey cattle

Traits related to calving have a significant impact on animal welfare and farm profitability in dairy production systems. Identifying genomic regions associated with calving traits could contribute to refining dairy cattle breeding programs and management practices in the dairy industry. Therefore, the primary objectives of this study were to estimate genetic parameters and perform genome-wide association studies (GWAS) and functional enrichment analyses for stillbirth, gestation length, calf size, and calving ease traits in North American Jersey cattle. A total of 40,503 animals with phenotypic records and 5,398 animals genotyped for 45,101 single nucleotide polymorphisms (SNPs) were included in the analyses. Genetic parameters were estimated based on animal models and Bayesian methods. The effects of SNPs were estimated using the Single-step Genomic Best Linear Unbiased Prediction (ssGBLUP) method. The heritability (standard error) estimates ranged from 0.01 (0.01) for stillbirths (SB) in heifers to 0.11 (0.01) for gestation length (GL) in cows. The genetic correlations ranged from -0.58 (0.11) between calving ease (CE) and SB in heifers to 0.44 (0.14) between calving ease and calf size (CZ) in cows. CE showed the highest genetic correlation between heifers and cows, 0.8 (0.22) respectively. The candidate genes identified, including MTHFR, SERPINA5, IGFBP3, and ZRANB1, are involved in key biological processes and metabolic pathways related to the studied traits. Reducing environmental variation and identifying novel indicators of reproduction traits in the Jersey breed are needed given the low heritability estimates for most traits evaluated in this study. In conclusion, this study provides a characterization of the genetic background of calving-related traits in Jersey cattle. The estimates obtained can be used to improve or build selection indexes in Jersey cattle breeding programs in North America.

Genetic Foundations of Nellore Traits: A Gene Prioritization and Functional Analyses of Genome-Wide Association Study Results

The main goal of this study was to pinpoint functional candidate genes associated with multiple economically important traits in Nellore cattle. After quality control, 1830 genomic regions sourced from 52 scientific peer-reviewed publications were used in this study. From these, a total of 8569 positional candidate genes were annotated for reproduction, 11,195 for carcass, 5239 for growth, and 3483 for morphological traits, and used in an over-representation analysis. The significant genes (adjusted p-values < 0.05) identified in the over-representation analysis underwent prioritization analyses, and enrichment analysis of the prioritized over-represented candidate genes was performed. The prioritized candidate genes were GFRA4, RFWD3, SERTAD2, KIZ, REM2, and ANKRD34B for reproduction; RFWD3, TMEM120A, MIEF2, FOXRED2, DUSP29, CARHSP1, OBI1, JOSD1, NOP58, and LOXL1-AS1 for the carcass; ANKRD34B and JOSD1 for growth traits; and no genes were prioritized for morphological traits. The functional analysis pinpointed the following genes: KIZ (plays a crucial role in spindle organization, which is essential in forming a robust mitotic centrosome), DUSP29 (involved in muscle cell differentiation), and JOSD1 (involved in protein deubiquitination, thereby improving growth). The enrichment of the functional candidate genes identified in this study highlights that these genes play an important role in the expression of reproduction, carcass, and growth traits in Nellore cattle.

Genome-wide association study considering genotype-by-environment interaction for productive and reproductive traits using whole-genome sequencing in Nellore cattle

BACKGROUND

The genotype-by-environment interaction (GxE) in beef cattle can be investigated using reaction norm models to assess environmental sensitivity and, combined with genome-wide association studies (GWAS), to map genomic regions related to animal adaptation. Including genetic markers from whole-genome sequencing in reaction norm (RN) models allows us to identify high-resolution candidate genes across environmental gradients through GWAS. Hence, we performed a GWAS via the RN approach using whole-genome sequencing data, focusing on mapping candidate genes associated with the expression of reproductive and growth traits in Nellore cattle. For this purpose, we used phenotypic data for age at first calving (AFC), scrotal circumference (SC), post-weaning weight gain (PWG), and yearling weight (YW). A total of 20,000 males and 7,159 females genotyped with 770k were imputed to the whole sequence (29 M). After quality control and linkage disequilibrium (LD) pruning, there remained ∼ 2.41 M SNPs for SC, PWG, and YW and ∼ 5.06 M SNPs for AFC.

RESULTS

Significant SNPs were identified on Bos taurus autosomes (BTA) 10, 11, 14, 18, 19, 20, 21, 24, 25 and 27 for AFC and on BTA 4, 5 and 8 for SC. For growth traits, significant SNP markers were identified on BTA 3, 5 and 20 for YW and PWG. A total of 56 positional candidate genes were identified for AFC, 9 for SC, 3 for PWG, and 24 for YW. The significant SNPs detected for the reaction norm coefficients in Nellore cattle were found to be associated with growth, adaptative, and reproductive traits. These candidate genes are involved in biological mechanisms related to lipid metabolism, immune response, mitogen-activated protein kinase (MAPK) signaling pathway, and energy and phosphate metabolism.

CONCLUSIONS

GWAS results highlighted differences in the physiological processes linked to lipid metabolism, immune response, MAPK signaling pathway, and energy and phosphate metabolism, providing insights into how different environmental conditions interact with specific genes affecting animal adaptation, productivity, and reproductive performance. The shared genomic regions between the intercept and slope are directly implicated in the regulation of growth and reproductive traits in Nellore cattle raised under different environmental conditions.

Identification of candidate regulatory genes for intramuscular fatty acid composition in pigs by transcriptome analysis

BACKGROUND

Intramuscular fat (IMF) content and its fatty acid (FA) composition are typically controlled by several genes, each with a small effect. In the current study, to pinpoint candidate genes and putative regulators involved in FA composition, we performed a multivariate integrative analysis between intramuscular FA and transcriptome profiles of porcine longissimus dorsi (LD) muscle. We also carried out a combination of network, regulatory impact factor (RIF), in silico prediction of putative target genes, and functional analyses to better support the biological relevance of our findings.

RESULTS

For this purpose, we used LD RNA-Seq and intramuscular FA composition profiles of 129 Iberian × Duroc backcrossed pigs. We identified 378 correlated variables (13 FA and 365 genes), including six FA (C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, C18:0, and C16:1n-7) that were among the most interconnected variables in the predicted network. The detected FA-correlated genes include genes involved in lipid and/or carbohydrate metabolism or in regulation of IMF deposition (e.g., ADIPOQ, CHUK, CYCS, CYP4B1, DLD, ELOVL6, FBP1, G0S2, GCLC, HMGCR, IDH3A, LEP, LGALS12, LPIN1, PLIN1, PNPLA8, PPP1R1B, SDR16C5, SFRP5, SOD3, SNW1, and TFRC), meat quality (GALNT15, GOT1, MDH1, NEU3, PDHA1, SDHD, and UNC93A), and transport (e.g., EXOC7 and SLC44A2). Functional analysis highlighted 54 over-represented gene ontology terms, including well-known biological processes and pathways that regulate lipid and carbohydrate metabolism. RIF analysis suggested a pivotal role for six transcription factors (CARHSP1, LBX1, MAFA, PAX7, SIX5, and TADA2A) as putative regulators of gene expression and intramuscular FA composition. Based on in silico prediction, we identified putative target genes for these six regulators. Among these, TADA2A and CARHSP1 had extreme RIF scores and present novel regulators in pigs. In addition, the expression of TADA2A correlated (either positively or negatively) with C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, and that of CARHSP1 correlated (positively) with the C16:1n-7 lipokine. We also found that these two transcription factors share target genes that are involved in lipid metabolism (e.g., GOT1, PLIN1, and TFRC).

CONCLUSIONS

This integrative analysis of muscle transcriptome and intramuscular FA profile revealed valuable information about key candidate genes and potential regulators for FA and lipid metabolism in pigs, among which some transcription factors are proposed to control gene expression and modulate FA composition differences.

Feet and legs malformation in Nellore cattle: genetic analysis and prioritization of GWAS results

Beef cattle affected by feet and legs malformations (FLM) cannot perform their productive and reproductive functions satisfactorily, resulting in significant economic losses. Accelerated weight gain in young animals due to increased fat deposition can lead to ligaments, tendon and joint strain and promote gene expression patterns that lead to changes in the normal architecture of the feet and legs. The possible correlated response in the FLM due to yearling weight (YW) selection suggest that this second trait could be used as an indirect selection criterion. Therefore, FLM breeding values and the genetic correlation between FLM and yearling weight (YW) were estimated for 295,031 Nellore animals by fitting a linear-threshold model in a Bayesian approach. A genome-wide association study was performed to identify genomic windows and positional candidate genes associated with FLM. The effects of single nucleotide polymorphisms (SNPs) on FLM phenotypes (affected or unaffected) were estimated using the weighted single-step genomic BLUP method, based on genotypes of 12,537 animals for 461,057 SNPs. Twelve non-overlapping windows of 20 adjacent SNPs explaining more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of candidate genes identified six genes (ATG7, EXT1, ITGA1, PPARD, SCUBE3, and SHOX) that may play a role in FLM expression due to their known role in skeletal muscle development, aberrant bone growth, lipid metabolism, intramuscular fat deposition and skeletogenesis. Identifying genes linked to foot and leg malformations enables selective breeding for healthier herds by reducing the occurrence of these conditions. Genetic markers can be used to develop tests that identify carriers of these mutations, assisting breeders in making informed breeding decisions to minimize the incidence of malformations in future generations, resulting in greater productivity and animal welfare.

Testicular hypoplasia in Nellore Cattle: Genetic analysis and functional analysis of genome-wide association study results

Characterized by the incomplete development of the germinal epithelium of the seminiferous tubules, Testicular hypoplasia (TH) leads to decreased sperm concentration, increased morphological changes in sperm and azoospermia. Economic losses resulting from the disposal of affected bulls reduce the efficiency of meat production systems. A genome-wide association study and functional analysis were performed to identify genomic windows and the underlying positional candidate genes associated with TH in Nellore cattle. Phenotypic and pedigree data from 207,195 animals and genotypes (461,057 single nucleotide polymorphism, SNP) from 17,326 sires were used in this study. TH was evaluated as a binary trait measured at 18 months of age. A possible correlated response on TH resulting from the selection for scrotal circumference was evaluated by using a two-trait analysis. Thus, estimated breeding values were calculated by fitting a linear-threshold animal model in a Bayesian approach. The SNP effects were estimated using the weighted single-step genomic BLUP method. Twelve non-overlapping windows of 20 adjacent SNP that explained more than 1% of the additive genetic variance were selected for candidate gene annotation. Functional and gene prioritization analysis of the candidate genes identified genes (KHDRBS3, GPX5, STAR, ERLIN2), which might play an important role in the expression of TH due to their known roles in the spermatogenesis process, synthesis of steroids and lipid metabolism.

Genome-wide association study using haplotype libraries and repeated measures model to identify candidate genomic regions for stillbirth in Holstein cattle

Reduced fertility is one of the main causes of economic losses on dairy farms, resulting in economic losses estimated at $938 per stillbirth case in Holstein herds. The identification of genomic regions associated with stillbirth could help to develop better management and breeding strategies aimed to reduce the frequency of undesirable gestation outcomes. Here, 10,570 cows and 50,541 birth records were used to perform a haplotype-based GWAS. A total of 41 significantly associated pseudo-SNPs (haplotypes within haplotype blocks converted to a binary classification) were identified after Bonferroni adjustment for multiple tests. A total of 117 positional candidate genes were annotated within or close (in a 200-kb interval) to significant pseudo-SNPs (haplotype blocks). The guilt-by-association functional prioritization identified 31 potential functional candidate genes for reproductive performance out of the 117 positional candidate genes annotated. These genes play crucial roles in biological processes associated with pregnancy persistence, fetus development, immune response, among others. These results helped us to better understand the genetic basis of stillbirth in dairy cattle and may be useful for the prediction of stillbirth in Holstein cattle, helping to reduce the related economic losses caused by this phenotype.

Determination of shared genetic etiology and possible causal relations between tobacco smoking and depression

BACKGROUNDS

Cigarette smoking is strongly associated with major depressive disorder (MDD). However, any genetic etiology of such comorbidity and causal relations is poorly understood, especially at the genome-wide level.

METHODS

In the present in silico research, we analyzed summary data from the genome-wide association study of the Psychiatric Genetic Consortium for MDD (n = 191 005) and UK Biobank for smoking (n = 337 030) by using various biostatistical methods including Bayesian colocalization analysis, LD score regression, variant effect size correlation analysis, and Mendelian randomization (MR).

RESULTS

By adopting a gene prioritization approach, we identified 43 genes shared by MDD and smoking, which were significantly enriched in membrane potential, gamma-aminobutyric acid receptor activity, and retrograde endocannabinoid signaling pathways, indicating that the comorbid mechanisms are involved in the neurotransmitter system. According to linkage disequilibrium score regression, we found a strong positive correlation between MDD and current smoking (rg = 0.365; p = 7.23 × 10-25) and a negative correlation between MDD and former smoking (rg = -0.298; p = 1.59 × 10-24). MR analysis suggested that genetic liability for depression increased smoking.

CONCLUSIONS

These findings inform the concomitant conditions of MDD and smoking and support the use of self-medication with smoking to counteract depression.

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.

Genome-wide association study to identify genomic regions and positional candidate genes associated with male fertility in beef cattle

Fertility plays a key role in the success of calf production, but there is evidence that reproductive efficiency in beef cattle has decreased during the past half-century worldwide. Therefore, identifying animals with superior fertility could significantly impact cow-calf production efficiency. The objective of this research was to identify candidate regions affecting bull fertility in beef cattle and positional candidate genes annotated within these regions. A GWAS using a weighted single-step genomic BLUP approach was performed on 265 crossbred beef bulls to identify markers associated with scrotal circumference (SC) and sperm motility (SM). Eight windows containing 32 positional candidate genes and five windows containing 28 positional candidate genes explained more than 1% of the genetic variance for SC and SM, respectively. These windows were selected to perform gene annotation, QTL enrichment, and functional analyses. Functional candidate gene prioritization analysis revealed 14 prioritized candidate genes for SC of which MAP3K1 and VIP were previously found to play roles in male fertility. A different set of 14 prioritized genes were identified for SM and five were previously identified as regulators of male fertility (SOD2, TCP1, PACRG, SPEF2, PRLR). Significant enrichment results were identified for fertility and body conformation QTLs within the candidate windows. Gene ontology enrichment analysis including biological processes, molecular functions, and cellular components revealed significant GO terms associated with male fertility. The identification of these regions contributes to a better understanding of fertility associated traits and facilitates the discovery of positional candidate genes for future investigation of causal mutations and their implications.

Galaxy InteractoMIX: An Integrated Computational Platform for the Study of Protein-Protein Interaction Data

Protein interactions play a crucial role among the different functions of a cell and are central to our understanding of cellular processes both in health and disease. Here we present Galaxy InteractoMIX (http://galaxy.interactomix.com), a platform composed of 13 different computational tools each addressing specific aspects of the study of protein-protein interactions, ranging from large-scale cross-species protein-wide interactomes to atomic resolution level of protein complexes. Galaxy InteractoMIX provides an intuitive interface where users can retrieve consolidated interactomics data distributed across several databases or uncover links between diseases and genes by analyzing the interactomes underlying these diseases. The platform makes possible large-scale prediction and curation protein interactions using the conservation of motifs, interology, or presence or absence of key sequence signatures. The range of structure-based tools includes modeling and analysis of protein complexes, delineation of interfaces and the modeling of peptides acting as inhibitors of protein-protein interactions. Galaxy InteractoMIX includes a range of ready-to-use workflows to run complex analyses requiring minimal intervention by users. The potential range of applications of the platform covers different aspects of life science, biomedicine, biotechnology and drug discovery where protein associations are studied.

Weighted Gene Correlation Network Meta-Analysis Reveals Functional Candidate Genes Associated with High- and Sub-Fertile Reproductive Performance in Beef Cattle

Improved reproductive efficiency could lead to economic benefits for the beef industry, once the intensive selection pressure has led to a decreased fertility. However, several factors limit our understanding of fertility traits, including genetic differences between populations and statistical limitations. In the present study, the RNA-sequencing data from uterine samples of high-fertile (HF) and sub-fertile (SF) animals was integrated using co-expression network meta-analysis, weighted gene correlation network analysis, identification of upstream regulators, variant calling, and network topology approaches. Using this pipeline, top hub-genes harboring fixed variants (HF × SF) were identified in differentially co-expressed gene modules (DcoExp). The functional prioritization analysis identified the genes with highest potential to be key-regulators of the DcoExp modules between HF and SF animals. Consequently, 32 functional candidate genes (10 upstream regulators and 22 top hub-genes of DcoExp modules) were identified. These genes were associated with the regulation of relevant biological processes for fertility, such as embryonic development, germ cell proliferation, and ovarian hormone regulation. Additionally, 100 candidate variants (single nucleotide polymorphisms (SNPs) and insertions and deletions (INDELs)) were identified within those genes. In the long-term, the results obtained here may help to reduce the frequency of subfertility in beef herds, reducing the associated economic losses caused by this condition.

Gene Prioritization through Consensus Strategy, Enrichment Methodologies Analysis, and Networking for Osteosarcoma Pathogenesis

Osteosarcoma is the most common subtype of primary bone cancer, affecting mostly adolescents. In recent years, several studies have focused on elucidating the molecular mechanisms of this sarcoma; however, its molecular etiology has still not been determined with precision. Therefore, we applied a consensus strategy with the use of several bioinformatics tools to prioritize genes involved in its pathogenesis. Subsequently, we assessed the physical interactions of the previously selected genes and applied a communality analysis to this protein-protein interaction network. The consensus strategy prioritized a total list of 553 genes. Our enrichment analysis validates several studies that describe the signaling pathways PI3K/AKT and MAPK/ERK as pathogenic. The gene ontology described TP53 as a principal signal transducer that chiefly mediates processes associated with cell cycle and DNA damage response It is interesting to note that the communality analysis clusters several members involved in metastasis events, such as MMP2 and MMP9, and genes associated with DNA repair complexes, like ATM, ATR, CHEK1, and RAD51. In this study, we have identified well-known pathogenic genes for osteosarcoma and prioritized genes that need to be further explored.

Embedding Alternative Conformations of Proteins in Protein-Protein Interaction Networks

While many proteins act alone, the majority of them interact with others and form molecular complexes to undertake biological functions at both cellular and systems levels. Two proteins should have complementary shapes to physically connect to each other. As proteins are dynamic and changing their conformations, it is vital to track in which conformation a specific interaction can happen. Here, we present a step-by-step guide to embedding the protein alternative conformations in each protein-protein interaction in a systems level. All external tools/websites used in each step are explained, and some notes and suggestions are provided to clear any ambiguous point.

PathFXweb: a web application for identifying drug safety and efficacy phenotypes

Summary

Limited efficacy and intolerable safety limit therapeutic development and identification of potential liabilities earlier in development could significantly improve this process. Computational approaches which aggregate data from multiple sources and consider the drug’s pathways effects could add to identification of these liabilities earlier. Such computational methods must be accessible to a variety of users beyond computational scientists, especially regulators and industry scientists, in order to impact the therapeutic development process. We have previously developed and published PathFX, an algorithm for identifying drug networks and phenotypes for understanding drug associations to safety and efficacy. Here we present a streamlined and easy-to-use PathFX web application that allows users to search for drug networks and associated phenotypes. We have also added visualization, and phenotype clustering to improve functionality and interpretability of PathFXweb.

Availability and implementation

https://www.pathfxweb.net/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Prioritizing target-disease associations with novel safety and efficacy scoring methods

Biological target (commonly genes or proteins) identification is still largely a manual process, where experts manually try to collect and combine information from hundreds of data sources, ranging from scientific publications to omics databases. Targeting the wrong gene or protein will lead to failure of the drug development process, as well as incur delays and costs. To improve this process, different software platforms are being developed. These platforms rely strongly on efficacy estimates based on target-disease association scores created by computational methods for drug target prioritization. Here novel computational methods are presented to more accurately evaluate the efficacy and safety of potential drug targets. The proposed efficacy scores utilize existing gene expression data and tissue/disease specific networks to improve the inference of target-disease associations. Conversely, safety scores enable the identification of genes that are essential, potentially susceptible to adverse effects or carcinogenic. Benchmark results demonstrate that our transcriptome-based methods for drug target prioritization can increase the true positive rate of target-disease associations. Additionally, the proposed safety evaluation system enables accurate predictions of targets of withdrawn drugs and targets of drug trials prematurely discontinued.

Genetic mechanisms underlying spermatic and testicular traits within and among cattle breeds: systematic review and prioritization of GWAS results

Reduced bull fertility imposes economic losses in bovine herds. Specifically, testicular and spermatic traits are important indicators of reproductive efficiency. Several genome-wide association studies (GWAS) have identified genomic regions associated with these fertility traits. The aims of this study were as follows: 1) to perform a systematic review of GWAS results for spermatic and testicular traits in cattle and 2) to identify key functional candidate genes for these traits. The identification of functional candidate genes was performed using a systems biology approach, where genes shared between traits and studies were evaluated by a guilt by association gene prioritization (GUILDify and ToppGene software) in order to identify the best functional candidates. These candidate genes were integrated and analyzed in order to identify overlapping patterns among traits and breeds. Results showed that GWAS for testicular-related traits have been developed for beef breeds only, whereas the majority of GWAS for spermatic-related traits were conducted using dairy breeds. When comparing traits measured within the same study, the highest number of genes shared between different traits was observed, indicating a high impact of the population genetic structure and environmental effects. Several chromosomal regions were enriched for functional candidate genes associated with fertility traits. Moreover, multiple functional candidate genes were enriched for markers in a species-specific basis, taurine (Bos taurus) or indicine (Bos indicus). For the different candidate regions identified in the GWAS in the literature, functional candidate genes were detected as follows: B. Taurus chromosome X (BTX) (TEX11, IRAK, CDK16, ATP7A, ATRX, HDAC6, FMR1, L1CAM, MECP2, etc.), BTA17 (TRPV4 and DYNLL1), and BTA14 (MOS, FABP5, ZFPM2). These genes are responsible for regulating important metabolic pathways or biological processes associated with fertility, such as progression of spermatogenesis, control of ciliary activity, development of Sertoli cells, DNA integrity in spermatozoa, and homeostasis of testicular cells. This study represents the first systematic review on male fertility traits in cattle using a system biology approach to identify key candidate genes for these traits.

Gene prioritization, communality analysis, networking and metabolic integrated pathway to better understand breast cancer pathogenesis

Consensus strategy was proved to be highly efficient in the recognition of gene-disease association. Therefore, the main objective of this study was to apply theoretical approaches to explore genes and communities directly involved in breast cancer (BC) pathogenesis. We evaluated the consensus between 8 prioritization strategies for the early recognition of pathogenic genes. A communality analysis in the protein-protein interaction (PPi) network of previously selected genes was enriched with gene ontology, metabolic pathways, as well as oncogenomics validation with the OncoPPi and DRIVE projects. The consensus genes were rationally filtered to 1842 genes. The communality analysis showed an enrichment of 14 communities specially connected with ERBB, PI3K-AKT, mTOR, FOXO, p53, HIF-1, VEGF, MAPK and prolactin signaling pathways. Genes with highest ranking were TP53, ESR1, BRCA2, BRCA1 and ERBB2. Genes with highest connectivity degree were TP53, AKT1, SRC, CREBBP and EP300. The connectivity degree allowed to establish a significant correlation between the OncoPPi network and our BC integrated network conformed by 51 genes and 62 PPi. In addition, CCND1, RAD51, CDC42, YAP1 and RPA1 were functional genes with significant sensitivity score in BC cell lines. In conclusion, the consensus strategy identifies both well-known pathogenic genes and prioritized genes that need to be further explored.

Context-enriched interactome powered by proteomics helps the identification of novel regulators of macrophage activation

The role of pro-inflammatory macrophage activation in cardiovascular disease (CVD) is a complex one amenable to network approaches. While an indispensible tool for elucidating the molecular underpinnings of complex diseases including CVD, the interactome is limited in its utility as it is not specific to any cell type, experimental condition or disease state. We introduced context-specificity to the interactome by combining it with co-abundance networks derived from unbiased proteomics measurements from activated macrophage-like cells. Each macrophage phenotype contributed to certain regions of the interactome. Using a network proximity-based prioritization method on the combined network, we predicted potential regulators of macrophage activation. Prediction performance significantly increased with the addition of co-abundance edges, and the prioritized candidates captured inflammation, immunity and CVD signatures. Integrating the novel network topology with transcriptomics and proteomics revealed top candidate drivers of inflammation. In vitro loss-of-function experiments demonstrated the regulatory role of these proteins in pro-inflammatory signaling.

When human cells or tissues are injured, the body triggers a response known as inflammation to repair the damage and protect itself from further harm. However, if the same issue keeps recurring, the tissues become inflamed for longer periods of time, which may ultimately lead to health problems. This is what could be happening in cardiovascular diseases, where long-term inflammation could damage the heart and blood vessels.

Many different proteins interact with each other to control inflammation; gaining an insight into the nature of these interactions could help to pinpoint the role of each molecular actor. Researchers have used a combination of unbiased, large-scale experimental and computational approaches to develop the interactome, a map of the known interactions between all proteins in humans. However, interactions between proteins can change between cell types, or during disease. Here, Halu et al. aimed to refine the human interactome and identify new proteins involved in inflammation, especially in the context of cardiovascular disease.

Cells called macrophages produce signals that trigger inflammation whey they detect damage in other cells or tissues. The experiments used a technique called proteomics to measure the amounts of all the proteins in human macrophages. Combining these data with the human interactome made it possible to predict new links between proteins known to have a role in inflammation and other proteins in the interactome. Further analysis using other sets of data from macrophages helped identify two new candidate proteins – GBP1 and WARS – that may promote inflammation. Halu et al. then used a genetic approach to deactivate the genes and decrease the levels of these two proteins in macrophages, which caused the signals that encourage inflammation to drop.

These findings suggest that GBP1 and WARS regulate the activity of macrophages to promote inflammation. The two proteins could therefore be used as drug targets to treat cardiovascular diseases and other disorders linked to inflammation, but further studies will be needed to precisely dissect how GBP1 and WARS work in humans.

Predictive Systems Toxicology

In this review we address to what extent computational techniques can augment our ability to predict toxicity. The first section provides a brief history of empirical observations on toxicity dating back to the dawn of Sumerian civilization. Interestingly, the concept of dose emerged very early on, leading up to the modern emphasis on kinetic properties, which in turn encodes the insight that toxicity is not solely a property of a compound but instead depends on the interaction with the host organism. The next logical step is the current conception of evaluating drugs from a personalized medicine point of view. We review recent work on integrating what could be referred to as classical pharmacokinetic analysis with emerging systems biology approaches incorporating multiple omics data. These systems approaches employ advanced statistical analytical data processing complemented with machine learning techniques and use both pharmacokinetic and omics data. We find that such integrated approaches not only provide improved predictions of toxicity but also enable mechanistic interpretations of the molecular mechanisms underpinning toxicity and drug resistance. We conclude the chapter by discussing some of the main challenges, such as how to balance the inherent tension between the predicitive capacity of models, which in practice amounts to constraining the number of features in the models versus allowing for rich mechanistic interpretability, i.e., equipping models with numerous molecular features. This challenge also requires patient-specific predictions on toxicity, which in turn requires proper stratification of patients as regards how they respond, with or without adverse toxic effects. In summary, the transformation of the ancient concept of dose is currently successfully operationalized using rich integrative data encoded in patient-specific models.

Consensus strategy in genes prioritization and combined bioinformatics analysis for preeclampsia pathogenesis

BACKGROUND

Preeclampsia is a multifactorial disease with unknown pathogenesis. Even when recent studies explored this disease using several bioinformatics tools, the main objective was not directed to pathogenesis. Additionally, consensus prioritization was proved to be highly efficient in the recognition of genes-disease association. However, not information is available about the consensus ability to early recognize genes directly involved in pathogenesis. Therefore our aim in this study is to apply several theoretical approaches to explore preeclampsia; specifically those genes directly involved in the pathogenesis.

METHODS

We firstly evaluated the consensus between 12 prioritization strategies to early recognize pathogenic genes related to preeclampsia. A communality analysis in the protein-protein interaction network of previously selected genes was done including further enrichment analysis. The enrichment analysis includes metabolic pathways as well as gene ontology. Microarray data was also collected and used in order to confirm our results or as a strategy to weight the previously enriched pathways.

RESULTS

The consensus prioritized gene list was rationally filtered to 476 genes using several criteria. The communality analysis showed an enrichment of communities connected with VEGF-signaling pathway. This pathway is also enriched considering the microarray data. Our result point to VEGF, FLT1 and KDR as relevant pathogenic genes, as well as those connected with NO metabolism.

CONCLUSION

Our results revealed that consensus strategy improve the detection and initial enrichment of pathogenic genes, at least in preeclampsia condition. Moreover the combination of the first percent of the prioritized genes with protein-protein interaction network followed by communality analysis reduces the gene space. This approach actually identifies well known genes related with pathogenesis. However, genes like HSP90, PAK2, CD247 and others included in the first 1% of the prioritized list need to be further explored in preeclampsia pathogenesis through experimental approaches.

Enriching Traditional Protein-protein Interaction Networks with Alternative Conformations of Proteins

Traditional Protein-Protein Interaction (PPI) networks, which use a node and edge representation, lack some valuable information about the mechanistic details of biological processes. Mapping protein structures to these PPI networks not only provides structural details of each interaction but also helps us to find the mutual exclusive interactions. Yet it is not a comprehensive representation as it neglects the conformational changes of proteins which may lead to different interactions, functions, and downstream signalling. In this study, we proposed a new representation for structural PPI networks inspecting the alternative conformations of proteins. We performed a large-scale study by creating breast cancer metastasis network and equipped it with different conformers of proteins. Our results showed that although 88% of proteins in our network has at least two structures in Protein Data Bank (PDB), only 22% of them have alternative conformations and the remaining proteins have different regions saved in PDB. However, using even this small set of alternative conformations we observed a considerable increase in our protein docking predictions. Our protein-protein interaction predictions increased from 54% to 76% using the alternative conformations. We also showed the benefits of investigating structural data and alternative conformations of proteins through three case studies.

Genetic and functional characterization of disease associations explains comorbidity

Understanding relationships between diseases, such as comorbidities, has important socio-economic implications, ranging from clinical study design to health care planning. Most studies characterize disease comorbidity using shared genetic origins, ignoring pathway-based commonalities between diseases. In this study, we define the disease pathways using an interactome-based extension of known disease-genes and introduce several measures of functional overlap. The analysis reveals 206 significant links among 94 diseases, giving rise to a highly clustered disease association network. We observe that around 95% of the links in the disease network, though not identified by genetic overlap, are discovered by functional overlap. This disease network portraits rheumatoid arthritis, asthma, atherosclerosis, pulmonary diseases and Crohn's disease as hubs and thus pointing to common inflammatory processes underlying disease pathophysiology. We identify several described associations such as the inverse comorbidity relationship between Alzheimer's disease and neoplasms. Furthermore, we investigate the disruptions in protein interactions by mapping mutations onto the domains involved in the interaction, suggesting hypotheses on the causal link between diseases. Finally, we provide several proof-of-principle examples in which we model the effect of the mutation and the change of the association strength, which could explain the observed comorbidity between diseases caused by the same genetic alterations.

Combined GWAS and 'guilt by association'-based prioritization analysis identifies functional candidate genes for body size in sheep

Background

Body size in sheep is an important indicator of productivity, growth and health as well as of environmental adaptation. It is a composite quantitative trait that has been studied with high-throughput genomic methods, i.e. genome-wide association studies (GWAS) in various mammalian species. Several genomic markers have been associated with body size traits and genes have been identified as causative candidates in humans, dog and cattle. A limited number of related GWAS have been performed in various sheep breeds and have identified genomic regions and candidate genes that partly account for body size variability. Here, we conducted a GWAS in Frizarta dairy sheep with phenotypic data from 10 body size measurements and genotypic data (from Illumina ovineSNP50 BeadChip) for 459 ewes.

Results

The 10 body size measurements were subjected to principal component analysis and three independent principal components (PC) were constructed, interpretable as width, height and length dimensions, respectively. The GWAS performed for each PC identified 11 significant SNPs, at the chromosome level, one on each of the chromosomes 3, 8, 9, 10, 11, 12, 19, 20, 23 and two on chromosome 25. Nine out of the 11 SNPs were located on previously identified quantitative trait loci for sheep meat, production or reproduction. One hundred and ninety-seven positional candidate genes within a 1-Mb distance from each significant SNP were found. A guilt-by-association-based (GBA) prioritization analysis (PA) was performed to identify the most plausible functional candidate genes. GBA-based PA identified 39 genes that were significantly associated with gene networks relevant to body size traits. Prioritized genes were identified in the vicinity of all significant SNPs except for those on chromosomes 10 and 12. The top five ranking genes were TP53, BMPR1A, PIK3R5, RPL26 and PRKDC.

Conclusions

The results of this GWAS provide evidence for 39 causative candidate genes across nine chromosomal regions for body size traits, some of which are novel and some are previously identified candidates from other studies (e.g. TP53, NTN1 and ZNF521). GBA-based PA has proved to be a useful tool to identify genes with increased biological relevance but it is subjected to certain limitations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0316-3) contains supplementary material, which is available to authorized users.

InteractoMIX: a suite of computational tools to exploit interactomes in biological and clinical research

Virtually all the biological processes that occur inside or outside cells are mediated by protein–protein interactions (PPIs). Hence, the charting and description of the PPI network, initially in organisms, the interactome, but more recently in specific tissues, is essential to fully understand cellular processes both in health and disease. The study of PPIs is also at the heart of renewed efforts in the medical and biotechnological arena in the quest of new therapeutic targets and drugs. Here, we present a mini review of 11 computational tools and resources tools developed by us to address different aspects of PPIs: from interactome level to their atomic 3D structural details. We provided details on each specific resource, aims and purpose and compare with equivalent tools in the literature. All the tools are presented in a centralized, one-stop, web site: InteractoMIX (http://interactomix.com).

Network-based in silico drug efficacy screening

The increasing cost of drug development together with a significant drop in the number of new drug approvals raises the need for innovative approaches for target identification and efficacy prediction. Here, we take advantage of our increasing understanding of the network-based origins of diseases to introduce a drug-disease proximity measure that quantifies the interplay between drugs targets and diseases. By correcting for the known biases of the interactome, proximity helps us uncover the therapeutic effect of drugs, as well as to distinguish palliative from effective treatments. Our analysis of 238 drugs used in 78 diseases indicates that the therapeutic effect of drugs is localized in a small network neighborhood of the disease genes and highlights efficacy issues for drugs used in Parkinson and several inflammatory disorders. Finally, network-based proximity allows us to predict novel drug-disease associations that offer unprecedented opportunities for drug repurposing and the detection of adverse effects.

Attempts to predict novel use for existing drugs rarely consider information on the impact on the genes perturbed in a given disease. Here, the authors present a novel network-based drug-disease proximity measure that provides insight on gene specific therapeutic effect of drugs and may facilitate drug repurposing.

FN14 and GRP94 expression are prognostic/predictive biomarkers of brain metastasis outcome that open up new therapeutic strategies

Brain metastasis is a devastating problem in patients with breast, lung and melanoma tumors. GRP94 and FN14 are predictive biomarkers over-expressed in primary breast carcinomas that metastasized in brain. To further validate these brain metastasis biomarkers, we performed a multicenter study including 318 patients with breast carcinomas. Among these patients, there were 138 patients with metastasis, of whom 84 had brain metastasis. The likelihood of developing brain metastasis increased by 5.24-fold (95%CI 2.83-9.71) and 2.55- (95%CI 1.52-4.3) in the presence of FN14 and GRP94, respectively. Moreover, FN14 was more sensitive than ErbB2 (38.27 vs. 24.68) with similar specificity (89.43 vs. 89.55) to predict brain metastasis and had identical prognostic value than triple negative patients (p < 0.0001). Furthermore, we used GRP94 and FN14 pathways and GUILD, a network-based disease-gene prioritization program, to pinpoint the genes likely to be therapeutic targets, which resulted in FN14 as the main modulator and thalidomide as the best scored drug. The treatment of mice with brain metastasis improves survival decreasing reactive astrocytes and angiogenesis, and down-regulate FN14 and its ligand TWEAK. In conclusion our results indicate that FN14 and GRP94 are prediction/prognosis markers which open up new possibilities for preventing/treating brain metastasis.

HyDRA: gene prioritization via hybrid distance-score rank aggregation

Gene prioritization refers to a family of computational techniques for inferring disease genes through a set of training genes and carefully chosen similarity criteria. Test genes are scored based on their average similarity to the training set, and the rankings of genes under various similarity criteria are aggregated via statistical methods. The contributions of our work are threefold: (i) first, based on the realization that there is no unique way to define an optimal aggregate for rankings, we investigate the predictive quality of a number of new aggregation methods and known fusion techniques from machine learning and social choice theory. Within this context, we quantify the influence of the number of training genes and similarity criteria on the diagnostic quality of the aggregate and perform in-depth cross-validation studies; (ii) second, we propose a new approach to genomic data aggregation, termed HyDRA (Hybrid Distance-score Rank Aggregation), which combines the advantages of score-based and combinatorial aggregation techniques. We also propose incorporating a new top-versus-bottom (TvB) weighting feature into the hybrid schemes. The TvB feature ensures that aggregates are more reliable at the top of the list, rather than at the bottom, since only top candidates are tested experimentally; (iii) third, we propose an iterative procedure for gene discovery that operates via successful augmentation of the set of training genes by genes discovered in previous rounds, checked for consistency.

MOTIVATION

Fundamental results from social choice theory, political and computer sciences, and statistics have shown that there exists no consistent, fair and unique way to aggregate rankings. Instead, one has to decide on an aggregation approach using predefined set of desirable properties for the aggregate. The aggregation methods fall into two categories, score- and distance-based approaches, each of which has its own drawbacks and advantages. This work is motivated by the observation that merging these two techniques in a computationally efficient manner, and by incorporating additional constraints, one can ensure that the predictive quality of the resulting aggregation algorithm is very high.

RESULTS

We tested HyDRA on a number of gene sets, including autism, breast cancer, colorectal cancer, endometriosis, ischaemic stroke, leukemia, lymphoma and osteoarthritis. Furthermore, we performed iterative gene discovery for glioblastoma, meningioma and breast cancer, using a sequentially augmented list of training genes related to the Turcot syndrome, Li-Fraumeni condition and other diseases. The methods outperform state-of-the-art software tools such as ToppGene and Endeavour. Despite this finding, we recommend as best practice to take the union of top-ranked items produced by different methods for the final aggregated list.

AVAILABILITY AND IMPLEMENTATION

The HyDRA software may be downloaded from: http://web.engr.illinois.edu/∼mkim158/HyDRA.zip.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.