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

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

<i>Schizosaccharomyces pombe</i> Pol II transcription elongation factor ELL functions as part of a rudimentary super elongation complex

ELL family transcription factors activate the overall rate of RNA polymerase II (Pol II) transcription elongation by binding directly to Pol II and suppressing its tendency to pause. In metazoa, ELL regulates Pol II transcription elongation as part of a large multisubunit complex referred to as the Super Elongation Complex (SEC), which includes P-TEFb and EAF, AF9 or ENL, and an AFF family protein. Although orthologs of ELL and EAF have been identified in lower eukaryotes including Schizosaccharomyces pombe, it has been unclear whether SEClike complexes function in lower eukaryotes. In this report, we describe isolation from S. pombe of an ELL-containing complex with features of a rudimentary SEC. This complex includes S. pombe Ell1, Eaf1, and a previously uncharacterized protein we designate Ell1 binding protein 1 (Ebp1), which is distantly related to metazoan AFF family members. Like the metazoan SEC, this S. pombe ELL complex appears to function broadly in Pol II transcription. Interestingly, it appears to have a particularly important role in regulating genes involved in cell separation

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

Ribosomal DNA, tri- and bi-partite pericentromeres in the permanent translocation heterozygote Rhoeo spathacea

High- and low-stringency FISH and base-specific fluorescence were performed on the permanent translocation heterozygote Rhoeo spathacea (2n = 12). Our results indicate that 45S rDNA arrays, rDNA-related sequences and other GC-rich DNA fraction(s) are located within the pericentromeric regions of all twelve chromosomes, usually colocalizing with the chromomycin A3-positive bands. Homogenization of the pericentromeric regions appears to result from the concerted spread of GC-rich sequences, with differential amplification likely. We found new 5S rDNA patterns, which suggest a variability in the breakpoints and in the consequent chromosome reorganizations. It was found that the large 5S rDNA locus residing on each of the 8E and 9E arms consisted of two smaller loci. On each of the two chromosome arms 3b and 4b, in addition to the major subtelomeric 5S rDNA locus, a new minor locus was found interstitially about 40% along the arm length. The arrangement of cytotogenetic landmarks and chromosome arm measurements are discussed with regard to genome repatterning in Rhoeo

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

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

Specific MRI abnormalities reveal severe perrault syndrome due to CLPP defects

In establishing a genetic diagnosis in heterogeneous neurological disease, clinical characterization and whole exome sequencing (WES) go hand-in-hand. Clinical data are essential, not only to guide WES variant selection and define the clinical severity of a genetic defect but also to identify other patients with defects in the same gene. In an infant patient with sensorineural hearing loss, psychomotor retardation, and epilepsy, WES resulted in identification of a novel homozygous CLPP frameshift mutation (c.21delA). Based on the gene defect and clinical symptoms, the diagnosis Perrault syndrome type 3 (PRLTS3) was established. The patient's brain-MRI revealed specific abnormalities of the subcortical and deep cerebral white matter and the middle blade of the corpus callosum, which was used to identify similar patients in the Amsterdam brain-MRI database, containing over 3000 unclassified leukoencephalopathy cases. In three unrelated patients with similar MRI abnormalities the CLPP gene was sequenced, and in two of them novel missense mutations were identified together with a large deletion that covered part of the CLPP gene on the other allele. The severe neurological and MRI abnormalities in these young patients were due to the drastic impact of the CLPP mutations, correlating with the variation in clinical manifestations among previously reported patients. Our data show that similarity in brain-MRI patterns can be used to identify novel PRLTS3 patients, especially during early disease stages, when only part of the disease manifestations are present. This seems especially applicable to the severely affected cases in which CLPP function is drastically affected and MRI abnormalities are pronounced

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