Cholinergic regulation of vascular endothelial function by human ChAT+ T cells.

Endothelial dysfunction and impaired vasodilation are linked with adverse cardiovascular events. T lymphocytes expressing choline acetyltransferase (ChAT), the enzyme catalyzing biosynthesis of the vasorelaxant acetylcholine (ACh), regulate vasodilation and are integral to the cholinergic antiinflammatory pathway in an inflammatory reflex in mice. Here, we found that human T cell ChAT mRNA expression was induced by T cell activation involving the PI3K signaling cascade. Mechanistically, we identified that ChAT mRNA expression was induced following the attenuation of RE-1 Silencing Transcription factor REST-mediated methylation of the ChAT promoter, and that ChAT mRNA expression levels were up-regulated by GATA3 in human T cells. In functional experiments, T cell-derived ACh increased endothelial nitric oxide-synthase activity, promoted vasorelaxation, and reduced vascular endothelial activation and promoted barrier integrity by a cholinergic mechanism. Further, we observed that survival in a cohort of patients with severe circulatory failure correlated with their relative frequency of ChAT +CD4+ T cells in blood. These findings on ChAT+ human T cells provide a mechanism for cholinergic immune regulation of vascular endothelial function in human inflammation.

NCAE: data-driven representations using a deep network-coherent DNA methylation autoencoder identify robust disease and risk factor signatures.

Precision medicine relies on the identification of robust disease and risk factor signatures from omics data. However, current knowledge-driven approaches may overlook novel or unexpected phenomena due to the inherent biases in biological knowledge. In this study, we present a data-driven signature discovery workflow for DNA methylation analysis utilizing network-coherent autoencoders (NCAEs) with biologically relevant latent embeddings. First, we explored the architecture space of autoencoders trained on a large-scale pan-tissue compendium (n = 75 272) of human epigenome-wide association studies. We observed the emergence of co-localized patterns in the deep autoencoder latent space representations that corresponded to biological network modules. We determined the NCAE configuration with the strongest co-localization and centrality signals in the human protein interactome. Leveraging the NCAE embeddings, we then trained interpretable deep neural networks for risk factor (aging, smoking) and disease (systemic lupus erythematosus) prediction and classification tasks. Remarkably, our NCAE embedding-based models outperformed existing predictors, revealing novel DNA methylation signatures enriched in gene sets and pathways associated with the studied condition in each case. Our data-driven biomarker discovery workflow provides a generally applicable pipeline to capture relevant risk factor and disease information. By surpassing the limitations of knowledge-driven methods, our approach enhances the understanding of complex epigenetic processes, facilitating the development of more effective diagnostic and therapeutic strategies.

A DNA Methylation Signature from Buccal Swabs to Identify Tuberculosis Infection.

BACKGROUND: Tuberculosis (TB) is amongst the largest infectious causes of death worldwide and there is a need for a time- and resource-effective diagnostic method. In this novel and exploratory study, we show the potential of using buccal swabs to collect human DNA and investigate the DNA methylation (DNAm) signatures as a diagnostic tool for TB. METHODS: Buccal swabs were collected from pulmonary TB patients (n= 7), TB exposed (n= 7), and controls (n= 9) in Sweden. Using Illumina MethylationEPIC array the DNAm status was determined. RESULTS: We identified 5644 significant differentially methylated CpG sites between the patients and controls. Performing the analysis on a validation cohort of samples collected in Kenya and Peru (patients, n=26; exposed, n=9; control, n=10) confirmed the DNAm signature. We identified a TB consensus disease module, significantly enriched in TB-associated genes. Lastly, we used machine learning to identify a panel of seven CpG sites discriminative for TB and developed a TB classifier. In the validation cohort the classifier performed with an AUC of 0.94, sensitivity of 0.92, and specificity of 1. CONCLUSION: In summary, the result from this study shows clinical implications of using DNAm signatures from buccal swabs to explore new diagnostic strategies for TB.

TET2 as a tumor suppressor and therapeutic target in T-cell acute lymphoblastic leukemia.

