Identification of bifurcation transitions in biological regulatory networks using Answer-Set Programming

Background

Numerous cellular differentiation processes can be captured using discrete qualitative models of biological regulatory networks. These models describe the temporal evolution of the state of the network subject to different competing transitions, potentially leading the system to different attractors. This paper focusses on the formal identification of states and transitions that are crucial for preserving or pre-empting the reachability of a given behaviour.

Methods

In the context of non-deterministic automata networks, we propose a static identification of so-called bifurcations, i.e., transitions after which a given goal is no longer reachable. Such transitions are naturally good candidates for controlling the occurrence of the goal, notably by modulating their propensity. Our method combines Answer-Set Programming with static analysis of reachability properties to provide an under-approximation of all the existing bifurcations.

Results

We illustrate our discrete bifurcation analysis on several models of biological systems, for which we identify transitions which impact the reachability of given long-term behaviour. In particular, we apply our implementation on a regulatory network among hundreds of biological species, supporting the scalability of our approach.

Conclusions

Our method allows a formal and scalable identification of transitions which are responsible for the lost of capability to reach a given state. It can be applied to any asynchronous automata networks, which encompass Boolean and multi-valued models. An implementation is provided as part of the Pint software, available at http://loicpauleve.name/pint.

     

An efficient algorithm for identifying primary phenotype attractors of a large-scale Boolean network

Background

Boolean network modeling has been widely used to model large-scale biomolecular regulatory networks as it can describe the essential dynamical characteristics of complicated networks in a relatively simple way. When we analyze such Boolean network models, we often need to find out attractor states to investigate the converging state features that represent particular cell phenotypes. This is, however, very difficult (often impossible) for a large network due to computational complexity.

Results

There have been some attempts to resolve this problem by partitioning the original network into smaller subnetworks and reconstructing the attractor states by integrating the local attractors obtained from each subnetwork. But, in many cases, the partitioned subnetworks are still too large and such an approach is no longer useful. So, we have investigated the fundamental reason underlying this problem and proposed a novel efficient way of hierarchically partitioning a given large network into smaller subnetworks by focusing on some attractors corresponding to a particular phenotype of interest instead of considering all attractors at the same time. Using the definition of attractors, we can have a simplified update rule with fixed state values for some nodes. The resulting subnetworks were small enough to find out the corresponding local attractors which can be integrated for reconstruction of the global attractor states of the original large network.

Conclusions

The proposed approach can substantially extend the current limit of Boolean network modeling for converging state analysis of biological networks.

     

High-fidelity discrete modeling of the HPA axis: a study of regulatory plasticity in biology

Background

The hypothalamic-pituitary-adrenal (HPA) axis is a central regulator of stress response and its dysfunction has been associated with a broad range of complex illnesses including Gulf War Illness (GWI) and Chronic Fatigue Syndrome (CFS). Though classical mathematical approaches have been used to model HPA function in isolation, its broad regulatory interactions with immune and central nervous function are such that the biological fidelity of simulations is undermined by the limited availability of reliable parameter estimates.

Method

Here we introduce and apply a generalized discrete formalism to recover multiple stable regulatory programs of the HPA axis using little more than connectivity between physiological components. This simple discrete model captures cyclic attractors such as the circadian rhythm by applying generic constraints to a minimal parameter set; this is distinct from Ordinary Differential Equation (ODE) models, which require broad and precise parameter sets. Parameter tuning is accomplished by decomposition of the overall regulatory network into isolated sub-networks that support cyclic attractors. Network behavior is simulated using a novel asynchronous updating scheme that enforces priority with memory within and between physiological compartments.

Results

Consistent with much more complex conventional models of the HPA axis, this parsimonious framework supports two cyclic attractors, governed by higher and lower levels of cortisol respectively. Importantly, results suggest that stress may remodel the stability landscape of this system, favoring migration from one stable circadian cycle to the other. Access to each regime is dependent on HPA axis tone, captured here by the tunable parameters of the multi-valued logic. Likewise, an idealized glucocorticoid receptor blocker alters the regulatory topology such that maintenance of persistently low cortisol levels is rendered unstable, favoring a return to normal circadian oscillation in both cortisol and glucocorticoid receptor expression.

