In skeletal muscle and neural crest cells, SMCHD1 regulates biological pathways relevant for Bosma syndrome and facioscapulohumeral dystrophy phenotype.

Many genetic syndromes are linked to mutations in genes encoding factors that guide chromatin organization. Among them, several distinct rare genetic diseases are linked to mutations in SMCHD1 that encodes the structural maintenance of chromosomes flexible hinge domain containing 1 chromatin-associated factor. In humans, its function as well as the impact of its mutations remains poorly defined. To fill this gap, we determined the episignature associated with heterozygous SMCHD1 variants in primary cells and cell lineages derived from induced pluripotent stem cells for Bosma arhinia and microphthalmia syndrome (BAMS) and type 2 facioscapulohumeral dystrophy (FSHD2). In human tissues, SMCHD1 regulates the distribution of methylated CpGs, H3K27 trimethylation and CTCF at repressed chromatin but also at euchromatin. Based on the exploration of tissues affected either in FSHD or in BAMS, i.e. skeletal muscle fibers and neural crest stem cells, respectively, our results emphasize multiple functions for SMCHD1, in chromatin compaction, chromatin insulation and gene regulation with variable targets or phenotypical outcomes. We concluded that in rare genetic diseases, SMCHD1 variants impact gene expression in two ways: (i) by changing the chromatin context at a number of euchromatin loci or (ii) by directly regulating some loci encoding master transcription factors required for cell fate determination and tissue differentiation.

Random walk with restart on multilayer networks: from node prioritisation to supervised link prediction and beyond.

BACKGROUND: Biological networks have proven invaluable ability for representing biological knowledge. Multilayer networks, which gather different types of nodes and edges in multiplex, heterogeneous and bipartite networks, provide a natural way to integrate diverse and multi-scale data sources into a common framework. Recently, we developed MultiXrank, a Random Walk with Restart algorithm able to explore such multilayer networks. MultiXrank outputs scores reflecting the proximity between an initial set of seed node(s) and all the other nodes in the multilayer network. We illustrate here the versatility of bioinformatics tasks that can be performed using MultiXrank. RESULTS: We first show that MultiXrank can be used to prioritise genes and drugs of interest by exploring multilayer networks containing interactions between genes, drugs, and diseases. In a second study, we illustrate how MultiXrank scores can also be used in a supervised strategy to train a binary classifier to predict gene-disease associations. The classifier performance are validated using outdated and novel gene-disease association for training and evaluation, respectively. Finally, we show that MultiXrank scores can be used to compute diffusion profiles and use them as disease signatures. We computed the diffusion profiles of more than 100 immune diseases using a multilayer network that includes cell-type specific genomic information. The clustering of the immune disease diffusion profiles reveals shared shared phenotypic characteristics. CONCLUSION: Overall, we illustrate here diverse applications of MultiXrank to showcase its versatility. We expect that this can lead to further and broader bioinformatics applications.

Tracking clusters of patients over time enables extracting information from medico-administrative databases.

CONTEXT: Identifying clusters (i.e., subgroups) of patients from the analysis of medico-administrative databases is particularly important to better understand disease heterogeneity. However, these databases contain different types of longitudinal variables which are measured over different follow-up periods, generating truncated data. It is therefore fundamental to develop clustering approaches that can handle this type of data. OBJECTIVE: We propose here cluster-tracking approaches to identify clusters of patients from truncated longitudinal data contained in medico-administrative databases. MATERIAL AND METHODS: We first cluster patients at each age. We then track the identified clusters over ages to construct cluster-trajectories. We compared our novel approaches with three classical longitudinal clustering approaches by calculating the silhouette score. As a use-case, we analyzed antithrombotic drugs used from 2008 to 2018 contained in the Échantillon Généraliste des Bénéficiaires (EGB), a French national cohort. RESULTS: Our cluster-tracking approaches allow us to identify several cluster-trajectories with clinical significance without any imputation of data. The comparison of the silhouette scores obtained with the different approaches highlights the better performances of the cluster-tracking approaches. CONCLUSION: The cluster-tracking approaches are a novel and efficient alternative to identify patient clusters from medico-administrative databases by taking into account their specificities.

orsum: a Python package for filtering and comparing enrichment analyses using a simple principle.