Pediatric T-cell acute lymphoblastic leukemia (T-ALL) is an aggressive malignancy resulting from overproduction of immature T-cells in the thymus and is typified by widespread alterations in DNA methylation. As survival rates for relapsed T-ALL remain dismal (10 to 25%), development of targeted therapies to prevent relapse is key to improving prognosis. Whereas mutations in the DNA demethylating enzyme TET2 are frequent in adult T-cell malignancies, TET2 mutations in T-ALL are rare. Here, we analyzed RNA-sequencing data of 321 primary T-ALLs, 20 T-ALL cell lines, and 25 normal human tissues, revealing that TET2 is transcriptionally repressed or silenced in 71% and 17% of T-ALL, respectively. Furthermore, we show that TET2 silencing is often associated with hypermethylation of the TET2 promoter in primary T-ALL. Importantly, treatment with the DNA demethylating agent, 5-azacytidine (5-aza), was significantly more toxic to TET2-silenced T-ALL cells and resulted in stable re-expression of the TET2 gene. Additionally, 5-aza led to up-regulation of methylated genes and human endogenous retroviruses (HERVs), which was further enhanced by the addition of physiological levels of vitamin C, a potent enhancer of TET activity. Together, our results clearly identify 5-aza as a potential targeted therapy for TET2-silenced T-ALL.

Prominent epigenetic and transcriptomic changes in CD4+ and CD8+ T cells during and after pregnancy in women with multiple sclerosis and controls.

BACKGROUND: Multiple sclerosis (MS) is a neuroinflammatory disease in which pregnancy leads to a temporary amelioration in disease activity as indicated by the profound decrease in relapses rate during the 3rd trimester of pregnancy. CD4+ and CD8+ T cells are implicated in MS pathogenesis as being key regulators of inflammation and brain lesion formation. Although Tcells are prime candidates for the pregnancy-associated improvement of MS, the precise mechanisms are yet unclear, and in particular, a deep characterization of the epigenetic and transcriptomic events that occur in peripheral T cells during pregnancy in MS is lacking. METHODS: Women with MS and healthy controls were longitudinally sampled before, during (1st, 2nd and 3rd trimesters) and after pregnancy. DNA methylation array and RNA sequencing were performed on paired CD4+ and CD8+ T cells samples. Differential analysis and network-based approaches were used to analyze the global dynamics of epigenetic and transcriptomic changes. RESULTS: Both DNA methylation and RNA sequencing revealed a prominent regulation, mostly peaking in the 3rd trimester and reversing post-partum, thus mirroring the clinical course with improvement followed by a worsening in disease activity. This rebound pattern was found to represent a general adaptation of the maternal immune system, with only minor differences between MS and controls. By using a network-based approach, we highlighted several genes at the core of this pregnancy-induced regulation, which were found to be enriched for genes and pathways previously reported to be involved in MS. Moreover, these pathways were enriched for in vitro stimulated genes and pregnancy hormones targets. CONCLUSION: This study represents, to our knowledge, the first in-depth investigation of the methylation and expression changes in peripheral CD4+ and CD8+ T cells during pregnancy in MS. Our findings indicate that pregnancy induces profound changes in peripheral T cells, in both MS and healthy controls, which are associated with the modulation of inflammation and MS activity.

Differential effects of estradiol and progesterone on human T cell activation in vitro.

Estradiol (E2) and progesterone (P4) are steroid hormones important for the regulation of immune responses during pregnancy. Their increasing levels coincide with an improvement of T cell-mediated diseases such as multiple sclerosis (MS). Although immune-endocrine interactions are involved in this phenomenon, the relative contribution of hormones is not known. We here report a direct comparison of E2- and P4-mediated effects on human CD4+ T cells, key cells in immune regulation. T cells were stimulated to obtain different activation levels and exposed to a broad range of hormone concentrations. Activation level was assessed by CD69/CD25 expression by flow cytometry, and secreted proteins (n = 196) were measured in culture supernatants using proximity extension assay and electrochemiluminescence immunoassay. We found that in low activated cells, pregnancy-relevant E2 concentrations increased activation and the secretion of several immune- and inflammation-related proteins. P4, on the other hand, showed a biphasic pattern, where serum-related concentrations upregulated activation and protein secretion while placenta-relevant concentrations induced a prominent dampening irrespective of the initial activation level. Our results demonstrate the importance of P4 as a major hormone in the immune modulation of T cells during pregnancy and emphasize the need to further evaluate its potency in the treatment of diseases like MS.

