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.

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.

Thermodynamic model of gene regulation for the Or59b olfactory receptor in Drosophila.

Complex eukaryotic promoters normally contain multiple cis-regulatory sequences for different transcription factors (TFs). The binding patterns of the TFs to these sites, as well as the way the TFs interact with each other and with the RNA polymerase (RNAp), lead to combinatorial problems rarely understood in detail, especially under varying epigenetic conditions. The aim of this paper is to build a model describing how the main regulatory cluster of the olfactory receptor Or59b drives transcription of this gene in Drosophila. The cluster-driven expression of this gene is represented as the equilibrium probability of RNAp being bound to the promoter region, using a statistical thermodynamic approach. The RNAp equilibrium probability is computed in terms of the occupancy probabilities of the single TFs of the cluster to the corresponding binding sites, and of the interaction rules among TFs and RNAp, using experimental data of Or59b expression to tune the model parameters. The model reproduces correctly the changes in RNAp binding probability induced by various mutation of specific sites and epigenetic modifications. Some of its predictions have also been validated in novel experiments.

The phototransduction machinery in the rod outer segment has a strong efficacy gradient.

Rod photoreceptors consist of an outer segment (OS) and an inner segment. Inside the OS a biochemical machinery transforms the rhodopsin photoisomerization into electrical signal. This machinery has been treated as and is thought to be homogenous with marginal inhomogeneities. To verify this assumption, we developed a methodology based on special tapered optical fibers (TOFs) to deliver highly localized light stimulations. By using these TOFs, specific regions of the rod OS could be stimulated with spots of light highly confined in space. As the TOF is moved from the OS base toward its tip, the amplitude of saturating and single photon responses decreases, demonstrating that the efficacy of the transduction machinery is not uniform and is 5-10 times higher at the base than at the tip. This gradient of efficacy of the transduction machinery is attributed to a progressive depletion of the phosphodiesterase along the rod OS. Moreover we demonstrate that, using restricted spots of light, the duration of the photoresponse along the OS does not increase linearly with the light intensity as with diffuse light.

Metabolic Adaptation Processes That Converge to Optimal Biomass Flux Distributions.

In simple organisms like E.coli, the metabolic response to an external perturbation passes through a transient phase in which the activation of a number of latent pathways can guarantee survival at the expenses of growth. Growth is gradually recovered as the organism adapts to the new condition. This adaptation can be modeled as a process of repeated metabolic adjustments obtained through the resilencings of the non-essential metabolic reactions, using growth rate as selection probability for the phenotypes obtained. The resulting metabolic adaptation process tends naturally to steer the metabolic fluxes towards high growth phenotypes. Quite remarkably, when applied to the central carbon metabolism of E.coli, it follows that nearly all flux distributions converge to the flux vector representing optimal growth, i.e., the solution of the biomass optimization problem turns out to be the dominant attractor of the metabolic adaptation process.

A rate-distortion theory for gene regulatory networks and its application to logic gate consistency.

MOTIVATION: A gene regulatory network in which the modes (activation/inhibition) of the transcriptional regulations are known and in which gene expression assumes boolean values can be treated as a system of linear equations over a binary field, i.e. as a constraint satisfaction problem for an information code. RESULTS: For currently available gene networks, we show in this article that the distortion associated with the corresponding information code is much lower than expected from null models, and that it is close to (when not lower than) the Shannon bound determined by the rate-distortion theorem. This corresponds to saying that the distribution of regulatory modes is highly atypical in the networks, and that this atypicality greatly helps in avoiding contradictory transcriptional actions. Choosing a boolean formalism to represent the gene networks, we also show how to formulate criteria for the selection of gates that maximize the compatibility with the empirical information available on the transcriptional regulatory modes. Proceeding in this way, we obtain in particular that non-canalizing gates are upper-bounded by the distortion, and hence that the boolean gene networks are more canalizing than expected from null models.

Detection of transcriptional triggers in the dynamics of microbial growth: application to the respiratorily versatile bacterium Shewanella oneidensis.

The capacity of microorganisms to respond to variable external conditions requires a coordination of environment-sensing mechanisms and decision-making regulatory circuits. Here, we seek to understand the interplay between these two processes by combining high-throughput measurement of time-dependent mRNA profiles with a novel computational approach that searches for key genetic triggers of transcriptional changes. Our approach helped us understand the regulatory strategies of a respiratorily versatile bacterium with promising bioenergy and bioremediation applications, Shewanella oneidensis, in minimal and rich media. By comparing expression profiles across these two conditions, we unveiled components of the transcriptional program that depend mainly on the growth phase. Conversely, by integrating our time-dependent data with a previously available large compendium of static perturbation responses, we identified transcriptional changes that cannot be explained solely by internal network dynamics, but are rather triggered by specific genes acting as key mediators of an environment-dependent response. These transcriptional triggers include known and novel regulators that respond to carbon, nitrogen and oxygen limitation. Our analysis suggests a sequence of physiological responses, including a coupling between nitrogen depletion and glycogen storage, partially recapitulated through dynamic flux balance analysis, and experimentally confirmed by metabolite measurements. Our approach is broadly applicable to other systems.

