Knowledge Graph Completion to Predict Polypharmacy Side Effects

The polypharmacy side effect prediction problem considers cases in which two drugs taken individually do not result in a particular side effect; however, when the two drugs are taken in combination, the side effect manifests. In this work, we demonstrate that multi-relational knowledge graph completion achieves state-of-the-art results on the polypharmacy side effect prediction problem. Empirical results show that our approach is particularly effective when the protein targets of the drugs are well-characterized. In contrast to prior work, our approach provides more interpretable predictions and hypotheses for wet lab validation.Comment: 13th International Conference on Data Integration in the Life Sciences (DILS2018

Gene expression drives the evolution of dominance.

Dominance is a fundamental concept in molecular genetics and has implications for understanding patterns of genetic variation, evolution, and complex traits. However, despite its importance, the degree of dominance in natural populations is poorly quantified. Here, we leverage multiple mating systems in natural populations of Arabidopsis to co-estimate the distribution of fitness effects and dominance coefficients of new amino acid changing mutations. We find that more deleterious mutations are more likely to be recessive than less deleterious mutations. Further, this pattern holds across gene categories, but varies with the connectivity and expression patterns of genes. Our work argues that dominance arises as a consequence of the functional importance of genes and their optimal expression levels

Multitask Hopfield Networks

Multitask algorithms typically use task similarity information as a bias to speed up and improve the performance of learning processes. Tasks are learned jointly, sharing information across them, in order to construct models more accurate than those learned separately over single tasks. In this contribution, we present the first multitask model, to our knowledge, based on Hopfield Networks (HNs), named HoMTask. We show that by appropriately building a unique HN embedding all tasks, a more robust and effective classification model can be learned. HoMTask is a transductive semi-supervised parametric HN, that minimizes an energy function extended to all nodes and to all tasks under study. We provide theoretical evidence that the optimal parameters automatically estimated by HoMTask make coherent the model itself with the prior knowledge (connection weights and node labels). The convergence properties of HNs are preserved, and the fixed point reached by the network dynamics gives rise to the prediction of unlabeled nodes. The proposed model improves the classification abilities of singletask HNs on a preliminary benchmark comparison, and achieves competitive performance with state-of-the-art semi-supervised graph-based algorithms.Comment: 16 pages, 1 figure

Functional Diversity and Structural Disorder in the Human Ubiquitination Pathway

The ubiquitin-proteasome system plays a central role in cellular regulation and protein quality control (PQC). The system is built as a pyramid of increasing complexity, with two E1 (ubiquitin activating), few dozen E2 (ubiquitin conjugating) and several hundred E3 (ubiquitin ligase) enzymes. By collecting and analyzing E3 sequences from the KEGG BRITE database and literature, we assembled a coherent dataset of 563 human E3s and analyzed their various physical features. We found an increase in structural disorder of the system with multiple disorder predictors (IUPred - E1: 5.97%, E2: 17.74%, E3: 20.03%). E3s that can bind E2 and substrate simultaneously (single subunit E3, ssE3) have significantly higher disorder (22.98%) than E3s in which E2 binding (multi RING-finger, mRF, 0.62%), scaffolding (6.01%) and substrate binding (adaptor/substrate recognition subunits, 17.33%) functions are separated. In ssE3s, the disorder was localized in the substrate/adaptor binding domains, whereas the E2-binding RING/HECT-domains were structured. To demonstrate the involvement of disorder in E3 function, we applied normal modes and molecular dynamics analyses to show how a disordered and highly flexible linker in human CBL (an E3 that acts as a regulator of several tyrosine kinase-mediated signalling pathways) facilitates long-range conformational changes bringing substrate and E2-binding domains towards each other and thus assisting in ubiquitin transfer. E3s with multiple interaction partners (as evidenced by data in STRING) also possess elevated levels of disorder (hubs, 22.90% vs. non-hubs, 18.36%). Furthermore, a search in PDB uncovered 21 distinct human E3 interactions, in 7 of which the disordered region of E3s undergoes induced folding (or mutual induced folding) in the presence of the partner. In conclusion, our data highlights the primary role of structural disorder in the functions of E3 ligases that manifests itself in the substrate/adaptor binding functions as well as the mechanism of ubiquitin transfer by long-range conformational transitions. © 2013 Bhowmick et al

Spina bifida-predisposing heterozygous mutations in Planar Cell Polarity genes and Zic2 reduce bone mass in young mice

Fractures are a common comorbidity in children with the neural tube defect (NTD) spina bifida. Mutations in the Wnt/planar cell polarity (PCP) pathway contribute to NTDs in humans and mice, but whether this pathway independently determines bone mass is poorly understood. Here, we first confirmed that core Wnt/PCP components are expressed in osteoblasts and osteoclasts in vitro. In vivo, we performed detailed µCT comparisons of bone structure in tibiae from young male mice heterozygous for NTD-associated mutations versus WT littermates. PCP signalling disruption caused by Vangl2 (Vangl2Lp/+) or Celsr1 (Celsr1Crsh/+) mutations significantly reduced trabecular bone mass and distal tibial cortical thickness. NTD-associated mutations in non-PCP transcription factors were also investigated. Pax3 mutation (Pax3Sp2H/+) had minimal effects on bone mass. Zic2 mutation (Zic2Ku/+) significantly altered the position of the tibia/fibula junction and diminished cortical bone in the proximal tibia. Beyond these genes, we bioinformatically documented the known extent of shared genetic networks between NTDs and bone properties. 46 genes involved in neural tube closure are annotated with bone-related ontologies. These findings document shared genetic networks between spina bifida risk and bone structure, including PCP components and Zic2. Genetic variants which predispose to spina bifida may therefore independently diminish bone mass

Multi-level evidence of an allelic hierarchy of USH2A variants in hearing, auditory processing and speech/language outcomes.

Language development builds upon a complex network of interacting subservient systems. It therefore follows that variations in, and subclinical disruptions of, these systems may have secondary effects on emergent language. In this paper, we consider the relationship between genetic variants, hearing, auditory processing and language development. We employ whole genome sequencing in a discovery family to target association and gene x environment interaction analyses in two large population cohorts; the Avon Longitudinal Study of Parents and Children (ALSPAC) and UK10K. These investigations indicate that USH2A variants are associated with altered low-frequency sound perception which, in turn, increases the risk of developmental language disorder. We further show that Ush2a heterozygote mice have low-level hearing impairments, persistent higher-order acoustic processing deficits and altered vocalizations. These findings provide new insights into the complexity of genetic mechanisms serving language development and disorders and the relationships between developmental auditory and neural systems

Timed inhibition of CDC7 increases CRISPR-Cas9 mediated templated repair.

Repair of double strand DNA breaks (DSBs) can result in gene disruption or gene modification via homology directed repair (HDR) from donor DNA. Altering cellular responses to DSBs may rebalance editing outcomes towards HDR and away from other repair outcomes. Here, we utilize a pooled CRISPR screen to define host cell involvement in HDR between a Cas9 DSB and a plasmid double stranded donor DNA (dsDonor). We find that the Fanconi Anemia (FA) pathway is required for dsDonor HDR and that other genes act to repress HDR. Small molecule inhibition of one of these repressors, CDC7, by XL413 and other inhibitors increases the efficiency of HDR by up to 3.5 fold in many contexts, including primary T cells. XL413 stimulates HDR during a reversible slowing of S-phase that is unexplored for Cas9-induced HDR. We anticipate that XL413 and other such rationally developed inhibitors will be useful tools for gene modification

Sarcoidosis activates diverse transcriptional programs in bronchoalveolar lavage cells

Abstract Background Sarcoidosis is a multisystem immuno-inflammatory disorder of unknown etiology that most commonly involves the lungs. We hypothesized that an unbiased approach to identify pathways activated in bronchoalveolar lavage (BAL) cells can shed light on the pathogenesis of this complex disease. Methods We recruited 15 patients with various stages of sarcoidosis and 12 healthy controls. All subjects underwent bronchoscopy with lavage. For each subject, total RNA was extracted from BAL cells and hybridized to an Affymetrix U133A microarray. Rigorous statistical methods were applied to identify differential gene expression between subjects with sarcoidosis vs. controls. To better elucidate pathways differentially activated between these groups, we integrated network and gene set enrichment analyses of BAL cell transcriptional profiles. Results Sarcoidosis patients were either non-smokers or former smokers, all had lung involvement and only two were on systemic prednisone. Healthy controls were all non-smokers. Comparison of BAL cell gene expression between sarcoidosis and healthy subjects revealed over 1500 differentially expressed genes. Several previously described immune mediators, such as interferon gamma, were upregulated in the sarcoidosis subjects. Using an integrative computational approach we constructed a modular network of over 80 gene sets that were highly enriched in patients with sarcoidosis. Many of these pathways mapped to inflammatory and immune-related processes including adaptive immunity, T-cell signaling, graft vs. host disease, interleukin 12, 23 and 17 signaling. Additionally, we uncovered a close association between the proteasome machinery and adaptive immunity, highlighting a potentially important and targetable relationship in the pathobiology of sarcoidosis. Conclusions BAL cells in sarcoidosis are characterized by enrichment of distinct transcriptional programs involved in immunity and proteasomal processes. Our findings add to the growing evidence implicating alveolar resident immune effector cells in the pathogenesis of sarcoidosis and identify specific pathways whose activation may modulate disease progression

Genome-wide association and HLA fine-mapping studies identify risk loci and genetic pathways underlying allergic rhinitis

Allergic rhinitis is the most common clinical presentation of allergy, affecting 400 million people worldwide, with increasing incidence in westernized countries1,2. To elucidate the genetic architecture and understand the underlying disease mechanisms, we carried out a meta-analysis of allergic rhinitis in 59,762 cases and 152,358 controls of European ancestry and identified a total of 41 risk loci for allergic rhinitis, including 20 loci not previously associated with allergic rhinitis, which were confirmed in a replication phase of 60,720 cases and 618,527 controls. Functional annotation implicated genes involved in various immune pathways, and fine mapping of the HLA region suggested amino acid variants important for antigen binding. We further performed genome-wide association study (GWAS) analyses of allergic sensitization against inhalant allergens and nonallergic rhinitis, which suggested shared genetic mechanisms across rhinitis-related traits. Future studies of the identified loci and genes might identify novel targets for treatment and prevention of allergic rhinitis

Identification of Hypoxia-Regulated Proteins Using MALDI-Mass Spectrometry Imaging Combined with Quantitative Proteomics

Hypoxia is present in most solid tumors and is clinically correlated with increased metastasis and poor patient survival. While studies have demonstrated the role of hypoxia and hypoxia-regulated proteins in cancer progression, no attempts have been made to identify hypoxia-regulated proteins using quantitative proteomics combined with MALDI-mass spectrometry imaging (MALDI-MSI). Here we present a comprehensive hypoxic proteome study and are the first to investigate changes in situ using tumor samples. In vitro quantitative mass spectrometry analysis of the hypoxic proteome was performed on breast cancer cells using stable isotope labeling with amino acids in cell culture (SILAC). MS analyses were performed on laser-capture microdissected samples isolated from normoxic and hypoxic regions from tumors derived from the same cells used in vitro. MALDI-MSI was used in combination to investigate hypoxia-regulated protein localization within tumor sections. Here we identified more than 100 proteins, both novel and previously reported, that were associated with hypoxia. Several proteins were localized in hypoxic regions, as identified by MALDI-MSI. Visualization and data extrapolation methods for the in vitro SILAC data were also developed, and computational mapping of MALDI-MSI data to IHC results was applied for data validation. The results and limitations of the methodologies described are discussed. 2014 American Chemical Societ