Multiscale networks in multiple sclerosis
Keith E. Kennedy,Nicole Kerlero de Rosbo,Antonio Uccelli,Maria Cellerino,Federico Ivaldi,Paola Contini,Raffaele De Palma,Hanne F. Harbo,Tone Berge,Steffan D. Bos,Einar A. Høgestøl,Synne Brune-Ingebretsen,Sigrid A. de Rodez Benavent,Friedemann Paul,Alexander U. Brandt,Priscilla Bäcker-Koduah,Janina Behrens,Joseph Kuchling,Susanna Asseyer,Michael Scheel,Claudia Chien,Hanna Zimmermann,Seyedamirhosein Motamedi,Josef Kauer-Bonin,Julio Saez-Rodriguez,Melanie Rinas,Leonidas G. Alexopoulos,Magi Andorra,Sara Llufriu,Albert Saiz,Yolanda Blanco,Eloy Martinez-Heras,Elisabeth Solana,Irene Pulido-Valdeolivas,Elena H. Martinez-Lapiscina,Jordi Garcia-Ojalvo,Pablo Villoslada
DOI: https://doi.org/10.1371/journal.pcbi.1010980
2024-02-09
PLoS Computational Biology
Abstract:Complex diseases such as Multiple Sclerosis (MS) cover a wide range of biological scales, from genes and proteins to cells and tissues, up to the full organism. In fact, any phenotype for an organism is dictated by the interplay among these scales. We conducted a multilayer network analysis and deep phenotyping with multi-omics data (genomics, phosphoproteomics and cytomics), brain and retinal imaging, and clinical data, obtained from a multicenter prospective cohort of 328 patients and 90 healthy controls. Multilayer networks were constructed using mutual information for topological analysis, and Boolean simulations were constructed using Pearson correlation to identified paths within and among all layers. The path more commonly found from the Boolean simulations connects protein MK03, with total T cells, the thickness of the retinal nerve fiber layer (RNFL), and the walking speed. This path contains nodes involved in protein phosphorylation, glial cell differentiation, and regulation of stress-activated MAPK cascade, among others. Specific paths identified were subsequently analyzed by flow cytometry at the single-cell level. Combinations of several proteins (GSK3AB, HSBP1 or RS6) and immune cells (Th17, Th1 non-classic, CD8, CD8 Treg, CD56 neg, and B memory) were part of the paths explaining the clinical phenotype. The advantage of the path identified from the Boolean simulations is that it connects information about these known biological pathways with the layers at higher scales (retina damage and disability). Overall, the identified paths provide a means to connect the molecular aspects of MS with the overall phenotype. Complex diseases such as Multiple Sclerosis (MS) involve the contribution of a wide range of biological processes. We conducted a systems biology study of MS based on network analysis and deep phenotyping in a prospective cohort of patients with clinical, imaging, genetics, and omics assessments. The genes, proteins and cell paths explained variation in central nervous system damage, and in metrics of disease severity. Such multilayer paths explain the different phenotypes of the disease and can be developed as biomarkers of MS.
biochemical research methods,mathematical & computational biology