BACKGROUND: Enrichment analyses are widely applied to investigate lists of genes of interest. However, such analyses often result in long lists of annotation terms with high redundancy, making the interpretation and reporting difficult. Long annotation lists and redundancy also complicate the comparison of results obtained from different enrichment analyses. An approach to overcome these issues is using down-sized annotation collections composed of non-redundant terms. However, down-sized collections are generic and the level of detail may not fit the user's study. Other available approaches include clustering and filtering tools, which are based on similarity measures and thresholds that can be complicated to comprehend and set. RESULT: We propose orsum, a Python package to filter enrichment results. orsum can filter multiple enrichment results collectively and highlight common and specific annotation terms. Filtering in orsum is based on a simple principle: a term is discarded if there is a more significant term that annotates at least the same genes; the remaining more significant term becomes the representative term for the discarded term. This principle ensures that the main biological information is preserved in the filtered results while reducing redundancy. In addition, as the representative terms are selected from the original enrichment results, orsum outputs filtered terms tailored to the study. As a use case, we applied orsum to the enrichment analyses of four lists of genes, each associated with a neurodegenerative disease. CONCLUSION: orsum provides a comprehensible and effective way of filtering and comparing enrichment results. It is available at https://anaconda.org/bioconda/orsum .

A multi-objective genetic algorithm to find active modules in multiplex biological networks.

The identification of subnetworks of interest-or active modules-by integrating biological networks with molecular profiles is a key resource to inform on the processes perturbed in different cellular conditions. We here propose MOGAMUN, a Multi-Objective Genetic Algorithm to identify active modules in MUltiplex biological Networks. MOGAMUN optimizes both the density of interactions and the scores of the nodes (e.g., their differential expression). We compare MOGAMUN with state-of-the-art methods, representative of different algorithms dedicated to the identification of active modules in single networks. MOGAMUN identifies dense and high-scoring modules that are also easier to interpret. In addition, to our knowledge, MOGAMUN is the first method able to use multiplex networks. Multiplex networks are composed of different layers of physical and functional relationships between genes and proteins. Each layer is associated to its own meaning, topology, and biases; the multiplex framework allows exploiting this diversity of biological networks. We applied MOGAMUN to identify cellular processes perturbed in Facio-Scapulo-Humeral muscular Dystrophy, by integrating RNA-seq expression data with a multiplex biological network. We identified different active modules of interest, thereby providing new angles for investigating the pathomechanisms of this disease. Availability: MOGAMUN is available at https://github.com/elvanov/MOGAMUN and as a Bioconductor package at https://bioconductor.org/packages/release/bioc/html/MOGAMUN.html. Contact: anais.baudot@univ-amu.fr.

Random walk with restart on multiplex and heterogeneous biological networks.

Motivation: Recent years have witnessed an exponential growth in the number of identified interactions between biological molecules. These interactions are usually represented as large and complex networks, calling for the development of appropriated tools to exploit the functional information they contain. Random walk with restart (RWR) is the state-of-the-art guilt-by-association approach. It explores the network vicinity of gene/protein seeds to study their functions, based on the premise that nodes related to similar functions tend to lie close to each other in the networks. Results: In this study, we extended the RWR algorithm to multiplex and heterogeneous networks. The walk can now explore different layers of physical and functional interactions between genes and proteins, such as protein-protein interactions and co-expression associations. In addition, the walk can also jump to a network containing different sets of edges and nodes, such as phenotype similarities between diseases. We devised a leave-one-out cross-validation strategy to evaluate the algorithms abilities to predict disease-associated genes. We demonstrate the increased performances of the multiplex-heterogeneous RWR as compared to several random walks on monoplex or heterogeneous networks. Overall, our framework is able to leverage the different interaction sources to outperform current approaches. Finally, we applied the algorithm to predict candidate genes for the Wiedemann-Rautenstrauch syndrome, and to explore the network vicinity of the SHORT syndrome. Availability and implementation: The source code is available on GitHub at: https://github.com/alberto-valdeolivas/RWR-MH. In addition, an R package is freely available through Bioconductor at: http://bioconductor.org/packages/RandomWalkRestartMH/. Supplementary information: Supplementary data are available at Bioinformatics online.