Computing global structural balance in large-scale signed social networks.

Structural balance theory affirms that signed social networks (i.e., graphs whose signed edges represent friendly/hostile interactions among individuals) tend to be organized so as to avoid conflictual situations, corresponding to cycles of negative parity. Using an algorithm for ground-state calculation in large-scale Ising spin glasses, in this paper we compute the global level of balance of very large online social networks and verify that currently available networks are indeed extremely balanced. This property is explainable in terms of the high degree of skewness of the sign distributions on the nodes of the graph. In particular, individuals linked by a large majority of negative edges create mostly "apparent disorder," rather than true "frustration."

From populations to networks: Relating diversity indices and frustration in signed graphs.

Diversity indices of quadratic type, such as fractionalization and Simpson index, are measures of heterogeneity in a population. Even though they are univariate, they have an intrinsic bivariate interpretation as encounters among the elements of the population. In the paper, it is shown that this leads naturally to associate populations to weakly balanced signed networks. In particular, the frustration of such signed networks is shown to be related to fractionalization by a closed-form expression. This expression allows to simplify drastically the calculation of frustration for weakly balanced signed graphs.

Partial inhibition and bilevel optimization in flux balance analysis.

MOTIVATION: Within Flux Balance Analysis, the investigation of complex subtasks, such as finding the optimal perturbation of the network or finding an optimal combination of drugs, often requires to set up a bilevel optimization problem. In order to keep the linearity and convexity of these nested optimization problems, an ON/OFF description of the effect of the perturbation (i.e. Boolean variable) is normally used. This restriction may not be realistic when one wants, for instance, to describe the partial inhibition of a reaction induced by a drug. RESULTS: In this paper we present a formulation of the bilevel optimization which overcomes the oversimplified ON/OFF modeling while preserving the linear nature of the problem. A case study is considered: the search of the best multi-drug treatment which modulates an objective reaction and has the minimal perturbation on the whole network. The drug inhibition is described and modulated through a convex combination of a fixed number of Boolean variables. The results obtained from the application of the algorithm to the core metabolism of E.coli highlight the possibility of finding a broader spectrum of drug combinations compared to a simple ON/OFF modeling. CONCLUSIONS: The method we have presented is capable of treating partial inhibition inside a bilevel optimization, without loosing the linearity property, and with reasonable computational performances also on large metabolic networks. The more fine-graded representation of the perturbation allows to enlarge the repertoire of synergistic combination of drugs for tasks such as selective perturbation of cellular metabolism. This may encourage the use of the approach also for other cases in which a more realistic modeling is required.

Decompositions of large-scale biological systems based on dynamical properties.

MOTIVATION: Given a large-scale biological network represented as an influence graph, in this article we investigate possible decompositions of the network aimed at highlighting specific dynamical properties. RESULTS: The first decomposition we study consists in finding a maximal directed acyclic subgraph of the network, which dynamically corresponds to searching for a maximal open-loop subsystem of the given system. Another dynamical property investigated is strong monotonicity. We propose two methods to deal with this property, both aimed at decomposing the system into strongly monotone subsystems, but with different structural characteristics: one method tends to produce a single large strongly monotone component, while the other typically generates a set of smaller disjoint strongly monotone subsystems. AVAILABILITY: Original heuristics for the methods investigated are described in the article. CONTACT: altafini@sissa.it

Quantifying leadership in climate negotiations: A social power game.

We consider complex multistage multiagent negotiation processes such as those occurring at climate conferences and ask ourselves how can an agent maximize its social power, intended as influence over the outcome of the negotiation. This question can be framed as a strategic game played over an opinion dynamics model, in which the action of an agent consists in stubbornly defending its own opinion. We show that for consensus-seeking opinion dynamics models in which the interaction weights are uniform, the optimal action obeys to an early mover advantage principle, i.e. the agents behaving stubbornly in the early phases of the negotiations achieve the highest social power. When looking at data collected from the climate change negotiations going on at the United Nations Framework Convention on Climate Change, we find evidence of the use of the early mover strategy. Furthermore, we show that the social powers computed through our model correlate very well with the perceived leadership roles assessed through independent survey data, especially when non-uniform weights incorporating economical and demographic factors are considered.

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.

A dynamical feedback model for adaptation in the olfactory transduction pathway.