Conclusion

These results emphasize the significance of regulatory connectivity alone and how regulatory plasticity may be explored using simple discrete logic and minimal data compared to conventional methods.

     

Optimization of fecal sample preparation for untargeted LC-HRMS based metabolomics

Introduction

Collecting feces is easy. It offers direct outcome to endogenous and microbial metabolites.

Objectives

In a context of lack of consensus about fecal sample preparation, especially in animal species, we developed a robust protocol allowing untargeted LC-HRMS fingerprinting.

Methods

The conditions of extraction (quantity, preparation, solvents, dilutions) were investigated in bovine feces.

Results

A rapid and simple protocol involving feces extraction with methanol (1/3, M/V) followed by centrifugation and a step filtration (10 kDa) was developed.

Conclusion

The workflow generated repeatable and informative fingerprints for robust metabolome characterization.

     

A lost opportunity for science: journals promote data sharing in metabolomics but do not enforce it

Introduction

Data sharing is being increasingly required by journals and has been heralded as a solution to the ‘replication crisis’.

Objectives

(i) Review data sharing policies of journals publishing the most metabolomics papers associated with open data and (ii) compare these journals’ policies to those that publish the most metabolomics papers.

Methods

A PubMed search was used to identify metabolomics papers. Metabolomics data repositories were manually searched for linked publications.

Results

Journals that support data sharing are not necessarily those with the most papers associated to open metabolomics data.

Conclusion

Further efforts are required to improve data sharing in metabolomics.

     

Synchronous early primary adenocarcinoma of both rectum and gallbladder. Report of a case

Background

Synchronous early primary cancers are rare and in addition synchronous adenocarcinoma of both rectum and gallbladder is extremely rare.

Case report

We report an unusual case of synchronous early primary adenocarcinoma of rectum and gallbladder. The patient was a 72-year-old woman with complaints of bloody stools and constipation. An endoscopy revealed adenocarcinoma of the lower rectum. A through preoperative investigation showed also cholelithiasis. The patient underwent abdominoperineal resection and cholecystectomy. The histopathological diagnosis was well to middle differentiate adenocarcinoma of the gallbladder (T2, N0, M0; stage II) and middle differentiate adenocarcinoma of the rectum (T2, N0, M0; stage II).

Conclusion

For the cases of extracolonic primary cancer associated with colorectal primary carcinoma, Warren and Gates' diagnostic criteria are used. All patients with colorectal carcinoma, should undergo a throughout preoperative examination to exclude the possibility of synchronous early primary cancers.

     

Quantitative analysis of tetrahydrofolate metabolites from <Emphasis Type="Italic">clostridium autoethanogenum</Emphasis>

Introduction

Quantification of tetrahydrofolates (THFs), important metabolites in the Wood–Ljungdahl pathway (WLP) of acetogens, is challenging given their sensitivity to oxygen.

Objective

To develop a simple anaerobic protocol to enable reliable THFs quantification from bioreactors.

Methods

Anaerobic cultures were mixed with anaerobic acetonitrile for extraction. Targeted LC–MS/MS was used for quantification.

Results

Tetrahydrofolates can only be quantified if sampled anaerobically. THF levels showed a strong correlation to acetyl-CoA, the end product of the WLP.

Conclusion

Our method is useful for relative quantification of THFs across different growth conditions. Absolute quantification of THFs requires the use of labelled standards.

     

A methodology for elucidating regulatory mechanisms leading to changes in lipid profiles

Introduction

It is difficult to elucidate the metabolic and regulatory factors causing lipidome perturbations.

Objectives

This work simplifies this process.

Methods

A method has been developed to query an online holistic lipid metabolic network (of 7923 metabolites) to extract the pathways that connect the input list of lipids.

Results

The output enables pathway visualisation and the querying of other databases to identify potential regulators. When used to a study a plasma lipidome dataset of polycystic ovary syndrome, 14 enzymes were identified, of which 3 are linked to ELAVL1—an mRNA stabiliser.