Molecular Inverse Comorbidity between Alzheimer's Disease and Lung Cancer: New Insights from Matrix Factorization.

Matrix factorization (MF) is an established paradigm for large-scale biological data analysis with tremendous potential in computational biology. Here, we challenge MF in depicting the molecular bases of epidemiologically described disease-disease (DD) relationships. As a use case, we focus on the inverse comorbidity association between Alzheimer's disease (AD) and lung cancer (LC), described as a lower than expected probability of developing LC in AD patients. To this day, the molecular mechanisms underlying DD relationships remain poorly explained and their better characterization might offer unprecedented clinical opportunities. To this goal, we extend our previously designed MF-based framework for the molecular characterization of DD relationships. Considering AD-LC inverse comorbidity as a case study, we highlight multiple molecular mechanisms, among which we confirm the involvement of processes related to the immune system and mitochondrial metabolism. We then distinguish mechanisms specific to LC from those shared with other cancers through a pan-cancer analysis. Additionally, new candidate molecular players, such as estrogen receptor (ER), cadherin 1 (CDH1) and histone deacetylase (HDAC), are pinpointed as factors that might underlie the inverse relationship, opening the way to new investigations. Finally, some lung cancer subtype-specific factors are also detected, also suggesting the existence of heterogeneity across patients in the context of inverse comorbidity.

Molecular evidence for the inverse comorbidity between central nervous system disorders and cancers detected by transcriptomic meta-analyses.

There is epidemiological evidence that patients with certain Central Nervous System (CNS) disorders have a lower than expected probability of developing some types of Cancer. We tested here the hypothesis that this inverse comorbidity is driven by molecular processes common to CNS disorders and Cancers, and that are deregulated in opposite directions. We conducted transcriptomic meta-analyses of three CNS disorders (Alzheimer's disease, Parkinson's disease and Schizophrenia) and three Cancer types (Lung, Prostate, Colorectal) previously described with inverse comorbidities. A significant overlap was observed between the genes upregulated in CNS disorders and downregulated in Cancers, as well as between the genes downregulated in CNS disorders and upregulated in Cancers. We also observed expression deregulations in opposite directions at the level of pathways. Our analysis points to specific genes and pathways, the upregulation of which could increase the incidence of CNS disorders and simultaneously lower the risk of developing Cancer, while the downregulation of another set of genes and pathways could contribute to a decrease in the incidence of CNS disorders while increasing the Cancer risk. These results reinforce the previously proposed involvement of the PIN1 gene, Wnt and P53 pathways, and reveal potential new candidates, in particular related with protein degradation processes.

Discovery of Drug Synergies in Gastric Cancer Cells Predicted by Logical Modeling.

Discovery of efficient anti-cancer drug combinations is a major challenge, since experimental testing of all possible combinations is clearly impossible. Recent efforts to computationally predict drug combination responses retain this experimental search space, as model definitions typically rely on extensive drug perturbation data. We developed a dynamical model representing a cell fate decision network in the AGS gastric cancer cell line, relying on background knowledge extracted from literature and databases. We defined a set of logical equations recapitulating AGS data observed in cells in their baseline proliferative state. Using the modeling software GINsim, model reduction and simulation compression techniques were applied to cope with the vast state space of large logical models and enable simulations of pairwise applications of specific signaling inhibitory chemical substances. Our simulations predicted synergistic growth inhibitory action of five combinations from a total of 21 possible pairs. Four of the predicted synergies were confirmed in AGS cell growth real-time assays, including known effects of combined MEK-AKT or MEK-PI3K inhibitions, along with novel synergistic effects of combined TAK1-AKT or TAK1-PI3K inhibitions. Our strategy reduces the dependence on a priori drug perturbation experimentation for well-characterized signaling networks, by demonstrating that a model predictive of combinatorial drug effects can be inferred from background knowledge on unperturbed and proliferating cancer cells. Our modeling approach can thus contribute to preclinical discovery of efficient anticancer drug combinations, and thereby to development of strategies to tailor treatment to individual cancer patients.

The functional landscape of Hsp27 reveals new cellular processes such as DNA repair and alternative splicing and proposes novel anticancer targets.