Using high-throughput multi-omics data to investigate structural balance in elementary gene regulatory network motifs.

MOTIVATION: The simultaneous availability of ATAC-seq and RNA-seq experiments allows to obtain a more in-depth knowledge on the regulatory mechanisms occurring in gene regulatory networks. In this article, we highlight and analyze two novel aspects that leverage on the possibility of pairing RNA-seq and ATAC-seq data. Namely we investigate the causality of the relationships between transcription factors, chromatin and target genes and the internal consistency between the two omics, here measured in terms of structural balance in the sample correlations along elementary length-3 cycles. RESULTS: We propose a framework that uses the a priori knowledge on the data to infer elementary causal regulatory motifs (namely chains and forks) in the network. It is based on the notions of conditional independence and partial correlation, and can be applied to both longitudinal and non-longitudinal data. Our analysis highlights a strong connection between the causal regulatory motifs that are selected by the data and the structural balance of the underlying sample correlation graphs: strikingly, >97% of the selected regulatory motifs belong to a balanced subgraph. This result shows that internal consistency, as measured by structural balance, is close to a necessary condition for 3-node regulatory motifs to satisfy causality rules. AVAILABILITY AND IMPLEMENTATION: The analysis was carried out in MATLAB and the code can be found at https://github.com/albertozenere/Multi-omics-elementary-regulatory-motifs. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ComHub: Community predictions of hubs in gene regulatory networks.

BACKGROUND: Hub transcription factors, regulating many target genes in gene regulatory networks (GRNs), play important roles as disease regulators and potential drug targets. However, while numerous methods have been developed to predict individual regulator-gene interactions from gene expression data, few methods focus on inferring these hubs. RESULTS: We have developed ComHub, a tool to predict hubs in GRNs. ComHub makes a community prediction of hubs by averaging over predictions by a compendium of network inference methods. Benchmarking ComHub against the DREAM5 challenge data and two independent gene expression datasets showed a robust performance of ComHub over all datasets. CONCLUSIONS: In contrast to other evaluated methods, ComHub consistently scored among the top performing methods on data from different sources. Lastly, we implemented ComHub to work with both predefined networks and to perform stand-alone network inference, which will make the method generally applicable.

A validated generally applicable approach using the systematic assessment of disease modules by GWAS reveals a multi-omic module strongly associated with risk factors in multiple sclerosis.

BACKGROUND: There exist few, if any, practical guidelines for predictive and falsifiable multi-omic data integration that systematically integrate existing knowledge. Disease modules are popular concepts for interpreting genome-wide studies in medicine but have so far not been systematically evaluated and may lead to corroborating multi-omic modules. RESULT: We assessed eight module identification methods in 57 previously published expression and methylation studies of 19 diseases using GWAS enrichment analysis. Next, we applied the same strategy for multi-omic integration of 20 datasets of multiple sclerosis (MS), and further validated the resulting module using both GWAS and risk-factor-associated genes from several independent cohorts. Our benchmark of modules showed that in immune-associated diseases modules inferred from clique-based methods were the most enriched for GWAS genes. The multi-omic case study using MS data revealed the robust identification of a module of 220 genes. Strikingly, most genes of the module were differentially methylated upon the action of one or several environmental risk factors in MS (n = 217, P = 10- 47) and were also independently validated for association with five different risk factors of MS, which further stressed the high genetic and epigenetic relevance of the module for MS. CONCLUSIONS: We believe our analysis provides a workflow for selecting modules and our benchmark study may help further improvement of disease module methods. Moreover, we also stress that our methodology is generally applicable for combining and assessing the performance of multi-omic approaches for complex diseases.

LiPLike: towards gene regulatory network predictions of high certainty.

