Regional transcriptomic profiling reveals immune system enrichment in nonfailing atria and all chambers of the failing human heart.

The different chambers of the human heart demonstrate regional physiological traits and may be differentially affected during pathological remodeling, resulting in heart failure. Few previous studies, however, have characterized the different chambers at a transcriptomic level. We, therefore, conducted whole tissue RNA sequencing and gene set enrichment analysis of biopsies collected from the four chambers of adult failing (n = 8) and nonfailing (n = 11) human hearts. Atria and ventricles demonstrated distinct transcriptional patterns. When compared with nonfailing ventricles, the transcriptional pattern of nonfailing atria was enriched for many gene sets associated with cardiogenesis, the immune system and bone morphogenetic protein (BMP), transforming growth factor-ß (TGF-ß), MAPK/JNK, and Wnt signaling. Differences between failing and nonfailing hearts were also determined. The transcriptional pattern of failing atria was distinct compared with that of nonfailing atria and enriched for gene sets associated with the innate and adaptive immune system, TGF-ß/SMAD signaling, and changes in endothelial, smooth muscle cell, and cardiomyocyte physiology. Failing ventricles were also enriched for gene sets associated with the immune system. Based on the transcriptomic patterns, upstream regulators associated with heart failure were identified. These included many immune response factors predicted to be similarly activated for all chambers of failing hearts. In summary, the heart chambers demonstrate distinct transcriptional patterns that differ between failing and nonfailing hearts. Immune system signaling may be a hallmark of all four heart chambers in failing hearts and could constitute a novel therapeutic target.NEW & NOTEWORTHY The transcriptomic patterns of the four heart chambers were characterized in failing and nonfailing human hearts. Both nonfailing atria had distinct transcriptomic patterns characterized by cardiogenesis, the immune system and BMP/TGF-ß, MAPK/JNK, and Wnt signaling. Failing atria and ventricles were enriched for gene sets associated with the innate and adaptive immune system. Key upstream regulators associated with heart failure were identified, including activated immune response elements, which may constitute novel therapeutic targets.

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.

RNA-sequencing and mass-spectrometry proteomic time-series analysis of T-cell differentiation identified multiple splice variants models that predicted validated protein biomarkers in inflammatory diseases.

Profiling of mRNA expression is an important method to identify biomarkers but complicated by limited correlations between mRNA expression and protein abundance. We hypothesised that these correlations could be improved by mathematical models based on measuring splice variants and time delay in protein translation. We characterised time-series of primary human naïve CD4+ T cells during early T helper type 1 differentiation with RNA-sequencing and mass-spectrometry proteomics. We performed computational time-series analysis in this system and in two other key human and murine immune cell types. Linear mathematical mixed time delayed splice variant models were used to predict protein abundances, and the models were validated using out-of-sample predictions. Lastly, we re-analysed RNA-seq datasets to evaluate biomarker discovery in five T-cell associated diseases, further validating the findings for multiple sclerosis (MS) and asthma. The new models significantly out-performing models not including the usage of multiple splice variants and time delays, as shown in cross-validation tests. Our mathematical models provided more differentially expressed proteins between patients and controls in all five diseases. Moreover, analysis of these proteins in asthma and MS supported their relevance. One marker, sCD27, was validated in MS using two independent cohorts for evaluating response to treatment and disease prognosis. In summary, our splice variant and time delay models substantially improved the prediction of protein abundance from mRNA expression in three different immune cell types. The models provided valuable biomarker candidates, which were further validated in MS and asthma.

Deep neural network prediction of genome-wide transcriptome signatures - beyond the Black-box.

Prediction algorithms for protein or gene structures, including transcription factor binding from sequence information, have been transformative in understanding gene regulation. Here we ask whether human transcriptomic profiles can be predicted solely from the expression of transcription factors (TFs). We find that the expression of 1600 TFs can explain >95% of the variance in 25,000 genes. Using the light-up technique to inspect the trained NN, we find an over-representation of known TF-gene regulations. Furthermore, the learned prediction network has a hierarchical organization. A smaller set of around 125 core TFs could explain close to 80% of the variance. Interestingly, reducing the number of TFs below 500 induces a rapid decline in prediction performance. Next, we evaluated the prediction model using transcriptional data from 22 human diseases. The TFs were sufficient to predict the dysregulation of the target genes (rho = 0.61, P < 10-216). By inspecting the model, key causative TFs could be extracted for subsequent validation using disease-associated genetic variants. We demonstrate a methodology for constructing an interpretable neural network predictor, where analyses of the predictors identified key TFs that were inducing transcriptional changes during disease.

TFTenricher: a python toolbox for annotation enrichment analysis of transcription factor target genes.