Previously, we identified the stress-induced chaperone, Hsp27, as highly overexpressed in castration-resistant prostate cancer and developed an Hsp27 inhibitor (OGX-427) currently tested in phase I/II clinical trials as a chemosensitizing agent in different cancers. To better understand the Hsp27 poorly-defined cytoprotective functions in cancers and increase the OGX-427 pharmacological safety, we established the Hsp27-protein interaction network using a yeast two-hybrid approach and identified 226 interaction partners. As an example, we showed that targeting Hsp27 interaction with TCTP, a partner protein identified in our screen increases therapy sensitivity, opening a new promising field of research for therapeutic approaches that could decrease or abolish toxicity for normal cells. Results of an in-depth bioinformatics network analysis allying the Hsp27 interaction map into the human interactome underlined the multifunctional character of this protein. We identified interactions of Hsp27 with proteins involved in eight well known functions previously related to Hsp27 and uncovered 17 potential new ones, such as DNA repair and RNA splicing. Validation of Hsp27 involvement in both processes in human prostate cancer cells supports our system biology-predicted functions and provides new insights into Hsp27 roles in cancer cells.

System-level analysis of genes mutated in muscular dystrophies reveals a functional pattern associated with muscle weakness distribution.

Muscular dystrophies (MDs) are inherited genetic diseases causing weakness and degeneration of muscles. The distribution of muscle weakness differs between MDs, involving distal muscles or proximal muscles. While the mutations in most of the MD-associated genes lead to either distal or proximal onset, there are also genes whose mutations can cause both types of onsets. We hypothesized that the genes associated with different MD onsets code proteins with distinct cellular functions. To investigate this, we collected the MD-associated genes and assigned them to three onset groups: genes mutated only in distal onset dystrophies, genes mutated only in proximal onset dystrophies, and genes mutated in both types of onsets. We then systematically evaluated the cellular functions of these gene sets with computational strategies based on functional enrichment analysis and biological network analysis. Our analyses demonstrate that genes mutated in either distal or proximal onset MDs code proteins linked with two distinct sets of cellular processes. Interestingly, these two sets of cellular processes are relevant for the genes that are associated with both onsets. Moreover, the genes associated with both onsets display high centrality and connectivity in the network of muscular dystrophy genes. Our findings support the hypothesis that the proteins associated with distal or proximal onsets have distinct functional characteristics, whereas the proteins associated with both onsets are multifunctional.

EnrichNet: network-based gene set enrichment analysis.

MOTIVATION: Assessing functional associations between an experimentally derived gene or protein set of interest and a database of known gene/protein sets is a common task in the analysis of large-scale functional genomics data. For this purpose, a frequently used approach is to apply an over-representation-based enrichment analysis. However, this approach has four drawbacks: (i) it can only score functional associations of overlapping gene/proteins sets; (ii) it disregards genes with missing annotations; (iii) it does not take into account the network structure of physical interactions between the gene/protein sets of interest and (iv) tissue-specific gene/protein set associations cannot be recognized. RESULTS: To address these limitations, we introduce an integrative analysis approach and web-application called EnrichNet. It combines a novel graph-based statistic with an interactive sub-network visualization to accomplish two complementary goals: improving the prioritization of putative functional gene/protein set associations by exploiting information from molecular interaction networks and tissue-specific gene expression data and enabling a direct biological interpretation of the results. By using the approach to analyse sets of genes with known involvement in human diseases, new pathway associations are identified, reflecting a dense sub-network of interactions between their corresponding proteins. AVAILABILITY: EnrichNet is freely available at http://www.enrichnet.org. CONTACT: Natalio.Krasnogor@nottingham.ac.uk, reinhard.schneider@uni.lu or avalencia@cnio.es SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics Online.

A multilevel screening pipeline in zebrafish identifies therapeutic drugs for GAN.

Giant axonal neuropathy (GAN) is a fatal neurodegenerative disorder for which there is currently no treatment. Affecting the nervous system, GAN starts in infancy with motor deficits that rapidly evolve toward total loss of ambulation. Using the gan zebrafish model that reproduces the loss of motility as seen in patients, we conducted the first pharmacological screening for the GAN pathology. Here, we established a multilevel pipeline to identify small molecules restoring both the physiological and the cellular deficits in GAN. We combined behavioral, in silico, and high-content imaging analyses to refine our Hits to five drugs restoring locomotion, axonal outgrowth, and stabilizing neuromuscular junctions in the gan zebrafish. The postsynaptic nature of the drug's cellular targets provides direct evidence for the pivotal role the neuromuscular junction holds in the restoration of motility. Our results identify the first drug candidates that can now be integrated in a repositioning approach to fasten therapy for the GAN disease. Moreover, we anticipate both our methodological development and the identified hits to be of benefit to other neuromuscular diseases.