MOTIVATION: High correlation in expression between regulatory elements is a persistent obstacle for the reverse-engineering of gene regulatory networks. If two potential regulators have matching expression patterns, it becomes challenging to differentiate between them, thus increasing the risk of false positive identifications. RESULTS: To allow for gene regulation predictions of high confidence, we propose a novel method, the Linear Profile Likelihood (LiPLike), that assumes a regression model and iteratively searches for interactions that cannot be replaced by a linear combination of other predictors. To compare the performance of LiPLike with other available inference methods, we benchmarked LiPLike using three independent datasets from the Dialogue on Reverse Engineering Assessment and Methods 5 (DREAM5) network inference challenge. We found that LiPLike could be used to stratify predictions of other inference tools, and when applied to the predictions of DREAM5 participants, we observed an average improvement in accuracy of >140% compared to individual methods. Furthermore, LiPLike was able to independently predict networks better than all DREAM5 participants when applied to biological data. When predicting the Escherichia coli network, LiPLike had an accuracy of 0.38 for the top-ranked 100 interactions, whereas the corresponding DREAM5 consensus model yielded an accuracy of 0.11. AVAILABILITY AND IMPLEMENTATION: We made LiPLike available to the community as a Python toolbox, available at https://gitlab.com/Gustafsson-lab/liplike. We believe that LiPLike will be used for high confidence predictions in studies where individual model interactions are of high importance, and to remove false positive predictions made by other state-of-the-art gene-gene regulation prediction tools. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

MODifieR: an Ensemble R Package for Inference of Disease Modules from Transcriptomics Networks.

MOTIVATION: Complex diseases are due to the dense interactions of many disease-associated factors that dysregulate genes that in turn form the so-called disease modules, which have shown to be a powerful concept for understanding pathological mechanisms. There exist many disease module inference methods that rely on somewhat different assumptions, but there is still no gold standard or best-performing method. Hence, there is a need for combining these methods to generate robust disease modules. RESULTS: We developed MODule IdentiFIER (MODifieR), an ensemble R package of nine disease module inference methods from transcriptomics networks. MODifieR uses standardized input and output allowing the possibility to combine individual modules generated from these methods into more robust disease-specific modules, contributing to a better understanding of complex diseases. AVAILABILITY AND IMPLEMENTATION: MODifieR is available under the GNU GPL license and can be freely downloaded from https://gitlab.com/Gustafsson-lab/MODifieR and as a Docker image from https://hub.docker.com/r/ddeweerd/modifier. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Pre- and postnatal Lactobacillus reuteri treatment alters DNA methylation of infant T helper cells.

BACKGROUND: Perinatal childhood exposures, including probiotic supplementation, may affect epigenetic modifications and impact on immune maturation and allergy development. The aim of this study was to assess the effects of pre- and postnatal Lactobacillus reuteri supplementation on DNA methylation in relation to immune maturation and allergy development. METHODS: DNA methylation patterns were investigated for allergy-related T helper subsets using a locus-specific method and at a genome-wide scale using the Illumina 450K array. From a randomised, double-blind, placebo-controlled allergy prevention trial with pre- and postnatal probiotic supplementation, CD4+ T helper cells were obtained at birth (from cord blood), and 12 and 24 months of age (total (placebo/probiotics); locus-specific method: CB = 32 (17/15), 12 months = 24 (9/15), 24 months = 35 (15/20); Illumina: CB = 19 (10/9), 12 months = 10 (6/4), 24 months = 19(11/8)). RESULTS: Comparing probiotics to placebo, the greatest genome-wide differential DNA methylation was observed at birth, where the majority of sites were hypomethylated, indicating transcriptional accessibility in the probiotic group. Bioinformatic analyses, including network analyses, revealed a module containing 91 genes, enriched for immune-related pathways such as chemotaxis, PI3K-Akt, MAPK and TGF-ß signalling. A majority of the module genes were associated with atopic manifestations (OR = 1.43, P = 2.4 × 10-6 ), and a classifier built on this model could predict allergy development (AUC = 0.78, P = 3.0 × 10e-3 ). Pathways such as IFN-γ signalling and T-cell activation were more hypermethylated at birth compared with later in life in both intervention groups over time, in line with DNA methylation patterns in the IFNG locus obtained by the locus-specific methodology. CONCLUSION: Maternal L. reuteri supplementation during pregnancy alters DNA methylation patterns in CD4+ T cells towards enhanced immune activation at birth, which may affect immune maturation and allergy development.

LASSIM-A network inference toolbox for genome-wide mechanistic modeling.

