Discovering dynamic patterns from infectious disease data using dynamic mode decomposition.
Ontology highlight
ABSTRACT: Background
The development and application of quantitative methods to understand disease dynamics and plan interventions is becoming increasingly important in the push toward eradication of human infectious diseases, exemplified by the ongoing effort to stop the spread of poliomyelitis.Methods
Dynamic mode decomposition (DMD) is a recently developed method focused on discovering coherent spatial-temporal modes in high-dimensional data collected from complex systems with time dynamics. The algorithm has a number of advantages including a rigorous connection to the analysis of nonlinear systems, an equation-free architecture, and the ability to efficiently handle high-dimensional data.Results
We demonstrate the method on three different infectious disease sets including Google Flu Trends data, pre-vaccination measles in the UK, and paralytic poliomyelitis wild type-1 cases in Nigeria. For each case, we describe the utility of the method for surveillance and resource allocation.Conclusions
We demonstrate how DMD can aid in the analysis of spatial-temporal disease data. DMD is poised to be an effective and efficient computational analysis tool for the study of infectious disease.
SUBMITTER: Proctor JL
PROVIDER: S-EPMC4379984 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Many demographic and ecological processes generate seasonal and other periodicities. Seasonality in infectious disease transmission can result from climatic forces such as temperature and humidity; variation in contact rates as a result of migration or school calendar; or temporary surges in birth rates. Seasonal drivers of acute immunizing infections can also drive longer-term fluctuations.Tensor decomposition has been used in many disciplines to uncover dominant trends in multi-dimensional data. We introduce tensors as a novel method for decomposing oscillatory infectious disease time series.We illustrate the reliability of the method by applying it to simulated data. We then present decompositions of measles data from England and Wales. This paper leverages simulations as well as much-studied data to illustrate the power of tensor decomposition to uncover dominant epidemic signals as well as variation in space and time. We then use tensor decomposition to uncover new findings and demonstrate the potential power of the method for disease incidence data. In particular, we are able to distinguish between annual and biennial signals across locations and shifts in these signals over time.Tensor decomposition is able to isolate variation in disease seasonality as a result of variation in demographic rates. The method allows us to discern variation in the strength of such signals by space and population size. Tensors provide an opportunity for a concise approach to uncovering heterogeneity in disease transmission across space and time in large datasets.
| S-EPMC7756762 | biostudies-literature
Project description:Most neurodegenerative disorders are characterized by progressive loss of neurons throughout the course of disease in the form of specific spatio-temporal patterns. To capture and quantify these coherent patterns across both space and time, traditionally one would either fit a pre-defined model with spatial and temporal parameters or apply analysis in the spatial and temporal domains separately. In this work, we introduce and validate the use of dynamic mode decomposition (DMD), a data-driven multivariate approach, to extract coupled spatio-temporal patterns simultaneously. We apply the method to examine progressive dopaminergic degeneration in 41 patients with Parkinson's disease (PD) using [11C](±)dihydrotetrabenazine (DTBZ) Positron Emission Tomography (PET). DMD decomposed the progressive dopaminergic changes in the putamen into two orthogonal temporal progression curves associated with distinct spatial patterns: 1) an anterior-posterior gradient, the expression of which decreased gradually with disease progression with a higher initial expression in the less affected side; 2) a dorsal-ventral gradient in the less affected side, which was present in early disease stage only. In the caudate, we found a head-tail gradient analogous to the anterior-posterior gradient seen in the putamen; as in the putamen, the expression of this gradient decreased gradually with disease progression with higher expression in the less affected side. Our results with DTBZ PET data show the applicability and relevance of the proposed method for extracting spatio-temporal patterns of neurotransmitter changes due to neurodegeneration. The method is able to decompose known PD-induced dopaminergic denervation into orthogonal (and thus loosely independent) temporal curves, which may be able to reflect and separate either different mechanisms underlying disease progression and disease initiation, or differential involvement of striatal sub-regions at different disease stages, in a completely data driven way. It is expected that this method can be easily extended to other PET tracers and neurodegenerative disorders and may help to elucidate disease mechanisms in more details compared to traditional approaches.
| S-EPMC6948364 | biostudies-literature
Project description:Resting state networks (RSNs) extracted from functional magnetic resonance imaging (fMRI) scans are believed to reflect the intrinsic organization and network structure of brain regions. Most traditional methods for computing RSNs typically assume these functional networks are static throughout the duration of a scan lasting 5-15 min. However, they are known to vary on timescales ranging from seconds to years; in addition, the dynamic properties of RSNs are affected in a wide variety of neurological disorders. Recently, there has been a proliferation of methods for characterizing RSN dynamics, yet it remains a challenge to extract reproducible time-resolved networks. In this paper, we develop a novel method based on dynamic mode decomposition (DMD) to extract networks from short windows of noisy, high-dimensional fMRI data, allowing RSNs from single scans to be resolved robustly at a temporal resolution of seconds. After validating the method on a synthetic dataset, we analyze data from 120 individuals from the Human Connectome Project and show that unsupervised clustering of DMD modes discovers RSNs at both the group (gDMD) and the single subject (sDMD) levels. The gDMD modes closely resemble canonical RSNs. Compared to established methods, sDMD modes capture individualized RSN structure that both better resembles the population RSN and better captures subject-level variation. We further leverage this time-resolved sDMD analysis to infer occupancy and transitions among RSNs with high reproducibility. This automated DMD-based method is a powerful tool to characterize spatial and temporal structures of RSNs in individual subjects.
| S-EPMC6834549 | biostudies-literature
Project description:Theoretical models are typically developed through a deductive process where a researcher formulates a system of dynamic equations from hypothesized mechanisms. Recent advances in algorithmic methods can discover dynamic models inductively-directly from data. Most previous research has tested these methods by rediscovering models from synthetic data generated by the already known model. Here we apply Sparse Identification of Nonlinear Dynamics (SINDy) to discover mechanistic equations for disease dynamics from case notification data for measles, chickenpox, and rubella. The discovered models provide a good qualitative fit to the observed dynamics for all three diseases, However, the SINDy chickenpox model appears to overfit the empirical data, and recovering qualitatively correct rubella dynamics requires using power spectral density in the goodness-of-fit criterion. When SINDy uses a library of second-order functions, the discovered models tend to include mass action incidence and a seasonally varying transmission rate-a common feature of existing epidemiological models for childhood infectious diseases. We also find that the SINDy measles model is capable of out-of-sample prediction of a dynamical regime shift in measles case notification data. These results demonstrate the potential for algorithmic model discovery to enrich scientific understanding by providing a complementary approach to developing theoretical models.
| S-EPMC7184751 | biostudies-literature
Project description:Inertial particles advected by a background flow can show complex structures. We consider inertial particles in a 2D Taylor-Green (TG) flow and characterize particle dynamics as a function of the particle's Stokes number using dynamic mode decomposition (DMD) method from particle image velocimetry (PIV) like-data. We observe the formation of caustic structures and analyze them using DMD to (a) determine the Stokes number of the particles, and (b) estimate the particle Stokes number composition. Our analysis in this idealized flow will provide useful insight to analyze inertial particles in more complex or turbulent flows. We propose that the DMD technique can be used to perform similar analysis on an experimental system.
| S-EPMC8128860 | biostudies-literature
Project description:Accelerating the design of synthetic biological circuits for biosensing requires expanding the currently available genetic toolkit. Although whole-cell biosensors have been successfully engineered and deployed, particularly in applications such as environmental and medical diagnostics, novel sensing applications necessitate the discovery and optimization of novel biosensors. Here we develop a data-driven approach that combines dynamic mode decomposition and observability analysis to extract salient, endogenous genetic sensors for analytes of interest from dynamic gene expression profiles. We show that relatively few dynamic modes are required to capture the transcriptome dynamics. The concept of system observability is then used to rank genes based on their ability to reconstruct the cell state. Genes which contribute to the system's observability are dubbed encoder genes. As a demonstration of our method, we extract, construct, and characterize 15 genetic reporters for the organophosphate malathion in the host bacterium \textit{Pseudomonas fluorescens} SBW25. Furthermore, we demonstrate how the distinct responses of each genetic reporter can be used to construct a single-input, single-output virtual sensor network, approximating a pair of reference trajectories. This library of living malathion sensors can be optimized for use in environmental diagnostics while the developed machine learning tool can be applied to discover genetic reporters in the compendium of host organisms and environmental conditions.
2022-04-21 | GSE200822 | GEO
Project description:Detecting periodicity signals from time-series microarray data is commonly used to facilitate the understanding of the critical roles and underlying mechanisms of regulatory transcriptomes. However, time-series microarray data are noisy. How the temporal data structure affects the performance of periodicity detection has remained elusive. We present a novel method based on empirical mode decomposition (EMD) to examine this effect. We applied EMD to a yeast microarray dataset and extracted a series of intrinsic mode function (IMF) oscillations from the time-series data. Our analysis indicated that many periodically expressed genes might have been under-detected in the original analysis because of interference between decomposed IMF oscillations. By validating a protein complex coexpression analysis, we revealed that 56 genes were newly determined as periodic. We demonstrated that EMD can be used incorporating with existing periodicity detection methods to improve their performance. This approach can be applied to other time-series microarray studies.
| S-EPMC4221108 | biostudies-literature
Project description:Images of the kidneys using dynamic contrast-enhanced magnetic resonance renography (DCE-MRR) contains unwanted complex organ motion due to respiration. This gives rise to motion artefacts that hinder the clinical assessment of kidney function. However, due to the rapid change in contrast agent within the DCE-MR image sequence, commonly used intensity-based image registration techniques are likely to fail. While semi-automated approaches involving human experts are a possible alternative, they pose significant drawbacks including inter-observer variability, and the bottleneck introduced through manual inspection of the multiplicity of images produced during a DCE-MRR study. To address this issue, we present a novel automated, registration-free movement correction approach based on windowed and reconstruction variants of dynamic mode decomposition (WR-DMD). Our proposed method is validated on ten different healthy volunteers' kidney DCE-MRI data sets. The results, using block-matching-block evaluation on the image sequence produced by WR-DMD, show the elimination of 99% of mean motion magnitude when compared to the original data sets, thereby demonstrating the viability of automatic movement correction using WR-DMD.
| S-EPMC7010382 | biostudies-literature
Project description:A new dynamic mode decomposition (DMD) method is introduced for simultaneous system identification and denoising in conjunction with the adoption of an extended Kalman filter algorithm. The present paper explains the extended-Kalman-filter-based DMD (EKFDMD) algorithm which is an online algorithm for dataset for a small number of degree of freedom (DoF). It also illustrates that EKFDMD requires significant numerical resources for many-degree-of-freedom (many-DoF) problems and that the combination with truncated proper orthogonal decomposition (trPOD) helps us to apply the EKFDMD algorithm to many-DoF problems, though it prevents the algorithm from being fully online. The numerical experiments of a noisy dataset with a small number of DoFs illustrate that EKFDMD can estimate eigenvalues better than or as well as the existing algorithms, whereas EKFDMD can also denoise the original dataset online. In particular, EKFDMD performs better than existing algorithms for the case in which system noise is present. The EKFDMD with trPOD, which unfortunately is not fully online, can be successfully applied to many-DoF problems, including a fluid-problem example, and the results reveal the superior performance of system identification and denoising.
| S-EPMC6383872 | biostudies-literature
Project description:Infectious disease occurs when a person is infected by a pathogen from another person or an animal. It is a problem that causes harm at both individual and macro scales. The Korea Center for Disease Control (KCDC) operates a surveillance system to minimize infectious disease contagions. However, in this system, it is difficult to immediately act against infectious disease because of missing and delayed reports. Moreover, infectious disease trends are not known, which means prediction is not easy. This study predicts infectious diseases by optimizing the parameters of deep learning algorithms while considering big data including social media data. The performance of the deep neural network (DNN) and long-short term memory (LSTM) learning models were compared with the autoregressive integrated moving average (ARIMA) when predicting three infectious diseases one week into the future. The results show that the DNN and LSTM models perform better than ARIMA. When predicting chickenpox, the top-10 DNN and LSTM models improved average performance by 24% and 19%, respectively. The DNN model performed stably and the LSTM model was more accurate when infectious disease was spreading. We believe that this study's models can help eliminate reporting delays in existing surveillance systems and, therefore, minimize costs to society.
| S-EPMC6121625 | biostudies-literature