BACKGROUND: Transcription factors (TFs) are the upstream regulators that orchestrate gene expression, and therefore a centrepiece in bioinformatics studies. While a core strategy to understand the biological context of genes and proteins includes annotation enrichment analysis, such as Gene Ontology term enrichment, these methods are not well suited for analysing groups of TFs. This is particularly true since such methods do not aim to include downstream processes, and given a set of TFs, the expected top ontologies would revolve around transcription processes. RESULTS: We present the TFTenricher, a Python toolbox that focuses specifically at identifying gene ontology terms, cellular pathways, and diseases that are over-represented among genes downstream of user-defined sets of human TFs. We evaluated the inference of downstream gene targets with respect to false positive annotations, and found an inference based on co-expression to best predict downstream processes. Based on these downstream genes, the TFTenricher uses some of the most common databases for gene functionalities, including GO, KEGG and Reactome, to calculate functional enrichments. By applying the TFTenricher to differential expression of TFs in 21 diseases, we found significant terms associated with disease mechanism, while the gene set enrichment analysis on the same dataset predominantly identified processes related to transcription. CONCLUSIONS AND AVAILABILITY: The TFTenricher package enables users to search for biological context in any set of TFs and their downstream genes. The TFTenricher is available as a Python 3 toolbox at https://github.com/rasma774/Tftenricher , under a GNU GPL license and with minimal dependencies.

Hybrid modelling for stroke care: Review and suggestions of new approaches for risk assessment and simulation of scenarios.

Stroke is an example of a complex and multi-factorial disease involving multiple organs, timescales, and disease mechanisms. To deal with this complexity, and to realize Precision Medicine of stroke, mathematical models are needed. Such approaches include: 1) machine learning, 2) bioinformatic network models, and 3) mechanistic models. Since these three approaches have complementary strengths and weaknesses, a hybrid modelling approach combining them would be the most beneficial. However, no concrete approach ready to be implemented for a specific disease has been presented to date. In this paper, we both review the strengths and weaknesses of the three approaches, and propose a roadmap for hybrid modelling in the case of stroke care. We focus on two main tasks needed for the clinical setting: a) For stroke risk calculation, we propose a new two-step approach, where non-linear mixed effects models and bioinformatic network models yield biomarkers which are used as input to a machine learning model and b) For simulation of care scenarios, we propose a new four-step approach, which revolves around iterations between simulations of the mechanistic models and imputations of non-modelled or non-measured variables. We illustrate and discuss the different approaches in the context of Precision Medicine for stroke.

Progesterone Dampens Immune Responses in In Vitro Activated CD4+ T Cells and Affects Genes Associated With Autoimmune Diseases That Improve During Pregnancy.

The changes in progesterone (P4) levels during and after pregnancy coincide with the temporary improvement and worsening of several autoimmune diseases like multiple sclerosis (MS) and rheumatoid arthritis (RA). Most likely immune-endocrine interactions play a major role in these pregnancy-induced effects. In this study, we used next generation sequencing to investigate the direct effects of P4 on CD4+ T cell activation, key event in pregnancy and disease. We report profound dampening effects of P4 on T cell activation, altering the gene and protein expression profile and reversing many of the changes induced during the activation. The transcriptomic changes induced by P4 were significantly enriched for genes associated with diseases known to be modulated during pregnancy such as MS, RA and psoriasis. STAT1 and STAT3 were significantly downregulated by P4 and their downstream targets were significantly enriched among the disease-associated genes. Several of these genes included well-known and disease-relevant cytokines, such as IL-12ß, CXCL10 and OSM, which were further validated also at the protein level using proximity extension assay. Our results extend the previous knowledge of P4 as an immune regulatory hormone and support its importance during pregnancy for regulating potentially detrimental immune responses towards the semi-allogenic fetus. Further, our results also point toward a potential role for P4 in the pregnancy-induced disease immunomodulation and highlight the need for further studies evaluating P4 as a future treatment option.

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.

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.

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.

Cross-talks via mTORC2 can explain enhanced activation in response to insulin in diabetic patients.

The molecular mechanisms of insulin resistance in Type 2 diabetes have been extensively studied in primary human adipocytes, and mathematical modelling has clarified the central role of attenuation of mammalian target of rapamycin (mTOR) complex 1 (mTORC1) activity in the diabetic state. Attenuation of mTORC1 in diabetes quells insulin-signalling network-wide, except for the mTOR in complex 2 (mTORC2)-catalysed phosphorylation of protein kinase B (PKB) at Ser473 (PKB-S473P), which is increased. This unique increase could potentially be explained by feedback and interbranch cross-talk signals. To examine if such mechanisms operate in adipocytes, we herein analysed data from an unbiased phosphoproteomic screen in 3T3-L1 adipocytes. Using a mathematical modelling approach, we showed that a negative signal from mTORC1-p70 S6 kinase (S6K) to rictor-mTORC2 in combination with a positive signal from PKB to SIN1-mTORC2 are compatible with the experimental data. This combined cross-branch signalling predicted an increased PKB-S473P in response to attenuation of mTORC1 - a distinguishing feature of the insulin resistant state in human adipocytes. This aspect of insulin signalling was then verified for our comprehensive model of insulin signalling in human adipocytes. Introduction of the cross-branch signals was compatible with all data for insulin signalling in human adipocytes, and the resulting model can explain all data network-wide, including the increased PKB-S473P in the diabetic state. Our approach was to first identify potential mechanisms in data from a phosphoproteomic screen in a cell line, and then verify such mechanisms in primary human cells, which demonstrates how an unbiased approach can support a direct knowledge-based study.