Recent technological advancements have made time-resolved, quantitative, multi-omics data available for many model systems, which could be integrated for systems pharmacokinetic use. Here, we present large-scale simulation modeling (LASSIM), which is a novel mathematical tool for performing large-scale inference using mechanistically defined ordinary differential equations (ODE) for gene regulatory networks (GRNs). LASSIM integrates structural knowledge about regulatory interactions and non-linear equations with multiple steady state and dynamic response expression datasets. The rationale behind LASSIM is that biological GRNs can be simplified using a limited subset of core genes that are assumed to regulate all other gene transcription events in the network. The LASSIM method is implemented as a general-purpose toolbox using the PyGMO Python package to make the most of multicore computers and high performance clusters, and is available at https://gitlab.com/Gustafsson-lab/lassim. As a method, LASSIM works in two steps, where it first infers a non-linear ODE system of the pre-specified core gene expression. Second, LASSIM in parallel optimizes the parameters that model the regulation of peripheral genes by core system genes. We showed the usefulness of this method by applying LASSIM to infer a large-scale non-linear model of naïve Th2 cell differentiation, made possible by integrating Th2 specific bindings, time-series together with six public and six novel siRNA-mediated knock-down experiments. ChIP-seq showed significant overlap for all tested transcription factors. Next, we performed novel time-series measurements of total T-cells during differentiation towards Th2 and verified that our LASSIM model could monitor those data significantly better than comparable models that used the same Th2 bindings. In summary, the LASSIM toolbox opens the door to a new type of model-based data analysis that combines the strengths of reliable mechanistic models with truly systems-level data. We demonstrate the power of this approach by inferring a mechanistically motivated, genome-wide model of the Th2 transcription regulatory system, which plays an important role in several immune related diseases.

DNA methylation changes separate allergic patients from healthy controls and may reflect altered CD4+ T-cell population structure.

Altered DNA methylation patterns in CD4(+) T-cells indicate the importance of epigenetic mechanisms in inflammatory diseases. However, the identification of these alterations is complicated by the heterogeneity of most inflammatory diseases. Seasonal allergic rhinitis (SAR) is an optimal disease model for the study of DNA methylation because of its well-defined phenotype and etiology. We generated genome-wide DNA methylation (N(patients) = 8, N(controls) = 8) and gene expression (N(patients) = 9, Ncontrols = 10) profiles of CD4(+) T-cells from SAR patients and healthy controls using Illumina's HumanMethylation450 and HT-12 microarrays, respectively. DNA methylation profiles clearly and robustly distinguished SAR patients from controls, during and outside the pollen season. In agreement with previously published studies, gene expression profiles of the same samples failed to separate patients and controls. Separation by methylation (N(patients) = 12, N(controls) = 12), but not by gene expression (N(patients) = 21, N(controls) = 21) was also observed in an in vitro model system in which purified PBMCs from patients and healthy controls were challenged with allergen. We observed changes in the proportions of memory T-cell populations between patients (N(patients) = 35) and controls (N(controls) = 12), which could explain the observed difference in DNA methylation. Our data highlight the potential of epigenomics in the stratification of immune disease and represents the first successful molecular classification of SAR using CD4(+) T cells.

scDrugPrio: a framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases.