Conclusion

This method provides a simplified approach to identifying potential regulators causing lipid-profile perturbations.

     

Pseudo current density maps of electrophysiological heart,nerve or brain function and their physical basis

Background

In recent years the visualization of biomagnetic measurement data by so-called pseudo current density maps or Hosaka-Cohen (HC) transformations became popular.

Methods

The physical basis of these intuitive maps is clarified by means of analytically solvable problems.

Results

Examples in magnetocardiography, magnetoencephalography and magnetoneurography demonstrate the usefulness of this method.

Conclusion

Hardware realizations of the HC-transformation and some similar transformations are discussed which could advantageously support cross-platform comparability of biomagnetic measurements.

     

Using an NMR metabolomics approach to investigate the pathogenicity of amyloid-beta and alpha-synuclein

Introduction

The pathogenicity at differing points along the aggregation pathway of many fibril-forming proteins associated with neurodegenerative diseases is unclear. Understanding the effect of different aggregation states of these proteins on cellular processes is essential to enhance understanding of diseases and provide future options for diagnosis and therapeutic intervention.

Objectives

To establish a robust method to probe the metabolic changes of neuronal cells and use it to monitor cellular response to challenge with three amyloidogenic proteins associated with neurodegenerative diseases in different aggregation states.

Method

Neuroblastoma SH-SY5Y cells were employed to design a robust routine system to perform a statistically rigorous NMR metabolomics study into cellular effects of sub-toxic levels of alpha-synuclein, amyloid-beta 40 and amyloid-beta 42 in monomeric, oligomeric and fibrillar conformations.

Results

This investigation developed a rigorous model to monitor intracellular metabolic profiles of neuronal cells through combination of existing methods. This model revealed eight key metabolites that are altered when neuroblastoma cells are challenged with proteins in different aggregation states. Metabolic pathways associated with lipid metabolism, neurotransmission and adaptation to oxidative stress and inflammation are the predominant contributors to the cellular variance and intracellular metabolite levels. The observed metabolite changes for monomer and oligomer challenge may represent cellular effort to counteract the pathogenicity of the challenge, whereas fibrillar challenge is indicative of system shutdown. This implies that although markers of stress are more prevalent under oligomeric challenge the fibrillar response suggests a more toxic environment.

Conclusion

This approach is applicable to any cell type that can be cultured in a laboratory (primary or cell line) as a method of investigating how protein challenge affects signalling pathways, providing additional understanding as to the role of protein aggregation in neurodegenerative disease initiation and progression.

     

A critical review of PASBio's argument structures for biomedical verbs

Background

Propositional representations of biomedical knowledge are a critical component of most aspects of semantic mining in biomedicine. However, the proper set of propositions has yet to be determined. Recently, the PASBio project proposed a set of propositions and argument structures for biomedical verbs. This initial set of representations presents an opportunity for evaluating the suitability of predicate-argument structures as a scheme for representing verbal semantics in the biomedical domain. Here, we quantitatively evaluate several dimensions of the initial PASBio propositional structure repository.

Results

We propose a number of metrics and heuristics related to arity, role labelling, argument realization, and corpus coverage for evaluating large-scale predicate-argument structure proposals. We evaluate the metrics and heuristics by applying them to PASBio 1.0.

Conclusion

PASBio demonstrates the suitability of predicate-argument structures for representing aspects of the semantics of biomedical verbs. Metrics related to theta-criterion violations and to the distribution of arguments are able to detect flaws in semantic representations, given a set of predicate-argument structures and a relatively small corpus annotated with them.

     

Metabolomic profiling of maternal hair suggests rapid development of intrahepatic cholestasis of pregnancy

Introduction

Intrahepatic cholestasis of pregnancy (ICP) is a common maternal liver disease; development can result in devastating consequences, including sudden fetal death and stillbirth. Currently, recognition of ICP only occurs following onset of clinical symptoms.

Objective

Investigate the maternal hair metabolome for predictive biomarkers of ICP.

Methods