Mutated genes, pathways and processes in tumours.

Integration of the many available sources of cancer gene information--such as large-scale tumour-resequencing studies--identifies the 'usual suspect' genes, mutated in many tumour types, as well as different sets of mutated genes according to the specific tumour type. Scaling-up the analysis reveals that this large collection of mutated genes cluster into a smaller number of signalling pathways and processes. From this, we draw a map of the altered processes, and their combinations, in more than 10 tumours types. Literature searches identify pathways and processes that are covered sparsely in the literature, and invite the proposal of new hypotheses to investigate cancer initiation and progression.

Interactome mapping of the phosphatidylinositol 3-kinase-mammalian target of rapamycin pathway identifies deformed epidermal autoregulatory factor-1 as a new glycogen synthase kinase-3 interactor.

The phosphatidylinositol 3-kinase-mammalian target of rapamycin (PI3K-mTOR) pathway plays pivotal roles in cell survival, growth, and proliferation downstream of growth factors. Its perturbations are associated with cancer progression, type 2 diabetes, and neurological disorders. To better understand the mechanisms of action and regulation of this pathway, we initiated a large scale yeast two-hybrid screen for 33 components of the PI3K-mTOR pathway. Identification of 67 new interactions was followed by validation by co-affinity purification and exhaustive literature curation of existing information. We provide a nearly complete, functionally annotated interactome of 802 interactions for the PI3K-mTOR pathway. Our screen revealed a predominant place for glycogen synthase kinase-3 (GSK3) A and B and the AMP-activated protein kinase. In particular, we identified the deformed epidermal autoregulatory factor-1 (DEAF1) transcription factor as an interactor and in vitro substrate of GSK3A and GSK3B. Moreover, GSK3 inhibitors increased DEAF1 transcriptional activity on the 5-HT1A serotonin receptor promoter. We propose that DEAF1 may represent a therapeutic target of lithium and other GSK3 inhibitors used in bipolar disease and depression.

No paradox, no progress: inverse cancer comorbidity in people with other complex diseases.

In the past 5 years, several leading groups have attempted to explain why individuals with Down's syndrome have a reduced risk of many solid tumours and an increased risk of leukaemia and testicular cancer. Niels Bohr, the Danish physicist, noted that a paradox could initiate progress. We think that the paradox of a medical disorder protecting against cancer could be formalised in a new model of inverse cancer morbidity in people with other serious diseases. In this Personal View, we review evidence from epidemiological and clinical studies that supports a consistently lower than expected occurrence of cancer in patients with Down's syndrome, Parkinson's disease, schizophrenia, diabetes, Alzheimer's disease, multiple sclerosis, and anorexia nervosa. Intriguingly, most comorbidities are neuropsychiatric or CNS disorders. We provide a brief overview of evidence indicating genetic and molecular connections between cancer and these complex diseases. Inverse comorbidity could be a valuable model to investigate common or related pathways or processes and test new therapies, but, most importantly, to understand why certain people are protected from the malignancy.

Tracking Temporal Clusters from Patient Networks.

Creating homogeneous groups (clusters) of patients from medico-administrative databases provides a better understanding of health determinants. But in these databases, patients have truncated care pathways. We developed an approach based on patient networks to construct care trajectories from such truncated data. We tested this approach on antithrombotic treatments prescribed from 2008 to 2018 contained in the échantillon généraliste des bénéficiaires (EGB). We constructed a patient network for each patients' age (years from birth). We then applied the Markov clustering algorithm in each network. The care trajectories were finally constructed by matching clusters identified in two consecutive networks. We calculated the silhouette score to assess the performance of this network approach compared to three existing approaches. We identified 12 care trajectories that we were able to associate with pathologies. The best silhouette score was obtained for the network approach. Our approach allowed to highlight care trajectories taking into account the longitudinal, multidimensional and truncated nature of data from medico-administrative databases.

