1
|
Liu R, Zhu L. Specification testing for ordinary differential equation models with fixed design and applications to COVID-19 epidemic models. Comput Stat Data Anal 2023; 180:107616. [PMID: 36128441 PMCID: PMC9479380 DOI: 10.1016/j.csda.2022.107616] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Revised: 09/05/2022] [Accepted: 09/05/2022] [Indexed: 01/25/2023]
Abstract
Checking the models about the ongoing Coronavirus Disease 2019 (COVID-19) pandemic is an important issue. Some famous ordinary differential equation (ODE) models, such as the SIR and SEIR models have been used to describe and predict the epidemic trend. Still, in many cases, only part of the equations can be observed. A test is suggested to check possibly partially observed ODE models with a fixed design sampling scheme. The asymptotic properties of the test under the null, global and local alternative hypotheses are presented. Two new propositions about U-statistics with varying kernels based on independent but non-identical data are derived as essential tools. Some simulation studies are conducted to examine the performances of the test. Based on the available public data, it is found that the SEIR model, for modeling the data of COVID-19 infective cases in certain periods in Japan and Algeria, respectively, maybe not be appropriate by applying the proposed test.
Collapse
Affiliation(s)
- Ran Liu
- School of Mathematics and Statistics, Beijing Jiaotong University, Beijing, China
- Department of Mathematics, Hong Kong Baptist University, Hong Kong, China
| | - Lixing Zhu
- Advanced Institute of Natural Sciences, Beijing Normal University at Zhuhai, Zhuhai, Guangdong, China
- Department of Mathematics, Hong Kong Baptist University, Hong Kong, China
| |
Collapse
|
2
|
Wang Y, Ghosh SK. Nonparametric estimation of isotropic covariance function. J Nonparametr Stat 2022. [DOI: 10.1080/10485252.2022.2146111] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
Affiliation(s)
- Yiming Wang
- Department of Statistics, North Carolina State University, Raleigh, NC, USA
| | - Sujit K. Ghosh
- Department of Statistics, North Carolina State University, Raleigh, NC, USA
| |
Collapse
|
3
|
Bhaumik P, Shi W, Ghosal S. Two-step Bayesian methods for generalized regression driven by partial differential equations. BERNOULLI 2022. [DOI: 10.3150/21-bej1363] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Affiliation(s)
- Prithwish Bhaumik
- Department of Statistics, North Carolina State University, 2311 Stinson Drive, Raleigh, North Carolina 27695-8203, USA
| | - Wenli Shi
- Department of Statistics, North Carolina State University, 2311 Stinson Drive, Raleigh, North Carolina 27695-8203, USA
| | - Subhashis Ghosal
- Department of Statistics, North Carolina State University, 2311 Stinson Drive, Raleigh, North Carolina 27695-8203, USA
| |
Collapse
|
4
|
Bradley W, Boukouvala F. Two-Stage Approach to Parameter Estimation of Differential Equations Using Neural ODEs. Ind Eng Chem Res 2021. [DOI: 10.1021/acs.iecr.1c00552] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- William Bradley
- School of Chemical & Biomolecular Engineering, Georgia Institute of Technology, 311 Ferst Dr., N.W., Atlanta, Georgia 30332-0100, United States
| | - Fani Boukouvala
- School of Chemical & Biomolecular Engineering, Georgia Institute of Technology, 311 Ferst Dr., N.W., Atlanta, Georgia 30332-0100, United States
| |
Collapse
|
5
|
Regazzoni F, Chapelle D, Moireau P. Combining data assimilation and machine learning to build data-driven models for unknown long time dynamics-Applications in cardiovascular modeling. INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN BIOMEDICAL ENGINEERING 2021; 37:e3471. [PMID: 33913623 PMCID: PMC8365699 DOI: 10.1002/cnm.3471] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/07/2021] [Revised: 03/12/2021] [Accepted: 04/23/2021] [Indexed: 06/12/2023]
Abstract
We propose a method to discover differential equations describing the long-term dynamics of phenomena featuring a multiscale behavior in time, starting from measurements taken at the fast-scale. Our methodology is based on a synergetic combination of data assimilation (DA), used to estimate the parameters associated with the known fast-scale dynamics, and machine learning (ML), used to infer the laws underlying the slow-scale dynamics. Specifically, by exploiting the scale separation between the fast and the slow dynamics, we propose a decoupling of time scales that allows to drastically lower the computational burden. Then, we propose a ML algorithm that learns a parametric mathematical model from a collection of time series coming from the phenomenon to be modeled. Moreover, we study the interpretability of the data-driven models obtained within the black-box learning framework proposed in this paper. In particular, we show that every model can be rewritten in infinitely many different equivalent ways, thus making intrinsically ill-posed the problem of learning a parametric differential equation starting from time series. Hence, we propose a strategy that allows to select a unique representative model in each equivalence class, thus enhancing the interpretability of the results. We demonstrate the effectiveness and noise-robustness of the proposed methods through several test cases, in which we reconstruct several differential models starting from time series generated through the models themselves. Finally, we show the results obtained for a test case in the cardiovascular modeling context, which sheds light on a promising field of application of the proposed methods.
Collapse
Affiliation(s)
- Francesco Regazzoni
- MOX—Mathematics DepartmentPolitecnico di MilanoMilanoItaly
- M3DISIMInstitut National de Recherche en Informatique et en AutomatiquePalaiseauFrance
- LMSEcole Polytechnique, CNRS, Institut Polytechnique de ParisPalaiseauFrance
| | - Dominique Chapelle
- M3DISIMInstitut National de Recherche en Informatique et en AutomatiquePalaiseauFrance
- LMSEcole Polytechnique, CNRS, Institut Polytechnique de ParisPalaiseauFrance
| | - Philippe Moireau
- M3DISIMInstitut National de Recherche en Informatique et en AutomatiquePalaiseauFrance
- LMSEcole Polytechnique, CNRS, Institut Polytechnique de ParisPalaiseauFrance
| |
Collapse
|
6
|
Meng L, Zhang J, Zhang X, Feng G. Bayesian estimation of time-varying parameters in ordinary differential equation models with noisy time-varying covariates. COMMUN STAT-SIMUL C 2021. [DOI: 10.1080/03610918.2019.1565584] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Affiliation(s)
- Lixin Meng
- School of Statistics, Jilin University of Finance and Economics, Changchun, P R China
| | - Jiwei Zhang
- School of Mathematics and Statistics, Yunnan University, Kunming, P R China
| | - Xue Zhang
- China Institute of Rural Education Development, Northeast Normal University, Changchun, P R China
| | - Guozhong Feng
- School of Information Science and Technology, Northeast Normal University, Changchun, P R China
| |
Collapse
|
7
|
Clairon Q. A regularization method for the parameter estimation problem in ordinary differential equations via discrete optimal control theory. J Stat Plan Inference 2021. [DOI: 10.1016/j.jspi.2020.04.007] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
|
8
|
Hu X, Hu Y, Wu F, Leung RWT, Qin J. Integration of single-cell multi-omics for gene regulatory network inference. Comput Struct Biotechnol J 2020; 18:1925-1938. [PMID: 32774787 PMCID: PMC7385034 DOI: 10.1016/j.csbj.2020.06.033] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2020] [Revised: 06/17/2020] [Accepted: 06/20/2020] [Indexed: 12/20/2022] Open
Abstract
The advancement of single-cell sequencing technology in recent years has provided an opportunity to reconstruct gene regulatory networks (GRNs) with the data from thousands of single cells in one sample. This uncovers regulatory interactions in cells and speeds up the discoveries of regulatory mechanisms in diseases and biological processes. Therefore, more methods have been proposed to reconstruct GRNs using single-cell sequencing data. In this review, we introduce technologies for sequencing single-cell genome, transcriptome, and epigenome. At the same time, we present an overview of current GRN reconstruction strategies utilizing different single-cell sequencing data. Bioinformatics tools were grouped by their input data type and mathematical principles for reader's convenience, and the fundamental mathematics inherent in each group will be discussed. Furthermore, the adaptabilities and limitations of these different methods will also be summarized and compared, with the hope to facilitate researchers recognizing the most suitable tools for them.
Collapse
Affiliation(s)
- Xinlin Hu
- Shenzhen Key Laboratory of Advanced Machine Learning and Applications, College of Mathematics and Statistics, Shenzhen University, Shenzhen 518060, China
| | - Yaohua Hu
- Shenzhen Key Laboratory of Advanced Machine Learning and Applications, College of Mathematics and Statistics, Shenzhen University, Shenzhen 518060, China
| | - Fanjie Wu
- School of Pharmaceutical Sciences (Shenzhen), Sun Yat-sen University, Shenzhen 518107, China
| | - Ricky Wai Tak Leung
- School of Pharmaceutical Sciences (Shenzhen), Sun Yat-sen University, Shenzhen 518107, China
| | - Jing Qin
- School of Pharmaceutical Sciences (Shenzhen), Sun Yat-sen University, Shenzhen 518107, China
| |
Collapse
|
9
|
Ranciati S, Wit EC, Viroli C. Bayesian smooth‐and‐match inference for ordinary differential equations models linear in the parameters. STAT NEERL 2020. [DOI: 10.1111/stan.12192] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Saverio Ranciati
- Department of Statistical SciencesUniversity of Bologna Bologna Italy
| | - Ernst C. Wit
- Institute of Computational ScienceUniversità della Svizzera Italiana Lugano Switzerland
| | - Cinzia Viroli
- Department of Statistical SciencesUniversity of Bologna Bologna Italy
| |
Collapse
|
10
|
Zhang T, Sun Y, Li H, Yan G, Tanabe S, Miao R, Wang Y, Caffo BS, Quigg MS. Bayesian inference of a directional brain network model for intracranial EEG data. Comput Stat Data Anal 2020. [DOI: 10.1016/j.csda.2019.106847] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
11
|
Li H, Wang Y, Tanabe S, Sun Y, Yan G, Quigg MS, Zhang T. Mapping epileptic directional brain networks using intracranial EEG data. Biostatistics 2019; 22:613-628. [PMID: 31879751 DOI: 10.1093/biostatistics/kxz056] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2019] [Revised: 11/26/2019] [Accepted: 11/29/2019] [Indexed: 11/13/2022] Open
Abstract
The human brain is a directional network system, in which brain regions are network nodes and the influence exerted by one region on another is a network edge. We refer to this directional information flow from one region to another as directional connectivity. Seizures arise from an epileptic directional network; abnormal neuronal activities start from a seizure onset zone and propagate via a network to otherwise healthy brain regions. As such, effective epilepsy diagnosis and treatment require accurate identification of directional connections among regions, i.e., mapping of epileptic patients' brain networks. This article aims to understand the epileptic brain network using intracranial electroencephalographic data-recordings of epileptic patients' brain activities in many regions. The most popular models for directional connectivity use ordinary differential equations (ODE). However, ODE models are sensitive to data noise and computationally costly. To address these issues, we propose a high-dimensional state-space multivariate autoregression (SSMAR) model for the brain's directional connectivity. Different from standard multivariate autoregression and SSMAR models, the proposed SSMAR features a cluster structure, where the brain network consists of several clusters of densely connected brain regions. We develop an expectation-maximization algorithm to estimate the proposed model and use it to map the interregional networks of epileptic patients in different seizure stages. Our method reveals the evolution of brain networks during seizure development.
Collapse
Affiliation(s)
- Huazhang Li
- Department of Statistics, University of Virginia 148 Amphitheater Way, Charlottesville, VA 22904-4135, USA
| | - Yaotian Wang
- Department of Statistics, University of Virginia 148 Amphitheater Way, Charlottesville, VA 22904-4135, USA
| | - Seiji Tanabe
- Department of Statistics, University of Virginia 148 Amphitheater Way, Charlottesville, VA 22904-4135, USA
| | - Yinge Sun
- Department of Statistics, University of Virginia 148 Amphitheater Way, Charlottesville, VA 22904-4135, USA
| | - Guofen Yan
- Department of Statistics, University of Virginia 148 Amphitheater Way, Charlottesville, VA 22904-4135, USA
| | - Mark S Quigg
- Department of Statistics, University of Virginia 148 Amphitheater Way, Charlottesville, VA 22904-4135, USA
| | - Tingting Zhang
- Department of Statistics, University of Virginia 148 Amphitheater Way, Charlottesville, VA 22904-4135, USA
| |
Collapse
|
12
|
Dattner I, Huppert A. Modern statistical tools for inference and prediction of infectious diseases using mathematical models. Stat Methods Med Res 2019; 27:1927-1929. [PMID: 29846149 DOI: 10.1177/0962280217746456] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
Affiliation(s)
- Itai Dattner
- 1 Department of Statistics, University of Haifa, Haifa, Israel
| | - Amit Huppert
- 2 The Biostatistics & BIomathematics Unit, The Gertner Institute for Epidemiology and Health Policy Research, Chaim Sheba Medical Center, Tel HaShomer, Israel.,3 School of Public Health, the Sackler Faculty of Medicine, Tel-Aviv University, Tel Aviv, Israel
| |
Collapse
|
13
|
Torres-Cerna CE, Morales JA, Hernandez-Vargas EA. Modeling Quorum Sensing Dynamics and Interference on Escherichia coli. Front Microbiol 2019; 10:1835. [PMID: 31481938 PMCID: PMC6710385 DOI: 10.3389/fmicb.2019.01835] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/20/2019] [Accepted: 07/25/2019] [Indexed: 01/16/2023] Open
Abstract
Bacteria control the expression of specific genes by Quorum Sensing (QS). This works using small signaling molecules called Autoinducers (AIs), for example, the Autoinducer-2 (AI-2). In this work, we present a mathematical model that represents the AI-2 dynamics on Escherichia coli, which is linked to the cell growth and the lsr operon expression. The model is adjusted using experimental data. Our results suggest that the extracellular AI-2 activity level depends on the cell growth rate, and this activity depends on the cell exponential growth phase. The model was adapted to simulate the interference of QS mechanisms in a co-culture of two E. coli strains: a wild type strain and a knock out strain that detects AI-2 but does not produce it. Co-culture simulations unveiled two conditions to avoid the QS on the wild strain: when the knock out takes control of the growth medium and overcomes the wild strain, or when is pre-cultured to its mid-exponential phase and then added to the wild strain culture. Model simulations unveiled new insights about the interference of bacterial communication and offer new tools for QS control.
Collapse
Affiliation(s)
| | - J Alejandro Morales
- Computer Science Department, Universidad de Guadalajara, Guadalajara, Mexico
| | | |
Collapse
|
14
|
Tian X, Pang W, Wang Y, Guo K, Zhou Y. LatinPSO: An algorithm for simultaneously inferring structure and parameters of ordinary differential equations models. Biosystems 2019; 182:8-16. [PMID: 31167112 DOI: 10.1016/j.biosystems.2019.05.006] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2018] [Revised: 05/01/2019] [Accepted: 05/14/2019] [Indexed: 10/26/2022]
Abstract
Simultaneously inferring both the structure and parameters of Ordinary Differential Equations (ODEs) for a complex dynamic system is more practical in many systems identification problems, but it remains challenging due to the complexity of the underlying search space. In this research, we propose a novel algorithm based on Particle Swarm Optimization (PSO) and Latin Hypercube Sampling (LHS) to address the above problem. The proposed algorithm is termed LatinPSO, and it can be effectively used for inferring the structure and parameters of ODE models through time course data. To start with, the real Human Immunodeficiency Virus (HIV) model and several synthetic models are used for evaluating the performance of LatinPSO. Experimental results demonstrated that LatinPSO could find satisfactory candidate ODE models with appropriate structure and parameters.
Collapse
Affiliation(s)
- Xinliang Tian
- College of Computer Science and Technology, Jilin University, Changchun 130012, China; Key Laboratory of Symbolic Computation and Knowledge Engineering of Ministry of Education, China
| | - Wei Pang
- School of Natural and Computing Sciences, University of Aberdeen, AB24 3UE, UK
| | - Yizhang Wang
- College of Computer Science and Technology, Jilin University, Changchun 130012, China; Key Laboratory of Symbolic Computation and Knowledge Engineering of Ministry of Education, China
| | - Kaimin Guo
- College of Computer Science and Technology, Jilin University, Changchun 130012, China; Key Laboratory of Symbolic Computation and Knowledge Engineering of Ministry of Education, China
| | - You Zhou
- College of Computer Science and Technology, Jilin University, Changchun 130012, China; Key Laboratory of Symbolic Computation and Knowledge Engineering of Ministry of Education, China.
| |
Collapse
|
15
|
Nardini JT, Bortz DM. The influence of numerical error on parameter estimation and uncertainty quantification for advective PDE models. INVERSE PROBLEMS 2019; 35:065003. [PMID: 34121793 PMCID: PMC8191598 DOI: 10.1088/1361-6420/ab10bb] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
Advective partial differential equations can be used to describe many scientific processes. Two significant sources of error that can cause difficulties in inferring parameters from experimental data on these processes include (i) noise from the measurement and collection of experimental data and (ii) numerical error in approximating the forward solution to the advection equation. How this second source of error alters parameter estimation and uncertainty quantification during an inverse problem methodology is not well understood. As a step towards a better understanding of this problem, we present both analytical and computational results concerning how a least squares cost function and parameter estimator behave in the presence of numerical error in approximating solutions to the underlying advection equation. We investigate residual patterns to derive an autocorrelative statistical model that can improve parameter estimation and confidence interval computation for first order methods. Building on our results and their general nature, we provide guidelines for practitioners to determine when numerical or experimental error is the main source of error in their inference, along with suggestions of how to efficiently improve their results.
Collapse
Affiliation(s)
- John T Nardini
- Statistical and Applied Mathematical Sciences Institute, 4501 Research Commons, Suite 300 79 T.W. Alexander Drive, PO Box 110207 Durham, NC 27709, United States of America
- Department of Mathematics, North Carolina State University, Box 8205, Raleigh, NC 27695, United States of America
| | - D M Bortz
- Department of Applied Mathematics, University of Colorado, 526 UCB, Boulder, CO 80309-0526, United States of America
| |
Collapse
|
16
|
Cao X, Sandstede B, Luo X. A Functional Data Method for Causal Dynamic Network Modeling of Task-Related fMRI. Front Neurosci 2019; 13:127. [PMID: 30872989 PMCID: PMC6402339 DOI: 10.3389/fnins.2019.00127] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2018] [Accepted: 02/05/2019] [Indexed: 01/15/2023] Open
Abstract
Functional MRI (fMRI) is a popular approach to investigate brain connections and activations when human subjects perform tasks. Because fMRI measures the indirect and convoluted signals of brain activities at a lower temporal resolution, complex differential equation modeling methods (e.g., Dynamic Causal Modeling) are usually employed to infer the neuronal processes and to fit the resulting fMRI signals. However, this modeling strategy is computationally expensive and remains to be mostly a confirmatory or hypothesis-driven approach. One major statistical challenge here is to infer, in a data-driven fashion, the underlying differential equation models from fMRI data. In this paper, we propose a causal dynamic network (CDN) method to estimate brain activations and connections simultaneously. Our method links the observed fMRI data with the latent neuronal states modeled by an ordinary differential equation (ODE) model. Using the basis function expansion approach in functional data analysis, we develop an optimization-based criterion that combines data-fitting errors and ODE fitting errors. We also develop and implement a block coordinate-descent algorithm to compute the ODE parameters efficiently. We illustrate the numerical advantages of our approach using data from realistic simulations and two task-related fMRI experiments. Compared with various effective connectivity methods, our method achieves higher estimation accuracy while improving the computational speed by from tens to thousands of times. Though our method is developed for task-related fMRI, we also demonstrate the potential applicability of our method (with a simple modification) to resting-state fMRI, by analyzing both simulated and real data from medium-sized networks.
Collapse
Affiliation(s)
- Xuefei Cao
- Division of Applied Mathematics, Brown University, Providence, RI, United States
| | - Björn Sandstede
- Division of Applied Mathematics, Brown University, Providence, RI, United States
| | - Xi Luo
- Department of Biostatistics and Data Science, School of Public Health, The University of Texas Health Science Center at Houston, Houston, TX, United States
| |
Collapse
|
17
|
Wu L, Qiu X, Yuan YX, Wu H. Parameter Estimation and Variable Selection for Big Systems of Linear Ordinary Differential Equations: A Matrix-Based Approach. J Am Stat Assoc 2019; 114:657-667. [PMID: 34385718 DOI: 10.1080/01621459.2017.1423074] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Ordinary differential equations (ODEs) are widely used to model the dynamic behavior of a complex system. Parameter estimation and variable selection for a "Big System" with linear ODEs are very challenging due to the need of nonlinear optimization in an ultra-high dimensional parameter space. In this article, we develop a parameter estimation and variable selection method based on the ideas of similarity transformation and separable least squares (SLS). Simulation studies demonstrate that the proposed matrix-based SLS method could be used to estimate the coefficient matrix more accurately and perform variable selection for a linear ODE system with thousands of dimensions and millions of parameters much better than the direct least squares (LS) method and the vector-based two-stage method that are currently available. We applied this new method to two real data sets: a yeast cell cycle gene expression data set with 30 dimensions and 930 unknown parameters and the Standard & Poor 1500 index stock price data with 1250 dimensions and 1,563,750 unknown parameters, to illustrate the utility and numerical performance of the proposed parameter estimation and variable selection method for big systems in practice.
Collapse
Affiliation(s)
- Leqin Wu
- Department of Mathematics, Jinan University, Guangzhou, China
| | - Xing Qiu
- Department of Biostatistics and Computational Biology University of Rochester, Rochester, New York, U.S.A
| | - Ya-Xiang Yuan
- Academy of Mathematics and System Sciences Chinese Academy of Sciences, Beijing, China
| | - Hulin Wu
- Department of Biostatistics, University of Texas Health Science Center at Houston, Houston, TX, U.S.A
| |
Collapse
|
18
|
Xue H, Kumar A, Wu H. Parameter Estimation for Semiparametric Ordinary Differential Equation Models. COMMUN STAT-THEOR M 2019; 48:5985-6004. [PMID: 32952273 DOI: 10.1080/03610926.2018.1523433] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
Abstract
We propose a new class of two-stage parameter estimation methods for semiparametric ordinary differential equation (ODE) models. In the first stage, state variables are estimated using a penalized spline approach; In the second stage, form of numerical discretization algorithms for an ODE solver is used to formulate estimating equations. Estimated state variables from the first stage are used to obtain more data points for the second stage. Asymptotic properties for the proposed estimators are established. Simulation studies show that the method performs well, especially for small sample. Real life use of the method is illustrated using Influenza specific cell-trafficking study.
Collapse
Affiliation(s)
- Hongqi Xue
- iCardiac Technologies, 150 Allens Creek Rd, Rochester, NY 14618
| | - Arun Kumar
- Livanova, 100 Cyberonics Blvd, Houston, TX 77058
| | - Hulin Wu
- The university of texas health science center of Houston, 7000 Fannin St 1200, Houston, TX 77030
| |
Collapse
|
19
|
Lee K, Lee J, Dass SC. Inference for differential equation models using relaxation via dynamical systems. Comput Stat Data Anal 2018. [DOI: 10.1016/j.csda.2018.05.014] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/14/2022]
|
20
|
Xue H, Wu S, Wu Y, Idarraga JCR, Wu H. Independence screening for high dimensional nonlinear additive ODE models with applications to dynamic gene regulatory networks. Stat Med 2018; 37:2630-2644. [PMID: 29722041 PMCID: PMC6940146 DOI: 10.1002/sim.7669] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/08/2016] [Revised: 01/18/2018] [Accepted: 03/08/2018] [Indexed: 11/12/2022]
Abstract
Mechanism-driven low-dimensional ordinary differential equation (ODE) models are often used to model viral dynamics at cellular levels and epidemics of infectious diseases. However, low-dimensional mechanism-based ODE models are limited for modeling infectious diseases at molecular levels such as transcriptomic or proteomic levels, which is critical to understand pathogenesis of diseases. Although linear ODE models have been proposed for gene regulatory networks (GRNs), nonlinear regulations are common in GRNs. The reconstruction of large-scale nonlinear networks from time-course gene expression data remains an unresolved issue. Here, we use high-dimensional nonlinear additive ODEs to model GRNs and propose a 4-step procedure to efficiently perform variable selection for nonlinear ODEs. To tackle the challenge of high dimensionality, we couple the 2-stage smoothing-based estimation method for ODEs and a nonlinear independence screening method to perform variable selection for the nonlinear ODE models. We have shown that our method possesses the sure screening property and it can handle problems with non-polynomial dimensionality. Numerical performance of the proposed method is illustrated with simulated data and a real data example for identifying the dynamic GRN of Saccharomyces cerevisiae.
Collapse
Affiliation(s)
- Hongqi Xue
- iCardiac Technologies, 150 Allens Creek Road, Rochester, NY 14618, USA
| | - Shuang Wu
- Biogen, 300 Binney Street, Cambridge, MA 02142, USA
| | - Yichao Wu
- Department of Mathematics, Statistics and Computer Science, University of Illinois at Chicago, Chicago, IL 60607-7045, USA
| | | | - Hulin Wu
- Department of Biostatistics and Data Science, School of Public Health, University of Texas Health Science Center at Houston, 1200 Pressler Street, RAS E833, Houston, TX 77030, USA
| |
Collapse
|
21
|
Clairon Q, Brunel NJB. Optimal Control and Additive Perturbations Help in Estimating Ill-Posed and Uncertain Dynamical Systems. J Am Stat Assoc 2018. [DOI: 10.1080/01621459.2017.1319841] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Affiliation(s)
- Quentin Clairon
- School of Mathematics and Statistics, University of Newcastle, UK
| | - Nicolas J.-B. Brunel
- ENSIIE and Laboratoire de Mathématiques et Modélisation d’Evry, Université d’Evry Val d’Essonne, UMR CNRS, France
| |
Collapse
|
22
|
Niu M, Macdonald B, Rogers S, Filippone M, Husmeier D. Statistical inference in mechanistic models: time warping for improved gradient matching. Comput Stat 2018; 33:1091-1123. [PMID: 31258254 PMCID: PMC6560940 DOI: 10.1007/s00180-017-0753-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2016] [Accepted: 07/19/2017] [Indexed: 10/31/2022]
Abstract
Inference in mechanistic models of non-linear differential equations is a challenging problem in current computational statistics. Due to the high computational costs of numerically solving the differential equations in every step of an iterative parameter adaptation scheme, approximate methods based on gradient matching have become popular. However, these methods critically depend on the smoothing scheme for function interpolation. The present article adapts an idea from manifold learning and demonstrates that a time warping approach aiming to homogenize intrinsic length scales can lead to a significant improvement in parameter estimation accuracy. We demonstrate the effectiveness of this scheme on noisy data from two dynamical systems with periodic limit cycle, a biopathway, and an application from soft-tissue mechanics. Our study also provides a comparative evaluation on a wide range of signal-to-noise ratios.
Collapse
Affiliation(s)
- Mu Niu
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Benn Macdonald
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| | - Simon Rogers
- Department of Computer Science, University of Glasgow, Glasgow, UK
| | | | - Dirk Husmeier
- School of Mathematics and Statistics, University of Glasgow, Glasgow, UK
| |
Collapse
|
23
|
Bhaumik P, Ghosal S. Efficient Bayesian estimation and uncertainty quantification in ordinary differential equation models. BERNOULLI 2017. [DOI: 10.3150/16-bej856] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
24
|
Zhang T, Yin Q, Caffo B, Sun Y, Boatman-Reich D. Bayesian inference of high-dimensional, cluster-structured ordinary differential equation models with applications to brain connectivity studies. Ann Appl Stat 2017. [DOI: 10.1214/17-aoas1021] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
25
|
|
26
|
Zhang X, Cao J, Carroll RJ. Estimating varying coefficients for partial differential equation models. Biometrics 2017; 73:949-959. [PMID: 28076654 DOI: 10.1111/biom.12646] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2015] [Revised: 11/01/2016] [Accepted: 12/01/2016] [Indexed: 11/29/2022]
Abstract
Partial differential equations (PDEs) are used to model complex dynamical systems in multiple dimensions, and their parameters often have important scientific interpretations. In some applications, PDE parameters are not constant but can change depending on the values of covariates, a feature that we call varying coefficients. We propose a parameter cascading method to estimate varying coefficients in PDE models from noisy data. Our estimates of the varying coefficients are shown to be consistent and asymptotically normally distributed. The performance of our method is evaluated by a simulation study and by an empirical study estimating three varying coefficients in a PDE model arising from LIDAR data.
Collapse
Affiliation(s)
- Xinyu Zhang
- Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China and Statistics and Mathematics College, Yunnan University of Finance and Economics, Kunming 650221, China
| | - Jiguo Cao
- Department of Statistics and Actuarial Science, Simon Fraser University, Burnaby, British Columbia, V5A1S6, Canada
| | - Raymond J Carroll
- Department of Statistics, Texas A&M University, 3143 TAMU, College Station, Texas 77843-3143, U.S.A.,School of Mathematical and Physical Sciences, University of Technology, Sydney, PO Box 123, Broadway, New South Wales 2007, Australia
| |
Collapse
|
27
|
Chen S, Shojaie A, Witten DM. Network Reconstruction From High-Dimensional Ordinary Differential Equations. J Am Stat Assoc 2017; 112:1697-1707. [PMID: 29618851 DOI: 10.1080/01621459.2016.1229197] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
We consider the task of learning a dynamical system from high-dimensional time-course data. For instance, we might wish to estimate a gene regulatory network from gene expression data measured at discrete time points. We model the dynamical system nonparametrically as a system of additive ordinary differential equations. Most existing methods for parameter estimation in ordinary differential equations estimate the derivatives from noisy observations. This is known to be challenging and inefficient. We propose a novel approach that does not involve derivative estimation. We show that the proposed method can consistently recover the true network structure even in high dimensions, and we demonstrate empirical improvement over competing approaches. Supplementary materials for this article are available online.
Collapse
Affiliation(s)
- Shizhe Chen
- Department of Biostatistics, University of Washington, WA
| | - Ali Shojaie
- Departments of Biostatistics and Statistics, University of Washington, WA
| | - Daniela M Witten
- Departments of Biostatistics and Statistics, University of Washington, WA
| |
Collapse
|
28
|
Paul D, Peng J, Burman P. Nonparametric estimation of dynamics of monotone trajectories. Ann Stat 2016. [DOI: 10.1214/15-aos1409] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
29
|
Duvigneau S, Sharma-Chawla N, Boianelli A, Stegemann-Koniszewski S, Nguyen VK, Bruder D, Hernandez-Vargas EA. Hierarchical effects of pro-inflammatory cytokines on the post-influenza susceptibility to pneumococcal coinfection. Sci Rep 2016; 6:37045. [PMID: 27872472 PMCID: PMC5181841 DOI: 10.1038/srep37045] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2016] [Accepted: 10/20/2016] [Indexed: 02/07/2023] Open
Abstract
In the course of influenza A virus (IAV) infections, a secondary bacterial infection frequently leads to serious respiratory conditions provoking high hospitalization and death tolls. Although abundant pro-inflammatory responses have been reported as key contributing factors for these severe dual infections, the relative contributions of cytokines remain largely unclear. In the current study, mathematical modelling based on murine experimental data dissects IFN-γ as a cytokine candidate responsible for impaired bacterial clearance, thereby promoting bacterial growth and systemic dissemination during acute IAV infection. We also found a time-dependent detrimental role of IL-6 in curtailing bacterial outgrowth which was not as distinct as for IFN-γ. Our numerical simulations suggested a detrimental effect of IFN-γ alone and in synergism with IL-6 but no conclusive pathogenic effect of IL-6 and TNF-α alone. This work provides a rationale to understand the potential impact of how to manipulate temporal immune components, facilitating the formulation of hypotheses about potential therapeutic strategies to treat coinfections.
Collapse
Affiliation(s)
- Stefanie Duvigneau
- Infection Immunology Group, Institute of Medical Microbiology, Disease Prevention and Control, Otto-von-Guericke University Magdeburg, Germany.,Systems Medicine of Infectious Disease Group, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig, Germany.,Immune Regulation Group, Helmholtz Centre for Infection Research, Braunschweig, Germany.,Chair for Automation/Modeling, Institute for Automation Engineering, Otto-von-Guericke University Magdeburg, Germany
| | - Niharika Sharma-Chawla
- Infection Immunology Group, Institute of Medical Microbiology, Disease Prevention and Control, Otto-von-Guericke University Magdeburg, Germany.,Immune Regulation Group, Helmholtz Centre for Infection Research, Braunschweig, Germany
| | - Alessandro Boianelli
- Systems Medicine of Infectious Disease Group, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig, Germany
| | - Sabine Stegemann-Koniszewski
- Infection Immunology Group, Institute of Medical Microbiology, Disease Prevention and Control, Otto-von-Guericke University Magdeburg, Germany.,Immune Regulation Group, Helmholtz Centre for Infection Research, Braunschweig, Germany
| | - Van Kinh Nguyen
- Systems Medicine of Infectious Disease Group, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig, Germany
| | - Dunja Bruder
- Infection Immunology Group, Institute of Medical Microbiology, Disease Prevention and Control, Otto-von-Guericke University Magdeburg, Germany.,Immune Regulation Group, Helmholtz Centre for Infection Research, Braunschweig, Germany
| | - Esteban A Hernandez-Vargas
- Systems Medicine of Infectious Disease Group, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig, Germany
| |
Collapse
|
30
|
Sun X, Hu F, Wu S, Qiu X, Linel P, Wu H. Controllability and stability analysis of large transcriptomic dynamic systems for host response to influenza infection in human. Infect Dis Model 2016; 1:52-70. [PMID: 29928721 PMCID: PMC5963324 DOI: 10.1016/j.idm.2016.07.002] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2016] [Accepted: 07/08/2016] [Indexed: 12/20/2022] Open
Abstract
Background Gene regulatory networks are complex dynamic systems and the reverse-engineering of such networks from high-dimensional time course transcriptomic data have attracted researchers from various fields. It is also interesting and important to study the behavior of the reconstructed networks on the basis of dynamic models and the biological mechanisms. We focus on the gene regulatory networks reconstructed using the ordinary differential equation (ODE) modelling approach and investigate the properties of these networks. Results Controllability and stability analyses are conducted for the reconstructed gene response networks of 17 influenza infected subjects based on ODE models. Symptomatic subjects tend to have larger numbers of driver nodes, higher proportions of critical links and lower proportions of redundant links than asymptomatic subjects. We also show that the degree distribution, rather than the structure of networks, plays an important role in controlling the network in response to influenza infection. In addition, we find that the stability of high-dimensional networks is very sensitive to randomness in the reconstructed systems brought by errors in measurements and parameter estimation. Conclusions The gene response networks of asymptomatic subjects are easier to be controlled than those of symptomatic subjects. This may indicate that the regulatory systems of asymptomatic subjects are easier to recover from disease stimulations, so these subjects are less likely to develop symptoms. Our results also suggest that stability constraint should be considered in the modelling of high-dimensional networks and the estimation of network parameters.
Collapse
Affiliation(s)
- Xiaodian Sun
- Biostatistics and Bioinformatics Core, Sylvester Comprehensive Cancer Center, University of Miami, Miami, USA
| | - Fang Hu
- Department of Biostatistics and Computational Biology, University of Rochester School of Medicine and Dentistry, Rochester, NY, USA
| | - Shuang Wu
- Genus PLC, ABS Global, Deforest, WI, USA
| | - Xing Qiu
- Department of Biostatistics and Computational Biology, University of Rochester School of Medicine and Dentistry, Rochester, NY, USA
| | | | - Hulin Wu
- Department of Biostatistics, School of Public Health, University of Texas Health Science Center at Houston, Houston, TX, USA
| |
Collapse
|
31
|
Conrad PR, Girolami M, Särkkä S, Stuart A, Zygalakis K. Statistical analysis of differential equations: introducing probability measures on numerical solutions. STATISTICS AND COMPUTING 2016; 27:1065-1082. [PMID: 32226237 PMCID: PMC7089645 DOI: 10.1007/s11222-016-9671-0] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/05/2015] [Accepted: 05/05/2016] [Indexed: 06/10/2023]
Abstract
In this paper, we present a formal quantification of uncertainty induced by numerical solutions of ordinary and partial differential equation models. Numerical solutions of differential equations contain inherent uncertainties due to the finite-dimensional approximation of an unknown and implicitly defined function. When statistically analysing models based on differential equations describing physical, or other naturally occurring, phenomena, it can be important to explicitly account for the uncertainty introduced by the numerical method. Doing so enables objective determination of this source of uncertainty, relative to other uncertainties, such as those caused by data contaminated with noise or model error induced by missing physical or inadequate descriptors. As ever larger scale mathematical models are being used in the sciences, often sacrificing complete resolution of the differential equation on the grids used, formally accounting for the uncertainty in the numerical method is becoming increasingly more important. This paper provides the formal means to incorporate this uncertainty in a statistical model and its subsequent analysis. We show that a wide variety of existing solvers can be randomised, inducing a probability measure over the solutions of such differential equations. These measures exhibit contraction to a Dirac measure around the true unknown solution, where the rates of convergence are consistent with the underlying deterministic numerical method. Furthermore, we employ the method of modified equations to demonstrate enhanced rates of convergence to stochastic perturbations of the original deterministic problem. Ordinary differential equations and elliptic partial differential equations are used to illustrate the approach to quantify uncertainty in both the statistical analysis of the forward and inverse problems.
Collapse
Affiliation(s)
| | - Mark Girolami
- Department of Statistics, University of Warwick, Coventry, UK
- Present Address: Alan Turing Institute, London, UK
| | - Simo Särkkä
- Department of Electrical Engineering and Automation, Aalto University, Espoo, Finland
| | - Andrew Stuart
- Department of Mathematics, University of Warwick, Coventry, UK
| | | |
Collapse
|
32
|
Vujačić I, Mahmoudi SM, Wit E. Generalized Tikhonov regularization in estimation of ordinary differential equations models. Stat (Int Stat Inst) 2016. [DOI: 10.1002/sta4.111] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Ivan Vujačić
- Department of Mathematics; VU University Amsterdam; Room S3.30, De Boelelaan 1081a Amsterdam 1081 HV The Netherlands
| | - Seyed Mahdi Mahmoudi
- Statistics and Probability; Johann Bernoulli Institute; Nijenborgh 9 Groningen 9747 AG The Netherlands
| | - Ernst Wit
- Statistics and Probability; Johann Bernoulli Institute; Nijenborgh 9 Groningen 9747 AG The Netherlands
| |
Collapse
|
33
|
Frasso G, Jaeger J, Lambert P. Inference in dynamic systems using B-splines and quasilinearized ODE penalties. Biom J 2015; 58:691-714. [PMID: 26602190 DOI: 10.1002/bimj.201500082] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2015] [Revised: 05/13/2015] [Accepted: 08/20/2015] [Indexed: 11/09/2022]
Abstract
Nonlinear (systems of) ordinary differential equations (ODEs) are common tools in the analysis of complex one-dimensional dynamic systems. We propose a smoothing approach regularized by a quasilinearized ODE-based penalty. Within the quasilinearized spline-based framework, the estimation reduces to a conditionally linear problem for the optimization of the spline coefficients. Furthermore, standard ODE compliance parameter(s) selection criteria are applicable. We evaluate the performances of the proposed strategy through simulated and real data examples. Simulation studies suggest that the proposed procedure ensures more accurate estimates than standard nonlinear least squares approaches when the state (initial and/or boundary) conditions are not known.
Collapse
Affiliation(s)
- Gianluca Frasso
- Faculté des Sciences Sociales, Méthodes Quantitatives en Sciences Sociales, Université de Liège, Boulevard du Rectorat 7, B-4000, Liège, Belgium
| | | | - Philippe Lambert
- Faculté des Sciences Sociales, Méthodes Quantitatives en Sciences Sociales, Université de Liège, Boulevard du Rectorat 7, B-4000, Liège, Belgium.,Institut de Statistique, Biostatistique et Sciences Actuarielles, Université Catholique de Louvain, Belgium
| |
Collapse
|
34
|
Boianelli A, Nguyen VK, Ebensen T, Schulze K, Wilk E, Sharma N, Stegemann-Koniszewski S, Bruder D, Toapanta FR, Guzmán CA, Meyer-Hermann M, Hernandez-Vargas EA. Modeling Influenza Virus Infection: A Roadmap for Influenza Research. Viruses 2015; 7:5274-304. [PMID: 26473911 PMCID: PMC4632383 DOI: 10.3390/v7102875] [Citation(s) in RCA: 80] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2015] [Revised: 09/28/2015] [Accepted: 09/28/2015] [Indexed: 12/24/2022] Open
Abstract
Influenza A virus (IAV) infection represents a global threat causing seasonal outbreaks and pandemics. Additionally, secondary bacterial infections, caused mainly by Streptococcus pneumoniae, are one of the main complications and responsible for the enhanced morbidity and mortality associated with IAV infections. In spite of the significant advances in our knowledge of IAV infections, holistic comprehension of the interplay between IAV and the host immune response (IR) remains largely fragmented. During the last decade, mathematical modeling has been instrumental to explain and quantify IAV dynamics. In this paper, we review not only the state of the art of mathematical models of IAV infection but also the methodologies exploited for parameter estimation. We focus on the adaptive IR control of IAV infection and the possible mechanisms that could promote a secondary bacterial coinfection. To exemplify IAV dynamics and identifiability issues, a mathematical model to explain the interactions between adaptive IR and IAV infection is considered. Furthermore, in this paper we propose a roadmap for future influenza research. The development of a mathematical modeling framework with a secondary bacterial coinfection, immunosenescence, host genetic factors and responsiveness to vaccination will be pivotal to advance IAV infection understanding and treatment optimization.
Collapse
Affiliation(s)
- Alessandro Boianelli
- Systems Medicine of Infectious Diseases, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
| | - Van Kinh Nguyen
- Systems Medicine of Infectious Diseases, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
| | - Thomas Ebensen
- Department of Vaccinology and Applied Microbiology, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
| | - Kai Schulze
- Department of Vaccinology and Applied Microbiology, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
| | - Esther Wilk
- Department of Infection Genetics, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
| | - Niharika Sharma
- Immune Regulation, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
| | | | - Dunja Bruder
- Immune Regulation, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
- Infection Immunology, Institute of Medical Microbiology, Infection Control and Prevention, Otto-von-Guericke-University, Magdeburg 39106, Germany.
| | - Franklin R Toapanta
- Center for Vaccine Development, University of Maryland School of Medicine, Baltimore, Maryland 21201, USA.
| | - Carlos A Guzmán
- Department of Vaccinology and Applied Microbiology, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
| | - Michael Meyer-Hermann
- Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
- Institute for Biochemistry, Biotechnology and Bioinformatics, Technische Universität Braunschweig, Braunschweig 38106, Germany.
| | - Esteban A Hernandez-Vargas
- Systems Medicine of Infectious Diseases, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research, Braunschweig 38124, Germany.
| |
Collapse
|
35
|
Qiu X, Wu S, Hilchey SP, Thakar J, Liu ZP, Welle SL, Henn AD, Wu H, Zand MS. Diversity in Compartmental Dynamics of Gene Regulatory Networks: The Immune Response in Primary Influenza A Infection in Mice. PLoS One 2015; 10:e0138110. [PMID: 26413862 PMCID: PMC4586376 DOI: 10.1371/journal.pone.0138110] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2015] [Accepted: 08/26/2015] [Indexed: 01/23/2023] Open
Abstract
Current approaches to study transcriptional profiles post influenza infection typically rely on tissue sampling from one or two sites at a few time points, such as spleen and lung in murine models. In this study, we infected female C57/BL6 mice intranasally with mouse-adapted H3N2/Hong Kong/X31 avian influenza A virus, and then analyzed the gene expression profiles in four different compartments (blood, lung, mediastinal lymph nodes, and spleen) over 11 consecutive days post infection. These data were analyzed by an advanced statistical procedure based on ordinary differential equation (ODE) modeling. Vastly different lists of significant genes were identified by the same statistical procedure in each compartment. Only 11 of them are significant in all four compartments. We classified significant genes in each compartment into co-expressed modules based on temporal expression patterns. We then performed functional enrichment analysis on these co-expression modules and identified significant pathway and functional motifs. Finally, we used an ODE based model to reconstruct gene regulatory network (GRN) for each compartment and studied their network properties.
Collapse
Affiliation(s)
- Xing Qiu
- Department of Biostatistics and Computational Biology, University of Rochester Medical Center, Rochester, NY, 14642, United States of America
| | - Shuang Wu
- Department of Biostatistics and Computational Biology, University of Rochester Medical Center, Rochester, NY, 14642, United States of America
| | - Shannon P. Hilchey
- Department of Medicine, University of Rochester Medical Center, Rochester, NY, 14642, United States of America
| | - Juilee Thakar
- Department of Biostatistics and Computational Biology, University of Rochester Medical Center, Rochester, NY, 14642, United States of America
- Department of Microbiology and Immunology, University of Rochester, Rochester, NY, 14642 United States of America
| | - Zhi-Ping Liu
- Department of Biostatistics and Computational Biology, University of Rochester Medical Center, Rochester, NY, 14642, United States of America
- Department of Biomedical Engineering, Shandong University, Jinan, Shandong, China
| | - Stephen L. Welle
- Functional Genomics Center, University of Rochester, Rochester, NY, 14642, United States of America
| | - Alicia D. Henn
- Department of Medicine, University of Rochester Medical Center, Rochester, NY, 14642, United States of America
- Department of Microbiology and Immunology, University of Rochester, Rochester, NY, 14642 United States of America
| | - Hulin Wu
- Department of Biostatistics, School of Public Health, University of Texas Health Science Center at Houston, Houston, TX, 77030, United States of America
- * E-mail: (HW); (MSZ)
| | - Martin S. Zand
- Department of Medicine, University of Rochester Medical Center, Rochester, NY, 14642, United States of America
- * E-mail: (HW); (MSZ)
| |
Collapse
|
36
|
Linel P, Wu S, Deng N, Wu H. Dynamic transcriptional signatures and network responses for clinical symptoms in influenza-infected human subjects using systems biology approaches. J Pharmacokinet Pharmacodyn 2015; 41:509-21. [PMID: 25015847 DOI: 10.1007/s10928-014-9365-1] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2014] [Accepted: 06/14/2014] [Indexed: 01/10/2023]
Abstract
Recent studies demonstrate that human blood transcriptional signatures may be used to support diagnosis and clinical decisions for acute respiratory viral infections such as influenza. In this article, we propose to use a newly developed systems biology approach for time course gene expression data to identify significant dynamically response genes and dynamic gene network responses to viral infection. We illustrate the methodological pipeline by reanalyzing the time course gene expression data from a study with healthy human subjects challenged by live influenza virus. We observed clear differences in the number of significant dynamic response genes (DRGs) between the symptomatic and asymptomatic subjects and also identified DRG signatures for symptomatic subjects with influenza infection. The 505 common DRGs shared by the symptomatic subjects have high consistency with the signature genes for predicting viral infection identified in previous works. The temporal response patterns and network response features were carefully analyzed and investigated.
Collapse
|
37
|
Dattner I. A model-based initial guess for estimating parameters in systems of ordinary differential equations. Biometrics 2015; 71:1176-84. [PMID: 26172865 DOI: 10.1111/biom.12348] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2014] [Revised: 04/01/2015] [Accepted: 05/01/2015] [Indexed: 11/28/2022]
Abstract
The inverse problem of parameter estimation from noisy observations is a major challenge in statistical inference for dynamical systems. Parameter estimation is usually carried out by optimizing some criterion function over the parameter space. Unless the optimization process starts with a good initial guess, the estimation may take an unreasonable amount of time, and may converge to local solutions, if at all. In this article, we introduce a novel technique for generating good initial guesses that can be used by any estimation method. We focus on the fairly general and often applied class of systems linear in the parameters. The new methodology bypasses numerical integration and can handle partially observed systems. We illustrate the performance of the method using simulations and apply it to real data.
Collapse
Affiliation(s)
- Itai Dattner
- Department of Statistics, University of Haifa, 199 Aba Khoushy Ave. Mount Carmel, Haifa, 3498838, Israel
| |
Collapse
|
38
|
Hu T, Qiu Y, Cui H. Robust estimation of constant and time-varying parameters in nonlinear ordinary differential equation models. J Nonparametr Stat 2015. [DOI: 10.1080/10485252.2015.1042377] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
39
|
Wu H, Miao H, Xue H, Topham DJ, Zand M. Quantifying Immune Response to Influenza Virus Infection via Multivariate Nonlinear ODE Models with Partially Observed State Variables and Time-Varying Parameters. STATISTICS IN BIOSCIENCES 2015; 7:147-166. [PMID: 26085850 PMCID: PMC4465846 DOI: 10.1007/s12561-014-9108-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2013] [Accepted: 01/12/2014] [Indexed: 01/22/2023]
Abstract
Influenza A virus (IAV) infection continues to be a global health threat, as evidenced by the outbreak of the novel A/California/7/2009 IAV strain. Previous flu vaccines have proven less effective than hoped for emerging IAV strains, indicating a more thorough understanding of immune responses to primary infection is needed. One issue is the difficulty in directly measuring many key parameters and variables of the immune response. To address these issues, we considered a comprehensive workflow for statistical inference for ordinary differential question (ODE) models with partially observed variables and time-varying parameters, including identifiability analysis, two-stage and NLS estimation, and model selection etc‥ In particular, we proposed a novel one-step method to verify parameter identifiability and formulate estimating equations simultaneously. Thus, the pseudo-LS method can now deal with general ODE models with partially observed state variables for the first time. Using this workflow, we verified the relative significance of various immune factors to virus control, including target epithelial cells, cytotoxic T-lymphocyte (CD8+) cells and IAV specific antibodies (IgG and IgM). Factors other than cytotoxic T-lymphocyte (CTL) killing contributed the most to the loss of infected epithelial cells, though the effects of CTL are still significant. IgM antibody was found to be the major contributor to neutralization of free infectious viral particles. Also, the maximum viral load, which correlates well with mortality, was found to depend more on viral replication rates than infectivity. In contrast to current hypotheses, the results obtained via our methods suggest that IgM antibody and viral replication rates may be worth of further explorations in vaccine development.
Collapse
Affiliation(s)
- Hulin Wu
- Department of Biostatistics and Computational Biology, University of Rochester School of Medicine and Dentistry, 601 Elmwood Avenue, Box 630, Rochester, New York 14642
| | - Hongyu Miao
- Department of Biostatistics and Computational Biology, University of Rochester School of Medicine and Dentistry, 601 Elmwood Avenue, Box 630, Rochester, New York 14642
| | - Hongqi Xue
- Department of Biostatistics and Computational Biology, University of Rochester School of Medicine and Dentistry, 601 Elmwood Avenue, Box 630, Rochester, New York 14642
| | - David J. Topham
- David H. Smith Center for Vaccine Biology & Immunology, Department of Microbiology and Immunology, University of Rochester School of Medicine and Dentistry, NY, 14642
| | - Martin Zand
- Department of Medicine, Division of Nephrology, University of Rochester School of Medicine and Dentistry, 601 Elmwood Avenue, Box 675, Rochester, New York 14642
| |
Collapse
|
40
|
Nguyen VK, Binder SC, Boianelli A, Meyer-Hermann M, Hernandez-Vargas EA. Ebola virus infection modeling and identifiability problems. Front Microbiol 2015; 6:257. [PMID: 25914675 PMCID: PMC4391033 DOI: 10.3389/fmicb.2015.00257] [Citation(s) in RCA: 42] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2014] [Accepted: 03/16/2015] [Indexed: 12/11/2022] Open
Abstract
The recent outbreaks of Ebola virus (EBOV) infections have underlined the impact of the virus as a major threat for human health. Due to the high biosafety classification of EBOV (level 4), basic research is very limited. Therefore, the development of new avenues of thinking to advance quantitative comprehension of the virus and its interaction with the host cells is urgently needed to tackle this lethal disease. Mathematical modeling of the EBOV dynamics can be instrumental to interpret Ebola infection kinetics on quantitative grounds. To the best of our knowledge, a mathematical modeling approach to unravel the interaction between EBOV and the host cells is still missing. In this paper, a mathematical model based on differential equations is used to represent the basic interactions between EBOV and wild-type Vero cells in vitro. Parameter sets that represent infectivity of pathogens are estimated for EBOV infection and compared with influenza virus infection kinetics. The average infecting time of wild-type Vero cells by EBOV is slower than in influenza infection. Simulation results suggest that the slow infecting time of EBOV could be compensated by its efficient replication. This study reveals several identifiability problems and what kind of experiments are necessary to advance the quantification of EBOV infection. A first mathematical approach of EBOV dynamics and the estimation of standard parameters in viral infections kinetics is the key contribution of this work, paving the way for future modeling works on EBOV infection.
Collapse
Affiliation(s)
- Van Kinh Nguyen
- Systems Medicine of Infectious Diseases, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research Braunschweig, Germany
| | - Sebastian C Binder
- Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research Braunschweig, Germany
| | - Alessandro Boianelli
- Systems Medicine of Infectious Diseases, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research Braunschweig, Germany
| | - Michael Meyer-Hermann
- Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research Braunschweig, Germany ; Institute for Biochemistry, Biotechnology and Bioinformatics, Technische Universität Braunschweig Braunschweig, Germany
| | - Esteban A Hernandez-Vargas
- Systems Medicine of Infectious Diseases, Department of Systems Immunology and Braunschweig Integrated Centre of Systems Biology, Helmholtz Centre for Infection Research Braunschweig, Germany
| |
Collapse
|
41
|
McGoff K, Mukherjee S, Pillai N. Statistical inference for dynamical systems: A review. STATISTICS SURVEYS 2015. [DOI: 10.1214/15-ss111] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
42
|
Dattner I, Klaassen CAJ. Optimal rate of direct estimators in systems of ordinary differential equations linear in functions of the parameters. Electron J Stat 2015. [DOI: 10.1214/15-ejs1053] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
43
|
Bhaumik P, Ghosal S. Bayesian two-step estimation in differential equation models. Electron J Stat 2015. [DOI: 10.1214/15-ejs1099] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
|
44
|
Abstract
Existing estimation methods for ordinary differential equation (ODE) models are not applicable to discrete data. The generalized ODE (GODE) model is therefore proposed and investigated for the first time. We develop the likelihood-based parameter estimation and inference methods for GODE models. We propose robust computing algorithms and rigorously investigate the asymptotic properties of the proposed estimator by considering both measurement errors and numerical errors in solving ODEs. The simulation study and application of our methods to an influenza viral dynamics study suggest that the proposed methods have a superior performance in terms of accuracy over the existing ODE model estimation approach and the extended smoothing-based (ESB) method.
Collapse
|
45
|
Progressive contraction of the latent HIV reservoir around a core of less-differentiated CD4⁺ memory T Cells. Nat Commun 2014; 5:5407. [PMID: 25382623 PMCID: PMC4241984 DOI: 10.1038/ncomms6407] [Citation(s) in RCA: 103] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2014] [Accepted: 09/29/2014] [Indexed: 12/26/2022] Open
Abstract
In patients who are receiving prolonged antiretroviral treatment (ART), HIV can persist within a small pool of long-lived resting memory CD4+ T cells infected with integrated latent virus. This latent reservoir involves distinct memory subsets. Here we provide results that suggest a progressive reduction of the size of the blood latent reservoir around a core of less-differentiated memory subsets (central memory and stem cell-like memory (TSCM) CD4+ T cells). This process appears to be driven by the differences in initial sizes and decay rates between latently infected memory subsets. Our results also suggest an extreme stability of the TSCM sub-reservoir, the size of which is directly related to cumulative plasma virus exposure before the onset of ART, stressing the importance of early initiation of effective ART. The presence of these intrinsic dynamics within the latent reservoir may have implications for the design of optimal HIV therapeutic purging strategies. HIV can persist in CD4+ T cells of patients receiving long-term antiretroviral therapy. Here the authors show the presence of intrinsic dynamics that progressively contract the latent HIV reservoir around a core of less-differentiated CD4 T-cell memory subsets.
Collapse
|
46
|
Wu S, Xue H, Wu Y, Wu H. Variable Selection for Sparse High-Dimensional Nonlinear Regression Models by Combining Nonnegative Garrote and Sure Independence Screening. Stat Sin 2014; 24:1365-1387. [PMID: 25170239 PMCID: PMC4142445 DOI: 10.5705/ss.2012.316] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
In many regression problems, the relations between the covariates and the response may be nonlinear. Motivated by the application of reconstructing a gene regulatory network, we consider a sparse high-dimensional additive model with the additive components being some known nonlinear functions with unknown parameters. To identify the subset of important covariates, we propose a new method for simultaneous variable selection and parameter estimation by iteratively combining a large-scale variable screening (the nonlinear independence screening, NLIS) and a moderate-scale model selection (the nonnegative garrote, NNG) for the nonlinear additive regressions. We have shown that the NLIS procedure possesses the sure screening property and it is able to handle problems with non-polynomial dimensionality; and for finite dimension problems, the NNG for the nonlinear additive regressions has selection consistency for the unimportant covariates and also estimation consistency for the parameter estimates of the important covariates. The proposed method is applied to simulated data and a real data example for identifying gene regulations to illustrate its numerical performance.
Collapse
Affiliation(s)
- Shuang Wu
- Department of Biostatistics and Computational Biology, University of Rochester, Rochester, NY
| | - Hongqi Xue
- Department of Biostatistics and Computational Biology, University of Rochester, Rochester, NY
| | - Yichao Wu
- Department of Statistics, North Carolina State University, Raleigh, NC
| | - Hulin Wu
- Department of Biostatistics and Computational Biology, University of Rochester, Rochester, NY
| |
Collapse
|
47
|
Wu H, Lu T, Xue H, Liang H. Sparse Additive Ordinary Differential Equations for Dynamic Gene Regulatory Network Modeling. J Am Stat Assoc 2014; 109:700-716. [PMID: 25061254 DOI: 10.1080/01621459.2013.859617] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The gene regulation network (GRN) is a high-dimensional complex system, which can be represented by various mathematical or statistical models. The ordinary differential equation (ODE) model is one of the popular dynamic GRN models. High-dimensional linear ODE models have been proposed to identify GRNs, but with a limitation of the linear regulation effect assumption. In this article, we propose a sparse additive ODE (SA-ODE) model, coupled with ODE estimation methods and adaptive group LASSO techniques, to model dynamic GRNs that could flexibly deal with nonlinear regulation effects. The asymptotic properties of the proposed method are established and simulation studies are performed to validate the proposed approach. An application example for identifying the nonlinear dynamic GRN of T-cell activation is used to illustrate the usefulness of the proposed method.
Collapse
Affiliation(s)
- Hulin Wu
- Department of Biostatistics and Computational Biology, University of Rochester, School of Medicine and Dentistry, 601 Elmwood Avenue, Box 630, Rochester, NY 14642
| | - Tao Lu
- Department of Epidemiology and Biostatistics, State University of New York, Albany, NY 12144
| | - Hongqi Xue
- Department of Biostatistics and Computational Biology, University of Rochester, School of Medicine and Dentistry, 601 Elmwood Avenue, Box 630, Rochester, NY 14642
| | - Hua Liang
- Department of Statistics, George Washington University, 801 22nd St. NW, Washington, D.C. 20052
| |
Collapse
|
48
|
Chervoneva I, Freydin B, Hipszer B, Apanasovich TV, Joseph JI. Estimation of nonlinear differential equation model for glucose-insulin dynamics in type I diabetic patients using generalized smoothing. Ann Appl Stat 2014; 8:886-904. [PMID: 33833847 DOI: 10.1214/13-aoas706] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/19/2022]
Abstract
In this work, we develop an ordinary differential equations (ODE) model of physiological regulation of glycemia in type 1 diabetes mellitus (T1DM) patients in response to meals and intravenous insulin infusion. Unlike for majority of existing mathematical models of glucose-insulin dynamics, parameters in our model are estimable from a relatively small number of noisy observations of plasma glucose and insulin concentrations. For estimation, we adopt the generalized smoothing estimation of nonlinear dynamic systems of Ramsay et al. (2007). In this framework, the ODE solution is approximated with a penalized spline, where the ODE model is incorporated in the penalty. We propose to optimize the generalized smoothing by using penalty weights that minimize the covariance penalties criterion (Efron, 2004). The covariance penalties criterion provides an estimate of the prediction error for nonlinear estimation rules resulting from nonlinear and/or non-homogeneous ODE models, such as our model of glucose-insulin dynamics. We also propose to select the optimal number and location of knots for B-spline bases used to represent the ODE solution. The results of the small simulation study demonstrate advantages of optimized generalized smoothing in terms of smaller estimation errors for ODE parameters and smaller prediction errors for solutions of differential equations. Using the proposed approach to analyze the glucose and insulin concentration data in T1DM patients we obtained good approximation of global glucose-insulin dynamics and physiologically meaningful parameter estimates.
Collapse
|
49
|
Wu S, Liu ZP, Qiu X, Wu H. Modeling genome-wide dynamic regulatory network in mouse lungs with influenza infection using high-dimensional ordinary differential equations. PLoS One 2014; 9:e95276. [PMID: 24802016 PMCID: PMC4011728 DOI: 10.1371/journal.pone.0095276] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2013] [Accepted: 03/26/2014] [Indexed: 12/20/2022] Open
Abstract
The immune response to viral infection is regulated by an intricate network of many genes and their products. The reverse engineering of gene regulatory networks (GRNs) using mathematical models from time course gene expression data collected after influenza infection is key to our understanding of the mechanisms involved in controlling influenza infection within a host. A five-step pipeline: detection of temporally differentially expressed genes, clustering genes into co-expressed modules, identification of network structure, parameter estimate refinement, and functional enrichment analysis, is developed for reconstructing high-dimensional dynamic GRNs from genome-wide time course gene expression data. Applying the pipeline to the time course gene expression data from influenza-infected mouse lungs, we have identified 20 distinct temporal expression patterns in the differentially expressed genes and constructed a module-based dynamic network using a linear ODE model. Both intra-module and inter-module annotations and regulatory relationships of our inferred network show some interesting findings and are highly consistent with existing knowledge about the immune response in mice after influenza infection. The proposed method is a computationally efficient, data-driven pipeline bridging experimental data, mathematical modeling, and statistical analysis. The application to the influenza infection data elucidates the potentials of our pipeline in providing valuable insights into systematic modeling of complicated biological processes.
Collapse
Affiliation(s)
- Shuang Wu
- Department of Biostatistics and Computational Biology, University of Rochester, Rochester, New York, United States of America
| | - Zhi-Ping Liu
- Department of Biostatistics and Computational Biology, University of Rochester, Rochester, New York, United States of America
| | - Xing Qiu
- Department of Biostatistics and Computational Biology, University of Rochester, Rochester, New York, United States of America
| | - Hulin Wu
- Department of Biostatistics and Computational Biology, University of Rochester, Rochester, New York, United States of America
| |
Collapse
|
50
|
Brunel NJB, Clairon Q, d’Alché-Buc F. Parametric Estimation of Ordinary Differential Equations With Orthogonality Conditions. J Am Stat Assoc 2014. [DOI: 10.1080/01621459.2013.841583] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|