The maternal hair metabolome (gestational age of sampling between 17 and 41 weeks) of 38 Chinese women with ICP and 46 pregnant controls was analysed using gas chromatography–mass spectrometry.

Results

Of 105 metabolites detected in hair, none were significantly associated with ICP.

Conclusion

Hair samples represent accumulative environmental exposure over time. Samples collected at the onset of ICP did not reveal any metabolic shifts, suggesting rapid development of the disease.

     

Extraction of natural products from bark of <Emphasis Type="Italic">Betula pendula</Emphasis> using ionic liquids

Introduction

Aqueous–methanol mixtures have successfully been applied to extract a broad range of metabolites from plant tissue. However, a certain amount of material remains insoluble.

Objectives

To enlarge the metabolic compendium, two ionic liquids were selected to extract the methanol insoluble part of trunk from Betula pendula.

Methods

The extracted compounds were analyzed by LC/MS and GC/MS.

Results

The results show that 1-butyl-3-methylimidazolium acetate (IL-Ac) predominantly resulted in fatty acids, whereas 1-ethyl-3-methylimidazolium tosylate (IL-Tos) mostly yielded phenolic structures. Interestingly, bark yielded more ionic liquid soluble metabolites compared to interior wood.

Conclusion

From this one can conclude that the application of ionic liquids may expand the metabolic snapshot.

     

Untargeted metabolomics suffers from incomplete raw data processing

Introduction

Untargeted metabolomics is a powerful tool for biological discoveries. To analyze the complex raw data, significant advances in computational approaches have been made, yet it is not clear how exhaustive and reliable the data analysis results are.

Objectives

Assessment of the quality of raw data processing in untargeted metabolomics.

Methods

Five published untargeted metabolomics studies, were reanalyzed.

Results

Omissions of at least 50 relevant compounds from the original results as well as examples of representative mistakes were reported for each study.

Conclusion

Incomplete raw data processing shows unexplored potential of current and legacy data.

     

Preoperative detection of pleural adhesions by respiratory dynamic computed tomography

Background

Video-assisted thoracic surgery (VATS) plays an important role in thoracic surgery because it is less invasive. However, the existence of severe pleural adhesions may make VATS difficult and complicated. The aim of this study was to assess the utility of inspiration and expiration computed tomography (respiratory dynamic CT (RD-CT)) in evaluation of pleural adhesions preoperatively.

Methods

RD-CT was performed on 107 patients undergoing thoracotomies (both VATS and open). We assessed synchronous motion during respiration on RD-CT. Comparing the results of RD-CT and intraoperative findings, we assessed the utility of preoperative evaluation.

Results

A negative correlation between sliding score and adhesion grade was revealed. Sliding score in adhesion negative patients was significantly higher than that in adhesion positive patients (P?<?0.0001). The sensitivity of RD-CT was 63.6%, specificity was 74.1%, and accuracy was 72%. Among 62 patients with a CT-Respiration Ratio of less than 0.65, the sensitivity of RD-CT was 77.8%, specificity was 86.8%, and accuracy was 85.5%.

Conclusions

RD-CT may be clinically useful for detecting the presence of pleural adhesions. It can be adopted as one of the criteria for deciding the surgical approach.

     

NMRProcFlow: a graphical and interactive tool dedicated to 1D spectra processing for NMR-based metabolomics

Introduction

Concerning NMR-based metabolomics, 1D spectra processing often requires an expert eye for disentangling the intertwined peaks.

Objectives

The objective of NMRProcFlow is to assist the expert in this task in the best way without requirement of programming skills.

Methods

NMRProcFlow was developed to be a graphical and interactive 1D NMR (¹H & ¹³C) spectra processing tool.

Results