Facioscapulohumeral dystrophy weakened sarcomeric contractility is mimicked in induced pluripotent stem cells-derived innervated muscle fibres.

BACKGROUND: Facioscapulohumeral dystrophy (FSHD) is a late-onset autosomal dominant form of muscular dystrophy involving specific groups of muscles with variable weakness that precedes inflammatory response, fat infiltration, and muscle atrophy. As there is currently no cure for this disease, understanding and modelling the typical muscle weakness in FSHD remains a major milestone towards deciphering the disease pathogenesis as it will pave the way to therapeutic strategies aimed at correcting the functional muscular defect in patients. METHODS: To gain further insights into the specificity of the muscle alteration in this disease, we derived induced pluripotent stem cells from patients affected with Types 1 and 2 FSHD but also from patients affected with Bosma arhinia and microphthalmia. We differentiated these cells into contractile innervated muscle fibres and analysed their transcriptome by RNA Seq in comparison with cells derived from healthy donors. To uncover biological pathways altered in the disease, we applied MOGAMUN, a multi-objective genetic algorithm that integrates multiplex complex networks of biological interactions (protein-protein interactions, co-expression, and biological pathways) and RNA Seq expression data to identify active modules. RESULTS: We identified 132 differentially expressed genes that are specific to FSHD cells (false discovery rate < 0.05). In FSHD, the vast majority of active modules retrieved with MOGAMUN converges towards a decreased expression of genes encoding proteins involved in sarcomere organization (P value 2.63e-12 ), actin cytoskeleton (P value 9.4e-5 ), myofibril (P value 2.19e-12 ), actin-myosin sliding, and calcium handling (with P values ranging from 7.9e-35 to 7.9e-21 ). Combined with in vivo validations and functional investigations, our data emphasize a reduction in fibre contraction (P value < 0.0001) indicating that the muscle weakness that is typical of FSHD clinical spectrum might be associated with dysfunction of calcium release (P value < 0.0001), actin-myosin interactions, motor activity, mechano-transduction, and dysfunctional sarcomere contractility. CONCLUSIONS: Identification of biomarkers of FSHD muscle remain critical for understanding the process leading to the pathology but also for the definition of readouts to be used for drug design, outcome measures, and monitoring of therapies. The different pathways identified through a system biology approach have been largely overlooked in the disease. Overall, our work opens new perspectives in the definition of biomarkers able to define the muscle alteration but also in the development of novel strategies to improve muscle function as it provides functional parameters for active molecule screening.

miR-376a-3p and miR-376b-3p overexpression in Hutchinson-Gilford progeria fibroblasts inhibits cell proliferation and induces premature senescence.

Hutchinson-Gilford progeria syndrome (HGPS) is a rare genetic disorder, in which an abnormal and toxic protein called progerin, accumulates in cell nuclei, leading to major cellular defects. Among them, chromatin remodeling drives gene expression changes, including miRNA dysregulation. In our study, we evaluated miRNA expression profiles in HGPS and control fibroblasts. We identified an enrichment of overexpressed miRNAs belonging to the 14q32.2-14q32.3 miRNA cluster. Using 3D FISH, we demonstrated that overexpression of these miRNAs is associated with chromatin remodeling at this specific locus in HGPS fibroblasts. We then focused on miR-376b-3p and miR-376a-3p, both overexpressed in HGPS fibroblasts. We demonstrated that their induced overexpression in control fibroblasts decreases cell proliferation and increases senescence, whereas their inhibition in HGPS fibroblasts rescues proliferation defects and senescence and decreases progerin accumulation. By targeting these major processes linked to premature aging, these two miRNAs may play a pivotal role in the pathophysiology of HGPS.

TopoGSA: network topological gene set analysis.

UNLABELLED: TopoGSA (Topology-based Gene Set Analysis) is a web-application dedicated to the computation and visualization of network topological properties for gene and protein sets in molecular interaction networks. Different topological characteristics, such as the centrality of nodes in the network or their tendency to form clusters, can be computed and compared with those of known cellular pathways and processes. AVAILABILITY: Freely available at http://www.infobiotics.net/topogsa.