Olfactory transduction exhibits two distinct types of adaptation, which we denote multipulse and step adaptation. In terms of measured transduction current, multipulse adaptation appears as a decrease in the amplitude of the second of two consecutive responses when the olfactory neuron is stimulated with two brief pulses. Step adaptation occurs in response to a sustained steplike stimulation and is characterized by a return to a steady-state current amplitude close to the prestimulus value, after a transient peak. In this article, we formulate a dynamical model of the olfactory transduction pathway, which includes the kinetics of the CNG channels, the concentration of Ca ions flowing through them, and the Ca-complexes responsible for the regulation. Based on this model, a common dynamical explanation for the two types of adaptation is suggested. We show that both forms of adaptation can be well described using different time constants for the kinetics of Ca ions (faster) and the kinetics of the feedback mechanisms (slower). The model is validated on experimental data collected in voltage-clamp conditions using different techniques and animal species.

A system-level approach for deciphering the transcriptional response to prion infection.

MOTIVATION: Deciphering the response of a complex biological system to an insulting event, at the gene expression level, requires adopting theoretical models that are more sophisticated than a one-to-one comparison (i.e. t-test). Here, we investigate the ability of a novel reverse engineering approach (System Response Inference) to unveil non-obvious transcriptional signatures of the system response induced by prion infection. RESULTS: To this end, we analyze previously published gene expression data, from which we extrapolate a putative full-scale model of transcriptional gene-gene dependencies in the mouse central nervous system. Then, we use this nominal model to interpret the gene expression changes caused by prion replication, aiming at selecting the genes primarily influenced by this perturbation. Our method sheds light on the mode of action of prions by identifying key transcripts that are the most likely to be responsible for the overall transcriptional rearrangement from a nominal regulatory network. As a first result of our inference, we have been able to predict known targets of prions (i.e. PrP(C)) and to unveil the potential role of previously unsuspected genes. CONTACT: altafini@sissa.it SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multi-omics protein-coding units as massively parallel Bayesian networks: Empirical validation of causality structure.

In this article we use high-throughput epigenomics, transcriptomics, and proteomics data to construct fine-graded models of the "protein-coding units" gathering all transcript isoforms and chromatin accessibility peaks associated with more than 4000 genes in humans. Each protein-coding unit has the structure of a directed acyclic graph (DAG) and can be represented as a Bayesian network. The factorization of the joint probability distribution induced by the DAGs imposes a number of conditional independence relationships among the variables forming a protein-coding unit, corresponding to the missing edges in the DAGs. We show that a large fraction of these conditional independencies are indeed verified by the data. Factors driving this verification appear to be the structural and functional annotation of the transcript isoforms, as well as a notion of structural balance (or frustration-free) of the corresponding sample correlation graph, which naturally leads to reduction of correlation (and hence to independence) upon conditioning.

Relating balance and conditional independence in graphical models.

When data are available for all nodes of a Gaussian graphical model, then, it is possible to use sample correlations and partial correlations to test to what extent the conditional independencies that encode the structure of the model are indeed verified by the data. In this paper, we give a heuristic rule useful in such a validation process: When the correlation subgraph involved in a conditional independence is balanced (i.e., all its cycles have an even number of negative edges), then a partial correlation is usually a contraction of the corresponding correlation, which often leads to conditional independence. In particular, the contraction rule can be made rigorous if we look at concentration subgraphs rather than correlation subgraphs. The rule is applied to real data for elementary gene regulatory motifs.

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.

A signed network perspective on the government formation process in parliamentary democracies.

In parliamentary democracies, government negotiations talks following a general election can sometimes be a long and laborious process. In order to explain this phenomenon, in this paper we use structural balance theory to represent a multiparty parliament as a signed network, with edge signs representing alliances and rivalries among parties. We show that the notion of frustration, which quantifies the amount of "disorder" encoded in the signed graph, correlates very well with the duration of the government negotiation talks. For the 29 European countries considered in this study, the average correlation between frustration and government negotiation talks ranges between 0.42 and 0.69, depending on what information is included in the edges of the signed network. Dynamical models of collective decision-making over signed networks with varying frustration are proposed to explain this correlation.

Achieving consensus in multilateral international negotiations: The case study of the 2015 Paris Agreement on climate change.

The purpose of this paper is to propose a dynamical model describing the achievement of the 2015 Paris Agreement on climate change. To represent the complex, decade-long, multiparty negotiation process that led to the accord, we use a two time scale dynamical model. The short time scale corresponds to the discussion process occurring at each meeting and is represented as a Friedkin-Johnsen model, a dynamical multiparty model in which the parties show stubbornness, i.e., tend to defend their positions during the discussion. The long time scale behavior is determined by concatenating multiple Friedkin-Johnsen models (one for each meeting). The proposed model, tuned on real data extracted from the Paris Agreement meetings, achieves consensus on a time horizon similar to that of the real negotiations. Remarkably, the model is also able to identify a series of parties that exerted a key leadership role in the Paris Agreement negotiation process.