NMRProcFlow (http://nmrprocflow.org), dedicated to metabolic fingerprinting and targeted metabolomics, covers all spectra processing steps including baseline correction, chemical shift calibration and alignment.

Conclusion

Biologists and NMR spectroscopists can easily interact and develop synergies by visualizing the NMR spectra along with their corresponding experimental-factor levels, thus setting a bridge between experimental design and subsequent statistical analyses.

     

KniMet: a pipeline for the processing of chromatography–mass spectrometry metabolomics data

Introduction

Data processing is one of the biggest problems in metabolomics, given the high number of samples analyzed and the need of multiple software packages for each step of the processing workflow.

Objectives

Merge in the same platform the steps required for metabolomics data processing.

Methods

KniMet is a workflow for the processing of mass spectrometry-metabolomics data based on the KNIME Analytics platform.

Results

The approach includes key steps to follow in metabolomics data processing: feature filtering, missing value imputation, normalization, batch correction and annotation.

Conclusion

KniMet provides the user with a local, modular and customizable workflow for the processing of both GC–MS and LC–MS open profiling data.

     

A Boolean network model of human gonadal sex determination

Background

Gonadal sex determination (GSD) in humans is a complex biological process that takes place in early stages of embryonic development when the bipotential gonadal primordium (BGP) differentiates towards testes or ovaries. This decision is directed by one of two distinct pathways embedded in a GSD network activated in a population of coelomic epithelial cells, the Sertoli progenitor cells (SPC) and the granulosa progenitor cells (GPC). In males, the pathway is activated when the Sex-Determining Region Y (SRY) gene starts to be expressed, whereas in females the WNT4/ β-catenin pathway promotes the differentiation of the GPCs towards ovaries. The interactions and dynamics of the elements that constitute the GSD network are poorly understood, thus our group is interested in inferring the general architecture of this network as well as modeling the dynamic behavior of a set of genes associated to this process under wild-type and mutant conditions.

Methods

We reconstructed the regulatory network of GSD with a set of genes directly associated with the process of differentiation from SPC and GPC towards Sertoli and granulosa cells, respectively. These genes are experimentally well-characterized and the effects of their deficiency have been clinically reported. We modeled this GSD network as a synchronous Boolean network model (BNM) and characterized its attractors under wild-type and mutant conditions.

Results

Three attractors with a clear biological meaning were found; one of them corresponding to the currently known gene expression pattern of Sertoli cells, the second correlating to the granulosa cells and, the third resembling a disgenetic gonad.

Conclusions

The BNM of GSD that we present summarizes the experimental data on the pathways for Sertoli and granulosa establishment and sheds light on the overall behavior of a population of cells that differentiate within the developing gonad. With this model we propose a set of regulatory interactions needed to activate either the SRY or the WNT4/ β-catenin pathway as well as their downstream targets, which are critical for further sex differentiation. In addition, we observed a pattern of altered regulatory interactions and their dynamics that lead to some disorders of sex development (DSD).

     

Dynamical crises,multistability and the influence of the duration of immunity in a seasonally-forced model of disease transmission

Background

Highly successful strategies to make populations more resilient to infectious diseases, such as childhood vaccinations programs, may nonetheless lead to unpredictable outcomes due to the interplay between seasonal variations in transmission and a population’s immune status.

Methods

Motivated by the study of diseases such as pertussis we introduce a seasonally-forced susceptible-infectious-recovered model of disease transmission with waning and boosting of immunity. We study the system’s dynamical properties using a combination of numerical simulations and bifurcation techniques, paying particular attention to the properties of the initial condition space.

Results

We find that highly unpredictable behaviour can be triggered by changes in biologically relevant model parameters such as the duration of immunity. In the particular system we analyse — previously used in the literature to study pertussis dynamics — we identify the presence of an initial-condition landscape containing three coexisting attractors. The system’s response to interventions which perturb population immunity (e.g. vaccination "catch-up" campaigns) is therefore difficult to predict.

Conclusion

Given the increasing use of models to inform policy decisions regarding vaccine introduction and scheduling and infectious diseases intervention policy more generally, our findings highlight the importance of thoroughly investigating the dynamical properties of those models to identify key areas of uncertainty. Our findings suggest that the often stated tension between capturing biological complexity and utilising mathematically simple models is perhaps more nuanced than generally suggested. Simple dynamical models, particularly those which include forcing terms, can give rise to incredibly complex behaviour.