BACKGROUND: Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs. METHODS: Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs. RESULTS: scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment. CONCLUSIONS: We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package ( https://github.com/SDTC-CPMed/scDrugPrio ).

A genetic-epigenetic interplay at 1q21.1 locus underlies CHD1L-mediated vulnerability to primary progressive multiple sclerosis.

Multiple Sclerosis (MS) is a heterogeneous inflammatory and neurodegenerative disease with an unpredictable course towards progressive disability. Treating progressive MS is challenging due to limited insights into the underlying mechanisms. We examined the molecular changes associated with primary progressive MS (PPMS) using a cross-tissue (blood and post-mortem brain) and multilayered data (genetic, epigenetic, transcriptomic) from independent cohorts. In PPMS, we found hypermethylation of the 1q21.1 locus, controlled by PPMS-specific genetic variations and influencing the expression of proximal genes (CHD1L, PRKAB2) in the brain. Evidence from reporter assay and CRISPR/dCas9 experiments supports a causal link between methylation and expression and correlation network analysis further implicates these genes in PPMS brain processes. Knock-down of CHD1L in human iPSC-derived neurons and knock-out of chd1l in zebrafish led to developmental and functional deficits of neurons. Thus, several lines of evidence suggest a distinct genetic-epigenetic-transcriptional interplay in the 1q21.1 locus potentially contributing to PPMS pathogenesis.

A comprehensive mechanistic model of adipocyte signaling with layers of confidence.

Adipocyte signaling, normally and in type 2 diabetes, is far from fully understood. We have earlier developed detailed dynamic mathematical models for several well-studied, partially overlapping, signaling pathways in adipocytes. Still, these models only cover a fraction of the total cellular response. For a broader coverage of the response, large-scale phosphoproteomic data and systems level knowledge on protein interactions are key. However, methods to combine detailed dynamic models with large-scale data, using information about the confidence of included interactions, are lacking. We have developed a method to first establish a core model by connecting existing models of adipocyte cellular signaling for: (1) lipolysis and fatty acid release, (2) glucose uptake, and (3) the release of adiponectin. Next, we use publicly available phosphoproteome data for the insulin response in adipocytes together with prior knowledge on protein interactions, to identify phosphosites downstream of the core model. In a parallel pairwise approach with low computation time, we test whether identified phosphosites can be added to the model. We iteratively collect accepted additions into layers and continue the search for phosphosites downstream of these added layers. For the first 30 layers with the highest confidence (311 added phosphosites), the model predicts independent data well (70-90% correct), and the predictive capability gradually decreases when we add layers of decreasing confidence. In total, 57 layers (3059 phosphosites) can be added to the model with predictive ability kept. Finally, our large-scale, layered model enables dynamic simulations of systems-wide alterations in adipocytes in type 2 diabetes.

scDrugPrio: A framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases.

Background: Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs. Methods: Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs. Results: scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment. Conclusion: We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package (https://github.com/SDTC-CPMed/scDrugPrio).

Proteomics reveal biomarkers for diagnosis, disease activity and long-term disability outcomes in multiple sclerosis.

Sensitive and reliable protein biomarkers are needed to predict disease trajectory and personalize treatment strategies for multiple sclerosis (MS). Here, we use the highly sensitive proximity-extension assay combined with next-generation sequencing (Olink Explore) to quantify 1463 proteins in cerebrospinal fluid (CSF) and plasma from 143 people with early-stage MS and 43 healthy controls. With longitudinally followed discovery and replication cohorts, we identify CSF proteins that consistently predicted both short- and long-term disease progression. Lower levels of neurofilament light chain (NfL) in CSF is superior in predicting the absence of disease activity two years after sampling (replication AUC = 0.77) compared to all other tested proteins. Importantly, we also identify a combination of 11 CSF proteins (CXCL13, LTA, FCN2, ICAM3, LY9, SLAMF7, TYMP, CHI3L1, FYB1, TNFRSF1B and NfL) that predict the severity of disability worsening according to the normalized age-related MS severity score (replication AUC = 0.90). The identification of these proteins may help elucidate pathogenetic processes and might aid decisions on treatment strategies for persons with MS.

Combining gene expression microarray- and cluster analysis with sequence-based predictions to identify regulators of IL-13 in allergy.

The Th2 cytokine IL-13 plays a key role in allergy, by regulating IgE, airway hyper secretion, eosinophils and mast cells. In this study, we aimed to identify novel transcription factors (TFs) that potentially regulated IL-13. We analyzed Th2 polarized naïve T cells from four different blood donors with gene expression microarrays to find clusters of genes that were correlated or anti-correlated with IL13. These clusters were further filtered, by selecting genes that were functionally related. In these clusters, we identified three transcription factors (TFs) that were predicted to regulate the expression of IL13, namely CEBPB, E2F6 and AHR. siRNA mediated knockdowns of these TFs in naïve polarized T cells showed significant increases of IL13, following knockdown of CEBPB and E2F6, but not AHR. This suggested an inhibitory role of CEBPB and E2F6 in the regulation of IL13 and allergy. This was supported by analysis of E2F6, but not CEBPB, in allergen-challenged CD4+ T cells from six allergic patients and six healthy controls, which showed decreased expression of E2F6 in patients. In summary, our findings indicate an inhibitory role of E2F6 in the regulation of IL-13 and allergy. The analytical approach may be generally applicable to elucidate the complex regulatory patterns in Th2 cell polarization and allergy.