Click here to close
Hello! We notice that you are using Internet Explorer, which is not supported by Xenbase and may cause the site to display incorrectly.
We suggest using a current version of Chrome,
FireFox, or Safari.
BMC Bioinformatics
2022 Nov 24;231:503. doi: 10.1186/s12859-022-05055-5.
Show Gene links
Show Anatomy links
Reconstructing gene regulatory networks of biological function using differential equations of multilayer perceptrons.
Mao G
,
Zeng R
,
Peng J
,
Zuo K
,
Pang Z
,
Liu J
.
???displayArticle.abstract???
BACKGROUND: Building biological networks with a certain function is a challenge in systems biology. For the functionality of small (less than ten nodes) biological networks, most methods are implemented by exhausting all possible network topological spaces. This exhaustive approach is difficult to scale to large-scale biological networks. And regulatory relationships are complex and often nonlinear or non-monotonic, which makes inference using linear models challenging.
RESULTS: In this paper, we propose a multi-layer perceptron-based differential equation method, which operates by training a fully connected neural network (NN) to simulate the transcription rate of genes in traditional differential equations. We verify whether the regulatory network constructed by the NN method can continue to achieve the expected biological function by verifying the degree of overlap between the regulatory network discovered by NN and the regulatory network constructed by the Hill function. And we validate our approach by adapting to noise signals, regulator knockout, and constructing large-scale gene regulatory networks using link-knockout techniques. We apply a real dataset (the mesoderm inducer Xenopus Brachyury expression) to construct the core topology of the gene regulatory network and find that Xbra is only strongly expressed at moderate levels of activin signaling.
CONCLUSION: We have demonstrated from the results that this method has the ability to identify the underlying network topology and functional mechanisms, and can also be applied to larger and more complex gene network topologies.
Fig. 1
Fully linked neural network. a Schematic diagram of the fully connected neural network training synthetic term f. b Three genes are used as an example to illustrate. The synthetic term π2
(π3
) of π2
(π3
) is evaluated by a fully connected network. π2
and π3
(wheat and light blue) can depend on all three variables: π3
, π2
and the input signal π1
Fig. 2
Time evolution process under noisy conditions. (1) Under the stimulus after adding Gaussian white noise to the input signal (the red line of Input), without any constraints on g2, the time evolution curves g3 and g2 obtained after training the NN, and the expression level of g3 (blue line) is the same as The target time progress value (targetβs blue dotted line) basically matches. (2) Cross-section information obtained by training a NN under noisy conditions. Three panels show the dependence of π3
, π2
on I, π3
or π2
with the other two variables fixed. (3) The regulatory network obtained from (2)
Fig. 3
Simulate regulator knockout. (1) The perturbed f function can be iterated to simulate the effect of mutants in which specific regulatory chains are deleted. For example, deletion of π3
βs modulating effect on π3
leads to an increase in π3
(from darker to brighter solid green lines), indicating self-inhibition (shown in the third panel). The difference in π2
levels is not important here (dotted line). A similar argument applies to the other three panel. (2) Describes the sensitivity of network sequence and adapt to the error. Pane shows the knockout technology through links, step 1 to 4 of the evolution process of the network topology. The minimum incoherent feed forward motif appears naturally (topology #4), before the network has too few links to adapt
Fig. 4
The AUROC and AUPR of GENIE3, BiXGBoost, SIGNET, GNIPLP, PoLoBag and our methods on DREAM4 InSilico_Size100 five networks
Fig. 5
E. coli network including 1484 genes. Each bar represents the performance of one method in which the abscissas are the corresponding AUROC (right) and AUPR (left) values
Fig. 6
the activin/gsc/Xbra system. The Activin gene was activated by the input signal of morphogenetic gradient (Bcd), so it began to imitate its gradient mode. The Activin gene activated Xbra gene and opened the positive feedback of Xbra gene at a certain threshold. The Activin gene activates the Goosecold gene, and when the concentration of the two genes accumulates high enough, it forces the Xbra gene down. However, the concentration is highest only on the left side, when the concentration of Goosecold gene is low and its inhibitory effect is low, so that Xbra gene reaches a stable state
Aderhold,
Approximate Bayesian inference in semi-mechanistic models.
2017, Pubmed
Aderhold,
Approximate Bayesian inference in semi-mechanistic models.
2017,
Pubmed
Aubin-Frankowski,
Gene regulation inference from single-cell RNA-seq data with linear differential equations and velocity inference.
2020,
Pubmed
Brunton,
Discovering governing equations from data by sparse identification of nonlinear dynamical systems.
2016,
Pubmed
Ehsan Elahi,
A method for estimating Hill function-based dynamic models of gene regulatory networks.
2018,
Pubmed
Ferrell,
Perfect and Near-Perfect Adaptation in Cell Signaling.
2016,
Pubmed
Gama-Castro,
RegulonDB version 9.0: high-level integration of gene regulation, coexpression, motif clustering and beyond.
2016,
Pubmed
Gardner,
Construction of a genetic toggle switch in Escherichia coli.
2000,
Pubmed
Ghosh Roy,
PoLoBag: Polynomial Lasso Bagging for signed gene regulatory network inference from expression data.
2021,
Pubmed
Green,
Morphogen gradients, positional information, and Xenopus: interplay of theory and experiment.
2002,
Pubmed
,
Xenbase
Hsiao,
Inferring robust gene networks from expression data by a sensitivity-based incremental evolution method.
2012,
Pubmed
Hsiao,
Reverse engineering gene regulatory networks: coupling an optimization algorithm with a parameter identification technique.
2014,
Pubmed
Huynh-Thu,
Inferring regulatory networks from expression data using tree-based methods.
2010,
Pubmed
Jozefczuk,
Metabolomic and transcriptomic stress response of Escherichia coli.
2010,
Pubmed
Karlebach,
Modelling and analysis of gene regulatory networks.
2008,
Pubmed
Kimura,
Inference of S-system models of genetic networks using a cooperative coevolutionary algorithm.
2005,
Pubmed
Luo,
SIGNET: single-cell RNA-seq-based gene regulatory network prediction using multiple-layer perceptron bagging.
2022,
Pubmed
Ma,
Defining network topologies that can achieve biochemical adaptation.
2009,
Pubmed
Mandal,
Reverse engineering of gene regulatory networks based on S-systems and Bat algorithm.
2016,
Pubmed
Marbach,
Revealing strengths and weaknesses of methods for gene network inference.
2010,
Pubmed
Marbach,
Wisdom of crowds for robust gene network inference.
2012,
Pubmed
Matsumoto,
SCODE: an efficient regulatory network inference algorithm from single-cell RNA-Seq during differentiation.
2017,
Pubmed
Muzzey,
A systems-level analysis of perfect adaptation in yeast osmoregulation.
2009,
Pubmed
Nakayama,
Inference of S-system models of gene regulatory networks using immune algorithm.
2011,
Pubmed
Oates,
Causal network inference using biochemical kinetics.
2014,
Pubmed
Penfold,
CSI: a nonparametric Bayesian approach to network inference from multiple perturbed time series gene expression data.
2015,
Pubmed
Qiao,
Network Topologies That Can Achieve Dual Function of Adaptation and Noise Attenuation.
2019,
Pubmed
Ren,
Finding Robust Adaptation Gene Regulatory Networks Using Multi-Objective Genetic Algorithm.
2016,
Pubmed
Shen,
Finding gene network topologies for given biological function with recurrent neural network.
2021,
Pubmed
Wahde,
Coarse-grained reverse engineering of genetic regulatory networks.
2000,
Pubmed
Wu,
Computational optimization for S-type biological systems: cockroach genetic algorithm.
2013,
Pubmed
Zhang,
Inference of gene regulatory networks using pseudo-time series data.
2021,
Pubmed
Zheng,
BiXGBoost: a scalable, flexible boosting-based method for reconstructing gene regulatory networks.
2019,
Pubmed
de Jong,
Modeling and simulation of genetic regulatory systems: a literature review.
2002,
Pubmed